环境科学  2024, Vol. 45 Issue (9): 5298-5307   PDF    
结合生态网络解析白洋淀浮游生物生态位和种间联结性特征
霍笑康1,2, 王永刚2, 周灵同1,2, 王书航1, 姜霞1, 陈珂3, 王鹏飞1     
1. 中国环境科学研究院,湖泊水污染治理与生态修复技术国家工程实验室,国家环境保护湖泊污染控制重点实验室,北京 100012;
2. 兰州理工大学生命科学与工程学院,兰州 730050;
3. 青海省生态环境监测中心,西宁 810002
摘要: 为了解白洋淀浮游生物群落结构及其优势种的生态位特征, 于2022年春季(3月)、夏季(7月)和秋季(9月)对白洋淀浮游生物进行采样调查. 通过构建生态网络图并结合非度量多维尺度分析(NMDS)解析3个季节浮游生物群落结构变化特征, 并运用改进的Levins公式和Petraitis指数评估了浮游生物优势种的生态位宽度和生态位重叠度, 通过卡方检验和种间联结性系数判断总体和优势种种间联结度. 结果表明, 浮游生物在整个淀区的生态位宽度较窄. 浮游动物以轮虫为主, 浮游植物主要是硅藻门、蓝藻门和绿藻门. 浮游生物的群落变化具有明显的季节性特征. 与夏秋季相比, 春季浮游生物种类较少, 种间联通度较低. 夏季浮游动物优势种重叠度高, 种间竞争加剧, 而浮游植物的种间重叠在3个季节均处于较低水平. 夏秋季节浮游植物各优势种种间的联结性均呈显著正联结关系(W > χ0.052), 群落结构稳定. 浮游动物和浮游植物的跨域生态网络在秋季表现出较高的负相关性比例, 尤其在浮游动物的桡足类和枝角类与浮游植物的绿藻门和蓝藻门之间. 白洋淀浮游生物的物种丰富, 季节性差异明显, 优势种以窄生态位为主, 浮游生物群落总体处于较稳定的状态, 桡足类和枝角类与绿藻和蓝藻之间存在较强的捕食关系.
关键词: 白洋淀      浮游生物      生态网络      生态位      种间联结     
Characterization of the Ecological Niche and Interspecific Connectivity of Plankton in Baiyangdian Lake by Combining Ecological Networks
HUO Xiao-kang1,2 , WANG Yong-gang2 , ZHOU Ling-tong1,2 , WANG Shu-hang1 , JIANG Xia1 , CHEN Ke3 , WANG Peng-fei1     
1. National Engineering Laboratory for Lake Pollution Control and Ecological Restoration, State Environment Protection Key Laboratory for Lake Pollution Control, Chinese Research Academy of Environmental Science, Beijing 100012, China;
2. School of Life Science and Engineering, Lanzhou University of Technology, Lanzhou 730050, China;
3. Monitoring Center of Ecology and Environment of Qinghai Province, Xining 810002, China
Abstract: To understand the structure of the plankton community and the ecological niche characteristics of their dominant species, sampling surveys of plankton were conducted in Baiyangdian Lake in the spring (March), summer (July), and autumn (September) of 2022. The changes in the plankton community during the three seasons were analyzed by constructing ecological network diagrams, non-metric multidimensional scaling analysis (NMDS), and the ecological niche width. The niche overlap of zooplankton dominant species was evaluated by the improved Levins' formula and Petraitis' index. The interspecific connectivity of dominant species was judged using the chi-square test and interspecies connectivity coefficients. The results showed that the niche width of plankton in the whole area was low. Zooplankton was dominated by rotifers, and phytoplankton was dominated by diatoms, cyanobacteria, and green algae. There were significant seasonal changes in the community structures of plankton. Compared with that in summer and autumn, there were fewer species of plankton in spring and lower interspecies connectivity. The overlap of dominant species of zooplankton was high in summer, and the interspecific competition was intensified, whereas the interspecific overlap of phytoplankton was at a low level in all three seasons. There was a significant positive correlation (W > χ0.052) between phytoplankton in summer and autumn, and the community structure was stable. The interdomain ecological network of zooplankton and phytoplankton showed a high negative correlation ratio in autumn, especially between copepods and cladoceras of zooplankton and chlorophyta and cyanophyta of phytoplankton. The plankton species in Baiyangdian Lake were abundant, with obvious seasonal differences. The dominant species were mainly a narrow ecological niche. The plankton community was generally in a stable state, and there was a strong predation relationship between copepods and cladoceras and green algae and cyanobacteria.
Key words: Baiyangdian Lake      plankton      ecological network      ecological niche      interspecies connectivity     

浮游生物通常指浮游植物和浮游动物, 是湖泊食物链中不可或缺的存在, 其中浮游植物是湿地生态系统的初级生产者之一, 为消费者如浮游动物, 提供了直接的食物来源[1, 2]. 浮游生物的群落动态变化是指示湖泊生态状况的重要指标[3~5]. 生态位概念描述群落中一个物种与其他物种在特定时间、空间和功能地位的关联性[6], 可表征各种群在群落中的功能地位和生态适应性, 反映物种对资源的竞争强度[7], 已被广泛应用于植物[8]、动物[9]和微生物[10]的群落特征研究.

