阿蓬江黔江段干支流水體微生物群落特征分析

摘 要 以阿蓬江（黔江段）為研究區域，通過對研究區域水體微生物的組成成分和群落特征的研究，為長江上游水環境保護和長江流域微生物群落特征和多樣性研究提供數據支持。結果表明：通過微生物16S片段高通量測序數據，可知阿蓬江（黔江段）水樣樣品中包含27門、46綱、98目、158科。在門水平分類上，以變形菌門占絕對優勢，在綱分類上則是以γ-變形菌綱和放線菌綱為優勢類群，在目分類上以伯克氏菌目為優勢類群，在科分類上以伯克氏菌科占據絕對優勢?；诟咄繙y序數據分析，干支流空間差異較為明顯。

關鍵詞 微生物； 群落特征；阿蓬江

中圖分類號：Q89 文獻標志碼：A DOI：10.19415/j.cnki.1673-890x.2024.03.002

我國水體微生物研究主要集中在長江[1]、黃河[2]等主要河流地，以及鄱陽湖[3]等典型湖泊，對于支流等小型河流、湖泊及其他小型水域的微生物缺乏研究。河流生態系統中，浮游細菌是重要組成部分，在水體中浮游細菌數量龐大，代謝活動強烈，是高級生物的主要食物，對水中污染物的降解和水循環過程有重要作用[4]。在水生態系統中，微生物對水體外的環境因子的變化十分敏感。近年來，眾多學者研究了浮游細菌與環境因子之間的相關性 [5]，但大都集中在海洋和湖泊[6]。河流作為一個動態系統，它們的環境條件在時間和空間尺度上具有高度差異性[7]。相關研究[8]表明其在濕地、養殖廢水和深海沉積物等環境中優勢類群均是變形菌門，且大多數在生物脫氮除磷和污染物降解過程中起重要作用的微生物也是變形菌門[9]。在綱水平上，β-變形菌綱占據的比例最大，達到了11.01%到31.50%不等，其大多數存在于河流和湖泊中（松花江[10]、海寧長山河[11]、東江[12]、Prealpine湖[13]等），對氮、磷等污染物的降解起到了重要的作用[14]。

本研究選擇阿蓬江（黔江段）為研究區域，對其水體中的微生物種類、多樣性和群落分布特征進行分析，以期為觀察阿蓬江水污染狀況，為長江上游水環境保護和長江流域微生物群落特征和多樣性研究提供數據支持。

1" 材料與方法

1.1" 研究區概況及樣品采集

阿蓬江（黔江段）位于重慶市黔江區南部，全長約60 km。主要支流有太極河、細沙河、段溪河及黔江河等河流。于2019年在阿蓬江（黔江段）干支流中設定11個樣點并取樣，其中支流設有4個樣點，分別是QW1、QW7、QW9、QW11樣點，其余7個全是干流樣點。

1.2" 水體微生物群落結構分析

1.2.1" 水體微生物總DNA提取

水樣通過0.22 μm的濾膜進行抽濾，過濾好后的濾膜放置在無菌離心管中于-20 ℃保存，并使用15 kg干冰送測。

1.2.2" 水體微生物高通量DNA測序

用特定的DNA提取試劑盒進行基因組DNA提取后，再用0.8%瓊脂糖凝膠電泳檢測DNA。16S引物使用：341F（5′-CCTAYGGGRBGCASCAG-3′） 和 806R（5′-GGACTACHVGGGTWTCTAAT-3′）以稀釋后的基因組 DNA 為模板，根據測序區域的選擇，使用帶 Barcode 的特異引物進行PCR。 每一個25 μL的體系包括 1×PCR buffer、1.5 mM MgCl2、0.4 μM dNTPs、正向和反向引物各1.0 μM、0.5 U KOD-Plus-Neo 酶（TOYOBO） 和 10 ng 模板。 PCR程序包括起始 94 ℃ 1 min， 然后 30 循環 （變性 94 ℃ 20 s，退火 54 ℃ 30 s 和 延伸 72 ℃ 30 s），最后 72 ℃ 5 min。 每個樣本進行3個PCR技術重復。PCR 與 1/6 體積的 6×loading buffer 混合，用 2% 瓊脂糖凝膠電泳檢測。取目的條帶用來回收，回收使用 QIAquick Gel Extraction Kit （QIAGEN）。用Qubit@ 2.0 Fluorometer （Thermo Scientific）定量。最后等摩爾量混合。建庫使用 TruSeq DNA PCR-Free Sample Prep Kit，構建好的文庫經過定量和文庫檢測合格后，使用羅寧生物的Miseq平臺PE300模式測序。測序工作在上海羅寧生物有限公司進行。

