Marine Ecology

Phytoplankton community structure and its relationship with environmental factors in the spring coastal region of Nan’ao based on morphology and high-throughput sequencing

  • ZHOU Zhixi , 1 ,
  • TANG Huijuan , 1 ,
  • KE Zhixin 2, 3 ,
  • LIU Jiaxing 2 ,
  • ZHOU Weihua 2
Expand
  • 1. South China Agricultural University, Guangzhou 510640, China
  • 2. Key Laboratory of Tropical Marine Bio-resources and Ecology, Chinese Academy of Sciences, Guangzhou 510301, China
  • 3. Comprehensive Observatory of Marine Ecosystems in the Upwelling Regions of Eastern Guangdong, Shantou 515041, China
TANG Huijuan. email:

Copy editor: SUN Cuici

Received date: 2024-02-28

  Revised date: 2024-03-25

  Online published: 2024-05-17

Supported by

Special Fund for Science and Technology Planning Project of Guangdong Province of China(2021B1212050023)

Guangdong Basic and Applied Basic Research Foundation(2022A1515010656)

Guangzhou Municipal Key Science and Technology Project(2023B03J1328)

Abstract

The structure of phytoplankton community was investigated based on morphological identification and high-throughput sequencing in the Nan’ao waters in May 2022. The spatial distribution of phytoplankton in relation to environmental factors was explored and results from two different methods were also compared and discussed. According to the cluster analysis of environmental factors, the survey area can be divided into nearshore area and offshore area. Totally, 105 species belong to 52 genera, 39 families, 27 orders, 8 orders, 5 phyla were morphological identified, while 543 species in 6 phyla, 32 orders, 97 families, 155 families, 272 genera through high-throughput sequencing were identified. Bacillariophyta and Dinophyta were the main groups, composing 74.54% and 24.78% through microscopic identification, and 17.52% and 67.84% through high-throughput sequencing, respectively. The abundance of phytoplankton identified by morphology ranged from 0.11×105 to 6.85×105 cells·L-1, which was significantly correlated with the distribution of chlorophyll a concentration, and the phytoplankton diversity index was lower than that of the high-throughput sequencing method. 7 dominant species were obtained by both methods, of which the absolutely dominant specie in the morphological identification results was Chaetoceros compressus, and the absolute dominant specie in the high-throughput sequencing identification results was Heterocapsa rotundata. Correlation analysis showed that the main environmental factor affecting phytoplankton community structure in Nan’ao in spring was pH, salinity and DIP (dissolved inorganic phosphorus). Using the combination of morphological identification and high-throughput sequencing technology, the structure and diversity of the phytoplankton community can be described more comprehensively and accurately.

Cite this article

ZHOU Zhixi , TANG Huijuan , KE Zhixin , LIU Jiaxing , ZHOU Weihua . Phytoplankton community structure and its relationship with environmental factors in the spring coastal region of Nan’ao based on morphology and high-throughput sequencing[J]. Journal of Tropical Oceanography, 2025 , 44(1) : 53 -65 . DOI: 10.11978/2024046

海洋中浮游植物作为初级生产者, 决定着生态系统的生物多样性, 为维持海洋生态系统的物质循环和能量传递提供了重要保障(Irigoien et al, 2004)。浮游植物分布范围广、生活周期短、对环境变化敏感(唐涛 等, 2002), 其群落结构组成和分布格局主要受环境因子调控, 可反映水质的状况变化, 是评价环境健康的重要指示生物之一(谭香 等, 2011)。虽然许多浮游植物的生长对生态系统健康是有益的, 然而其过量生长(如水华和赤潮)已成为世界范围内的一个主要环境问题。Hallegraeff 等(2021)发现, 赤潮事件的数量与水产养殖规模的扩大之间存在显著联系。众多研究显示赤潮的发生频率和分布范围将随着未来气候变化及人类活动影响而增加(Schläpfer et al, 1999; Barton et al, 2016; Gobler, 2020)。开展浮游植物调查工作对了解海域生态系统稳定性, 监测和防控赤潮藻种的暴发具有重要意义(Reynolds et al, 2002; 孙军 等, 2004)。
南澳岛位于广东省东部沿海, 总面积约为114.74km2, 是广东省唯一的海岛县级旅游城市, 也是中国重要的海洋牧场示范区之一, 养殖方式包括网箱养鱼、贝类底播和海藻养殖等(彭璇 等, 2014)。南澳海域位于南海北部, 是中国乃至世界海洋物种多样性最丰富的地区之一(Ma et al, 2008), 它受到韩江、榕江、莲阳河、黄冈河等径流以及台湾暖流、浙闽沿岸流、南海暖流等大洋环流共同影响(苏纪兰 等, 1990), 地理位置特殊, 生态环境比较脆弱(马华栋 等, 2021)。近几十年在人类活动的影响下, 本海域富营养化状态呈现出越来越严重的趋势(朱小山 等, 2005; 蔡德华 等, 2020; 黎素菊 等, 2022), 且出现过严重的赤潮现象, 多发种为球形棕囊藻(Phaeocystis globosa)(黄长江 等, 1999; 陈炜婷 等, 2023)。赤潮的爆发原因一般是解除了某种营养盐对浮游植物生长的限制作用, 这些营养盐的来源包括河流径流输入、沉积物的释放、上升流及沿岸流和海水养殖等(暨卫东 等, 1990; Lv et al, 2011; Xu et al, 2015)。研究显示, 在南澳周边海域浮游植物生长一般受到磷营养盐的限制(周凯 等, 2002; 陈丹婷 等, 2020; Yu et al, 2022)。
浮游植物群落结构常用于评价生态系统的健康程度。根据以往南澳岛附近海域基于形态学鉴定的调查, 浮游植物物种数量可达322种, 其中硅藻是主要类群(杜虹 等, 2011)。高通量测序技术相比传统方法更加便捷, 并且在低分类单元上有更强的分辨能力, 被广泛认为可以作为形态学检测结果的良好补充(Gilbert et al, 2008; Vandeputte et al, 2017)。鉴于部分浮游植物是微小的单细胞生物, 显微镜观察时难免会忽略 (Field et al, 1998; Esenkulova et al, 2020), 本次研究使用形态学鉴定方法结合高通量测序技术, 对南澳海域浮游植物群落结构及多样性进行分析, 并结合环境因子阐释浮游植物群落时空变化的驱动因素, 可为该海域环境保护和管理提供基础数据。

