《Poultry Science》：Transcriptomic analysis identifies candidate genes associated with striated muscle contraction during breast muscle development in meat-type ducks
编辑推荐：
摘要专业翻译
胸肌的发育对肉用鸭的产肉量和肉质至关重要。然而，肉用鸭胸肌发育背后的转录变化仍知之甚少。在本研究中，研究人员整合了表型、组织学和肌源性调控因子（myogenic regulatory factor, MRF）分析，以刻画胸肌的发育轨迹，确定胚胎第
摘要专业翻译
胸肌的发育对肉用鸭的产肉量和肉质至关重要。然而，肉用鸭胸肌发育背后的转录变化仍知之甚少。在本研究中，研究人员整合了表型、组织学和肌源性调控因子（myogenic regulatory factor, MRF）分析，以刻画胸肌的发育轨迹，确定胚胎第19天（E19）为一个关键转折点，其特征是肌肉生长速率达到峰值，且MYOD、MYOG和MYF5的表达达到最大值。随后，研究人员利用15个cDNA文库进行RNA测序（RNA-seq），在三个代表性发育阶段——胚胎第11天（E11）、E19和出壳日（D1）——对胸大肌进行了转录组学分析。两两比较分别在E11 vs. E19、E19 vs. D1和E11 vs. D1中鉴定出5,300、5,167和8,393个差异表达基因（differentially expressed genes, DEGs）。这些DEGs随后进行了基因本体论（Gene Ontology, GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）富集分析。在鉴定出三个发育阶段共有的DEGs后，研究人员将蛋白质-蛋白质相互作用（protein-protein interaction, PPI）网络分析与基因集富集分析（gene set enrichment analysis, GSEA）相结合。该分析揭示了24个参与肉用鸭生长过程中胸肌发育的共有关键基因，并将研究重点聚焦于横纹肌收缩基因集。此外，通过整合CytoHubba排名和GSEA结果，研究人员鉴定出一组与鸭胸肌发育过程中横纹肌收缩相关的六个核心候选基因，即ACTN2、MYH11、MYL1、TNNC1、TNNI1和TNNI2。总之，研究结果描绘了鸭胸肌发育的动态转录组图谱，并鉴定了与横纹肌收缩相关的候选基因，这将有助于未来对这些基因的功能研究以及鸭分子标记的开发。
论文解读
一、研究背景与科学问题
中国是全球最大的肉用鸭生产国和消费国，年出栏量超过35亿只，占全球肉用鸭产量的75%。骨骼肌，尤其是胸肌，是家禽最具经济价值的组成部分之一，肌肉量的积累直接决定产肉量和肉质。在鸟类和哺乳动物中，肌纤维的数量主要在胚胎期确定，出壳或出生后基本保持恒定，出生后的肌肉生长主要依赖于卫星细胞与现有肌纤维融合所驱动的肥大过程。因此，阐明胚胎期骨骼肌发育的分子机制，对于理解肉类生产的生物学基础至关重要。
胚胎肌生成是一个高度有序的生物学过程，涉及成肌细胞融合形成多核肌纤维。肌源性调控因子（MRF）家族是调控肌源性前体细胞分化的核心转录因子，由MYOD、MYOG、MRF4和MYF5四个基因编码。其中，MYOD和MYF5负责胚胎发育早期成肌细胞的定向分化和增殖，而MYOG对终末分化和肌管形成至关重要。值得注意的是，功能性肌肉的成熟不仅依赖于肌纤维的融合，还需要收缩装置的协调组装，这一过程伴随着横纹肌收缩和肌节组织相关基因表达的协调变化。
近年来，RNA测序（RNA-seq）技术极大地推进了畜禽骨骼肌生长和成熟相关转录谱的研究。然而，已有研究主要聚焦于剩余采食量、肉质参数和品种间生长差异等生产性状，对胚胎发育晚期与横纹肌收缩装置功能性组装相关的阶段性特异转录程序仍缺乏系统表征。基于此，本研究旨在鉴定肉用鸭胸肌发育相关的调控基因和信号通路。
二、研究设计与主要技术方法
研究人员以黄山强英鸭业有限公司提供的随机交配群体肉用鸭为实验材料，该群体经过12代闭锁选育。收集E11、E15、E19、E23、E27和D1六个时间点的胚胎样本，分离双侧胸肌进行表型和组织学分析。通过分子性别鉴定（扩增性连锁CHD1基因）筛选雄性胚胎，最终选取E11、E19和D1三个阶段的15个雄性胚胎（每阶段5个生物学重复）进行RNA-seq。主要技术方法包括：利用Illumina HiSeq 4000平台进行转录组测序；采用edgeR包进行差异表达基因（DEGs）鉴定（阈值：|log?FC| ≥ 1且FDR < 0.05）；通过OmicSmart云平台进行GO和KEGG富集分析；利用STRING数据库构建蛋白质-蛋白质相互作用（PPI）网络并经Cytoscape可视化；采用基因集富集分析（GSEA）评估肌肉相关生物学过程；运用CytoHubba的六种拓扑算法（MCC、MNC、Degree、Bottleneck、Closeness和EPC）进行核心基因排序；最后通过定量实时PCR（qPCR）验证RNA-seq结果的可靠性。
三、研究结果
1. 不同胚胎发育阶段胸肌特征及基因表达变化
表型分析显示，胚胎体重在整个发育期持续稳定增长；胸肌重量从E11到E23持续增加，E23至出壳期间保持相对稳定；胸肌比重随时间先升后降，在E19达到峰值。组织学分析表明，肌纤维直径随发育逐渐增大，而纤维密度随时间呈下降趋势。MYOD、MYOG和MYF5三个肌源性调控因子的相对表达量均在E19达到峰值。综合表型、组织学和分子数据，E11、E19和D1被选为RNA-seq分析的三个代表性阶段，分别代表早期发育阶段、胸肌生长速率峰值和肌源性标志物表达最大化的过渡阶段，以及出壳阶段。
2. RNA-seq数据描述性统计
所有样本的clean reads比例超过97%，Q20和Q30质量分数分别大于96.31%和89.96%，GC含量稳定在48.02%至52.91%之间，总比对率超过82.30%，82.77%至86.36%的比对reads位于外显子区域，表明转录组测序具有高特异性和高质量。
3. 差异表达基因分析
主成分分析（PCA）显示组内生物学重复性高、组间分离清晰。E11 vs. E19比较鉴定出5,300个DEGs（1,189个上调，4,111个下调）；E19 vs. D1比较鉴定出5,167个DEGs（1,520个上调，3,647个下调）；E11 vs. D1比较鉴定出8,393个DEGs（6,854个上调，1,539个下调）。三个比较组共有1,549个基因重叠。
4. 差异表达基因的富集分析
GO富集分析显示，E11 vs. E19的DEGs主要富集于系统发育、发育过程、解剖结构发育和多细胞生物过程等生物学过程；E19 vs. D1的DEGs在细胞组分方面主要涉及细胞外区域、细胞外基质、超分子纤维、突触膜和细胞骨架；E11 vs. D1的DEGs在分子功能方面显著富集于结合、细胞骨架蛋白结合和蛋白质结合。KEGG通路分析表明，Wnt、PI3K-Akt和Hippo信号通路在胸大肌发育早期和中期显著差异；ECM-受体相互作用和细胞周期通路在所有三个阶段均显著富集；黏附连接通路在胸大肌发育早期特异性富集。
5. 蛋白质-蛋白质相互作用网络分析
E11 vs. E19的PPI网络包含78个节点和748条边，关键基因包括CAV3、TNNC2和TNNI2；E19 vs. D1的PPI网络包含74个节点和406条边，关键基因为NCK2、RAC3和MAPK11；E11 vs. D1的PPI网络包含88个节点和674条边，关键基因为CTNNB1、TPM2、EGF、TRRAP和TGFB3。三个网络的节点重叠分析鉴定出24个共有基因，包括ACTN2、MYH11、MYL1、TNNC1、TNNI1、TNNI2、MYOZ1、NEB和TNNT2等。功能富集分析显示这些共有基因显著关联于横纹肌收缩、TGF-β信号传导和成纤维细胞生长因子受体信号传导，其中11个基因与横纹肌收缩基因集相关。
6. 基因集富集分析
GSEA鉴定出三个显著富集的基因集：肌肉细胞发育（关键贡献基因为ACTN2和MYH11）、肌肉结构组成（主要由ACTN2、MYH11和MYL1驱动）以及骨骼肌收缩（关键贡献基因包括TNNC1、TNNI1和TNNI2）。这六个基因均映射到横纹肌收缩基因集，进一步支持收缩装置相关过程在胚胎胸肌发育中的重要性。
7. 横纹肌收缩基因集关键候选基因的鉴定
基于24个共有基因构建PPI网络，11个基因与横纹肌收缩基因集显著关联。CytoHubba六种拓扑算法的前10名基因交集产生9个候选核心基因。qPCR验证结果显示这9个基因的表达谱与RNA-seq数据总体一致。表达模式分析表明，ACTN2、MYH11、MYL1和TNNI2的表达水平逐渐升高并在D1达到峰值；TNNC1和TNNI1在E19显著上调后在D1明显下降；SHC4、TNNT2和TGFB2从E11到D1逐渐下调。最终，整合CytoHubba排名和GSEA结果，ACTN2、MYH11、MYL1、TNNC1、TNNI1和TNNI2被确定为核心候选基因，PPI网络显示这六个基因高度互联，形成密集的相互作用模块。
四、讨论与结论
本研究发现E19是肉用鸭胸肌发育的关键转折点，标志着从活跃的成肌细胞增殖阶段向肌管形成和成熟阶段的转变。功能富集分析揭示Wnt、PI3K-Akt、ECM-受体相互作用和Hippo等多条信号通路协同调控胚胎期肌肉发育。六个候选基因中，ACTN2编码α-辅肌动蛋白-2，是骨骼肌和心肌中主要表达的细胞骨架锚定蛋白，通过调节黏着斑信号通路调控肌肉干细胞活化和分化；MYH11编码肌球蛋白重链11，参与肌纤维组成调控；MYL1编码肌球蛋白轻链1，在骨骼肌损伤修复中促进卫星细胞活化和肌纤维重塑；TNNC1编码肌钙蛋白C，是横纹肌收缩的核心调控蛋白；TNNI1（慢肌）和TNNI2（快肌）编码肌钙蛋白I的两种亚型，在无钙条件下抑制肌动蛋白-肌球蛋白相互作用。这些基因共同涵盖了肌节支架构建（ACTN2）、肌球蛋白驱动的力产生（MYH11、MYL1）和肌钙蛋白复合体组装（TNNC1、TNNI1、TNNI2）等多个功能层面，其阶段性表达模式表明收缩装置的重塑与成肌细胞增殖向肌纤维成熟的转变密切相关。
研究结论：本研究描绘了肉用鸭胸肌发育的动态转录组图谱，确定E19为与胸肌生长和肌源性标志物表达变化相关的关键转折点。整合转录组学分析、PPI网络分析和GSEA，突出了横纹肌收缩基因集作为胸肌发育相关分子特征的重要性。ACTN2、MYH11、MYL1、TNNC1、TNNI1和TNNI2六个核心基因被优先确定为与胸肌发育过程中横纹肌收缩相关的候选基因。这些发现为鸭胚胎胸肌发育的生物学基础提供了新见解，未来需进一步研究这些基因在调控肉用鸭胸肌发育中的功能作用。