1.3" 數據分析與處理

高通量測序產出的原始圖像數據文件經CASAVA軟件進行堿基識別后轉化為原始的測序序列，稱為Raw Data或Raw Reads，結果以Fastq格式儲存。對序列進行質控，得到高質量序列，質控后得到的序列稱為Clean Tag。使用UCHIME7算法去除嵌合體，得到Effective Tag。為了方便研究樣本的物種組成及多樣性信息，對序列進行聚類產生操作時的分類單元，即OTU（Operational Taxonomic Units）。在97%的相似水平下進行OTU聚類分析， 每個OTU代表的是一類相似序列的集合，通過OTU聚類分析，得到OTU豐度表，此豐度表是后續分析的基礎文件。為防止測序深度對后續多樣性分析產生影響，對數據進行重抽樣，使得各樣本保持一致的序列數。用Circos網站在線繪圖功能（http：//circos.ca/）對阿蓬江水體微生物種群門、綱、目、科水平的相對豐度變化情況進行可視化表達。

2" 結果與分析

2.1" 阿蓬江干支流水體微生物組成成分

在阿蓬江（黔江段）區域所取樣的11個樣點所得的所有細菌16 SrDNA片段中，將得到的序列按門、綱、目、科分類，可以知道研究區域所取樣的11個樣點共有27門、46綱、98目、158科。

2.2" 阿蓬江干支流水體微生物群落多樣性

通過采用α多樣性分析，對阿蓬江干支流水體微生物群落多樣性進行表征，包括物種豐富度指數（Chao 1）、Shannon多樣性指數（H’）和Simpson多樣性指數（D）（見表1）。對阿蓬江干支流水體樣品測序后，獲得97%相似水平及以上的OTUs分別有約360個和337個，且對微生物群落覆蓋率達到99%，OTU值干流明顯大于支流。Chao 1豐富度指數表現為干流水體大于支流水體，干支流微生物群落OTU值均低于Chao 1豐富度指數， Shannon和Simpson指數也均表現為干流水體大于支流水體。總體而言，阿蓬江水體微生物的群落豐富度和多樣性均表現為干流大于支流。

2.3" 門水平上的微生物群落分布特征

通過高通量測序分析，發現細菌群落在門水平上有27個細菌門，其中主要為變形菌門、放線菌門和擬桿菌門等細菌門，而在這些細菌門中占據絕對優勢的是變形菌門，說明阿蓬江（黔江段）是以浮游細菌組成的典型的淡水水體，符合我國典型的淡水水體微生物組成的群落特征。

從表2可知，在干流中的優勢細菌門是放線菌門、擬桿菌門、厚壁菌門、浮霉菌門和裝甲菌門，支流中的優勢細菌門則是變形菌門、疣微菌門、異常球菌門、芽單胞菌門和奇古菌門，從標準偏差來看，放線菌門、擬桿菌門和異常球菌門在干支流的優劣勢并不明顯，擬桿菌門和厚壁菌門雖然在占比上是干流優于支流，但是其干流的標準偏差是遠大于支流的，相反的是，異常球菌門在占比上是支流優于干流，其標準偏差是支流略大于干流，由此可看出細菌群落在空間的分布上是有一定的空間差異，但其波動幅度不大。

2.4" 綱水平上的微生物群落結構分析

通過高通量測序分析，在綱水平上一共發現46個細菌綱。 從表3的干支流占比分析來看，干流中的優勢細菌綱有放線菌綱和α-變形菌綱，在支流中占據優勢的細菌綱是γ-變形菌綱和擬桿菌綱。從標準差來看，首先γ-變形菌綱在支流中的標準差是要小于干流，其次是放線菌綱、擬桿菌綱和α-變形菌綱，其標準差均是干流大于支流，最后是疣微菌綱和Acidimicrobiia，這兩大細菌綱在干支流中的標準差相距太大，疣微菌綱在支流中的標準差是干流中的4倍，Acidimicrobiia在干流的標準差是支流的3倍，所以這兩大細菌綱在干支流中的空間分布有一定的波動幅度。

