利用小亞基核糖體RNA技術分析溫室黃瓜近根土壤古菌和真菌多樣性

趙志祥，蘆曉飛，陳國華，茆振川，楊宇紅，劉二明，謝丙炎

1 中國農業科學院蔬菜花卉研究所，北京 100081

2 湖南農業大學生物安全科技學院，長沙 410128

利用小亞基核糖體RNA技術分析溫室黃瓜近根土壤古菌和真菌多樣性

趙志祥1,2，蘆曉飛1，陳國華1，茆振川1，楊宇紅1，劉二明2，謝丙炎1

1 中國農業科學院蔬菜花卉研究所，北京 100081

2 湖南農業大學生物安全科技學院，長沙 410128

土壤古菌和真菌在溫室生態系統是僅次于細菌的微生物，具有類似于細菌的重要生態功能。通過構建古菌16S rRNA和真菌18S rRNA基因克隆文庫，分析溫室黃瓜近根土壤古菌和真菌群落結構組成，為開發利用溫室這一特殊的生態環境中豐富的微生物資源以及理解微生物與植物間的互作提供參考依據。采用研磨-凍融-溶菌酶-蛋白酶K-SDS熱處理以及CTAB處理等理化方法，提取和純化微生物總DNA，構建古菌16S rRNA和真菌18S rRNA基因克隆文庫。利用DOTUR軟件將古菌和真菌序列按照相似性97%的標準分成若干個可操作分類單元 (OTUs)。土壤古菌克隆文庫主要包括泉古菌門和未分類的古菌兩大類，并有少部分廣域古菌類群，所有泉古菌均屬于熱變形菌綱，共45個OTUs；真菌克隆文庫包括真菌門的大多數亞門真菌，共24個OTUs，未發現擔子菌亞門真菌。古菌多樣性比較豐富，且發現少量的廣域古菌 (甲烷菌)，這一情況可能與溫室長期高溫高濕，高有機質含量，土壤處于缺氧環境有關；土壤真菌的優勢種群為子囊菌，占到土壤真菌的80%以上，這可能與絕大多數植物真菌性病害屬于土傳病害，通過菌絲體、菌核或子囊殼在土壤病殘體中越冬有一定的關系。

土壤微生物群落，古菌16S rRNA基因，真菌18S rRNA基因，系統進化發育分析

1 Institute of Vegetables and Flowers, Chinese Academy of Agricultural Science, Beijing 100081, China

Received: March 2, 2010; Accepted: April 30, 2010

Supported by: National Basic Research Program of China (No. 2009CB119000), Public Service Sectors (agriculture) Special Funding for Research (No. nyhyzx07-050), Key Projects in the National Science and Technology Pillar Program during the Eleventh Five-Year Plan Period (Nos.2006BAD07B03, 2006BAD08A08, 2006BAD17B08).

Corresponding author: Bingyan Xie. Tel: +86-10-82109545; E-mail: xiebingyan2003@yahoo.com.cn

國家重點基礎研究發展計劃 (973計劃) (No. 2009CB119000)，公益性行業 (農業) 科研專項經費 (No. nyhyzx07-050)，國家“十一五”國家科技支撐計劃 (Nos. 2006BAD07B03, 2006BAD08A08, 2006BAD17B08) 資助。

2 College of Biosafety Science and Technology, Hunan Agricultural University, Changsha 410128, China

Abstract:Soil archaea and fungi play important roles in the greenhouse soil ecosystem. To develop and apply rich microbial resources in greenhouse ecological environment, and to understand the interaction between microbes and plants, we constructed archaeal 16S rRNA and fungal 18S rRNA gene libraries to analyze the compositions of archaeal and fungal communityies. Total greenhouse soil DNA was directly extracted and purified by skiving-thawing-lysozyme-proteinase K-SDS hot treatment and treatment of cetyltriethylammnonium bromide (CTAB). After PCR amplification, retrieving, ligating, transforming, screening of white clones, archaeal 16S rRNA and fungal 18S rRNA gene libraries were constructed. The sequences of archaea and fungi were defined into operational taxonomic units (OTUs) when 97% similarity threshold for OTU assignment was performed by using the software DOTUR. Phylogenetic analysis showed that crenarchaeota and unidentified-archaea were the two major sub-groups and only a few of euryarchaeota existed in the archaeal clone library, total 45 OTUs. All the crenarchaeota belonged to thermoprotei;except for Basidiomycotina, the other four sub-group fungi were discovered in the fungal library, total 24 OTUs. The diversities of archaea were very abundant and a few euryarchaeota (methanebacteria) existed in the archaeal clone library, it might be directly related to the long-term high temperature, high humidity, and high content of organic matter. The limitation of oxygen was the other reason for causing this phenomenon; Ascomycotina (over 80%) was the dominant sub-groups in fungal library. It was because most of the plant fungal diseases belonged to soil-borne diseases which gone through the winter by the ways of scierotium or perithecium and became the sources of primary infection.

