近交是保护生物学与动物育种学中的关键议题,可显著影响家养物种的存活率、繁殖成功率及遗传多样性(GD)。布朗瑞士牛于约一个世纪前引入南美洲,但育种体系衍生的遗传后果在多代中被育种者忽视。本研究基于系谱数据库,旨在阐明引入南美洲的布朗瑞士牛的种群演化、近交水平与GD动态。研究纳入截至2023年的8686头官方登记个体的多代记录,评估了系谱完整性、近交系数(F)、有效种群大小(Ne)、世代间隔(GI)等种群结构参数,同时解析了基因起源概率(有效创始人数量[fe]、有效祖先数量[fa])及GD丧失特征。结果显示,研究群体最大世代数(GMax)为12.33,完整世代数(GCom)为2.37,等效世代数(GEqu)为5.43;平均GI为7.11年,历史群体与近期群体的平均F分别为1.27%和5.60%,近期降至2.78%;F每10年上升0.10%–0.50%,呈持续增加趋势。平均普查有效种群大小(Ne-census)、基于近交增量(ΔFp)的有效种群大小(Ne-ΔFp)及基于等效世代数(GEqu)的有效种群大小(Ne-GEqu)分别为339.50、143.80和95.40,均随时间推移呈下降趋势。基因起源概率衍生参数为fe=55、fa=32,fe/fa比值为1.72,提示该品种可能曾经历偶发种群瓶颈。GD分析显示,总杂合度随时间累计丧失3.35%,但不同时间段存在显著变异。综上,布朗瑞士牛群体内存在显著的GD丧失,主要源于较低的创始人效应及维持中等遗传多样性库所需的个体数量不足,亟需优化管理策略。研究结果强调了监测GD与降低近交影响的必要性,成功的育种规划对完善南美布朗瑞士牛的保护策略与长期可持续生产实践具有重要意义。
《Veterinary Medicine International》论文解读:南美布朗瑞士牛近交、遗传多样性与种群进化研究
研究背景与意义
布朗瑞士牛原产于瑞士阿尔卑斯山区,以极强的环境适应性和 longevity 著称,20世纪上半叶被引入南美洲以提升本地牛种生产性能,目前在厄瓜多尔已成为第二大乳用专门化品种。长期以来,该品种的育种管理忽视了近交与遗传多样性(GD)的动态变化,尤其是20世纪70年代土地改革后,育种群体规模缩减,后续依赖从美国、加拿大大量引进活体、胚胎与冻精,进一步加剧了潜在的遗传风险。欧洲与北美的研究已表明,高强度人工选配导致布朗瑞士牛GD持续下降、有效种群大小(Ne)缩减,但南美洲缺乏系统的系谱层面评估。在此背景下,研究人员针对厄瓜多尔官方登记的布朗瑞士牛群体开展多维度遗传分析,旨在量化近交趋势、GD丧失程度及种群演化规律,为制定可持续育种策略提供依据。
关键技术方法
研究基于厄瓜多尔布朗瑞士协会(ABSE)、美国奶牛育种委员会(CDCB)、瑞士Braunvieh协会及加拿大奶牛网络(CDN)提供的系谱数据,纳入截至2023年出生的8686头个体(含3317头公牛、5369头母牛)。按出生年份划分为5个时段群体:历史全群(1984–2023)、1984–1993、1994–2003、2004–2013、2014–2023,并定义双亲信息完整的参考群体。采用ENDOG v4.8、POPREP、CFC v1.0、GRAIN v2.2等专业软件,计算系谱完整性指标(最大世代数[GMax]、完整世代数[GCom]、等效世代数[GEqu])、世代间隔(GI)、近交系数(F)(Meuwissen法、递归法等)、平均亲缘关系(AR)、共祖系数(C)、非随机交配系数(α)、遗传保护指数(GCI)、有效种群大小(Ne)(普查法、ΔFp法、GEqu法)、基因起源概率(有效创始人数量[fe]、有效祖先数量[fa]、创始人基因组当量[fg])及GD丧失量,并通过Nei遗传距离与Wright固定指数(FIS、FST、FIT)分析种群遗传结构。
研究结果
系谱完整性衍生参数
研究群体的GMax达12.33,GCom为2.37,GEqu为5.43;近期的系谱完整性显著提升,第1至第4代的完整性达0.6以上,表明系谱记录质量随世代推进逐步改善,可支持可靠的GD评估。
世代间隔(GI)
全群平均GI为7.11年,不同时段呈下降趋势:1984–2013年平均GI维持在7.20–7.36年,2014–2023年降至5.60年。分通径看,父系-女儿通径GI始终高于父系-儿子通径,母系-儿子与母系-女儿通径GI相近;近期所有通径GI均显著缩短,介于4.76–7.03年,可能与基因组选择技术的应用缩短了种用个体的世代更替周期有关。
近交、平均亲缘关系、共祖与非随机交配
不同方法估算的近交系数呈现分化:递归法、祖先历史系数(AHC)与Ballou祖先近交系数(fa_Ballou)值最高,Kalinowski法最低,经典Meuwissen法介于两者之间。时序分析显示,F从1984–1993年的0.28%升至2014–2023年的2.78%,近交增量(ΔF)同期从0.10%增至0.50%;2014–2023年群体中92.38%的个体存在近交,高度近交个体占比2.14%。平均亲缘关系(AR)从0.99%升至2.54%,非随机交配系数(α)在近期转为正值,提示存在偏离随机交配的趋势。
基于近交的GD参数
有效种群大小(Ne)呈持续下降趋势:2014–2023年Ne-census约为200,Ne-ΔFp与Ne-GEqu均降至约100,低于FAO建议的长期可持续阈值(Ne≥500)。遗传保护指数(GCI)从2.37升至5.23,反映群体遗传贡献集中度提升。
基因起源概率、祖先贡献与GD
fe/fa比值从历史的1.43升至2014–2023年的1.72,fg/fe比值从0.41降至0.27,提示种群瓶颈与遗传漂变导致的GD丧失加剧。解释50% GD所需的祖先数量从21个减至11个,11个核心祖先贡献了50%的GD,其中公牛TOP ACRES ELEGANT SIMON单独贡献9.90%,母牛IDYL WILD IMPROVER JINX贡献3.82%,反映出少数精英个体的过度使用。
GD丧失
全群GD从1984–1993年的99.01%降至2014–2023年的96.65%,累计丧失3.35%;其中2.44%由遗传漂变导致,0.91%由创始人贡献不均导致,表明遗传漂变是GD丧失的主导因素。
国家间遗传关系
基于15个遗传来源国的分析显示,FIS=0.000992、FST=0.001551、FIT=0.002541,提示亚群内近交水平极低、国家间遗传分化微弱。Nei遗传距离聚类显示,厄瓜多尔群体与美国、瑞士、加拿大的遗传关系最近,印证了其遗传物质主要来源于北美与欧洲的核心育种群体。
讨论与结论总结
讨论部分指出,厄瓜多尔布朗瑞士牛的系谱完整性优于多数南美本地品种,但等效世代数仍低于欧洲核心群,反映出引种后的系谱深度不足。GI的缩短与全球奶牛育种的基因组化趋势一致,但也导致Ne加速下降。近交系数的方法差异源于对未知亲本的处理假设,递归法与祖先近交指标更能反映历史累积效应。fe/fa比值升高验证了引种过程中少数精英个体的“奠基者效应”被放大,而fg/fe比值降低则表明当前群体的实际GD远低于理论预期。GD丧失中遗传漂变占比超70%,与厄瓜多尔近年系谱登记规模缩减直接相关。国家间低分化则源于全球布朗瑞士牛育种的遗传交流高度频繁。
研究结论明确:厄瓜多尔布朗瑞士牛的Ne偏低、F持续上升、GD呈下降趋势,创始人贡献不均与种群规模偏小威胁遗传资源保护。研究人员建议引入低亲缘系数(AR)的新血统、优化配种策略以维持GD,未来需结合基因组工具(如纯合片段[ROH]分析)进一步提升育种管理的精准性。该研究为全球热带与亚热带地区的布朗瑞士牛可持续育种提供了关键的系谱基准数据。