張爽，晏磊，陳志寶

（黑龍江八一農(nóng)墾大學生命科學技術(shù)學院，大慶 163319）

硫氧化菌的地域分布特征

張爽，晏磊，陳志寶

（黑龍江八一農(nóng)墾大學生命科學技術(shù)學院，大慶 163319）

通過對NCBI數(shù)據(jù)庫中硫氧化菌的現(xiàn)有16S rDNA序列分析，并構(gòu)建硫氧化菌系統(tǒng)進化樹，進而探討不同地域硫氧化菌分布特性。NCBI中共有硫氧化菌16S rDNA序列298條，經(jīng)相似性分析后發(fā)現(xiàn)共有53條不同的序列，分別源自溫泉湖、礦區(qū)、瀉湖等區(qū)域。全長序列的系統(tǒng)進化分析顯示，全球硫氧化菌的分布具有一定的地域性，即溫泉湖與礦區(qū)、瀉湖的硫氧化菌差別明顯。不同地域硫氧化菌的系統(tǒng)進化分析顯示，鉛鋅礦區(qū)與銅礦區(qū)的硫氧化菌區(qū)分明顯；而同為溫泉湖或礦區(qū)、瀉湖生境中的硫氧化菌相似性極高，即相同生境中的硫氧化菌相似性較高，表明硫氧化菌的類型與生境條件有較大的關(guān)聯(lián)性。基于對硫氧化菌16S rDNA的序列分析，推測硫氧化菌的起源可能是多源的，生境條件作為重要的進化壓力在硫氧化菌長期演化過程中具有舉足輕重的作用。

硫氧化菌；16S rDNA序列分析；系統(tǒng)發(fā)育；地域分布

硫氧化菌（Sulfur oxidizing bacteria，SOB）是一類通過氧化還原性硫化物獲取能源生長的特殊微生物，在pH 1.0～9.0的環(huán)境中均能生長[1]。SOB在自然界的分布極為廣泛，從火山熔巖到礦山廢水[2]，從海洋底泥到鹽堿湖泊[3]，從陸地熱泉到油藏鹵水[4-5]，從根際土壤到地表河流[6]，即使在金屬離子濃度極高的區(qū)域[7]，都分布了SOB的一些種類，是極端環(huán)境中占優(yōu)勢的古老微生物之一。由于SOB在硫的地球化學循環(huán)中的重要作用，以及在生物浸礦和環(huán)境治理方面的應用，使其成為當前地球科學、環(huán)境科學以及生物科學等領(lǐng)域的研究熱點。

目前為止，從NCBI信息庫中搜索到的硫氧化菌序列共298條，但由于硫氧化菌在培養(yǎng)過程中對營養(yǎng)及氧的要求比較苛刻，實驗室環(huán)境難以重現(xiàn)其真實生存環(huán)境，因此獲得的純培養(yǎng)硫氧化菌種類較少。16S rDNA在結(jié)構(gòu)與功能上具有極高的保守性，因此常被用作原核生物系統(tǒng)發(fā)育以及多樣性研究的標志性基因，用表型方法難以鑒定的細菌，用16S rDNA序列分析均能進行準確的鑒定[8]。由于目前獲得的硫氧化菌純培養(yǎng)菌株種類較少，使得多數(shù)未純培養(yǎng)的硫氧化菌通過分析16S rDNA序列來確定它們的系統(tǒng)進化位置。近年來，由于硫氧化菌在分子形態(tài)、生理活性以及系統(tǒng)發(fā)育等方面具有多樣性，硫氧化菌的生態(tài)學意義被廣泛關(guān)注。研究旨在通過對來自NCBI中298條不同地域和生境及相同地域和生境的硫氧化菌16S rDNA序列進行的分析，及對溫泉湖、礦區(qū)、瀉湖等已發(fā)現(xiàn)的硫氧化菌系統(tǒng)進化樹的構(gòu)建，初步探討硫氧化菌的地域分布特征。

