海洋保护区是全球使用最广泛的海洋管理工具之一,它对于保护生物多样性至关重要,它可以提供保护,使其免受人为活动的威胁。微生物群落在有机污染物生物地球化学循环中起着十分关键的作用。微生物群落结构和多样性及其变化在一定程度上反映了土壤的质量。为准确掌握我国黄海沿海潮间带沉积物微生物群落,了解海洋保护区与相邻区域微生物群落结构差异,于2018年9月对辽宁省大连长海珍稀海洋生物保护区及附近区域潮间带沉积物微生物开展了详细调查。应用Illumina MiSeq宏基因组高通量测序技术,在调查区域潮间带共鉴定沉积物微生物43门类,变形菌门为最主要的优势门类,其序列数占总数的45.7%。其它较为丰富的门类包括酸杆菌门、拟杆菌门、放线菌门、浮霉菌门、蓝细菌门和疣微菌门。采用Mothur软件计算Chao1、Shannon和Simpson多样性指数,结果发现:菌群丰度指数Chao1的均值为5493,Shannon菌群多样性的指数为9.71,Simpson指数0.994,显示出海域沉积物中微生物多样性较高,通过主成份分析发现,保护区和大连入海河流潮间带微生物群落结构差异较大。 Marine protected areas (MPAs) are one of the most widely utilized marine management tools globally. Marine protected areas are essential for conservation of biodiversity by offering protection from major anthropogenic threats. Microbial community plays a key role in biogeochemical cycles of organic pollutants. Their community structure and diversity, and the changes reflect certain aspects of sediment quality. In order to study the microbial community structure in intertidal sediments of coastal area in Yellow Sea and to understand the differences in the structure of microbial communities between marine protected areas and adjacent areas, metagenomic analysis based on the Illumina Miseq high throughput sequencing technology was applied in Liaoning Dalian Changhai Marine Precious Biological Provincial Nature Reserve intertidal zone in September 2018. A total of 43 phyla of microoganisms were detected from the intertidal sediment samples. Proteobacteria was the dominant group, comprising 45.7% of the total sequences. Other abundant groups include phylum Acidobacteria, Bacteroidetes, Actinobacteria, Planctomycetes, Cyanobacteria and Verrucomicrobia. The Chao1, Shannon and Simpson diversity indexes were calculated using Mothur software. The average value of Chao1 index was 5493, Shannon diversity index was 9.71 for all the samples, and the Simpson index was 0.994. These indices indicated high microbial community diversity in the area. Through principal component analysis, it was found that the structure of microbial communities in the intertidal zone of MPAs and Dalian's estuarine rivers was quite different.
——以大连长海珍稀海洋生物保护区为例
唐嘉威1,黄备2*,孙沛雯3,孙明霞4,刘子成3,李翔宇3,俞存根1
1浙江海洋大学水产学院,浙江 舟山
2浙江省舟山海洋生态环境监测站,浙江 舟山
3辽宁省大连生态环境监测中心,辽宁 大连
4辽宁省大连市长海县环境监测站,辽宁 大连
收稿日期:2020年2月19日;录用日期:2020年3月6日;发布日期:2020年3月13日
海洋保护区是全球使用最广泛的海洋管理工具之一,它对于保护生物多样性至关重要,它可以提供保护,使其免受人为活动的威胁。微生物群落在有机污染物生物地球化学循环中起着十分关键的作用。微生物群落结构和多样性及其变化在一定程度上反映了土壤的质量。为准确掌握我国黄海沿海潮间带沉积物微生物群落,了解海洋保护区与相邻区域微生物群落结构差异,于2018年9月对辽宁省大连长海珍稀海洋生物保护区及附近区域潮间带沉积物微生物开展了详细调查。应用Illumina MiSeq宏基因组高通量测序技术,在调查区域潮间带共鉴定沉积物微生物43门类,变形菌门为最主要的优势门类,其序列数占总数的45.7%。其它较为丰富的门类包括酸杆菌门、拟杆菌门、放线菌门、浮霉菌门、蓝细菌门和疣微菌门。采用Mothur软件计算Chao1、Shannon和Simpson多样性指数,结果发现:菌群丰度指数Chao1的均值为5493,Shannon菌群多样性的指数为9.71,Simpson指数0.994,显示出海域沉积物中微生物多样性较高,通过主成份分析发现,保护区和大连入海河流潮间带微生物群落结构差异较大。
关键词 :海洋保护区,微生物群落,潮间带,高通量测序
Copyright © 2020 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
由于人类活动导致的生境丧失等因素的影响,全球范围内的生物多样性持续减少。由于在生物多样性保护上可投入的资源有限,加上生物资源的分布本身存在地区差异,集中力量优先保护一些重要的地区,是目前较为现实和高效的途径 [
尽管MPA已在世界各地建立,但它们与其他生态系统之间的联系非常复杂,需要进行大量研究。许多海洋保护区靠近污染源,例如工业园区、港口码头、水产养殖场、石油码头和平台,市区和排污口等。由于海洋和河口的流动性使污染物可以长距离传输,因此陆上和海洋活动会持续对海洋生态系统产生负面影响,距污染源不同距离的MPA也可能受到威胁 [
微生物群落是地球生命的基础,并通过催化驱动全球营养循环的生物地球化学反应以及与生物群落直接相互作用,提供基本的生态系统服务 [
本文采用基于Illumina MiSeq高通量测序平台,对辽宁省大连市长海珍稀海洋生物保护区及相邻区域潮间带沉积物微生物多样性进行了深入的调查,通过分析沉积物微生物群落16S rRNA基因测序数据,探索并区分保护区与邻近区域之间微生物组成和多样性特征,讨论微生物群落与所处生境的关系,进一步评价人类活动对海洋保护区微生物群落的影响,以期为我国海洋生态和海洋生物多样性保护提供技术支撑。
大连长海珍稀海洋生物省级自然保护区位于辽宁省大连市长海县小长山镇东南部(122˚44'02''~122˚45'40''E, 39˚12'49''~39˚13'34''N)。保护区面积220公顷,其中海域面积192.25公顷,土地面积27.75公顷。长海县位于中纬度地带,属暖温带半湿润季风性气候。四季分明,冬暖夏凉,日照充足,温差较小,年平均气温10℃左右。夏季平均气温25℃,是消夏避暑胜地。保护区位于中国黄海北部,是目前我国唯一的温带海洋珍稀生物自然保护区。保护区位于著名的海洋岛渔场,有近百种鱼类,二十多种贝类,数十种藻类。海参,鲍鱼,海胆,扇贝和其他海珍品在国内外都广为人知。保护区所在的长海县隶属于辽宁省大连市,位于辽东半岛东侧的黄海北部海域,全县共由252个海岛组成,海域面积10,324平方公里,海岸线总长359公里,全县植被茂密,绿树成荫,森林覆盖率达44.4%,是国家级海岛森林公园。
在大连长海珍稀海洋生物省级自然保护区及附近区域设置潮间带采样点10个(图1、表1),其中P01和P02测点位于保护区核心区域核大坨子岛,C01~C04位于相邻保护区的长海县潮间带,D01~D04为大连市区与长海县相邻区域入海河流入海口。2018年9月12~13日,完成上述10个点位采样。每个测点取表层泥样于已灭菌的50 mL离心中,冷冻保存,带回实验室后保存于−80℃超低温冰箱中。
图1. 监测点位图
采样点位 | 经度 | 纬度 |
---|---|---|
C01 | 122.74 | 39.23 |
C02 | 122.73 | 39.25 |
C03 | 122.64 | 39.27 |
C04 | 122.49 | 39.30 |
D01 | 122.33 | 39.38 |
D02 | 122.29 | 39.38 |
D03 | 122.23 | 39.34 |
D04 | 122.10 | 39.18 |
P01 | 122.76 | 39.23 |
P02 | 122.76 | 39.23 |
表1. 各采样点位置经纬度
使用TIANGEN TIANamp Soil DNA Kit试剂盒提取泥样的总DNA,根据规定实验步骤操作 [
微生物OUT的划分、Chao1、Shannon和Simpson多样性指数的计算均使用Mothur软件进行;通过R语言绘制种类组成图和维恩图等。
潮间带沉积物微生物多样性通过Chao1、Shannon和Simpson多样性指数来表示,并根据Coverage计算结果评估测序深度。各指数的计算公式如下所示:
Schao1 = Sobs + n1(n1 − 1)/2(n2 + 1);其中Schao1为估计的OUT数量,Sobs为实际观察的OUT数量,n1为只含有一条序列的OUT数量,n2为只含有二条序列的OUT数量。
H ′ = − ∑ i = 1 s ( n i N ) L o g 2 ( n i N )
式中:S为观察到的OUT数量,ni为只含有i条序列的OUT数量,N为样品中所有的序列数。
C = 1 − n 1 / N
其中C为Coverage指数,n1为只含有一条序列的OUT数量,N为样品中所有的序列数。
D s i m p s o n = ∑ i = 1 S o b s n i ( n i − 1 ) N ( N − 1 )
其中Sobs为实际观察的OUT数量,ni为只含有i条序列的OUT数量,N为样品中所有的序列数。
本调查共获得大连长海珍稀海洋生物保护区及相邻海域10个站位潮间带微生物原始序列344,985对,按照2.2步骤,筛选出高质量序列330,369条。原始序列通过预处理,按97%序列相似度进行OTU划分,其详细分布见表2。
站位 | 原始序列 | 有效序列 | 优化序列 | 可操作分类单位(OTUs) |
---|---|---|---|---|
C01 | 31,369 | 31,099 | 30,072 | 3815 |
C02 | 32,549 | 32,304 | 30,874 | 2041 |
C03 | 34,205 | 33,936 | 32,623 | 3549 |
C04 | 37,058 | 36,728 | 35,452 | 5469 |
D01 | 36,881 | 36,395 | 35,997 | 4659 |
D02 | 31,646 | 31,288 | 29,808 | 4393 |
D03 | 37,919 | 37,455 | 36,183 | 3307 |
D04 | 31,618 | 31,251 | 30,050 | 4434 |
P01 | 33,614 | 33,221 | 32,558 | 4667 |
P02 | 38,126 | 37,744 | 36,752 | 5623 |
表2. 沉积物微生物序列信息统计
随机抽取经过优化的序列,用抽到的序列数与它们的Chao1指数构建稀释曲线。由图2可知,各站位曲线均趋向平坦,说明本实验测序数量合理,OTU取样已达到或接近饱和,调查结果基本反映了研究海域潮间带微生物的多样性。
图2. 微生物稀释曲线
本调查在大连市长海县附近区域潮间带沉积物中共鉴定微生物43门302科647属,微生物种类较多。从门水平上看,变形菌门(Proteobacteria)是第一优势菌,其相对丰度占总丰度的45.7%。其它较为丰富的门类包括酸杆菌门(Acidobacteria),拟杆菌门(Bacteroidetes)、放线菌门(Actinobacteria)、浮霉菌门(Planctomycetes)、蓝细菌(Cyanobacteria)和疣微菌门(Verrucomicrobia)。调查区域微生物门类组成详见图3。从属水平上看,Gammar变形菌在各站位中是最主要的优势菌,另外常见种类还有黄杆菌属(Flavobacterium)、黏球菌属(Myxococcus)和脱硫杆菌科(Desulfobacteraceae)、固氮菌科(Rhodobacteraceae)的未明确分类地位的属以及某些变形菌、酸杆菌、浮霉菌的属等。各站位主要种属组成详见图4。
海洋沉积物是地球上最复杂的微生物栖息地,海洋沉积物中微生物在沉积环境和上覆水生态系统的能量与物质流动以及营养循环中都扮演着重要的角色 [
本研究中不同区域微生物群落组成的差异情况见图5 (门水平)和图6 (属水平)。将样品分成保护区(P组)、长海县(C组)和大连入海河流(D组)三组,从门水平上看,各区域微生物群落非常相似,总共43个
图3. 潮间带微生物种类组成(门水平)
图4. 潮间带微生物种类组成(属水平)
门类中有32个是三组共有物种,占74.4%,其中D组的物种完全包含了P组,而C组也只有一个门类为其独有物种。从属水平上看,各区域微生物总共647个属中有268个是三组共有物种,占41.4%,其中P组的独有物种最少,只有6个属;C组次之,有47属;D组最多有145属(图3~7)。由此可见,在较高的分类等阶上(门水平)保护区和相邻海域的微生物群落差异并不大。
图5. 微生物群落差异分析(门水平)
图6. 微生物群落差异分析(属水平)
大连市长海县附近区域潮间带沉积物微生物多样性指数见表3。从表3中可以看出,各站位Chao1指数较高,均值为5493,不过各站位间数值波动很大(最大值为7227,最小值为2806)。Chao1在生态学中常用来估计物种总数,由Chao (1984)最早提出。而反映菌群多样性的Shannon指数平均为9.71,Simpson指数平均为0.994,显示微生物多样性非常高。各站位测序深度指数coverage平均为0.916,说明各样品文库的覆盖率较高,样本中序列没有被测出的概率极低。
Chao1指数 | Coverage指数 | Shannon指数 | Simpson指数 | |
---|---|---|---|---|
C01 | 5296 | 0.92 | 9.35 | 0.99 |
C02 | 2806 | 0.96 | 7.69 | 0.98 |
C03 | 4615 | 0.92 | 9.2 | 0.99 |
C04 | 7227 | 0.89 | 10.52 | 1 |
D01 | 6313 | 0.91 | 10.11 | 1 |
D02 | 5539 | 0.91 | 9.81 | 0.99 |
D03 | 4427 | 0.93 | 8.74 | 0.99 |
D04 | 5489 | 0.92 | 10.41 | 1 |
P01 | 5956 | 0.91 | 10.41 | 1 |
P02 | 7262 | 0.89 | 10.82 | 1 |
表3. 沉积物微生物多样性指数
PCA分析(主成分分析),是一种使用最广泛的数据降维算法。它能将大量相关变量转化为一组很少的不相关变量,这些无关变量称为主成分。在数据压缩消除冗余和数据噪音消除等领域都有广泛的应用。通过分析不同样品OTU组成可以反映样品间的差异和距离。利用R语言Scatterplot 3D软件包绘制微生物群落PCA3D图(图7),各个测点在图中的位置较好地反映了采样点地理环境差异。其中C组和P组各测点非常接近,D组与它们相距较远且D组内各点之间距离也较远。说明大连入海河流(D组)的潮间带微生物群落与保护区(P组)及长海县(C组)与差异较大,且组内各点也相距较远,说明各入海河流土壤微生物群落差异较大。沉积物的分布具有区域性特点,沉积环境的不同可能是导致微生物多样性差异的重要原因。沉积物的环境和菌群的类型之间是互相选择的,沉积物样品的微生物群落组成存在差别,这种差别可能与地理位置相关。很多实验证明,由于沉积物样品和土壤样品的复杂性,其中的菌群多样性远远高于水体等其他环境样品中微生物的多样性,沉积物和土壤样品中菌群微生物多样性所能够达到的范围一直以来备受争议 [
1) 大连市长海县及附近区域潮间带沉积物中共鉴定微生物43门302科647属,微生物种类较多。变形菌门(Proteobacteria)是第一优势种类,其相对丰度占总数的45.7%。其它主要门类包括酸杆菌门,拟杆菌门、放线菌门、浮霉菌门、蓝细菌和疣微菌门等。
图7. 潮间带沉积物微生物群落PCA图(属水平)
2) 各站微生物多样性非常高,其中Chao1指数较高,均值为5493,Shannon指数平均为9.71,Simpson指数平均接近1。
3) 通过主成份分析发现,保护区(P组)和长海县(C组)各测点微生物群落结构非常接近,而大连入海河流(D组)则差异较大,且D组内各点之间也各不相同。
浙江省生态环境厅科研项目(2016A012;2017A011;2018A022及2019A001)。
唐嘉威,黄 备,孙沛雯,孙明霞,刘子成,李翔宇,俞存根. 海洋保护区潮间带沉积物微生物群落的研究——以大连长海珍稀海洋生物保护区为例Study on Microbial Communities in the Intertidal Sediments of Marine Protected Areas—A Case Study on Dalian Changhai Marine Precious Biological Nature[J]. 海洋科学前沿, 2020, 07(01): 13-23. https://doi.org/10.12677/AMS.2020.71003