《Journal of Biogeography》：Links Between Lagrangian Connectivity, Thermal Thresholds and Genetic Structure in the Floating Marine Siphonophore Physalia
编辑推荐：
**目的**：研究人员旨在探究遗传上不同的僧帽水母属（*Physalia*）物种是否在扩散行为和环境暴露方面存在差异，并评估大尺度风－海洋输运在塑造其遗传结构中的作用。
**位置**：热带和亚热带海域，即*Physalia*搁浅事件与分布范围重叠的区域。
**目的**：研究人员旨在探究遗传上不同的僧帽水母属（*Physalia*）物种是否在扩散行为和环境暴露方面存在差异，并评估大尺度风－海洋输运在塑造其遗传结构中的作用。
**位置**：热带和亚热带海域，即*Physalia*搁浅事件与分布范围重叠的区域。
**时间跨度**：基于2006年至2023年间的近期遗传采样数据及再分析资料，开展为期一年的漂移模拟。
**研究的主要类群**：*Physalia*（葡萄牙战舰/蓝瓶水母），一个由五种遗传上不同的物种组成的群体水螅纲（hydrozoan）属。
**方法**：研究人员将来自五个谱系的125个个体的遗传鉴定与数值拉格朗日粒子回溯相结合。漂移轨迹利用海洋和大气再分析场进行模拟，并采用一套实验推导的漂移参数化方案，该方案考虑了*Physalia*的非对称形态和风驱运动特征。由模拟输运得到的连通性度量指标与观测到的遗传分化进行了比较。
**结果**：漂移模拟揭示了谱系特异性的扩散和温度暴露差异。沿漂移路径的平均海表温度范围从*Physalia minuta*的15°C到*Physalia mikazuki*的28°C不等。谱系内的遗传分化与扩散连通性显著相关：地理距离与拉格朗日输运的结合可解释高达57%的遗传距离方差。
**主要结论**：风－海洋输运是驱动*Physalia*环境暴露和遗传结构的关键因素。扩散路径的差异使不同谱系暴露于不同的热力环境，这对生态学和生物地理学具有深远意义。随着气候变化改变风场、洋流和海洋温度，连通性和热力栖息地将以谱系特异性的方式发生改变，可能重塑*Physalia*物种的分布格局和演化轨迹。
**僧帽水母属（*Physalia*）海洋连通性、热力阈值与遗传结构关联研究解读**
**一、研究背景与科学问题**
海洋物种的遗传结构模式受到生活史特征与海洋学变率之间相互作用的影响。洋流、涡旋和锋面等海洋特征既可作为扩散屏障减少基因流从而隔离种群，又可作为长距离输运通道促进偏远地区间的扩散事件。此外，海洋环流还调节海洋温度，进而塑造热力栖息地并影响物种的生理耐受阈值。僧帽水母属（*Physalia*），俗称葡萄牙战舰或蓝瓶水母，是一种缺乏自主运动能力的群体水螅纲动物，几乎完全被动地漂浮于海气界面，主要依赖风、表层流和波浪驱动的输运进行扩散。这种缺乏强自主推进运动的特性为分离大尺度物理强迫在塑造跨洋盆遗传连通性中的作用提供了独特机遇。
*Physalia*频繁搁浅并与渔业和游泳者发生交互，在热带和亚热带地区常成群冲上岸滩，蜇伤游客、堵塞渔网并威胁养殖鱼类。然而，尽管其高可见度和显著影响，学界对其种群动态、栖息地偏好或开阔大洋分布知之甚少。目前缺乏其栖息地偏好信息——特别是其耐受和最适温度范围——使得预测*Physalia*如何响应变化中的海洋条件（如洋流和温度）变得困难。鉴于其扩散路径几乎完全由大气和海洋动力决定，全球风场和洋流模式的预期变化及已观测变化预计将强烈影响*Physalia*种群动态，要么增强先前隔离种群间的混合，要么将其输运至通常不栖息的区域。
*Physalia*的漂移和扩散形态高度复杂。每个个体由一个充满气体的"浮囊"维持浮力、一个容纳消化和生殖zooids的水下"膨体"以及用于捕食的长拖曳触手组成，这使其同时暴露于风和海洋强迫。此外，种群存在二态性，个体根据水下膨体相对于浮囊的位置而分为左旋型或右旋型，导致其与风向成一定角度漂移——分别偏向风的右侧或左侧。该漂移角并非固定不变，而是随风速变化：弱风时较大（约40°），强风时减小（约0°）。Church et al. (2025)基于遗传和形态分析对*Physalia*多样性进行了全球重新评估，鉴定出五个遗传上截然不同的物种：*Physalia physalis*（谱系A）、*Physalia utriculus*（谱系B1）、*Physalia mikazuki*（谱系B2）、*Physalia megalista*（谱系C1）和*Physalia minuta*（谱系C2）。全基因组数据的固定指数（FST）显示物种间遗传分化极高，范围达0.3至0.6。
本研究旨在弥合遗传结果与环境暴露和物理连通性之间的鸿沟。研究核心问题包括：个体在一年内能漂移多远？在此期间它们暴露于何种温度范围？通过回溯每个遗传采样个体的搁浅时间和位置，计算遗传研究中采样个体间的拉格朗日连通性，并评估漂移衍生的连通性模式是否与定义谱系内亚群的观测遗传结构一致。
**二、研究的技术方法**
研究人员利用欧洲中期天气预报中心（ECMWF）第五代再分析资料ERA5提供的10米高度逐时风场数据（水平分辨率0.25°）和Bluelink再分析BRAN2020提供的海洋环流数据（最浅层约2.5米，分辨率0.1°）作为环境强迫场。基于OpenDrift (v1.12.0)开源软件包开展拉格朗日粒子回溯实验，从Church et al. (2025)报告的125个遗传采样个体的搁浅位置和日期释放1000个虚拟粒子，回溯模拟一年。漂移参数化基于Bourg et al. (2024)利用3D打印*Physalia*复制品在水槽中开展的实验室实验结果，包含随风速变化的漂移角指数关系及平均为10米风速1.7%的漂移速度，并对已知旋向性样本赋予对应的漂移方向。遗传距离基于全基因组测序数据的主成分分析（PCA）计算，物种特异性PCA按谱系分别进行。统计分析采用普通最小二乘（OLS）回归模型，以配对遗传距离为响应变量，分别测试水上距离、拉格朗日相遇概率、拉格朗日最短距离及组合预测因子的解释力。
**三、主要研究结果**
**3.1 不同*Physalia*谱系在分布、海岸邻近度和温度暴露方面的差异**
基于Church et al. (2025)中*Physalia*样本搁浅位置回溯的一年拉格朗日轨迹显示，所有谱系的轨迹大体局限于其采集所在洋盆，但个体常穿越整个洋盆。*P. physalis*样本局限于北大西洋；*P. utriculus*分布最广，遍布太平洋、南大西洋和南印度洋；*P. mikazuki*仅见于北太平洋；*P. megalista*限于南半球，覆盖南印度洋和西南太平洋大部；*P. minuta*空间上最受限制，主要集中在西南太平洋。漂移能力方面，*Physalia*每年可旅行数千公里，但谱系间存在明显差异。180天后，*P. mikazuki*搁浅率最高（61%），其次为*P. minuta*（49%）；*P. physalis*（37%）、*P. megalista*（40%）和*P. utriculus*（26%）则保持更偏外海的轨迹。*P. physalis*平均位移最远（890公里），其后为*P. utriculus*（750公里）、*P. megalista*（730公里）、*P. minuta*（650公里）和*P. mikazuki*（420公里）。沿轨迹海表温度（SST）方面，各谱系平均暴露温度介于15°C至27°C之间：*P. physalis*中位数为22°C，*P. utriculus*为25°C，*P. mikazuki*为28°C，*P. megalista*为21°C，*P. minuta*为15°C。*P. mikazuki*限于最温暖条件（24°C–29°C），而*P. minuta*占据最冷范围（11°C–19°C）。子采样分析确认这些差异并非仅由采样量不均导致。2023年SST空间概率图显示，*P. mikazuki*和*P. minuta*的热力栖息地范围极为有限，且二者SST范围完全不重叠。
**3.2 海洋－风输运能否解释*Physalia*谱系内的遗传相似性**
对同属一个谱系的每对遗传样本，研究人员基于从每个样本位置释放的1000个粒子的一年回溯轨迹计算了拉格朗日相遇概率（LMP），并与样本采集位置间的水上距离（D）进行比较。所有谱系均表现为距离与连通性之间的负相关关系：相距小于1680公里的样本对均显示出非零相遇概率，而相距超过8570公里的样本对在一年模拟中从不连接。然而也存在"近而不连"和"远而相连"的例外。*P. physalis*的距离-连通性关系最清晰；*P. utriculus*变异性最大——相距不足1000公里的一些样本对几乎无相遇概率，而相距超过6000公里的其他样本对可达50%的相遇概率。将拉格朗日连通性结果与Church et al. (2025)的PCA遗传聚类比较显示：在*P. physalis*中，大西洋两侧（美洲沿岸与非洲-欧洲沿岸）的样本与各自大陆边缘显示出更高相遇概率，来自百慕大和亚速尔群岛的岛屿样本与最近大陆海岸及其彼此之间均显示强连通性。对于*P. utriculus*和*P. megalista*，因样本分布极远直接相遇概率极低，研究人员引入拉格朗日最短距离（LSD）指标——即回溯轨迹起点位置间的平均距离。*P. utriculus*中遗传相近的区域对（如日本-关岛、日本-东澳大利亚、东澳大利亚-北太平洋岛屿、南非-西澳大利亚）显示出最短的起点间距离，而遗传距离远的区域对（如西澳大利亚-北太平洋、南美洲-日本、南非-太平洋岛屿）对应最长的起点间距离。OLS回归模型显示，将所有谱系内样本对合并分析时，模型解释了27%至57%的遗传距离方差，所有p值低于0.05。组合预测因子D×LSD解释的方差（57%）高于单独使用D（42%）、LMP（27%）或LSD（51%）。按谱系分别分析时，*P. utriculus*和*P. minuta*的模型解释力最高（分别达55%和46%），最佳预测因子均为D×LSD组合；*P. physalis*的水上距离单独解释力与组合指标相同（均为29%）；*P. megalista*的最佳预测因子为LMP（解释26%方差）。
**四、讨论与结论**
本研究表明，所有谱系在模型中遵循相同的漂移行为，谱系间在扩散能力、海岸邻近度和温度暴露方面的显著差异反映了不同谱系所处环境的差异——即热带与高纬度分布、沿岸与岛屿样本、高动力与平静海域。*P. physalis*和*P. utriculus*展示出最高的位移速率、最低的搁浅率和最宽的SST暴露范围；*P. mikazuki*和*P. minuta*则漂移速度较慢、搁浅更频繁、SST暴露范围更窄。特别值得注意的是，*P. utriculus*和*P. mikazuki*作为遗传分化最低的两个谱系（Church et al. 2025），其SST暴露范围存在显著差异——*P. mikazuki*的SST下限（第10百分位数）几乎等于*P. utriculus*所遭遇的平均SST。SST空间概率图与iNaturalist公民科学平台观测记录的对比显示，*P. physalis*和*P. megalista*的纬度范围估计与观测吻合良好，*P. minuta*的观测点均落在预测的狭窄纬度带内，而*P. utriculus*在东南亚海域的缺失也与其SST范围估计一致。
物理输运与遗传结构的对比揭示了多重一致性。跨物种连通性分析显示，遗传相似度低的谱系（如*P. utriculus*、*P. megalista*和*P. minuta*）即使空间位置接近也表现出有限的拉格朗日连通性；而遗传分化最低的*P. utriculus*和*P. mikazuki*对在模拟中显示最高的长距离连通性——相距超过6000公里的部分样本对仍有高达50%的粒子轨迹重叠。谱系内模式亦部分反映Church et al. (2025)发现的空间遗传聚类：*P. physalis*中北大西洋两侧的样本在各自洋盆边缘内比跨洋盆更具连通性，岛屿样本（百慕大、亚速尔群岛）作为潜在连接枢纽；*P. utriculus*中遗传相近的区域往往具有较短的推断起点间距离。统计回归进一步支持拉格朗日连通性在塑造遗传分化中的作用——结合距离与连通性指标可将遗传距离方差的解释率提升至57%。
研究存在若干局限性：样本量在不同物种间存在差异且全球仅不足200个个体，采样偏向北半球；拉格朗日模拟未考虑波浪效应及风-浪-流的非线性相互作用；漂移参数化可能过度简化了形态多样性和行为变异性。尽管如此，本研究整合了目前关于*Physalia*遗传结构和漂移参数化的最佳可用知识。研究结论指出，不同谱系在热暴露和搁浅趋势方面存在明显差异，这些差异应推动对每个谱系生理耐受性的进一步研究。由于*Physalia*通过风和洋流被动扩散，其追踪适宜环境条件的能力可能有限。研究还发现拉格朗日追踪得出的物理连通性与Church et al. (2025)报告的遗传分化之间存在强烈相似性，这表明大尺度风-洋流输运在塑造*Physalia*种群结构中发挥着重要作用。随着全球风场和洋流模式的持续变化，种群间的连通性可能发生改变，从而影响*Physalia*物种的栖息地适宜性、分布模式和基因流，这构成未来研究的重要方向。
**结论**：风－海洋输运是驱动*Physalia*环境暴露和遗传结构的关键因素。不同谱系的扩散路径暴露于截然不同的热力环境，这对生态学和生物地理学具有深远意义。随着气候变化改变风场、洋流和海洋温度，连通性和热力栖息地将以谱系特异性的方式发生改变，可能重塑*Physalia*物种的分布格局和演化轨迹。
