全基因组重测序揭示 Rhabdophis nuchalis 的遗传多样性、群体结构与适应性进化

《Ecology and Evolution》：Whole-Genome Resequencing Reveals the Genetic Diversity and Adaptive Evolution of Rhabdophis nuchalis

蛇类因其形态特化、广泛分布和强适应性，对研究适应性进化具有重要价值。然而，针对 Rhabdophis nuchalis 的基因组研究有限。研究人员对来自4个地区的19个个体进行了重测序，以评估遗传多样性和适应性。研究人员获得了39,226,324个高质量单核苷酸多态性（single nucleotide polymorphisms, SNPs）和15,838,348个高质量插入缺失（insertions and deletions, INDELs）。遗传结构分析显示，四川和云南群体形成一个聚类，湖北与陕西群体形成另一个聚类。在四个采样群体（四川、陕西、湖北和云南）中，云南表现出较高的核苷酸多样性和观测杂合度，表明其具有较高的遗传多样性。陕西表现出较低连锁不平衡，提示其有效种群大小较大，与历史估计一致。基因本体（Gene Ontology, GO）富集揭示了与繁殖和钙依赖性细胞黏附相关的分化基因。京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）仅鉴定出细胞黏附分子通路。基础局部比对搜索工具（basic local alignment search tool, BLAST）将这些序列与原钙黏蛋白β和γ亚簇比对。总体而言，四个群体之间存在清晰的遗传聚类和分化，分化相关基因涉及繁殖和细胞黏附，为遗传多样性保护提供了见解。

论文解读：全基因组重测序揭示 Rhabdophis nuchalis 的遗传多样性、群体结构与适应性进化

蛇类具有高度形态特化，广泛分布于除极端极地环境外的大陆，并在生态系统中作为高效捕食者维持生态平衡。自其穴居祖先分化以来，蛇类发生了多种适应性进化改变，部分类群还形成了独特形态与生理特征。因此，蛇类被视为研究遗传多样性和适应性选择的重要对象。Rhabdophis nuchalis 属爬行纲有鳞目蛇亚目游蛇科颈槽蛇属，因颈部中线两行背鳞形成浅槽而得名。该物种命名和分类历经多次修订，Boulenger 曾以 Tropidonotus nuchalis 描述，后归入 Natrix sensu lato，Malnate 又将其分立并建立 R. nuchalis。现有研究多集中于形态、系统发育和生物地理，而基因组研究很少，迄今仅有单样本基因组组装，缺乏多样本深入研究。全基因组重测序（whole-genome resequencing）是下一代测序（next-generation sequencing, NGS）的重要应用，可通过比对参考基因组识别遗传变异，研究物种生物多样性与进化。此前仅少数研究将重测序用于蛇类。尽管已知 R. nuchalis 在中国分布较广，但不同地理群体间遗传结构和分化仍不清楚。因此，研究人员开展本研究，旨在利用全基因组重测序探讨 R. nuchalis 的遗传多样性和适应性进化特征。研究人员对来自中国四川（SC）、云南（YN）、陕西（SX）和湖北（HB）四个地区的19个个体进行全基因组重测序，分析遗传多样性和群体结构，并结合全基因组选择信号分析，初步鉴定不同地区群体间显著遗传分化的关键通路及受选择同源基因片段。该研究为制定 R. nuchalis 保护策略、维持区域生物多样性提供理论基础和数据支持，也为非模式生物基因组研究提供案例。

关键技术方法方面，研究人员采集湖北神农架3只、陕西安康8只、四川攀枝花5只和云南大理白族自治州3只，共19个 R. nuchalis 个体，海拔范围为1200–2750 m，采集肝组织并送至上海Personal Biotechnology Co. Ltd.进行重测序。文库插入片段为400 bp，在BGISEQ DNBSEQ-T7平台进行双端2×150 bp测序，目标深度10×。质控后，用Burrows-Wheeler Aligner（BWA）比对参考基因组GCA_039707465.1，用GATK进行变异检测并筛选高质量SNP和INDEL。随后采用PLINK、FastTree、Admixture、VCFtools、PopLDdecay、SMC++、TBtools和R v4.3.3进行群体结构、遗传多样性、历史动态、选择信号和功能富集分析，并以Tajima's D进行中性验证。

结果部分按原文小标题总结如下。

3.1 Whole-Genome Resequencing and SNP Calling：研究人员对19个 R. nuchalis 样本进行全基因组重测序，获得277.2 Gb原始数据，质控后得到268.1 Gb高质量数据，平均Q20为98.75%，平均Q30为96.51%。高质量数据成功比对到参考基因组，平均覆盖深度11.5×，平均比对率99.71%。GATK变异检测共鉴定84,612,445个SNP和16,002,200个INDEL，质控后获得39,226,324个高质量SNP和15,838,348个高质量INDEL。

3.2 Population Genetic Structure Analysis：主成分分析（principal component analysis, PCA）显示，第一主成分（PC1）解释28.5%的总方差，有效区分SC和YN群体与HB和SX群体；第二主成分（PC2）解释15.19%的总方差，进一步揭示SC与YN的分化，而HB与SX紧密聚类。Admixture分析显示，K=2为最佳支持，YN和SC所有个体及部分SX个体共享同一遗传成分，其余SX个体与HB个体形成另一聚类；K=3和K=4时SX群体出现亚簇。基于SNP的系统发育树显示HB、YN和SC各自成簇，YN与SC聚为一大谱系，SX未形成单系，提示其内部可能存在分化，与Admixture结果一致。

