Analysis of changes in intestinal structure and microbial composition in Elentheronema tetradactylum juvenile at different days of age
-
摘要:
为研究四指马鲅 (Elentheronema tetradactylum) 幼鱼的肠道发育规律和菌群变化特征,为其健康养殖和饵料投喂管理提供科学依据,采集了养殖实验过程中不同日龄 (37、44、58、72、86和114日龄)的四指马鲅,并结合生长性能、肠道组织结构、肠道菌群组成变化等指标进行研究。结果显示:各日龄幼鱼间的生长性能存在显著性差异 (p<0.05)。随着鱼苗日龄增加,其肠道绒毛长度、宽度、肌层厚度和杯状细胞数量等逐渐增加,在58日龄时肠道结构相对完善。在幼鱼肠道中共鉴定出40个门和521个属的微生物,在门水平上,变形菌门、放线菌门和厚壁菌门为优势菌门,占所有肠道细菌的78.69%。在属水平上,不同日龄的优势菌属存在显著性差异 (p<0.05),丰度变化较大。伯克氏菌属(Burkholderia)、弧菌属(Vibrio)、罗尔斯通氏菌属(Ralstonia)与绒毛长度和肌层厚度均呈显著正相关性 (p<0.05),芽殖杆菌属(Gemmobacter)、硝酸盐还原污物球菌属 (Defluviimonas)与绒毛长度、宽度和肌层厚度均呈显著负相关性 (p<0.05)。转换颗粒饲料后,肠道内病原菌属的相对丰度上升,建议不宜在幼鱼早期投喂过多的颗粒饲料,且应加强换料期间的水质监管和病原菌防控。
Abstract:In order to explore the intestinal development and microbial communities of Elentheronema tetradactylum juvenile, and to provide a scientific basis for its healthy breeding and feeding management, we investigated the growth performance, intestinal structure and gut microbiota of the juveniles at different days of age (37, 44, 58, 72, 86, 114). The results show that there were significant differences in the growth performance among different days of age (p<0.05). As the age of the juveniles increased, the length and width of the villi, the thickness of the myenteric layer, and the number of cup-shaped cells increased gradually and became more completed. On Day 58, the morphological structure of the intestinal tract was relatively well-developed. A total of 40 microbial phylum and 521 microbial genus were identified in the intestine. Proteobacteria, Actinobacteria and Firmicutes were the dominant phylum, which accounted for 78.69% of the total intestinal microbiota. On genus level, there were significant differences in the relative abundances of microbial genus at different days of age (p<0.05), with large variations in abundance. Pearson correlation analysis reveals that the villus length, width, and muscular thickness were negatively correlated with Gemmobacter and Defluviimonas (p<0.05), but positively correlated with Burkholderia, Vibrio and Ralstonia (p<0.05). The relative abundances of pathogenic genus (Ralstonia, Vibrio) exhibited an increase after being fed with the pellet feed. It is not advisable to feed E. tetradactylum juvenile with an excessive amount of pellet feed at early stage, and it should enhance the water quality monitoring and pathogen prevention during the feed transition period.
-
四指马鲅 (Elentheronema tetradactylum) 隶属于辐鳍鱼纲、鲈形目、马鲅亚目、马鲅科、四指马鲅属,俗称午鱼、马友,是一种暖温、广盐性鱼类,主要分布于太平洋、印度洋西部以及澳大利亚西海岸,在中国东海、南海均有分布[1-2]。四指马鲅具有肉质鲜美、生长速度快、营养价值高等特点[3],是我国重要的海洋经济养殖鱼类[4]。近年来,随着四指马鲅人工繁育技术的突破,已在广东、广西、福建、海南、浙江等沿海地区推广养殖,现已发展成为水产养殖的新品种[5]。国内对于四指马鲅的研究主要集中在组织结构[6-7]、环境胁迫[5,8]、营养成分[9]及养殖技术[10]等方面。目前四指马鲅的养殖存在苗种培育阶段死亡率较高的问题,严重制约了其养殖生产规模的扩大。
肠道作为鱼类消化吸收的主要器官[11],其结构与功能对鱼类的生长、发育至关重要。幼鱼的肠道发育情况直接关系到鱼类对营养物质的吸收和利用效率[12],肠道内存在的免疫细胞可以帮助幼鱼抵抗病原体入侵[13],肠黏膜所形成的特殊构造如皱襞高度、绒毛高度及隐窝深度等是评估肠道发育及其消化吸收功能的重要组织形态学指标[14]。此外,肠道绒毛上的杯状细胞也是反映肠道消化能力的重要指标[15]。幼鱼在生长过程中,其肠道结构逐步完善,饵料来源也趋于多样化,明显影响到肠道菌群的变化[16]。幼鱼摄食不同的饵料也会对肠道结构发育和肠道菌群组成产生影响[17]。此外,健康的肠道菌群可促进鱼类生长和减少肠道炎症,肠道菌群与宿主之间存在着密切的相互作用关系,它们不仅参与了营养物质代谢和能量转化[18],还可通过激活免疫细胞调节宿主的免疫系统[19],增强疾病抵抗能力[20]。以往研究表明,多种鱼类的肠道菌群结构在生长阶段处于动态平衡状态,3和6月龄花鲈 (Lateolabrax japonicus) 肠道微生物中的鲸杆菌属 (Cetobacterium) 相对丰度较高,9月龄的梭状芽孢杆菌属 (Clostridium) 相对丰度较高[13];大口黑鲈 (Micropterus salmoides) 肠道微生物中,气单胞菌属 (Aeromonas) 在第5天占据优势,而第25天假单胞菌属 (Pseudomonas)、乳球菌属 (Lactococcus)和不动杆菌属 (Acinetobacter) 相对丰度增加;第45天假单胞菌属和芽孢杆菌属(Bacillus)相对丰度较高,第180天成鱼肠道中检测出以邻单胞菌属 (Plssiomonas) 为代表的新优势菌群[21]。因此,研究四指马鲅不同生长发育阶段肠道结构形态及肠道菌群组成,可掌握其肠道对外源饵料的适应性及肠道菌群的变化规律,对其幼鱼期间的饵料投喂策略和病害防治具有重要意义。
本研究选用四指马鲅幼鱼进行养殖生长实验,结合肠道组织苏木素-伊红 (HE)染色切片观察、肠道菌群检测,比较分析了四指马鲅幼鱼在不同阶段的生长性能、肠道组织发育和结构形态变化以及肠道菌群结构;旨在揭示四指马鲅在幼鱼阶段的摄食生长与肠道发育、菌群变化之间的内在联系,为制定适宜的饵料投喂制度、病害防控等日常管理提供理论支撑。
1. 材料与方法
1.1 实验材料与管理
在广东省中山市四指马鲅养殖基地开展实验。2023年5月,采用同一批次的鱼苗6 000尾开展实验,平均全长为 (1.26±0.07) cm,平均体质量为 (1.5±0.1)×10−3 g。在室内育苗池 (4 m×5 m×1.2 m) 中暂养7 d,水体盐度为10‰,pH为7.6~8.0,溶解氧 (DO) 质量浓度>5 mg·L−1。暂养期间,投喂鳗鱼粉 (质量分数:蛋白质≥45%,脂肪≥10%),逐渐降低盐度至3‰。暂养结束后,转移到池塘内的3口网箱 (长3 m×宽4 m×高1.7 m) 中,养殖密度为100尾·m−3,池塘水体盐度为3‰~5‰,pH为7.5~8.1,养殖周期共12周。实验开始时,鱼为37日龄,在37~44日龄投喂鳗鱼粉,44~58日龄投喂粒径为0.6 mm的饲料,58~72日龄和72~86日龄均投喂粒径为1 mm的饲料,86~114日龄投喂粒径为2 mm的饲料,饲料粗蛋白质量分数为48%,粗脂肪质量分数为5%,粗灰分质量分数为6%。每天饱食投喂2次,日投喂量为体质量的5%。
1.2 肠道组织学分析
分别在37、44、58、72、86和114日龄采样,从每个网箱中随机挑选5尾鱼 (采样前停食1 d),记录全长、体长、体质量,计算特定生长率 (SGR)、平均日增质量 (ADG) 和肥满度 (CF);取3尾鱼中肠1 cm置于4% (w)多聚甲醛溶液固定,经过乙醇脱水、二甲苯透明、石蜡包埋、切片、HE染色和中性树脂封片,采用光学显微镜 (ZEISS Axio ScopeA1型)进行观察和记录图像。使用Image J在每个切片选取5根肠道绒毛,测量绒毛长度、宽度、肌层厚度,计算肠道结构变化情况。
1.3 肠道微生物16S rDNA基因高通量测序与分析
取3尾鱼,用75% (φ) 乙醇消毒体表,用无菌镊子取肠道组织,置于2 mL冻存管内并于 −80 ℃保存,以测定肠道菌群组成。采用HiPure Stool DNA mini Kit试剂盒 (Magen,中国)提取肠道内容物中微生物的总基因组DNA。DNA浓度使用NanoDrop 2000检测。使用正向引物338F (5'-ACTCCTACGAGGCAGCA-3') 和反向引物806R (5'-GGACTACHVGGGTWTCTAAT-3'),对16S rDNA基因的可变V3—V4区域进行扩增:94 ℃ 5 min;94 ℃ 30 s,53 ℃ 30 s,72 ℃ 30 s,共29个循环;72 ℃ 8 min。扩增后的PCR产物通过AMPure XP Beads进行纯化,并经过建库检测,将符合条件的样本文库使用Illumina NovaSeq平台对测序样本进行双端测序。
使用QIIME2 (V 2020.2)软件[22]对原始数据进行切除Barcods和引物、质量过滤、去噪和去嵌合体,将相似性100%的序列归为扩增子测序变体 (Amplicon sequence variants, ASV),得到包含细菌特征及相对丰度的表格和细菌代表性序列文件。使用Silva123数据库对代表性序列进行注释确定物种分类信息。通过QIIME命令summarize_taxa_through_plots.py生成各分类水平的细菌群落组成信息。通过QIIME命令alpha_diversity.py和collate_alpha.py 计算细菌α-多样性 (香农多样性指数Shannon和物种丰富度指数Richness)。通过QIIME命令beta_ diversity_through_plots.py计算细菌群落的β-多样性,利用基于Bray-Curtis距离的主坐标分析 (Principal coordinates analysis, PCoA) 表征不同组间肠道菌群结构的差异。使用线性判别分析LEfSe (Linear discriminant analysis effect Size)[23]确定不同日龄幼鱼肠道细菌组成的显著性差异,阈值对数评分为4.0。使用SPSS 22.0软件对肠道菌群属水平与肠道组织结构进行Pearson相关性分析,采用Origin 2024软件作图。
1.4 统计与分析
四指马鲅幼鱼的生长性能指标包括特定生长率 (SGR)、平均日增质量 (ADG) 和肥满度 (CF),计算公式如下:
$$ \mathrm{SGR=(ln} {W} _{ {t} } \mathrm{−ln} {W} _{ \mathrm{0}} \mathrm{)/} {t} \mathrm{\times 100\% } $$ (1) $$ \mathrm{ADG=(}{W}_{{t}}\mathrm{−}{W}_{\mathrm{0}}\mathrm{)/}{t} $$ (2) $$ \mathrm{CF=} {W} _{ {t} } \mathrm{/} {L} ^{ \mathrm{3}} $$ (3) 式中:W0为初始体质量 (g);Wt为终末体质量 (g);t为试验周期 (d);L为平均体长 (cm)。
实验数据均以“平均值±标准误 ($\bar x $±$ s_{\bar x}$, n=5)”表示,通过SPSS 20.0软件进行单因素方差分析 (One-way ANOVA),通过Tukey's 检验方法进行多重比较。p<0.05表示具有显著性差异。
2. 结果
2.1 四指马鲅生长性能测定
实验期间,不同日龄的四指马鲅在特定生长率、平均日增质量、肥满度上均存在显著性差异 (p<0.05,表1)。其中,44日龄的特定生长率显著高于其他日龄 (p<0.05),37日龄的显著高于58和72日龄,其余各日龄之间差异不显著;86日龄的平均日增质量显著高于其他日龄 (p<0.05),其余日龄差异不显著;58~114日龄的肥满度随幼鱼生长而逐渐升高,且72、86和114日龄的肥满度显著高于58日龄 (p<0.05)。37与44日龄的肥满度显著高于其他日龄 (p<0.05),原因可能是其消化道未发育完全,存在卵黄囊,导致体型较为肥满。
表 1 不同日龄四指马鲅生长性能Table 1. Growth performance of E. tetradactylum at different days of age指标
Indicator组别 Group A B C D E F 体长 Body length/cm 2.37±0.25 3.47±0.21 6.20±.032 7.00±0.30 10.08±0.36 13.20±1.19 体质量 Body mass/g 0.24±0.02 0.71±0.07 1.60±0.34 3.46±0.65 8.62±1.21 44.72±11.37 特定生长率SGR/(%·d−1) 8.58±0.23b 15.36±1.73a 5.72±1.45c 5.38±1.62c 6.43±1.16bc 6.52±2.02bc 平均日增质量ADG/(g·d−1) 0.01±0.00a 0.07±0.01a 0.06±0.02a 0.13±0.05a 0.52±0.29b 1.15±0.22a 肥满度CF/(g·cm−3) 2.88±0.95a 2.47±0.32a 0.85±0.03c 1.44±0.17b 1.67±0.09b 1.77±0.42b 注:同行中不同小写字母间存在显著性差异 (p<0.05)。A、B、C、D、E、F分别表示37、44、58、72、86、114日龄的幼鱼。 Note: Values with different lowercase letters within the same row are significantly different (p<0.05). A, B, C, D, E and F represent the E. tetradactylum juveniles of 37, 44, 58, 72, 86 and 114 days of age, respectively. 2.2 四指马鲅肠道组织结构变化
随着养殖时间的增加,四指马鲅肠道绒毛的长度、宽度和肌层厚度不断增加,且各日龄间存在显著性差异 (p<0.05),绒毛数量和杯状细胞数量也有所增加 (表2,图1)。绒毛长度由91.69 μm增长到312.78 μm,绒毛宽度由19.39 μm增长到42.76 μm,肌层厚度由12.31 μm增加到56.72 μm。在44日龄时肠道壁结构界限模糊,杯状细胞数量稀少,每个绒毛上分布1~4个,主要分布于肠道绒毛前端或近前端;58日龄时,肠道组织结构清晰,杯状细胞增多,每个绒毛上分布5~10个。72、86和114日龄时,肠道组织结构清晰,黏膜层、黏膜下层、肌层和浆膜结构明显,杯状细胞数量显著增多,每个绒毛分布11~34个,且均匀分布在绒毛的前、中、后部。
表 2 不同日龄四指马鲅肠道组织结构变化Table 2. Intestinal tissue structure of E. tetradactylum at different days of age指标
Indicator组别 Group B C D E F 绒毛长度 Villus length/μm 91.69±8.42d 124.33±22.00c 168.11±42.35b 189.60±25.91b 312.78±45.92a 绒毛宽度 Villus width/μm 19.39±4.39c 20.14±3.29c 28.40±4.33b 38.28±8.84a 42.76±10.51a 肌层厚度 Muscular thickness/μm 12.31±1.57d 16.57±6.22bc 20.03±5.86c 31.58±5.07b 56.72±15.13a 注:同行中不同小写字母间存在显著性差异 (p<0.05)。B、C、D、E、F分别表示44、58、72、86、114日龄的幼鱼。 Note: Values with different lowercase letters within the same row are significantly different (p<0.05). B, C, D, E and F represent the E. tetradactylum juveniles of 44, 58, 72, 86 and 114 days of age, respectively. 2.3 肠道菌群结构变化
2.3.1 Alpha多样性指数
本研究共获得原始序列数2 534 257条,经质控和过滤后共得到有效序列2 388 796条,序列长度为236~249 bp,平均长度243 bp。按照97%的序列相似性聚类为7 543个OTUs,属于39个细菌门类、73个纲、128个目、235个科和520个属。不同日龄四指马鲅肠道内共有OTUs为228个,占总数的2.67%,37、44、58、72和86日龄的特有OTUs分别为1 567、1 330、1 354、1 073、1 005和986个 (图2)。随着日龄的增加,肠道内特有OTUs总数逐渐降低。不同日龄四指马鲅鱼苗的肠道菌群Alpha多样性情况见表3。不同日龄的Chao 1和ACE指数差异均不显著,在37日龄具有最大值。从Shannon指数来看,37日龄群落多样性最高,显著高于86日龄 (p<0.05)。不同日龄间的Simpson指数无显著性差异 (p>0.05)。
图 2 不同日龄四指马鲅肠道内细菌韦恩图分析注:每个椭圆代表一个分组,不同椭圆交集部分为不同分组间的共有物种。其中A、B、C、D、E、F分别表示37、44、58、72、86、114日龄的幼鱼。Figure 2. VENN analysis of intestinal microbiota in E. tetradactylum at different days of ageNote: Each ellipse represents a member group, and the intersection of different ellipses represents the common species between different groups. A, B, C, D, E and F represent the E. tetradactylum juveniles of 37, 44, 58, 72, 86 and 114 days of age, respectively.表 3 不同日龄四指马鲅鱼肠道菌群Alpha多样性统计Table 3. Alpha diversity of intestinal microbiota in E. tetradactylum at different days of age指标
Indicator组别 Group A B C D E F Chao 1 指数 Chao 1 index 1 019.70±247.37 919.92±24.98 795.38±167.93 818.00±32.87 739.96±231.98 805.94±80.17 ACE 指数 ACE index 986.88±211.51 938.39±29.54 810.60±163.07 817.51±20.64 774.20±258.03 866.18±93.10 香农指数Shannon index 4.65±0.01a 4.21±0.08ab 4.72±0.61a 3.86±0.47ab 3.16±1.58b 4.52±0.37ab 辛普森指数Simpson index 0.97±0.01 0.95±0.01 0.96±0.02 0.93±0.04 0.94±0.02 0.97±0.02 注:同行中不同小写字母间存在显著性差异 (p<0.05)。A、B、C、D、E、F分别表示37、44、58、72、86、114日龄的幼鱼。 Note: Values with different lowercase letters within the same row are significantly different (p<0.05). A, B, C, D, E and F represent the E. tetradactylum juveniles of 37, 44, 58, 72, 86 and 114 days of age, respectively. 2.3.2 基于门水平肠道菌群分析
在门水平上,四指马鲅肠道细菌种类及丰度从高到低 (取丰度前5类细菌)依次为变形菌门、放线菌门、厚壁菌门、梭杆菌门和蓝细菌门 (图3)。在幼鱼生长过程中,变形菌门为绝对优势菌群,其相对丰度在86日龄阶段有所下降,随后又明显回升。在多数采样阶段,放线菌门和厚壁菌门也是主要的细菌门类,变形菌门、放线菌门和厚壁菌门占肠道细菌总丰度的75%以上。在86日龄时,四指马鲅肠道内的梭杆菌门和蓝细菌门明显增加,伴随着变形菌门相对丰度的显著降低 (p<0.05)。
2.3.3 基于属水平的肠道菌群分析
在属水平上,四指马鲅肠道细菌种类及丰度按从高到低 (取丰度前5类细菌) 顺序为伯克氏菌属 (Burkholderia)、鲸杆菌属 (Cetobacterium)、申氏杆菌属 (Shinella)、芽殖杆菌属 (Gemmobacter) 和硝酸盐还原污物球菌属 (Defluviimonas) (图4)。其中,37日龄的优势属为芽殖杆菌属、副球菌属 (Paracocccus) 和红细菌属 (Rhodobacter);44日龄为硝酸盐还原污物球菌属、芽殖杆菌属和申氏杆菌属;58日龄为伯克氏菌属、弧菌属 (Vibrio) 和罗尔斯通氏菌属 (Ralstonia);72日龄为申氏杆菌属、博斯氏菌属 (Bosea) 和红细菌属;86日龄为鲸杆菌属、副球菌属和伯克氏菌属;114日龄为伯克氏菌属、罗尔斯通氏菌属和弧菌属。其中,硝酸盐还原污物球菌属、鲸杆菌属和伯克氏菌属在44、58、86和114日龄显著增加,而申氏杆菌、红细菌属和博斯氏菌属在72日龄显著增加;副球菌属和伯克氏菌属在58和72日龄显著减少,而红细菌属在58和86日龄显著减少。另外,弧菌属的相对丰度与罗尔斯通氏菌属丰度呈正相关,红细菌属的相对丰度与弧菌属呈负相关。
2.3.4 不同日龄四指马鲅肠道内菌群变化
采用LEfSe分析方法比较了不同日龄四指马鲅肠道中所有细菌的相对丰度差异性,在6个阶段 (A—F) 中有38个细菌分类群的相对丰度存在显著性差异 (图5)。在属水平上表现出显著差异的种类有16个:37日龄时的Nesiotabacter、Nautella、聚球藻属(Synechococcus)和杨氏菌属 (Yangia),44日龄时的硝酸盐还原污物球菌属、芽殖杆菌属和红球菌属的相对丰度明显增加,58日龄时的不动杆菌属和弧菌属,72日龄时的申氏菌属、红细菌属和剑菌属 (Ensifer)的相对丰度明显增加,86日龄时的鲸杆菌属和斯塔普氏菌属,14日龄时的伯克霍尔德氏菌属和罗尔斯通属的相对丰度明显增加。
在四指马鲅幼鱼肠道中变形菌门所包含的特征微生物范围最广,其次为蓝细菌门和梭杆菌门 (图6)。变形菌门中的α-变形菌纲、γ-变形菌纲和β-变形菌纲的丰度最高。进一步分析发现,α-变形菌纲中的红杆菌科在37日龄时显著高于其他阶段,但其下的属在所有日龄中均表现出显著性表现;莫拉菌科和弧菌科在58日龄时的丰度最高,根瘤菌科在72日龄时最高,梭杆菌科在86日龄时最高,伯克氏菌科在114日龄时最高。
图 6 物种分类学分支图注:Cladogram由内到外辐射的圆圈分别对应界、门、纲、目、科、属不同的分类层级,层级间的连线代表所属关系。每个圆圈节点代表1个物种,黄色节点代表分组间差异不显著,差异物种跟随组进行着色。有颜色的扇形区域标注了特征微生物的下属分类区间。A、B、C、D、E、F 分别表示37、44、58、72、86、114 日龄的幼鱼。Figure 6. Cladogram of speciesNote: Cladogram consists of circles radiating from the inside out, each corresponding to different levels of classification such as kingdom, phylum, class, order, family, and genus. The lines connecting the circles represent their relationships. Each circle node represents a species, with yellow nodes indicating insignificant differences between groups, while significantly different species are colored according to their group. Colored sector regions represent the subclassification intervals of characteristic microorganisms.A, B, C, D, E and F represent the E. tetradactylum juveniles of 37, 44, 58, 72, 86 and 114 days of age, respectively.2.3.5 肠道菌群和肠道组织结构相关性分析
对四指马鲅幼鱼肠道组织结构和肠道菌属的相对丰度进行相关性分析发现 (图7),伯克氏菌属、弧菌属和罗尔斯通氏菌属与绒毛长度、肌层厚度呈显著正相关性 (p<0.05),芽殖杆菌属和硝酸盐还原污物球菌属与绒毛长度、宽度和肌层厚度呈显著负相关性 (p<0.05)。
图 7 相关性分析热图注:色块的颜色表示相关系数,其中红色代表正相关,蓝色代表负相关;颜色深浅代表相关系数大小 (*. p<0.05)。Figure 7. Heat map of correlation analysisNote: The colors of the color blocks represent the correlation coefficients. Red represents positive correlation and blue represents negative correlation; the color shades represent the size of the correlation coefficients (*. p<0.05).3. 讨论
已有研究表明稚鱼、幼鱼阶段的消化道发育尚不完善,消化能力较弱,进食频率过高或过低都可能造成肠道受损,进而导致免疫力下降甚至感染疾病[24]。随着稚鱼成长为成鱼,其肠道内的组织、结构和功能逐渐发育成熟[25]。一般情况下,海水鱼类的仔稚鱼期在30日龄内[26-28],已有研究发现四指马鲅幼鱼43日龄已具完整的黏膜褶皱,消化系统结构与功能的发育已能基本满足机体的生长代谢需要[7]。肠道绒毛数量越多,长度和宽度越大,肠道吸收面积越大,机体的消化吸收能力越强[29]。杯状细胞可分泌黏液和消化酶,黏液有保护黏膜和润滑作用,消化酶可提高消化能力[30],杯状细胞数量的增多,也反映了消化能力的增强[31]。卵形鲳鲹 (Trachinotus ovatus) 肠道杯状细胞最早于孵化后的第3天出现在前肠,且数量随着仔鱼的生长而增加[32]。在本研究中,通过肠道组织切片观察发现37日龄四指马鲅肠道结构基本成型,但是绒毛数量、杯状细胞仍较少,肠道壁与黏膜层界限模糊,表明37~44日龄的四指马鲅仍处于稚鱼期,可投喂少量粉状饲料,但投喂量不宜过多。在58~114日龄,四指马鲅肠道结构发育明显,与成鱼基本一致,肠道绒毛长度、宽度逐渐增加,肠道绒毛数量和杯状细胞数量明显增加,表明其消化能力不断增强,因此在58日龄后逐渐增加人工饵料的投喂量可能更符合四指马鲅幼鱼的消化机能发育规律。此外,通过肠道属水平细菌与肠道结构进行Pearson相关性分析发现,伯克氏菌属和罗尔斯通氏菌属的相对丰度随着绒毛长度和肌层厚度的增加而升高,芽殖杆菌属和硝酸盐还原污物球菌属的相对丰度随着绒毛长度、宽度和肌层厚度的增加而降低。已有研究表明,伯克氏菌属和罗尔斯通氏菌属广泛存在于动物肠道和水环境中[33],随着幼鱼的生长发育,肠道绒毛增长,细菌的相对丰度也会增加。以上结果表明,随着鱼类的不断生长,其肠道细菌多样性在属水平上具有较大差异,且肠道内优势菌属与肠道结构发育存在相关性。优势菌属动态变化与组织结构相互联系,使得水生动物能够更好地获得生长所需的营养。
鱼类不同生长发育阶段的肠道菌群处于动态变化状态。肠道菌群易受生物和非生物因素影响,如水环境理化指标[34]、发育阶段[13]、个体差异[35]、不同区域[36]和生活习性[37]等因素均会影响肠道的菌群组成。肠道菌群在自然选择下与宿主形成寄生或共生关系[38],主要参与消化分解、营养代谢并影响宿主免疫系统[39]。已有研究表明,海水鱼肠道内的细菌门类主要以变形菌门、厚壁菌门和拟杆菌门等为主[37]。在本研究中,四指马鲅幼鱼肠道内的优势菌门为变形菌门、放线菌门、厚壁菌门、梭杆菌门和蓝细菌门,与虹鳟 (Oncorhynchus mykiss)[40]、赤点石斑鱼 (Epinephelus akaara)[41]、军曹鱼 (Rachycentron canadum)[42]肠道内的优势菌门相近。研究表明,大口黑鲈肠道内的气单胞菌属、不动杆菌属、肠杆菌属分别在第5、第25和第45天成为优势菌属,在第180 天时增加了多个新的优势菌属,例如邻单胞菌属、大肠埃希菌属-志贺氏菌属和鲸杆菌属等;随着幼鱼生长和饲料转换,七彩神仙鱼 (Symphysodon haraldi) 肠道内的内优势菌属由鲸杆菌属、邻单胞菌属和未识别的紫单胞菌科转化为鲸杆菌属、未分类消化链球菌科和Breznakia;不同月龄间花鲈的肠道微生物群落结构也表现出显著性差异,3和6月龄以鲸杆菌属、邻单胞菌属和罗姆布茨菌属为主,9月龄以梭菌属为优势菌属。鲸杆菌属作为一种共生菌,普遍存在于多种鱼类肠道中,具有调节肠道菌群组成、改善肝脏健康的功能[43]。在集约化养殖水库中的6种健康鱼类肠道中均检测到鲸杆菌属,该菌在杂交条纹鲈 (Morone saxatilis)、欧洲鲈 (Dicentrarchus labrax)、美国红鱼 (Sciaenops ocellatus)、奥尼罗非鱼 (Oreochromis)、平头鲻 (Mugil cephalus) 中均属于优势菌属[44]。本研究中,鲸杆菌属为四指马鲅在86日龄时的主要优势菌属,而在其余阶段丰度较低,这可能与其消化能力增强有关,也可能与宿主的健康状况密切相关。有研究表明,鲸杆菌属是潜在的益生菌,可产生维生素B12和抗菌肽,促进营养吸收,提高机体抗病力[43,45]。此外,随着鱼类生长过程中肠道发育的逐渐完善,其肠道菌群变得更加复杂,不同时期的菌属具有较大差异,但菌群多样性却呈下降趋势[35,46];与本研究结果一致,即随着幼鱼生长,肠道菌群的Alpha多样性指数呈降低趋势,其肠道微生物在不同时期处于动态平衡状态,与大口黑鲈[21]、七彩神仙鱼[47]、花鲈[13]等肠道菌群的变化情况相近。
由于幼鱼的肠道发育尚未完善,饵料来源和转换较为频繁,易造成菌群组成不稳定[48]。有研究表明,不同开口饵料、转食过程及饲料种类均会对鱼类肠道菌群造成影响。本研究中,对37~44日龄的四指马鲅幼鱼投喂鳗鱼粉,其余日龄投喂同一品牌的营养成分一致但粒径大小不同的饲料。在幼鱼转食颗粒饲料后的58日龄,肠道菌群的多样性增加、丰富度降低,与大口黑鲈转食饲料[16]前后肠道菌群的Alpha多样性指数变化一致;而在转食颗粒饲料后,在门水平上,优势菌门未发生改变,说明变形菌门、厚壁菌门和放线菌门是维系四指马鲅肠道菌群稳态的主要门类;而在属水平上,弧菌属和罗尔斯通氏菌属明显增加,在58和114日龄时为优势菌属。罗尔斯通氏菌属和弧菌属中的大多数菌被认为是水产养殖中的机会致病菌[49],危害水产动物健康。本研究中,肠道罗尔斯通氏菌属和弧菌属的相对丰度在转食颗粒饲料后上升,在养殖中应加强对这一阶段的防控措施。添加益生菌、中草药制剂可以有效抑制弧菌属的生长和繁殖[50]。
水环境与鱼类的健康息息相关,良好的水环境不仅能为其提供适宜的生存环境,还能减少疾病发生。与成鱼相比,幼鱼肠道微生物更易受水环境影响[35]。伯克氏菌属[51]、申氏菌属[52]、红细菌属[53]和博斯氏菌属[54]均属于反硝化菌,广泛存在于水环境和土壤中,可促进养殖环境的氮循环,在牙鲆 (Paralichthys olivaceus) 幼鱼肠道中亦有存在[55]。本研究中,幼鱼肠道中的伯克氏菌属、申氏菌属、红细菌属和博斯氏菌属丰度在44和72日龄时出现显著变化,且红细菌属的相对丰度较高,弧菌的相对丰度较低,二者可能存在一定程度的相互抑制作用。已有研究表明,光合细菌可以抑制弧菌[56],光合细菌与弧菌有一定的拮抗作用,通过添加光合菌制剂可使弧菌数量最多下降23.5%[57]。本研究中,不同日龄的四指马鲅幼鱼肠道中均存在不同比例的弧菌属,但均未发生明显病害,可能是因为弧菌属于原籍菌群或具有条件致病性[58]。
4. 结论
58日龄的四指马鲅肠道结构发育较为完善,在此阶段适合增加人工颗粒饵料投喂。变形菌门、放线菌门和厚壁菌门是四指马鲅幼鱼肠道中的优势菌门,不同日龄间的优势菌属具有显著性差异,且肠道内优势菌属的相对丰度与肠道结构发育状况具有相关性。在四指马鲅的生长过程中,由鳗鱼粉转食颗粒饲料阶段,有害菌丰度增加,应格外注意饵料转换时期的病害防治。
-
图 2 不同日龄四指马鲅肠道内细菌韦恩图分析
注:每个椭圆代表一个分组,不同椭圆交集部分为不同分组间的共有物种。其中A、B、C、D、E、F分别表示37、44、58、72、86、114日龄的幼鱼。
Figure 2. VENN analysis of intestinal microbiota in E. tetradactylum at different days of age
Note: Each ellipse represents a member group, and the intersection of different ellipses represents the common species between different groups. A, B, C, D, E and F represent the E. tetradactylum juveniles of 37, 44, 58, 72, 86 and 114 days of age, respectively.
图 6 物种分类学分支图
注:Cladogram由内到外辐射的圆圈分别对应界、门、纲、目、科、属不同的分类层级,层级间的连线代表所属关系。每个圆圈节点代表1个物种,黄色节点代表分组间差异不显著,差异物种跟随组进行着色。有颜色的扇形区域标注了特征微生物的下属分类区间。A、B、C、D、E、F 分别表示37、44、58、72、86、114 日龄的幼鱼。
Figure 6. Cladogram of species
Note: Cladogram consists of circles radiating from the inside out, each corresponding to different levels of classification such as kingdom, phylum, class, order, family, and genus. The lines connecting the circles represent their relationships. Each circle node represents a species, with yellow nodes indicating insignificant differences between groups, while significantly different species are colored according to their group. Colored sector regions represent the subclassification intervals of characteristic microorganisms.A, B, C, D, E and F represent the E. tetradactylum juveniles of 37, 44, 58, 72, 86 and 114 days of age, respectively.
图 7 相关性分析热图
注:色块的颜色表示相关系数,其中红色代表正相关,蓝色代表负相关;颜色深浅代表相关系数大小 (*. p<0.05)。
Figure 7. Heat map of correlation analysis
Note: The colors of the color blocks represent the correlation coefficients. Red represents positive correlation and blue represents negative correlation; the color shades represent the size of the correlation coefficients (*. p<0.05).
表 1 不同日龄四指马鲅生长性能
Table 1 Growth performance of E. tetradactylum at different days of age
指标
Indicator组别 Group A B C D E F 体长 Body length/cm 2.37±0.25 3.47±0.21 6.20±.032 7.00±0.30 10.08±0.36 13.20±1.19 体质量 Body mass/g 0.24±0.02 0.71±0.07 1.60±0.34 3.46±0.65 8.62±1.21 44.72±11.37 特定生长率SGR/(%·d−1) 8.58±0.23b 15.36±1.73a 5.72±1.45c 5.38±1.62c 6.43±1.16bc 6.52±2.02bc 平均日增质量ADG/(g·d−1) 0.01±0.00a 0.07±0.01a 0.06±0.02a 0.13±0.05a 0.52±0.29b 1.15±0.22a 肥满度CF/(g·cm−3) 2.88±0.95a 2.47±0.32a 0.85±0.03c 1.44±0.17b 1.67±0.09b 1.77±0.42b 注:同行中不同小写字母间存在显著性差异 (p<0.05)。A、B、C、D、E、F分别表示37、44、58、72、86、114日龄的幼鱼。 Note: Values with different lowercase letters within the same row are significantly different (p<0.05). A, B, C, D, E and F represent the E. tetradactylum juveniles of 37, 44, 58, 72, 86 and 114 days of age, respectively. 表 2 不同日龄四指马鲅肠道组织结构变化
Table 2 Intestinal tissue structure of E. tetradactylum at different days of age
指标
Indicator组别 Group B C D E F 绒毛长度 Villus length/μm 91.69±8.42d 124.33±22.00c 168.11±42.35b 189.60±25.91b 312.78±45.92a 绒毛宽度 Villus width/μm 19.39±4.39c 20.14±3.29c 28.40±4.33b 38.28±8.84a 42.76±10.51a 肌层厚度 Muscular thickness/μm 12.31±1.57d 16.57±6.22bc 20.03±5.86c 31.58±5.07b 56.72±15.13a 注:同行中不同小写字母间存在显著性差异 (p<0.05)。B、C、D、E、F分别表示44、58、72、86、114日龄的幼鱼。 Note: Values with different lowercase letters within the same row are significantly different (p<0.05). B, C, D, E and F represent the E. tetradactylum juveniles of 44, 58, 72, 86 and 114 days of age, respectively. 表 3 不同日龄四指马鲅鱼肠道菌群Alpha多样性统计
Table 3 Alpha diversity of intestinal microbiota in E. tetradactylum at different days of age
指标
Indicator组别 Group A B C D E F Chao 1 指数 Chao 1 index 1 019.70±247.37 919.92±24.98 795.38±167.93 818.00±32.87 739.96±231.98 805.94±80.17 ACE 指数 ACE index 986.88±211.51 938.39±29.54 810.60±163.07 817.51±20.64 774.20±258.03 866.18±93.10 香农指数Shannon index 4.65±0.01a 4.21±0.08ab 4.72±0.61a 3.86±0.47ab 3.16±1.58b 4.52±0.37ab 辛普森指数Simpson index 0.97±0.01 0.95±0.01 0.96±0.02 0.93±0.04 0.94±0.02 0.97±0.02 注:同行中不同小写字母间存在显著性差异 (p<0.05)。A、B、C、D、E、F分别表示37、44、58、72、86、114日龄的幼鱼。 Note: Values with different lowercase letters within the same row are significantly different (p<0.05). A, B, C, D, E and F represent the E. tetradactylum juveniles of 37, 44, 58, 72, 86 and 114 days of age, respectively. -
[1] WANG J J, PENG S, FEI Y, et al. Low mtDNA Cytb diversity and shallow population structure of Eleutheronema tetradactylum in the East China Sea and the South China Sea[J]. Biochem Syst Ecol, 2014, 55: 268-274. doi: 10.1016/j.bse.2014.03.026
[2] MOORE B R, STAPLEY J, ALLLSOP Q, et al. Stock structure of blue threadfin Eleutheronema tetradactylum across northern Australia, as indicated by parasites[J]. J Fish Biol, 2011, 78(3): 923-936. doi: 10.1111/j.1095-8649.2011.02917.x
[3] 郭海波, 吴益春, 罗海军, 等. 四指马鲅鱼的营养成分分析[J]. 食品安全质量检测学报, 2017, 8(1): 88-92. [4] BRUIN G D, RUSSELL B, BOGUSCH A. FAO species identification field guide for fishery purposes. Fishes of the Cambodian Mekong[J]. 1996, 2: 156-158.
[5] 牛莹月, 罗志平, 区又君, 等. 不同盐度对四指马鲅幼鱼存活及鳃组织形态变化的影响[J]. 南方农业学报, 2021, 52(6): 1719-1726. [6] 蓝军南, 区又君, 温久福, 等. 四指马鲅卵巢发育和卵子发生的组织学观察[J]. 中国水产科学, 2020, 27(12): 1415-1426. [7] 谢木娇, 区又君, 李加儿, 等. 四指马鲅 (Eleutheronema tetradactylum) 消化系统胚后发育组织学观察[J]. 渔业科学进展, 2017, 38(2): 50-58. [8] 区又君, 刘奇奇, 温久福, 等. 急性低温胁迫对四指马鲅幼鱼肝脏、肌肉以及鳃组织结构的影响[J]. 生态科学, 2018, 37(5): 53-59. [9] 江梁正, 王珺, 区又君, 等. 2种盐度池塘养殖四指马鲅的肌肉营养组成分析[J]. 中国渔业质量与标准, 2023, 13(2): 11-17. [10] 李俊伟, 区又君, 温久福, 等. 室内循环水和池塘养殖四指马鲅的生长性能及形态性状与体质量的相关性研究[J]. 南方水产科学, 2020, 16(1): 27-35. [11] 高淼, 周秋白, 王自蕊, 等. 饲料原料粉碎粒度对大鳞副泥鳅生长性能的影响[J]. 饲料工业, 2018, 39(22): 21-25. [12] 李芳源, 王常安, 刘红柏. 丁酸梭菌对水产动物营养及免疫研究进展[J]. 水产学杂志, 2021, 34(5): 104-110. [13] 梁祖銮, 赵吉臣, 廖敏泽, 等. 不同生长阶段中国花鲈肠道和环境微生物群落分析[J]. 大连海洋大学学报, 2024, 39(2): 215-224. [14] GAO Y J, YANG H J, LIU Y J, et al. Effects of graded levels of threonine on growth performance, biochemical parameters and intestine morphology of juvenile grass carp Ctenopharyngodon idella[J]. Aquaculture, 2014, 424: 113-119.
[15] 郭鑫伟, 张洋, 迟淑艳, 等. 三种铁源对珍珠龙胆石斑鱼幼鱼生长性能、肝脏抗氧化酶活性及肠道发育形态的影响[J]. 中国海洋大学学报, 2020, 50(11): 53-61. [16] 马兴宇, 唐忠林, 陈树桥, 等. 转食饲料对大口黑鲈幼鱼的存活率、抗氧化酶和消化酶活性及肠道菌群的影响[J]. 中国水产科学, 2024, 31(4): 403-415. [17] 范梓健, 张紫玥, 曹建萌, 等. 不同开口饵料对罗非鱼发育早期肠道菌群结构的影响[J]. 农业生物技术学报, 2023, 31(5): 1032-1042. [18] 孟晓林, 李文均, 聂国兴. 鱼类肠道菌群影响因子研究进展[J]. 水产学报, 2019, 43(1): 143-155. [19] 张碧云, 杨红玲, 汪攀, 等. 鱼类肠道微生物与宿主免疫系统相互作用研究进展[J]. 微生物学报, 2021, 61(10): 3046-3058. [20] 郁维娜, 戴文芳, 陶震, 等. 健康与患病凡纳滨对虾肠道菌群结构及功能差异研究[J]. 水产学报, 2018, 42(3): 399-409. [21] 李鸣霄, 强俊, 徐钢春, 等. 不同养殖阶段的大口黑鲈肠道结构和肠道微生物组成变化的比较[J]. 动物营养学报, 2023, 35(9): 5886-5903. [22] BOILEN E, RIDEOUT J R, DILLON M R, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2[J]. Nat Biotechnol, 2019, 37(8): 852-857. doi: 10.1038/s41587-019-0209-9
[23] SEGATA N, LZARD J, WALDRON L, et al. Metagenomic biomarker discovery and explanation[J]. Genome Biol, 2011, 12(6): 60. doi: 10.1186/gb-2011-12-6-r60
[24] LI X F, TIAN H Y, ZHANG D Z, et al. Feeding frequency affects stress, innate immunity and disease resistance of juvenile blunt snout bream Megalobrama amblycephala[J]. Fish Shellfish Immunol, 2014, 38(1): 80-87. doi: 10.1016/j.fsi.2014.03.005
[25] BUDDINGTON R K. Digestive secretions of lake sturgeon, Acipenser fulvescens, during early development[J]. Fish Shellfish Immunol, 1985, 26(6): 715-723.
[26] 谭泽宇, 李涛, 姜敬哲, 等. 喷点雪印小丑鱼的胚胎及仔、稚、幼鱼形态发育观察[J]. 南方水产科学, 2024, 20(2): 73-82. doi: 10.12131/20230119 [27] 黄伟卿, 宋炜, 刘铮, 等. 大黄鱼仔、稚鱼期的脊椎骨及附肢骨生长发育研究[J]. 海洋渔业, 2024, 3: 1-10. [28] 兰真强, 郑纪涛, 陈芸, 等. 丝鳍圆天竺鲷的繁殖习性、胚胎发育和胚后发育观察[J]. 热带海洋学报, 2024, 43(1): 116-125. [29] HUANG B C, ZHANG S, DONG X H, et al. Effects of fishmeal replacement by black soldier fly on growth performance, digestive enzyme activity, intestine morphology, intestinal flora and immune response of pearl gentian grouper (Epinephelus fuscoguttatus♀ × Epinephelus lanceolatus♂)[J]. Fish Shellfish Immunol, 2022, 120: 497-506. doi: 10.1016/j.fsi.2021.12.027
[30] 赵彦花, 区又君, 李加儿, 等. 黄唇鱼消化系统组织结构及黏液细胞分布特征[J]. 渔业科学进展, 2019, 40(3): 80-86. [31] 曹剑香, 简纪常, 黄洋, 等. 多鳞鱚消化系统的形态学和组织学研究[J]. 水产科学, 2010, 29(6): 339-343. [32] 区又君, 何永亮, 李加儿. 卵形鲳鲹消化系统的胚后发育[J]. 台湾海峡, 2011, 30(4): 533-539. [33] 向佳丽.南极鱼类肠道微生物群落结构特征及其与环境因子关系的研究[D]. 上海: 上海海洋大学, 2023: 44-46. [34] SYLVAIN F E, CHEAIB B, LLEWELLYN M, et al. PH drop impacts differentially skin and gut microbiota of the Amazonian fish tambaqui (Colossoma macropomum)[J]. Sci Rep, 2016, 6(1): 32032. doi: 10.1038/srep32032
[35] STEPHENS W Z, BURNS A R, STAGAMAN K, et al. The composition of the zebrafish intestinal microbial community varies across development[J]. ISME J, 2016, 10(3): 644-654. doi: 10.1038/ismej.2015.140
[36] NIELSEN S, WILKES W J, VERGES A, et al. Microbiome patterns across the gastrointestinal tract of the rabbitfish Siganus fuscescens[J]. PeerJ, 2017, 5: e3317. doi: 10.7717/peerj.3317
[37] HUANG Q, SHAM R C, DENG Y, et al. Diversity of gut microbiomes in marine fishes is shaped by host–related factors[J]. Mol Ecol, 2020, 29(24): 5019-5034. doi: 10.1111/mec.15699
[38] 高权新, 吴天星, 王进波. 肠道微生物与寄主的共生关系研究进展[J]. 动物营养学报, 2010, 22(3): 519-526. [39] LI W H, ZHOU Z X, LI Q H, et al. Successional changes of microbial communities and host-microbiota interactions contribute to dietary adaptation in allodiploid hybrid fish[J]. Microb Ecol, 2023, 85(4): 1190-1201. doi: 10.1007/s00248-022-01993-y
[40] ZHAO C Y, WANG J, REN W H, et al. Histological, immune, and intestine microbiota responses of the intestine of rainbow trout (Oncorhynchus mykiss) to high temperature stress[J]. Aquaculture, 2024, 582: 740465. doi: 10.1016/j.aquaculture.2023.740465
[41] SHANG Y Q, ZHONG H M, LIU G, et al. Characteristics of microbiota in different segments of the digestive tract of Lycodon rufozonatus[J]. Animals, 2023, 13(4): 731. doi: 10.3390/ani13040731
[42] WANG W Z, HUANG J S, ZHANG J D, et al. Effects of hypoxia stress on the intestinal microflora of juvenile of cobia (Rachycentron canadum)[J]. Aquaculture, 2021, 536: 736419. doi: 10.1016/j.aquaculture.2021.736419
[43] ZHAO Y, LI S P, LESSING D J, et al. The attenuating effects of synbiotic containing Cetobacterium somerae and Astragalus polysaccharide against trichlorfon-induced hepatotoxicity in crucian carp (Carassius carassius)[J]. J Hazard Mater, 2024, 461: 132621. doi: 10.1016/j.jhazmat.2023.132621
[44] OFEK T, LALZAR M, LAVIAD S, et al. Comparative study of intestinal microbiota composition of six edible fish species[J]. Front Microbiol, 2021, 12: 760266. doi: 10.3389/fmicb.2021.760266
[45] 谷雪勤, 王庆奎, 王洋. 海水鱼循环水养殖系统中细菌群落结构研究进展[J]. 天津农学院学报, 2022, 29(2): 76-84. [46] WEI N, WANG C C, XIAO S J, et al. Intestinal microbiota in large yellow croaker, Larimichthys crocea, at different ages[J]. J World Aquacult Soc, 2018, 49(1): 256-267. doi: 10.1111/jwas.12463
[47] ZHANG Y, WEN B, MENG L J, et al. Dynamic changes of gut microbiota of discus fish (Symphysodon haraldi) at different feeding stages[J]. Aquaculture, 2021, 531: 735912. doi: 10.1016/j.aquaculture.2020.735912
[48] CICALA F, LAGO L A, GOMEZ G B, et al. Gut microbiota shifts in the giant tiger shrimp, Penaeus monodon, during the postlarvae, juvenile, and adult stages[J]. Aquac Int, 2020, 28(4): 1421-1433. doi: 10.1007/s10499-020-00532-1
[49] 李路宽. 四种新型蛋白源替代鱼粉对大口黑鲈肠道健康的影响[D]. 武汉: 华中农业大学, 2023: 39-40. [50] 汤菊芬, 黄瑜, 蔡佳, 等. 中草药复合微生态制剂对吉富罗非鱼 (Oreochromis niloticus) 生长、肠道菌群及抗病力的影响[J]. 渔业科学进展, 2016, 37(4): 104-109. [51] 李艳丽, 杨垒, 张志昊, 等. 好氧反硝化细菌Burkholderia sp. ZH8的脱氮特性与生物强化作用研究[J]. 中国环境科学, 2024, 44(8): 4282-4291. [52] 潘丹, 黄巧云, 陈雯莉. 两株异养硝化细菌的分离鉴定及其脱氮特性[J]. 微生物学报, 2011, 51(10): 1382-1389. [53] 杨浩, 张国珍, 杨晓妮, 等. 16S rRNA高通量测序研究集雨窖水中微生物群落结构及多样性[J]. 环境科学, 2017, 38(4): 1704-1716. [54] DACEWICA E, LENART B A. Waste polyurethane foams as biomass carriers in the treatment process of domestic sewage with increased ammonium nitrogen content[J]. Materials, 2023, 16(2): 619. doi: 10.3390/ma16020619
[55] 刘增新, 柳学周, 史宝, 等. 牙鲆 (Paralichthys olivaceus) 仔稚幼鱼肠道菌群结构比较分析[J]. 渔业科学进展, 2017, 38(1): 111-119. [56] 王雨婷, 周荣翔, 李霁虹, 等. 抗弧菌光合细菌的分离鉴定及对氨氮和亚硝态氮的降解特性[J]. 南方水产科学, 2021, 17(5): 26-33. doi: 10.12131/20210016 [57] 付保荣, 曹向宇, 冷阳, 等. 光合细菌对水产养殖水质和水生生物的影响[J]. 生态科学, 2008(2): 102-106. doi: 10.3969/j.issn.1008-8873.2008.02.007 [58] 陈明霞, 李和阳, 李刚, 等. 深圳海域弧菌种类组成、数量分布及其与环境因子的关系研究[J]. 海洋学报, 2010, 32(5): 117-126.