一种恢复力强的珊瑚物种沿富营养化梯度的弱遗传结构及环境关联位点

《Ecology and Evolution》：Weak Structure and Environment-Associated Loci Across a Eutrophication Gradient in a Resilient Coral Species

　　

在海洋物种中，大的有效种群规模和高基因流（connectivity）可最小化遗传漂变（drift）并维持遗传多样性。然而，这种特性如何影响其适应人为效应所致环境变化的能力尚不明确。石珊瑚（Scleractinian corals）是面临沿海污染和气候变化日益严峻威胁的关键物种（keystone species）。在本研究中，研究人员利用限制性位点关联DNA测序（RADseq）检查了扁枝滨珊瑚（Oulastrea crispata）的种群结构。该物种为壳状珊瑚，原产于印度-太平洋地区，分布范围广泛，且表现出高度的环境胁迫耐受性。为评估扁枝滨珊瑚的潜在适应性差异，研究人员在香港范围内沿一个空间范围较小但显著的富营养化梯度（溶解无机氮：1.6–9.0 μM）进行了个体采样。一组39,527个基因组范围内的单核苷酸多态性位点（SNPs）显示，尽管环境条件存在变异，但90个个体数据集中未检测到广泛的的中性种群结构。多样性水平估计和低近交系数表明，尽管环境具有挑战性，珊瑚种群并未出现多样性丧失。虽然F^ST离群分析未能揭示跨位点的任何明确选择信号，但基因型-环境关联分析鉴定出137个与多个水质参数相关的离群SNPs。对这些离群位点的基于转录组的注释表明其可能在环境响应中发挥作用。研究人员还发现，在所采样的富营养化梯度范围内，扁枝滨珊瑚始终与环境胁迫耐受性共生藻属Durusdinium sp.相关联。尽管核位点呈现泛群交配（panmixia），环境关联等位的存在与由当地环境条件驱动的适应性反应相容。胁迫耐受性Symbiodiniaceae物种的检出提示了一种互补的全息生物（holobiont）水平机制。扁枝滨珊瑚的表型可塑性（phenotypic plasticity）可能是其胁迫恢复力的另一机制。未来需要额外的全基因组测序和转录组分析研究来阐明扁枝滨珊瑚宿主恢复力的遗传结构，以及宿主遗传学和全息生物介导效应的相对贡献。

**研究背景与意义**

海洋物种通常具有大的有效种群规模和高基因流，这有助于维持遗传多样性，但这也引发了一个关键问题：这种遗传特征如何影响它们适应人为环境变化的能力。石珊瑚作为海洋生态系统中的关键物种（keystone species），正面临气候变化和沿海污染的双重威胁，包括极端温度、海洋酸化、盐度变化和沿海富营养化等。这些胁迫因素协同作用，加剧了珊瑚生存的严峻形势。与此同时，全球约有9.97亿人口居住在距珊瑚种群100公里以内的沿海地区，珊瑚礁生态系统为人类聚居地提供了至关重要的生态系统服务。然而，沿海城市化带来的污染、栖息地破坏和过度捕捞正严重威胁着这些脆弱的生态系统。

边际珊瑚礁生境（marginal coral reef habitats）指存在于次优环境条件下的生境，为研究多重胁迫因素对珊瑚生存的影响提供了宝贵的天然实验室。香港的水域受到城市化、富营养化和自然河口影响，水质指标在数十公里范围内呈现显著梯度变化：西部水域因人为污染源和珠江三角洲径流影响，溶解无机氮水平和浊度升高，盐度波动较大；而东部水域则因持续的海洋输入保持相对寡营养状态。这种溶解氮、盐度和浊度的局部差异可能在短距离内施加不同的选择压力，为研究人员研究高基因流环境下的本地适应提供了理想条件。扁枝滨珊瑚（Oulastrea crispata）是一种环境耐受性强的斑马珊瑚，分布范围从热带印度-太平洋延伸至温带区域的济州岛和日本南部，甚至作为引入种出现在地中海。在香港，该物种主导着浑浊的潮间带和潮下带、海堤和遮阴缝隙，包括少数其他珊瑚物种难以生存的区域。除了环境耐受性，扁枝滨珊瑚还具有生殖多样性：作为雌雄同体珊瑚，它既能进行有性繁殖（ broadcast spawning/体外排放和broading/体内保育），也能进行无性繁殖（触手内和触手外出芽），并且具有营养灵活性，能够从依赖内共生藻获取营养转变为完全异养。这些生活史特征与其环境耐受性相结合，使其成为探索边际条件下珊瑚恢复力的潜在模式生物。基于此，研究人员假设扁枝滨珊瑚在香港强环境梯度中的成功反映了其对特定胁迫因素（包括污染和其他人为影响）的本地适应。

该研究发表于《Ecology and Evolution》。

