基因流下的局部适应在大西洋海滨歌雀中具有多基因基础

《Ecology and Evolution》：Local Adaptation With Gene Flow Has a Polygenic Basis in a Threatened Coastal Songbird

　　

跨越异质性环境的物种通常表现出局部适应，尽管种群之间存在基因流。在这种情况下，局部适应的维持受若干因素控制，包括选择强度、迁移、选型交配和基因组结构。歌带鹀（Melospiza melodia）是一种分布广泛的鸣禽，具有多个亚种，其中包括适应海岸环境的大西洋歌带鹀（M. m. atlantica）。尽管与一个广布型亚种之间存在大量基因流，大西洋歌带鹀种群仍被发现能维持其海岸适应性表型。研究人员为该物种构建了一个新的染色体水平参考组装，并利用全基因组数据探索该亚种中基因流存在下的分化基因组结构及适应的遗传基础。大西洋歌带鹀的适应具有多基因基础，处于分化选择下的位点散布于整个基因组，而非集中在低重组连锁区块。精细尺度种群结构、有限的基因组分化以及一小部分基因上的分化选择信号提示，自然选择在维持该亚种边界中发挥作用。研究结果突显了选择、基因流和重组在塑造基因流下适应过程中的复杂相互作用，并为广泛分布、形态多样的物种中生态分化的基因组结构提供了见解。

# 论文解读：大西洋海滨歌雀在基因流下的局部适应具有多基因基础

## 一、研究背景与科学问题

物种分布范围跨越异质环境时，不同选择压力常导致种群间出现可遗传性状分化，即局部适应。当相邻的适应性分化种群间存在广泛基因流时，中性基因组背景趋于同质化，而差异分化通常仅局限于支撑局部适应的基因组区域。因此，高分辨率基因组数据是检测此类种群结构的关键，同时也为识别局部适应的遗传基础提供了独特窗口——在低分化的基因组背景下，受分化选择的位点更容易被检测。维持基因流下适应性分化的能力受选择强度、迁移程度、选型交配以及基因组物理连锁与结构互作等因素控制。单一强效应位点可能仅靠强选择即可维持，但多基因适应常与重组冷点（如着丝粒、倒位）相关，使共适应等位基因以连锁块形式遗传。然而，理论预测多基因适应在高重组条件下可更高效地筛选有益等位基因组合。歌带鹀（*Melospiza melodia*）是研究这一问题的理想系统：它包含25个已识别亚种，其中大西洋海滨歌雀（*M. m. atlantica*）特化于北美东海岸盐沼与海岸沙丘，虽与广布型亚种东部歌雀（*M. m. melodia*）存在广泛基因流，仍维持着大喙、深色条纹等“盐沼综合征”表型。此前基于简化基因组的研究虽能区分两个亚种，但分辨率不足以揭示精细结构或确定选择目标。本研究旨在通过染色体水平参考基因组与全基因组重测序，刻画基因流下分化的基因组结构，并定位海岸适应的候选基因。

## 二、研究内容与主要结论

研究人员构建了歌带鹀新的染色体水平参考组装，并对覆盖两个亚种分布区的个体进行全基因组重测序，综合运用群体结构分析、全基因组分化统计与单倍型选择扫描，系统探究了局部适应在基因流下的维持机制。结果发现，大西洋海滨歌雀的适应具有明显的多基因基础，受选择位点散布于整个基因组而非集中于低重组区块；同时检测到亚种内精细的种群结构，但基因组总体分化程度很低，表明自然选择在有限位点上维持亚种边界。这一发现强调了选择、基因流和重组在塑造适应性分化中的复杂相互作用，也为生态分化早期过程的基因组机制提供了实证。

## 三、主要技术方法

研究人员于2021—2022年繁殖季在缅因州至北卡罗来纳州的11个地点采集样本，最终纳入53只个体，分为内陆（*M. m. melodia*）、沿海大陆和南部屏障岛三类。技术流程包括：使用PacBio HiFi长读段与Hi-C数据构建染色体水平参考基因组；对每个个体进行约15×覆盖度的Illumina双端重测序；经比对、变异检测与严格过滤获得高质量SNP数据集；基于连锁不平衡剪枝和统计定相构建单倍型；采用主成分分析（PCA）和判别主成分分析（DAPC）解析群体结构；用pixy计算核苷酸多样性（π）、固定指数（FST）和绝对分化度（dxy）；利用跨群体扩展单倍型纯合度（xpEHH）统计量扫描分化选择信号，并在显著峰附近10 kbp范围内注释候选基因。

