硫磺泉鱼类种群性状表达的共享与独特模式——以哥斯达黎加瓜纳卡斯特为例

《Environmental Biology of Fishes》：Shared and unique patterns of trait expression in a sulfide spring fish population in Guanacaste, Costa Rica

　　

摘要 当生物体暴露于相似的环境压力时，往往表现出相似的特征。在富含有毒硫化氢（H2S）且严重缺氧的硫磺泉中，花鳉科（Poeciliidae）成员已多次演化出相似的生活史、生理和形态特征。本研究量化了来自哥斯达黎加瓜纳卡斯特一处硫磺泉的吉利花鳉（Poecilia gillii）种群中与硫化氢耐受相关的形态和生理特征，并与邻近的非硫磺种群以及墨西哥南部研究较为深入的墨西哥花鳉（P. mexicana）种群进行了比较。尽管在P. gillii中发现的大多数差异表达基因在硫磺环境中的P. mexicana种群中并未出现差异调控，但确实存在若干与硫化氢解毒、硫转运及细胞缺氧反应相关的差异表达基因与P. mexicana共享。此外，与硫磺环境中的P. mexicana相似，来自硫磺泉生境的P. gillii相较于非硫磺鱼类具有更长的头部、更延长的腹部以及减小的标准体长（standard length）。此前在墨西哥南部的研究表明，现存遗传变异（standing genetic variation）是硫磺泉鱼类演化的重要驱动因素，但本研究强调，即使由于地理隔离导致祖先现存遗传变异不太可能共享，相似的选择梯度仍可导致相似特征的产生。未来利用共同园（common-garden）饲养环境对P. gillii开展的研究，将需要探讨这些结果是否代表对硫磺泉的趋同适应。

《环境生物学中的趋同适应：哥斯达黎加硫磺泉吉利花鳉的多维度性状研究》解读

一、研究背景与科学问题

极端环境为研究生物适应性演化提供了天然实验场。硫磺泉是一类同时存在两种主要环境胁迫因子的极端生境：高浓度有毒硫化氢（H₂S）和严重缺氧。H₂S在微摩尔浓度下即可对多数后生动物产生致死效应，其毒性机制在于直接抑制电子传递链中的细胞色素c氧化酶（cytochrome c oxidase, COX），干扰氧化磷酸化过程。与此同时，H₂S的解毒过程需要消耗氧气，且H₂S在水溶液中会自发氧化消耗溶解氧，进一步加剧缺氧状态。尽管环境极为严苛，花鳉科鱼类中墨西哥花鳉（Poecilia mexicana）物种复合群已在墨西哥南部多个硫磺泉中成功定殖并适应。

此前对墨西哥南部P. mexicana种群的研究揭示了硫磺泉鱼类在基因表达、形态、生理等多层面的适应性特征：硫磺生境种群普遍上调硫化物解毒通路相关基因（如sqr、ethe1、slc25a10），表现出更小的体型、更大的头部、以硫细菌和无脊椎动物为主的食性、较少但较大的后代以及降低的代谢率。然而，这些种群间的相似性可能部分归因于基因流和现存遗传变异（standing genetic variation）上的选择，这些因素在地理隔离的种群间不太可能发挥作用。

哥斯达黎加瓜纳卡斯特地区硫磺泉中的吉利花鳉（Poecilia gillii）种群提供了一个独特的研究机会：由于与墨西哥种群之间存在巨大的地理隔离，共享祖先变异和基因流的可能性极低，因此可以检验相似环境是否独立驱动相似性状的演化。本研究旨在回答这一核心科学问题：在地理隔离、祖先变异不共享的条件下，硫磺泉环境是否仍能导致趋同的适应性特征？

二、主要技术方法

研究人员采用了多维度整合的研究策略。样本采集自哥斯达黎加Río Tempisque水系内的一处硫磺位点（10.838°N, 85.516°W）和邻近非硫磺位点（10.851°N, 85.485°W），分别于2015年（形态数据）和2024年（基因表达、耐受性数据）采集。关键技术方法包括：几何形态测量学分析（geometric morphometrics），基于16个地标点量化体形变异；急性H₂S和缺氧耐受性测定，采用生存分析（Cox回归）评估耐受差异；鳃组织转录组测序（RNA-seq），并与已发表的墨西哥P. mexicana三个水系（Pichucalco、Puyacatengo、Tacotalpa）的转录组数据（PRJNA290391）进行整合比较；差异表达分析使用edgeR包，基因本体（Gene Ontology, GO）富集分析使用topGO包。实验鱼在实验室常氧非硫磺条件下随机繁殖约20个月后进行耐受性测试。

三、研究结果

形态变异

