盐度与混合营养型驱动长江口有害藻华物种分布的关键机制

《Diversity》：Salinity and Mixotrophy as Key Drivers of HAB Species Distributions in the Changjiang Estuary

　　

长江口是中国有害藻华（HABs）的主要热点区域。研究人员于2025年3月、8月、10月和2026年1月四个季节航次，在14个站位开展宏条形码（metabarcoding）分析，揭示该区域有害藻华物种的组成、动态及环境驱动因子。共检出63种有害藻华物种，其中硅藻全年占优势，但甲藻从春季（21.6%）至夏秋季（34.8–44.6%）显著增加。盐度成为塑造空间格局的首要因子：广盐性硅藻中肋骨条藻（Skeletonema costatum）主导低盐近岸水域，而多种甲藻包括卡尔藻（Karlodinium veneficum）、微小原甲藻（Prorocentrum micans）、圆形异甲藻（Heterocapsa rotundata）和链状亚历山大藻（Alexandrium catenella）则在离岸海域占优。季节尺度上，硅藻与甲藻的比例从春季的10.6下降至秋季的1.6，这与温度升高和氮磷比（N:P ratio）降低同步发生。值得注意的是，大多数高丰度甲藻属于混合营养型（mixotrophic），这一特征很可能赋予其竞争优势。研究结果表明，有害藻华物种的动态受到盐度驱动的栖息地过滤与营养功能性状之间相互作用的塑造，这为将混合营养型作为关键指标纳入未来有害藻华监测与预警框架提供了依据。未来的监测工作应将这些分子趋势与传统定量方法相结合进行验证。

盐度与混合营养型驱动长江口有害藻华物种分布的关键机制——基于四季度宏条形码调查的研究解读

一、研究背景与科学问题

长江口是中国近海有害藻华（HABs）危害最严重的区域之一。过去五十年来，受人为富营养化、气候变暖和水文变化驱动，全球范围内有害藻华的频率、强度和地理范围显著增加。长江口作为连接中国大陆最大河流与西太平洋的重要过渡带，接纳了大量农业、工业和生活源营养盐输入，为浮游植物水华创造了极为有利的条件。历史上，该区域曾多次发生由鱼毒性甲藻引发的水产养殖重大经济损失事件，以及贝类毒素导致的人体中毒事件。

然而，当前对长江口有害藻华的研究存在若干关键空白：已有的宏条形码研究多限于单季采样，无法捕捉完整的季节变异；有害藻华物种的全年动态、其季节性高丰度窗口期、触发增殖的环境因子，以及确定性环境驱动与随机过程在群落组装中的相对重要性均缺乏系统评估。此外，传统形态学鉴定方法劳动强度大、易出错，难以区分形态相似但遗传不同的隐存种，小型浮游植物容易被遗漏或误判，导致对有害藻华物种丰富度的系统性低估。

因此，研究人员开展了一项四季宏条形码调查，旨在：（1）全面记录长江口各季节浮游植物和有害藻华物种多样性；（2）刻画优势类群的季节演替模式；（3）解析有害藻华物种沿近岸—离岸梯度的空间分布格局及其与环境驱动因子的关系。

二、研究方法概述

研究人员于2025年3月（春季）、8月（夏季）、10月（秋季）和2026年1月（冬季）在长江口14个站位采集表层海水样品。采用CTD剖面仪记录原位温度和盐度，按国家标准方法测定硝酸盐、亚硝酸盐、氨氮、磷酸盐和硅酸盐浓度。提取滤膜样品总基因组DNA后，扩增真核生物18S rRNA基因V4区，在Illumina NovaSeq 6000平台进行双端测序。原始数据经质量过滤、引物去除、拼接和嵌合体剔除后，利用DADA2流程降噪生成扩增子序列变异体（ASVs），以PR2数据库进行分类注释（≥99%相似度定为种级）。基于国际海洋考察理事会—联合国教科文组织政府间海洋学委员会（IOC-UNESCO）有害藻华物种参考列表及长江口历史记录确定有害藻华物种。所有统计分析在R软件中进行，包括α多样性指数计算、基于Bray-Curtis相异性的非度量多维尺度分析（NMDS）和置换多元方差分析（PERMANOVA），以及Spearman相关性分析。

三、研究结果

3.1 浮游植物组成与相对丰度

测序获得2,461,747条原始读长，经处理后得到12,312个ASVs，其中5,239个（42.6%）归属浮游植物。这些ASV主要隶属于6个门，以褐藻门（Ochrophyta，51.0%）和甲藻门（Dinophyta，33.5%）最为丰富。属水平上，中肋骨条藻属（Skeletonema）丰度最高（328个ASVs），其次为小环藻属（Cyclotella，260个ASVs）。

