海洋与湖沼  2022, Vol. 53 Issue (6): 1441-1454   PDF    
http://dx.doi.org/10.11693/hyhz20220400113
中国海洋湖沼学会主办。
0

文章信息

李雯璐, 徐大鹏, 季道德. 2022.
LI Wen-Lu, XU Da-Peng, JI Dao-De. 2022.
莱州湾近岸海域浮游纤毛虫群落组成和季节变化
COMMUNITY COMPOSITION AND SEASONAL VARIATION OF PLANKTONIC CILIATES IN SANSHAN ISLAND, LAIZHOU BAY
海洋与湖沼, 53(6): 1441-1454
Oceanologia et Limnologia Sinica, 53(6): 1441-1454.
http://dx.doi.org/10.11693/hyhz20220400113

文章历史

收稿日期:2022-04-28
收修改稿日期:2022-05-12
莱州湾近岸海域浮游纤毛虫群落组成和季节变化
李雯璐1, 徐大鹏2, 季道德1     
1. 烟台大学海洋学院 山东烟台 264005;
2. 厦门大学近海海洋环境与科学国家重点实验室 海洋与地球学院 海洋微型生物与地球圈层研究所 福建厦门 361104
摘要:浮游纤毛虫是真核微型生物的重要类群, 在浮游生态系统中发挥着重要作用。利用传统形态学与高通量测序技术相结合的方法, 于2020年9月至2021年8月对莱州湾三山岛近岸浮游纤毛虫进行每月两次的调查。对环境DNA的18S rRNA基因的V4区进行高通量测序, 结果显示无壳寡毛类(占纤毛虫序列丰度的40.14%)和砂壳纤毛虫(10.67%)是优势类群, 并且具有明显的季节变化。无壳寡毛类在冬季序列相对丰度最高, 主要由急游虫属(Strombidium)和急游虫科未定类群(Strombidiidae_X)组成; 砂壳类在春夏季序列相对丰度最高, 主要由砂壳目未定类群(Tintinnida_XX)和拟铃虫属(Tintinnopsis)组成。浮游纤毛虫群落的β多样性呈现显著的季节变化, 而α和β多样性在站点间差异均不显著。与环境因子的RDA分析显示, 总有机碳和溶解氧含量是影响纤毛虫群落结构的主要环境驱动因子。形态学方法共检获鉴定17种纤毛虫, 分别隶属于7个属, 其中砂壳类在多样性和检获频率上均为最优势类群。研究填补了莱州湾三山岛海域浮游纤毛虫的多样性、群落结构和季节动态研究的空白, 同时为深入研究温带海湾近岸浮游纤毛虫群落的长期变化规律以及对环境因子的响应机制提供本底资料。
关键词浮游纤毛虫    多样性    近岸海域    高通量测序    形态学    
COMMUNITY COMPOSITION AND SEASONAL VARIATION OF PLANKTONIC CILIATES IN SANSHAN ISLAND, LAIZHOU BAY
LI Wen-Lu1, XU Da-Peng2, JI Dao-De1     
1. School of Ocean, Yantai University, Yantai 264005, China;
2. State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Science, Institute of Marine Microbes and Ecospheres, Xiamen University, Xiamen 361104, China
Abstract: Planktonic ciliates are important eukaryotic microorganisms and play a crucial role in the planktonic ecosystem. Semimonthly ecological research was conducted on planktonic ciliates in Sanshan Island, Laizhou Bay from September 2020 to August 2021 combining traditional morphology and high-throughput sequencing technology (HTS). HTS for the V4 region of the 18S rRNA gene showed that aloricate oligotrichs (accounting for 40.14% of the abundance of ciliate sequences) and tintinnids (10.67%) were dominant groups with obvious seasonal variation. The relative abundance of oligotrichs consisted of mainly Strombidium and Strombidiidae_X, and was the highest in winter. The relative abundance of tintinnids was the highest in spring and summer, and they were composed of mainly Tintinnida_XX and Tintinnopsis. Beta diversity showed significant seasonal variation, but both alpha diversity and beta diversity showed no considerable difference in spatial distribution. Redundancy analysis demonstrated that total organic carbon and dissolved oxygen content were the significant environmental factors on the seasonal variation of planktonic ciliates in the coastal water of Sanshan Island. Seventeen species were identified in microscopy method, and they belong to 7 genera. Tintinnids were the most dominant group in diversity and occurrence. This work filled the gaps in species diversity, community structure, and seasonal dynamic of planktonic ciliates in Sanshan Island, Laizhou Bay. Meanwhile, this study provided essential data for future studies on the long-term variation of planktonic ciliates community in the coastal waters of temperate bays and their response mechanisms to environmental factors.
Key words: planktonic ciliates    diversity    coastal water    High-throughput sequencing    morphology    

浮游纤毛虫是一类广泛分布于海洋中上层、营浮游生活的纤毛虫, 粒级大小一般在5~200 μm之间。它们是pico-级(0.2~2 μm)和nano-级(2~20 μm)浮游生物的主要消费者(Suzuki et al, 2007), 同时也是后生动物和鱼类幼虫的食物来源(Stoecker et al, 1987; Gómez, 2007), 在连接微食物网和传统食物链的物质循环和能量流动中发挥着重要作用(Pierce et al, 1992; Weisse et al, 2022)。

海洋浮游纤毛虫主要包括砂壳纤毛虫和无壳寡毛类, 也包括一些前口类、盾纤类等类群。依照Lynn (2008)的分类系统, 砂壳类纤毛虫隶属于旋毛纲(Spirotrichea)、环毛亚纲(Choreotrichia)、砂壳目(Tintinnida), 其重要特征是具有管状或瓶状的外壳, 因其近岸丰度高、采样易获得的特点, 目前世界各海区已经鉴定约1 000种(Dolan, 2010)。壳的口径和长度、壳的形状、黏着颗粒的有无是砂壳类重要的分类学依据(Capriulo et al, 1982; Dolan et al, 2002)。无壳寡毛类纤毛虫隶属于旋毛纲(Spirotrichea)、寡毛亚纲(Oligotrichia)、寡毛目(Oligotrichida)。目前已报道无壳寡毛类约200余种(Liu et al, 2017), 然而据Foissner等人(2008)Agatha (2011)估计, 无壳寡毛类在世界范围内约有830~1 280种, 其多样性可能有83%~89%以上是未知的。