生态网络分析以随机矩阵理论为基础, 构建生物群落数据的模块化网络[11], 具有较高的可靠性和敏感性. 目前生态网络分析被广泛应用于微生物群落、植物微生物互作等领域. Kruk等[12]通过构建浮游生物生态网络研究了潟湖的浮游动物和浮游植物类群之间的生物结构关系. Goździejewska等[13]通过构建浮游动物生态网络研究了波兰中部的一个矿坑湖在浊度梯度下的浮游动物的结构变化. Yan等[14]通过Illumina MiSeq测序技术研究了陕西金彭水库浮游植物的生态网络.

白洋淀位于雄安新区腹地, 是华北平原上最大的草型浅水湖泊[15], 淀区内台田遍布且沟壕纵横, 开阔水面被分割, 下垫面类型十分复杂[16], 浮游动植物的生境结构多样. 张子良等[17]和陈佳林等[18]分别研究了白洋淀典型区域的浮游动物和浮游植物多样性并进行了水质评价. 目前国内针对浮游动植物的研究主要集中在群落结构调查水平, 缺少对于浮游动植物的生态位和种间联结特征研究.

本研究在白洋淀全淀区春、夏和秋这3个季节浮游生物采样调查基础上, 构建浮游生物生态网络, 计算生态位宽度和生态位重叠度, 并结合非度量多维尺度分析白洋淀浮游生物在季节间的变化情况以及浮游生物的生态位和种间联结度, 通过解析浮游生物群落结构和生态位特征, 以期对湖泊生态系统评估和生物多样性保护提供基础数据支持.

1 材料与方法 1.1 区域概况及样点布设

白洋淀(38°43′~39°02′N, 115°38′~116°07′E)位于河北省中部, 属海河流域大清河南支水系湖泊, 总面积366 km2, 平均年份蓄水量13.2亿m3. 为了全面地获取淀区浮游生物信息, 于2022年春季(3月)、夏季(7月)和秋季(9月)在白洋淀布设36个采样点(图 1), 采样点覆盖了白洋淀各种水域类型, 包括开阔水面、荷塘、沼塘、入淀河口、出淀河口、航道和沟壕等.

图 1 白洋淀浮游生物采样点位示意 Fig. 1 Sampling sites of plankton in Baiyangdian Lake

1.2 样品采集及处理

浮游植物与轮虫样品使用柱状采水器(容积2.5 L)取表层湖水4次, 在10 L容器中充分混合后, 取出1 L并立即加入15 mL鲁哥氏固定液, 静置48 h, 采用虹吸法将上清液缓慢吸出, 剩余30 mL浓缩水样倒入样品瓶中保存. 枝角类和桡足类样品使用柱状采水器采集水样10 L, 经25号浮游生物网过滤后, 浓缩水样加入4%的甲醛溶液保存. 样品带回实验室, 镜检、计数和鉴定[19, 20].

1.3 数据处理 1.3.1 浮游动植物生态网络构建

使用统计分析软件R(Version 4.03 for windows)的“igraph”程序包[21], 基于Spearman相关性系数处理浮游生物群落组成数据, 设定阈值|r|0.4, 去除degree(度)为0的节点, degree较高的节点代表关键物种, 使用Gephi 0.10.1和Cytoscape3.10.1软件构建浮游生物的共现生态网络.

1.3.2 非度量多维尺度分析(NMDS)

使用“dplyr”程序包合并3次采样的浮游生物物种数目数据, 基于Bray-Curtis相似性系数使用“vegan”和“ggplots2”程序包进行了NMDS分析并绘制了图形.

1.3.3 优势度指数

依据优势度指数(Y)确定优势种[22], Y > 0.02则该物种为优势种.

1.3.4 生态位宽度

生态位宽度(Bi)采用加权修正Levins指数[23]计算. 生态位宽度一般分为广生态位(Bi > 0.60)、中生态位(0.30 < Bi ≤ 0.60)和窄生态位(Bi ≤ 0.30)这3种类型[24].

1.3.5 生态位重叠

生态位重叠值(Oik)选用Pianka指数[25]计算. 生态位重叠值越高, 表明种间占用空间和资源利用越相似, 竞争较大;生态位重叠值越低, 表明种间没有太大交集, 也不会争夺相同的资源, 竞争较小.

1.3.6 总体联结性检验

