海拔梯度下藏香鸡线粒体基因组变异及其与肉品质性状的关联分析

《Animals》：Acute Hyperkalaemia During General Anaesthesia in Two Cows

　　

本研究旨在表征藏香鸡沿海拔梯度的线粒体（MT）基因组特征，并探讨其与肉品质的关联。研究人员选取了来自三个海拔组的18只150日龄雄性鸡（LX组，1500 m；MX组，2800 m；HX组，3500 m；每组n = 6；每组从30只饲养个体中随机抽取6只），进行全基因组重测序（参考基因组：Gallus gallus GCA_000002315.5）。研究共鉴定出90个MT变异位点（其中85个为多态性位点；转换/颠换比83:7），涵盖58个编码区位点（42个同义突变、16个非同义突变）、20个D环/非编码区位点、10个rRNA位点和2个tRNA位点；未检测到MT插入/缺失（insertions/deletions，InDels）或结构变异（structural variants，SVs）。18只鸡聚为三个母系谱系（A、B、C），其海拔结构仅部分分明（A谱系在低海拔更常见；B谱系在中、高海拔占优）。基于12个Gallus属线粒体基因组构建的邻接树（neighbour-joining tree）将该群体定位于接近Gallus gallus gallus和中国地方鸡种的位置。研究共检测到15个异质性（heteroplasmy）位点，其中16S rRNA的nt3946位点（G/C/A）在所有个体中均呈异质性。相对mtDNA拷贝数在各组间无显著差异（p = 0.52）。肉品质随海拔而变化（2800 m时蛋白质含量最高、脂肪含量最低；3500 m时肌苷单磷酸（inosine monophosphate，IMP）含量最高），但不随母系谱系变化，提示所观察到的表型差异主要与饲养海拔相关，而非与母系线粒体谱系相关；MT基因组对这种差异贡献的是已有的遗传变异（standing genetic variation），而非其直接的遗传学基础。

**研究背景与科学问题**

藏香鸡是长期自然选择与农户驯化在青藏高原形成的地方鸡种，是少数能在约2200–4100 m海拔完成生命周期的鸡种之一，其肌肉富含氨基酸和肌苷单磷酸（inosine monophosphate，IMP），具有浓郁鲜味和高营养价值。木里藏族自治县位于横断山脉东缘，海拔跨度超过4400 m，长期地理隔离可能造就与其他藏香鸡不同的特征。线粒体是真核细胞有氧呼吸的核心细胞器，其基因组（mtDNA）为母系遗传、进化速率快、拷贝数高且频繁携带异质性（heteroplasmy），是母系结构与适应性进化研究的重要标志。由于呼吸链复合体由核基因和线粒体基因共同编码，线粒体编码亚基的变异可能参与复合体对低氧环境的适应，因此被视为缺氧适应的候选位点，既往研究已将藏香鸡MT-CO3单倍型及MT-COI、ATP6变异与高海拔适应相关联。然而，现有线粒体研究多聚焦单一基因或D环，缺乏沿海拔梯度整合序列变异、母系结构、异质性和拷贝数的系统性研究，木里群体的线粒体特征及其与肉品质的关联尚未被表征。同时，高原地方鸡种正面临与商业品种杂交的遗传侵蚀风险，记录母系谱系可为遗传侵蚀监测提供基线；肌肉线粒体特征可能间接影响肉品质性状，这对藏香鸡高端市场具有生产意义。

**研究内容、结论与意义**

研究人员以木里县三个海拔（1500 m、2800 m、3500 m）饲养的18只150日龄雄性藏香鸡为研究对象，采集胸肌进行全基因组重测序，围绕五个目标开展研究：表征MT变异谱、界定并在系统发育上定位母系谱系、描述异质性、估算相对mtDNA拷贝数，并将上述特征与表型和肉品质变异相关联以检验其是否随海拔变化。研究得出：MT基因组以转换变异为主、结构高度保守；母系谱系仅部分沿海拔分化；相对拷贝数在各海拔组间无显著差异；而蛋白质、脂肪、IMP等肉品质性状主要随海拔而非母系谱系变化。结论认为，该品种的高海拔适应主要由核基因和环境驱动，线粒体基因组主要贡献已有的遗传多样性，同时可作为品种鉴定、溯源与保护规划的母系标记。这对高原鸡种遗传资源保护和高端藏香鸡产品生产均具有实践价值。

**关键技术方法**

样本队列为四川省木里藏族自治县示范农场150日龄雄性藏香鸡18只，分低海拔（LX，1500 m）、中海拔（MX，2800 m）、高海拔（HX，3500 m）三组，每组6只（每组从30只中随机抽取）。关键技术方法包括：Illumina NovaSeq 6000平台PE150双端全基因组重测序（约4×深度），BWA比对至参考基因组Gallus gallus GCA_000002315.5；SAMtools mpileup鉴定SNP与小InDel，BreakDancer检测结构变异，CNVnator检测拷贝数变异，ANNOVAR功能注释；基于90个MT变异位点构建UPGMA矩阵及邻接树（neighbour-joining tree）并与12个Gallus属线粒体基因组比较定位；依据等位基因深度（allele depth，AD）与变异等位基因频率（variant allele frequency，VAF）评估异质性；以MT与核基因组平均深度之比估算相对拷贝数；对132项肉品质及表型性状进行方差分析或Kruskal–Wallis检验，并以Benjamini–Hochberg校正控制错误发现率。