传统的形态学由于分辨率低、费时费力等弊端, 难以避免地对形态上不可区分的物种(即隐性多样性)产生遗漏。高通量测序(high-throughput sequencing, HTS)因其通量高、耗时短、相对成本低等特点, 对微型生物生态学研究产生了革命性影响。对环境DNA的高通量测序可以揭示纤毛虫的多样性、分布模式和长时间尺度上的动态变化, 从而获得更为全面的群落信息(Slapeta et al, 2006; De Schepper et al, 2019)。然而研究显示基于HTS得出的纤毛虫多样性大大超过已知形态物种数量(Bachy et al, 2013), 甚至高达6倍之多(De Vargas et al, 2017), 因此基于形态和分子的分布模式是否一致尚未可知。近年来, 越来越多的研究者尝试将多种方法相结合对纤毛虫进行生态调查: Santoferrara等人(2016)对Lugol’s固定样品镜检的同时结合HTS, 对大西洋砂壳类纤毛虫的群落组成和分布模式进行研究。Zhao等(2017)利用HTS比较了环境DNA和RNA两种方法在纤毛虫群落垂直分布上的差异, 显示出RNA获得的多样性高于DNA。在调查南海浮游纤毛虫多样性和垂直分布时, Sun等(2019)将HTS与定量蛋白银染色(QPS)相结合, 结果显示两种方法在揭示较低分类等级的群落组成时存在差异。另外, 有研究者将QPS的形态学方法与DNA及其cDNA的HTS三者相结合, 调查沉积物中纤毛虫的多样性(黄平平等, 2017)。综合来看, 利用分子测序与传统形态学方法所揭示的群落组成及其变化趋势较为一致, 尤其是常见/优势物种, 而高通量测序因其测序深度高可发现大量稀有物种序列。因此将两种方法相结合能够在一定程度上弥补两者的弊端, 从而获得更为全面的纤毛虫群落信息。

近岸环境与人类社会的关系最为密切, 且同时受到人为因素和自然因素的双重影响而具有一定的特殊性。一方面, 三山岛近岸海域受到多重人类生产活动的影响, 海水富营养化严重、水质下降; 另外, 受到风暴潮等自然因素的影响, 水体扰动频繁使得浑浊度较高。这些因素可能对该生境中微食物网成员(如纤毛虫)的丰度、物种多样性和群落结构产生影响。本研究首次利用传统形态学和分子测序技术相结合的手段, 对莱州湾三山岛近岸海域的浮游纤毛虫进行周年生态调查, 旨在了解该海域的浮游纤毛虫的多样性、群落组成和时空变化, 以及对环境因子的响应。

1 材料与方法 1.1 调查站位、样品采集和预处理

本研究于2020年9月至2021年8月在莱州湾三山岛近岸海域设置3个采样点, A、B、C站点分别位于码头、离岸300 m处和养殖厂排水口(图 1), 使用渔船进行每月两次的浮游样品采集。由于海况原因, 缺失3个站点的1、2、7、8月上旬和5月下旬的样品及数据, 以及离岸的B站位11、1和4月下旬样品及数据。因此, 秋季(9~11月)共采集样品17份, 冬季(12~2月)样品11份, 春季(3~5月)样品14份, 夏季(6~8月)样品12份, 共计54个样品。

图 1 莱州湾三山岛近岸采样点示意图 Fig. 1 Locations of the sampling sites in Sanshan Island, Laizhou Bay

使用便携式水质分析仪(YSI, USA)测量原位水温、盐度、溶解氧和pH值。使用卡扣式采水器采集表层水样(深度20 cm), 每一个站位采集两个水样各1 L, 混合为一个2 L的样品。用孔径为200 μm的筛绢进行预过滤, 再使用聚碳酸酯滤膜(Millipore, Ireland) (直径45 mm、孔径0.22 μm)对1 500 mL水样进行抽滤收集, 设置三组重复(每张膜500 mL)。将抽滤得到的含有环境DNA的滤膜保存于冻存管中, 立即置于−80 ℃冰箱中保存。取100 mL表层原位海水, 经滤膜孔径为0.45 μm的针头过滤器过滤后收集在聚乙烯瓶中, 并置于−20 ℃冰箱中保存, 用于营养盐分析。氨态氮(NH4+), 硝态氮(NO3-), 亚硝态氮(N)2-), 磷酸盐(PO43-)和硅酸盐含量(SiO32-)使用营养盐自动分析仪(Seal, Germany)进行分析。水样的总碳(TC)和总无机碳(TIC)含量使用Vario微型立方元素分析仪(Elementar, Germany)进行测定。使用400目浮游生物网(孔径37 μm)在A、B两个站点水平拖网采集镜检样品并用广口采样瓶收集, 获得35个样品用于纤毛虫的形态学分类鉴定。

1.2 高通量测序方法

环境总DNA使用Fast DNA SPIN Kit for Soil DNA快速提取试剂盒(Mpbio, USA) 进行提取。使用NanoDrop One (Thermo-Fisher, USA)对DNA的浓度和纯度进行测定后, 于−80 ℃保存。

利用真核生物通用引物528F (5′-GCGGTAATT CCAGCTCCAA-3′)和706R (5′-AATCCRAGAATTT CACCTCT-3′)(Cheung et al, 2010) 对18S rRNA基因的V4区进行PCR扩增, 于北京诺禾致源生物信息科技有限公司利用Illumina NovaSeq测序平台对文库进行双末端PE250测序。

1.3 数据处理及分析