物种间总体的关联性检验采用方差比率法, 分别计算方差比率(VR)和统计量(W[26]. 如果VR > 1, 正相关;如果VR < 1, 负相关;VR=1, 符合所有物种无关联的零假设. 统计量W用来检验VR值偏离1的程度:根据卡方分布, 若W < χ0.95N2W > χ0.05N2, 则物种间总体关联显著(P < 0.05);若χ0.95N2 < W < χ0.05N2, 则优势种种间总体关联不显著(P > 0.05).

1.3.7 优势种种间联结性

种间联结性检验采用Yates的连续校正公式计算卡方值(χ2), 进行卡方检验[27]χ2 < 3.841表示种间联结不显著(P > 0.05);3.841 < χ2 < 6.635表示种间联结显著(0.01 < P < 0.05);χ2 > 6.635表示种间联结极显著(P < 0.01). χ2值不能定量区分联结强度, 需要结合联结系数(AC)对种间联结程度进一步分析, AC值域为[-1, 1], 越接近1, 说明物种间的正联结性越强;越接近-1, 说明负联结性越强;AC=0说明物种间完全独立[28].

生态位宽度、生态位重叠计算和种间联结分析及作图使用“spaa”程序包[29]和Excel 2021. 使用ArcGIS 10.8进行采样点的绘制.

2 结果与分析 2.1 浮游动植物群落组成及优势种

白洋淀全淀区共检出浮游动物133种, 其中轮虫88种, 占总种数的66.16%;枝角类21种, 桡足类23种. 秋季浮游动物种类最多(81种), 春季最少(41种). 白洋淀浮游动物优势种共12种, 除桡足类的广布中剑水蚤外均为轮虫(表 1). 桡足类的无节幼体和桡足幼体在夏季也具有较高的优势度, 因优势种一般指成体, 未计入优势种.

表 1 白洋淀浮游生物优势种的优势度和生态位宽度1) Table 1 Dominance and niche width of plankton dominant species in Baiyangdian Lake

全淀区共检出浮游植物共计8门197种, 其中绿藻门76种, 占比最高(38.57%);硅藻门66种, 占33.50%;蓝藻门27种, 占13.70%. 白洋淀浮游植物优势种共12种, 其中春季4种, 以硅藻为主, 夏季和秋季优势种分别为4种和6种, 均以蓝藻门为主.

2.2 浮游动植物群落生态网络与拓扑特性

分别构建了春夏秋3个季节浮游动物、浮游植物和浮游生物的共现网络. 单独就浮游动物或浮游植物而言, 物种间的相互作用皆以正相关为主, 占比在96.43%及以上, 夏季的正相关百分率均为100%(图 2表 2). 浮游动物和浮游植物的关键物种和多样性存在季节性变化. 浮游动物共现网络图的网络节点数、边数和平均度均表现为秋季 > 夏季 > 春季;而浮游植物为夏季 > 秋季 > 春季. 浮游植物的节点数和边数均比对应季节浮游动物多, 表明浮游植物具有更深层的网络规模(表 2).

节点的大小与节点度的大小有关,针对度较大的节点添加物种信息;绿色连线为正相关,红色连线为负相关,连线的粗细表示相关性系数绝对值的大小 图 2 白洋淀浮游动物和浮游植物群落生态网络 Fig. 2 Ecological network of phytoplankton community in Baiyangdian Lake

表 2 白洋淀浮游生物共现网络拓扑指标 Table 2 Topological index of zooplankton co-occurrence network in Baiyangdian Lake

浮游动物和浮游植物的跨域生态网络在秋季表现出较高的负相关性比例(4.01%), 高于春夏季以及单独的浮游植物或浮游动物, 主要集中在浮游动物的桡足类(剑水蚤幼体)和枝角类(简弧象鼻溞)与浮游植物的绿藻(四足十字藻)和蓝藻(束丝藻)之间(图 3). 3个季节的关键物种均以浮游动物的轮虫类和浮游植物的硅藻门和绿藻门为主. 由春季到秋季, 浮游生物群落的关键物种由浮游植物占比较高逐渐过渡到浮游动物, 以degree排序前15的物种代表关键物种, 春季仅3种为浮游动物, 夏季增加至4种, 秋季为8种且degree最高(图 3表 3).

节点的大小与节点度的大小有关, 仅对度较大的节点添加物种信息;绿色连线为正相关, 红色连线为负相关, 连线的粗细表示相关性系数绝对值的大小, 菱形表示浮游动物, 圆形表示浮游植物 图 3 白洋淀浮游生物跨域生态网络 Fig. 3 Interdomain ecological network of plankton in Baiyangdian Lake

表 3 3个季节跨域生态网中的关键种所属的门(类)和度 Table 3 Phylum and degree of key species in three seasons interdomain ecological network

NMDS分析表明(图 4), 白洋淀不同点位浮游植物种属分布特征春季与夏秋季明显不同, 春季的各采样点表现出明显的聚类效果, 各采样点浮游植物种属相似性较高, 而夏秋季的相似度较低. 与浮游植物相反, 各采样点浮游动物的属种分布在春季比夏秋季更分散, 相似性较低. 整体而言, 与浮游动物相比, 浮游植物在不同季节的分布差异更为显著, 浮游植物种属更易受到季节性变化的影响.

图 4 基于浮游生物季节间物种数目变化的NMDS图 Fig. 4 NMDS plot based on the number of plankton species in seasonal variation

2.3 浮游动植物优势种的生态位 2.3.1 生态位宽度

白洋淀浮游动物优势种的生态位宽度(Bi)在0.083~0.464之间(表 1), 夏季的曲腿龟甲轮虫(Keratella valga)和秋季的广布中剑水蚤(Mesocyclops leuckarti)处于中生态位, 其余优势种均为窄生态位种(类)群. 浮游植物优势种的生态位宽度(Bi)变化范围为0.098~0.445(表 1), 夏季的小环藻(Cyclotella sp.)和秋季的小球藻(Chlorella vulgaris)处于中生态位, 其它物种均处于窄生态位.

2.3.2 生态位重叠值

