漂变与选择在近期殖民海域中的博弈：波罗的海港湾鼠海豚低盐适应性的基因组基础

《Molecular Ecology》：Drift and Selection in a Recently Colonized Sea: The Genomic Basis of Adaptation to Low-Salinity Waters in Baltic Harbour Porpoises

　　

专业术语翻译 理解种群动态（demography）与自然选择如何在新的环境殖民过程中相互作用，是进化生物学中的一个核心问题。波罗的海区域结合了世界上最强的盐度梯度之一以及特征明确的、不足约12,000年的冰后期历史。港湾鼠海豚（Phocoena phocoena）沿该梯度分布为三个种群——北海（North Sea）、贝尔特海（Belt Sea）和本波罗的海（Proper Baltic Sea）种群。利用来自59个个体的全基因组重测序（whole genome resequencing）数据，研究人员重建了波罗的海殖民的种群动态历史，并鉴定了局部适应的基因组特征。种群动态建模支持一个序贯殖民模型，其中贝尔特海种群首先从一个祖先的北海种群分化出来，随后本波罗的海种群独立分化，有两个同等支持的定年假说将初始殖民置于约尔迪亚海（Yoldia Sea）时期（约11,800年前）或波罗的海时期（约8,000年前）。本波罗的海种群表现出与奠基者效应（founder effect）一致的漂变主导的分化，而在北海与贝尔特海种群之间发现了两个离散的分化基因组岛（genomic islands of differentiation），具有选择性清除（selective sweep）信号，指向近期的正选择。这两个区域均包含渗透压调节候选基因，最引人注目的是SCNN1D，它编码上皮钠离子通道（epithelial Na+ channels）的一部分，这是钠稳态的一个关键调节因子。SCNN1D中的一个非同义突变（non-synonymous mutation）沿盐度梯度显示出显著的等位基因频率分化，并在约400个个体中得到确认。这些结果共同揭示了遗传漂变（genetic drift）与分歧选择（divergent selection）如何相互作用，驱动对半咸水环境的快速适应，并将SCNN1D确定为鲸类渗透压调节适应的一个有力候选基因。

论文解读文章

一、研究背景与研究目的

理解物种如何扩张其分布范围并适应新环境是进化生物学的核心问题。当殖民新地区涉及全新的环境条件时，新的选择压力可能驱动快速的局部适应，尤其是当选择作用于已在祖先群体中分离的现存遗传变异时。然而，殖民过程很少是纯粹的选择性过程；伴随范围扩张的种群瓶颈和收缩等种群动态事件，会严重降低有效种群大小（effective population size, Ne）并侵蚀选择所依赖的现存遗传变异。这种选择与种群动态约束之间的张力对于理解种群多样化至关重要：尽管物种范围内的环境异质性可通过将种群暴露于不同的选择压力下促进局部适应，但与种群动态历史相关的遗传漂变可能会抵消或掩盖这些适应性信号。因此，厘清这些过程需要在种群近期殖民了环境差异显著的生境的系统中进行研究，从而重建殖民的种群动态历史和快速适应的基因组印记。

波罗的海是一个理想的此类系统。它是一个地质年轻、环境动态变化的盆地，在末次盛冰期后经历了两次连续的湖-海循环，才形成了现代的半咸水波罗的海。该区域拥有世界上最强的盐度梯度之一，从北海到波罗的海内部盐度急剧下降。港湾鼠海豚（Phocoena phocoena）是该海域唯一的定居鲸类物种，沿此盐度梯度分布为三个遗传和管理上不同的种群：北海（North Sea, NOS）种群、贝尔特海（Belt Sea, BES）种群和本波罗的海（Proper Baltic Sea, PBS）种群。此前的研究表明，盐度是北大西洋港湾鼠海豚基因组变异的主要驱动力，暗示波罗的海种群经历了对低盐条件的局部适应。然而，关于这些种群如何殖民波罗的海、其种群动态历史如何，以及快速适应低盐环境的基因组基础是什么，这些问题尚不清楚。本研究旨在填补这一空白，通过全基因组测序数据，重建殖民历史并鉴定局部适应的分子足迹。

二、研究方法概述

研究人员对来自波罗的海区域的59个港湾鼠海豚样本进行了全基因组重测序（whole genome sequencing, WGS）。利用高质量SNP数据集，研究人员首先通过主成分分析和群体遗传结构分析确定了三个遗传种群。随后，使用fastsimcoal2进行种群动态历史推断，测试了13种不同的分化场景。为了量化基因组多样性和分化水平，研究人员计算了核苷酸多样性（π）、相对分化指数（Fst）和绝对分化指数（dxy）。通过三种互补的方法（Pcadapt、基于位点的Fst和环境关联分析LFMM）鉴定受选择的异常SNP，并使用跨种群扩展单倍型纯合度（XP-EHH）检测选择性清除信号。最后，研究人员开发了一个由192个高度分化的SNP组成的简化SNP分型面板，利用Fluidigm技术对更大规模的约400个样本进行基因分型和种群分配验证。