Keywords:soil microbial communities, archaeal 16S rRNA gene, fungal 18S rRNA gene, phylogenetic analysis

土壤微生物代表著地球環境中相當大部分的微生物量[1]。據估計，每公頃的表層土壤中包含著103~104kg的微生物生物量[2]。它的活動影響著地球環境中的營養循環、有機質的分解、土壤肥力、以及物質與能量的交換[3]。盡管具有結構和功能多樣性，長期以來，人們對其的研究極其有限，多數情況下，局限于傳統的可純培養的方法[4]。近年來，隨著分子生物學技術在微生物生態學研究中的應用，小亞基核糖體RNA (Small subunit rRNA，SSU rRNA)序列分析作為微生物分類系統的主要依據也得到了廣泛認同，隨著微生物核糖體數據庫的日益完善，該技術成為微生物分類和鑒定的一個有力工具。

在土壤微生物群體中，人們研究得最多的是細菌群體。據估計，1 g土壤樣品中，含有103~107個細菌種 (Species)[5-6]。然而古菌和真菌多樣性同樣豐富。到目前為止，僅有美國學者Fierer等[7]系統地研究過同一環境樣品中細菌、古菌、真菌和病毒多樣性，國內尚未見報道。

溫室常年種植單一作物，高溫高濕度 (年平均氣溫＞28 ℃，相對濕度＞80%)，有機質含量相當豐富。從而造成了該特定環境下微生物種類的特異性和物種豐度的穩定性，以及單一微生物物種對該環境的偏好性。因而，研究溫室土壤微生物多樣性，對了解微生物與主栽作物之間的互作關系有十分重要的意義，對開發利用溫室土壤中豐富的微生物資源以及理解土壤溫度、濕度、酸堿度等物理化學性質的變化對土壤微生物群落結構多樣性的影響提供參考依據。在之前的研究中，我們結合細菌16S rRNA基因克隆文庫和宏基因組末端測序技術，分析了溫室黃瓜近根土壤細菌多樣性，初步了解了溫室黃瓜近根土壤細菌群落結構組成。但細菌并不是唯一的土壤微生物，古菌和真菌也占有相當大的比例。

因此，為全面弄清溫室黃瓜近根土壤微生物的多樣性和群落結構組成情況，本試驗繼續對溫室黃瓜近根土壤古菌和真菌多樣性進行了研究。從中國農業科學院黃瓜溫室采集土壤樣品，直接提取和純化土壤微生物總DNA，構建古菌16S rRNA和真菌18S rRNA基因克隆文庫，分析了溫室黃瓜近根土壤古菌和真菌多樣性。對開發利用溫室這一特殊的生態環境中豐富的微生物資源以及理解微生物與植物間的互作提供了參考依據。

1 材料與方法

1.1 材料

1.1.1 溫室土壤類型及栽培管理

土壤類型為潮土，前茬作物為黃瓜 (種植年限10年以上)，品種為中農118。種植前，先將地塊深翻，并且撒上13.3 kg/hm2的蛭石，與土壤混勻，增加土壤疏松度和保濕能力。用平板鋤頭或者鍬挖出245 mm×175 mm大小的坑，播種黃瓜前，每坑補上120 g的干雞糞，與周邊土壤稍混勻。人工管理，中途不施加其他化學肥料和農藥。從苗子長出到 6葉期，少量噴霧澆水的方式，每日1次，每次10 min，6葉期后，加大噴水次數，每日分早、中、晚3次，每次10 min。每年3月份，溫室最后一季黃瓜收完之后，取近根土壤。其中一份土壤經測量濕重為200.004 g，25 ℃、100 r/min的烘箱中，烘至土壤不再粘結成塊為止，土壤干重為168.1 g。因此，估計土壤濕度在15%以上。

1.1.2 土壤樣品

