《Animals》：Mitochondrial Genome Variation in Breast Muscle of Tibetan Fragrant Chickens Across an Altitudinal Gradient and Its Association with Meat Quality
编辑推荐：
该研究旨在表征藏香鸡沿海拔梯度的线粒体（mitochondrial，MT）基因组，并将其与肉品质关联。选取18只150日龄雄性鸡，来自三组（LX，1500 m；MX，2800 m；HX，3500 m；每组n = 6；每组饲养30只中随机选6只），进行全基因组重
该研究旨在表征藏香鸡沿海拔梯度的线粒体（mitochondrial，MT）基因组，并将其与肉品质关联。选取18只150日龄雄性鸡，来自三组（LX，1500 m；MX，2800 m；HX，3500 m；每组n = 6；每组饲养30只中随机选6只），进行全基因组重测序（whole-genome resequencing；参考：Gallus gallus GCA_000002315.5）。鉴定90个MT变异（85个多态性；转换/颠换（transition/transversion）83:7），包括58个编码区（42个同义、16个非同义）、20个D-loop/非编码区、10个核糖体RNA（ribosomal RNA，rRNA）和2个转运RNA（transfer RNA，tRNA）位点；未发现MT插入/缺失（insertion/deletion，InDel）或结构变异（structural variant，SV）。18只鸡聚类为三个母系世系（maternal lineage；A、B、C），仅部分按海拔分层（A在低海拔更常见；B在中、高海拔）。与12个Gallus线粒体基因组构建的邻接（neighbour-joining，NJ）树将该群体置于接近Gallus gallus gallus和中国地方品种的位置。检测到15个异质性（heteroplasmy）位点，其中16S rRNA nt3946（G/C/A）在所有鸡中均异质。相对线粒体DNA（mtDNA）拷贝数组间无差异（p = 0.52）。肉品质随海拔变化（2800 m蛋白最高、脂肪最低；3500 m肌苷酸（inosine monophosphate，IMP）最高），但不随母系世系变化，提示观察到的表型差异主要与饲养海拔而非母系线粒体世系相关，且MT基因组贡献的是现存遗传变异，而非这些差异的直接功能基础。
研究背景：藏香鸡是在青藏高原经长期自然选择和家庭驯化形成的地方品种，可在约2200–4100 m完成生活周期；木里藏族自治县位于横断山东缘，海拔跨度超过4400 m，长期地理隔离可能形成不同于其他藏鸡的特征。线粒体是需氧呼吸核心细胞器，线粒体DNA（mitochondrial DNA，mtDNA）母系遗传、拷贝数高且常见异质性，编码13个氧化磷酸化（oxidative phosphorylation，OXPHOS）亚基、2个rRNA、22个tRNA和控制区（D-loop）。低氧下呼吸链效率和能量供应直接影响适应能力，呼吸链基因是低氧适应候选位点。鸡mtDNA单倍群A–E反映母系起源、地理和人工选择；藏鸡MT-CO3单倍型及MT-COI、ATP6变异曾与高海拔适应关联。但多数藏鸡线粒体研究聚焦单基因或D-loop及藏鸡群体，尚缺乏沿海拔梯度系统整合序列变异、母系结构、异质性和拷贝数，木里群体线粒体特征及其与肉品质关系未表征。该研究分析三海拔（1500、2800、3500 m）藏香鸡胸肌全基因组重测序数据，目的为表征MT变异谱、定义并系统发育定位母系世系、描述异质性、估计相对mtDNA拷贝数、关联表型和肉品质并检验是否随海拔变化；同时为监测杂交导致的遗传侵蚀提供母系基线，并评估肌肉线粒体特征对氧化能力和肉品质的潜在影响。比较基因组研究提示低氧适应涉及缺氧诱导因子（hypoxia-inducible factor，HIF）通路和EPAS1等基因，结构变异如多态转座元件也可能参与。研究结论：木里群体MT变异以点突变为主、结构保守，三个母系仅弱海拔结构，异质性常见，拷贝数稳定；肉品质随海拔而非母系变化，提示高海拔适应主要由核基因和环境驱动，MT基因组主要贡献现存遗传多样性和母系标记。意义在于为品种认证、追溯和保护提供标记，并提示优质藏香鸡生产应考虑饲养海拔。
研究方法：18只150日龄雄性藏香鸡来自四川木里县示范农场；日龄雏鸡在2800 m孵化至30日龄后转至LX 1500 m、MX 2800 m、HX 3500 m饲养至150日龄，每组30只中随机选6只。采集胸大肌，全基因组重测序（Illumina NovaSeq 6000，PE150，平均深度约4×，参考Gallus gallus GCA_000002315.5），BWA比对、SAMtools变体调用并提取“MT”变异，ANNOVAR注释；UPGMA聚类定义母系，NJ树与12个Gallus线粒体基因组比较；以变异等位基因频率（variant allele frequency，VAF）判定异质性；以MT深度/核深度估计相对mtDNA拷贝数；测定体重、脏器及132项肉品质指标，采用方差分析或Kruskal–Wallis检验，Benjamini–Hochberg校正。
研究结果：3.1 测序质量与比对：18个样本Q30为95.96%–96.75%，GC为43.34%–45.09%，平均比对率95.03%，平均深度2.99×–4.69×，≥1×和≥4×覆盖分别为88.51%和32.23%；HX1比对率62.94%、≥1×覆盖71.95%明显偏低，解释需谨慎但不影响群体结论。通过重测序质控与比对得出数据总体高质量且组间可比。3.2 线粒体变异谱：以GRCg6a MT 16,775 bp为参考，90个MT变异中85个多态，转换/颠换比约11.9；58个编码（42同义、16非同义）、20个D-loop/非编码、10个rRNA、2个tRNA；无InDel或SV；非同义/同义约0.38，提示编码基因受纯化选择；多数位点由多只鸡共享。16个非同义替换分布于ND2（4个）、CO1（3个）、CO2（2个）、ND4L（1个）、ND5（5个）和CYTB（1个），集中于呼吸链复合体I、IV和III，ATP合酶及ND3/ND4未检出。通过变异提取和注释得出变异以点突变为主、结构保守，复合体I/III/IV存在候选非同义变异。3.3 母系单倍型与系统发育：基于90个位点将18只鸡分为A（LX1、LX3、LX4、LX6、MX4、HX6）、B（HX1、HX2、HX4、HX5、MX1、MX2、MX3、MX5、MX6、LX2）和C（LX5、HX3）；A在低海拔组较多（4/6），B在中高海拔占优，呈弱母系分化但非完全由海拔决定。NJ树中群体靠近Gallus gallus gallus（NC_007236）和中国地方品种；与NC_007236差异最小（21.3个位点），与PQ664589、PQ664590、PQ412228约26个，与NC_001323为37.2个；B含三海拔个体，C仅2只。通过聚类和NJ树得出母系部分结构化、母系起源与中国地方鸡接近。3.4 异质性：检测15个异质位点（0.05 < VAF < 0.95）；16S rRNA nt3946（G/C/A）在全部18只鸡中异质，非参考等位频率0.55–0.87且读段支持强；nt258、nt278、nt253、nt240、nt222在多只样本呈中间频率；ND5 nt13450仅HX2（VAF = 0.24）和HX5（VAF = 0.30）异质；低深度样本如MX2（DP = 15）需谨慎。异质位点富集于16S rRNA和D-loop，少数在编码区。通过VAF和读段水平证据得出异质集中非编码及翻译相关区域。3.5 mtDNA拷贝数：个体变异大（MX2约1989，MX5约11484），组均LX 7656、MX 6405、HX 8488，Kruskal–Wallis p = 0.52，无海拔趋势。通过MT深度/核深度估计得出相对拷贝数组间无显著差异。3.6 表型与肉品质：心、肺重量和IMP随海拔增加（心LX/MX/HX为8.54/8.23/11.02 g；IMP为775.2/778.5/828.3 mg/100 g），按母系无一致差异（心A/B/C为8.94/9.41/9.48 g），提示随海拔而非MT世系。蛋白最高在MX（94.28 g/100 g），最低在LX（91.28 g/100 g，p = 0.041）；脂肪最高在LX（2.10 g/100 g），最低在MX（1.73 g/100 g）；IMP最高在HX（828.33 mg/100 g，p = 0.025）；UMP最低在MX（8.71 mg/100 g，p = 0.011）；总游离氨基酸、脯氨酸、赖氨酸和精氨酸在中高海拔较高（p < 0.05），苯丙氨酸在MX较低（3.53 g/100 g，p = 0.008）；油酸在LX最高（0.39 g/100 g，p = 0.020），DHA在中高海拔较高（p = 0.024）；正己烷在中高海拔较高（p = 0.021）。BH校正后p ≤ 0.041仍显著（q ≈ 0.036），脂肪（p = 0.093）和水分（p = 0.899）不显著。通过表型测定和统计得出肉品质随海拔协调变化，而不与母系对应。
讨论总结：90个位点大多为转换，85个多态，编码变异以同义为主，非同义/同义约0.38，提示整体受纯化选择；16个非同义集中于复合体I（ND2、ND4L、ND5）、复合体IV（CO1、CO2）和复合体III（CYTB），这些复合体与高海拔及高能量需求适应有关；D-loop和rRNA多态丰富，符合鸡单倍群定义区域。多数变异可能中性或近中性，但少量非同义可能影响亚基折叠、质子转运及核-线粒体互作，候选如ND5 nt14697 p.His543Asn和CYTB nt15535 p.Asp215Asn需结构建模和功能验证；复合体IV为终端氧化酶，复合体I为活性氧（reactive oxygen species，ROS）主要产生位点，肌肉能量需求高，低氧下效应可能更明显，但本数据不能证明。三系A低海拔多、B中高海拔多，弱结构；NJ靠近G. gallus gallus和中国地方品种，可能含红原鸡祖先；同组多系共存、跨海拔共享反对强海拔特异选择，提示鸟类移动和历史母系多样性；村养群体交换、重叠世代和小有效群体造成马赛克，线粒体单独作为起源或适应代理受限，需结合核基因组。异质性15个位点，16S rRNA nt3946全群异质且支持强，可能反映复制-分离动态、松弛纯化选择、突变-选择平衡，也可能受深度和比对影响，需Sanger测序或靶向深度确认；ND5 nt13450仅少数个体；异质性受遗传瓶颈、Balbiani体和线粒体自噬影响，随年龄和代谢状态变化；本群体异质多低频且集中于16S rRNA和D-loop；检测深度依赖，15个为下限。拷贝数个体间>5倍变异，组间无差异无趋势；mtDNA拷贝数受组织能量需求和氧化应激调控，骨骼肌高拷贝数个体变异大；高海拔适应不是简单增加mtDNA量，可能调控OXPHOS效率和ROS控制等质量因素；相对值需qPCR或数字PCR绝对定量；纤维类型、年龄和代谢状态也影响，不应过度解释。心、肺重量随海拔增加，符合藏鸡适应；肉品质随海拔而非母系：中海拔高蛋白低脂，高海拔IMP最高（主要鲜味化合物），低海拔油酸最高但DHA较低，中高海拔DHA较高；与藏鸡代谢调整和慢生特性一致；转录组研究提示能量代谢和缺氧信号广泛变化、线粒体和代谢通路协调重塑，说明海拔差异可能由转录重塑而非MT序列变异单独驱动；无母系一致关系是重要阴性结果，线粒体对肉品质贡献可能相对核基因、年龄、日粮和饲养环境小；家禽经济性状多基因，单个MT变异影响弱；选择应依赖表型和核基因组，MT数据用于追溯和母系多样性监测。局限包括样本量18、深度约4×有限，异质性和拷贝数为初步，需靶向深度验证；MT SV/CNV检测不敏感，需高深度或长读长；MT变异-表型关联为探索性，候选非同义功能需实验；每海拔单一饲养点使海拔与地点环境不能完全分离，应解释为海拔关联而非直接适应证据；未做qPCR/数字PCR绝对定量及Sanger/靶向深度正交验证；需扩大样本、海拔范围和饲养环境数量，整合核基因组和功能证据。综合意义：木里藏香鸡MT基因组有大量现存变异，多同义或非编码，母系仅部分按海拔结构，异质常见，拷贝数稳定；肉品质随海拔显著但与母系无关；高海拔适应主要由核基因和环境驱动，MT贡献现存多样性和母系标记；这些标记和候选变异支持品种认证、追溯和保护规划，生产优质藏香鸡产品应考虑饲养海拔。
研究结论翻译：木里县藏香鸡胸肌线粒体基因组含丰富序列变异（90个位点，其中16个非同义）和异质性（15个位点，16S rRNA nt3946在每只鸡中均异质）。鸡形成三个母系（A、B、C），仅轻度海拔分化，母系起源与中国地方鸡相似。相对mtDNA拷贝数组间无差异。肉品质包括2800 m高蛋白低脂和3500 m最高IMP，主要与海拔而非母系相关。综上，这些发现初步表征该品种线粒体变异：观察到的表型差异主要与饲养海拔而非母系线粒体世系相关，MT基因组似乎贡献现存遗传多样性及可用于品种认证和保护的母系标记。这些解释为探索性，需更大样本和独立站点确认。
