Analysis of bacterial community and diversity in gill tissues of bony fishes in adjacent South China Sea
-
摘要: 围绕中国南海丰富的渔业资源,采用16S rDNA扩增子测序技术,分析了南海近岸硬骨鱼类鳃组织微生物群落分布特征,并探讨了7个不同站点细菌群落结构的差异。结果显示,测序共获得有效拼接片段 (Clean tags) 2 952 366条,平均每个文库56 776条。分别在门、属水平对其优势类群进行了分析,其中门水平上变形菌门 (Proteobacteria) 最高 (71.3%),属水平上变形菌门的不动杆菌属 (Acinetobacter) 最高 (17.2%)。不同站点的α多样性具有显著差异,其中I9和H8站点的物种丰富度 (Chao1指数) 最高,D3站点的多样性 (Shannon指数) 最高。不同站点来源样品之间的β多样性具有显著差异 (P<0.05),但在宿主分类的目水平上无显著差异 (P>0.05)。中国南海近岸硬骨鱼鳃组织中的细菌群落组成丰富,采样站位相比宿主分类对鳃组织上细菌群落具有更重要的影响,它们可能在辅助宿主营养物质转运及代谢方面发挥积极作用。
-
关键词:
- 硬骨鱼鳃 /
- 16S rDNA扩增子 /
- 细菌多样性 /
- 南海
Abstract: Focusing on the abundant fishery resources in the South China Sea, we analyzed the distribution characteristics of the microbial community in the gill tissue of bony fishes in the adjacent South China Sea by using 16S amplicon sequencing technology. In addition, we discussed the differences in the bacterial community structure at different sites. The results show that the sequencing obtained a total of 2 952 366 effective spliced fragments (Clean tags) with an average of 56 776 in each library. The dominant taxa were analyzed at phylum and genus levels. Proteobacteria was the highest (71.3%) at phylum level, and Acinetobacter was the highest at genus level (17.2%). The alpha diversity of different sites was significantly different, and the I9 and H8 sites had the highest species richness (Chao1 index), and the D3 site had the highest diversity (Shannon index). The beta diversity of samples from different sites was significantly different (P<0.05), but there was no significant difference between the host classification (Order level) groups (P>0.05). The composition of bacterial communities in the gill tissues of bony fishes in the South China Sea was abundant. Compared with the classification of the host, the sampling station has more important influence on the community of the bacterial in the gill tissues, and may play an active role in assisting the host's nutrient transport and metabolism.-
Key words:
- Bony fish gill /
- 16S rDNA amplicon /
- Bacterial diversity /
- South China Sea
-
图 5 Alpha多样性指数箱线图
注:a—e. 基于站点分组的Chao1、ACE、Observed、Simpson、Shannon 多样性指数箱线图,不同颜色代表不同站点的样本;f. 基于样本宿主(目水平) 分类分组的Chao1 指数箱线图。
Figure 5. Boxplot of Alpha diversity index
Note: a–e. Boxplot of diversity indexes of Chao1, ACE, Observed, Simpson and Shannon based on site grouping, and different colors represent samples of different sites; f. Boxplot of Chao1 index based on sample host (Order level) classification.
图 6 基于不同分组方式的NMDS分析 [非参数多元方差分析 (PerMANOVA)]
注:a. 基于站点分组的 NMDS 分析; b. 基于样本的目水平分组。图中两个样本点越接近,表示两个样本物种组成越相似。横纵坐标表示样本间的相对距离,无实际意义。图形中的每一个点代表一个样本,不同颜色代表不同的样本分组信息。NMDS 结果的优劣用胁迫系数 (stress) 来衡量,此值越小越好,当小于 0.2 时表示可以用 NMDS 的二维点图表示组间或组内差异。
Figure 6. NMDS analysis based on different grouping modes [Nonparametric multivariate analysis of variance (PerMANOVA)]
Note: a. NMDS analysis based on site grouping; b. Grouping based on the order level of the sample. The closer the two sample points are,the more similar the composition of the two sample species are. The abscissa and ordinate indicate the relative distance between samples, which is of no practical significance. Each point in the graph represents a sample, and different colors represent different sample grouping information. The advantages and disadvantages of testing NMDS results are measured by stress coefficients, the smaller, the better, and the value less than 0.2 indicates that the two-dimensional dot plot of NMDS can be used to represent differences between or within groups.
图 7 22个KEGG-L2功能层级相对丰度热图
注:横坐标代表不同的样本,纵坐标代表不同的功能层级,图片上方色条表示不同分组,并按照样本聚类。
Figure 7. Heatmap of relative abundance of 22 KEGG-L2 functional levels
Note: The abscissa represents different samples; the ordinate represents different functional levels; the color bars indicate different groups, which are clustered according to samples.
表 1 样品分类信息
Table 1. Taxonomy information of samples
站点
Sampling station样品编号
Sample ID目
Order物种
SpeciesD3 D3.1S 鲈形目 条纹䱨 Terapon theraps D3 D3.2S 鲈形目 日本金线鱼 Nemipterus japonicus D3 D3C.2S 鲈形目 列牙䱨 Pelates quadrilineatus F3 F3.2S 鲈形目 鹿斑鲾 Leiognathus ruconius F3 F3.4S 鲽形目 少牙斑鲆 Pseudorhombus oligodon F3 F3.5S 鲈形目 带鱼 Trichiurus lepturus D8 D8.1S 鲈形目 蓝圆鲹 Decapterus maruadsi D8 D8.3S 鳗鲡目 网纹裸胸鳝 Gymnothorax reticularis D8 D8.5S 鮟目 黑鮟 Lophiomus setigerus D8 D8.6S 鲉形目 环纹蓑鲥 Pterois lunulata D8 D8.8S 鲈形目 深水金线鱼 Nemipterus bathybius D8 D8.12S 鲉形目 深海红娘鱼 Lepidotrigla abyssalis F8 F8.1S 鲉形目 环纹蓑鲥 Pterois lunulata F8 F8.4S 鲉形目 单棘豹鲂鮄 Daicocus peterseni F8 F8.8S 鲈形目 二长棘鲷 Parargyrops edita F8 F8.11S 鮟目 棘茄鱼 Halieutaea stellata F8 F8.16S 鲉形目 日本红娘鱼 Lepidotrigla japonice F8 F8.17S 鲉形目 拟蓑鲉 Parapterois heterurus F8 F8.18S 海龙鱼目 鳞烟管鱼 Fistularia petimba H3 H3.3S 鲈形目 翼红娘鱼 Lepidotrigla alata H3 H3.4S 鲈形目 中华䲢 Uranoscopus chinensis H3 H3.5S 鳗鲡目 大鳞鳞头鲉 Sebastapistes megalepis H3 H3.7S 鮟目 无备虎鲉 Minous inermis H8 H8.1S 鲉形目 大鳞短额鲆 Engyprosopon grandisquama H8 H8.2S 鲈形目 鳞烟管鱼 Fistularia petimba H8 H8.3S 鲉形目 繁星鲆 Bothus myriaster H8 H8.7S 鲉形目 孔鰕虎鱼 Trypauchen vagina H8 H8.10S 鲉形目 斑鰶 Konosirus punctatus H8 H8.12S 鲈形目 皮氏叫姑鱼 Johnius belangerii H8 H8.15S 鮟目 短吻鲾 Leiognathus brevirostris H8 H8.24S 鲉形目 大头狗母鱼 Trachinocephalus myops H8 H8.25S 鲉形目 大头狗母鱼 Trachinocephalus myops I9 I9.1S 海龙鱼目 大鳞鳞头鲉 Sebastapistes megalepis I9 I9.5S 鲉形目 黄纹拟鲈 Parapercis xanthozona I9 I9.7S 鲈形目 寿鱼 Banjos banjos I9 I9.8S 鲉形目 棕斑宽吻鲀 Amblyrhynchotes rufopunctatus I9 I9.9S 鲉形目 单棘豹鲂鮄 Daicocus peterseni I9 I9.10S 鲽形目 松球鱼 Monocentris japonica I9 I9.11S 海龙鱼目 新平鲉 Neosebastes entaxis I9 I9.12S 鲽形目 二长棘鲷 Parargyrops edita I9 I9.13S 鲈形目 横带眶棘鲈 Scolopsis inermis I9 I9.14S 鲱形目 盔蓑鲉 Ebosia bleekeri I9 I9.15S 鲈形目 冠鲽 Samaris cristatus I9 I9.16S 鲈形目 叉斑狗母鱼 Synodus macrops I9 I9.17S 海龙鱼目 美尾䲗 Calliurichthys japonicus I9 I9.18S 鲉形目 瑞氏红鲂鮄 Satyrichthys rieffeli I9 I9.20S 鲉形目 虻鲉 Erisphex pottii I9 I9.21S 鲉形目 红鲬 Bembras japonicus I9 I9.22S 鲈形目 六带拟鲈 Parapercis sexfasciata I9 I9.23S 鲈形目 黄纹拟鲈 Parapercis xanthozona I9 I9.24S 鲈形目 尖牙鲈 Synagrops japonicus I9 I9.25S 鲈形目 水珍鱼 Argentina kagoshimae -
[1] 王正, 毛志华, 李晓娟. 气候变化对南海浮游植物藻华形成的影响研究进展[J]. 环境污染与防治, 2017, 39(12): 1384-1390. [2] 苏莹佳, 陈国宝, 周艳波, 等. 2015—2017年南海海域伏季休渔制度实施效果评价[J]. 南方水产科学, 2019, 15(2): 20-28. doi: 10.12131/20180149 [3] 史登福, 张魁, 蔡研聪, 等. 南海北部带鱼群体结构及生长、死亡和性成熟参数估计[J]. 南方水产科学, 2020, 16(5): 51-59. doi: 10.12131/20200055 [4] 陈大刚, 张美昭. 中国海洋鱼类[M]. 青岛: 中国海洋大学出版社, 2016: 1-2. [5] 张俊, 江艳娥, 陈作志, 等. 南海中南部中层鱼资源声学积分值及时空分布初探[J]. 中国水产科学, 2017(1): 120-135. [6] 肖军. 常见养殖鱼类鳃病的病理症状及防治方法[J]. 当代水产, 2019, 44(6): 97-101. doi: 10.3969/j.issn.1674-9049.2019.06.027 [7] TOENSHOFF E R, KVELLESTAD A, MITCHELL S O, et al. A novel betaproteobacterial agent of gill epitheliocystis in seawater farmed Atlantic salmon (Salmo salar)[J]. PLOS ONE, 2012, 7(3): e32696. doi: 10.1371/journal.pone.0032696 [8] WANG C Y, LI Z B. Growth performance, digestive enzyme activity and immune response of Japanese sea bass, Lateolabrax japonicus fed with fructooligosaccharide[J]. Aquac Nutr, 2020, 26(2): 296-305. doi: 10.1111/anu.12990 [9] DOUJET T L, SANTI C D, KLEMETSEN T, et al. Closely-related Photobacterium strains comprise the majority of bacteria in the gut of migrating Atlantic cod (Gadus morhua)[J]. Microbiome, 2019, 7(1): 64. doi: 10.1186/s40168-019-0681-y [10] FOGARTY C, BURGESS C M, COTTER P D, et al. Diversity and composition of the gut microbiota of Atlantic salmon (Salmo salar) farmed in Irish waters[J]. J Appl Microbiol, 2019, 127(3): 648-657. doi: 10.1111/jam.14291 [11] 樊瑞峰, 王印庚, 梁友, 等. 一株广温性大菱鲆肠道益生菌的筛选与鉴定[J]. 渔业科学进展, 2011, 32(1): 40-46. doi: 10.3969/j.issn.1000-7075.2011.01.006 [12] 郭国军, 覃映雪, 陈强, 等. 大黄鱼病原副溶血弧菌拮抗菌的筛选[J]. 海洋学报, 2008, 30(1): 127-134. [13] 王小彤, 许强华. 应用高通量测序法分析四种南极鱼鳃组织的微生物群落多样性[C]//2018年中国水产学会学术年会论文摘要集. 西安: 中国水产学会, 2018: 132. [14] PRATTE Z A, BESSON M, HOLLMAN R D, et al. The gills of reef fish support a distinct microbiome influenced by host-specific factors[J]. Appl Environ Microbiol, 2018, 84(9): e00063-18. [15] CAPORASO J G, LAUBER C L, WALTERS W A, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample[J]. Proc Nat Acad Sci, 2011, 108(Suppl 1): 4516-4522. [16] 杨翠, 李祖明. 分子生物学技术在肠道微生物研究中的应用进展[J]. 中国微生态学杂志, 2017, 29(2): 229-233,238. [17] CHEN S, ZHOU Y, CHEN Y, et al. fastp: an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics, 2018, 34(17): i884-i890. doi: 10.1093/bioinformatics/bty560 [18] EDGAR R C. UPARSE: highly accurate OTU sequences from microbial amplicon reads[J]. Nat Methods, 2013, 10(10): 996-998. doi: 10.1038/nmeth.2604 [19] CHAO A. Nonparametric estimation of the number of classes in a population[J]. Scand J Stat, 1984, 11(4): 265-270. [20] CHAO A, YANG M C K. Stopping rules and estimation for recapture debugging with unequal failure rates[J]. Biometrika, 1993, 80(1): 193-201. doi: 10.1093/biomet/80.1.193 [21] SHANNON C E. A mathematical theory of communication[J]. Bell Syst Tech J, 1948, 27(3): 379-423. doi: 10.1002/j.1538-7305.1948.tb01338.x [22] DOUGLAS G M, BEIKO R G, LANGILLE M G I. Predicting the functional potential of the microbiome from marker genes using PICRUSt[J]. Methods Mol Biol, 2018, 1849: 169-177. [23] CHONG J, LIU P, ZHOU G, et al. Using Microbiome analyst for comprehensive statistical, functional, and meta-analysis of microbiome data[J]. Nat Protoc, 2020, 15(3): 799-821. doi: 10.1038/s41596-019-0264-1 [24] 穆大帅, 卢德臣, 郑维爽, 等. 我国海洋细菌新物种鉴定与资源研发进展[J]. 生物资源, 2017, 39(6): 391-397. [25] 白洁, 刘小沙, 侯瑞, 等. 南海南部海域浮游细菌群落特征及影响因素研究[J]. 中国环境科学, 2014(11): 2950-2957. [26] 徐新亚, 杨宏, 宁小清, 等. 北部湾海洋微生物物种多样性与化学多样性研究进展[J]. 广西科学, 2020, 27(5): 433-450,461. [27] 王丹. 海洋细菌SAR86全球分布及其多样性研究[D]. 舟山: 浙江海洋大学, 2020: 11-16. [28] DELMONT T O, QUINCE C, SHAIBER A, et al. Nitrogen-fixing populations of Planctomycetes and Proteobacteria are abundant in surface ocean metagenomes[J]. Nat Microbiol, 2018, 3(7): 804-813. doi: 10.1038/s41564-018-0176-9 [29] CHEN Y, PENNER G B, LI M, et al. Changes in bacterial diversity associated with epithelial tissue in the beef cow rumen during the transition to a high-grain diet[J]. Appl Environ Microbiol, 2011, 77(16): 5770-5781. doi: 10.1128/AEM.00375-11 [30] EVANS N J, BROWN J M, MURRAY R D, et al. Characterization of novel bovine gastrointestinal tract treponema isolates and comparison with bovine digital dermatitis treponemes[J]. Appl Environ Microbiol, 2011, 77(1): 138-147. doi: 10.1128/AEM.00993-10 [31] LEGRAND T P R A, CATALANO S R, WOS-OXLEY M L, et al. The inner workings of the outer surface: skin and gill microbiota as indicators of changing gut health in yellowtail kingfish[J]. Front Microbiol, 2018, 8: 2664. doi: 10.3389/fmicb.2017.02664 [32] KESSEL M A H J V, MESMAN R J, ARSHAD A, et al. Branchial nitrogen cycle symbionts can remove ammonia in fish gills[J]. Environ Microbiol Rep, 2016, 8(5): 590-594. doi: 10.1111/1758-2229.12407 [33] 刘玉华, 王慧, 胡晓珂. 不动杆菌属 (Acinetobacter) 细菌降解石油烃的研究进展[J]. 微生物学通报, 2016, 43(7): 1579-1589. [34] 欧阳晴晴, 姚伯卿, 焦亚彬, 等. 南极中山站长城站土壤及海洋沉积物可培养细菌多样性[J]. 生物资源, 2021, 43(2): 133-142. [35] 翟万营, 汪倩, 王佳华, 等. 南极中山站附近海水微生物多样性和宏基因组分析[J]. 上海海洋大学学报, 2022, 31(1): 288-297. [36] ZHANG E, THIBAUT L M, TERAUDS A, et al. Lifting the veil on arid-to-hyperarid Antarctic soil microbiomes: a tale of two oases[J]. Microbiome, 2020, 8(1): 37. doi: 10.1186/s40168-020-00809-w -