3.3 Population Genetic Diversity Analysis：核苷酸多样性（nucleotide diversity, θπ）分析显示YN最高，SC最低。所有群体观测杂合度（observed heterozygosity, H_o）均低于期望杂合度（expected heterozygosity, H_e），YN的H_o最高，HB最低，表明存在杂合子缺陷。连锁不平衡（linkage disequilibrium, LD）衰减分析显示SX衰减最快，在100 Kb时其LD水平（r²）低于YN、HB和SC；YN、HB和SC的LD水平相近且高于SX。固定指数（fixation index, F_ST）显示HB与SX之间分化相对较低，提示二者间可能存在较强基因流。

3.4 Population Historical Dynamics Analysis：研究人员对四个群体的历史有效种群大小（effective population size, N_e）进行溯祖分析。总体上四个群体历史N_e均呈下降趋势，早期略有增加。HB、SC和YN在初期轻微增加后急剧下降，过程中虽有波动但总体持续下降。SX则表现为先波动下降、随后快速恢复、再下降并逐渐趋稳。

3.5 Selective Sweep Analysis：研究人员以F_ST前5%且θπ后5%的窗口为候选区域，进行选择清除分析及GO和KEGG富集。GO富集中，HB与SX比较无显著富集；HB与SC比较主要富集于多生物生殖过程和发育过程；SX与YN比较仅富集于卵子发生；HB与YN比较显著富集于钙依赖性细胞–细胞黏附；SC与SX比较显著富集于生殖细胞发育、多细胞生物生殖过程和配子发生；SC与YN比较显著富集于钙依赖性细胞–细胞黏附。KEGG分析中，HB与YN、HB与SC、SC与YN比较均显著富集于细胞黏附分子通路。结合Tajima's D验证，HB相关比较中HB本身无显著选择清除信号，提示分化信号可能主要来自其他三个群体；非HB比较中SC、SX和YN均表现出显著选择信号，说明SC、SX和YN在本研究观察到的分化过程中可能发挥更重要作用。

讨论部分总结：PC1区分SC-YN与HB-SX两个主要地理遗传群，支持地理隔离解释。Zhu et al. 认为四川盆地构成地理屏障，阻碍物种迁移，可能限制SC-YN与HB-SX之间的基因流。Admixture显示并非严格按地理隔离聚类，提示群体间可能存在基因流或历史扩散。Zhu et al. 发现 R. nuchalis 起源于横断山并向秦岭-大巴山扩张，四川、陕西和湖北群体约在1.54 Mya分化为两支，陕西个体被分配到不同簇，这与本研究的Admixture结果相似。有学者提出，滇西高山峡谷、南岭以及湘川鄂边境受第四纪冰期影响较小，可作为重要避难所；冰期避难所可促进种下分化。SC-YN与HB-SX的分化可能与末次冰期以来不同避难所隔离有关，SX的混合祖先可能反映后续扩散或二次接触中的基因流。YN具有最高核苷酸多样性，可能与其所在云南丰富的物种多样性和生态系统类型有关。所有群体H_oe，可能反映Wahlund效应或一定程度的近交。SX的LD衰减最快，提示其可能具有更高重组率或更大N_e，与历史N_e估计一致。历史N_e分析显示HB、SC、SX在0.1–0.01 Ma出现波动后快速下降，可能与大理冰期（Dali Glaciation，0.07–0.01 Mya）有关；YN急剧下降可能因缺乏合适避难所或避难所内强烈种内竞争；SX在约0.001 Mya快速增加，可能与太白冰期（Taibai Glaciation，0.01862–0.00073 Mya）有关。GO富集中，繁殖相关功能在多个群体对中反复显著富集，提示繁殖相关基因分化可能是种内群体分化的核心驱动之一。YN与HB、YN与SC比较中，钙依赖性细胞–细胞黏附显著富集，cadherin超家族成员可能参与其中。N-钙黏蛋白（N-cadherin）可介导细胞–细胞相互作用，并可与磷脂酰肌醇3-激酶（phosphatidylinositol 3-kinase, PI3K）和Rac家族小GTP酶（Rac family small GTPases, Rac）等信号分子相互作用，影响细胞定向迁移。KEGG仅鉴定出细胞黏附分子通路，BLAST显示相关序列与protocadherin基因家族β和γ亚簇高度同源。原钙黏蛋白主要参与脊椎动物神经元连接和细胞间通讯，是神经系统发育的关键基因家族，其分化提示不同地理群体可能在神经系统相关基因上发生适应性分化；γ亚簇具有复杂可变剪接特征，相关分化可能影响剪接模式或异构体表达。保护建议方面，研究人员建议优先保护遗传多样性最高且位于生物多样性热点地区的YN群体；SX群体个体呈现混合祖先，可能为不同遗传群之间的混合区，其遗传组成应持续监测；鉴于主要地理群体间存在遗传分化，建议将其作为独立管理单元，优先监测YN群体的遗传多样性。

局限性方面，研究人员指出样本量相对较小且分布不均，可能影响等位基因频率估计和遗传多样性评估；由于缺乏 Rhabdophis 属物种特异性突变率，研究人员采用 Hydrophis curtus 的突变率4.71×10^?9 per site per year推断种群历史，可能给历史动态准确估计带来不确定性，因此对上述结果应谨慎解释。未来研究应纳入更大群体样本和更精确的物种特异性参数，以进一步验证并深化对 Rhabdophis nuchalis 适应性进化机制的理解。

结论翻译：本研究阐明了来自不同地理区域的 R. nuchalis 群体的遗传结构和适应性分化模式。云南群体表现出高遗传多样性，陕西群体维持较大的历史有效种群大小，为该物种遗传资源保护提供了重要参考价值。繁殖过程和细胞黏附通路可能通过调控繁殖和细胞通讯功能在物种遗传分化中发挥关键作用。这些发现加深了对 R. nuchalis 适应性选择的理解，为制定以遗传资源保护为中心的群体管理策略提供了科学依据，也为爬行动物微进化提供了新见解。