纳米比亚猎豹表型转变相关的表观遗传变化

《Molecular Ecology》：Epigenetic Changes Associated With a Phenotypic Transition in Namibian Cheetahs

　　

在遗传多样性较低的种群中，适应环境变化的能力有限。虽然表观遗传机制可以产生表型可塑性，但其在介导野生种群适应性生活史转变中的作用在很大程度上尚不清楚。本研究表明，在成年纳米比亚猎豹（Acinonyx jubatus jubatus）——一种遗传多样性极低的物种中，从非领地流浪个体（floater）转变为占有领地的雄性个体，伴随着广泛的DNA甲基化变化，其中领地雄性以高甲基化为主。研究人员利用同一批个体在成为领地雄性之前和之后采集的配对血液样本，鉴定出差异甲基化基因，这些基因在解剖结构及头部发育、大脑发育与神经发育、cAMP生物合成及离子通道转运等相关通路中发挥作用，与观察到的形态、空间运动及攻击行为的表型转变一致。这些发现表明，DNA甲基化在分子水平上介导了这一生活史转变。更广泛而言，它们提供了罕见的证据，表明表观遗传调控与遗传多样性较低物种的表型多样性相关，并可能在种群水平上支持其适应潜力。

该论文发表于《Molecular Ecology》，聚焦于纳米比亚猎豹（Acinonyx jubatus jubatus）从非领地流浪雄性（floater）向领地雄性转变过程中伴随的表观遗传变化。当前全球生物多样性面临快速下降，遗传多样性降低使物种适应能力受限，但部分低多样性物种仍表现出一定韧性，提示可能存在遗传变异之外的可塑性机制。表观遗传修饰尤其是DNA甲基化可以在不改变DNA序列的情况下调控基因表达，但关于其介导内在生活史转变的研究仍十分有限。纳米比亚猎豹遗传多样性极低，其雄性个体在成年初期均为floater，部分个体在与领地雄性争斗后获得领地，进而表现出去领地雄性显著的形态和行为改变，包括体重增加、头骨增大、攻击性和大胆性增强，是研究野生哺乳动物表观遗传修饰的理想自然体系。研究人员基于长期追踪的配对样本，旨在解决三个问题：鉴定与空间策略转变相关的差异甲基化区域（DMR）；确定这些DMR相关基因是否与观察到的形态和行为变化相符，并校正年龄混杂因素；评估这些表观遗传信号是否可作为区分两种空间策略的预测性生物标志物。

在技术方法方面，研究人员从纳米比亚长期猎豹研究项目中获取样本队列：10只雄性个体先以floater身份被采集血样，待其转变为领地雄性后再次采集配对样本；另有8只始终为floater的个体用于体重比较。关键方法包括：从外周血单个核细胞（PBMC）提取DNA，进行甲基化DNA结合域测序（MBD-Seq）；将测序读段比对到猎豹参考基因组，以300 bp窗口为单位计数；采用嵌套“留一个体”交叉验证框架，并联合非参数检验、参数分类模型和参数连续模型三种统计方法，对各折间具有一致方向的DMR进行筛选；对差异甲基化基因（DMG）进行基因本体（GO）功能注释；利用ElasticNet和随机森林模型评估DMR区分floater与领地雄性的预测能力。

在结果部分，研究涵盖以下发现：

3.1 保持floater的雄性与转变为领地雄性的雄性之间的体重差异
配对体重数据显示，10只完成转变的雄性平均体重由floater阶段的46.81 kg（SD=5.61）增至领地雄性阶段的54.35 kg（SD=5.18），个体内平均增加7.54 kg（SD=3.49，配对t检验p<0.001）；而8只始终保持floater的雄性仅从46.66 kg（SD=4.53）增至47.59 kg（SD=3.87），平均增加0.93 kg（SD=1.39，p=0.102）。两组间体重增加量差异为6.62 kg（Welch t检验p<0.001）。在floater样本年龄分层分析中，年龄6类（>2.0至≤3.5岁）与年龄7类（>3.5至≤7.0岁）个体的体重差异不显著，表明年龄并非体重变化的主要原因。

3.2 测序结果
MBD-Seq共获得高质量配对末端读段，floater样本平均为29,001,044条clean reads，领地雄性为29,541,538条；平均比对率分别为94.3%和90.8%。除一个领地雄性样本比对率为77.0%外，其余样本均高于88%，整体数据质量良好，满足下游差异分析要求。