1 數(shù)據(jù)來源與分析方法

1.1 數(shù)據(jù)來源

研究所使用的數(shù)據(jù)均來自截至2014年3月在NCBI信息庫中搜索，獲得硫氧化菌298條16S rDNA。搜索關(guān)鍵詞包括Sulfobacillus acidophilus，Sulfobacillus thermotolerans，Sulfolobus tokodaii等，盡可能涵蓋所有硫氧化菌16S rDNA序列，并對獲得的硫氧化菌16S rDNA序列進行必要的篩選，刪除非硫氧化菌的16S rDNA序列。

1.2 數(shù)據(jù)整理與分析方法

將篩選后的硫氧化菌 16S rDNA序列進行BLAST比對（http：//www.ncbi.nlm.nih.gov／blast），每條序列即可獲得與其相似性最高的序列，再利用clustal1.83軟件對得到的硫氧化菌16S rDNA序列進行相似性分析。Stackebrandt等[9]認為，當16S rDNA相似性高于97%以上時，即可判定為同屬于一個種。如若同一個地點、同一次提交的序列中有2條序列的相似性均高于97%，即可判定兩者為同一種硫氧化菌的序列，則保留其中一個序列即可。將所有序列逐一整理，按生境（火山熔巖、礦山廢水、海洋底泥、鹽堿湖泊、陸地熱泉、油藏鹵水、根際土壤、地表河流湖泊和礦區(qū)等）、分類地位、發(fā)現(xiàn)地域和所屬國家、以及序列長度等方面進行歸類并逐一統(tǒng)計。

由于所提交的16S rDNA序列存在長度不等問題，無法直接構(gòu)建其系統(tǒng)進化樹。因此，研究利用硫氧化菌16S rDNA序列的前段（約750 bp）、中段（約600 bp）及后段（約750 bp）分別進行系統(tǒng)進化樹分析，并從中挑選全長序列進行系統(tǒng)進化分析。充分應用MEGA5.0軟件利用N-J法構(gòu)建系統(tǒng)進化樹，從而實現(xiàn)硫氧化菌的系統(tǒng)進化分析。同時，對不同生境（湖泊和礦區(qū)火山熔巖、礦山廢水、海洋底泥、鹽堿湖泊、陸地熱泉、油藏鹵水、根際土壤、地表河流湖泊和礦區(qū)等）的硫氧化菌進行分析，對不同地域的硫氧化菌進行分析，探索其規(guī)律性。

2 結(jié)果與分析

2.1 原始數(shù)據(jù)整理

從NCBI信息庫中共搜索到硫氧化菌16S rDNA序列298條，其中包括湖泊硫氧化菌序列185條，礦區(qū)硫氧化菌序列115條。利用clustal1.83軟件對獲得的硫氧化菌16S rDNA序列進行相似性分析，若相似性高于97%即判定為屬于同一序列[9]；若有2條序列相似性均高于97%，則保留其中一個序列即可。刪除重復序列后所獲得的湖泊、礦區(qū)硫氧化菌序列分別為36條、17條，共計53條，見表1。

表1 硫氧化菌種數(shù)統(tǒng)計（據(jù)提交NCBI的16S rDNA序列數(shù)統(tǒng)計）Table 1 Sulfur oxidizing bacteria statistics（according to the 16S rDNA submitted to the NCBI）

在獲得的53條不同硫氧化菌序列中，γ-變形菌綱序列15條，熱變形菌綱序列14條，綠菌綱序列13條，綠彎菌綱和β-變形菌綱的硫氧化菌序列較少。

研究發(fā)現(xiàn)，硫氧化菌的發(fā)現(xiàn)地點主要集中在中國、印度、智利、伊朗等。各地發(fā)現(xiàn)的硫氧化菌均以熱變形菌綱居多，而中國發(fā)現(xiàn)的以梭狀芽孢桿菌綱居多。

2.2 硫氧化菌系統(tǒng)進化分析