土樣于2008年3月采自海淀區中國農業科學院蔬菜花卉研究所試驗農場黃瓜種植地。選擇10個種植黃瓜的地塊，采用五點采樣法，用土壤采樣器采集5個黃瓜近根 (采集直徑10 cm，采集深度0~15 cm)的土壤樣品，混合后用密封袋帶回，土樣用10 mesh(2 mm) 篩網過篩，去除小顆粒的石頭和殘存的植物根系，存放于-20 ℃冰箱中。土壤性質見表1。

1.2 DNA的提取與純化

參考Bertrand等的方法并稍作改進[8]。用100 μL TE溶解，1%瓊脂糖電泳檢測，-20 ℃保存。獲得的粗DNA用上海申能博彩公司生產的DNA柱純化試劑盒進行純化。

1.3 古菌 16S rRNA和真菌 18S rRNA基因的PCR擴增

以純化后的DNA為模板，分別用古菌16S rRNA和真菌18S rRNA基因通用引物進行PCR擴增。PCR 擴增所用反應體系：10×Buffer 5 μL，Mg2+(25 mmol/L) 4 μL，dNTPs (5 mmol/L) 2 μL，上下游引物 (10 pmol/L) 各2 μL，模板 (純化后的土壤DNA)10 ng，Taq DNA 聚合酶 (5 U) 0.7 μL，用 ddH2O 補足50 μL體系。PCR反應條件：94 ℃預變性5 min；94 ℃變性30 s，60 ℃ (或 52 ℃，擴增真菌 18S rRNA基因) 退火30 s，72 ℃延伸1 min，共34個循環；最后72 ℃延伸10 min，12 ℃保溫1 h。退火溫度和延伸時間見表2。

1.4 基因克隆文庫的構建

PCR擴增產物經1%的瓊脂糖凝膠電泳檢測，膠回收900 bp的古菌和600 bp左右的真菌PCR擴增產物，連接至pGEM-T上，大腸桿菌TOP10轉化，在含有IPTG和X-Gal以及氨芐青霉素的LB平板上37 ℃過夜培養。挑白斑，37 ℃、220 r/min振蕩培養12 h，每克隆菌液取1 μL，T7和SP6進行菌液PCR，篩選陽性克隆，構建基因克隆文庫。T7：5′-TAATACGACTCACTATAGGG-3′，SP6：5′-ACGAT TTAGGTGACACTATAG-3′。

表1 土壤理化性質Table 1 The soil samples tested in this study

表2 擴增古菌16S rRNA基因和真菌18S rRNA基因所用引物Table 2 Primers for amplifying archaeal 16S rRNA and fungal 18S rRNA genes

1.5 序列測定及序列分析

1.5.1 古菌ARDRA分型和測序

將鑒定為陽性克隆的古菌菌液 PCR產物進行rDNA擴增片段限制性內切酶分析 (Amplified ribosomal DNA restriction analysis，ARDRA)，用HinfⅠ和Hae Ⅲ對菌液PCR擴增產物37 ℃酶切2 h，根據酶切圖譜，將陽性克隆分成若干個可操作分類單元 (OTU)，統計OTU的種類和各OTU所含陽性克隆的數量[9]，將每一個OTU類型送3~5個進行測序分析，測序工作由北京諾賽生物技術公司完成。

1.5.2 真菌克隆文庫的測序

真菌克隆文庫插入片段較小，小于580 bp，不利于 ARDRA分型后的電泳檢測，酶切條帶呈輻射狀，不利于觀察，因此，我們隨機挑取 200個真菌陽性克隆進行測序，測序工作同樣由北京諾賽生物技術公司完成。

1.6 系統進化發育樹的構建和物種豐度分析

1.6.1 古菌系統進化發育樹的構建

序列相似性≥97%的認為是一個同樣的類型并且用一個單獨的序列來代表[10]，即用單個的OUT來代表, 將每個OTU代表序列在NCBI進行序列比對，并下載一致性最高的序列作為參比序列構建系統進化樹。采用Mega4.0軟件包[11]，以鄰接法構建系統進化樹。

1.6.2 真菌界系統進化發育樹的構建

