1. 领域背景与文献
文献英文标题：Alternative splicing landscape of TGFβ-induced myofibroblast differentiation reveals a conserved 5 alternative splicing event signature；发表期刊：BMC Molecular and Cell Biology；影响因子：未公开；研究领域：细胞生物学-肌成纤维细胞分化调控。
肌成纤维细胞分化是发育生物学及组织修复、纤维化等多种病理生理过程的核心环节，领域发展关键节点包括：20世纪90年代鉴定出转化生长因子β（TGFβ）为肌成纤维细胞分化的核心调控因子，随后发现α-平滑肌肌动蛋白（α-SMA）、I/III型胶原等经典分化标志物，2000年后鉴定出纤连蛋白EDA+剪接体（EDA+FN1）参与分化调控。当前研究热点聚焦于肌成纤维细胞在纤维化疾病中的功能机制，以及寻找特异性调控靶点与诊断标志物。未解决的核心问题包括：肌成纤维细胞分化过程中可变剪接的整体调控图谱尚未系统解析，除EDA+FN1外是否存在其他功能性可变剪接事件仍不明确，可变剪接与基因表达调控的层级关系有待阐明。
针对上述研究空白，本研究通过高通量转录组测序结合双生物信息学分析流程，系统解析TGFβ诱导肌成纤维细胞分化过程中的可变剪接事件，旨在鉴定保守的、具有功能意义的可变剪接特征，为肌成纤维细胞分化的调控机制提供新的视角，同时为纤维化疾病的诊断与治疗提供潜在的生物标志物。
2. 文献综述解析
作者对领域内现有研究按“基因表达调控”和“可变剪接调控”两个核心维度分类评述，系统梳理了肌成纤维细胞分化的经典调控机制与已知可变剪接事件的研究现状，明确了现有研究的局限性与本研究的创新方向。
在基因表达调控方面，现有研究已明确TGFβ-SMAD信号通路是肌成纤维细胞分化的核心通路，通过过表达、沉默实验鉴定出α-SMA、骨膜蛋白（POSTN）、I/III型胶原等多个分化标志物与调控因子，技术方法上多采用细胞模型与实时荧光定量PCR（qRT-PCR）、免疫组化等经典检测手段，优势在于实验体系成熟、结果重复性较好，但存在局限性，如多聚焦于单个基因的功能解析，缺乏对转录组层面整体调控网络的系统研究，且部分研究缺乏多组织、多模型的验证。在可变剪接调控方面，仅EDA+FN1被证实为肌成纤维细胞分化的驱动性剪接体，现有研究多围绕该剪接体的功能展开，技术方法上多采用剪接体特异性抗体与PCR验证，优势在于靶点明确，但局限性在于未系统解析整个转录组的可变剪接图谱，无法全面揭示可变剪接在分化中的调控作用，也未明确可变剪接与基因表达调控的层级关系。
通过对比现有研究的未解决问题，本研究的创新点在于首次采用双生物信息学分析流程对TGFβ诱导的肌成纤维细胞分化进行全转录组可变剪接分析，鉴定出250余个相关剪接体，并验证了5个保守的可变剪接事件特征，突破了现有研究仅聚焦单个剪接体的局限；同时通过基因表达与可变剪接的关联分析，明确了可变剪接调控是独立于基因表达调控的另一层调控机制，为肌成纤维细胞分化的调控网络提供了新的补充，具有重要的学术价值。
3. 研究思路总结与详细解析
本研究的研究目标是解析TGFβ诱导肌成纤维细胞分化过程中的可变剪接调控图谱，鉴定保守的功能性可变剪接事件；核心科学问题是肌成纤维细胞分化过程中是否存在除EDA+FN1外的其他保守可变剪接事件，以及可变剪接与基因表达调控的层级关系；技术路线遵循“细胞模型构建与验证→转录组测序与可变剪接分析→候选事件实验验证→多数据集保守性验证→机制补充探究”的闭环逻辑。
3.1 肌成纤维细胞分化细胞模型构建与验证
实验目的：验证TGFβ处理可诱导原代皮肤成纤维细胞分化为肌成纤维细胞，为后续转录组分析提供可靠的细胞模型。方法细节：采用GM05386原代人皮肤成纤维细胞，培养至70-80%汇合度后，更换为无血清培养基并添加终浓度10 ng/mL的TGFβ处理2天；通过免疫荧光（IF）实验检测α-SMA的表达与应力纤维形成，采用qRT-PCR检测POSTN、ACTA2（α-SMA编码基因）的mRNA表达水平。结果解读：免疫荧光图像显示，TGFβ处理组细胞形成明显的α-SMA应力纤维，
，荧光信号定量结果显示处理组α-SMA表达显著上调（文献未明确提供具体数值，基于图表趋势推测）；qRT-PCR结果显示，处理组POSTN、ACTA2的mRNA表达均显著升高（文献未明确提供具体数值，基于图表趋势推测）。产品关联：实验所用关键产品：Abnova的抗α-平滑肌肌动蛋白抗体（货号H00000059-M02）、Abeomics的TGFβ（货号32-1763-10）、Qiagen的RNeasy mini RNA提取柱（货号74104）等。
3.2 转录组测序与基因表达分析
实验目的：解析TGFβ诱导肌成纤维细胞分化过程中的基因表达变化，验证细胞模型的分化状态，同时为可变剪接分析提供原始数据。方法细节：提取TGFβ处理组与未处理组细胞的总RNA，构建RNA-seq文库并进行Illumina NovaSeq测序；采用DESeq2软件分析差异表达基因，进行基因本体（GO）与KEGG通路富集分析。结果解读：主成分分析显示生物学重复聚类良好，表明实验重复性可靠；共鉴定出453个上调基因与1998个下调基因，上调基因包括多个胶原编码基因（COL4A2、COL11A1等）与细胞外基质（ECM）重塑相关基因（MMP2、LOXL2等），ACTA2基因表达显著上调；
，GO富集分析显示上调基因显著富集于ECM-受体相互作用、黏着斑等肌成纤维细胞分化相关生物学过程，进一步验证了细胞模型的分化状态。产品关联：实验所用关键产品：Illumina NovaSeq测序平台、Qiagen的RNeasy mini RNA提取柱（货号74104）、KAPA Library Quantification Kits（Universal）等。
3.3 全转录组可变剪接事件筛选
实验目的：系统鉴定TGFβ诱导肌成纤维细胞分化过程中的可变剪接事件。方法细节：采用RNomics UdeS pipeline与rMATS两个生物信息学分析工具，对RNA-seq数据进行可变剪接分析，筛选差异可变剪接事件（ASEs）；通过交叉匹配获得两个分析流程共有的ASEs。结果解读：rMATS分析鉴定出251个与肌成纤维细胞分化相关的显著可变剪接体，RNomics UdeS pipeline鉴定出8个相关ASEs；交叉匹配得到6个共有ASEs，包括ACTN1-19A/19B、COL5A1-64A/64B、COL6A3外显子4、FLNA外显子30、EDA+FN1、TPM1-6a/6b；
，
，其中COL6A3的剪接变化为外显子4的包含或缺失，且外显子3在所有样本中均未检测到。产品关联：文献未提及具体生物信息学分析工具的产品信息，领域常规使用rMATS、STAR等开源生物信息学软件。
3.4 候选可变剪接事件的实验验证
实验目的：验证RNA-seq鉴定的候选ASEs在肌成纤维细胞分化中的差异表达。方法细节：采用数字液滴PCR（ddPCR）对6个共有ASEs进行验证，其中FLNA外显子30因短剪接体无法被ddPCR检测，采用等位基因特异性PCR（AS-PCR）；以剪接指数（PSI值）表示剪接体的相对丰度。结果解读：5个ASEs（除EDA+FN1外）的剪接指数变化与RNA-seq结果一致，且在TGFβ处理组中呈现显著差异（P≤0.05，n=3）；EDA+FN1的剪接指数未呈现显著变化，与RNA-seq结果不符，提示其剪接调控可能具有细胞特异性或时间依赖性。产品关联：实验所用关键产品：Bio-Rad的QX200 ddPCR EvaGreen Supermix（货号1864034）、Bio-Rad的QX200液滴生成仪等。
3.5 可变剪接与基因表达调控的层级关系分析
实验目的：探究可变剪接事件对应的宿主基因的表达变化，明确可变剪接调控与基因表达调控的独立性。方法细节：通过DESeq2分析5个验证的ASEs对应的宿主基因（ACTN1、COL5A1、COL6A3、FLNA、TPM1）的表达水平，采用ddPCR验证宿主基因的整体mRNA表达；通过韦恩图分析251个ASEs对应的宿主基因与453个上调基因的重叠情况。结果解读：RNA-seq显示宿主基因有上调趋势，但ddPCR验证未发现显著差异（文献未明确提供具体数值，基于图表趋势推测）；
，251个ASEs对应的202个宿主基因中，仅9个与上调基因重叠，表明可变剪接调控是独立于基因表达调控的另一层调控机制，为肌成纤维细胞分化提供了额外的调控维度。产品关联：实验所用关键产品：Bio-Rad的QX200 ddPCR EvaGreen Supermix（货号1864034）等。
3.6 可变剪接事件的多数据集保守性验证
实验目的：验证5个ASEs特征在不同来源成纤维细胞中的保守性。方法细节：检索并分析3个公开的TGFβ诱导成纤维细胞分化的RNA-seq数据集（1个皮肤成纤维细胞、2个肺成纤维细胞），采用RNomics UdeS pipeline分析5个ASEs的剪接指数变化。结果解读：尽管存在个别例外（如COL5A1在1个肺成纤维细胞数据集、COL6A3在皮肤成纤维细胞数据集中无显著变化），但5个ASEs的剪接变化方向在所有数据集中均一致；
，表明该特征具有跨组织、跨实验条件（处理时间2-20天、TGFβ浓度2-10 ng/mL）的保守性，可作为潜在的通用肌成纤维细胞标志物。产品关联：文献未提及具体数据集分析的产品信息，领域常规使用GEO数据库获取公开转录组数据。
3.7 EDA+FN1剪接调控的补充探究
实验目的：探究TGFβ处理对EDA+FN1剪接的调控机制，解释其未被验证的原因。方法细节：采用免疫荧光实验检测EDA+FN1的蛋白表达，qRT-PCR检测FN1整体mRNA、EDA+FN1与EDA-FN1剪接体的mRNA表达，计算剪接指数（QPSI值）；采用AS-PCR验证剪接指数变化。结果解读：免疫荧光显示TGFβ处理组EDA+FN1蛋白表达上调；
，qRT-PCR显示FN1整体mRNA、EDA+FN1与EDA-FN1剪接体的mRNA均显著上调约5倍（n=3，P<0.01），但剪接指数无显著变化；AS-PCR结果与qRT-PCR一致，表明TGFβ通过上调FN1整体表达而非调控剪接来增加EDA+FN1的水平，解释了其剪接指数未呈现显著变化的原因。产品关联：实验所用关键产品：Santa Cruz的抗EDA+FN抗体（货号sc-59826）、Jackson ImmunoResearch的Donkey anti-Mouse AF488二抗（货号715-545-151）等。
4. Biomarker研究及发现成果
本研究鉴定的肌成纤维细胞分化相关生物标志物为5个保守的可变剪接事件特征（ACTN1-19A/19B、COL5A1-64A/64B、COL6A3外显子4、FLNA外显子30、TPM1-6a/6b），筛选与验证逻辑为：通过RNA-seq双生物信息学流程筛选候选事件→采用ddPCR/AS-PCR进行细胞水平实验验证→利用多组织公开数据集进行保守性验证，形成完整的筛选-验证链条。
该Biomarker来源于原代人皮肤成纤维细胞与公开数据集的肺/皮肤成纤维细胞的转录组数据，覆盖不同组织来源的成纤维细胞模型；验证方法包括数字液滴PCR（ddPCR）、等位基因特异性PCR（AS-PCR）、高通量转录组测序分析，其中ddPCR具有高灵敏度与特异性，可精准定量剪接体的相对丰度；特异性与敏感性方面，5个ASEs在TGFβ诱导的肌成纤维细胞中呈现一致的剪接变化方向，跨组织（皮肤、肺）、跨实验条件（处理时间2-20天、TGFβ浓度2-10 ng/mL）的保守性表明其具有较高的特异性，文献未明确提供敏感性相关数据。
该5个ASEs特征可作为潜在的通用肌成纤维细胞标志物，其功能关联包括：ACTN1、FLNA、TPM1均为肌动蛋白结合蛋白，参与肌成纤维细胞应力纤维的刚性与稳定性调控；COL5A1、COL6A3为ECM蛋白，参与细胞外基质重塑，提示该可变剪接特征可能通过调控细胞骨架与ECM来影响肌成纤维细胞的功能；创新性在于首次系统鉴定并验证了肌成纤维细胞分化过程中除EDA+FN1外的保守可变剪接特征，明确了可变剪接是独立于基因表达的调控层级，为肌成纤维细胞分化的调控网络提供了新的补充；统计学结果显示，5个ASEs在验证实验中均呈现显著的剪接变化（P≤0.05，n=3），多数据集验证中变化方向一致，表明结果具有可靠性。