3.2 浮游植物群落的组成动态

α多样性呈现显著的季节和空间变异：秋季平均丰富度最高（313），夏季最低（162），但冬季Shannon指数较高（3.38）源于高均匀度。空间上，夏季α多样性从近岸向离岸递减，秋季则在近岸和离岸均保持极高值。三个主要门类（褐藻门、甲藻门和隐藻门Cryptophyta）占总相对丰度的95%以上。硅藻（Coscinodiscophyceae）在夏季达到峰值（48.8%），冬季骤降至9.8%，伴随隐藻（Cryptophyceae）激增至26.9%。属水平上，中肋骨条藻属（Skeletonema）从春季到秋季占优（25.7–28.7%），但在冬季被隐藻属（Teleaulax，18.0%）和甲藻鳞甲藻属（Lepidodinium，7.9%）取代。NMDS分析证实季节解释了33.1%的群落变异。

3.3 有害藻华物种的季节与空间分布

共检出63种已知有害藻华物种。春季以中肋骨条藻（Skeletonema costatum）为主，在近岸站位占比14.1–36.8%；圆形异甲藻（Heterocapsa rotundata）随离岸距离增加而上升，在站位12达30.4%。夏季有害藻华丰度最高，中肋骨条藻在离岸站位14占92.5%，卡尔藻（Karlodinium veneficum）在多个站位达26.9–34.2%。秋季物种丰富度最高（48种），分布最均匀。冬季中肋骨条藻大幅下降（多数站位0.2–1.9%），而微小原甲藻（Prorocentrum micans）在站位3跃升至27.7%。

3.4 环境因子季节变化及其与有害藻华物种的关联

水温从春季9.69°C升至夏季29.45°C，冬季降至7.97°C。盐度冬季最高（15.35‰），其他季节低于9.80。硝酸盐和硅酸盐是最丰富的无机营养盐，秋季达峰值。氮磷比（N:P ratio）范围为64.76–118.30，春季偏高、冬季偏低。相关性分析显示：中肋骨条藻与温度正相关、与盐度负相关；卡尔藻与盐度最强正相关（r=0.79）；微小原甲藻与温度强负相关（r=?0.85）、与氨氮正相关（r=0.74）；圆形异甲藻与盐度和氨氮均正相关，但与磷酸盐和硅酸盐强负相关。

四、讨论与结论

讨论部分总结

研究人员指出，硅藻全年占据优势，但甲藻比例呈现显著季节和空间变异：硅藻与甲藻比例从春季10.6降至秋季1.6，伴随温度升高和氮磷比降低，符合经典框架中硅藻在凉爽、营养充足条件下占优而甲藻在温暖、营养受限条件下被偏好的规律。与历史记录的125种相比，本研究仅检出63种，差异可能源于许多历史记录种的稀有性或瞬时性。空间上，有害藻华物种的丰富度和相对丰度一致表现为离岸高于近岸，这与盐度梯度主导的生态格局一致，同时黑潮暖流可能在夏季为离岸甲藻提供外来种群来源。

盐度作为首要空间过滤器，划分出两类响应群体：低盐偏好组以中肋骨条藻为代表，与盐度负相关；高盐偏好组涵盖多种甲藻，如卡尔藻、微小原甲藻、链状亚历山大藻等。在此基础上，不同形式的无机氮进一步塑造了物种水平的生态位分化。尤为重要的是，高丰度甲藻普遍具有混合营养型（mixotrophic）特征，这种营养灵活性使其在长江口长期氮磷失衡背景下可能发挥越来越重要的作用，且多种产毒甲藻在磷限制下毒素产量升高，凸显了混合营养型有害藻华物种带来的更高生态和经济风险。

结论部分翻译

本研究通过四个季节航次的宏条形码调查，研究了长江口有害藻华物种的动态。结果显示，群落组成在空间和季节上均有显著变化：盐度充当主要的物理过滤器，而季节变暖和氮磷比下降驱动了从硅藻到甲藻的演替。一个关键发现是混合营养型在高丰度甲藻中的普遍存在，这很可能赋予其相对于专性自养者的竞争优势。这种营养可塑性，加上多种甲藻类群的产毒潜力，构成了常规监测可能低估的更高风险。然而，宏条形码无法提供定量数据，读长比例可能因引物结合和扩增效率差异以及参考数据库不完整而产生偏差。这些发现应通过将宏条形码与显微镜检查、细胞计数、叶绿素a和色素分析以及扩展的参考数据库相结合来进一步验证。