Trophic niche analysis of Collichthys lucidus in Pearl River Estuary
-
摘要:
棘头梅童鱼 (Collichthys lucidus) 是重要的小型经济鱼类，为了解其营养结构和食性等特征，于2021和2022年在珠江口水域采集样本，并基于碳、氮稳定同位素 (δ13C、δ15N) 技术从营养生态位、营养结构指标、食性等角度对其进行研究。结果表明：棘头梅童鱼δ13C为 (−18.74±1.41)‰，δ15N为 (12.35±0.57)‰，营养级为 (3.25±0.17)‰。δ13C在体长小于100.00 mm时随体长增加逐渐上升，之后随体长增加而下降；δ15N在体长 <110.00 mm时随体长增加而上升，在110.00~120.00 mm时出现骤降拐点，δ15N最小，之后随体长增加而上升。在体长 <100.00 mm和 >120.00 mm两阶段的营养多样性较高，营养生态位宽幅较大；在体长100.00~120.00 mm阶段，种群冗余度与均匀度较高，群落内营养生态位分布均匀，拥有相似营养特征的个体占比高。贝叶斯混合模型溯源表明，浮游动物对棘头梅童鱼的食物贡献比例最高。研究结果可为棘头梅童鱼的资源保护与开发提供理论依据。
Abstract:Collichthys lucidus is an important small economic fish, and in order to understand its nutritional structure and dietary characteristics, we collected samples from the waters of the Pearl River Estuary in 2021 and 2022, and studied them from the perspectives of trophic niche, trophic structure index and feeding habit by on stable isotope technology. The results indicate that the δ13C and δ15N values of C. lucidus were (−18.74±1.41)‰ and (12.35±0.57)‰, respectively, with a trophic level of (3.25±0.17)‰. The δ13C value gradually increased when the body length was less than 100.00 mm, then decreased with the increase of body length. The δ15N value increased with the increase of body length when the body length was less than 110.00 mm, and the inflection point occurred at body length of 110.00–120.00 mm. The δ15N value was the lowest, and then increased with the increase of body length. When the body length was less than 100.00 mm and more than 120.00 mm, the trophic diversity was high, and the trophic niche width was large. When the body length was 100.00–120.00 mm, the redundancy and evenness were high, and the trophic niche distribution was uniform in the community, and the individuals with similar trophic characteristics accounted for the majority. Bayesian mixing model tracing reveals that the highest food contribution proportion for C. lucidus was from zooplankton. The research provides a theoretical basis for the resource conservation and management of C. lucidus.
-
Keywords:
- Collichthys lucidus /
- Stable isotope /
- Niche /
- Trophic level /
- Pearl River Estuary
-
近年来，海水中以氮 (N)、磷 (P)为代表的污染物浓度呈波动增加的趋势，近岸海域富营养化问题日益突出[1-3]。N、P等污染物的来源，除陆源地表径流携带输入外，局部海域密集的海水养殖也成为富营养化的重要因素[4-8]。按照《水产养殖业污染源产排污系数手册》估算，2012年我国沿海水产养殖总氮 (Total nitrogen, TN) 和总磷 (Total phosphorus, TP) 排放量已分别达17 414 t和3 146 t，而同年监测的全国425个直排海污染源 (日排放量>100 m3) 氨氮和TP排放量分别为17 100 t和2 921 t，海水养殖总产量与沿海省份污染负荷总量存在显著的线性相关[9-10]。海水养殖污染物排放相关的研究表明，高位池海水养殖投喂饲料中仅约20%的N和30%的P转化为蛋白质，大量未经利用的残饵以及养殖排泄物最终进入海洋环境，进一步加剧了富营养化程度并由此引发了藻华等灾害事件[11-12]。以广东省海洋环境状况为例，海水养殖业分布与劣四类海水的主要分布区高度重叠[13]。
微藻作为海洋生态系统初级生产力的重要组成部分，对维持生态系统稳定和健康具有不可替代的作用。三角褐指藻 (Phaeodactylum tricornutum)、等鞭金藻 (Isochrysis sp.)、亚心形扁藻 (Platymonas subcordiformis) 等适应能力强，生长迅速的微藻被广泛作为海水养殖的天然饵料。然而，沿海频繁暴发的藻华等灾害也对海上生产、生活带来了极显著的影响。广东省近岸海域藻华记录种达60种以上，其中有毒有害种类主要有亚历山大藻 (Alexandrium sp.)、膝沟藻 (Gonyaulax sp.)、裸甲藻 (Gymnodinium sp.)、夜光藻 (Noctiluca scintillans)、中肋骨条藻 (Skeletonema costatum)、棕囊藻 (Phaeocystis sp.)以及旋沟藻 (Cochlodinium schütt)等。富营养化与藻华暴发之间存在显著的因果关系，但一定浓度的N、P又是微藻生长所必需的营养物质[14-15]，如何科学利用环境良好型微藻破解富营养化问题和赤潮灾害成为生物操纵领域的研究热点之一。目前，生态化海水养殖模式中对养殖排放N、P等污染物的综合利用途径之一就是培育优良的饵料微藻。一方面通过微藻、细菌等的共同作用，分解转化养殖代谢产物，并将N、P和二氧化碳 (CO2) 再次同化利用，降低养殖塘中氨氮和磷酸盐的浓度，同时释放出氧气；另一方面通过种间竞争作用，抑制其他有害藻类繁殖，且为自身提供优良饵料，进而改善和维持良好的养殖生态系统[16-20]。
为筛选适宜海水养殖尾水处理的优良饵料微藻，对比其生长特性和氮、磷营养物质利用效率，选取5种赤潮藻和2种可作为饵料的微藻为研究对象，开展不同生源要素[碳(C)、N、P]浓度对海洋微藻增长实验，氮源形态 (硝酸盐、铵盐和尿素) 对微藻细胞增长的影响实验以及微藻，对氮、磷营养盐的吸收利用特性研究，对比分析经济微藻和赤潮藻的营养盐利用特性，以期筛选出适宜的饵料微藻，开展海水养殖尾水的藻类生物去除，实现养殖尾水的循环利用。
1. 材料与方法
1.1 实验藻种和培养条件
藻华种类微藻包括甲藻类双胞旋沟藻 (Cochlodinium geminatum)、米氏凯伦藻 (Karenia mikimotoi)；硅藻类中肋骨条藻 (Skeletonema costatum)、牟氏角毛藻 (Chaetoceros muelleri)；金藻类球形棕囊藻 (Phaeocystis globosa)。饵料微藻为绿藻类亚心形扁藻 (Platymonas subcordiformis)；金藻类眼点拟微绿球藻 (Nannochloropsis oculata)。上述藻种均取自暨南大学水生生物研究中心 (生态学系)，培养基采用人工海水 [ 不添加碳酸氢钠 (NaHCO3)]，添加相应的组分配置成f/2培养液[21-22]。其中碳源 (NaHCO3)、磷源[ 磷酸二氢钠 (NaH2PO4·H2O)]和氮源 {[ 碳酸铵 (NH)2CO、氯化铵 (NH4Cl)、硝酸铵 (NaNO3)]}浓度分别为2.02、3.63 × 10−2和8.83 × 10−2 mmol·L−1。
微藻采用人工气候培养箱 (GXZ-160D，宁波江南仪器厂) 培养，光照强度为约为4 000 lx，光暗比为12 h∶12 h，温度为24 ℃。微藻培养所用器皿和培养基等均经高压蒸汽灭菌 (120 ℃，20 min)，冷却后使用，实验操作均在超净工作台 (SW-CJ-lF型，苏州) 内进行。
1.2 实验方法
1.2.1 碳氮磷生源要素对微藻种群增长的影响
预处理：将指数生长期的实验微藻，分别接种到C、N、P单一要素缺失的f/2培养基中饥饿培养4 d，浓缩收集预培养的微藻进行接种。
生源要素浓度梯度设置采用稀释法，浓度为正常培养基的1/500、1/50、1/5，对照组不添加相应碳/氮/磷。双胞旋沟藻起始密度为3 000个·mL−1，其余6种藻细胞起始密度约5 × 104 个·mL−1。实验采用50 mL玻璃管，接种培养液20 mL，每组设置3个平行，隔天测定细胞密度，并计算第12天的比生长速率，实验的光照、温度等培养条件同1.1。
1.2.2 微藻对硝酸盐和磷酸盐的吸收
时间序列实验中，藻细胞起始密度均设为3 000 个·mL−1，用100 mL锥形瓶，接种f/2培养液60 mL，分别在第1、第3、第5、第7、第9 天取样测定藻细胞密度和培养液中的硝酸盐、磷酸盐含量。
浓度序列实验中，硝酸盐/磷酸盐浓度设定为正常f/2培养基的1/10、1/20、1/40、1/80、1/160、1/320，双胞旋沟藻细胞起始密度为2 000 个·mL−1，其余微藻细胞密度约为8 × 103 个·mL−1，用50 mL玻璃管，接种培养液20 mL，培养时间5 h，实验的光照、温度等培养条件同1.1。
1.3 生长指标测定
1.3.1 藻细胞密度测定
利用酶标仪 (Infinite 200 Pro, Tecan, 瑞典) 测量藻的光密度 (OD680)。通过前期实验建立藻细胞密度 (N, 个·mL−1) 和光密度的回归方程计算藻细胞密度。
1.3.2 比生长速率 (μ, d−1) 计算
公式为μ= (lnNt−lnN0)/(Tt−T0)，式中Nt为在时间Tt时刻的细胞数，N0为接种时的细胞数。
1.3.3 硝酸盐和磷酸盐浓度的测定
藻类样品采集后经0.45 μm滤膜过滤，采用AA3-连续流动元素分析仪 (AutoAnalyzer3，BRAN+LUEBBE，德国) 测定N、P元素含量。
1.4 统计与分析
实验数据采用Origin 8.0整理并制图，采用SPSS 13.0软件包进行统计分析，单因素方差分析 (One-Way ANOVA) 用来检查处理之间的差异性，P<0.05和 P<0.01分别认为差异显著和极其显著。
2. 结果
2.1 微藻在不同C、N、P浓度下的生长
7种微藻在不同C、P浓度下培养11 d后，藻细胞比生长速率变化情况见图1。碳源浓度增加对双胞旋沟藻、米氏凯伦藻、中肋骨条藻、亚心形扁藻和眼点拟微绿球藻的生长速率表现出促进作用，与不添加C、P的对照组相比，最大比生长速率分别提高了35%、7%、19%、9%和19%，双胞旋沟藻比生长速率的增加比例高于其他种类微藻。磷源浓度增加对7种微藻比生长速率的最大促进作用范围介于25.01%~446.60%，且与对照组相比，其对7种微藻的生长速度均表现出显著的促进作用，最大促进作用在7种微藻间的大小顺序是双胞旋沟藻<中肋骨条藻<眼点拟微绿球藻<球形棕囊藻<米氏凯伦藻<亚心形扁藻<牟氏角毛藻。从7种微藻在磷处理对照组的比生长速率情况看，磷饥饿处理对角毛藻生长的抑制效应显著高于其他微藻，另外，米氏凯伦藻比生长速率在0.064 μmol·L−1磷处理组即表现出了显著的促进作用。
![]()
图 1 不同碳、磷起始浓度培养11 d 后7种微藻的生长速率*. 表示处理组与对照组有显著差异 (P<0.05) ，后图同此Figure 1. Specific growth rate of algae under different concentrations of bio-carbonate and phosphate after 11 days of culture*. Significant difference between treatments and control (P<0.05); the same case in the following figures3种形态的氮源培养，7种微藻比生长速率均随氮浓度增加而升高 (图2)。与不添加氮源的对照组相比，7种微藻最大比生长速率在硝酸盐处理组增加了16.54%~77.52%；在铵盐处理组增加了15.79%~88.82%；在尿素处理组增加了25.16%~71.43%。3种氮源含量增加，米氏凯伦藻最大比生长速率增加量最高，其次为牟氏角毛藻，球形棕囊藻细胞净增加量最低。从7种微藻空白对照组比生长速率情况来看，双胞旋沟藻、米氏凯伦藻、牟氏角毛藻对照组仍受到氮饥饿处理影响，抑制效应显著，但随着氮浓度增加，细胞生长速率迅速增加；球形棕囊藻对照组仍保持较高的比生长速率，也反映出在前期的氮饥饿处理后细胞内仍残留部分的氮源，也表明球形棕囊藻具有转化储存多余氮源的细胞调节机制。中肋骨条藻、亚心形扁藻和眼点拟微绿球藻对照组也表现出了较高水平的比生长速率，且随着浓度增加，比生长速率均表现出显著增加趋势，显示出较强的氮源竞争能力。
![]()
图 2 氮源形态和浓度对微藻比生长速率的影响Figure 2. Specific growth rate of algae under different concentrations of nitrate, ammonium and urea双胞旋沟藻最大比生长速率在硝酸盐和尿素处理组无显著差异，但显著高于铵盐处理组(P<0.05)。米氏凯伦藻、亚心形扁藻和眼点拟微绿球藻的最大比生长速率为铵盐>硝酸盐>尿素，前者在铵盐和硝酸盐处理组的最大比生长速率无显著差异 (P>0.05)，后者在硝酸盐和尿素处理组无显著差异 (P>0.05)。牟氏角毛藻最大比生长速率为硝酸盐>铵盐>尿素；球形棕囊藻和中肋骨条藻最大比生长速率在3种氮源下无显著差异 (P>0.05)。
2.2 微藻对培养液中N、P的吸收利用
相同细胞起始密度下培养9 d后，7种微藻的种群增长曲线和比生长速率见图3。比生长速率大小顺序为眼点拟微绿球藻>球形棕囊藻>牟氏角毛藻>中肋骨条藻>亚心形扁藻>米氏凯伦藻>双胞旋沟藻。7种微藻细胞增长的规律，从细胞大小看，细胞体积小的种类显著高于细胞体积大的。这一趋势在同一门类的微藻中也能观察到，细胞体积越小，比生长速率越高。
![]()
图 3 相同细胞起始密度培养9 d后微藻种群增长曲线和比生长速率不同字母表示处理组间有显著差异 (P<0.05)，后图同此Figure 3. Population growth curve and specific growth rate of same algae cell after 9 days of cultureValues with different letters indicate significant difference at P=0.05 level. The same cases in the following figures相同细胞起始密度下，微藻培养液中的N、P含量随培养时间延长呈显著降低的趋势(图4、图5)。根据第9天的细胞密度增量以及培养液中N、P含量计算出7种微藻对单位氮原子的利用效率为：眼点拟微绿球藻>球形棕囊藻>牟氏角毛藻>中肋骨条藻>米氏凯伦藻>亚心形扁藻>双胞旋沟藻；7种微藻对单位磷原子的利用效率为：眼点拟微绿球藻>球形棕囊藻>牟氏角毛藻>中肋骨条藻>亚心形扁藻>米氏凯伦藻>双胞旋沟藻。7种微藻在单位N、P原子下的种群增长与细胞体积呈一定的负相关关系。
![]()
图 4 相同细胞起始密度下培养液中氮含量变化和单位藻细胞的氮元素利用效率Figure 4. Variation of nitrogen concentrations in culture medium and nitrogen utilization efficiency at same algae cell density after 9 days of culture![]()
图 5 相同细胞起始密度下培养液中磷含量变化和单位藻细胞的磷元素利用效率Figure 5. Variation of phosphorus concentrations in culture medium and phosphorus utilization efficiency at same algae cell density培养液中N、P原子质量比变化情况见表1。7种微藻培养液中氮磷比范围介于1.10~19.56，硅藻、绿藻和金藻种类培养液中氮磷比变化情况显著高于甲藻，中肋骨条藻在第9天时培养液中氮磷比最大；7种微藻培养液中氮磷比均呈由小到大的增加趋势，这一现象主要与微藻对P的吸收速率相对较高有关。
表 1 微藻培养液中氮磷比变化Table 1. Variations of N/P ratio in culture medium培养时间/d
Culture time双胞旋沟藻
C. geminatum米氏凯伦藻
K. mikimotoi球形棕囊藻
P. globosa中肋骨条藻
S. costatum牟氏角毛藻
C. mulleri亚心形扁藻
P. subcordiformis眼点拟微绿球藻
N. oculata0 1.10 1.10 1.10 1.10 1.10 1.10 1.10 1 2.84 4.04 4.12 10.98 6.09 9.88 5.95 3 3.01 4.83 4.97 11.72 7.07 11.60 6.86 5 3.19 5.59 5.52 11.45 7.11 11.81 7.58 7 3.52 6.51 6.92 12.81 8.69 11.68 7.20 9 2.98 6.95 10.20 19.56 9.35 13.73 12.47 微藻对不同起始浓度N、P的吸收速率见图6、图7。7种微藻对不同起始浓度的N、P吸收速率存在一定的差异。双胞旋沟藻对N、P的吸收表现出饱和现象，最大吸收率分别为5%和4%，主要受起始细胞密度较低的限制。其余6种微藻对不同浓度N的吸收率均高于90%；米氏凯伦藻对磷的吸收率为60%，其余5种微藻对P的吸收率均达到90%。从单位细胞密度的氮磷吸收效率上看，7种微藻的N吸收效率顺序为亚心形扁藻>米氏凯伦藻>双胞旋沟藻>球形棕囊藻>中肋骨条藻>牟氏角毛藻>眼点拟微绿球藻；7种微藻P吸收效率顺序为亚心形扁藻>双胞旋沟藻>米氏凯伦藻>球形棕囊藻>中肋骨条藻>牟氏角毛藻>眼点拟微绿球藻。7种海洋微藻单位细胞密度的氮磷吸收效率在细胞体积较大的种类中高于细胞体积小的种类；而从微藻所属门类上看，绿藻门>甲藻门>金藻门>硅藻门。
![]()
图 6 不同氮起始浓度下培养5 h藻细胞氮吸收量和单位细胞吸收速率Figure 6. Nitrogen absorption of seven algae under different initial concentrations of nitrate after culture of 5 h![]()
图 7 不同磷起始浓度下培养5 h的藻细胞磷吸收量和单位细胞吸收效率Figure 7. Phosphorus absorption of seven algae under different initial concentrations of phosphate after culture of 5 h3. 讨论
一定浓度范围内，生源要素 (C、N、P) 含量的增加能够显著促进海洋微藻的生长，而浓度过高或过低时又成为某些海洋微藻种群增长的限制因子。与温室效应、富营养化相伴而生的海洋赤潮或有害藻类水华已成为我国近岸海域主要的环境问题之一[2, 15, 23-24]。米氏凯伦藻、球形棕囊藻、中肋骨条藻、牟氏角毛藻是南海海域常见的藻华种，而双胞旋沟藻在2006—2013年间的间断性暴发成为南海近岸海域新的环境热点问题[25-26]。海水养殖等生产活动直接排放的氮磷等营养物质逐渐成为赤潮灾害的重要诱发因素[4, 8, 11-13, 24]。
Wang等[27]的研究表明，海洋藻类Rubisco的Km (CO2) 约为40~60 μmol·L−1，但实际海水中CO2浓度约为10 μmol·L−1 (pH 8.0, 20 ℃)，因此海水中低溶解态的CO2也可能成为某些海洋微藻光合固碳的限制因子。本研究显示，在C浓度0~400 μmol·L−1条件下，双胞旋沟藻、眼点拟微绿球藻、中肋骨条藻比生长速率表现出显著的促进作用，最大增加量为35% (双胞旋沟藻)，其他微藻比生长速率未表现出显著增加，这一现象可能与藻细胞自身调整CCM (CO2 concentrating mechanism, CCM) 对CO2的亲和力有关[28]。本研究发现，在其他营养盐充足的条件下，碳源浓度增加，对双胞旋沟藻、中肋骨条藻的种群增长表现出一定的促进作用，但碳源浓度增加对7种微藻种群增长的总体促进作用影响较小，类似的研究结果在聚球藻 (Synechococcus sp.) 中也得到证实。聚球藻从高无机碳浓度培养条件转移到低浓度培养时，其胞外碳酸酐酶在4~6 h内即可恢复高CO2亲和力[29]。因此，温室效应导致的海洋溶解无机碳增加并非是有害藻类赤潮频发的关键诱导因素。
近岸海域藻华种微藻种类日益增多的现状，除部分种类随船舶压舱排放直接引发藻华外，大多数藻华种类微藻已逐步适应并呈现固定季节暴发的趋势[30]。因此，近岸海域氮、磷富营养化被认定为各类藻华暴发的基础。根据近几年的《中国近岸海域环境质量公报》，我国近岸海域四大海区的平均无机氮浓度已达11~35 μmol·L−1，平均活性磷酸盐浓度为0.6~1.8 μmol·L−1。本研究设定的N (0~176 μmol·L−1)、P (0~6.4 μmol·L−1) 浓度涵盖了我国绝大多数海域的水质状况，P浓度对照组比生长速率与最大浓度 (6.4 μmol·L−1)组相比下降了20%~81.71%。从门类上看，甲藻 (双胞旋沟藻、米氏凯伦藻)、硅藻 (中肋骨条藻、牟氏角毛藻)、金藻 (球形棕囊藻) 的平均比生长速率分别下降了28.58%、52.71%和35.66%。2种饵料微藻 (亚心形扁藻、眼点拟微绿球藻) 在设定的P浓度范围内，比生长速率仅明显低于球形棕囊藻，与其他藻华种的比生长速率相比无显著差别。微藻对磷源的的利用差异，在其他研究中也有证实，微藻对磷源的利用特性不仅与磷源形态 (正磷酸盐、有机磷等) 有关，而且受藻细胞磷酸酶活性的影响显著[15, 23-24, 31]。卵圆卡盾藻 (Chattonella ovate) 在无机磷源缺乏的条件下可以利用有机磷源，中肋骨条藻和东海原甲藻 (Prorocentrum donghaiense) 无机磷利用效率低于有机磷源[32]。
氮源浓度和形态对微藻细胞增长和生理生化组分的影响研究已有广泛开展，如东海原甲藻、多环旋沟藻 (Cochlodinium polykrikoides)、海洋卡盾藻 (C. marina)、赤潮异弯藻 (Heterosigma akashiwo)等[15, 23-24, 26, 33]。本研究设置的3种氮源形态培养下，硝酸盐对照组比生长速率与最大浓度组相比下降了16.54%~77.52%，甲藻 (双胞旋沟藻、米氏凯伦藻)、硅藻(中肋骨条藻、牟氏角毛藻)、金藻(球形棕囊藻)的平均抑制率分别为47.03%、58.13%、30.59%；铵盐对照组比生长速率与最大浓度组相比下降了15.79%~87.82%，甲藻 (双胞旋沟藻、米氏凯伦藻)、硅藻 (中肋骨条藻、牟氏角毛藻)、金藻 (球形棕囊藻)的平均抑制率分别为81.45%、43.96%、15.79；尿素对照组比生长速率与最大浓度组相比下降了25.16%~71.43%，甲藻 (双胞旋沟藻、米氏凯伦藻)、硅藻 (中肋骨条藻、牟氏角毛藻)、金藻 (球形棕囊藻) 的平均抑制率分别为69.31%、53.98%和29.55%。2种饵料微藻 (亚心形扁藻、眼点拟微绿球藻) 在3种氮源下的最大比生长速率仅明显低于球形棕囊藻，处于较高水平，也进一步表明，受试饵料微藻在氮源利用效率上具有较高水平。
对比3种氮源形态下藻华种类微藻的种群增长，甲藻 (双胞旋沟藻、米氏凯伦藻) 对铵盐和尿素限制最敏感；硅藻 (中肋骨条藻、牟氏角毛藻) 对硝酸盐和尿素限制敏感性更高；金藻 (球形棕囊藻) 对氮源限制敏感性最弱。上述结果在一定程度上解释了近岸海域藻华演变先硅藻后甲藻的趋势，在营养盐丰富的条件下，硅藻种群增长迅速；随着营养盐消耗，硅藻种群增长更易受到磷源限制而表现出抑制效应。同时，硝酸盐的大量消耗也进一步限制了硅藻门微藻的种群增长。而甲藻门种类可在低浓度硝酸盐和磷酸盐条件下维持较高的种群增长速率，并因此获得种间竞争优势。金藻门球形棕囊藻受到氮、磷限制的响应最低，也证实其种群增长具有明显的营养盐利用优势，这一优势有利于其在低浓度营养盐条件下快速增长并长时间维持高种群密度。甲藻门 (裸甲藻) 在低磷条件下仍能维持较高的种群增长，且与硅藻门 (假微型海链藻) 的种间竞争实验结果均呈现先硅藻后甲藻的群落演替规律[34]。然而，也有研究表明，受试微藻种类的不同导致其对氮、磷限制的响应也存在差异，甲藻 (赤潮异弯藻) 对磷限制的敏感性显著高于氮限制，而同属针胞藻纲的海洋卡盾藻却对相同条件下的氮、磷限制响应无显著差异[33]。本研究发现金藻 (球形棕囊藻) 和非赤潮藻亚心形扁藻、眼点拟微绿球藻在低氮、磷浓度下均保持较高水平的生长速率，这一特性在一定程度上为经济藻类处理利用养殖废水，并维持相对的种间竞争优势提供了生物学基础。
微藻在不同氮、磷浓度下的生长情况，受其对氮、磷吸收利用效率的影响。本研究发现，微藻比生长速率与细胞体积呈一定的负相关关系，而单位细胞密度的氮、磷营养盐吸收速率却与细胞体积呈正相关。从本研究采用的受试微藻门类上看，硅藻 (中肋骨条藻、牟氏角毛藻) 生长速率显著高于甲藻 (双胞旋沟藻、米氏凯伦藻)，研究结果与近岸海域浮游植物群落演替的调查结果相符合[25]，金藻 (球形棕囊藻) 细胞增长和氮、磷营养盐吸收利用均高于上述甲藻和硅藻，进一步揭示了富营养条件下球形棕囊藻种群暴发赤潮的潜在生理学基础。本研究选取的亚心形扁藻、眼点拟微绿球藻在生源要素的吸收利用以及细胞生长速率上均能够达到甚至超过部分典型藻华种，表现出较好的营养盐吸收利用特性，能够作为养殖尾水氮、磷去除的优良藻株，后续拟综合其他环境因素 (光照、温度、pH等) 构建微藻混合培养，进一步探索亚心形扁藻和眼点拟微绿球藻作为养殖尾水的原位生物处理能效。
-
![]()
图 1 棘头梅童鱼各体长组稳定同位素特征图
Figure 1. Stable isotope characteristics of each body length group of C. lucidus
![]()
图 2 棘头梅童鱼δ13C、δ15N与体长的相关性
Figure 2. Correlation between δ13C, δ15N values and body length of C. lucidus
![]()
图 4 棘头梅童鱼各体长组营养级变化图
Figure 4. Trophic level variation of each body length group of C. lucidus
![]()
图 5 棘头梅童鱼各体长组营养生态位特征图
Figure 5. Trophic niche characteristics of each body length group of C. lucidus
![]()
图 6 棘头梅童鱼各体长组的食物碳源贡献比
Figure 6. Contribution ratio of carbon sources to food of each body length group of C. lucidus
表 1 棘头梅童鱼体长组划分及其稳定同位素特征
Table 1 Classification of body length groups and stable isotope characteristics of C. lucidus
体长
Body length (L)/mm样本数
Number of samples/尾平均体长
Average length/mm营养级
Trophic levelδ15N/‰ δ13C/‰ L≤90 11 64.43±17.31 3.21±0.15 12.21±0.52ab −20.69±1.81a 90＜L≤100 23 96.36±2.19 3.28±0.14 12.45±0.47a −17.78±0.82d 100＜L≤110 34 104.71±2.34 3.28±0.12 12.46±0.40a −17.89±0.79d 110＜L≤120 40 114.97±3.17 3.19±0.12 12.15±0.40a −18.43±0.81c 120＜L≤130 35 123.84±2.51 3.25±0.19 12.34±0.65ab −19.36±1.23b 130＜L 15 137.04±12.39 3.34±0.28 12.65±0.94a −20.06±1.56a 注：同列不同字母表示组间差异显著 (P<0.05)。 Note: Values with different letters within the same column are significantly different (P<0.05). 
下载: 导出CSV
表 2 棘头梅童鱼各体长组的营养结构指标
Table 2 Trophic structure indexes of each body length group of C. lucidus
体长
Body length (L)/mmδ13C差值
CRδ15N差值
NR平均离心距离
CD平均最邻近距离
MNND平均最邻近距离标准差
SDNND多边形总面积
TA矫正后的标准椭圆面积
SEAcL≤90 5.79 1.87 1.59 0.61 0.52 5.93 2.92 90<L≤100 4.17 2.01 0.60 0.27 0.65 2.86 0.97 100<L≤110 3.53 1.75 0.65 0.20 0.27 3.18 0.78 110<L≤120 2.96 1.76 0.77 0.21 0.18 3.84 0.99 120<L≤130 5.33 3.71 1.09 0.26 0.35 8.35 1.88 130<L 5.13 3.41 1.52 0.55 0.42 6.45 2.73 Note: CR. δ13C range; NR. δ15N range; CD. Mean distance to centroid; MNND. Mean nearest neighbor distance; SDNND. Standard deviation of nearest neighbor distance; TA. Total area; SEAc. Corrected standard ellipse area.
下载: 导出CSV
表 3 棘头梅童鱼各体长组营养生态位椭圆重叠面积比率
Table 3 Ratio of overlap area of trophic niche ellipses in each body length group of C. lucidus
体长
Body length (L)/mmL≤90 90<L≤100 100<L≤110 110<L≤120 120<L≤130 130<L L≤90 2.92 0.00% 0.00% 4.76% 32.30% 34.23% 90<L≤100 0.00% 0.97 75.24% 13.21% 0.00% 0.00% 100<L≤110 0.00% 93.24% 0.78 22.54% 0.00% 0.00% 110<L≤120 14.09% 12.94% 17.82% 0.99 58.44% 33.56% 120<L≤130 50.12% 0.00% 0.00% 30.66% 1.88 83.42% 130<L 36.58% 0.00% 0.00% 12.13% 57.46% 2.73 注：加粗字表示该体长组标准椭圆面积(SEAc)，右上角数据和左下角数据均表示两个体长组SEAc占左列体长组SEAc的百分比。 Note: Bold represents the standard elliptical area (SEAc) of the body length group. The data in the upper right corner and lower left corner both represent the percentage of the two body length groups' SEAc to the left column body length group's SEAc.
下载: 导出CSV
-
[1] ZHANG S, LI M, ZHU J F, et al. An integrated approach to determine the stock structure of spinyhead croaker Collichthys lucidus (Sciaenidae) in Chinese coastal waters[J]. Front Mar Sci, 2021, 8: 693954. doi: 10.3389/fmars.2021.693954
[2] 熊朋莉, 陈作志, 侯刚, 等. 珠江河口棘头梅童鱼生物学特征的年代际变化[J]. 南方水产科学, 2021, 17(6): 31-38. [3] 区又君, 廖锐, 李加儿, 等. 利用耳石日轮研究珠江口棘头梅童鱼的产卵期及生长[J]. 台湾海峡, 2012, 31(1): 85-88. [4] SONG N, YIN L N, SUN D R, et al. Fine-scale population structure of Collichtys lucidus populations inferred from microsatellite markers[J]. J Appl Ichthyol, 2019, 35(3): 709-718. doi: 10.1111/jai.13902
[5] 黄良敏, 谢仰杰, 李军, 等. 闽江口及附近海域棘头梅童鱼的生物学特征[J]. 集美大学学报(自然科学版), 2010, 15(4): 8-13. [6] 农业农村部渔业渔政管理局, 全国水产技术推广总站, 中国水产学会. 2022中国渔业统计年鉴[M]. 北京: 中国农业出版社, 2022: 59-64. [7] 农业部渔业渔政管理局. 2014中国渔业统计年鉴[M]. 北京: 中国农业出版社, 2014: 59-63. [8] WANG Y Y, ZHANG H A, CHEN Y W, et al. Trophic niche width and overlap of three benthic living fish species in Poyang Lake: a stable isotope approach[J]. Wetlands, 2019, 39(1): 17-23. doi: 10.1007/s13157-018-1028-3
[9] WANG J, JIANG R J, XIAO Y, et al. Trophic niche partitioning of five Sciaenidae species sampled in Zhoushan Archipelago waters via stable isotope analysis[J]. Front Mar Sci, 2022, 9: 880123. doi: 10.3389/fmars.2022.880123
[10] WANG S Q, WANG X F, XU L X, et al. Feeding habits and trophic niche of rainbow runner Elagatis bipinnulata in the western and central Pacific Ocean[J]. Environ Biol Fishes, 2022, 105(1): 139-149. doi: 10.1007/s10641-021-01200-w
[11] KVAAVIK C, OSKARSSON G J, PÉTURSDÓTTIR H, et al. New insight into trophic niche partitioning and diet of mackerel (Scomber scombrus) and herring (Clupea harengus) in Icelandic waters[J]. ICES J Mar Sci, 2021, 78(4): 1485-1499. doi: 10.1093/icesjms/fsaa100
[12] CHENG J, MA G Q, MIAO Z Q, et al. Complete mitochondrial genome sequence of the spinyhead croaker Collichthys lucidus (Perciformes, Sciaenidae) with phylogenetic considerations[J]. Mol Biol Rep, 2012, 39(4): 4249-4259. doi: 10.1007/s11033-011-1211-6
[13] MA Q Y, TIAN S Q, HAN D Y, et al. Growth and maturity heterogeneity of three croaker species in the East China Sea[J]. Reg Stud Mar Sci, 2021, 41: 101483.
[14] LIU H B, JIANG T, HUANG H H, et al. Estuarine dependency in Collichthys lucidus of the Yangtze River Estuary as revealed by the environmental signature of otolith strontium and calcium[J]. Environ Biol Fish, 2015, 98(1): 165-172. doi: 10.1007/s10641-014-0246-7
[15] 高世科, 黄金玲, 于雯雯, 等. 吕泗渔场两种石首鱼科鱼类营养生态学特征: 来自稳定同位素的证据[J]. 应用海洋学学报, 2021, 40(3): 413-420. [16] 贺舟挺, 薛利建, 金海卫. 东海北部近海棘头梅童鱼食性及营养级的探讨[J]. 海洋渔业, 2011, 33(3): 265-273. [17] 王建锋, 赵峰, 宋超, 等. 长江口棘头梅童鱼食物组成和摄食习性的季节变化[J]. 应用生态学报, 2016, 27(1): 291-298. [18] 杨纪明. 渤海鱼类的食性和营养级研究[J]. 现代渔业信息, 2001(10): 10-19. [19] 宋业晖, 薛莹, 徐宾铎, 等. 海州湾3种石首鱼的食物组成和生态位重叠[J]. 水产学报, 2020, 44(12): 2017-2027. [20] 王军, 苏永全, 柳建英, 等. 罗源湾五种石首鱼类的食性研究[J]. 厦门水产学院学报, 1994(2): 34-39. [21] 王静. 舟山群岛海域四种经济鱼类的摄食生态研究[D]. 舟山: 浙江海洋大学, 2022: 25-39. [22] POST D M. Using stable isotopes to estimate trophic position: models, methods, and assumptions[J]. Ecology, 2002, 83(3): 703-718. doi: 10.1890/0012-9658(2002)083[0703:USITET]2.0.CO;2
[23] LAYMAN C A, ARRINGTON D A, MONTAA C G, et al. Can stable isotope ratios provide for community-wide measures of trophic structure?[J]. Ecology, 2007, 88(1): 42-48. doi: 10.1890/0012-9658(2007)88[42:CSIRPF]2.0.CO;2
[24] JACKSON A L, INGER R, PARNELL A C, et al. Comparing isotopic niche widths among and within communities: SIBER-Stable Isotope Bayesian Ellipses in R[J]. J Anim Ecol, 2011, 80(3): 595-602. doi: 10.1111/j.1365-2656.2011.01806.x
[25] 石焱. 基于碳氮稳定同位素的闽江口常见鱼类营养生态位季节性变化[D]. 厦门: 集美大学, 2018: 26-33. [26] 叶学瑶, 任泷, 匡箴, 等. 基于稳定同位素技术的阳澄湖鱼类群落营养结构研究[J]. 中国水产科学, 2021, 28(6): 703-714. [27] 赖丽华, 张申增, 陆丽仪, 等. 2017—2020年珠江口棘头梅童鱼的种群特征[J]. 应用生态学报, 2022, 33(5): 1413-1419. [28] PAULY D, PALOMARES M L. Fishing down marine food web: it is far more pervasive than we thought[J]. Bull Mar Sci, 2005, 76(2): 197-211.
[29] JENNINGS S, GREENSTREET S P R, HILL L, et al. Long-term trends in the trophic structure of the North Sea fish community: evidence from stable-isotope analysis, size-spectra and community metrics[J]. Mar Biol, 2002, 141(6): 1085-1097. doi: 10.1007/s00227-002-0905-7
[30] 郑德锋, 赵金良, 周文玉, 等. 我国沿海棘头梅童鱼(Collichthys lucidus)群体遗传结构的AFLP分析[J]. 海洋与湖沼, 2011, 42(3): 443-447. [31] QUILLFELDT P, EKSCHMITT K, BRICKLE P, et al. Variability of higher trophic level stable isotope data in space and time: a case study in a marine ecosystem[J]. Rapid Commun Mass Spectrom, 2015, 29(7): 667-674. doi: 10.1002/rcm.7145
[32] 廖建基, 郑新庆, 杜建国, 等. 基于氮稳定同位素的九龙江口鱼类营养级研究[J]. 海洋学报, 2015, 37(2): 93-103. [33] 高春霞, 戴小杰, 田思泉, 等. 基于稳定同位素技术的浙江南部近海主要渔业生物营养级[J]. 中国水产科学, 2020, 27(4): 438-453. [34] 何雄波, 李波, 王锦溪, 等. 不同时期北部湾日本带鱼营养生态位差异[J]. 应用生态学报, 2021, 32(2): 683-690. [35] 银利强, 孔业富, 吴忠鑫, 等. 南海中西部海域春季三种金枪鱼类的营养生态位比较[J]. 生态学杂志, 2020, 39(12): 4121-4130. [36] 闫光松, 张涛, 赵峰, 等. 基于稳定同位素技术对长江口主要渔业生物营养级的研究[J]. 生态学杂志, 2016, 35(11): 3131-3136. [37] 王淼, 徐开达, 梁君. 杭州湾北部棘头梅童鱼繁殖群体生物学特征初步分析[J]. 上海海洋大学学报, 2018, 27(5): 781-788. [38] GUL G, DEMIREL N. Ontogenetic shift in diet and trophic role of Raja clavata inferred by stable isotopes and stomach content analysis in the Sea of Marmara[J]. J Fish Biol, 2022, 101(3): 560-572. doi: 10.1111/jfb.15123
[39] 操亮亮, 刘必林, 李建华. 基于稳定同位素技术的东南太平洋公海茎柔鱼摄食生态分析[J]. 大连海洋大学学报, 2022, 37(1): 120-128. [40] ZHOU F, WU H P, JIA S S, et al. Ontogenetic variation of trophic habitat for sympatric benthic octopods in East China Sea derived from isotopic analysis on beaks[J]. Fish Res, 2021, 238: 105902. doi: 10.1016/j.fishres.2021.105902
[41] 郭家彤, 王腾, 陈得仿, 等. 大亚湾黑棘鲷的摄食习性[J]. 中国水产科学, 2021, 28(8): 1041-1050. [42] BARNES C L, BEAUDREAU A H, YAMADA R N. The role of size in trophic niche separation between two groundfish predators in Alaskan waters[J]. Mar Coast Fish, 2021, 13(1): 69-84. doi: 10.1002/mcf2.10141
[43] PARK H J, KWAK J H, LEE Y J, et al. Trophic structures of two contrasting estuarine ecosystems with and without a dike on the temperate coast of Korea as determined by stable isotopes[J]. Estuar Coast, 2020, 43(3): 560-577. doi: 10.1007/s12237-019-00522-4
[44] 杨蕊, 韩东燕, 高春霞, 等. 浙江南部近海前肛鳗营养生态位变化研究: 基于稳定同位素技术[J]. 生态学报, 2022, 42(23): 9796-9807. [45] 黄佳兴, 龚玉艳, 徐姗楠, 等. 南海中西部海域鸢乌贼中型群和微型群的营养生态位[J]. 应用生态学报, 2019, 30(8): 2822-2828. [46] 殷宝法, 淮虎银, 张镱锂, 等. 可可西里地区藏羚羊、藏原羚和藏野驴的营养生态位[J]. 应用生态学报, 2007, 18(4): 766-770. doi: 10.3321/j.issn:1001-9332.2007.04.010 [47] 徐超, 王思凯, 赵峰, 等. 长江口水生动物食物网营养结构及其变化[J]. 水生生物学报, 2019, 43(1): 155-164. [48] 曾艳艺, 赖子尼, 杨婉玲, 等. 珠江河口渔业生物稳定同位素营养级分析[J]. 生态学杂志, 2018, 37(1): 194-202. [49] ZHANG Y L, ZHANG C L, XU B D, et al. Impacts of trophic interactions on the prediction of spatio-temporal distribution of mid-trophic level fishes[J]. Ecol Indic, 2022, 138(2/3): 108826.
[50] AYELEN T, FRANCO C, ANÍBAL G N, et al. Trophic niche partitioning of five skate species of genus Bathyraja in northern and central Patagonia, Argentina[J]. J Fish Biol, 2020, 97(3): 656-667. doi: 10.1111/jfb.14416
[51] 李振华, 徐开达, 蒋日进, 等. 东海中北部小眼绿鳍鱼的食物组成及摄食习性的体长变化[J]. 中国水产科学, 2011, 18(1): 185-193. -
期刊类型引用(10)
1. 宋广军，宋伦，刘印，王昆，李楠. 滤食性贝类对饵料微藻摄食选择性研究进展. 水产科学. 2025(02): 320-326 . 
百度学术
2. 张聪颖，薛瑞萍，彭瑞冰，韩庆喜，蒋霞敏. 4株微藻悬浮和固定培养生长和氮磷去除效果的比较. 生态科学. 2023(01): 76-82 . 百度学术
3. 阮烨烨，林雄平，翁晴，何子泓，林浩，张芷滢. 不同氮磷浓度对莱布新月藻生长与脂溶性物质含量的影响. 亚热带植物科学. 2023(01): 18-22 . 百度学术
4. 马东怡，郭永军，李永仁，石八军，张中华. 盐田虾蛤混养模式下硬壳蛤的生长、水质及浮游植物特征分析. 天津农林科技. 2023(06): 6-12 . 百度学术
5. 朱昔恩，林原有，项桂德，熊建华，陈田聪，陈晓汉，张彬. 盐度、温度和光照对广西沿海牟氏角毛藻生长的影响. 广东农业科学. 2022(01): 136-141 . 百度学术
6. 叶舟，姚有智. 水体表层沉积物对磷的吸收及释放研究进展. 河北环境工程学院学报. 2021(01): 60-64 . 百度学术
7. 赵东会，赵丽达，尤宏，秦松，王寅初，焦绪栋. 耐高温球等鞭金藻(Isochrysis galbana)藻种的选育与饵料特性初步分析. 海洋与湖沼. 2021(01): 206-212 . 百度学术
8. 叶舟，吕丹亚，姚有智，吴康，华飞，靳维维. 河流表层沉积物磷的生物吸收、吸附机制及影响因素. 西昌学院学报(自然科学版). 2021(01): 70-74 . 百度学术
9. 曾茹，李亚军，何金曼，黄晓晴，邓晓东. 清澜港近岸海湾春秋两季浮游植物群落结构特征及其与环境因子的关系. 安徽农业科学. 2021(10): 77-82+100 . 百度学术
10. 宋高飞，朱宇轩，Anila P.Ajayan，杨丽华，贾云璐，姜传奇，毕永红. 重金属污染区秋季浮游植物群落特征及其与环境因子的关系. 水生生物学报. 2021(05): 1082-1092 . 百度学术
其他类型引用(11)
计量
- 文章访问数: 117
- HTML全文浏览量: 18
- PDF下载量: 40
- 被引次数: 21
粤公网安备 44010502001741号
