绕桶赛夸特马表现韧性的遗传指标

《Animal Genetics》：Genetic Indicators of Performance Resilience in Barrel Racing Quarter Horses

　　

韧性（Resilience），定义为在可变条件下维持一致表现的能力，已成为畜牧系统中的一个相关育种目标，但在马属动物群体中仍 largely unexplored。本研究的目的是：(1) 从巴西夸特马（Brazilian Quarter Horses）的纵向绕桶赛时间（Barrel Racing Time, BRT）记录中推导韧性指标；(2) 估计这些指标的遗传参数；(3) 量化其与BRT和累积职业收入（Accumulated Career Earnings, ACE）的遗传关系。数据集包含2010年至2024年间举行的绕桶赛事的BRT记录，涉及10,665匹马。韧性指标是在去除系统效应和持久性马匹与骑手特定表现水平后，围绕预期表现轨迹的短期偏差推导而来。评估了三个指标：中位数绝对偏差（Median Absolute Deviation, MAD）、偏差方差的对数（Logarithm of the Variance of Deviations, LNVAR）和正偏差比例（Proportion of Positive Deviations, PPOS）。韧性指标的后验平均遗传力估计值较低（0.06至0.09），但仍表明存在可利用的加性遗传变异。遗传相关性显示MAD和LNVAR之间近乎完全重叠（0.987），而PPOS捕获了部分独特的遗传信息。所有韧性指标均显示出与ACE的有利遗传关联（平均-0.59）和与BRT的正向关联（平均0.44），表明遗传上更具韧性的马倾向于表现更快并取得更高的经济回报。因此，从纵向表现记录中推导的韧性指标捕获了与表现一致性和经济成功相关的遗传上有意义的变异。这些指标补充了传统表现性状，并提供了一个可扩展的框架，用于将韧性和长期规律性纳入运动马的遗传评估。

**论文解读：绕桶赛夸特马表现韧性的遗传指标**

**研究背景与问题**
在现代动物育种中，韧性（Resilience）被定义为动物在环境挑战下维持生产功能或从应激事件后迅速恢复的能力。该概念在畜牧系统中已被视为重要的育种目标，因其与健康、繁殖和长寿等关键功能性状密切相关，有助于提升系统可持续性和盈利能力。然而，在马属动物群体，尤其是运动马中，韧性相关的遗传评估仍 largely unexplored。虽然已有研究通过异方差或通道化模型（heteroscedastic or canalization models）探索了马匹竞赛等级和形态的变异性遗传基础，但并未关注从纵向运动表现中推导韧性指标。绕桶赛（Barrel Racing）在夸特马（Quarter Horses）中提供了一个理想框架，因其产生大量重复表现记录，且成功取决于速度和规律性，共同影响职业收入。为此，研究人员开展了此项研究，旨在：(1) 从巴西夸特马的纵向绕桶赛时间（Barrel Racing Time, BRT）记录中推导韧性指标；(2) 估计这些指标的遗传参数；(3) 量化韧性指标与BRT和累积职业收入（Accumulated Career Earnings, ACE）的遗传关系。该论文发表于《Animal Genetics》。

**主要技术方法**
研究人员基于2010至2024年间巴西23个州的绕桶赛事数据，样本队列来源为“Sistema de Gerenciamento de Provas”（SGP），经质量控后保留334,508条记录，涉及10,665匹夸特马。关键方法包括：(1) 两步建模策略：先用线性模型预调整BRT数据以去除系统非遗传效应（如当代群体、性别、年龄），再用线性混合模型去除持久性马匹和骑手特定表现水平，获取短期偏差残差；(2) 基于残差计算三个韧性指标：中位数绝对偏差（Median Absolute Deviation, MAD）、偏差方差的对数（Logarithm of the Variance of Deviations, LNVAR）和正偏差比例（Proportion of Positive Deviations, PPOS）；(3) 采用贝叶斯框架下的四性状动物模型估计遗传参数，并计算育种值（EBV）相关系数及偏遗传相关。所有分析在BLUPF90系列程序中进行。

**研究结果**

**3.1 数据**
韧性指标表现出中等离散度。MAD（范围0.019-0.976，变异系数CV=37.1%）、LNVAR（范围-5.810-0.456，CV=32.1%）和PPOS（范围0.100-0.857，CV=19.6%）。ACE原始尺度高度偏斜（范围1-784,169 BRL，CV=200.5%），对数转换后CV降至26.9%。

**3.2 遗传参数**
韧性指标后验平均遗传力较低：MAD为0.072（HPD95%: 0.045-0.106），LNVAR为0.092（0.053-0.143），PPOS为0.063（0.031-0.100）。ACE遗传力为0.081（0.049-0.124）。BRT骑手效应占表型方差0.245（0.235-0.255），遗传力0.207（0.186-0.230），重复率0.397。这表明加性遗传变异虽有限但可被利用。

**3.3 韧性指标与ACE的遗传相关**
MAD与LNVAR遗传相关极强（0.987），近乎完全重叠。PPOS与MAD（0.546）和LNVAR（0.453）呈中等相关。韧性指标与ACE均呈中等负向有利相关：MAD为-0.671，LNVAR为-0.564，PPOS为-0.539。控制BRT后，偏遗传相关略有降低（MAD：-0.604，LNVAR：-0.491，PPOS：-0.376），表明韧性指标与ACE的遗传关联部分独立于竞速能力。

**3.4 EBV相关性**
韧性指标与BRT的EBV相关均为有利正向中等：PPOS最强（0.594），其次为MAD（0.386）和LNVAR（0.336）。ACE与BRT呈负相关（-0.452），说明遗传上更快的马匹收入更高。所有相关均显著（p<0.0001）。

**3.5 表型模式**
根据EBV将马匹分为最韧性（top 200）和最不耐性（bottom 200）两组。韧性组BRT记录分散更窄，平均BRT更低（约18.0s或更低），ACE更高（>15,000 BRL）；不耐性组BRT分散更宽，平均BRT接近或超过18.9s，ACE低于2,500 BRL。表型差异与遗传结果一致。

**讨论与结论**
讨论部分指出，该研究首次从定量遗传学角度证明，表现韧性可作为运动马的评估性状。与畜牧系统类似，韧性指标（MAD、LNVAR、PPOS）虽遗传力低，但存在可选择的加性方差。MAD与LNVAR近乎重合，建议在实际评估中仅选其一；PPOS则提供补充信息，反映不利偏差的频率而非总体变异。韧性指标与ACE和BRT的有利遗传关系表明，遗传上更韧性的马匹兼具更快速度和更高经济回报。偏遗传相关分析进一步确认韧性对经济成功的独立贡献。该方法为将韧性和长期规律性纳入运动马遗传评估提供了可扩展框架，拓展了育种目标从峰值表现向稳定性的转变。翻译结论部分：本研究证明，运动马的表现韧性可通过纵向竞争记录推导的遗传指标进行量化。基于巴西夸特马绕桶赛时间轨迹开发的韧性指标捕获了可检测的加性遗传成分，且MAD与LNVAR高度重叠，PPOS提供独特信息。与竞速表现和职业收入的有利遗传关系表明，遗传上更韧性的马趋向于更快并取得更大经济成功。这些指标提供了补充传统表现性状的遗传上有意义的信息，并为将表现韧性纳入运动马遗传评估提供了可扩展的方法。研究结果为马匹育种计划开辟了新视角，将选择目标从峰值表现扩展到改善韧性与长期稳定性。