标准体长分析显示，生境和性别均显著影响体长。硫磺种群鱼类平均比非硫磺种群短49%，雄性平均比雌性短17%。几何形态测量学分析表明，PC1轴解释了62.4%的变异，主要反映性别二态性；PC2轴解释了12.2%的变异，反映生境类型差异。Procrustes ANOVA显示生境对体形有显著影响，但该效应依赖于质心大小（centroid size），且两个种群体长范围几乎不重叠，因此体形差异可能受强烈的异速生长效应影响。与墨西哥硫磺泉鱼类一致，硫磺生境的P. gillii雌雄个体均表现出更长的头部和更延长的体形。

急性H₂S和缺氧耐受性

H₂S耐受性分析显示生境与性别之间存在显著交互作用：硫磺种群雄性的硫化物耐受性显著高于非硫磺种群雄性，但雌性间差异不显著。缺氧耐受性分析发现生境、性别和标准体长之间存在显著的三方交互作用，其中较大的硫磺种群雌性表现出增强的缺氧耐受性。总体而言，种群间理化耐受性差异不如墨西哥P. mexicana种群间差异显著，这可能表明表型可塑性在P. gillii的胁迫耐受性塑造中发挥更强作用，或非硫磺种群本身具有较高的基础耐受水平。

硫磺与非硫磺P. gillii的基因表达**

鳃转录组测序共获得2.076亿条读段，过滤后保留18,303个基因用于下游分析。硫磺种群上调的3,621个基因富集于134个GO生物学过程条目，包括网格蛋白依赖性内吞作用、硫化合物代谢过程、细胞对热反应的调控、三羧酸循环、糖异生和经典糖酵解等。下调的3,190个基因富集于105个GO条目，主要涉及细胞质翻译、炎症反应和免疫功能。

跨种群基因表达比较

多维标度（MDS）分析显示，基因表达总体变异沿MDS1轴按物种分离（哥斯达黎加P. gillii vs. 墨西哥P. mexicana），沿MDS2轴按生境类型分离。差异表达分析显示，P. gillii具有最多的独有差异表达基因（2,300个上调和1,822个下调）。四个水系共有的上调基因富集于40个GO条目，包括H₂S解毒（GO:0070813硫化氢代谢过程、GO:0006749谷胱甘肽代谢过程）、硫化合物转运与加工、氧化磷酸化（GO:0006123线粒体电子传递，细胞色素c到氧）、糖异生、经典糖酵解以及细胞对缺氧的反应（GO:0071456）。共有的下调基因富集于14个GO条目，涉及角质化和脂质代谢。P. gillii独有上调基因富集于131个GO类别，包括染色质重塑、应激颗粒组装、蛋白质稳定化、内吞作用和泛素依赖性蛋白分解过程；独有下调基因富集于29个GO类别，涉及细胞质翻译和突触前胞质钙浓度调节等。

四、讨论与结论

本研究表明，哥斯达黎加硫磺泉中的P. gillii在多个性状维度上与墨西哥硫磺泉P. mexicana种群表现出趋同特征：上调硫化物解毒、氧化磷酸化和细胞缺氧反应相关关键基因，同时表现出更长的头部和更小的体型。这些相似性出现在地理隔离、祖先共享变异和基因流均不太可能的情况下，表明相似的选择压力梯度可以独立驱动相似的适应性演化。

形态差异方面，硫磺泉鱼类较小的体型可能是一种降低能量需求的适应策略，而较大的头部可能有助于增加口腔容积以提高通气效率。然而，由于种群间体长范围几乎不重叠且存在显著的异速生长效应，体形差异的解释需保持谨慎。

理化耐受性方面，实验室饲养条件下种群间差异不如墨西哥种群显著，可能反映可塑性在P. gillii耐受性塑造中的重要作用，或非硫磺种群本身具有较高的基础耐受水平。

基因表达方面，四个硫磺种群共享的差异表达基因集中于H₂S解毒、硫转运、缺氧反应和能量代谢调整，表明这些是硫磺泉适应的核心通路。P. gillii大量独有差异表达基因主要涉及保守的细胞应激反应（如蛋白质修复与稳定、染色质重塑、应激颗粒组装），这可能表明该种群的定殖历史较墨西哥种群更近，适应尚未完全推进；也可能反映不同硫磺泉间环境特征的差异，包括温度、pH、电导率以及潜在的人为 stressors 如农业污染和重金属等。

本研究的核心结论是：即使祖先遗传变异不共享，相似的环境选择梯度仍可导致趋同的适应性特征。未来的共同园饲养实验将有助于区分遗传适应与表型可塑性的相对贡献，进一步验证这些特征是否代表对硫磺泉环境的真正趋同适应。该研究发表于《Environmental Biology of Fishes》，为理解极端环境下的重复性演化提供了跨地理尺度的新证据。