**关键技术方法**

研究人员于2018年8月在香港西-东水质梯度沿线的8个地点采集了90个扁枝滨珊瑚个体样本，样本量每点4–17个不等。环境数据来源于香港特别行政区环境保护署（EPD）的月度监测数据，采用最近环境监测站的地表水参数（1米深度）作为各采样点环境条件的替代指标，涵盖九项水质参数。基因组DNA经提取后，送至Floragenex公司构建PstI酶切的RADseq文库，随后在美国俄勒冈大学GC3F设施的NovaSeq 6000平台进行单端122 bp测序。数据分析采用de novo Stacks流程进行RADseq组装，通过Kraken2去除非珊瑚来源的污染物后，最终获得39,527个双等位SNPs用于后续分析。种群结构分析采用主成分分析（PCA)、ADMIXTURE分析和分子方差分析（AMOVA），并计算观察杂合度、期望杂合度、核苷酸多样性(π)和近交系数(F_is)等多样性指标。选择压力检测结合使用BayeScan、OutFLANK、潜在因子混合模型(LFMM)和冗余分析(RDA)四种方法，仅保留多种方法共同鉴定的离群位点。RDA离群位点通过比对扁枝滨珊瑚转录组进行功能注释。此外，研究人员还将RADseq读段映射至11种Symbiodiniaceae参考基因组，通过Bray-Curtis相异矩阵和置换多元方差分析(PERMANOVA)评估共生藻组成与水质类别的关联。

**研究结果**

**同质性种群结构伴细微局部分歧**

基于39,527个不连锁SNPs和88个个体（移除两个高亲缘关系个体后）的主成分分析显示，扁枝滨珊瑚在香港范围内呈遗传同质性分布，前两个主成分分别仅解释2%和1%的总变异。ADMIXTURE分析进一步支持这一遗传均一性，最优K值为1（交叉验证误差=0.55），表明种群内无明显的遗传结构分化。

**种群多样性估计**

各采样点的异质性和核苷酸多样性差异极小。观察杂合度(0.00589–0.00664)始终低于期望杂合度(0.00865–0.0134)，核苷酸多样性(π)从最小样本量地点的0.01095到最大样本量地点的0.01417不等。AMOVA分析表明，数据集中71.16%的遗传变异存在于个体间，28.76%存在于种群内个体间，而种群间比较仅占总变异的0.08%。近交系数(F_is)在各点一致，范围从0.01342至0.02405，表明采样点内近交证据有限。成对F_ST值极低（范围：-0.0070至0.0060），进一步支持采样点间不存在种群结构的结论。

**离群分析**

BayeScan和OutFLANK未检测到任何离群位点，LFMM仅检出1个离群位点，但因三种方法间无重叠，该位点未作为候选位点报告以减少假阳性。基于RDA基因组扫描方法，研究人员检出137个离群SNPs，经马氏距离分布校正基因组膨胀因子和错误发现率（FDR）调整后具有统计显著性（q值阈值为0.05）。这些离群SNPs主要与pH（43个）、悬浮固体（39个）、温度（25个）、溶解氧（20个）和叶绿素（10个）相关联，尽管变量间存在共线性。置换检验显示整体RDA模型显著不同于零分布（F=1.0351，p=0.001），解释总变异的5%；第一约束轴(RDA1)显著（F=1.1768，p=0.001），RDA2-5轴则不显著。

在137个离群位点中，95个成功映射至扁枝滨珊瑚转录组的contigs。其中，93个(97.9%)返回blast hit，76个contigs(80%)获得GO和InterPro GO术语注释。进一步分析确定54个独特GO术语，涵盖21个生物学过程、20个分子功能和13个细胞组分GO术语。功能注释表明，候选离群位点主要关联于：（i）外部信号感知与转导（如转录调控、蛋白质磷酸化、信号转导）；（ii）膜运输和离子稳态，特别是Ca²⁺相关运输；（iii）内质网/高尔基体相关蛋白质转运和细胞代谢。

**扁枝滨珊瑚内的Symbiodiniaceae内生藻**

在调查的11种Symbiodiniaceae物种中，Durusdinium trenchii是整个水质梯度（低、中、高富营养化水平）中 representation 最高的种类。映射到Symbiodiniaceae物种的读段计数相对比例在三类水质间存在显著差异（PERMANOVA：F=3.324，p=0.013）。然而，鉴于Durusdinium与其他Symbiodiniaceae之间的读段数量存在巨大差异，这种宿主内物种组成变异可能反映稀有物种读段数量的微小变化，功能相关性尚不明确。

**讨论总结**