2.5" 目水平上的微生物群落結構分析

目水平上，一共發現98個細菌目，如表4所示，最優勢的細菌目是伯克氏菌目，從干支流占比對比分析，伯克氏菌目在支流中的比例大于干流中所占比例，法蘭克氏菌目干流中所占比例大于支流，從兩者標準差來看，干支流中的空間分布沒有明顯差異。擬桿菌目和鞘脂單胞菌目在干流中的占比略優于支流，而其在干流中的空間分布波動相比于支流是比較大的，而剩下的假單胞菌目、噬幾丁質菌目、熱纖梭菌目等細菌目的相對比例差距不大，空間波動幅度不明顯。

2.6" 科水平上的微生物群落結構分析

科水平上，細菌科較多，所有樣品中一共有158個細菌科，如表5所示，最優勢的細菌科是伯克氏菌科。從干支流的占比來看，支流中占比較大的是伯克氏菌科，而在干流中所占比例較大的是魚孢菌科。從標準差的對比來看，伯克氏菌科和鞘脂單胞菌科在干流中的空間波動幅度遠大于支流，魚孢菌科、紅桿菌科和假單胞菌科在干支流中的空間分布沒有明顯差異，而未列出的細菌科占據比例較大，在空間分布上也沒有明顯差異。

3" 討論與結論

3.1" 討論

利用Illumina測序技術的16S擴增子測序，檢測了阿蓬江（黔江段）水體微生物群落的組成成分，發現所取樣的11個樣點共有27門、46綱、98目、158科。在門分類水平上，變形菌門是11個樣品中最優勢細菌門，如吳曉冰等[15]在研究長江干流浮游細菌群落結構中變形菌門（7.49%～85.63%）和張旭[16]在研究海寧市不同類型區域河流底泥微生物群落結構特征研究時，其中變形菌門為細菌門中最優勢類群，均符合典型的淡水水體微生物的組成特征。變形菌門在研究區域呈現出支流含量高、干流含量低的特征，變形菌門在所有樣點中都是絕對優勢類群，說明變形菌門相比其他細菌更具有適應性；放線菌門和擬桿菌門作為次優勢類群，在研究區域呈現出干流含量高、支流含量低的特征；而剩下的厚壁菌門、疣微菌門及浮霉菌門等細菌門的占比相對較低，在10%以下。在綱水平上，γ-變形菌綱是最優勢類群，在研究區域呈現出支流含量高、干流含量低的特征；其次次優勢類群是放線菌綱、擬桿菌綱及α-變形菌綱，在研究區域呈現出干流含量高、支流含量低的特征；其余的疣微菌綱、厚壁菌綱及β-變形菌綱等細菌綱的占比均不超過10%。γ-變形菌綱、放線菌綱、擬桿菌綱及α-變形菌綱在研究區域是優勢類群，這與我國大多數淡水生態系統的調查結果是一樣的[17-19]。在目和科水平中伯克氏菌均處于優勢地位，其干支流的占比差距明顯。初步分析伯克氏菌在支流上含量高的原因可能是支流采樣點都是河流，其周圍城鎮發展迅速，其發展產生的各種污水導致水體污染加劇，其水體內營養物質增多，細菌含量增高。

3.2" 結論

1）通過微生物16S片段高通量測序分析，阿蓬江（黔江段）所取得的11樣品在界的水平上總共發現27門、46綱、98目、158科。

2）阿蓬江水體微生物的群落豐富度和多樣性均表現為干流大于支流。

3）在門水平分類上，以變形菌門為絕對優勢類群，在綱水平分類上則是以γ-變形菌綱和放線菌綱為優勢類群，在目分類上以伯克氏菌目為優勢類群，在科分類上以伯克氏菌科占據絕對優勢。

