《Ecology and Evolution》：Genetic Structure and Signs of Local Adaptation in Icelandic Brown Trout Populations
编辑推荐：
鲑科鱼类溯河洄游（anadromy）的丧失通常发生在冰川后定殖与栖息地变化之后，这一过程受选择压力驱动，此时留居（residency）的收益超过迁移成本。约一万年前末次冰期结束时，溯河洄游褐鳟抵达冰岛并入侵新形成的栖息地。一些群体在冰层退缩后不久即被隔离于河源
鲑科鱼类溯河洄游（anadromy）的丧失通常发生在冰川后定殖与栖息地变化之后，这一过程受选择压力驱动，此时留居（residency）的收益超过迁移成本。约一万年前末次冰期结束时，溯河洄游褐鳟抵达冰岛并入侵新形成的栖息地。一些群体在冰层退缩后不久即被隔离于河源区，推测加速了局部适应和留居习性的演化。在本研究中，研究人员在全国尺度上评估了褐鳟的遗传亲缘关系与多样性，以解析独特群体的数量及其连通性。研究人员检测了来自河源区的基因流，并评估了隔离与适应的影响。研究人员利用ddRADseq技术对来自40个采样点的600个个体进行基因分型，鉴定出2946个推定中性SNP。研究人员识别出16组遗传相似个体（遗传簇），其中5组主要由溯河洄游个体组成。在大尺度上，褐鳟群体的遗传变异与亲缘关系反映了水系的地理格局，而在流域内部，地形（如瀑布和湖泊的存在）则是关键因素。与直觉相反，来自隔离河源群体的下游基因流仅在小河群体中观察到，而来自大型湖泊群体的基因流似乎可忽略不计。利用离群SNP，研究人员鉴定出数百个与温度和栖息地类型（即湖泊与河流）相关的位点。与温度相关的标记定位到与免疫应答、脂质代谢及跨膜蛋白基因相关的位点上。研究结果凸显了隔离对河源群体遗传多样性的影响，以及识别独特群体在保护工作中的重要性。
**论文解读：冰岛褐鳟群体的遗传结构、局部适应及其保护意义**
**一、研究背景与科学问题**
鲑科鱼类普遍存在溯河洄游（anadromy）与淡水留居（residency）两种生活史策略的共存现象。溯河洄游个体需承担长距离迁移、盐度转换生理压力及海洋捕食风险等高昂成本，但可通过海洋环境中的快速生长和更高繁殖力获得补偿；留居策略虽规避了迁移风险，却受限于淡水栖息地的生长潜力和繁殖输出。理论上，两种策略的维持取决于其成本与收益的平衡。当自然或人为障碍（如不可逾越的瀑布）形成后，海-淡水迁移受阻，留居表型将受到正选择，伴随早熟、迁移行为减弱及渗透调节能力改变等适应性演化。
冰岛位于大西洋中脊，末次冰期后经历了频繁的火山活动与地壳均衡运动。约10,000 1?C年前，海平面较当前高约40 m，溯河洄游褐鳟（*Salmo trutta*）与北极红点鲑（*Salvelinus alpinus*）随冰川消退而入侵新形成的水系。随后，地壳均衡回弹、冰川河流流量减少、构造活动及火山喷发塑造了大量不可逾越的瀑布，将众多群体隔离于河源区。尽管冰岛褐鳟线粒体DNA多样性较低，但其表型变异极为丰富，涵盖地热溪流中的侏儒型个体（最大体长仅30 cm）、大型捕食性个体以及不同洄游习性的群体。
以往研究多聚焦于局部地理尺度的褐鳟遗传结构，全国尺度的系统性研究仍属空白。本研究旨在回答三个核心问题：（i）地理格局与水系物理演化（如瀑布形成）如何塑造冰岛褐鳟的冰川后演化与分化；（ii）群体隔离如何通过遗传漂变与局部适应导致遗传多样性丧失；（iii）是否存在与温度及栖息地类型（湖泊与河流）相关的局部适应基因组信号，以预测鲑科鱼类对未来气候变化的响应。
**二、主要技术方法**
研究人员利用双酶切限制性位点关联DNA测序（ddRADseq）技术，对来自冰岛全国58个采样点（13个湖泊、1个潟湖、25条河流）的1,120尾褐鳟进行基因分型，并以苏格兰斯凯岛Loch Slapin的溯河褐鳟作为参考群体。经过严格的数据筛选——包括去除旁系同源位点、一级亲缘关系个体、偏离哈迪-温伯格平衡（Hardy-Weinberg equilibrium）的位点以及缺失数据过多的个体——最终保留600个个体的2,946个推定中性SNP用于多样性与亲缘关系分析，另保留5,614个标记用于环境关联分析。
数据分析流程包括：利用Stacks v2.4进行 reads 处理与位点组装，参考基因组为*S. trutta*基因组（RefSeq assembly accession: GCF_901001165.1）；通过主成分分析（PCA）、最大似然系统发育树（RAxML-NG）、祖先成分分析（snmf）及F_ST统计评估遗传结构；采用潜在因子混合模型（LFMM）检测等位基因频率与温度的关联；通过比较湖泊与河流群体及成对湖泊与溯河群体的F_ST值鉴定与留居习性相关的离群位点。
**三、研究结果**
**3.1 冰岛褐鳟分为16个遗传簇，河源群体与下游群体显著分化**
祖先成分分析（K=17，包括苏格兰参考群体）揭示16个冰岛褐鳟遗传簇（C1–C16），其中5个主要由溯河洄游个体组成。聚类模式在流域尺度上反映水系地理格局，如?lfusá河流域形成7个独立簇（C1–C7）；南部推定溯河群体聚为单簇（C10），Faxaflói湾西北侧河流群体聚为C11；北部Svínavatn和Myvatn分别形成独立簇（C12、C15）；西峡湾区形成南北两个簇（C13、C14）。河源群体与下游群体在所有分析中均表现出极强分化，且大型湖泊（Tingvallavatn、Myvatn）向下游的基因流贡献可忽略不计，而小河群体（如Tverá向Tingvallavatn、Forv?e向Skjálftavatn）存在下游基因流证据。
**3.2 系统发育树支持河源群体的独立演化地位**
最大似然系统发育树显示，苏格兰褐鳟与冰岛褐鳟位于最深的分支；冰岛分支中，Veieiv?tn/Tórisvatn群体构成高支持率的单系群。多个河源隔离群体（如Tverá、Leirvogsvatn、Svínavatn、Myvatn、Forv?e）位于高支持率（85%–100%）的末端分支上，而推定溯河群体多位于支持率较低的内部支系，与其广泛的基因交流历史一致。
**3.3 多数河源群体遗传多样性显著降低**
核苷酸多样性（π）分析显示，大多数河源褐鳟群体的遗传多样性显著低于溯河群体，后者与苏格兰参考群体水平相当。多样性沿河流向上游递减，如Varmá至Hengladalsá（Innstidalur）π值从6.82×10??降至4.68×10??。长期隔离的湖泊群体（Tingvallavatn/úlfljótsvatn、Myvatn、Stóra Fossvatn）同样呈现低多样性（π=3.34–4.38×10??），而栖居于瀑布上方小生境中的Tverá和Leirvogsvatn群体则表现出全冰岛最低值。等位基因丰富度（AR）方面，溯河群体AR=1.31–1.87，而隔离群体低至1.13–1.39。值得注意的是，有记录进行人工放流的湖泊（除Stóra Fossvatn外）显示出与溯河群体相当的AR值，反映放流活动的遗传效应。
**3.4 多个位点等位基因频率与栖息地温度相关**
以7个具有夏季温度记录（5–10月平均值，6.5°C–15.7°C）的采样点为解释变量，LFMM分析在假发现率q=0.001水平鉴定出40个与温度显著相关的标记（q=0.05时为461个），分布于22个连锁群，其中19个位于编码区、21个邻近（<5000 bp）编码区。候选基因功能涉及跨膜蛋白（如12号染色体富亮氨酸重复跨膜神经元蛋白4类似物）、脂质代谢（如15号染色体长链脂肪酸-CoA连接酶6类似物）及免疫应答（如2号染色体锌指CCHC域蛋白2类似物）。多数位点的稀有等位基因频率随温度升高而增加，呈现多基因适应的特征。
**3.5 与溯河洄游丧失相关的基因组分化**
将5个湖泊群体（?xará、Stóra Fossvatn、Myvatn、Svínavatn、Hestvatn）分别与推定溯河群体（Sog河）进行成对F_ST比较（阈值>0.5），鉴定出的离群标记数量为：?xará 214个、Stóra Fossvatn 155个、Myvatn 132个、Svínavatn 53个、Hestvatn 10个。60%–66%的推定适应标记为3个隔离湖泊所独有，仅2个标记在所有5个湖泊中均出现分化。功能注释揭示这些标记涉及渗透与热应激响应（3号染色体双特异性激酶）、细胞骨架组织（7号染色体Coronin肌动蛋白结合蛋白）、细胞分裂（14号染色体stathmin-3类似物）及血管重构（血管生成素2b）等与洄游生理相关的基因。将参考群体替换为苏格兰Loch Slapin时得到相似模式。
**四、讨论与结论**
本研究揭示了冰岛褐鳟遗传结构的两级驱动机制：在大尺度上，水系地理格局决定群体间亲缘关系；在流域尺度内，地形障碍（瀑布、湖泊）塑造分化模式。16个遗传簇的存在反映了冰岛复杂的冰川后定殖历史，支持"一次延长迁移波"假说——早期到达个体被隔离于河源区后经历独立演化。河源群体普遍低遗传多样性与高遗传分化的特征，是遗传漂变与局部适应长期作用的必然结果。
关于下游基因流，研究发现其可能性取决于隔离栖息地的大小与稳定性。小河群体（Tverá、Forv?e）存在向下游的基因流贡献，而大型湖泊群体（Tingvallavatn、Myvatn）则未检测到明显下游输出。这一不对称格局可能源于：湖泊生态系统提供充足的生长机会，促使被困个体转向淡水迁移模式（如在湖中摄食、在支流繁殖），从而降低向海迁移的倾向；而小河群体因资源有限，仍保留部分向海迁移的遗传变异。
温度关联分析鉴定的候选基因具有明确生物学合理性：免疫基因对温度敏感已获广泛证实，脂质代谢与膜流动性调节是变温动物温度适应的经典机制。但研究者也指出，温度与地理（海拔、纬度）的共线性可能混淆因果关系，且有限的温度采样站点降低了统计效力。
关于溯河洄游丧失的遗传基础，研究表明不同湖泊群体的离群标记重叠度有限，暗示遗传平行性程度较低，这与北极红点鲑、三刺鱼和白鲑的研究发现相似。?xará（Tingvallavatn产卵群体）的高分化标记数可能反映该群体早期定殖后的长期独立演化历史。
本研究对于保护生物学具有重要启示：鉴定出的河源独特群体（如Tverá、Leirvogsvatn、Stóra Fossvatn）遗传多样性极低且基因库独特，应优先保护以维持其作为演化储库的功能；人工放流活动虽然短期内提高了部分湖泊的遗传多样性，但也引入了混合祖先成分（如Veieiv?tn系统除Stóra Fossvatn外的湖泊均含未知来源的C9成分），可能掩盖自然遗传模式并造成基因渗入风险。未来应利用全基因组重测序技术，结合更多群体的复制研究，进一步阐明温度适应与洄游丧失的遗传机制，为鲑科鱼类保护提供科学依据。
**研究结论**
本研究描述了40个采样点中的16个独特遗传簇，揭示了隔离分化与局部适应的证据。河源独特群体（?xará、Tverá、上Hengladalsá、Leirvogsvatn、Stóra Fossvatn和Forv?e）相比溯河群体表现出显著降低的遗传多样性，这一模式可能由数千年的适应与漂变所驱动。在水系内部，群体结构反映了水道拓扑格局，下游基因流证据仅局限于小型隔离河流群体（如从Tverá到Tingvallavatn）。出乎意料的是，来自大型河源湖泊（Tingvallavatn和Myvatn）的下游基因流似乎可忽略不计，凸显了栖息地类型对连通性的差异影响。人类干预在Veieiv?tn系统中的影响清晰可见——经放流的湖泊显示出混合祖先群体，而未放流的Stóra Fossvatn则保持独特遗传特征。此外，研究结果还揭示了局部适应信号，部分标记与温度或留居程度（湖沼型褐鳟与溯河河川型褐鳟）显著相关。全基因组测序与更多群体的重复研究将有助于进一步阐明温度制度和溯河洄游变异相关局部适应的遗传基础。利用如本研究所采用的遗传分析识别独特群体，对其管理具有重要意义，并可促进其保护。显然，应优先保护遗传多样性有限的河源群体（如Tverá、Leirvogsvatn），以维持其作为演化储库的功能，同时对人工放流湖泊进行基因渗入特征监测，以避免本地基因库的丧失。