从下机数据中拆分出各样本数据, 截去barcode和引物序列后使用FLASH软件(V 1.2.11, http://ccb.jhu.edu/software/FLASH/)(Magoč et al, 2011)对样本的测序结果进行拼接, 得到低质量原始测序数据。随后使用fastp软件对得到的低质量原始测序数据进行质控, 得到高质量数据。利用QIIME和Mothur软件包(Schloss et al, 2009; Caporaso et al, 2010)去除嵌合体, 利用UCLUST以97%的相似性水平划分可操作分类单元(operational taxonomic units, OTUs), 去除只有一条代表序列的OTU, 利用原生生物核糖体参考数据库(Protist Ribosomal Reference, PR2)数据库进行物种分类注释, 得到最终的OTU表格。抽取纤毛门数据用于纤毛虫群落的分析, 将各样品的序列相对丰度数据按照季节划分整合。

为了保证样品间α多样性比较的可靠性和准确性, 以所有样品中最小序列数为标准, 对样品数据进行标准化, 然后计算α多样性指数(Chao 1、Richness、Shannon和Simpson指数)。通过单因素方差分析(ANOVA)对浮游原生生物α多样性的时空差异进行统计检验。

基于Bray-Curtis距离矩阵进行浮游原生生物群落的β多样性的分析。使用PRIMER 6.0软件对浮游原生生物群落进行非计量多维尺度分析(non-metric multidimensional scaling, NMDS); 利用R中的ggplot2包进行主坐标分析(principal co-ordinates analysis, PCoA), 并利用ANOSIM分析检验三个站点不同季节的浮游原生生物群落结构的差异情况; 使用SPSS软件对浮游原生生物主要类群和环境因子进行Spearman相关分析; 使用R中的vegan包进行冗余分析(redundancy analysis, RDA), 寻找影响浮游原生生物群落结构的主要环境因子。

1.4 形态学方法

将拖网样品充分混匀后利用解剖镜进行镜检观察, 依据文献对浮游纤毛虫进行鉴定和计数(Kofoid et al, 1929; 宋微波等, 2009; 张武昌等, 2012), 并拍照记录生物学特征及出现频次, 以获得浮游纤毛虫多样性信息。将镜检鉴定后的纤毛虫用灭菌水清洗数次后转移至灭菌后的1.5 mL离心管中, 用于提取DNA。引物为: EukA (5′-AACCTGGTTGATCCTGCCA GT-3′)、EukB (5′-TGATCCTTCTGCAGGTTCACCT AC-3′)、Euk82F (5′-GAAACTGCGAATGGCTC-3′)和U1492R (5′-GGTTACCTTGTTACGACTT-3′)(Medlin et al, 1988; López-García et al, 2001; Doty et al, 2005)。PCR扩增程序依照(Wang et al, 2019)进行。使用TIANgel Midi Purification Kit (天根生化科技有限公司, 上海) 试剂盒纯化PCR产物。将纯化后的扩增产物插入质粒载体pTZ57R/T (InsTAcloneTM PCR Cloning Kit, Fermentas, Canada), 并将重组质粒转移至感受态细胞Trans1-T1 (全式金生物技术股份有限公司, 中国北京)中。摇菌培养1 h后接种至琼脂糖培养基上过夜培养。随机挑选6个重组克隆, 利用引物M13-20和M13-26对菌样进行PCR验证。将满足条件的菌液送至青岛生工测序公司进行测序。将测序后的DNA序列与数据库中的序列进行比对, 相似度99%以上则认为是同一物种。

2 结果 2.1 环境因子

12个环境因子(包括温度、盐度、溶解氧、pH、总有机碳、总无机碳、总碳、磷酸盐、氨态氮、亚硝态氮、硝态氮、硅酸盐含量)时空变化如图 2所示。三个站点间的所有环境因子均无显著差异, 有6个环境因子表现出明显的季节模式: 水温(ANOVA, P < 0.05)、总碳(ANOVA, P < 0.01)、总无机碳(ANOVA, P < 0.01)和总有机碳含量(ANOVA, P < 0.01)一致呈现为夏秋季 > 冬春季的趋势; 盐度(ANOVA, P < 0.01)和溶解氧(ANOVA, P < 0.05)呈现冬春季 > 夏秋季的趋势。三个站点间盐度和有机碳含量的季节变化折线图变化趋势相似。

图 2 不同站点与季节的环境因子 Fig. 2 Environmental factors of different sites and seasons. 注: Ps为季节对环境因素的单因素方差分析(ANOVA)的显著性, P < 0.05表示差异显著; Pr是站点对环境因素的显著性, P < 0.05表示差异显著。

A站点水温年变化为1.7~27.9℃, 平均15.4 ℃, 其中最高温出现在2020年9月上旬; 水温最低出现在2021年1月下旬; B站点水温年变化为4~27 ℃, 平均14 ℃, 其中最高温出现在2021年8月下旬; 水温最低出现在2020年12月和2021年2月下旬; C站点水温年变化为10.9~30.4 ℃, 平均22.2 ℃, 其中最高温出现在2021年7月下旬, 水温最低出现在2021年2月下旬。三个站点的溶解氧在夏秋两季含量较为一致且最低值出现在夏季, 而冬季A站点溶解氧含量最高(27.14 mg/L), 其次是B站点(21.52 mg/L)。A、B两个站点的pH值在秋冬两季的变化幅度较大。A、B站点的磷酸盐含量在春夏两季较高, 氨态氮和亚硝态氮含量在春夏季交替时升高, 硝态氮则相反, 在春夏季交替时含量降低, 硅酸盐含量在夏秋两季最高, 而在春季达到最低值。C站点的营养盐含量季节变化较大, 氨态氮、硝态氮和亚硝态氮含量均呈现出夏季最低, 秋冬季高的特点, 最低值出现在7月上旬和8月下旬, 最高值出现在9月上旬和12月下旬。磷酸盐含量呈现夏季最低, 春秋季最高, 最低值出现在6月下旬, 最高值出现在4月上旬。

在三个站点间, A、B站点环境因子的全年变化拟合程度较高, 但pH在秋冬季变化趋势相反, 磷酸盐浓度在春夏季呈现明显不同的变化。总的来看, C站点海水理化因子变化较为平缓而营养盐浓度变化浮动较大, 与A、B站点变化拟合程度较差。

2.2 群落组成与相对丰度

高通量测序结果表明, 旋毛纲(Spirotrichea)是最优势类群, 占纤毛虫序列的64.10%, 其次是寡膜纲(Oligohymenophorea, 2.17%) (图 3a)。将旋毛纲(Spirotrichea)在目(order)水平上进一步细分如图 3b所示, 无壳寡毛类(急游虫目和环毛目共占序列总丰度的40.14%)和砂壳类(砂壳目, 10.67%)是浮游纤毛虫的主要优势类群, 并且显示出明显的季节模式。其中急游虫目(Strombidiida)序列的相对丰度为27.27%, 环毛目(Choreotrichida)为12.88%。急游虫目在全年的变化较为稳定并且具有明显的季节变化模式和站点差异。季节变化表现为秋冬 > 春夏, 且在春夏季样品中呈现A > B > C, 在秋冬季样品B > A > C。环毛目在站点间无明显差异, 而显示出秋 > 夏 > 冬 > 春的季节模式, 春季相对丰度明显降低。砂壳目在站点间差异显示出A(4.34%) > B(4.10%) > C(2.23%), 且A和B站点砂壳目在夏秋两季的相对丰度相近, 而在冬春两季的相对丰度显示出相反的变化模式。

图 3 浮游纤毛虫在纲(a)、目(b)阶元上的相对丰度, 以及主要优势类群在科(c)、属(d)阶元上的相对丰度 Fig. 3 Relative sequence abundance of planktonic ciliates at class (a), order (b), family (c), and genus (d) level

寡毛类和砂壳类在科(family)级水平上的群落组成如图 3c所示, 寡毛类中相对丰度较高的有Strombidiidae_M (16.62%)、Strombidiidae (14.63%)、Strombidiidae_R (9.77%)和Strobilidiidae_G (8.00%); 砂壳类中相对丰度较高的是Tintinnida_X (7.59%)、Favellidae (2.53%)和TIN_07 (2.32%)。

在属阶元上对寡毛类和砂壳类进行进一步划分如图 3d, 寡毛类的相对丰度主要由急游虫属(Strombidium, 19.01%)和急游虫科未定类群(Strombidiidae_X, 12.24%)组成, 相对丰度显示出先下降后上升的趋势, 并在冬季达到最大值。砂壳类主要由砂壳目未定类群(Tintinnida_XX, 3.84%)、拟铃虫属(Tintinnopsis, 3.27%)和网纹虫属(Favella, 1.32%)组成, 其相对丰度在春夏季达到最大值后下降, 冬季的相对丰度最低。

站点间的物种组成存在明显差异, 值得注意的是无论处于哪个分类阶元, C站点的夏季样品在物种组成和相对丰度上均与其他分组样品有明显的不同, 主要表现在总的物种多样性低, 以及除寡毛类和砂壳类外的纤毛虫相对丰度与其他样品有明显差别。在纲水平上, 旋毛纲相对丰度明显降低, 而叶咽纲(Phyllopharyngea)、叶口纲(Litostomatea)、寡膜纲(Oligohymenophore)具有明显高的相对丰度。在旋毛纲中, 管口目(Cyrtophoria)、前口目(Haptoria)、缘毛目(Peritrichia)的相对丰度明显高于其他站点, 而CONThreeP_X、Tintinnida和Strombidiida明显低于其他站点。A、B站点物种组成的季节变化较为相似(图 3b), 不同之处主要在于A站点的寡毛类纤毛虫相对丰度在春季明显大于B站点, 随着季节推移二者差距减小, 直至秋季相对丰度相当。冬季则呈现出B站点的急游虫目纤毛虫相对丰度大于A站点。而砂壳类纤毛虫的相对丰度变化趋势相反, 在春季B站点大于A站点, 随后差距缩小, 直至冬季A站点砂壳类纤毛虫的相对丰度大于B站点。

2.3 α多样性

计算54个样品的α多样性指数, 并进行站点和季节间的ANOVA方差分析检验, 结果如表 1所示, 浮游纤毛虫群落的丰富度和多样性指数均未显示出站点间和季节间差异(P > 0.05)。Chao 1指数显示出秋季 > 春季 > 冬季 > 夏季的季节分布, 而除Chao 1指数外的α多样性指数均显示出夏季 > 秋季 > 春季 > 冬季, 其中春秋两季的α多样性指数均最为接近。

表 1 各季节和各站点的浮游纤毛虫α-多样性指数 Tab. 1 α-diversity of planktonic ciliates in different samples
Chao 1 Richness Shannon Simpson
A站点 29.22±9.25 13.77±2.99 3.30±0.58 0.84±0.09
B站点 32.73±11.71 14.39±3.32 3.41±0.61 0.86±0.10
C站点 33.49±10.40 14.57±2.74 3.45±0.48 0.86±0.07
F 1.26 1.04 0.92 0.87
P 0.30 0.38 0.44 0.46
春季 33.71±9.43 14.73±2.48 3.46±0.44 0.87±0.06
夏季 27.99±13.10 13.38±3.43 3.26±0.57 0.84±0.08
秋季 34.37±8.15 14.90±2.42 3.50±0.44 0.87±0.07
冬季 29.39±10.78 13.50±3.67 3.24±0.77 0.83±0.14
F 0.81 0.37 0.39 0.34
P 0.45 0.69 0.68 0.71

为了进一步明确α多样性指数和环境因子之间的相关性, 我们进行了Spearman相关性分析, 结果如表 2所示。Chao 1指数与总碳和总无机碳含量显著正相关(R > 0, P < 0.05)。Richness指数与所有环境因子均不呈现相关性。Shannon指数和Simpson指数均与盐度呈现显著相关关系(R < 0, P < 0.05)。

表 2 浮游纤毛虫α多样性指数与环境因子的Spearman相关性分析 Tab. 2 Spearman's correlation coefficients between α diversity indices and environmental factors
Chao 1 Richness Shannon Simpson
温度/(℃) 0.009 0.051 0.048 0.075
盐度 0.169 0.268 0.290* 0.314*
溶解氧/(mg/L) 0.111 0.051 0.047 0.013
pH −0.027 −0.055 −0.063 −0.072
总碳/(mg/L) 0.297* −0.105 −0.066 −0.018
总无机碳/(mg/L) 0.282* −0.122 −0.090 −0.046
总有机碳/(mg/L) −0.266 −0.053 −0.011 0.031
磷酸盐/(μg/L) 0.041 0.054 0.049 0.047
氨态氮(μg/L) −0.116 −0.095 −0.091 −0.081
亚硝态氮/(μg/L) −0.007 0.038 0.019 0.031
硝态氮/(μg/L) −0.029 0.004 0.009 0.016
硅酸盐/(mg/L) −0.173 −0.086 −0.051 −0.028
注: * 表示P≤0.05, 并以加粗标识。
2.4 β多样性

基于Bray-Curtis距离矩阵对分组后的浮游纤毛虫群落进行非度量多维尺度分析(non-metric multidimensional scaling, NMDS)。如图 4所示, 每个季节的浮游纤毛虫样品聚合性良好且和其他季节样品有明显分离, 说明纤毛虫群落具有季节性。秋冬两季样品有部分交叉并且和春夏季样品显著分开。而不同站点的样品混合散乱排布且重叠明显, 纤毛虫群落未显示出明显的水平分布差异。

图 4 不同季节(a)和不同站点(b)的浮游纤毛虫群落结构的非度量多维尺度图 Fig. 4 Nonmetric multidimensional scaling plots showing the community structure of planktonic ciliates in season (a) and in region (b)

对不同组间浮游纤毛虫群落结构差异性进行ANOSIM检验可知(表 3), 不同站点的纤毛虫群落结构没有显著性差异(ANOSIM, P > 0.05)。而不同季节的纤毛虫群落结构差异显著(ANOSIM, P < 0.01)。

表 3 站点和季节变化对浮游纤毛虫群落β-多样性的ANOSIM分析 Tab. 3 ANOSIM results on changes of plankton ciliates β-diversity in season and in region
分组 R P
站点间 −0.033 0.907
A, B −0.053 0.937
A, C −0.021 0.715
B, C −0.023 0.703
季节间 0.509 0.001
秋季, 冬季 0.206 0.011
秋季, 春季 0.711 0.001
秋季, 夏季 0.717 0.001
冬季, 春季 0.305 0.002
冬季, 夏季 0.604 0.001
春季, 夏季 0.333 0.002
注: 加粗字体表示具有显著差异(P < 0.05)
2.5 与环境因子间的相关关系

为了解影响浮游纤毛虫群落多样性差异的环境因子, 对12个环境因子以及54个样品进行了冗余分析(RDA), 结果如图 5所示。总有机碳和溶解氧是影响三山岛近岸海域浮游纤毛虫群落结构的主要环境因子。总有机碳含量与B站点的部分样品呈现正相关性; 与C站点的夏秋季样品的群落结构呈现正相关性。溶解氧含量与A站点样品和B站点绝大多数样品的群落结构呈现正相关, 且总有机碳与溶解氧含量之间显示出负相关关系。

图 5 浮游纤毛虫群落与环境因子的冗余分析 Fig. 5 The redundancy analyses (RDA) plot showing correlations between ciliate communities and environmental factors

为了解纤毛虫主要类群与环境因子的相关性, 对属(genus)阶元下的无壳寡毛类和砂壳纤毛虫与环境因子进行了Spearman相关性分析(表 4), 与环境因子无相关关系的类群未呈现在表格中。结果显示无壳寡毛类和砂壳类纤毛虫与氨态氮和磷酸盐均无相关性。可以看到在无壳寡毛类中相对序列丰度最高的Choreotrichida_XX与盐度(R < 0, P < 0.05)和碳含量(R < 0, P < 0.01)呈现显著负相关关系; Leegaardiella与硅酸盐呈显著负相关(R < 0, P < 0.01)。砂壳类中相对序列丰度最高的Tintinnopsis与温度(R < 0, P < 0.01)呈显著负相关, 而与溶解氧(R > 0, P < 0.01)、pH (R > 0, P < 0.05)和硅酸盐(R > 0, P < 0.01)呈显著正相关关系。

表 4 属(genus)水平上的浮游纤毛虫主要类群与环境因子的Spearman相关性分析 Tab. 4 Spearman's correlation coefficients (R) between major groups in genus level and environmental factors
2.6 形态学分类鉴定

对A、B两个站点的35个拖网样品进行活体观察和单细胞18S rRNA基因测序, 累计分离鉴定17种纤毛虫(图 6; 表 5), 分别隶属于7个属: 拟铃虫属(Tintinnopsis)、类铃虫属(Codonellopsis)、网纹虫属(Favella)、真铃虫属(Eutintinnus)、前口虫属(Frontonia)、栉毛虫属(Didinium)、旋游虫属(Spirostrombidium), 有9种已测序定种。其中砂壳类有13种, 11种属于黏着壳, 仅有两种透明壳物种。如表 5所示, 砂壳类不仅物种多样性较高, 在检获频率上也远大于其他类群。由于少数物种检获个数和频次少, 未能获得清晰的显微照片或单细胞测序结果, 因此仅能够鉴定到属级阶元或未能呈现。

图 6 浮游纤毛虫的活体显微照片 Fig. 6 Microscopic photograph of planktonic ciliates 注: a. 妥肯丁拟铃虫; b. 卡拉直克拟铃虫; c. 穗缘拟铃虫; d. 罗氏拟铃虫; e. 直颈拟铃虫; f. 圆桶拟铃虫; g. 微细拟铃虫; h. 拟铃虫属; i. 费拉克拟铃虫; j. 清兰拟铃虫; k. 运动类铃虫; l. 巴拿马网纹虫; m. 微小真铃虫; n. 特氏前口虫; o. 前口虫属; p. 栉毛虫属; q. 旋游虫属。比例尺= 50 μm (a~o, q), 25 μm (p)

表 5 三山岛近岸海域浮游纤毛虫种名录 Tab. 5 List of planktonic ciliates species in Sanshan Island
物种名 拉丁文名称 壳长或体长/μm 口径/μm 出现频次
妥肯丁拟铃虫 Tintinnopsis tocantinensis 110~130 30~45 11
卡拉直克拟铃虫 Tintinnopsis karajacensis 90~100 40~45 10
穗缘拟铃虫 Tintinnopsis fimbriata 74~79 45~49 2
罗氏拟铃虫 Tintinnopsis lohmanni 81~91 40~44 3
直颈拟铃虫 Tintinnopsis directa 85~91 35~40 2
圆桶拟铃虫 Tintinnopsis cylindrica 92~136 34~49 8
微细拟铃虫 Tintinnopsis parva 50 23 2
拟铃虫属未知种 Tintinnopsis sp. 73 42 2
费拉克拟铃虫 Tintinnopsis failakkaensis 140 105 1
清兰拟铃虫 Tintinnopsis chinglanensis 109 40 1
运动类铃虫 Codonellopsis mobilis 126~138 72~80 6
巴拿马网纹虫 Favella panamensis 180~230 78~89 6
微小真铃虫 Eutintinnus exigua 140 35 2
特氏前口虫 Frontonia tchibisovae 183 2
前口虫属未知种 Frontonia sp. 133 1
栉毛虫属未知种 Didinium sp. 41 1
旋游虫属未知种 Spirostrombidium sp. 65 1

与高通量测序结果相比, 形态学结果无论是α多样性还是β多样性均较低, 而砂壳类的季节变化趋势在两种方法的结果中较为一致。形态学结果显示拟铃虫属在春夏秋三季均大量出现, 出现频率约为纤毛虫总检获频次的69%, 其中主要在秋季大量检获。这与分子测序结果相似: 拟铃虫属春夏秋三季的序列相对丰度占全年总丰度的80%, 春季的相对丰度最高。形态和分子结果均显示透明壳种类仅在夏秋季高水温时出现。另外, 站点间差异在两种方法的结果中也显示出一致性, 二者均显示A站点寡毛类和砂壳类多样性大于B站点。

3 讨论 3.1 浮游纤毛虫群落组成的时空变动

迄今为止在莱州湾进行的浮游纤毛虫生态研究很少, 仅有陈雪等人(2014, 2015)针对砂壳类纤毛虫丰度的时空分布调查。本研究发现无壳寡毛类和砂壳类纤毛虫是浮游纤毛虫的优势类群, 均呈现明显的季节变化。其中无壳寡毛类在秋冬季的序列相对丰度最高, 主导了冬季纤毛虫的多样性和相对丰度, 因此对季节模式具有较高的贡献度。无壳寡毛类纤毛虫在冬季的生存优势可能与它们的营养模式有关。研究发现有些纤毛虫不仅可以进行吞噬营养, 还可以利用内共生藻类或滞留叶绿体进行光合营养(Bernard et al, 1994; Pettigrosso et al, 2016), 据Stoecker等人(2017)估计, 约30%~40%的寡毛类能够进行混合营养。光合作用可以为混合营养纤毛虫提供大量的有机碳, 因此相比于专性异养的纤毛虫, 混合营养纤毛虫能够通过营养模式的转变, 提高摄食效率、增加碳源, 并降低对氧气、营养盐和食物浓度的依赖程度, 从而更好的适应多变的环境(Esteban et al, 2010)。本研究发现冬季氨态氮含量显著提高, 这种结果在一定程度上与Schoener (2013)的推测相符: 混合营养纤毛虫能够为藻类提供氨态氮, 而非与藻类竞争无机氮。因此推测大多数寡毛类纤毛虫能够降低冬季营养盐变动和食物来源减少带来的不利影响, 从而在冬季样品中占据优势地位, 并能够通过固碳和固氮效应参与浮游生态系统的物质循环。

在其他温带海湾, 如渤海湾(于莹等, 2018, 2022)、Hiroshima湾(Kamiyama et al, 1996)、胶州湾(Chen et al, 2009; Jiang et al, 2011; 于莹等, 2011)、桑沟湾(于莹等, 2013)等的周年调查与本研究中利用HTS得到的结果呈现较高的一致性: 砂壳类纤毛虫的丰度和物种组成呈现明显的季节变化, 且在春夏季达到峰值。另外, 本研究通过镜检观察发现大多为黏着壳类物种, 少数为透明壳物种, 这可能与近岸海水浑浊度高有关(Dolan et al, 2013)。

在C站点的夏季样品中, 不仅物种多样性低而且相对丰度与其他两站点显著不同。其中旋毛纲(Spirotrichea)相对丰度全年最低, 而寡膜纲(Oligohymenophorea)、叶口纲(Litostomatea)和叶咽纲(Phyllopharyngea)的相对丰度显著高于其他两站点。对环境因子的全年监测显示C站点的氮磷元素浓度在夏季出现骤降, 推测氮磷等营养盐含量能够影响浮游纤毛虫的物种组成和优势种的相对丰度。

无论是对浮游纤毛虫种群的α多样性分析还是群落的β多样性分析均未显示出站点间差异, 原因可能是站点设置距离太近。然而在群落组成上, 砂壳目在A站点的相对丰度最高。浮游纤毛虫的空间分布同样会受到环境因子的驱动(Yang et al, 2020), 本研究中溶解氧含量可能是A站点浮游纤毛虫群落的重要控制因素, 前人研究也显示出溶解氧对浮游纤毛虫总生物量的影响(Zhang et al, 2015)。Capriulo等(1982)的研究表明在近岸和上升流区域, 水体扰动能够使一部分沉降的砂壳类克服重力影响而分散在表层水体中。另外, 不同水团也会对近岸浮游纤毛虫群落生态特性的显著影响(Gu et al, 2021)。因此我们推测由于A站点位于码头, 船只靠岸和海浪冲击导致的水体扰动会明显影响近岸环境中砂壳类纤毛虫的水平分布模式。

除环境因子对浮游纤毛虫群落结构产生直接影响外, 生物间相互作用(捕食、寄生、共生等) (Dhib et al, 2013; Elloumi et al, 2014; Fu et al, 2020)或与环境因子的叠加效应能够共同作用于纤毛虫群落(Huang et al, 2021), 从而产生一定的时空分布模式。

3.2 传统形态学与分子生态学方法相结合

在以往对于莱州湾纤毛虫的研究中, 大多利用传统形态学的方法对特定类群进行的丰度和生物量调查。传统的镜检方法较为费时费力, 且镜检样品的体积受到时间和精力的限制(约50 mL), 往往体型较大物种(> 40 μm)容易被检获并鉴定到种水平, 并且鉴定过程需要观察者具备一定的分类学基础和经验。与传统形态学研究方法相比, 基于基因测序的生态调查能够揭示更高的生物多样性(隐性种), 目前对18S V4区进行高通量测序已经广泛应用于真核微生物的生态调查中(Majaneva et al, 2018; Gran-Stadniczeñko et al, 2019; Sun et al, 2019)。

本研究使用高通量测序和显微镜检方法所揭示的季节变动和群落组成基本一致, 且砂壳类纤毛虫群落的β多样性并未受到两种方法样本量不同的影响, 这与之前利用两种方法的调查研究相似(Santoferrara et al, 2016)。另外, 高通量测序结果显示全年寡毛类(40.14%)与砂壳类(10.67%)相对丰度的比例大致为4︰1, 这与胶州湾的调查结果十分接近(寡毛类占总纤毛虫丰度的59.91%, 砂壳类占21.81%) (Jiang et al, 2011), 可能暗示了温带海湾浮游纤毛虫群落具有相似的群落组成模式。

对两种方法得到的结果进行比较, 我们发现砂壳类序列的相对丰度较低(10.67%), 而镜检中的出现频率远高于寡毛类, 导致这种差异可能的原因主要有: (1) 拖网采样对较小个体的遗漏以及皮膜脆弱个体的损坏; (2) 分类学经验不足, 对体型小于40 μm个体的忽略; (3) 高通量测序仅能得到相对丰度数据, 由于类群间rRNA基因拷贝数的差异(寡毛类的rRNA基因拷贝数可能较高), 从而放大了某些物种的相对丰度(Medinger et al., 2010)。

另外, 我们对物种的分类注释结果中显示出大量未知种序列, 如寡毛类中相对丰度最高的Choreotrichida_XX和砂壳类中的Tintinnida_XX, 由于参考序列数据库中基于形态鉴定的寡毛类和砂壳类的序列较少, 数据库不完善, 因此可能会导致OTU分配出现错误, 降低了分类的广度和质量。另外, 高通量测序不能区分活体和胞囊信息, 因此会大大高估环境中纤毛虫多样性。同时结合高通量测序和显微镜检的研究显示, 砂壳类纤毛虫群落的β多样性并未受到两种方法样本量不同的影响(Santoferrara et al, 2016), 但由于OTU信息很可能无法完全对应形态学鉴定的物种(Santoferrara et al, 2020), 该结果仅适用于对纤毛虫群落整体的研究, 而不针对某一类群。

目前我们仍迫切需要环境中的OTU数据来补充纤毛虫生物多样性信息(Warren et al, 2017), 但由于分子测序方法具有一定的局限性, 例如仅能获得物种相对丰度数据、难以确定获得的多样性是真实存在的还是被人为放大, 同时测序结果也会受到测序过程和数据处理方法的影响, 因此应该与形态学方法结合使用、相互补充。未来应利用形态学方法补充完善数据库中参考序列信息, 连接基因种和形态种, 以实现更精确的分类识别。另外由于该海域相关研究的欠缺, 未来也需要进行长时间尺度和更大范围的重复采样, 以获得较为全面的纤毛虫多样性和分布信息, 并通过一定程度的标准化来连接形态学和高通量测序结果。

4 结论

本研究首次利用传统形态学和分子生态学方法相结合的手段, 对莱州湾三山岛近岸海域浮游纤毛虫进行周年生态调查。高通量测序结果显示无壳的寡毛类和有壳的砂壳类是浮游纤毛虫的主要优势类群, 且具有显著的季节变化模式。寡毛类占总纤毛虫序列的40.14%, 在冬季多样性最高, 其中主要优势类群是急游虫(Strombidium)。砂壳类相对丰度为10.67%, 并在春夏季出现峰值。A站点(码头)黏着壳砂壳类的相对丰度明显多于其他站点, 推测水体扰动会影响近岸环境中砂壳类纤毛虫的物种组成和水平分布模式。然而站点间水平分布差异不显著, 原因可能是站点设置距离过近导致的。C站点(养殖厂排水口)夏季样品中寡毛类和砂壳类纤毛虫多样性和相对丰度均低于其他站点, 可能的原因是养殖厂排放的低溶氧尾水会对浮游纤毛虫的群落结构产生影响。在全年12个环境因子中, 总有机碳和溶解氧含量是驱动浮游纤毛虫群落结构变化的主要环境因子。通过传统形态方法共检获17种纤毛虫, 分别隶属于7个属, 多数为砂壳类, 其中拟铃虫属(Tintinnopsis)检获次数最高。本研究主要基于相对丰度数据, 因此某些类群的相对丰度可能会受到其他分类单元绝对丰度变化的影响。目前高通量测序和显微镜检两种方法均有不足, 需要结合使用以获得更为全面的环境生态信息。

参考文献
于莹, 王宇, 张博伦, 等, 2022. 渤海湾浮游纤毛虫丰度和生物量的周年变化. 生态学报, 42(9): 3822-3831
于莹, 张武昌, 赵楠, 等, 2011. 胶州湾浮游纤毛虫丰度和生物量的周年变化. 海洋与湖沼, 42(5): 690-701
于莹, 张武昌, 蒋增杰, 等, 2013. 桑沟湾浮游纤毛虫丰度和生物量分布的季节变化. 海洋学报, 35(3): 215-224 DOI:10.3969/j.issn.0253-4193.2013.03.025
于莹, 周锋, 房恩军, 等, 2018. 渤海湾近岸海域浮游纤毛虫丰度和生物量的季节变化. 渔业科学进展, 39(1): 37-45 DOI:10.11758/yykxjz.20161222001
宋微波, 沃 伦A, 胡晓钟, 2009. 中国黄渤海的自由生纤毛虫. 北京: 科学出版社, 1-515
张武昌, 丰美萍, 于莹, 等, 2012. 砂壳纤毛虫图谱. 北京: 科学出版社, 1-499
陈雪, 吴强, 栾青杉, 等, 2015. 莱州湾大型砂壳纤毛虫丰度的水平分布. 渔业科学进展, 36(2): 11-22 DOI:10.11758/yykxjz.20150202
陈雪, 张武昌, 吴强, 等, 2014. 莱州湾大型砂壳纤毛虫群落季节变化. 生物多样性, 22(5): 649-657 DOI:10.3724/SP.J.1003.2014.13257
黄平平, 赵峰, 徐奎栋, 2017. 基于形态学与核糖体DNA及其cDNA高通量测序的海洋沉积物中纤毛虫多样性比较. 海洋与湖沼, 48(2): 285-296
AGATHA S, 2011. Global diversity of aloricate oligotrichea (Protista, Ciliophora, Spirotricha) in marine and brackish sea water. PLoS One, 6(8): e22466 DOI:10.1371/journal.pone.0022466
BACHY C, DOLAN J R, LÓPEZ-GARCÍA P, et al, 2013. Accuracy of protist diversity assessments: morphology compared with cloning and direct pyrosequencing of 18s rrna genes and its regions using the conspicuous tintinnid ciliates as a case study. The ISME Journal, 7(2): 244-255 DOI:10.1038/ismej.2012.106
BERNARD C, RASSOULZADEGAN F, 1994. Seasonal variations of mixotrophic ciliates in the northwest Mediterranean Sea. Marine Ecology Progress Series, 108: 295-301 DOI:10.3354/meps108295
CAPORASO J G, KUCZYNSKI J, STOMBAUGH J, et al, 2010. Qiime allows analysis of High-throughput community sequencing data. Nature Methods, 7(5): 335-336 DOI:10.1038/nmeth.f.303
CAPRIULO G M, GOLD K, OKUBO A, 1982. Evolution of the lorica in tintinnids: a possible selective advantage. Annales de L'Institut Océanographique, 58: 319-323
CHEN Y H, YANG Y F, 2009. Characteristics of the microzooplankton community in Jiaozhou Bay, Qingdao, China. Chinese Journal of Oceanology and Limnology, 27(3): 435-442 DOI:10.1007/s00343-009-9153-8
CHEUNG M K, AU C H, CHU K H, et al, 2010. Composition and genetic diversity of picoeukaryotes in subtropical coastal waters as revealed by 454 pyrosequencing. The ISME Journal, 4(8): 1053-1059 DOI:10.1038/ismej.2010.26
DE SCHEPPER S, RAY J L, SKAAR K S, et al, 2019. The potential of sedimentary ancient DNA for reconstructing past sea ice evolution. The ISME Journal, 13(10): 2566-2577 DOI:10.1038/s41396-019-0457-1
DE VARGAS C, AUDIC S, Tara Oceans Consortium, et al, 2017. Total V9 rDNA information organized at the metabarcode level for the tara oceans expedition (2009-2012). Pangaea DOI:10.1594/PANGAEA.873277
DHIB A, BRAHIM M B, ZIADI B, et al, 2013. Factors driving the seasonal distribution of planktonic and epiphytic ciliates in a eutrophicated mediterranean lagoon. Marine Pollution Bulletin, 74(1): 383-395 DOI:10.1016/j.marpolbul.2013.06.021
DOLAN J R, 2010. Morphology and ecology in tintinnid ciliates of the marine plankton: correlates of lorica dimensions. Acta Protozoologica, 49(3): 235-244
DOLAN J R, CLAUSTRE H, CARLOTTI F, et al, 2002. Microzooplankton diversity: relationships of tintinnid ciliates with resources, competitors and predators from the atlantic coast of morocco to the eastern mediterranean. Deep Sea Research Part I: Oceanographic Research Papers, 49(7): 1217-1232 DOI:10.1016/S0967-0637(02)00021-3
DOLAN J R, PIERCE R W, 2013. Diversity and distributions of tintinnid ciliates [M] // DOLAN J R, MONTAGNES D J S, AGATHA S, et al. Biology and Ecology of Tintinnid Ciliates: Models for Marine Plankton. Oxford: John Wiley & Sons, Ltd: 214-243.
DOTY S L, DOSHER M R, SINGLETON G L, et al, 2005. Identification of an endophytic rhizobium in stems of populus. Symbiosis, 39(1): 27-35
ELLOUMI J, DRIRA Z, HAMZA A, et al, 2014. Space-time variation of ciliates related to environmental factors in 15 nearshore stations of the Gulf of Gabes. Mediterranean Marine Science, 16(1): 162-179 DOI:10.12681/mms.792
ESTEBAN G F, FENCHEL T, FINLAY B J, 2010. Mixotrophy in ciliates. Protist, 161(5): 621-641 DOI:10.1016/j.protis.2010.08.002
FOISSNER W, CHAO A N, KATZ L A, 2008. Diversity and geographic distribution of ciliates (Protista: Ciliophora). Biodiversity and Conservation, 17(2): 345-363 DOI:10.1007/s10531-007-9254-7
FU Y J, ZHENG P F, ZHANG X L, et al, 2020. Protist interactions and seasonal dynamics in the coast of Yantai, northern Yellow Sea of China as revealed by metabarcoding. Journal of Ocean University of China, 19(4): 961-974 DOI:10.1007/s11802-020-4461-x
GÓMEZ F, 2007. Trends on the distribution of ciliates in the open Pacific Ocean. Acta Oecologica, 32(2): 188-202 DOI:10.1016/j.actao.2007.04.002
GRAN-STADNICZEÑKO S, EGGE E, HOSTYEVA V, et al, 2019. Protist diversity and seasonal dynamics in Skagerrak plankton communities as revealed by metabarcoding and microscopy. Journal of Eukaryotic Microbiology, 66(3): 494-513 DOI:10.1111/jeu.12700
GU B W, WANG Y, XU J, et al, 2021. Water mass shapes the distribution patterns of planktonic ciliates (Alveolata, Ciliophora) in the subtropical Pearl River Estuary. Marine Pollution Bulletin, 167: 112341 DOI:10.1016/j.marpolbul.2021.112341
HUANG H, YANG J P, HUANG S X, et al, 2021. Spatial distribution of planktonic ciliates in the western Pacific Ocean: along the transect from Shenzhen (China) to Pohnpei (Micronesia). Marine Life Science & Technology, 3(1): 103-115 DOI:10.1007/s42995-020-00075-7
JIANG Y, XU H L, AL-RASHEID K A S, et al, 2011. Planktonic ciliate communities in a semi-enclosed bay of Yellow Sea, northern China: annual cycle. Journal of the Marine Biological Association of the United Kingdom, 91(1): 97-105 DOI:10.1017/S002531541000175X
KAMIYAMA T, TSUJINO M, 1996. Seasonal variation in the species composition of tintinnid cilates in Hiroshima Bay, the Seto Inland Sea of Japan. Journal of Plankton Research, 18(12): 2313-2327 DOI:10.1093/plankt/18.12.2313
KOFOID C A, CAMPBELL A S, 1929. A conspectus of the Marine and Freshwater Ciliata Belonging to the Suborder Tintinnoinea, with Descriptions of New Species Principally from the Agassiz Expedition to the Eastern Tropical Pacific 1904-1905. California: University of California Press, 1-403
LÓPEZ-GARCÍA P, RODRÍGUEZ-VALERA F, PEDRÓS-ALIÓ C, et al, 2001. Unexpected diversity of small eukaryotes in deep-sea Antarctic plankton. Nature, 409(6820): 603-607 DOI:10.1038/35054537
LIU W W, JIANG J M, XU Y, et al, 2017. Diversity of free-living marine ciliates (Alveolata, Ciliophora): faunal studies in coastal waters of China during the years 2011-2016. European Journal of Protistology, 61: 424-438 DOI:10.1016/j.ejop.2017.04.007
LYNN D H, 2008. The Ciliated Protozoa: Characterization, Classification, and Guide to the Literature. 3rd ed. Dordrecht, Netherlands: Springer
MAGOČ T, SALZBERG S L, 2011. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics, 27(21): 2957-2963 DOI:10.1093/bioinformatics/btr507
MAJANEVA M, DISERUD O H, EAGLE S H C, et al, 2018. Environmental DNA filtration techniques affect recovered biodiversity. Scientific Reports, 8(1): 4682 DOI:10.1038/s41598-018-23052-8
MEDINGER R, NOLTE V, PANDEY R V, et al, 2010. Diversity in a hidden world: potential and limitation of next-generation sequencing for surveys of molecular diversity of eukaryotic microorganisms. Molecular Ecology, 19: 32-40 DOI:10.1111/j.1365-294X.2009.04478.x
MEDLIN L, ELWOOD H J, STICKEL S, et al, 1988. The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. Gene, 71(2): 491-499 DOI:10.1016/0378-1119(88)90066-2
PETTIGROSSO R E, GARCIA M D, UIBRIG R, et al, 2016. Mixotrophic ciliate dynamics in two zones of a temperate and highly turbid estuary in south America, Argentina. Ecología Austral, 26(2): 107-119
PIERCE R W, TURNER J T, 1992. Ecology of plankton ciliates in marine food webs. Reviews in Aquatic Sciences, 6(2): 139-181
SANTOFERRARA L, BURKI F, FILKER S, et al, 2020. Perspectives from Ten Years of Protist Studies by high-throughput metabarcoding. Journal of Eukaryotic Microbiology, 67(5): 612-622 DOI:10.1111/jeu.12813
SANTOFERRARA L F, GRATTEPANCHE J D, KATZ L A, et al, 2016. Patterns and processes in microbial biogeography: do molecules and morphologies give the same answers?. The ISME Journal, 10(7): 1779-1790 DOI:10.1038/ismej.2015.224
SCHLOSS P D, WESTCOTT S L, RYABIN T, et al, 2009. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology, 75(23): 7537-7541 DOI:10.1128/aem.01541-09
SCHOENER D M, 2013. Inorganic carbon and nitrogen utilization in mixotrophic ciliates [D]. Connecticut: University of Connecticut: 59-60.
SLAPETA J, LÓPEZ-GARCÍA P, MOREIRA D, 2006. Global dispersal and ancient cryptic species in the smallest marine eukaryotes. Molecular Biology and Evolution, 23(1): 23-29 DOI:10.1093/molbev/msj001
STOECKER D K, HANSEN P J, CARON D A, et al, 2017. Mixotrophy in the marine plankton. Annual Review of Marine Science, 9: 311-335 DOI:10.1146/annurev-marine-010816-060617
STOECKER D K, MICHAELS A E, DAVIS L H, 1987. Grazing by the jellyfish, Aurelia aurita, on microzooplankton. Journal of Plankton Research, 9(5): 901-915 DOI:10.1093/plankt/9.5.901
SUN P, HUANG L Y, XU D P, et al, 2019. Integrated space-time dataset reveals high diversity and distinct community structure of ciliates in mesopelagic waters of the northern south China sea. Frontiers in Microbiology, 10: 2178 DOI:10.3389/fmicb.2019.02178
SUZUKI T, MIYABE C, 2007. Ecological balance between ciliate plankton and its prey candidates, pico- and nanoplankton, in the East China Sea. Hydrobiologia, 586(1): 403-410 DOI:10.1007/s10750-007-0715-5
WANG Y F, JI D D, YIN J H, 2019. Morphology and phylogeny of two Phialina species (Ciliophora, Haptoria) from northern China. European Journal of Protistology, 67: 46-58 DOI:10.1016/j.ejop.2018.10.002
WARREN A, PATTERSON D J, DUNTHORN M, et al, 2017. Beyond the "code": a guide to the description and documentation of biodiversity in ciliated protists (Alveolata, Ciliophora). Journal of Eukaryotic Microbiology, 64(4): 539-554 DOI:10.1111/jeu.12391
WEISSE T, MONTAGNES D J S, 2022. Ecology of planktonic ciliates in a changing world: concepts, methods, and challenges. Journal of Eukaryotic Microbiology, 69(5): e12879 DOI:10.1111/jeu.12879
YANG J P, HUANG S X, FAN W X, et al, 2020. Spatial distribution patterns of planktonic ciliate communities in the East China Sea: potential indicators of water masses. Marine Pollution Bulletin, 156: 111253 DOI:10.1016/j.marpolbul.2020.111253
ZHANG C X, ZHANG W C, NI X B, et al, 2015. Influence of different water masses on planktonic ciliate distribution on the East China Sea shelf. Journal of Marine Systems, 141: 98-111 DOI:10.1016/j.jmarsys.2014.09.003
ZHAO F, FILKER S, XU K D, et al, 2017. Patterns and drivers of vertical distribution of the ciliate community from the surface to the abyssopelagic zone in the Western Pacific Ocean. Frontiers in Microbiology, 8: 02559 DOI:10.3389/fmicb.2017.02559