將真菌克隆文庫中已測得的 176條序列，在NCBI (http://www.ncbi.nlm.nih.gov/VecScreen/Vecscreen.html) 中在線去除載體序列，并截取500 bp的片段后，在GenBank數據庫進行比對分析，由Blastn去掉同源性比對相同的結果，選取文庫中代表單一OTU的基因序列，采用Mega4.0軟件包[11]，以鄰接法構建系統進化樹。

2 結果與分析

2.1 直接提取和純化土壤微生物總DNA

直接法提取的微生物總 DNA更能代表樣品的微生物群落，能夠比較全面地了解、鑒定土壤中微生物的多樣性。但是提取DNA的同時也提取了其他有機土壤成分，比如腐質酸和褐菌酸等。本試驗采用上海申能博彩公司生產的 DNA柱純化試劑盒進行純化后可獲得高質量的土壤微生物總DNA，且片段大小＞23 kb (圖1)，適合進行后續反應。

圖1 純化后的土壤微生物總DNAFig. 1 Extraction and purification of total soil microbial DNA.M: DNA marker λHind ; 1Ⅲ?7: DNA template.

2.2 PCR擴增古菌16S rRNA和真菌18S rRNA基因

經探索，擴增兩種基因時，DNA濃度稀釋倍數在30~100倍之間均可以達到較理想的擴增效果。本試驗將7份DNA原液稀釋50倍作為模板，擴增古菌16S rRNA和真菌18S rRNA基因，每次PCR反應設置3次重復，將21次重復的PCR產物混合后，利用PCR產物回收試劑盒回收。結果分別得到900 bp和580 bp大小的片段 (圖2A、2B)。將各類回收的PCR產物分別連接在PGEM-T載體上，用于構建古菌16S rRNA和真菌18S rRNA基因克隆文庫。

2.3 基因克隆文庫的構建

2.3.1 古菌16S rRNA基因克隆文庫的構建

通過連接，轉化，藍白斑篩選，挑白斑，搖菌，挑取672個古菌克隆，用引物T7和SP6進行陽性克隆菌液PCR反應檢測，603個插入目的片段的古菌克隆，陽性率分別為90%。通過ARDRA分析，得到40個不同的可操作分類單元 (OTUs)，其中有15個OTU類型代表單克隆，最多的RFLP類型包括86個克隆。反映古菌豐度的稀缺性曲線見圖3A。

通過ARDRA獲得的CoverageC[12]值相對較高，為 97.5%。CoverageC理論上表示克隆文庫中包含的微生物種類 (OTU) 占樣品中全部微生物種類的比例，它的計算式如下：(a) C=1?n1/N，N代表文庫的庫容，n1代表在克隆文庫中僅出現一次的 OTU數量。如果 CoverageC很高，甚至達到 100%，說明庫容已經飽和。

然而，ARDRA分型有其局限性，且通過對酶切后的電泳膠圖的觀察，誤差較大。為盡可能地覆蓋所有古菌類型，挑取了310個古菌陽性克隆進行測序，測序共得到297條古菌序列。

用 DOTUR軟件分析，把序列相似性≥97%的定義為同一個OTU[10]。這樣測得的297條古菌序列歸屬為45個OTU。其中有19個OTU只有1個克隆，最多的 5個 OTU類型分別含有 39、28、24、20和18個克隆。利用DOTUR軟件生成的數據，繪制稀缺性曲線 (圖3B)，分析表明，這些克隆代表了文庫中大多數古菌的多樣性。

古菌文庫中部分克隆尚難確定其分類地位，可能代表新的屬和種。這些序列已在GenBank中登錄，并獲得核酸序列登錄號。已獲得的古菌序列登錄號有：FJ455763-FJ455766、FJ490632-FJ490675、GU394949-GU394952、GU434679-GU434688。

圖2 古菌16S rRNA基因 (A) 和真菌18S rRNA基因 (B) 的PCR擴增結果Fig. 2 PCR amplification results of total DNA with archaeal 16S rRNA universal primers (A) and PCR amplification of total DNA with fungal 18S rRNA universal primers (B). M1, M2: DNA marker DL2000; 1–7: PCR amplication of archaeal 16S rRNA; 8?12: PCR amplication of fungal 18S rRNA.

圖3 ARDRA分析繪制的古菌稀缺性曲線 (A) 和DOTUR軟件分析繪制的古菌稀缺性曲線 (B) (DNA水平上有3%的序列差異，則認為不同的OTU)Fig. 3 Rarefation curves generated for 16S rRNA genes in the archeal clone library by analysis of ARDRA (A) and rarefation curves generated for 16S rRNA genes in the archeal clone library by DOTUR software (B). OTUs were defined as groups of sequences differing by 3% at the DNA level.

2.3.2 真菌18S rRNA基因克隆文庫的構建

隨機挑選400個克隆，通過鑒定，共有386個插入目的片段的克隆，陽性率達到 96.5%。為盡可能地覆蓋環境樣品中真菌類型，挑取196個陽性克隆進行測序。測序共得到176 個真菌18S rRNA基因序列。

運用 DOTUR軟件進行分析，176 個真菌 18S rRNA基因序列歸屬25個OTUs，其中11 個OTU只有1個克隆，最多的3個OTU類型分別含有71、59和 10個克隆。利用 DOTUR軟件繪制真菌18S rRNA基因克隆文庫稀缺性曲線 (圖4)，分析表明，文庫的CoverageC為94%，這些克隆代表了大多數真菌的多樣性。

圖4 DOTUR軟件分析繪制的真菌18S基因克隆文庫稀缺性曲線Fig. 4 Rarefation curves generated for 18S rRNA genes in the clone library by DOTUR software. OTUs were defined as groups of sequences differing by 3% at the DNA level.

同樣地，少部分真菌克隆尚難確定其分類地位，可能代表新的屬和種。這些序列已在GenBank中登錄，并獲得核酸序列登錄號：GU395091-GU395093，GU394948。

2.4 系統進化發育分析

2.4.1 古菌16S rRNA基因系統進化發育分析

利用Clustalx對297條古菌16S rRNA基因序列進行全局比對，并將產生的phylip文件，輸入phylip軟件中進行排序分析，產生相應的輸出文件(Outfile)，然后將該輸出文件 (Outfile) 變成DOTUR軟件中的輸入文件(Infile)，序列相似性≥97%的定義為一個古菌序列，共得到45種不同類型的古菌。

通過對古菌16S rRNA基因克隆文庫進行系統進化發育樹分析 (圖5) 發現，該文庫主要覆蓋了嗜泉古菌 (Crenarchaeota)、未分類的古菌 (Unclassified_Archaea) 和廣域古菌 (Euryarchaeota)。其中 232個克隆屬于嗜泉古菌 (Crenarchaeota)，占克隆文庫的78.1%，56個克隆屬于未分類的古菌 (Unclassified_Archaea)，占克隆文庫的19.2%，另外，還有9個克隆屬于未純培養的廣域古菌 (Uncultured Euryarchaeota)，占克隆文庫的 2.7%，并且，這 9個克隆序列與NCBI中廣域古菌 (Euryarchaeota) 甲烷菌 (Methanococci) 的序列相似性在 91%~98 %之間，克隆HBA287 (FJ490663) 與未培養的廣域古菌克隆 (Uncultured euryarchaeote clone CFS-F12)序列相似性為96%，克隆HBA29 (GU434686) 相似性更低，僅為94%。

從系統發育樹來看 (圖 5)，嗜泉古菌(Crenarchaeota) 占了 3大部分，其中，有一部分與廣域古菌 (Euryarchaeota) 親緣關系較近，未鑒定的古菌與嗜泉古菌、廣域古菌的親緣關系較遠，且與從 NCBI上下載的極端環境 (火山口) 未鑒定的古菌序列相似性較高。

另外，所獲得的全部古菌序列在核糖體數據庫(Ribosomal Database Project，RDP) 中進行比對，除了9條序列屬于廣域古菌 (Euryarchaeota) 外，其他288條古菌序列都屬于熱變形菌綱 (Thermoprotei)，其中包括4個目，分別是脫琉球菌目 (Desulfurococales)、硫化葉菌目 (Sulfolobales)、暖球形菌目 (Caldisphaerales)和熱變形菌目 (Thermoproteales)，分別占文庫中嗜泉古菌 (Crenarchaeota) 的比例分別為48.5%、20.5%、17.8%和13.2%。

2.4.2 真菌18S rRNA基因系統進化發育分析

利用DOTUR軟件分析測序得到的176個真菌18S rRNA基因序列，把序列相似性≥97%的定義為一個真菌序列[10]。選取文庫中代表不同OUT的序列與環境樣品中的序列共同構建系統發育樹 (圖6)。這些環境樣品中的序列來自森林土壤[13]、mojave沙漠土壤[7]、草地土壤[14]、火炬松山地土壤[15]等。

圖5 古菌16S rRNA基因克隆文庫系統進化發育樹 (自舉數據集1 000次)Fig. 5 Phylogenetic tree based on archaeal 16S rRNA sequences of the clones obtained from greenhouse soil samples. The tree was constructed via the neighbor-joining method, bootstrap values above 1 000 are shown as percentage.

圖6 溫室土壤真菌18S rRNA基因克隆文庫系統進化發育樹 (自舉數據集1 000次)Fig. 6 Phylogenetic tree based on fungal 18S rRNA sequences of the clones obtained from Beijing greenhouse soil samples. The tree was constructed via the neighbor-joining method, bootstrap values above 1 000 are shown as percentage.

通過對真菌 18S rRNA基因克隆文庫進行系統進化發育 (圖6) 分析發現，文庫包含了子囊菌亞門 (Ascomycotina)、接合菌亞門 (Zygomycotina)、鞭毛菌亞門 (Mastigomycotina)、半知菌亞門(Deuteromycotina) 和未分類的不可純培養真菌(Uncultured fungi)，未發現擔子菌亞門 (Basidiomycotina)真菌。子囊菌亞門 (Ascomycotina) 是土壤真菌的優勢菌群，在測定的176條真菌18S rRNA基因序列中，有 155條序列屬于該菌群，占克隆文庫的88.1%，12個克隆屬于半知菌亞門 (Deuteromycotina)，占克隆文庫的8.5%，鞭毛菌亞門 (Mastigomycotina)、接合菌亞門 (Zygomycotina) 和不可純培養真菌(Uncultured fungi) 克隆分別為5個、4個和4個，分別占克隆文庫的2.8%、2.3%和2.3%。

從系統進化發育樹 (圖6) 上可以看出，部分未純培養的真菌克隆 (Uncultured fungus clone) 的親緣關系比較接近鞭毛菌亞門 (Mastigomycotina)，部分離子囊菌亞門 (Ascomycotina) 較近。子囊菌亞門(Ascomycotina) 中，盤菌 Pezizomycotina和毛殼菌Chaetomium是其優勢菌，分別占到子囊菌亞門(Ascomycotina) 的52.3%和33.5%。青霉菌Penicillium和鐮孢菌Fusarium是半知菌亞門 (Deuteromycotina)的優勢菌，分別占54.8%和39.7%。

3 討論

本研究結合 ARDRA分析和基因克隆文庫序列測定的方法來研究溫室黃瓜近根土壤古菌多樣性，能夠快速而準確地對克隆文庫中古菌種類進行分型，將相同酶切圖譜的克隆劃分為一種類型 (OTU)，然后，在每類 OTU中選取有代表性的克隆進行測序，這樣就避免了克隆測序的盲目性，減少了測序數量以及人力、物力和財力的使用。兩者在分析OTUs時存在一定的誤差，這可能因為構建的庫容較大，利用 ARDRA分型后，觀察結果時存在誤差，需要借助計算機才能更快更準確地完成分型工作，減少誤差。

古菌的鑒定，是我們了解古菌在其小生境中的生態功能重要的一步，對全面地分析微生物群落結構和功能多樣性有十分重要的作用。因此，古菌多樣性的研究越來越受到人們的重視。近年來，人們從極地海洋[16]、淡水湖[17]、紅樹林土壤[18]和西藏米拉山口[9]等不同的環境樣品取樣，研究古菌群落結構多樣性，同時，鑒定了大量新奇的古菌類型。此外，通過比較系統進化發育分析，揭示了古菌群體典型地系統發育特征。因此，為評價溫室古菌土壤多樣性，本研究直接以純化后的溫室土壤總DNA為模板，通過PCR擴增，構建古菌16S rRNA基因克隆文庫。

通過對文庫進行分析發現，嗜泉古菌(Crenarchaeota) 占78.1%，未分類的古菌 (Unclassified_Archaea) 占19.2%，古菌多樣性非常豐富，且有少部分未培養的廣域古菌 (甲烷菌)。這一情況的可能與溫室長期高溫高濕，有機質含量相對豐富，整個土壤長期處于高溫發酵，缺氧的狀態有關。這一結果與Leclerc等[19]的報道是相符合的。Sekiguchi等[20]也報道高溫發酵缺氧的消化池中的古菌大部分是廣域古菌 (甲烷菌)。

在嗜泉古菌 (Crenarchaeota) 中，脫琉球菌(Desulfurococales) 和硫化葉菌 (Sulfolobales) 是其優勢菌亞群。在環境樣品中，硫化還原型細菌和古菌均能將有機碳礦化成CH4和CO2，且大多數可培養的硫化還原細菌對醋酸纖維素和 H2有高度地親和力[21]。因此，本研究中脫琉球菌 (Desulfurococales)和硫化葉菌 (Sulfolobales) 與溫室土壤碳和硫的循環有關。至于以何種方式參與碳和硫的循環有待于進一步驗證。

子囊菌 (Ascomycotina) 是土壤真菌的優勢菌群，占克隆文庫的 88.1%，其次為半知菌(Deuteromycotina)，占克隆文庫的 8.5%，鞭毛菌(Mastigomycotina)、接合菌 (Zygomycotina) 和不可純培養真菌 (Uncultured fungus) 分別占克隆文庫的2.8%、2.3%和2.3%，未發現擔子菌 (Basidiomycotina)。這與近幾年來其他研究人員利用18S rRNA 基因克隆文庫的方法對土壤樣品進行分析所獲得的結果有所區別。Hunt等[14]利用18S rDNA途徑研究草地土壤真菌多樣性，分析發現該環境樣品中土壤真菌包括接合菌 (Zygomycotina)、擔子菌 (Basidiomycotina) 和子囊菌 (Ascomycotina)，而無鞭毛菌 (Mastigomycotina)和半知菌 (Deuteromycotina)。O’Brien等[15]構建了真菌ITS和18S rRNA基因克隆文庫，分析土壤樣品真菌多樣性，結果發現兩種文庫中的子囊菌(Ascomycotina) 和擔子菌 (Basidiomycotina) 所占的比例基本相當，而無其他3種類型的真菌。

該環境樣品中未發現擔子菌亞門真菌(Basidiomycotina)，可能有以下兩方面的原因：其一，擔子菌亞門真菌 (Basidiomycotina) 在該環境樣品中的豐度極低，200條測序樣品，甚至增加測序樣品到400或1 000條都難找到一個擔子菌亞門真菌(Basidiomycotina) 序列，因此，可以認為擔子菌對該基因克隆文庫中總的真菌物種豐度不會有影響；其二，前茬作物一直為黃瓜，少數真菌種群在近根選擇性富集，另一部分種群由于不適應此微環境變化而數量逐漸減少。作為該環境中弱勢種群的擔子菌亞門真菌 (Basidiomycotina)，其豐度逐漸減少，甚至檢測不出來，很可能與不適應近根微環境的變化有關。隨連作茬次增加，選擇性富集于此土壤微環境的真菌類群開始形成優勢菌群，并對其他種群產生抑制作用，導致真菌多樣性水平降低，并趨于平穩。這與胡元森等[22]的報道是一致的。

到目前為止，雖無文獻報道農事操作會影響土壤真菌群落結構，但是，土壤真菌群落結構會隨季節的變化而變化[23]。大多數植物真菌性病害是由子囊菌亞門真菌引起的，冬季，這些病原物已經成熟，以子囊殼和菌核的形式在土壤中越冬，成為來年的初侵染源。因而，子囊菌亞門的物種數相對較多，土壤真菌的群落結構組成會相應地發生變化。這也是該克隆文庫中子囊菌是優勢菌群的重要原因。

本研究應用SSU rRNA基因克隆文庫技術分別構建了溫室黃瓜近根土壤古菌16S和真菌18S rRNA基因克隆文庫，弄清了溫室黃瓜近根土壤古菌和真菌群落的組成情況。但由于病毒無高度保守的序列結構，還不能全面認識，需要進一步研究。同時，為徹底弄清高強度利用條件下溫室土壤微生物的多樣性和動態變化情況，還需考慮土壤理化性質、作物種類、施肥情況以及人事操作對微生物的影響等。

REFERENCES

[1] Whitman WB, Coleman DC, Wiebe WJ. Prokaryotes:the unseen majority. Proc Natl Acad Sci USA, 1998, 95(12):6578?6583.

[2] Brady NC, Weil RR. The nature and properties of soils.Pearson Education, Inc., Upper Saddle River, NJ, 2002.

[3] Kennedy AC. Bacterial diversity in agroecosystems. Agr Ecosyst Environ, 1999, 74(1/3): 65?76.

[4] Pace NR. A molecular view of microbial diversity and the biosphere. Science, 1997, 276(5313): 734?740.

[5] Gans J, Wolinsky M, Dunbar J. Computational improvements reveal great bacterial diversity and high metal toxicity in soil. Science, 2005, 309: 1387?1390.

[6] Tringe SG, von Mering C, Salamov AA, et al.Comparative metagenomics of microbial communities.Science, 2005, 308(5721): 554?557.

[7] Fierer N, Breitbart M, Nulton J, et al. Metagenomic and small-subunit rRNA analyses reveal the genetic diversity of bacteria, archaea, fungi, and viruses in soil. Appl Environ Microbiol, 2007, 73(21): 7059?7066.

[8] Bertrand H, Poly F, Van VT, et al. High molecular weight DNA recovery from soils prerequisite for biotechnological metagenomic library construction. J Microbial Meth,2005, 62(1): 1?11.

[9] Meng XW, Mao ZC, Chen GH, et al. Diversity of soil Archaea in Tibetan Mila Mountains. Acta Microbiol Sin,2009, 49(8): 994?1002.

孟祥偉, 茆振川, 陳國華, 等. 西藏米拉山土壤古菌16S rRNA及amoA基因多樣性分析. 微生物學報, 2009,49(8): 994?1002.

[10] Mccaig AE, Glover LA, Prosser JI. Molecular analysis of bacterial community structure and diversity in unimproved and improved upland grass pastures. Appl Environ Microbiol, 1999, 65(4): 1721?1730.

[11] Tamura K, Dudley J, Nei M, et al. MEGA4:molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol, 2007, 24(8): 1596?1599.

[12] Good IJ. The population frequencies of species and the estimation of population parameters. Biometrika, 1953,40(3/4): 237?264.

[13] Bailly J, Fraissinet-Tachet L, Verner MC, et al. Soil eukaryotic functional diversity, a metatranscriptomic approach. Interna Soc Micro Ecol, 2007, 1(7): 632?642.

[14] Hunt J, Boddy L, Randerson PF, et al. An evaluation of 18S rDNA approaches for the study of fungal diversity in grassland soils. Micro Ecol, 2004, 47(4): 385?395.

[15] O′Brien HE, Parrent JL, Jackson JA, et al. Fungal community analysis by large-scale sequencing of environmental samples. Appl Environ Microbiol, 2005,71(9): 5544?5550.

[16] Vetriani C, Reysenbach AL, Doré J. Recovery and phylogenetic analysis of archaeal rRNA sequences from continental shelf sediments. FEMS Microbiol Lett, 1998,161(1): 83?88.

[17] Jurgens G, Gl?ckner F, Amann R, et al. Identification of novel archaea in bacterioplankton of a boreal forest lake by phylogenetic analysis and fluorescent in situ hybridization. FEMS Microbiol Ecol, 2000, 34(1): 45?56.

[18] Yan B, Hong K, Yu ZN. Archaeal communities in mangrove soil characterized by 16S rRNA gene clones. J Microbiol, 2006, 44(5): 566?571.

[19] Leclerc M, Delgènes JP, Godon JJ. Diversity of the archaeal community in 44 anaerobic digesters as determined by single strand conformation polymorphism analysis and 16S rDNA sequencing. Environ Microbiol,2004, 6(8): 809?819.

[20] Sekiguchi Y, Kamagata Y, Nakamura K, et al. Fluorescence in situ hybridization using 16S rRNA-targeted oligonucleoides reveals localization of methanogens and selected uncultured bacteria in mesophilic and thermophilic sludge granules. Appl Environ Microbiol,1999, 65(3): 1280?1288.

[21] Lovely DR, Dwyer DF, Klug MJ. Kinetic analysis of competition between sulfate reducers and methanogens for hydrogen in sediments. Appl Environ Microbiol, 1982,43(6): 1373?1379.

[22] Hu YS, Liu YF, Wu K, et al. Variation of microbial community structure in relation to successive cucumber cropping soil. Chin J Soil Sci, 2006, 37(1): 126?129.

胡元森, 劉亞峰, 吳坤, 等. 黃瓜連作土壤微生物區系變化研究. 土壤通報, 2006, 37(1): 126?129.

[23] Daniell TJ, Husband R, Fitter AH, et al. Molecular diversity of arbuscular mycorrhizal fungi colonising arable crops. FEMS Microbiol Ecol, 2001, 36(2/3): 203?209.

Diversity analysis of archaeal and fungal communities in adjacent cucumber root soil samples in greenhouse by small-subunit rRNA gene cloning

Zhixiang Zhao1,2, Xiaofei Lu1, Guohua Chen1, Zhenchuan Mao1, Yuhong Yang1, Erming Liu2,and Bingyan Xie1

Received: April 19, 2010; Accepted: July 27, 2010

Supported by: National Basic Research Program of China (973 Program) (No. 2010CB833603), National Natural Science Foundation of China (No.30872404).

Corresponding author: Jiye Cai. Tel/Fax: +86-20-85223569; E-mail: tjycai@jnu.edu.cn

國家重點基礎研究發展計劃 (973計劃) (No. 2010CB833603)，國家自然科學基金 (No. 30872404) 資助。