研究篩選后的16S rDNA序列共53條，因獲得的16S rDNA序列長度不等，無法直接構(gòu)建其系統(tǒng)進化樹。故分別選取硫氧化菌16S rDNA序列的前段（約750 bp）、中段（約600 bp）及后段（約750 bp），構(gòu)建較為精確的系統(tǒng)進化樹，并進行比較分析（圖1）。從系統(tǒng)進化樹中不難發(fā)現(xiàn)湖泊（溫泉湖、瀉湖）硫氧化菌與礦區(qū)硫氧化菌聚類分區(qū)明顯，各自分區(qū)聚類在一起，但是由于所屬綱不同，硫氧化菌的分布有所不同。

2.3 不同區(qū)域的硫氧化菌比較

2.3.1 湖泊與礦區(qū)硫氧化菌比較

湖泊與礦區(qū)硫氧化菌的比較（見圖1），由圖1可見，湖泊中以熱變形菌綱和綠菌綱硫氧化菌為主，約為75.0%；而礦區(qū)中的硫氧化菌以γ-變形菌綱為主，約為88.2%。同時，由圖1的進化樹也可以發(fā)現(xiàn)，湖泊與礦區(qū)硫氧化菌相似程度較低，表明系統(tǒng)進化位置具有明顯差別。

圖1 湖泊與礦區(qū)中硫氧化菌系統(tǒng)進化樹Fig.1 Phylogenetic tree based on 16S rDNA of SOB in lake and mine

2.3.2 溫泉湖與瀉湖硫氧化菌比較

溫泉湖（中國廈門海）與瀉湖中硫氧化菌的比較（見圖2），由圖2可見，瀉湖中發(fā)現(xiàn)的硫氧化菌以綠菌綱為主，相較于溫泉湖中發(fā)現(xiàn)的熱變形菌綱和梭狀芽孢綱，其同源性約為78%，表明溫泉湖中的硫氧化菌與瀉湖中的硫氧化菌相似程度較低。溫泉湖中的硫氧化菌具有嗜熱嗜酸的特性，而瀉湖中發(fā)現(xiàn)的硫氧化菌通常能夠耐受較高濃度，推測溫泉湖與瀉湖硫氧化菌區(qū)分明顯可能是其自身特性不同產(chǎn)生的結(jié)果。

圖2 溫泉湖與瀉湖中發(fā)現(xiàn)的硫氧化菌系統(tǒng)進化樹Fig.2 Phylogenetic tree based on 16S rDNA of SOB in hot spring and lagoon

2.3.3 不同地域礦區(qū)硫氧化菌比較

不同地域礦區(qū)硫氧化菌的比較，僅針對中國與智利、秘魯?shù)V區(qū)生境硫氧化菌進行比較（見圖3），由圖3可見，在智利與秘魯?shù)V區(qū)中發(fā)現(xiàn)的硫氧化菌屬于梭狀芽孢桿菌綱，而中國礦區(qū)發(fā)現(xiàn)的硫氧化菌多屬于γ-變形菌綱。因不同地域礦區(qū)硫氧化菌所屬分類地位不同，使得硫氧化菌的相似性極低，各自形成獨立分支。伴隨著研究的不斷深入，不同地域硫氧化菌的多樣性特性將不斷明晰。

圖3 中國與智利發(fā)現(xiàn)的礦區(qū)中硫氧化菌系統(tǒng)進化樹Fig.3 Phylogenetic tree based on 16S rDNA of SOB in mine from China and Chile

3 討論

通過對不同地域硫氧化菌系統(tǒng)進化的分析，不難發(fā)現(xiàn)相近生境的硫氧化菌相似性較高。此外，目前發(fā)現(xiàn)的硫氧化菌多存在于湖泊，并以溫泉湖為主，且多為溫泉湖出口。由此推測硫氧化菌的分布可能具有區(qū)域性特征，其在世界各地的生境中可能以“土著”的方式存在。為了適應生境，不同地域類型硫氧化菌經(jīng)過長期進化被逐漸形成。