3.3 甲基化分析
在三种统计方法中，与floater相比，领地雄性的多数DMR表现为高甲基化。参数分类分析（PCa）鉴定到1675个高甲基化DMR，参数连续分析（PCo）为624个，非参数分析（nP）为1544个，均显著高于随机窗口中的高甲基化比例（95%引导置信区间53%–58%）。三种方法共同检测到98个一致性DMR，至少两种方法共同检测到的DMR有554个，这些方法无关的DMR对应481个差异甲基化基因，其中97个基因被三种方法同时检出。年龄分析鉴定了745个年龄相关DMR，与空间策略转变相关的DMR在DMR水平上完全没有重叠；但在基因水平上，年龄相关DMG与策略相关DMG存在少量但高于随机预期的重叠，表明年龄信号与策略转变信号在基因注释层面存在一定汇聚。

3.4 生物学过程
对DMG进行GO注释后，按证据得分排序最高的80个GO生物过程术语聚类为七组：神经发育、解剖结构发育、头部与脑发育、细胞运动、学习与认知、离子通道转运和cAMP生物合成。这些功能类别与领地雄性在形态（体重、头骨）、行为（攻击性、大胆性）和空间运动（家域缩小）方面的变化一致。相比之下，年龄相关GO术语分布更为松散，主要涉及淋巴系统发育、甾醇代谢和细胞突起等一般性过程，但其中也包括部分与策略信号共享的术语，如解剖结构形态发生和脑发育。

3.5 与空间策略转变相关的启动子重叠DMR
启动子区域甲基化对基因调控具有重要作用。研究共在61个策略转变相关基因中检测到启动子重叠DMR，其中55个在领地雄性中表现为高甲基化，仅6个为低甲基化，且这些基因与年龄相关启动子集无重叠。方法无关的18个基因涉及发育过程、肌肉收缩、细胞信号传导和代谢等，例如MYO18B、DTNA和ADAM12参与肌节装配、膜稳定和肌母细胞融合；PLXNA1和PTPRS参与神经元导向和轴突生长；TRPM1和MFSD10参与离子通道和能量代谢；SMC3参与染色体黏连。这些启动子甲基化改变可能为领地转变过程中的体细胞和神经肌肉变化提供表观遗传基础。

3.6 一致性DMR作为区分floater与领地雄性的生物标志物
在交叉验证框架下，基于非参数方法得到的DMR分类性能最高，平均准确率为78.2%（95%引导置信区间75%–85%）；参数分类和参数连续方法的准确率分别为76.2%（70%–85%）和74.0%（65%–80%）。非参数方法在折间稳定性也最佳，有19个DMR在所有10折中被ElasticNet一致选择。综合三种方法，最终筛选出21个候选生物标志物DMR，其中7个重叠于已知基因（GLI2、ZFYVE28、PARVB、MUSK、ZNF423、SHISA6、EFCC1），2个位于未知基因区，12个位于基因间区。这些标志物需在独立队列中进一步验证。

在讨论部分，研究人员指出，尽管年龄因素与策略转变存在部分基因水平重叠，但策略转变信号在功能上高度凝聚，而年龄信号较为弥散，说明表型转变可能依赖于多个微弱甲基化信号汇聚为一致性调控模式，而非单一DMR的显著变化。交叉验证设计能够提供样本外预测评估，但最终生物标志物列表仍需独立验证。研究还强调，甲基化信号来自PBMC，不一定代表脑、肌肉或颅骨组织的直接功能改变，而可能反映与空间策略转变相关的系统性生理状态。淋巴细胞比例在两类雄性中相似且较低，排除了急性感染或主要细胞亚群比例差异的干扰。此外，配对设计控制了固定遗传差异，但由于缺乏转变过程中的中间时间点样本，尚无法判断甲基化变化是转变的先因、伴随还是结果。

研究结论如下：本研究证明了在自由放养的雄性猎豹中，DNA甲基化与一次生活史转变及其伴随的表型变化之间存在一致性关联。所鉴定的差异甲基化区域重叠的基因在解剖发育、神经发育和能量转换等生物学过程中发挥功能，这与floater和领地雄性之间的形态与行为差异相呼应。通过利用野生哺乳动物的配对样本设计，研究人员的结果表明，个体内部从一个表型到另一个表型的转变与独特的DNA甲基化模式相关，无需个体间的固定遗传差异。这种调控灵活性可能有助于种群水平所表达的表型响应范围，这在遗传多样性有限的物种中可能尤为相关。未来的研究应探索类似机制是否在其他受威胁物种中运行，以及是否存在跨类群的通用表观遗传特征可作为适应性能力的生物标志物。预测性表观遗传标记的鉴定也使保护生物学中的实际应用成为可能，例如对空间策略分布或种群状态进行非侵入性监测。最终，将表观遗传视角整合到保护研究中，可以在生物多样性加速丧失的时代完善灭绝风险预测并指导管理策略。