《Journal of Dairy Science》：Genomic Dissection of Genetic Correlation Between Stillbirth and Gestation Length in German Holstein Cows
编辑推荐：
本研究针对奶牛死胎（SB）这一关乎动物福利和经济的重要问题，探索其与妊娠期长度（GL）的复杂遗传关系。研究人员利用高密度和全基因组测序基因型数据，通过区域遗传力定位（RHM）及局部遗传相关性分析，在基因组水平揭示了SB与GL之间的拮抗性局部遗传关联，并鉴定出ZNF415、SLC43A2等关键候选基因及相关生物学通路。该研究为深入理解奶牛死胎的遗传机制、优化育种策略以改善动物福利和生产效益提供了新的见解。
在奶牛养殖业中，小牛出生时或出生后48小时内死亡，即死胎（Stillbirth, SB），不仅是一个严峻的动物福利问题，也造成了巨大的经济损失。尽管通过育种和管理手段，德国荷斯坦头胎牛的SB率在过去20年间从11.5%下降至6.9%，但进一步降低的空间依然存在。SB的风险受到犊牛性别、母牛胎次、难产和双胎等多种因素影响。其中，妊娠期长度（Gestation Length, GL）与SB之间存在着一种非线性的有趣关系：无论是缩短还是延长的妊娠期，似乎都与更高的死胎率相关。然而，这种关联背后的遗传机制却像一团迷雾，尚未被完全揭开。理解SB与GL之间共享的遗传效应，对于通过基因组选择进行更精准的育种、从根本上改善奶牛福利和生产效率至关重要。
传统上，研究人员通过估算全基因组遗传相关性（global genetic correlation）来评估性状间的共同遗传背景。但这种方法就像取一个平均数，可能会掩盖基因组不同区域（local）存在的、方向相反的遗传效应。例如，在基因组某个区域，延长GL的基因可能同时增加SB风险（负相关）；而在另一个区域，同样的基因效应却可能降低SB风险（正相关）。这种“拮抗”效应在全基因组层面相互抵消，从而无法被察觉。为了深入“解剖”这种复杂的遗传关系，一项发表在《Journal of Dairy Science》上的研究应运而生，它采用了更为精细的基因组分析工具，旨在绘制出SB与GL之间局部遗传相关性的详细图谱。
该研究由德国霍恩海姆大学动物科学研究所的团队主导。为了开展这项研究，研究人员整合了大规模基因型与表型数据。他们获得了34,497头德国荷斯坦奶牛的基因型数据，这些数据最初来源于低密度或50K SNP芯片，并经过基因型填充（imputation）统一至50K水平，随后进一步填充至高密度（High-Density, HD）和全基因组测序（Whole-Genome Sequencing, WGS）水平，确保了分析的广度和精度。表型数据则采用了官方育种值估计中获得的去回归育种值（Deregressed Proofs, DRP），涉及死胎的直接效应（SBd）和母体效应（SBm）、妊娠期长度的直接效应（第一胎GL_d1，第二胎GL_d2）和母体效应（GL_m1, GL_m2），以及作为生产性能参照的乳蛋白产量（PY）。研究主要运用了三种核心分析方法：首先，利用区域遗传力定位（Regional Heritability Mapping, RHM）来识别对性状方差贡献显著的基因组区域；其次，通过双基因组限制性最大似然法（bivariate GREML）估算性状间的全基因组遗传相关性；最后，采用相关性扫描（Correlation Scan, CS）方法，在滑动窗口框架下精细估算局部遗传相关性。在找到显著区域后，研究还利用变异效应预测器（Variant Effect Predictor, VEP）进行基因注释，并通过基因本体（Gene Ontology, GO）富集分析探索潜在的生物学通路。
研究结果
区域遗传力定位
RHM分析成功鉴定出多个对目标性状方差有显著贡献的基因组区域。对于妊娠期长度的直接效应，在第一胎（GL_d1）和第二胎（GL_d2）中均发现了两个相同的显著区域，分别位于Bos taurusautosome 18（BTA18）和19（BTA19）上。对于妊娠期长度的母体效应，仅在第二胎（GL_m2）于BTA3上发现一个显著区域。对于死胎的母体效应（SBm），在BTA9上发现一个显著区域；而死胎的直接效应（SBd）未发现通过严格校正的显著区域。作为对照，乳蛋白产量（PY）在BTA5、6、14、19和20上发现了五个显著区域，这与该性状较高的遗传力相符，也验证了分析方法的效力。
全局遗传相关性
全基因组水平的遗传相关性分析显示，大多数性状对之间的相关性接近于零且不显著。唯一一个显著的相关性出现在死胎直接效应（SBd）与第一胎妊娠期长度直接效应（GL_d1）之间，为-0.17（标准误=0.04），表明两者存在轻微的负向遗传关联，即从全基因组平均水平看，较短的GL可能与较低的SB风险略相关。其他组合，如SB与PY之间，未检测到显著相关性，提示针对产奶量的选择在遗传上并不必然对抗死胎率的降低。
相关性扫描（局部遗传相关性）
局部遗传相关性分析揭示了在全基因组分析中被掩盖的细节。对于SBd和GL_d1这一对性状，在RHM发现的显著区域内部，CS分析检测到了方向相反的局部相关性：在BTA18上，两者呈负相关；而在BTA19上，则呈正相关。这种“拮抗”模式在SBd和GL_d2中也存在，但信号较弱。对于SBm与GL_m1/GL_m2，在SBm的显著区域BTA9上，检测到了强烈的正局部相关性。
变异注释与候选基因
研究人员对RHM和CS均显著的基因组区域内的遗传变异进行了注释，以寻找潜在的候选基因。在SBd与GL_d1呈负相关的BTA18区域，发现了ZNF415（锌指蛋白415）基因，该基因此前已被报道与公牛的女儿难产率和死胎率相关。在呈正相关的BTA19区域，则发现了SLC43A2基因，它编码一个氨基酸转运蛋白。在SBm的显著区域BTA9，发现了AKAP12和SCAF8等基因，其中SCAF8的高影响移码变异曾与导致严重难产的先天性畸形相关联。
基因本体富集分析
GO富集分析为理解SBd与GL_d1的遗传互作提供了功能层面的见解。最显著的富集集群与“有机酸和氨基酸转运”通路相关，涉及32个基因，其中包括SLC43A2和另一个氨基酸转运蛋白基因SLC1A5。这提示营养物质（特别是氨基酸）的跨膜运输在调节妊娠期长度和影响犊牛存活方面可能扮演关键角色。另一个较小的集群与“（后）突触结构”相关。
研究结论与意义
这项研究通过整合区域遗传力定位、局部遗传相关性扫描以及生物信息学注释，深入解析了德国荷斯坦牛死胎与妊娠期长度之间复杂的遗传架构。主要结论与意义包括：
- 1.
揭示了拮抗性局部遗传相关的存在：研究首次在基因组尺度上证实，SB与GL之间的遗传关联并非均匀一致。在BTA18上，延长GL的等位基因可能增加SB风险（负相关）；而在BTA19上，同样的等位基因却可能降低SB风险（正相关）。这种方向相反的局部效应在全基因组相关性中被平均化而隐藏，凸显了局部分析对于理解复杂性状互作的重要性。
- 2.
鉴定了新的候选基因和通路：研究将ZNF415和SLC43A2锁定为影响SB和GL互作的关键候选基因。GO富集分析进一步将“氨基酸转运”通路置于核心位置，这不仅为已知的SLC43A2功能提供了支持，也指向了围产期营养（如蛋氨酸）补充可能通过影响这些转运体来优化妊娠期长度、改善犊牛存活率的潜在干预途径。
- 3.
阐明了性状效应的特异性：研究发现显著的遗传相关性主要局限于死胎的直接效应（SBd）与第一胎妊娠期长度的直接效应（GL_d1）之间。母体效应之间的相关性很弱，且生产性状（PY）与SB之间无显著遗传冲突。这为育种实践提供了精细的指导：针对降低SB的育种选择，可能需要特别关注影响首胎妊娠期长度的犊牛自身基因，而无需过分担忧其对产奶量的负面影响。
- 4.
提供了方法论范式和未来研究方向：本研究成功展示了结合RHM（定位性状相关区域）与CS（精细绘制局部相关性）的分析流程的有效性。研究也指出，许多已鉴定的候选基因功能未知，未来的基因表达分析和功能实验将有助于阐明其具体机制，从而最终为通过基因组选择、管理干预或营养策略综合降低奶牛死胎率，提升动物福利和牧场经济效益提供坚实的科学基础。