生态位重叠值(Oik)可以有效地反映浮游生物种类之间对于资源的竞争程度. 浮游动物优势种春、夏、秋季的Oik范围分别为0.078~1.4、0.325~0.903、0.163~0.832, 平均值分别为0.434、0.605、0.452, 夏季 > 秋季 > 春季(图 5);浮游植物优势种春、夏、秋季的Oik范围分别为0.094~0.373、0.002~1.039、0.105~1.035, 平均值分别为0.199、0.394、0.415, 秋季 > 夏季 > 春季(图 5). 3个季节浮游动物生态位重叠值均大于浮游植物, 种间竞争更为激烈.

颜色越深表示生态位重叠值越高;Z1~Z12和P1~P12表示的物种种类见表 1, 下同 图 5 白洋淀浮游动植物优势种生态位重叠值 Fig. 5 Niche overlap of dominant zooplankton species in Baiyangdian Lake

2.3.3 优势种总体关联性

白洋淀浮游生物优势种种间关联性的方差比率VR值均大于1, 表明浮游生物优势种间呈正关联性(表 4). 统计量W值检验结果表明, 夏季和秋季浮游植物优势种总体关联性显著(W < χ0.95N2W > χ0.05N2), 优势种总体联结较为紧密, 群落结构稳定, 外界干扰对夏季和秋季浮游植物群落稳定性影响较小. 而白洋淀浮游动物优势种在春夏秋3个季节以及春季的浮游植物优势种总体正关联性不显著(χ0.95N2 < W < χ0.05N2), 浮游动物和春季的浮游植物优势种种间总体联结较为松散, 趋于随机性, 群落结构不稳定, 外界干扰对群落稳定性影响较大.

表 4 白洋淀优势种总体关联性分析 Table 4 Overall association analysis of dominant plankton species in Baiyangdian Lake

2.3.4 优势种种间联结性

白洋淀浮游生物优势种种间联结性χ2检验结果和种间联结性系数如图 6所示. 浮游动物优势种在3个季节共有22种对, 优势种种间联结关系均不显著(χ2 < 3.841). 浮游植物优势种在3个季节共有36对, 其中有6个种对表现出显著联结关系(3.841 < χ2 < 6.635), 分别是春季的小环藻(P1)和针杆藻(P2)负联结, 夏季的点形平裂藻(P7)和小席藻(P8)、微小平裂藻(P6)和粘球藻(P11)、粘球藻(P11)和鱼腥藻(P12), 秋季的点形平裂藻(P7)和微囊藻(P10)、小席藻(P8)和微囊藻(P10)呈正联结关系. 联结系数(AC)表明, 浮游动物3个季节共计22个种对中共有16个种对表现出正联结关系(AC≥0), 占总数的72.73%;浮游植物3个季节共计36个种对中共有27个种对表现出正联结关系(AC≥0), 占总数的75.00%, 高于浮游动物.

数值表示χ2检验结果, 图形表示种间联结系数半矩阵 图 6 白洋淀浮游生物优势种种间联结 Fig. 6 Interspecific association of plankton dominance in Baiyangdian Lake

3 讨论 3.1 浮游动植物群落特征及优势种演变

白洋淀浮游生物种类数多, 接近黑龙江西部的扎龙湿地[27, 30]和贵州的草海[31], 丰度高且复杂多变, 该类草型湖泊对于生物多样性的保护具有极高的价值.

浮游动物在3个季节均以小型轮虫为主, 轮虫主要为孤雌生殖且世代周期短, 可以在短时间内大量繁殖, 迅速占据优势的生态位[32, 33]. 夏季白洋淀水温适宜, 是桡足类幼虫出生发育的重要时间[34], 因而夏季桡足类无节幼虫和桡足幼虫丰度较高. 白洋淀浮游植物优势种由春季硅藻门为主向夏秋季蓝藻门为主转变, 水温变化可能是一个重要驱动因素[35], 在东平湖和大汶河也观察到同样的演变特征[36]. 春季温度较低, 硅藻门作为狭冷性物种能够在低温环境下通过调节生理代谢活动适应外部环境, 从而在低温且光照强度较弱的条件下取得竞争优势[37], 且硅藻可以分散到冰的底部, 从而固定它们在光区内的位置, 以维持光合作用[38, 39]. 春季白洋淀浮游植物聚类效应明显, 不同采样点的属种分布相似性较高, 而夏秋季分布较分散, 各采样点之间的属种差异较大. 由春季到夏季, 水温的升高影响着浮游植物群落的结构, 小型物种的增殖增加了初级生产力, 加快了群落演替的速度[40, 41].

3.2 浮游动植物的生态位特征

生物所处的生态环境的异质性和种群对于环境的适应和利用能力决定了生态位宽度的大小[24]. 相比于扎龙湿地和白马湖[27, 42], 白洋淀浮游动物优势种的生态位较窄, 一定程度上反映出白洋淀生态环境异质性高于扎龙湿地和白马湖.

两个物种之间的生态位重叠通常被解释为衡量它们对于有限资源的竞争[43]. 浮游动物的生态位重叠值在夏季最高, 而浮游植物的生态位重叠值在夏季最低, 表明在夏季浮游动物对于资源(浮游植物)的竞争激烈, 而浮游植物对资源的竞争并不激烈, 可能与白洋淀夏季水位低、碳氮磷等生源要素含量高有关, 加之不同区域之间的较大异质性, 充足的营养元素及其特殊的环境特征为夏季浮游植物繁殖提供了极好的条件[44, 45].