研究人员指出，尽管扁枝滨珊瑚在记录详尽的环境梯度上广泛分布，但未检测到遗传变异的空间结构。然而，通过RDA分析鉴定出137个环境关联的候选离群位点。鉴于西部水域以较高溶解无机氮、叶绿素a、磷酸盐和较低盐度为特征，这些应对珊瑚生存施加选择压力的因素，任何适应性反应（如果存在）仅限于少数位点而非广泛的基因组分歧。各点观察杂合度低于期望杂合度的模式与遗传漂变或扁枝滨珊瑚历史种群规模缩减一致，而非近交或种群结构所致，尽管RADseq的方法特异性偏差可能导致杂合度低估，结果解释需谨慎。

极低的中性种群结构程度与长浮游幼体期和泛化特征介导的高基因流物种模式相似。扁枝滨珊瑚既能进行体外排放又能进行幼体保育，但低F_is值表明异交很可能是主要繁殖模式。香港西南季风期间，沿海水流可在不到一周内输送水体超过50公里，使幼虫能够到达并定殖于环境迥异的位置。这种被动扩散可能是所有采样点间广泛遗传混合的基础。尽管存在这种同质性，RDA模型捕获了5%总变异的潜在离群位点子集，其注释功能涉及外部刺激检测和基因表达变化，指向扁枝滨珊瑚广泛环境耐受性的潜在机制。

方法间的不一致性可能反映多基因（多位点）效应，因为RDA等多元方法可能比LFMM等单变量方法更好地捕捉弥散的多位点信号。然而，研究人员强调需谨慎解释：首先，环境值在站点水平分配，有效重复数为站点数而非个体数，这减少有效自由度并增加假重复风险；其次，理论研究表明高连接性系统中多基因本地适应不太可能；第三，样本量和站点重复有限，研究可能缺乏解析多基因适应微妙协调等位基因频率变化所需的统计效力。

离群位点的功能注释提示其与扁枝滨珊瑚环境响应细胞通路的关联。信号转导、离子运输、膜转运和代谢等过程已知有助于细胞稳态和钙化相关过程。蛋白质磷酸化和碳水化合物代谢涉及刺胞动物-藻类共生调控，可能反映低光照或富营养/浑浊条件下的能量预算变化。Ca²⁺稳态和阳离子运输术语与珊瑚钙化过程一致，而Ca²⁺稳态破坏与珊瑚应激反应相关。

关于未来遗传研究的关键问题，是这些适应性状的遗传结构。尽管高连接性可能抑制本地适应，但理论上大效应位点或降低重组的遗传结构可支持适应。当前RADseq数据集未能在检视的空间尺度上检测到大效应位点或这种结构的明确证据。适应性变异可能由"软"选择清扫（soft selective sweeps）驱动，但这难以通过简化基因组测序分辨。

扁枝滨珊瑚 tolerance 宽广环境条件的适应性可能部分源于其与Durusdinium属的稳定关联。Durusdinium trenchii是映射读段中占比最高的Symbiodiniaceae物种，但研究承认仅有一种Durusdinium代表物种限制了种级分辨率。根据先前研究，热带地区O. crispata的优势共生菌株倾向于D. eurythalpos，温带地区则为D. boreaum。除Durusdinium外，其他Symbiodiniaceae物种组成的细微但显著差异提示需进一步实验探究。细菌微生物组在珊瑚健康和营养处理中也发挥关键作用，但先前16S宏条形码研究显示扁枝滨珊瑚细菌群落跨香港采样点相对稳定。

有限的生理学研究已证明扁枝滨珊瑚个体可表现出对广泛环境条件的高度耐受性，可能通过表型可塑性实现。在其分布极北边缘，该物种可耐受7-10°C低温，并在实验室中展现出对短期温度波动的acclimate能力。香港的营养研究显示其高度异养性，实验室中个体可在黑暗中长时间存活，暗示了浑浊低光环境下的生存途径。这种生理灵活性允许种群持续存在，可直接促进适应，也可间接增加自然选择作用于适应性状的时间，但也可能通过减少可遗传变异来限制其气候变化适应能力。

研究最终结论表明：尽管扁枝滨珊瑚在其分布范围内表现出高连接性和弱中性结构，该研究仍通过RDA检测到与环境变异相关的候选离群位点。除宿主基因组响应外，其他机制可能贡献其恢复力，包括与胁迫耐受性内生藻的稳定关联、泛化灵活的营养策略，以及允许个体在不同胁迫 regime 中存活的微生物群落。这种策略在快速变化或边际礁环境中可能提供选择优势，因为本地适应在此类环境中不太可行或过于缓慢。这些生活史特征和生理灵活性凸显了其作为边际条件下珊瑚恢复力研究潜在模式生物的价值。未来研究应采用全基因组分辨率检测和验证适应性位点，并结合幼虫运输的生物物理模型，以阐明海洋学过程如何塑造种群连接性，从而将扁枝滨珊瑚定位为研究遗传和非遗传机制如何贡献珊瑚面对气候变化和沿海干扰恢复力的有用系统。