對比礦區(qū)與湖泊硫氧化菌，不難發(fā)現(xiàn)溫泉湖硫氧化菌16S rDNA序列大多聚集在一起，與礦區(qū)或瀉湖中發(fā)現(xiàn)的硫氧化菌形成明顯分支差異。存在于相似生境的硫氧化菌相似性程度較高，究其原因可能是由于相似生境具有相似條件（如溫度、酸度、鹽度及其他種類元素等），而生境又是推動硫氧化菌進化的主動力，因此促使硫氧化菌向同一方向進化，最終形成相似性較高的硫氧化菌。反之，不同生境的硫氧菌經(jīng)長期進化形成不同類型的硫氧化菌。

縱使相似生境中發(fā)現(xiàn)的硫氧化菌相似性較高，但也存在特例。研究發(fā)現(xiàn)礦石浸出物中的硫氧化菌（EF101930.1）與在溫泉湖發(fā)現(xiàn)的硫氧化菌（EU419196.1），經(jīng)軟件分析發(fā)現(xiàn)兩者相似性高達92%。換言之，在不同地域生境中發(fā)現(xiàn)了2株16S rDNA序列比較接近的硫氧化菌?？赡苁怯捎谠惊毩⒋嬖诘牧蜓趸谶M化過程中被趨同進化，進而形成相似性極高的硫氧化菌類型。這2株相似性極高的硫氧化菌可能是進化過程的偶然現(xiàn)象，也有可能是生境被作為進化壓力的某些環(huán)境因子（如溫度、酸度、鹽度及其他種類元素等）發(fā)揮了作用，導致硫氧化菌朝同一方向進化。僅通過數(shù)據(jù)分析并不能完全認定這2株相似性較高的硫氧化菌為功能相同的一類原核微生物。據(jù)張文燕等[8，10]報道，引起痢疾的病原菌與普通大腸桿菌的16S rDNA相似性高于97%，系統(tǒng)進化分析認為二者屬于同一類細菌，但實則二者生理功能不同，且功能迥異。因此，鑒于具有特定功能的一類原核微生物而言，如硫氧化菌，功能基因或全基因組的分類分析較16S rDNA構(gòu)建的系統(tǒng)進化樹意義更為深遠。

研究受NCBI數(shù)據(jù)庫硫氧化菌16S rDNA序列數(shù)量的限制，故有些分析結(jié)果還需更深入的研究來證實。伴隨著硫氧化菌研究的不斷深入及微生物分離及培養(yǎng)技術(shù)的不斷完善，更多、更純的硫氧化菌將被獲得。同時，隨著種類繁多、數(shù)量龐大的硫氧化菌的不斷發(fā)現(xiàn)，輔以分子和基因水平的深入研究，必將獲得不同地域有關(guān)于硫氧化菌更多、更全的生態(tài)信息，為更深入的認識和理解硫氧化菌生境意義和應用前景打下更為殷實的基礎(chǔ)。

4 結(jié)論

利用NCBI數(shù)據(jù)庫分析硫氧化菌16S rDNA序列，系統(tǒng)進化分析發(fā)現(xiàn)處于不同地域的硫氧化菌具有明顯的地域分布特性，如湖泊（溫泉湖、瀉湖）與礦區(qū)之間的硫氧化菌區(qū)分明顯。而處于相同生境的硫氧化菌，如同處于湖泊（溫泉湖、瀉湖）或礦區(qū)，則相似性極高，充分說明硫氧化菌類型與生境之間具有密切聯(lián)系。這些發(fā)現(xiàn)均驗證生境作為進化壓力在硫氧化菌長期進化過程中決定其向同一方向進化（相似環(huán)生境形成相似硫氧化菌類型）或多元化方向進化（不同生境形成不同硫氧化菌類型），由此推測硫氧化菌的起源可能是多源的。