4）基于高通量測序數據分析，干支流空間差異較為明顯，在門水平上，擬桿菌門和厚壁菌門在干流中的空間分布波動幅度較大，疣微菌門和芽單胞菌門在支流中的空間分布有較大差異。綱水平上，疣微菌綱在支流的標準差是干流的4倍，Acidimicrobiia在干流的標準差是支流的3倍，所以這兩大細菌綱在干支流中的空間分布有一定的差異。在目分類上，細菌目的空間分布有一定的差異，但熱纖梭菌目、Microtrichales和Selenomonadales在分布上差異較大。在科分類上，細菌科的空間分布沒有明顯的差異，其中的鞘脂單胞菌科、Ilumatobacteraceae和氨基酸球菌科的空間分布有些波動。

參考文獻：

[1] WOESE C R，FOX G E.Phylogenetic structure of the prokaryotic domain： The primary kingdoms[J]. Proceedings of the National Academy of Sciences， 1977， 74：5088-5090.

[2] KUCZYNSK J. Experimental and analytical tools for studying the human microbiome[J]. Nat Rev Genet， 2012，13：47-58.

[3] 余晨笛.長江河口沉積物全程氨氧化微生物的富集及多樣性分析[D].上海：華東師范大學，2018.

[4] 李明聰.黃河入?？谖⑸矧寗拥牡h過程及其影響因素的研究[D].泰安：山東農業大學，2019.

[5] 張婷.鄱陽湖河湖交錯帶水體微生物群落特征研究[D].南昌：江西師范大學，2018.

[6] 滿江紅，王先良，楊永堅，等.我國地表水環境微生物基準研究現狀[J].環境與健康雜志，2012，29（1）：82-84.

[7] STUART FINDLAY. Stream microbial ecology[J]. Journal of the North American Benthological Society，2010， 29（1）：170-181.

[8] POMEROY L R， LEB WILLIAMS P J， AZAM F， et al. The microbial loop[J]. Oceanography， 2007， 20（2）：28-33.

[9] R JAN STEVENSON， JOHN P SMOL. Use of Algae in Environmental Assessments[P]. Department of Zoology Michigan State University East Lansing， Michigan 48824， USA ， 2003.

[10] JIANG J G， SHEN Y F. Development of the microbial communities in Lake Donghu in relation to water quality[J]. Environmental Monitoring and Assessment， 2007， 127（1-3）：227-236.

[11] LAQUE T， FARJALLA V F， RODADO A S， et al. Spatiotemporal variation of bacterial community composition and possible controlling factors in tropical shallow lagoons[J]. Microbial Ecology， 2010， 59（4）：819-829

[12] DEL Giorgio P A， Cole J J. Bacterial growth efficiency in natural aquatic systems[J]. Annual Review of Ecology and Systematics， 1998， 29：503-541.

[13] PING XU， LAURA G LEFF. Longitudinal changes in the benthic bacterial community of the Mahoning River （Ohio， USA.）[J]. Hydrobiologia， 2004， 522（1-3）：329-335.

[14] LEFF L G. Stream bacterial ecology： A neglected field[J]. Am Soc Microbiol News， 1994， 60（3）：135-138.

[15] 吳曉冰，葉飛，姜毅，等.長江干流浮游細菌群落結構及影響因素[J].長江流域資源與環境，2019， 28（7）：1652-1662.

[16] 張旭.淡水水體底泥微生物群落特征及微生物指標評價研究[D].上海：上海大學，2016.

[17] 朱珍香，高肖飛，彭鳳，等.廈門后溪水質與流域景觀特征沿城鄉梯度的變化分析[J].生態學報， 2019，39（6）：2021-2033.

[18] 王鵬，陳波，李傳瓊，等.贛江南昌段豐水期細菌群落特征[J].中國環境科學，2016，36（8）：2453-2462.

[19] JORDAAN K， BEZUIDENHOUT C C. Bacterial community composition of an urban river in the North West Province， South Africa， in relation to physico-chemical water quality[J]. Environmental Science and Pollution Research， 2016， 23（6）：5868-5880.

（責任編輯：敬廷桃）

收稿日期：2023-07-31

作者簡介：許田甜（1983—），本科，工程師，主要研究方向為生態學、測繪測量。E-mail：29462830@qq.com。

*為通信作者，E-mail：7426997@qq.com。