三、研究结果

3.1 遗传结构

群体遗传结构分析证实了波罗的海区域存在三个港湾鼠海豚种群：北海（NOS）、贝尔特海（BES）和本波罗的海（PBS）种群。主成分分析（PCA）和群体遗传结构分析（K=3）清晰地分离了这三个集群。研究人员依据K=3下的祖源比例（≥75%）将个体分配到各群体，最终得到12个NOS、30个BES和5个PBS个体用于后续分析。

3.2 种群动态推断

通过对13种种群动态模型的比较，研究支持一个序贯殖民模型：贝尔特海（BES）种群首先从一个祖先的北海（NOS）种群分化出来，随后本波罗的海（PBS）种群也从北海（NOS）种群独立分化，而非从地理上相邻的贝尔特海（BES）种群分化。关于殖民时间，存在两个同等支持的假说：一个将首次分化置于波罗的海时期（约8044年前），另一个则回溯到更早的约尔迪亚海时期（约11,805年前）。无论哪种假说，第二次分化（PBS从NOS分化）都发生在现代半咸水波罗的海形成之后。

3.3 基因组分化与多样性景观

全基因组范围内的Fst分析显示，NOS-BES比较中存在两个突出的高Fst峰（位于伪染色体1和15），而BES-PBS和NOS-PBS的比较中则呈现相对均匀的分化模式。尽管缺乏突出的Fst峰，PBS与其他两个群体的整体分化水平却更高，这与其较低的核苷酸多样性（π）和Watterson's θ值一致，反映了强烈的遗传漂变影响。所有种群间的绝对分化（dxy）水平相近，且大部分基因组窗口的π与dxy值接近1:1线，符合近期共祖的预期。

3.4 选择扫描

三种异常SNP扫描方法（Pcadapt、Fst和LFMM）共同鉴定出78个推定的适应性SNP。这些SNP的基因型模式显示，生活在低盐环境中的BES和PBS个体共享相同的基因型，而与高盐环境中的NOS个体截然不同。XP-EHH分析在NOS-BES比较中检测到最强的选择性清除信号，其中12个推定适应性SNP与选择性清除区域重叠。这些信号主要集中在伪染色体1和15上的两个基因组岛上。

3.5 大规模基因分型

开发的192个SNP分型面板在质量控制后保留了149个可靠SNP，并对376个独特个体成功进行了基因分型。该面板能够完美地重现全基因组数据的三个种群结构，并以高准确性分配个体。重要的是，位于SCNN1D基因内的推定适应性非同义SNP的等位基因频率沿盐度梯度呈现出显著的分化：NOS个体中C等位基因频率为0.70，而BES和PBS个体中仅为0.20，表明该位点受到持续的正选择维持分化。

四、讨论与结论

本研究通过全基因组数据揭示了波罗的海港湾鼠海豚的殖民历史和局部适应的基因组基础。种群动态建模支持一个非直观的序贯殖民模型，即贝尔特海（BES）和本波罗的海（PBS）种群分别独立地从北海（NOS）祖先分化而来，这可能反映了历史上绕过贝尔特海的扩散走廊的存在。殖民时间的两种假说均将整个分化过程置于不到12,000年的冰后期内，使得波罗的海鼠海豚成为快速进化的一个范例。

在基因组层面，研究发现两种截然不同的分化模式：NOS-BES之间存在由选择驱动的基因组分化岛，而涉及PBS的比较则表现为漂变主导的全基因组均匀分化。前者反映了在两个较大且仍有基因流的种群间，分歧选择克服了漂变的影响；后者则揭示了PBS因极小的有效种群大小（Ne）和强烈的奠基者效应，导致遗传漂变压倒选择，即使存在适应性进化也难以产生可检测的经典选择性清除信号。

本研究最重要的发现是鉴定出两个与渗透压调节相关的强候选基因组区域，特别是位于伪染色体1上的SCNN1D基因。该基因编码上皮钠离子通道（ENaC）的δ亚基，是脊椎动物维持钠稳态的关键分子。在该基因外显子8中发现的一个非同义突变，在高低盐环境种群间呈现出极端的等位基因频率差异和强烈的选择性清除信号，使其成为鲸类快速渗透压调节适应的最强候选基因。此外，该区域还富集了神经发生相关的候选基因，暗示可能存在肾脏和神经系统渗透压调节通路的协同重组。

总之，这项研究得出结论：港湾鼠海豚对波罗的海低盐环境的适应是遗传漂变与分歧选择相互作用的产物。在较大的贝尔特海（BES）种群中，选择在SCNN1D等关键渗透压调节基因上留下了清晰的印记，推动了快速的局部适应；而在较小的本波罗的海（PBS）种群中，强烈的遗传漂变掩盖了适应性信号，但其共享的低盐适应等位基因表明选择同样在起作用。SCNN1D被确定为一个极具说服力的鲸类渗透压调节适应候选基因，为理解海洋哺乳动物如何快速适应极端环境变化提供了新的视角。