［1］ Biderre-Petit C，Boucher D，Kuever J，et al.Identification of sulfur-cycle prokaryotes in a low-sulfate lake（Lake Pavin）using aprA and 16S rRNA gene markers［J］.Microbial Ecology，2010，61（2）：313-327.

［2］ Ghosh W，Dam B.Biochemistry and molecular biology of lithotrophic sulfur oxidation by taxonomically and ecologically diverse bacteria and archaea［J］.FEMS Microbiology Reviews，2009，33（6）：999-1043.

［3］ Jung M Y，Pham V H，Park S J，et al.Metagenomic assessment of a sulfur-oxidizing enrichment culture derived from marine sediment［J］.The Journal of Microbiology，2010，48（6）：739-747.

［4］ Wei L，Ma F，Zhao G.Composition and dynamics of sulfate-reducing bacteria during the waterflooding process in the oil field application［J］.Bioresource Technology，2010，101（8）：2643-2650.

［5］ Ogg C D，Patel B K C.Desulfotomaculum varum sp.nov.，a moderately thermophilic sulfate-reducing bacterium isolated from a microbial mat colonizing a Great Artesian Basin bore well runoff channel.3［J］.Biotech，2011（3）：139-149.

［6］ Liu X Z，Zhang L M，Prosser J I，et al.Abundance and community structure of sulfate reducing prokaryotes in a paddy soil of southern China under different fertilization regimes［J］.Soil Biology and Biochemistry，2009，41（4）：687-694.

［7］ Ghosh W，Dam B.Biochemistry and molecular biology of lithotrophic sulfur oxidation by taxonomically and ecologically diverse bacteria and archaea［J］.FEMS Microbiology Reviews，2009，33（6）：999-1043.

［8］ 張文燕，張圣妲，肖天，等.趨磁細菌的地域分布特性［J］.環(huán)境科學，2010，31（2）：450-458.

［9］ 石艷玲，鄭偉.景觀水體微生態(tài)調(diào)控試驗研究——以日照職業(yè)技術(shù)學院樂學河生態(tài)修復為例［J］.長春大學學報，2014，24（4）：468-472.

［10］ 王長虹，王國興，晏磊，等.寒區(qū)條垛式和槽式堆肥工藝的比較研究［J］.黑龍江八一農(nóng)墾大學學報，2016，28（1）：68-72.

Geographical Distribution of Sulfur-oxidizing Bacteria

Zhang Shuang，Yan Lei，Chen Zhibao
（College of Life Science and Technology，Heilongjiang Bayi Agricultural University，Daqing 163319）

Sulfur-oxidizing bacteria（SOB）were widely distributed in aquatic environments.To investigate the correlation between their evolutionary relatedness and geographical distribution，298 16S rDNA sequence available in the NCBI were analyzed，and phylogenic trees based on the sequences were constructed.Total 53 16S rDNA sequences of SOB were analyzed in detail，which derived from hot spring，mine and lagoon.Phylogeny analysis based on those sequences suggested that the geographical distribution of SOB had certain regional distribution character：SOB was different from hot spring，mine and lagoon.Regional analysis based on phylogenic analysis also suggested that the difference from lead-zinc mine and copper mine.In contrast，the SOB from similar habitats were closely related，such as hot spring，mine and lagoon.It was found that similar species had a large regional distribution and tried to adopt the similar habitats，morphotypes of SOB and their living conditions had a significant relevance.The observation suggested that SOB might have multiple evolutionary origins.And it also suggested that the environmental conditions，as an important evolutionary pressure，played an important role in the long-term evolution of SOB.

Sulfur oxidizing bacteria（SOB）；16S rDNA sequence analysis；phylogeny analysis；geographical distribution

Q938

A

1002-2090（2017）02-0068-06

10.3969/j.issn.1002-2090.2017.02.014

2016-03-29

張爽（1981-），女，講師，黑龍江八一農(nóng)墾大學畢業(yè)，現(xiàn)主要從事微生物藥物分析方面的教學與研究工作。