牛科动物蒙大拿（Montana）合成群体（composite）品种组成与种群结构的基因组分析

《Mammalian Genome》：Genomic analysis of breed composition and population structure in Montana composite cattle

　　

蒙大拿（Montana）合成牛于巴西由瘤牛（Bos indicus）与欧洲牛（Bos taurus）杂交培育，并被划分为四个生物学类型：瘤牛型（Zebu，N）、适应型陶牛（adapted taurine，A）、英国陶牛（British taurine，B）和大陆陶牛（continental taurine，C）。研究人员旨在通过主成分分析（principal component analysis，PCA）、混合分析（admixture analysis）和赖特固定指数（Wright’s FST statistic）利用基因组数据分析该合成群体的遗传多样性与种群结构。PCA显示瘤牛与陶牛群体间存在明显分离，蒙大拿个体分布于中间区域；前两个主成分分别解释总变异的69.48%与3.45%。监督式混合估计表明陶牛贡献占优势，A型在K=4、9、11时分别占34.47%、52.64%与51.71%。提高祖先分辨率可细化各创始品种贡献，但不改变陶牛祖先总体占优。系谱与基因组所得品种比例多数生物学类型和祖先模型下差异显著（P<0.001），说明实现品种组成偏离理论预期。遗传分化估计证实瘤牛与陶牛组间分歧更大，共享祖先群体间距离更小。研究人员识别出该合成群体与部分创始品种（尤指Belmont Red、Senepol、Tuli）间特定关系。结果表明蒙大拿合成牛具有复杂基因组结构，基因组祖先随分辨率变化且不同于系谱预期。

研究背景方面，生产动物长期经受自然与人工选择，杂交可组合不同品种互补性状。蒙大拿合成牛1994年于巴西由四类生物学类型牛种杂交形成，目标为聚合抗逆性、性早熟、繁殖效率、胴体质量及热带亚热带适应性。该计划采用开放体系，各类型内具体品种不固定，且认定蒙大拿合成牛需满足至少三品种、A型最低12.5%、N+A合计最低25%、N最高37.5%等系谱比例要求。传统系谱方法不能反映孟德尔分离带来的变异，而基因组数据可更准确估计遗传参数、降低系谱误差并体现孟德尔抽样。因此，研究人员认为有必要在品种水平用基因组手段分析蒙大拿合成牛的种群结构、品种祖先与遗传分化，以指导育种规划与遗传资源保护。该研究论文发表于《Mammalian Genome》。

主要技术方法方面，研究人员基于371,596头系谱数据库识别主要贡献品种，并整合3,892头蒙大拿合成牛与15个创始品种共19,947头动物基因型；样本来自巴西戈亚斯、南马托格罗索、米纳斯吉拉斯、南里奥格兰德、圣保罗及乌拉圭共18个商品群。基因型经跨平台协调至ARS-UCD1.2参考基因组，取两芯片共有17,619个常染色体SNP，再用PREGSF90做质量控制和PLINK相关流程，最终保留约16,081个SNP。研究人员采用VanRaden法建基因组关系矩阵做PCA；用ADMIXTURE 1.3做无监督与监督混合分析，K=4、9、11分别对应生物学类型、部分细分与英国陶牛单品种分辨率，并将系谱与基因组比例用配对Wilcoxon检验与Benjamini–Hochberg校正比较；用PLINK 1.9按Weir and Cockerham法估计两两加权F_ST。

研究结果方面，主成分分析（Principal Component Analysis，PCA）研究人员通过基于SNP的基因组关系矩阵做三种情景PCA：全群体、去Nellore、去Nellore与蒙大拿。结果显示PC1解释69.48%变异并分离Bos indicus与Bos taurus，蒙大拿个体居中间且更近陶牛；去大样本Nellore后PC1、PC2解释为30.80%与11.87%；再去蒙大拿后欧洲陶牛、热带适应陶牛与瘤牛形成更清晰聚类，证明群体结构随参照集改变但陶牛亲和趋势稳定。

混合分析（admixture analysis）研究人员用监督ADMIXTURE估计祖先比例。K=4时A型34.47%、C型31.19%、B型21.93%、N型12.40%；K=9时Bonsmara 15.21%、Senepol 14.44%、英国陶牛组13.23%、Tuli 12.05%、Limousin 11.13%、归组瘤牛（N1）7.61%、Boran 6.54%；K=11时Belmont Red 17.22%、Senepol 13.55%、Red Angus 11.72%、Bonsmara 11.01%、Limousin 10.78%，Aberdeen Angus与Hereford仅3.00%与2.00%。地区分层显示南里奥格兰德与乌拉圭英国祖因中Red Angus占优，南马托格罗索以Limousin为大陆陶牛主源，戈亚斯瘤牛贡献最高。结论为K值决定分辨率，宽类群会掩盖创始品种真实贡献。

系谱对基因组品种组成研究人员对3,892头配对比较。K=4时N、A、C差异显著（错误发现率校正P<0.001），B无差异；K=9与K=11时四型均显著差异。系谱高估瘤牛贡献，低估Tuli与Belmont Red等适应陶牛贡献。结论为实现基因组组成受重组、选配与模型分辨率共同影响，不等于系谱理论比例。

遗传关系与分化（Genetic relationship and differentiation）研究人员用两两F_ST评估。生物学类型间N与陶牛型F_ST为0.168–0.246；蒙大拿与A型仅0.022、与N型0.115。品种水平Sindi–Nellore为0.013，Nellore–Hereford达0.316；蒙大拿与Belmont Red 0.034、与Bonsmara 0.147。结论为同生物学类型内分化低、瘤牛—陶牛间为高分化，蒙大拿与Belmont Red、Senepol、Tuli共享祖先最近。

讨论部分总结方面，研究人员指出PCA中热带适应陶牛近蒙大拿源于Sanga组历史混合，Senepol近蒙大拿源于热带适应选择，Brahman近蒙大拿源于陶牛渐渗。混合分析显示低K把异源创始品种归并到大陆陶牛等宽类群，高K揭示Belmont Red、Bonsmara、Senepol、Tuli真实占比；系谱—基因组差异既来自实现遗传也来自祖先模型分辨率。F_ST与PCA口径不同：前者量化等位基因频率分歧，后者保协方差结构，故Bonsmara在PCA近而在F_ST较远可由其合成起源与亚结构解释。研究人员强调复合群体祖先推断具模型依赖性，未来应评估不同分辨率基因组比例对遗传参数与基因组预测精度影响。

结论部分翻译方面，综合研究表明，蒙大拿合成牛呈现复杂遗传结构，瘤牛与陶牛基因素材间存在清晰分化，当前基因组组成以陶牛祖先占优。创始品种贡献不均反映其多品种起源与不同遗传背景。推断基因组组成随祖先模型变化，说明不同分辨率提供该合成群体遗传架构的补充信息。系谱与基因组估计差异表明实现品种组成可能偏离预期值。总体上，结果对蒙大拿合成牛的基因组组成与种群结构提供了详细表征。