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
-
大洋中上层鲨鱼多为顶级捕食者,可通过下行效应影响被捕食者的数量丰度和种群结构,对生态系统的稳定和物种多样性起关键作用[1]。鲨鱼具有生长缓慢、性成熟晚、繁殖率和种群密度低等特点[2-4],易受威胁从而导致数量减少,已成为保护生物学研究的重点关注对象。研究鲨鱼种间的食性差异、营养关系及资源利用模式有助于揭示其物种共存机制。
生物对其饵料资源实际和潜在的占据与利用部分即营养生态位,可反映物种间的营养关系及资源竞争状况[5],脂肪酸组成分析 (Fatty acid profiles analysis, FAs) 是研究大型海洋捕食者摄食的有力工具[6],摄食所得的外源性脂肪酸常以保守形态沉积于组织中,可反映捕食者一段时间内的摄食信息[7-8]。某些脂肪酸对生物的食性与营养关系具有特定的指示作用,例如:二十碳五烯酸 (C20:5n3) 的含量可指示以硅藻为基础的食物网,二十二碳酸和二十四碳酸 (C22:0+C24:0) 含量可追溯陆源植物食性[7]。因此,脂肪酸组成可用于反映鲨鱼与其食物间的动态关系和营养流动,进而指示其营养生态位[9]。不同组织的脂肪酸周转率存在差异,适用于探究不同时间尺度的摄食生态学信息[10-11]。肝脏是鲨鱼的一种能量储存器官,其脂肪酸组成随摄食变化周转较快,可反映鲨鱼短期内的摄食生态学信息;而鲨鱼肌肉组织脂肪酸随摄食变化周转较慢,可提供综合的长期摄食生态学信息[12]。Eduardo等[13]利用肝脏和肌肉脂肪酸组成的差异分析,探究了秘鲁北部3个海域锤头双髻鲨 (Sphyrna zygaena) 短期和长期营养生态位的差异。Xu等[10]也应用该技术探究了东太平洋同海域中上层鲨鱼的种间资源分配及利用情况,其中大眼长尾鲨 (Alopias superciliosus) 和镰状真鲨 (Carcharhinus falciformis) 的肝脏生态位重叠较高,而肌肉却不重叠,表明其短期内食性相似,而在较长时间尺度上存在摄食隔离。
目前,对于大西洋中上层鲨鱼摄食生态学的研究仍相对匮乏,不利于鲨鱼种群监测和资源保护。因此,本研究以热带大西洋大青鲨 (Prionace glauca)、长鳍鲭鲨 (Isurus paucus)、拟锥齿鲨 (Pseudocarcharias kamoharai) 和尖吻鲭鲨 (I. oxyrinchus) 为研究对象,通过分析其肌肉和肝脏脂肪酸组成,探讨其营养生态位的种间分化,以进一步掌握热带大西洋中上层鲨鱼的种间关系和生活史,并为相应的保护生物学研究提供参考。
1. 材料与方法
1.1 样品采集
4种鲨鱼样本于2016年5月及2017年5月在大西洋中部热带海域 (1°43'W-41°25'W, 12°33'N-8°88'S) 进行延绳钓渔业资源调查时兼捕所得,捕获点见图1,基础生物学信息见表1。鲨鱼肌肉和肝脏样本分别取自背鳍附近肌肉组织和肝脏任一叶前部,叉长在船上直接测定。将采集的样本真空包装后低温冷冻保存,运回实验室后放入 −80 ℃超低温冰箱保存,以便进行脂肪酸分析。
表 1 鲨鱼样本信息Table 1. Sample information of sharks
种名 Species样本 Sample 叉长 Fork length/mm 体质量 Body mass/kg 雄性
Male雌性
Female范围
Range$ \overline { x}\pm s$ 范围
Range$ \overline { x}\pm s$ 大青鲨
P. glauca5 9 790~2 360 1 994.30±374.90 6~80 53.50±18.57 长鳍鲭鲨
I. paucus2 5 1 100~1 970 1 522.70±271.70 13~68 37.80±18.82 拟锥齿鲨
P. kamoharai6 2 780~980 902.90±69.00 3~7 4.91±1.36 尖吻鲭鲨
I. oxyrinchus1 4 1 540~3 150 2 056.00±633.00 34~88 64.20±19.33 1.2 样品处理
肌肉样本去除表皮,肌肉和肝脏用超纯水冲洗干净,随后装入5 mL 冻存管,肝脏无需冷冻干燥和研磨。肌肉在冷冻干燥机 (Christ Alpha 1-4) 中−55 ℃干燥24 h,而后用冷冻混合球磨仪 (Mixer mill MM440) 磨碎成粉末 (直径5 µm)。粗脂肪提取参考Folch等[14],依据GB/T 17376-2008《动植物油脂 脂肪酸甲酯制备》三氟化硼法进行脂肪酸甲酯化。加入4 mL氢氧化钠-甲醇溶液 (0.5 mol·L−1) 复溶后,80 ℃水浴回流加热20 min;随后加入4 mL三氟化硼-甲醇溶液 (0.14 kg·L−1),水浴加热20 min;冷却至室温,加入4 ml正己烷,震荡30 s;再加入饱和氯化钠溶液定容,静置分层。取正己烷层过0.20 μm有机相滤膜后,移入配有聚四氟乙烯瓶盖的厚壁玻璃进样瓶中;后将正己烷层与50 mg·L−1的十九烷酸甲酯内标以1∶1 (体积比) 混合,移至进样瓶中待测。以对比37种脂肪酸混标的保留时间进行定性分析,并以十九烷酸甲酯标品为内标进行定量分析,脂肪酸采用气相色谱-质谱联用仪 (GC-MS 7890B, Agilent) 进行测定。
1.3 数据处理
每种脂肪酸以占各脂肪酸总含量的百分比表示,为避免可能的分析误差造成影响,在统计分析中剔除百分比小于0.5%的脂肪酸[15]。使用SPSS 26.0和Excel 2016软件对数据进行统计分析。先将两年样本进行时空差异分析,发现其无差异,故只进行种间差异分析。为检验鲨鱼种间、组织间的脂肪酸差异,利用Shapiro-Wilk检验分析脂肪酸含量数据是否符合正态分布并进行方差齐性检验,符合则采用单因素方差分析 (One-way ANOVA) 检验各脂肪酸含量是否存在显著的种间差异。若不符合,则使用Mann-Whitney U非参数检验。
鲨鱼脂肪酸组成测定结果以“平均值±标准差 ($\overline { x}\pm s $)”表示。依据Meyer等[16]研究结果,使用C16:0、C18:0、C16:1n7、C18:1n9、C22:1n9、C24:1n9、C18:2n6、C20:4n6、C20:5n3和C22:6n3共10种脂肪酸进行主成分分析,以检验鲨鱼不同组织的种间脂肪酸组成差异。并利用R语言中的SIAR 软件包计算标准椭圆面积 (SEAc) 和重叠百分比,以表征各鲨鱼物种的营养生态位[17]。使用Origin 2023和R 4.0.2软件绘图。
2. 结果
2.1 两种组织脂肪酸组成结果
4种鲨鱼肝脏中总脂含量均大于肌肉,2种组织中脂肪酸含量差异显著。4种鲨鱼的2种组织中各检测出35种脂肪酸,分别为15种饱和脂肪酸、9种单不饱和脂肪酸和11种多不饱和脂肪酸。C16:0、C18:0、C16:1n7、C18:1n9、C20:4n6、C20:5n3和C22:6n3在脂肪酸中占优势比例。
肌肉和肝脏中含量大于0.5%的脂肪酸各有14种 (表2、表3)。肌肉中的饱和脂肪酸含量介于35.55%~41.97%。其中,C16:0和C18:0含量最高,分别为17.07%~21.71%和8.57%~12.95%,大青鲨与拟锥齿鲨、长鳍鲭鲨的2种饱和脂肪酸的种间差异显著 (P<0.05)。单不饱和脂肪酸含量介于17.01%~27.40%,其中,C18:1n9含量较高 (9.88%~17.87%),长鳍鲭鲨肌肉中C18:1n9含量最高,拟锥齿鲨最低。长鳍鲭鲨与大青鲨、拟锥齿鲨、尖吻鲭鲨的C18:1n9含量种间差异显著 (P<0.05)。
表 2 鲨鱼肌肉脂肪酸组成Table 2. Fatty acid profiles of muscle tissue of sharks% 脂肪酸
Fatty acid大青鲨
P. glauca长鳍鲭鲨
I. paucus拟锥齿鲨
P. kamoharai尖吻鲭鲨
I. oxyrinchus肉豆蔻酸 C14:0 0.97±0.25 0.97±0.29 0.78 ±0.04 1.04±0.17 十五碳酸 C15:0 0.92±031a 0.44±0.16b 0.71±0.03ab 0.96±0.27a 棕榈酸 C16:0 17.07±3.73ac 21.71±2.17b 21.16±1.23b 17.61±4.12bc 十七碳酸 C17:0 1.10±0.33a 0.69±0.20b 0.88±0.05ab 1.19±0.32a 硬脂酸 C18:0 12.95±2.64a 8.57±2.25c 9.85±2.27bc 12.84±2.01ab 花生酸 C20:0 0.95±0.32a 0.40±0.17b 0.71±0.05ab 0.99±0.31a 二十一碳酸 C21:0 1.00±0.36a 0.39±0.20b 0.74±0.05ab 1.03±0.33a 二十二碳酸 C22:0 1.03±0.36a 0.40±0.20b 0.76±0.05ab 1.06±0.34a 二十三碳酸 C23:0 0.99±0.47a 0.42±0.21b 0.78±0.05ab 1.10±0.36a 二十四碳酸 C24:0 1.18±0.42ab 0.46±0.24c 0.75±0.33bc 1.22±0.39a 饱和脂肪酸 ∑SFA 40.99±2.02a 35.55±4.61b 39.22±1.64ab 41.97±4.15a 十六碳一烯酸 C16:1n7 1.89±1.22 4.70±4.74 1.80±0.81 1.70±0.38 十七碳一烯酸 C17:1n7 0.97±0.32a 0.55±0.16b 0.74±0.04ab 1.01±0.30a 油酸 C18:1n9 12.49±3.84b 17.87±6.30a 10.68±2.16b 9.88±1.90b 二十碳一烯酸 C20:1 1.72±0.60 2.27±1.32 1.39±0.21 1.63±0.19 芥酸 C22:1n9 1.70±0.78 0.93±0.73 0.98±0.30 0.95±0.62 二十四碳一烯酸 C24:1n9 1.11±0.37ab 0.48±0.19c 0.81±0.05bc 1.33±0.65a 单不饱和脂肪酸 ∑MUFA 21.33±5.14ab 27.40±11.10a 17.50±2.01b 17.01±1.96b 亚油酸 C18:2n6 2.34±1.36a 0.72±0.25b 1.39±0.83ab 1.70±1.51ab 二十碳二烯酸 C20:2 1.12±0.33a 0.60±0.15b 0.84±0.07ab 1.13±0.30a 十八碳三烯酸 C18:3n3 1.11±0.39a 0.45±0.22b 0.82±0.06ab 1.15±0.37a γ-亚麻酸 C18:3n6 0.90±0.48 0.38±0.20 0.61±0.27 0.83±0.57 二十碳三烯酸 C20:3n6 1.12±0.35a 0.44±0.20b 0.79±0.08ab 1.12±0.33a 花生四烯酸 C20:4n6 (ARA) 4.82±0.91a 4.18±1.18ab 3.61±0.85b 5.01±0.60a 二十碳五烯酸 C20:5n3 (EPA) 3.20±1.30a 0.82±0.32b 1.75±0.57b 1.70±0.32b 二十二碳二烯酸 C22:2n6 0.71±0.64ab 0.44±0.22b 0.69±0.31ab 1.14±0.37a 二十二碳六烯酸 C22:6n3 (DHA) 22.30±5.71b 28.98±6.56ab 32.70±4.30a 27.18±6.29ab 多不饱和脂肪酸 ∑PUFA 37.68±5.41a 37.05±8.29ab 43.28±3.22b 41.02±4.42ab 二十二碳酸和二十四碳酸 C22:0+C24:0 2.20±0.78a 0.86±0.44b 1.50±0.35c 2.27±0.73a 二十二碳六烯酸/二十碳五烯酸 DHA/EPA 8.31±5.21c 39.75±14.58a 20.94±8.44b 16.52±5.11bc 总脂 Total lipid/(mg·g−1) 129.84±61.97 181.50±92.48 151.89±55.95 190.20±40.39 注:同一行不同上标字母代表差异显著 (P<0.05)。 Note: Different superscript letters within the same row represent significant differences (P<0.05). 表 3 鲨鱼肝脏脂肪酸组成Table 3. Fatty acid profiles of liver tissue of sharks% 脂肪酸
Fatty acid大青鲨
P. glauca长鳍鲭鲨
I. paucus拟锥齿鲨
P. kamoharai尖吻鲭鲨
I. oxyrinchus肉豆蔻酸 C14:0 4.23±1.32a 1.30±0.45b 1.82±0.65b 5.18±1.82a 十五碳酸 C15:0 0.93±0.24b 0.66±0.18bc 0.51±0.07c 1.32±0.33a 棕榈酸 C16:0 20.48±3.52bc 22.55±3.09ab 25.75±2.72a 17.24±2.52c 十七碳酸 C17:0 1.22±0.35b 1.49±0.33ab 1.17±0.31b 1.89±0.40a 硬脂酸 C18:0 8.28±1.32b 11.86±1.59a 12.01±2.36a 12.38±1.32a 饱和脂肪酸 ∑SFA 34.46±5.69b 38.87±5.12ab 42.45±5.22a 39.17±1.94ab 十六碳一烯酸 C16:1n7 5.94±2.17b 3.55±1.04b 5.74±2.18b 11.22±3.85a 十七碳一烯酸 C17:1n7 0.49±0.20a 0.58±0.18b 0.65±0.07b 1.24±0.46b 油酸 C18:1n9 12.07±7.77b 17.98±3.42b 25.81±2.75a 14.39±6.13b 二十碳一烯酸 C20:1 2.95±2.69 1.24±1.56 1.00±1.47 1.97±2.53 芥酸 C22:1n9 3.91±3.30 1.76±1.13 1.22±0.88 2.78±2.94 二十四碳一烯酸 C24:1n9 0.61±0.55b 1.22±0.70ab 0.99±0.78ab 1.58±1.08a 单不饱和脂肪酸 ∑MUFA 26.16±7.15b 26.42±4.21b 35.57±3.74a 34.01±2.49a 亚油酸 C18:2n6 1.24±0.52b 1.03±0.33b 1.00±0.20b 2.07±0.67a 二十碳二烯酸 C20:2 0.90±0.28a 0.51±0.09b 0.49±0.13b 0.68±0.10ab 二十碳三烯酸 C20:3n3 0.86±1.12b 3.24±2.79a 0.62±0.76b 0.26±0.06b 花生四烯酸 C20:4n6 (ARA) 0.85±1.05b 3.62±2.45a 0.80±0.77b 0.27±0.05b 二十碳五烯酸 C20:5n3 (EPA) 6.01±2.77a 3.04±1.84ab 3.06±2.66ab 1.99±0.56b 二十二碳二烯酸 C22:2n6 2.73±1.37ab 1.29±0.85ab 4.35±5.17a 0.58±0.59b 二十二碳六烯酸 C22:6n3 (DHA) 26.21±4.47a 21.64±5.51a 11.31±6.15b 20.31±2.66a 多不饱和脂肪酸 ∑PUFA 39.38±6.68a 34.71±7.45ab 21.98±6.57bc 26.82±2.68c 二十二碳酸和二十四碳酸 C22:0+C24:0 0.19±0.11b 0.40±0.15a 0.38±0.10a 0.38±0.14a 二十二碳六烯酸/二十碳五烯酸 DHA/EPA 5.12±2.00 17.08±27.79 36.89±63.68 11.26±5.34 总脂 Total lipid/(mg·g−1) 332.10±99.59 295.95±66.59 321.03±83.73 273.33±25.39 注:同一行不同上标字母代表差异显著(P<0.05)。 Note: Different superscript letters within the same row represent significant differences (P<0.05). 多不饱和脂肪酸含量介于37.05%~43.28%,仅有大青鲨与拟锥齿鲨种间差异显著 (P<0.05)。其中,C20:4n6、C20:5n3和C22:6n3含量较高,分别介于3.61%~5.01%、0.82%~3.20%和22.30%~32.70%。C20:5n3、C20:4n6和C22:6n3分别在大青鲨肌肉、尖吻鲭鲨和拟锥齿鲨中含量最高,其中大青鲨与拟锥齿鲨中的C20:4n6、C20:5n3和C22:6n3含量种间差异显著 (P<0.05)。拟锥齿鲨与尖吻鲭鲨中C20:4n6含量种间差异显著 (P<0.05),大青鲨与长鳍鲭鲨、尖吻鲭鲨中C20:5n3含量种间差异显著 (P<0.05)。此外,C22:0+C24:0含量介于0.86%~2.27%,种间差异显著 (P<0.05),在长鳍鲭鲨中含量最低。鲨鱼肌肉中的DHA/EPA介于8.31~39.75,最大值和最小值分别出现在长鳍鲭鲨和大青鲨中。
肝脏中饱和脂肪酸含量介于34.46%~42.45%,其中C16:0和C18:0含量最高,分别介于17.24%~25.75%和8.28%~12.38%。大青鲨与拟锥齿鲨的2种饱和脂肪酸含量存在显著种间差异(P<0.05)。单不饱和脂肪酸含量介于26.16%~35.57%,其中C18:1n9含量较高 (12.07%~25.81%)。拟锥齿鲨肝脏中C18:1n9含量最高,大青鲨中最低,拟锥齿鲨与大青鲨、长鳍鲭鲨、尖吻鲭鲨中C18:1n9含量种间差异显著 (P<0.05)。
多不饱和脂肪酸含量介于21.98%~39.38%,大青鲨与拟锥齿鲨、尖吻鲭鲨,长鳍鲭鲨与尖吻鲭鲨种间差异显著 (P<0.05),其中,C20:5n3和C22:6n3含量较高,分别介于1.99%~6.01%和11.31%~26.21%。大青鲨肝脏中C20:5n3和C22:6n3含量最高,尖吻鲭鲨中C20:5n3含量最低,拟锥齿鲨中C22:6n3含量最低。大青鲨与尖吻鲭鲨中C20:5n3种间差异显著 (P<0.05),拟锥齿鲨与大青鲨、长鳍鲭鲨、尖吻鲭鲨中C22:6n3种间差异显著 (P<0.05)。此外,C22:0+C24:0含量介于0.19%~0.40%,大青鲨与长鳍鲭鲨、拟锥齿鲨、尖吻鲭鲨种间差异显著 (P<0.05),其中,大青鲨中含量最低。鲨鱼肝脏中DHA/EPA介于5.12~36.89,最大值和最小值分别出现在拟锥齿鲨和大青鲨中。
2.2 脂肪酸主成分分析
由主成分分析图可见造成差异的主要脂肪酸种类,肝脏相对于肌肉点的分布更分散 (图2)。尖吻鲭鲨与大青鲨的肌肉组织点基本散布在相同区域,表明其脂肪酸组成相似,两者肌肉中C20:5n3和C20:4n6含量较高;而拟锥齿鲨与长鳍鲭鲨的2种组织点均散布在相同区域,均具有较高的C18:1n9含量。2种组织点中,大青鲨与拟锥齿鲨和长鳍鲭鲨散布在不同区域,表明其脂肪酸组成存在差异,大青鲨的C18:2n6、C20:5n3和C22:1n9含量较高,其余2种鲨鱼的C16:0和C18:1n9含量较高。
2.3 营养生态位
营养生态位宽度及重叠见图3,尖吻鲭鲨肌肉组织中的生态位宽度较小 (1.15),而在肝脏中则较大 (4.71)。大青鲨2种组织营养生态位宽度均较大 (3.43~5.26),而拟锥齿鲨2种组织营养生态位宽度均相对较小 (0.47~2.71),长鳍鲭鲨肌肉生态位宽度较大 (5.37)。长鳍鲭鲨与拟锥齿鲨2种组织生态位图谱均存在一定重叠 (39.19%~42.13%)。在2种组织生态位图谱中,大青鲨、尖吻鲭鲨与长鳍鲭鲨、拟锥齿鲨均分离,而尖吻鲭鲨与大青鲨肌肉生态位图谱重叠度较高 (63.90%),肝脏中则完全分离。
3. 讨论
3.1 种间脂肪酸组成差异
脂肪酸是鱼类代谢的主要能量来源,是构成鱼体其他化合物的主要原料,影响鱼类的神经和视觉发育[18],可促进卵泡成熟,精子和类固醇的产生[19],因此鲨鱼不同组织对脂肪酸的吸收以及利用存在差异。4种鲨鱼肝脏中总脂含量大于肌肉,鲨鱼倾向于在肝脏中存储并合成脂质,之后通过血液、淋巴系统等将脂肪酸转移到其他组织[20]。与鲨鱼肌肉相比,肝脏中的单不饱和脂肪酸含量相对较高,而多不饱和脂肪酸含量较低,这可能由于食物中的多不饱和脂肪酸被优先输送到性腺和肌肉等组织中,供应生长和发育需求[21]。饱和脂肪酸和单不饱和脂肪酸可为机体提供能量,其中饱和脂肪酸是构成生物细胞膜的重要成分,由于肌肉中稳定的膜结构和功能需要,从而导致肌肉中的饱和脂肪酸含量高。2种组织的饱和脂肪酸中C16:0均含量最高,尤其在肝脏中;单不饱和脂肪酸中均以C18:1n9最为丰富,作为主要的供能物质,C18:1n9与生长、发育、基础代谢和运动关系密切[22]。
鲨鱼不能依靠自身从头合成C22:6n3、C20:5n3和C20:4n6,只能靠摄食来满足其生长发育所需的必需脂肪酸,可提供鲨鱼的摄食信息。C20:5n3主要由硅藻自身合成,且在头足类体内含量较高[23],大青鲨2种组织中含有相对较高水平的C20:5n3,表明大青鲨偏好捕食头足类,所处食物链底端初级生产者可能为硅藻,这与该种鲨鱼的胃含物分析结果一致[24]。
C20:4n6多由双壳类和甲壳类选择性的保留在体内[25],本研究中,长鳍鲭鲨2种组织中C20:4n6含量均较高,表明其倾向于捕食底栖甲壳类。长鳍鲭鲨个体较小且存在垂直洄游模式[26],可能捕食垂直迁移的甲壳类[27]。食物链传递过程中,DHA/EPA会随营养级的上升而升高[7],可指示生物营养级。本研究中大青鲨2种组织的DHA/EPA低于其他3种鲨鱼,表明其营养级相对较低,这可能与其较其他3种鲨鱼摄食区域范围更广[28]、食性更泛化,主要摄食中小型鱼类和头足类有关[29],与Mendez等[30]利用稳定同位素分析加利福尼亚西南海域大青鲨的营养级结果一致。深海物种一般C18:1n9含量丰富,可用来指示生物对深水高压环境的适应性[31],本研究中4种鲨鱼C18:1n9含量 (9.88%~25.81%) 高于浅海长尾鲨 (7.20%~13.00%) 与镰状真鲨 (8.80%~12.50%)[32],表明其栖息水层较深[23]。尖吻鲭鲨2种组织的C22:0+C24:0和C18:2n6含量相对较高,这两种脂肪酸可反映食物网主要碳源来自陆源植物[33]。本研究中尖吻鲭鲨多为未成熟幼鲨,尖吻幼鲨为获得足够的营养会洄游至大陆架等近岸海域[34-35],从近岸食物网获取营养物质。
3.2 鲨鱼种间营养生态位分化
营养生态位重叠反映了种间食物资源利用的相似程度及潜在竞争关系[36],同域共存的物种可通过生态位分化避免或减少对有限资源的竞争[37]。拟锥齿鲨和长鳍鲭鲨的肌肉和肝脏生态位均存在重叠,拟锥齿鲨体型较小,喜好捕食小型深海鱼类、底栖甲壳类与头足类[38],而本研究中长鳍鲭鲨均为未成熟个体,同样偏好捕食底栖甲壳类和小型硬骨鱼类[39],食性相似。由于大青鲨终身大洋性的生活史特征,摄食范围广泛且栖息水层相对较浅[28],与其他3种鲨鱼生态位分离,表明种间食物资源及栖息地利用的竞争程度较弱;与应用稳定同位素分析东北太平洋中部海域大青鲨与其他鲨鱼营养生态位分化的结果一致[40]。尖吻鲭鲨肝脏生态位差异较大且与其他3种鲨鱼完全分离,这可能是因其幼鲨短期内会频繁游向更深、温度更低的海域[41]的,食物组成与栖息环境较为复杂[11]。而尖吻鲭鲨肌肉生态位与大青鲨重叠较高可能是因为本研究中大青鲨数量相对尖吻鲭鲨较多所致,在较长时间尺度上两者存在栖息地和摄食竞争。Klarian等[42]胃含物分析结果显示,大西洋尖吻鲭鲨与大青鲨摄食相似,竞争激烈且有研究表明大青鲨主要以头足类、硬骨鱼类和甲壳类为食[24],而尖吻鲭鲨可捕食硬骨鱼类[43],两者有一定程度的重合。长鳍鲭鲨与尖吻鲭鲨虽同属,但营养生态位分离,与陈子昂等[44]应用稳定同位素分析探究其营养生态位的结果一致;这可能因为本研究中长鳍鲭鲨多为雄鲨且体型较小,而尖吻鲭鲨则多为雌鲨且体型较大所致,两种鲨鱼不同体型间存在食性差异[34,39],但性间差异仍有待探究。
营养生态位宽度可用于表征不同物种对食物资源的利用和栖息地竞争能力的强弱[36]。尖吻鲭鲨肌肉营养生态位宽度较小,可能是由于本研究中的尖吻鲭鲨多为未成熟个体,饵料相对单一,Maia等[45]利用胃含物分析发现尖吻鲭鲨幼鲨偏好捕食硬骨鱼类,其食物的均匀度与多样性[44]可能低于大青鲨和长鳍鲭鲨。而尖吻鲭鲨肝脏生态位宽度较大,可表明其短期内摄食的食物种类丰富和对资源的利用率较其他鲨鱼高。大青鲨肌肉和肝脏营养生态位宽度相对较大,表明其大洋性生活史特征及摄食饵料种类丰富。而拟锥齿鲨2种组织生态位相对较小,表明其对资源的利用能力不及其他3种鲨鱼。研究显示,体型较大的鲨鱼比较小的鲨鱼有更高的代谢需求,对栖息地利用范围更广,能更有效地捕获食物[46]。长鳍鲭鲨具有内温性,对环境的适应性较强[47],因此其肌肉生态位宽度较大。4种鲨鱼肝脏和肌肉生态位宽度结果存在差异,是由于肝脏内复杂的脂肪酸代谢过程[10],导致其脂肪酸组成变异系数较大。
4. 结论
本研究应用多组织脂肪酸组成对中上层鲨鱼在不同时间尺度上的种间摄食、栖息地利用及营养生态位关系变化进行了较为深入的研究。4种鲨鱼存在不同程度的食性或栖息地差异,大青鲨营养生态位宽度较大,表明其对环境的适应性较强,而拟锥齿鲨营养生态位宽度较小,表现为食性特化。鲨鱼在不同时期的资源分配模式存在差异,尖吻鲭鲨与大青鲨肌肉的营养生态位重叠最高,在较长时间尺度上存在摄食竞争。肝脏脂肪酸组成可反映鲨鱼短期摄食信息,而肌肉组织可提供其综合的长期摄食信息;而鱼类的个体大小、性别、成熟度、栖息水层、食物来源等均会影响其脂肪酸组成。为更全面地了解中上层鲨鱼的营养生态学,应将脂肪酸组成分析与稳定同位素分析、DNA条形码技术等多种方法结合,进一步加深对顶级捕食者同域共存关键机制的研究。
-
表 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). 表 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. 表 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. -
[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.