种间联结是指不同种类物种在时间和空间分布上的相互关联性, 单一的检验方法可能会导致结果产生误差[27]. 本研究结合卡方检验和种间联结系数(AC)两种方法进行检验, 种间联结性判断结果更为客观. 浮游生物优势种的总体关联分析和生态网络图分析均表明春季浮游生物优势种种间总体联结较为松散, 而夏秋较为紧密, 这和小兴凯湖和金彭水库的研究结果一致[46, 47]. 春季浮游生物的稳定性较差可能是由于温度较低, 冬季过后浮游生物种群处于恢复阶段, 物种的丰度较低[48]. 优势种总体检验结果和种间联结系数均表明白洋淀浮游生物种间联结关系以正联结为主, 说明白洋淀在春夏秋这3个季节的环境总体上适合浮游生物生长. 前人研究表明, 负联结比较多的群落处于演替阶段, 随着演替的进行物种间会趋于正联结, 正联结比例越高, 形成稳定群落后, 各种群之间达到动态平衡状态, 种间会呈现总体上的显著正联结, 结合白洋淀的正联结比例说明白洋淀的浮游动植物群落总体处于较稳定的状态[27, 49].

3.3 浮游动植物的生态网络特征

白洋淀浮游动物生态网络节点数、边数和平均度均表现出:秋季 > 夏季 > 春季, 说明秋季浮游动物物种间的连通度更高, 相互作用更复杂, 浮游植物则表现出:夏季 > 秋季 > 春季, 浮游植物在夏季的连通度更高[42, 43]. 结合浮游动物和浮游植物的跨域生态网络, 夏季浮游生物群落种间相互作用紧密, 连通度更高, 对外界扰动响应迅速, 抗干扰能力强. 群落中优势类群间联系密切, 存在模块化现象, 具有相似的生态功能, 尤其体现在秋季浮游动物网络中, 这和渭河基于eDNA得出的结论是一致的[50]. 网络图中的节点之间的作用包括正相互作用和负相互作用, 正相互作用在群落中主要表现为合作关系, 负相互作用主要表现为竞争关系[51]. 本研究表明白洋淀浮游动植物网络拓扑结构均呈现正相关关系(96.43%~100%). 相关研究表明表层水体的细菌网络和真核微生物网络拓扑结构同样呈现较高比例(92.86%~100%)的正相关关系[52], 而土壤、植物微生物的生态网络中正相关性比例(52.6%~98.3%)相对较低[53, 54]. 浮游动植物生活在表层水这一独特的生存环境, 位于食物链金字塔的底端, 缺少水柱、土壤和植物组织等介质的保护缓冲作用[55], 极易受到外界环境因素的影响, 故呈现出近乎百分百的正相关关系.

在春季的群落增长阶段, 关键种(类)群主要为浮游植物, 秋季浮游植物的生长能力变弱, 浮游生物的关键种转变以浮游动物为主. 本研究中, 作为春季优势种的萼花臂尾轮虫在夏秋季不再是优势种(表 1), 而萼花臂尾轮虫可被桡足类和枝角类捕食, 且桡足类和枝角类在跨域网络图中的种类数目逐季增加, 可以推断在白洋淀湖泊生态系统中萼花臂尾轮虫和桡足、枝角类之间存在重要的交流关系[12]. 在由浮游动物和浮游植物构成的跨域生态网络, 两者间的负相关关系表现为捕食关系或寄生关系, 浮游动物的生长在消耗浮游植物[56]. 秋季的浮游生物网络中, 浮游动物的枝角类和桡足类与浮游植物的绿藻门和蓝藻门之间表现出明显的负相关关系, 秋季枝角类和桡足类与绿藻门和蓝藻门之间可能也存在较强的捕食关系.

4 结论

(1)白洋淀浮游生物群落季节性差异显著, 春季种类数目较少, 物种间连通度较低, 生态网络稳定性较差. 浮游植物从春季到夏季的属种变化最为明显, 优势种由春季硅藻向夏秋季蓝藻转变. 浮游动物在3个季节的优势种主要为不同属种的轮虫.

(2)白洋淀由于特殊且复杂的地理环境, 不同区域间的异质性较强, 浮游生物主要为中窄生态位种群, 浮游植物种间竞争程度较低.

(3)浮游生物优势种的总联结性和种间联结性在3个季节均呈现正联结, 生态网络拓扑参数总体呈正相关, 白洋淀浮游生物群落总体处于较稳定的状态. 跨域生态网络分析表明, 秋季枝角类和桡足类与绿藻门和蓝藻门之间可能存在较强的捕食关系.

致谢: 感谢茅台学院的夏伟同学在关于生态位联结方面的问题探讨和中国环境科学研究院邢建宇和李贺在采样过程中的帮助.

