《Journal of Dairy Science》：Investigating mechanisms of divergent feed efficiency in dairy cows
编辑推荐：
本研究旨在探究残余采食量（RFI，即观测干物质采食量与预测值之差）与瘤胃微生物组、消化、行为及代谢之间的关联，从而阐明泌乳奶牛RFI差异的潜在机制。研究对象为3个队列共150头经基因分型的荷斯坦奶牛，其中114头奶牛进行了瘤胃微生物群测序，30头奶牛用于肝脏线
本研究旨在探究残余采食量（RFI，即观测干物质采食量与预测值之差）与瘤胃微生物组、消化、行为及代谢之间的关联，从而阐明泌乳奶牛RFI差异的潜在机制。研究对象为3个队列共150头经基因分型的荷斯坦奶牛，其中114头奶牛进行了瘤胃微生物群测序，30头奶牛用于肝脏线粒体呼吸分析。研究人员根据表型RFI（pRFI）或基因组RFI（gRFI）将奶牛按RFI排序并分为四分位组（Q1为最高效组，Q4为最低效组）。统计模型将线性和二次项RFI作为连续解释变量进行分析。结果表明，无论基于表型还是基因组分类，效率提高（即RFI值从大到小）均与干物质采食量（DMI）降低相关——pRFI的Q4与Q1相差3.0 kg/d，gRFI相差1.9 kg/d，且未影响能量校正乳产量（ECM）或体组织储备。这些差异使pRFI和gRFI极端组的饲料转化率分别提高了200 g ECM/kg DMI和100 g ECM/kg DMI。pRFI和gRFI均与乳脂中脂肪酸（FA）谱相关：随着效率提高，从头合成及混合来源FA比例下降，而预先形成FA（尤其是单不饱和FA）比例上升。此外，pRFI与gRFI呈中等程度相关（r=0.48），两种分组方法的奶牛排名一致性较高（ρ=0.44）。RFI降低与总反刍时间减少有关，但每千克DMI的反刍时间增加——pRFI和gRFI极端组间分别相差2.0 min/kg和1.7 min/kg。表型和基因组高效奶牛的瘤胃微生物α多样性（基于逆辛普森指数）较低。基于pRFI和gRFI的Q1与Q4比较发现，共有57个扩增子序列变异组（ASVG）丰度存在差异，其中普雷沃氏菌属（Prevotella）和琥珀酸弧菌科（Succinivibrionaceae）为两类RFI分类所共有。表型和基因组效率提高均与瘤胃氨氮（NH3-N）浓度升高相关。基因组高效奶牛尽管采食较少，但其瘤胃pH值倾向于降低（gQ1至gQ4：6.42 vs. 6.47 vs. 6.43 vs. 6.53）。pRFI降低与微生物氮产量减少相关，但微生物氮占氮摄入量的比例有增加趋势。综上，表型RFI与基因组RFI具有中等程度的一致性，与该性状的估计遗传力相符；改善饲料效率的机制与瘤胃微生物群和发酵的差异，以及每千克DM反刍时间的增加有关，而与全消化道消化率或肝脏线粒体呼吸无关。
论文解读文章
1. 研究背景与科学问题
饲料成本是奶牛养殖业最大的经济支出，同时饲料利用效率也直接影响畜牧业的可持续发展。残余采食量（RFI）作为衡量饲料效率的关键指标，定义为实际干物质采食量（DMI）与基于产奶能量、代谢体重及体组织变化预测DMI之间的差值。低RFI（负值）代表更高的饲料效率，已被纳入美国等多个国家荷斯坦牛的育种计划。然而，尽管RFI在遗传选择中广泛应用，其背后的生物学机制仍不完全清楚。已有研究表明，RFI差异可能与瘤胃微生物组成、营养物质消化率、采食行为及线粒体能量代谢等多种生理过程有关，但这些因素的相对贡献以及表型RFI（pRFI）与基因组RFI（gRFI）所反映的生物学机制是否一致尚缺乏系统研究。因此，本研究的目的是通过整合分析DMI、瘤胃发酵、细胞能量代谢指标以及基因组预测，全面探究泌乳奶牛饲料效率差异的生物学基础。
2. 研究方法与技术路线
本研究为前瞻性队列研究，共纳入3个队列150头泌乳荷斯坦奶牛。所有奶牛均进行基因分型（中密度SNP芯片，约4万个标记），并由美国奶牛育种委员会进行基因型填充。研究人员首先通过线性混合效应模型计算每头奶牛的预测DMI，进而获得pRFI（实测值与预测值之差）；gRFI则直接采用基因组育种值。根据RFI值将奶牛分为四分位组（Q1最有效至Q4最无效）。关键技术方法包括：（1）瘤胃液采集与16S rRNA基因（V1-V3区）扩增子测序，用于分析微生物α/β多样性和差异丰度；（2）气相色谱法测定瘤胃短链脂肪酸（SCFA）浓度；（3）酚-次氯酸盐比色法测定瘤胃氨氮（NH3-N）浓度；（4）使用尿嘌呤衍生物法估算微生物氮流量；（5）以未消化中性洗涤纤维（uNDF）为内源标记测定全消化道表观消化率；（6）高分辨率呼吸测量法（Oroboros O2K）评估肝脏分离线粒体的氧消耗速率（OCR）及偶联效率；（7）自动监测设备记录反刍时间和活动量。
3. 研究结果
表型与基因组RFI的相关性：pRFI与gRFI呈中等程度正相关（r=0.48，ρ=0.44），47%的pQ1奶牛同时属于gQ1，56%的pQ4奶牛属于gQ4，表明两种分类方法具有一定一致性。
DMI、产奶量与饲料转化率：pRFI与DMI呈二次相关（P=0.05），gRFI呈线性相关（P=0.002）。pQ1比pQ4少采食3.0 kg/d，gQ1比gQ4少采食1.9 kg/d，且不影响ECM产量或体组织储备。pRFI极端组饲料转化率（ECM/DMI）提高0.20 kg/kg，gRFI极端组提高0.10 kg/kg。
乳脂肪酸组成：pRFI和gRFI降低均与乳脂中从头合成FA（<16碳）和混合FA（16碳）比例下降、预先形成FA（>16碳）比例上升相关，尤其单不饱和FA（MUFA）显著增加。gRFI高效组还表现出牛奶脂肪抑制型共轭亚油酸（CLA）产量增加的趋势。
反刍行为与活动量：pRFI降低与每日总反刍时间减少（P=0.01）、但每千克DMI反刍时间增加（P=0.002）相关，极端组相差2.0 min/kg DMI。gRFI也呈现类似趋势。两种RFI分类均与奶牛活动量（步数/小时）无显著关联。
瘤胃微生物群落：pRFI降低与Chao1丰富度（P=0.04）和逆辛普森多样性指数（P=0.05）下降相关；gRFI降低同样导致逆辛普森指数下降（P=0.03）。β多样性分析显示，pRFI和gRFI的Q1与Q4组间瘤胃微生物组成均存在显著差异（PERMANOVA，P<0.001）。共鉴定出57个差异丰度ASVG，其中普雷沃氏菌属（Prevotella）和琥珀酸弧菌科（Succinivibrionaceae）为pRFI和gRFI共同差异类群。
瘤胃发酵参数：pRFI和gRFI均与单个SCFA摩尔比例或浓度无显著关联，但pRFI与总SCFA呈二次趋势（P=0.09）。pRFI和gRFI降低均与瘤胃NH3-N浓度显著升高相关（P<0.01）。gRFI高效组瘤胃pH值有降低趋势（P=0.09）。pRFI降低与微生物氮产量减少（P=0.03）相关，但微生物氮占氮摄入量的比例有增加趋势（P=0.07）。
全消化道消化率：pRFI和gRFI均与DM、OM、NDF、CP和淀粉的全消化道表观消化率无显著关联。
肝脏线粒体呼吸：pRFI和gRFI均与肝脏线粒体的漏呼吸、磷酸化呼吸（PCIc、PCIc+II）、ATP偶联呼吸、偶联效率及最大解偶联呼吸（ECI+II）等指标无显著关联。Pearson相关性分析也未发现线粒体能量代谢指标与RFI之间存在显著关系。
4. 讨论与结论
本研究表明，表型RFI与基因组RFI具有中等程度的一致性，这与该性状的估计遗传力相符。高效奶牛（低RFI）在不牺牲产奶量和体况的前提下显著降低了采食量，其机制主要与瘤胃微生物群落结构和功能的变化有关，而非全消化道消化率或肝脏线粒体呼吸效率的差异。具体而言，高效奶牛瘤胃微生物α多样性更低，提示可能存在更特化的微生物群落；瘤胃NH3-N浓度升高可能反映了更活跃的蛋白质降解与氮循环；每千克DMI反刍时间增加有助于更充分的物理消化和唾液缓冲。乳脂中MUFA比例的上升进一步支持了瘤胃氢化途径的改变。这些发现表明，改善饲料效率的核心机制在于瘤胃层面的养分转化效率提升，而非宿主组织的能量代谢优化。研究结果为通过微生物标志物间接筛选高效奶牛提供了理论基础，并强调了在遗传选育中整合表型和基因组信息的价值。
