《Mammalian Genome》：Genome-wide association study of mature cow size traits in American Angus cattle
编辑推荐：
为解决安格斯牛成年母牛体重（MWT）、体高（MHT）和体况评分（BCS）等体型性状遗传基础不明确的问题，研究人员开展了针对美国安格斯牛的大规模全基因组关联研究（ssGWAS）。该研究在BTA7、BTA14和BTA20等染色体上发现了与体型性状显著关联的基因组区域，鉴定出包括FBXO32、STC2、CREBRF在内的多个候选基因，并揭示了胰岛素/IGF1R信号、成纤维细胞生长因子受体（FGFR）级联、丝裂原活化蛋白激酶（MAPK）等通路在调控体型和能量平衡中的关键作用。研究结果加深了对安格斯牛体型性状遗传机制的理解，为平衡生长效率与维持成本的育种选择提供了重要依据。
在牛肉生产体系中，提高效率与可持续性是全球日益关注的重点。成年母牛的体重、身高和体况评分不仅是衡量其生产性能的关键指标，更直接影响着牛群的维持饲料成本、繁殖效率和整体经济效益。随着动物蛋白需求的增长，如何在追求更高生长速度的同时，避免因母牛体型过大而导致维持成本显著增加，成为育种实践中亟待平衡的难题。传统上，遗传选育多聚焦于断奶重、周岁重等早期生长性状，而成年母牛的体型性状虽对长期牛群盈利性至关重要，但其背后的具体遗传机制和关键效应基因在很大程度上仍是未知的。
为了深入解析美国安格斯牛成年母牛体型性状的遗传基础，一项题为“Genome-wide association study of mature cow size traits in American Angus cattle”的研究在《Mammalian Genome》发表。这项研究利用美国安格斯协会（AAA）提供的大规模表型与基因型数据，旨在识别与成年母牛体重、身高和体况评分相关的基因组区域、候选基因和多效性变异，并阐明其背后的生物学过程和代谢通路，以期为优化安格斯牛的育种策略提供科学依据。
为开展研究，研究人员主要运用了单步法全基因组关联研究（ssGWAS）和多性状荟萃分析等关键技术。研究数据源自美国安格斯协会的数据库，包含来自数十万头牛的庞大表型记录（如MWT、MHT、BCS）以及超过4.5万头牛的基因型数据。通过严格的质量控制，最终使用了来自45,452头牛的51,410个单核苷酸多态性（SNP）进行关联分析。为检测多效性位点，研究采用了Bolormaa等人提出的多性状荟萃分析方法，对经乔列斯基（Cholesky）变换校正后的性状进行分析。功能分析则利用GALLO等软件包进行基因注释、QTL富集以及通路富集分析，以阐释显著关联区域的生物学意义。
研究结果
成年母牛体型性状的全基因组关联研究
通过ssGWAS，研究在多个染色体上鉴定出与各体型性状显著关联的SNP。
- •
成年母牛体重：在12条染色体上发现59个显著SNP，解释了MWT 6.26%的遗传方差，其中BTA20、BTA7和BTA14贡献最大。这些区域涉及122个基因，包括HAUS4、MMP14、NBEA等。
- •
成年母牛身高：在12条染色体上发现47个显著SNP，解释了MHT 3.07%的遗传方差，BTA20、BTA14和BTA7仍是主要贡献染色体。关联区域包含111个基因，如ADAMTSL3、LAP3、PENK。
- •
体况评分：在14条染色体上发现29个显著SNP，解释了BCS 1.75%的遗传方差，关键区域位于BTA20、BTA7和BTA4。关联区域包含82个基因，如AKAP12、MTHFD1L。
值得注意的是，BTA7、BTA14和BTA20是三个性状共有的关键染色体，存在大量共享的SNP和基因，如位于BTA20的STC2、CREBRF、ERGIC1等。
多效性效应分析
运用多性状荟萃分析方法，研究发现了39个对MWT、MHT和BCS均具有多效性影响的显著SNP，涉及81个基因。BTA20、BTA14和BTA7同样是多效性基因富集的区域。与单性状分析重叠的基因包括FBXO32、DUSP1、STC2、ARRDC3等，同时该方法也识别出一些独特基因，如PECAM1、SMURF2，显示了多性状分析在发现共享遗传机制方面的优势。
QTL鉴定与功能分析
研究将显著关联的基因组区域与牛QTL数据库进行比对，发现它们与大量已报道的QTL存在重叠，这些QTL涉及生产、肉用与胴体、健康、繁殖等多个性状类别。
- •
对于MWT，富集的QTL与体重、平均日增重、代谢体重、胴体重等性状相关。
- •
对于MHT，富集的QTL与体型、体高、体深、臀部宽度等结构性状相关。
- •
对于BCS，富集的QTL与胴体重、眼肌面积、体长等性状相关。
- •
多效性区域富集的QTL同样主要与生产性状（如代谢体重、平均日增重、干物质采食量、胴体重）相关。
通路富集分析进一步揭示了调控这些性状的核心生物学过程：MWT和MHT主要富集于胰岛素/胰岛素样生长因子1受体（IGF1R）信号、成纤维细胞生长因子受体（FGFR）信号级联（如下游的PI3K/AKT、RAF/丝裂原活化蛋白激酶（MAPK）通路）等与生长、代谢和骨骼发育密切相关的通路；而BCS仅富集于RHOD GTP酶循环通路，这可能与细胞骨架重塑和能量代谢有关；多效性分析则突出了G蛋白信号通路作为连接能量平衡、肥胖和生长的共同调节机制。
结论与意义
本研究通过大规模的全基因组关联分析，系统描绘了美国安格斯牛成年母牛体型性状的遗传图谱。研究不仅确认了BTA7、BTA14和BTA20等染色体是影响体型性状的关键基因组热点，与既往关于生长、饲料效率和胴体性状的QTL报告高度重叠，还鉴定出一系列具有多效性潜力的核心候选基因，如FBXO32、STC2、DUSP1、CREBRF等。这些基因在肌肉发育、骨骼生长、脂质代谢和能量调节中扮演重要角色。
更重要的是，通路分析揭示了胰岛素/IGF1R、FGFR/MAPK等高度保守的生长调控通路是影响安格斯牛体型性状的核心生物学途径，而G蛋白信号则可能是协调不同性状间平衡的共同节点。这些发现从分子机制层面解释了为何这些体型性状彼此间存在中等到高度的遗传相关，以及它们如何与饲料效率、胴体品质、繁殖性能等其他经济性状产生内在联系。
此项研究加深了我们对安格斯牛体型性状复杂遗传结构的理解，所识别出的关键基因组区域、候选基因和通路为未来开发功能性分子标记、实施精准育种提供了重要靶点。在实践意义上，这些成果有助于育种者在选育过程中更有效地权衡生长速度与成年母牛维持成本之间的关系，从而培育出在生长效率、生产性能和长期牛群盈利性之间达到更优平衡的安格斯牛品种，最终促进牛肉产业朝着更高效、更可持续的方向发展。