**研究结果**

3.1 测序质量与比对：18个样本Q30为95.96%–96.75%，GC含量43.34%–45.09%，平均比对率95.03%，平均测序深度2.99×–4.69×；仅高海拔样本HX1比对率（62.94%）和≥1×覆盖度（71.95%）较低，其变异解读需谨慎，但不影响群体水平结论。

3.2 线粒体变异谱：基于GRCg6a MT序列（16,775 bp）鉴定出90个变异位点，其中85个为群体内多态；转换/颠换比为83:7（约11.9）；58个位点位于编码区（42个同义、16个非同义）、20个位于D环/非编码区、10个位于rRNA、2个位于tRNA；非同义/同义比约0.38，提示纯化选择作用；未检测到MT插入/缺失或结构变异，表明该群体MT基因组以点突变为变异主体且结构保守。16个非同义替换分布于复合体I（ND₂、ND₄L、ND₅）、复合体III（CYTB）和复合体IV（CO₁、CO₂），ND₅位点最多（5个），提示编码区变异并非局限于单一呼吸链复合体。

3.3 母系单倍型与系统发育定位：基于90个变异位点，18只鸡分为三个母系谱系（A、B、C），仅部分按海拔分化——A谱系在低海拔组更常见（6只中4只），B谱系在中、高海拔占优；邻接树（含12个Gallus属线粒体基因组）将该群体定位于接近Gallus gallus gallus和中国地方鸡种的位置，与NC_007236的平均差异最小（21.3个位点）。

3.4 异质性：检测到15个异质性位点（0.05
3.5 mtDNA拷贝数：相对拷贝数个体间差异超过5倍（MX2约1989至MX5约11484），但三组间无显著差异（Kruskal–Wallis p = 0.52），亦无一致海拔趋势，提示高海拔适应并非通过简单增加mtDNA拷贝数实现。

3.6 表型与肉品质：心脏和肺重量及IMP含量随海拔升高；肌肉蛋白质在2800 m最高（94.28 g/100 g）、脂肪在2800 m最低（1.73 g/100 g），呈"高蛋白、低脂肪"特征；3500 m时IMP最高（828.33 mg/100 g，p = 0.025）；苯丙氨酸在2800 m较低（3.53 g/100 g，p = 0.008）；油酸在低海拔最高（0.39 g/100 g，p = 0.020），DHA在中、高海拔较高（p = 0.024）；己烷在中、高海拔显著较高（p = 0.021）。经Benjamini–Hochberg校正后，蛋白质（p = 0.041）等关联仍显著（q≈0.036），而脂肪（p = 0.093）和水分（p = 0.899）不显著。按母系谱系分组未解析出一致模式，表明肉品质差异追踪海拔而非母系祖先。

**讨论部分总结**

讨论指出，90个MT变异位点以转换为主，大多数变异为同义或非编码变异，多数呈中性或近中性；16个非同义替换集中于复合体I、III、IV，与既往高海拔适应研究相符，但具体氨基酸替换与既往报告不同，其后果需结构建模与功能验证。母系结构仅部分与海拔相关，反映母系群体统计学（历史母系多样性与鸟类流动）而非强海拔特异性选择，线粒体谱系保留已有多样性。nt3946位点在所有个体中异质并不寻常，可能反映特异的复制–分离动态、松弛的纯化选择或突变–选择平衡，亦可能部分源于测序深度与比对特性，需正交方法验证；15个异质性位点应视为下限。拷贝数无海拔趋势提示适应可能通过调控线粒体功能（氧化磷酸化（oxidative phosphorylation，OXPHOS）效率与活性氧（reactive oxygen species，ROS）控制）实现，即"调整质量而非数量"。肉品质随海拔而非母系谱系变化是较有信息量的阴性结果，表明母系基因组本身不决定所测肌肉表型，对肉品质的选择应依赖表型和核基因组信息，线粒体数据最适用于溯源与母系多样性监测。研究局限性包括：样本量小（n = 18）、测序深度低（约4×）、每个海拔仅单一饲养点（海拔与地点效应不能完全分离）、关联分析属探索性、异质性与拷贝数结果需正交验证。综合而言，藏香鸡线粒体基因组携带大量已有变异，母系谱系仅部分沿海拔分化，异质性常见而拷贝数稳定；肉品质随海拔变化但与母系谱系无关；高海拔适应主要由核基因和环境驱动，线粒体基因组提供的母系标记可支持品种鉴定、溯源与保护规划，高端藏香鸡产品生产应考虑饲养海拔。

**研究结论部分翻译**

藏香鸡胸肌线粒体基因组具有丰富的序列变异（90个位点，其中16个为非同义突变）和异质性（15个位点，其中16S rRNA nt3946位点在每个个体中均呈异质性）。这些鸡只形成三个母系谱系（A、B和C），仅表现出轻微的海拔分化，其母系起源与中国地方鸡种相似。相对mtDNA拷贝数在各海拔组间无差异。肉品质，包括在2800 m处的高蛋白、低脂肪以及在3500 m处最高的IMP，主要与海拔而非母系谱系相关。综合来看，这些发现初步表征了这一品种的线粒体变异：所观察到的表型差异主要与饲养海拔相关，而非与母系线粒体谱系相关，且线粒体基因组似乎贡献的是已有的遗传多样性以及可能可用于品种鉴定和保护的母系标记。这些解释是探索性的，需要在更大样本量和独立地点的研究中得到证实。