参考文献
[1] Majumder A, Adak D, Bairagi N. Phytoplankton-zooplankton interaction under environmental stochasticity: survival, extinction and stability[J]. Applied Mathematical Modelling, 2021, 89: 1382-1404. DOI:10.1016/j.apm.2020.06.076
[2] Schaffner L R, Govaert L, De Meester L, et al. Consumer-resource dynamics is an eco-evolutionary process in a natural plankton community[J]. Nature Ecology & Evolution, 2019, 3(9): 1351-1358.
[3] García-Chicote J, Armengol X, Rojo C. Zooplankton abundance: a neglected key element in the evaluation of reservoir water quality[J]. Limnologica, 2018, 69: 46-54. DOI:10.1016/j.limno.2017.11.004
[4] Ding Y T, Pan B Z, Zhao G N, et al. Geo-climatic factors weaken the effectiveness of phytoplankton diversity as a water quality indicator in a large sediment-laden river[J]. Science of the Total Environment, 2021, 792. DOI:10.1016/j.scitotenv.2021.148346
[5] Oh H J, Jeong H G, Nam G S, et al. Comparison of taxon-based and trophi-based response patterns of rotifer community to water quality: applicability of the rotifer functional group as an indicator of water quality[J]. Animal Cells and Systems, 2017, 21(2): 133-140. DOI:10.1080/19768354.2017.1292952
[6] Pocheville A. The ecological niche: history and recent controversies[A]. In: Heams T, Huneman P, Lecointre G, et al (Eds.). Handbook of Evolutionary Thinking in the Sciences[M]. Dordrecht: Springer, 2015. 547-586.
[7] Parrish J A D, Bazzaz F A. Competitive interactions in plant communities of different successional ages[J]. Ecology, 1982, 63(2): 314-320. DOI:10.2307/1938948
[8] 林岩, 丁晓宇, 吕航, 等. 冰封期乌梁素海浮游植物生态位和种间联结性研究[J]. 水生态学杂志, 2023, 44(3): 102-109.
Lin Y, Ding X Y, Lv H, et al. Ecological niche and interspecific association of phytoplankton during the ice-on period of ulansuhai lake[J]. Journal of Hydroecology, 2023, 44(3): 102-109.
[9] 刘惠, 俞存根, 郑基, 等. 杭州湾灰鳖洋海域秋季主要游泳动物生态位与种间联结[J]. 生态学报, 2019, 39(11): 3828-3836.
Liu H, Yu C G, Zheng J, et al. Niche and interspecific association of major nekton in the Hui-bie-yang sea area of Hangzhou Bay[J]. Acta Ecologica Sinica, 2019, 39(11): 3828-3836.
[10] Luo W, Wang J L, Li Y L, et al. Local domestication of soybean leads to strong root selection and diverse filtration of root-associated bacterial communities[J]. Plant and Soil, 2022, 480(1): 439-455.
[11] Fuhrman J A. Microbial community structure and its functional implications[J]. Nature, 2009, 459(7244): 193-199. DOI:10.1038/nature08058
[12] Kruk M, Paturej E. Indices of trophic and competitive relations in a planktonic network of a shallow, temperate lagoon. A graph and structural equation modeling approach[J]. Ecological Indicators, 2020, 112. DOI:10.1016/j.ecolind.2019.106007
[13] Goździejewska A M, Kruk M. Zooplankton network conditioned by turbidity gradient in small anthropogenic reservoirs[J]. Scientific Reports, 2022, 12(1): 6329. DOI:10.1038/s41598-022-10387-6
[14] Yan M M, Chen S N, Huang T L, et al. Community compositions of phytoplankton and eukaryotes during the mixing periods of a drinking water reservoir: dynamics and interactions[J]. International Journal of Environmental Research and Public Health, 2020, 17(4). DOI:10.3390/ijerph17041128
[15] Li L Q, Chen X H, Zhang M Y, et al. The spatial variations of water quality and effects of water landscape in Baiyangdian Lake, North China[J]. Environmental Science and Pollution Research, 2022, 29(11): 16716-16726. DOI:10.1007/s11356-021-16938-0
[16] Zhao C P, Gong J G, Zeng Q H, et al. Landscape pattern evolution processes and the driving forces in the wetlands of Lake Baiyangdian[J]. Sustainability, 2021, 13(17). DOI:10.3390/su13179747
[17] 张子良, 朱弘阳, 杨志昭, 等. 白洋淀湿地浮游动物群落结构与水质评价[J]. 农业技术与装备, 2023(5): 65-67.
Zhang Z L, Zhu H Y, Yang Z Z, et al. Zooplankton community structure and water quality evaluation in Baiyangdian wetland[J]. Agricultural Technology and Equipment, 2023(5): 65-67. DOI:10.3969/j.issn.1673-887X.2023.05.024
[18] 陈佳林, 孔凡青, 余海军, 等. 白洋淀湿地浮游植物多样性分析及水质评价[J]. 湖北农业科学, 2021, 60(19): 59-64, 74.
Chen J L, Kong F Q, Yu H J, et al. Biodiversity analysis of phytoplankton and water quality evaluation in Baiyangdian Wetland[J]. Hubei Agricultural Sciences, 2021, 60(19): 59-64, 74.
[19] 胡鸿钧, 魏印心. 中国淡水藻类-系统、分类及生态[M]. 北京: 科学出版社, 2006.
[20] 赵文. 水生生物学[M]. 北京: 中国农业出版社, 2005.
[21] Csárdi G, Nepusz T. The igraph software package for complex network research[J]. InterJournal, Complex Systems, 2006, 1695(5): 1-9.
[22] Mcnaughton S J. Relationships among functional properties of californian grassland[J]. Nature, 1967, 216(5111): 168-169.
[23] Colwell R K, Futuyma D J. On the measurement of niche breadth and overlap[J]. Ecology, 1971, 52(4): 567-576. DOI:10.2307/1934144
[24] 张靖泽, 杨波, 肖晶, 等. 贵州南盘江和北盘江后生浮游动物群落结构及优势种生态位特征[J]. 生态学杂志, 2023, 42(9): 2184-2192.
Zhang J Z, Yang B, Xiao J, et al. The niche characteristics of dominant species and community structure of metazooplankton in North Pan River and South Pan River in Guizhou[J]. Chinese Journal of Ecology, 2023, 42(9): 2184-2192.
[25] Pianka E R. The structure of lizard communities[J]. Annual Review of Ecology and Systematics, 1973, 4: 53-74. DOI:10.1146/annurev.es.04.110173.000413
[26] Schluter D. A variance test for detecting species associations, with some example applications[J]. Ecology, 1984, 65(3): 998-1005. DOI:10.2307/1938071
[27] 李雪, 夏伟, 范亚文, 等. 扎龙湿地藻类优势种生态位及其种间联结性动态分析[J]. 生态学报, 2023, 43(10): 4098-4108.
Li X, Xia W, Fan Y W, et al. Dynamic analysis of niche and interspecific association of dominant algae species in Zhalong wetland[J]. Acta Ecologica Sinica, 2023, 43(10): 4098-4108.
[28] 张金屯. 数量生态学[M]. 北京: 科学出版社, 2004.
[29] Zhang J L. spaa: SPecies association analysis[EB/OL]. https://CRAN.R-project.org/package=spaa, 2020-06-09.
[30] 陈庆, 柴一涵, 王天鹅, 等. 基于浮游动物完整性的扎龙湿地水生态状况评价[J]. 中国环境监测, 2022, 38(5): 87-95.
Chen Q, Chai Y H, Wang T E, et al. Evaluation of water ecological status of Zhalong wetland based on zooplankton integrity[J]. Environmental Monitoring in China, 2022, 38(5): 87-95.
[31] 代亮亮, 李莉杰, 何梅, 等. 贵州草海秋季浮游植物群落结构与水质因子的关系[J]. 水生态学杂志, 2020, 41(2): 62-67.
Dai L L, Li L J, He M, et al. Autumn phytoplankton community structure and its relationship with water quality parameters in Caohai Wetland of Guizhou Province[J]. Journal of Hydroecology, 2020, 41(2): 62-67.
[32] Gómez A, Carvalho G R. Sex, parthenogenesis and genetic structure of rotifers: microsatellite analysis of contemporary and resting egg bank populations[J]. Molecular Ecology, 2000, 9(2): 203-214. DOI:10.1046/j.1365-294x.2000.00849.x
[33] García-Roger E M, Lubzens E, Fontaneto D, et al. Facing adversity: dormant embryos in rotifers[J]. The Biological Bulletin, 2019, 237(2): 119-144. DOI:10.1086/705701
[34] da Cruz Lopes da Rosa J, da Silva Matos T, da Silva D C B, et al. Seasonal changes in the size distribution of copepods is affected by coastal upwelling[J]. Diversity, 2023, 15(5). DOI:10.3390/d15050637
[35] Gogoi P, Das S K, Das Sarkar S, et al. Environmental factors driving phytoplankton assemblage pattern and diversity: insights from Sundarban eco-region, India[J]. Ecohydrology & Hydrobiology, 2021, 21(2): 354-367.
[36] 段萍萍, 陈明帅, 宋宁宁, 等. 东平湖与大汶河浮游植物群落结构季节变化及水质评价[J]. 水产学杂志, 2023, 36(4): 89-98.
Duan P P, Chen M S, Song N N, et al. Seasonal variation in phytoplankton community structure and water quality evaluation in Dongping Lake and Dawen River[J]. Chinese Journal of Fisheries, 2023, 36(4): 89-98. DOI:10.3969/j.issn.1005-3832.2023.04.013
[37] 王司阳, 张笑欣, 田世民, 等. 黄河流域干旱区湖泊冰封期浮游植物群落结构特征及影响因子研究[J]. 水利学报, 2020, 51(9): 1070-1079.
Wang S Y, Zhang X X, Tian S M, et al. Study on phytoplankton community structure characteristics and its influencing factors of lakes in arid regions of the Yellow River basin during ice-sealing period[J]. Journal of Hydraulic Engineering, 2020, 51(9): 1070-1079.
[38] Zepernick B N, Chase E E, Denison E R, et al. Declines in ice cover induce light limitation in freshwater diatoms[J]. bioRxiv, 2023. DOI:10.1101/2023.08.02.551672
[39] D'souza N A, Kawarasaki Y, Gantz J D, et al. Diatom assemblages promote ice formation in large lakes[J]. The ISME Journal, 2013, 7(8): 1632-1640. DOI:10.1038/ismej.2013.49
[40] Rasconi S, Gall A, Winter K, et al. Increasing water temperature triggers dominance of small freshwater plankton[J]. PLoS One, 2015, 10(10). DOI:10.1371/journal.Pone.0140449
[41] 王雅雯, 李迎鹤, 张博, 等. 嘉兴南湖不同湖区浮游动植物群落结构特征与环境因子关系[J]. 环境科学, 2022, 43(6): 3106-3117.
Wang Y W, Li Y H, Zhang B, et al. Structural characteristics of zooplankton and phytoplankton communities and its relationship with environmental factors in different regions of Nanhu Lake in Jiaxing City[J]. Environmental Science, 2022, 43(6): 3106-3117.
[42] 范林洁, 胡晓东, 陈文猛, 等. 白马湖浮游动物生态位及其生态分化影响因子[J]. 水生态学杂志, 2022, 43(5): 59-66.
Fan L J, Hu X D, Chen W M, et al. Zooplankton ecological niche and factors affecting ecological differentiation in Baima Lake[J]. Journal of Hydroecology, 2022, 43(5): 59-66.
[43] Macarthur R, Levins R. The limiting similarity, convergence, and divergence of coexisting species[J]. The American Naturalist, 1967, 101(921): 377-385. DOI:10.1086/282505
[44] Wei M J, Gao C, Zhou Y J, et al. Variation in spectral characteristics of dissolved organic matter in inland rivers in various trophic states, and their relationship with phytoplankton[J]. Ecological Indicators, 2019, 104: 321-332. DOI:10.1016/j.ecolind.2019.05.020
[45] Li C C, Feng W Y, Chen H Y, et al. Temporal variation in zooplankton and phytoplankton community species composition and the affecting factors in Lake Taihu-a large freshwater lake in China[J]. Environmental Pollution, 2019, 245: 1050-1057. DOI:10.1016/j.envpol.2018.11.007
[46] Ma C X, Mwagona P C, Yu H X, et al. Seasonal dynamics of zooplankton functional group and its relationship with physic-chemical variables in high turbid nutrient-rich Small Xingkai Wetland Lake, Northeast China[J]. Journal of Freshwater Ecology, 2019, 34(1): 65-79. DOI:10.1080/02705060.2018.1443847
[47] Zhang H H, Yang Y S, Liu X, et al. Novel insights in seasonal dynamics and co-existence patterns of phytoplankton and micro-eukaryotes in drinking water reservoir, Northwest China: DNA data and ecological model[J]. Science of the Total Environment, 2023, 857. DOI:10.1016/j.scitotenv.2022.159160
[48] Ding X, Liu J X, Liu W W, et al. Phytoplankton communities miniaturization driven by extreme weather in subtropical estuary under climate changes[J]. Water Research, 2023, 245. DOI:10.1016/j.watres.2023.120588
[49] 邓贤兰, 刘玉成, 吴杨. 井冈山自然保护区栲属群落优势种群的种间联结关系研究[J]. 植物生态学报, 2003, 27(4): 531-536.
Deng X L, Liu Y C, Wu Y. Interconnection among dominant plant populations of Castanopsis comimunity in Jinggang mountain nature reserve[J]. Acta Phytoecologica Sinica, 2003, 27(4): 531-536. DOI:10.3321/j.issn:1005-264X.2003.04.016
[50] 梁东, 夏军, 宋进喜, 等. 基于eDNA技术的渭河浮游动物多样性及关键种生态位特征[J]. 环境科学, 2021, 42(10): 4708-4716.
Liang D, Xia J, Song J X, et al. Diversity of zooplankton and niche characteristics of keystone species in the Weihe River based on eDNA[J]. Environmental Science, 2021, 42(10): 4708-4716.
[51] Zhang B G, Zhang J, Liu Y, et al. Co-occurrence patterns of soybean rhizosphere microbiome at a continental scale[J]. Soil Biology and Biochemistry, 2018, 118: 178-186. DOI:10.1016/j.soilbio.2017.12.011
[52] 王洵, 廖琴, 王沛芳, 等. 水库水深变化对不同浮游微生物群落及网络互作关键种的影响[J]. 环境科学, 2023, 44(7): 3881-3891.
Wang X, Liao Q, Wang P F, et al. Effects of reservoir water depth on different plankton communities and keystone species of network interaction[J]. Environmental Science, 2023, 44(7): 3881-3891.
[53] Zhang H Q, Wu C F, Wang F Y, et al. Wheat yellow mosaic enhances bacterial deterministic processes in a plant-soil system[J]. Science of the Total Environment, 2022, 812. DOI:10.1016/j.scitotenv.2021.151430
[54] Lv X, Wang Q K, Zhang X Y, et al. Community structure and associated networks of endophytic bacteria in pea roots throughout plant life cycle[J]. Plant and Soil, 2021, 468(1-2): 225-238. DOI:10.1007/s11104-021-05124-3
[55] Yang Y Z, Gao Y C, Huang X N, et al. Adaptive shifts of bacterioplankton communities in response to nitrogen enrichment in a highly polluted river[J]. Environmental Pollution, 2019, 245: 290-299. DOI:10.1016/j.envpol.2018.11.002
[56] Krebs C J. Ecology: the experimental analysis of distribution and abundance[M]. Harlow: Pearson, 2014.