1 材料与方法

1.1 调查区域

本研究于2022年5月(春季)在汕头南澳海域进行样品采集, 共设置12个站位(图1), 其中S1站位于柘林湾网箱养鱼集中区, S2、S4站附近分布有密集的牡蛎和大型海藻的养殖, S9站周边是海上风电建设区。
图1 南澳海域地图及调查站位分布

该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1609号的标准地图制作, 底图无修改

Fig.1 Map of the Nan’ao coastal region and location of sampling stations

1.2 样品采集和处理

使用水质多参数仪(YSI-6600, 美国)现场记录各站位的水温(T)、盐度(S)、pH和浊度(turbidity, Tur); 表层水样采用有机玻璃采水器从水表以下50cm处采集。取表层海水1L装入样品瓶, 立即加入1%体积的鲁哥试剂固定, 再将样品自然沉降24h后浓缩至50mL, 用于浮游植物形态学鉴定与计数。取1mL浓缩后的样品, 放入Sedgwick计数框内在显微镜下以200倍或400倍下进行形态鉴定和计数, 浮游植物形态学分类方法参照郭玉洁(2003)和林更铭 等(2021)专业书籍。将0.5~2L水样通过0.2μm滤膜过滤, 滤膜在液氮中快速冷冻, 然后保存于-20℃的冰箱, 以用于浮游植物的高通量测序分析(Jiang et al, 2015)。
采取水样0.25~1L, 用孔径为0.7µm的玻璃纤维滤膜(GF/F)过滤。过滤后的滤膜用中性滤纸吸干水分后再用锡箔纸包裹, 将滤膜放入-20℃冰箱保存用于叶绿素a浓度测定, 测定方法为萃取荧光法。过滤后的滤液放入-20℃冰箱保存用于营养盐分析, 溶解态无机营养盐[硝酸盐NO- 3 -N、亚硝酸盐NO- 2-N、氨氮NH+ 4-N、溶解态无机磷酸盐 (dissolved inorganic phosphorus, DIP)、溶解态无机硅酸盐 (dissolved inorganic silicon, DSi)]使用AA3连续流动分析仪(Seal Analytical, 德国)测定; 溶解态无机氮盐(dissolved inorganic nitrogen, DIN)为NO- 3 -N、NO- 2-N和NH+ 4-N的总和。以上操作参照《海洋调查规范 第6部分:海洋生物调查》(GB/T12763.6-2007) (中华人民共和国国家质量监督检验检疫总局 等, 2007)和《近岸海域环境监测规范》(HJ 442.3-2020)(中华人民共和国生态环境部, 2020)执行。

1.3 浮游植物高通量测序

参考刘卫东 等(2017), 采用CTAB法提取样品的宏基因组。对于真核浮游植物18S rDNA, 使用通用V4区引物V4-F(5ʹ-CCAGCASCYGCGGTAATTCC-3ʹ)和V4-R(5ʹ-ACTTTCGTTCTTGATYRATGA-3ʹ)进行扩增。PCR扩增采用TransStart Fastpfu DNA Polymerse (TransGen AP221-02), 反应体系为50μL, 扩增反应在ABI Gene Amp 9700型PCR仪上进行, 反应条件为: 94℃预变形3min; 94℃变形30s, 46℃下退火复性30s, 72℃延伸45s, 共30个循环; 最后于72℃下再延伸5min。每个样本按照上述步骤进行3个重复, 将同一样本的PCR产物混合后用2%琼脂糖凝胶电泳监测, 使用AxyPrepDNA凝胶回收试剂盒(AXYGEN公司)切胶回收PCR产物, 2%琼脂糖电泳检测。将检测合格的产物交给上海凌恩生物公司测序, 利用Illumina PE250平台进行文库构建, 使用Hiseq2500 PE250进行测序。

1.4 数据分析

使用mothur(version 1.30.1)软件对浮游植物物种多样性指数进行分析, 多样性指数有: Shannon-Wiener指数( H ')、Simpson指数(D)、Pielou指数( J ')。采用优势度指数(Y)描述浮游植物群落优势种, 以Y≥0.02确定为优势种。高通量测序结果计算时以质控