## 四、研究结果

### 4.1 参考基因组（Reference Genome）

新的歌带鹀参考组装（bMelMel2.pri.cur.20230821；GCF_035770615.1）包含353个支架、31条完整常染色体及完整的Z和W染色体，contig N50为23.12 Mbp，scaffold N50为82.36 Mbp，共注释20,371个基因。组装总长1.55 Gbp，大于近期其他歌带鹀亚种的组装（0.98–1.40 Gbp），这一差异可能源于组装技术或生物学因素，尚需进一步评估。

### 4.2 基因组数据（Genomic Data）

最终VCF包含20个屏障岛个体、15个沿海大陆个体和18个内陆个体。平均个体覆盖度为14.6×，平均基因型缺失率为2.5%。过滤后保留11,496,474个双等位SNP。连锁不平衡极低，r2在310 bp距离内降至0.10以下，300 kbp时降至0.03，表明重组率高。

### 4.3 群体结构（Population Structure）

PCA第一主成分（3.4%方差）清晰分离两个亚种，内陆个体紧密聚集，屏障岛个体位于另一极端，沿海大陆个体居中。PC2（2.4%）主要区分北卡罗来纳州皮岛与其他屏障岛种群。大西洋海滨歌雀内部存在精细结构：特拉华湾以北至纽约的沿海大陆为一个亚群；切萨皮克湾及特拉华/马里兰屏障岛为另一亚群；弗吉尼亚南部屏障岛为第三亚群；北卡罗来纳皮岛为最独特的种群。相比之下，内陆东部歌雀在整个1500 km采样范围内呈泛交配状态，未检测到地理结构。

### 4.4 全基因组分化（Genome-Wide Divergence）

亚种间平均FST仅为0.027，绝对序列分化度dxy为6.00e-03，且与两亚种的π高度相关（R2分别为0.983和0.988），但与FST几乎不相关（R2=0.011）。这些模式与近期分化及高基因流一致。大多数FST较高的窗口变异位点稀少，可能是低多样性区域而非选择驱动的分化峰。

### 4.5 选择扫描与候选基因（Selection Scan and Candidate Genes）

xpEHH比较（内陆对屏障岛、内陆对沿海大陆）共发现11个共享的显著峰，涉及20个基因。其中7个区域在大西洋海滨歌雀中显示出正选择信号，3个区域邻近着丝粒。部分基因与渗透调节功能相关，如ANO3（氯离子跨膜转运体家族成员）、S100A16、TRIP13（肾脏相关），以及预测基因LOC134425495（与CYP11A1直系同源，参与醛固酮合成和盐水平衡）。还有PCSK2（与黑色素合成相关）和COL15A1（与体型变异相关）。这些基因中的多数此前未被其他海岸-内陆鸣禽比较研究报道，表明平行适应可能通过谱系特异机制实现。

## 五、讨论总结与结论

讨论部分指出，全基因组数据揭示了此前简化基因组方法无法分辨的精细种群结构：大西洋海滨歌雀在约1000 km的繁殖范围内形成多个亚群，而东部歌雀在更大范围内保持泛交配，这与其高扩散能力和不可预测的迁徙行为一致。分化位点散布于整个基因组，而非聚集在倒位等连锁块中，提示强选择在缺乏选型交配的情况下足以维持多基因局部适应。少数选择峰近着丝粒表明重组景观仅起次要作用。FST峰与选择信号不一致的原因可能包括极低连锁不平衡、不完全扫荡及基因流阻碍固定。此外，两个亚种可能处于二次接触后的早期阶段，尚未达到迁移-选择平衡。遗传基础方面，尽管多数候选基因为谱系特异性，但TRIP13和PCSK2在沼泽麻雀的海岸-内陆分化中同样被检测到，表明共享与特异机制并存。

研究结论可概括为：大西洋海滨歌雀在持续基因流下维持局部适应的遗传基础是多基因的，受选择位点分散于高重组基因组背景中，自然选择作用于少量基因而非大型连锁区，从而维持亚种界限；该结果为理解广泛分布、形态多样物种中生态分化的基因组架构提供了关键证据。