程立君，王 磊，吳銀梅 ，高興國，師 睿 ，胡志芳

（1.昭通學院農學與生命科學學院，云南 昭通 657000；2.昭通高原特色農產品檢測平臺，云南 昭通 657000；3.昭通蘋果產業(yè)研究所，云南 昭通 657000）

【研究意義】珙桐（Davidia involucrateBaill）又稱為鴿子樹，屬珙桐科（Nyssaceae）珙桐屬（Davidia）雙子葉植物，是我國特有的第三紀古熱帶植物的孑遺樹種，被列為國家一級保護植物［1-2］，具有很高的觀賞價值和科研價值。但珙桐自然分布不連續(xù)，資源量很少，外界自然環(huán)境條件是影響其分布和繁育的主要因素，其中干旱是主要生態(tài)限制因子之一。因此可以嘗試篩選珙桐根際微生物優(yōu)勢菌群，利用土壤微生物的活動能力，提升珙桐的逆境條件下的生存能力?！厩叭搜芯窟M展】自20世紀80年代以來，前人對珙桐的生物學、生態(tài)學、人工繁殖技術及引種栽培、生理生化以及分子生物學和遺傳多樣性等多個方面開展了大量的研究，并取得了一些進展，但在珙桐引種繁殖、植物化學成分等研究方向上還需更深層次的發(fā)展，尤其是在異地引種和開花結實方面研究不夠成熟，限制了珙桐作為園藝觀賞植物進行廣泛栽培。有研究者對其繁育和引種開展了許多工作，但效果均不理想［3-5］，原因是珙桐對土壤干旱脅迫極為敏感［6］，耐旱能力較差［7］。土壤微生物是土壤生態(tài)系統(tǒng)中最活躍的成分［8］，其參與物質循環(huán)代謝過程，在改善土壤理化性質、調節(jié)土壤養(yǎng)分的貯存和釋放等方面起關鍵作用，并影響植物生長［9-10］，根際微生物多樣性能緩解干旱和鹽堿脅迫對植物生理的影響［11］，尤其是菌根菌能與植物根微系形成菌根共生體，菌根的物理生理特性能提高植物在逆境條件下的生存能力，如氣候干旱、土壤貧瘠等障礙因素［12］，如胡楊根際微生物對胡楊的逆境抗性有重要作用，王珊珊［13］從胡楊根際土壤中篩選得到10株高效促生抗逆菌；有研究表明，叢枝菌根真菌（AMF）可提高青楊雌雄株抗旱性［14］。【本研究切入點】本研究嘗試采用Illumina HiSeq高通量測序技術，對永善三江口、彝良縣小草壩和鎮(zhèn)雄塘房3個地區(qū)的光葉珙桐根際土壤中細菌和真菌的多樣性進行研究，分析光葉珙桐自然條件下根際土壤微生物的共同特征與差異?！緮M解決的關鍵問題】通過對比不同地區(qū)光葉珙桐根際土壤微生物多樣性和主要菌群的差異，為今后篩選提高其抗旱能力的根際微生物和探索引種繁育新途徑提供基礎數(shù)據(jù)和理論參考。

1 材料與方法

1.1 樣品采集

2017年8月，從昭通地區(qū)3個光葉珙桐自然分布區(qū)（永善三江口、彝良小草壩、鎮(zhèn)雄塘房）進行土樣采集（避免雨季），具體操作是在每個地點用土壤取樣器采集光葉珙桐根基土壤各5份，深度為0～40 cm，將在多個地點采集的土壤分別進行充分混合，并除去較大的石塊、雜草、根系等雜物后，對獲得的3個土樣分別采用四分法留取土樣，于-80℃冷凍保存，用于土壤微生物多樣性分析。

1.2 試驗方法

1.2.1 土壤微生物總DNA提取 準確稱取0.30 g土壤樣品，采用OMEGA試劑盒法提取土壤樣品中微生物總DNA，用1%瓊脂糖凝膠電泳檢測DNA完整性，用Mini Dorp測定DNA純度和濃度。

1.2.2 PCR擴增及測序 樣品總DNA提取后，采用Illumina HiSeq 2500測序平臺的通用引物對細菌16S rRNA V3+V4區(qū)域和真菌ITS1區(qū)域分別進行PCR擴增。反應程序：95℃預變性5 min；95℃變性35 s、57℃退火30 s、72℃延伸40 s，30個循環(huán)；72℃延伸5 min。對擴增產物進行純化、定量和均一化，形成測序文庫，由北京百邁客生物科技有限公司通過Illumina HiSeq 2500進行測序。PCR擴增通用引物包括V3+V4區(qū)域（428 bp）引物（上游引物5'-CCTACGGGNGGCWGCAG-3'，下游引 物 5'-GACTACHVGGGTATCTAATCC-3'） 和真菌ITS1區(qū)域（187～298 bp）引物（上游引物5'-CTTGGTCATTTAGAGGAAGTAA-3'，下游引物5'-GCTGCGTTCTTCATCGATGC-3'）。

1.3 數(shù)據(jù)分析與處理

1.3.1 生物信息學分析 原始數(shù)據(jù)存在一定的干擾數(shù)據(jù)，先將測得的原始序列截去Barcode序列和引物序列，使用FLASH v1.2.7軟件，通過overlap對每個樣品的reads進行拼接，得到的拼接序列即原始Tags數(shù)據(jù)（Raw Tags）；使用Trimmomatic v0.33軟件，對拼接得到的Raw Tags進行過濾，得到高質量的Tags數(shù)據(jù)（Clean Tags）；使用UCHIME v4.2軟件鑒定并去除嵌合體序列，得到最終有效數(shù)據(jù)（Effective Tags）。使用QIIME軟件中的UCLUST對Tags在97%的相似度水平下進行聚類，獲得OTU，統(tǒng)計各個土壤樣品每個OTU中的豐度信息。采用隨機抽樣的方法抽取數(shù)據(jù)，以抽到的序列數(shù)與它們所能代表OTU的數(shù)目構建曲線，即稀釋性曲線，利用Mothur軟件和R語言工具繪制稀釋曲線圖。

1.3.2 物種組成分析 基于Silva（細菌）和UNITE（真菌）分類學數(shù)據(jù)庫對OTU進行分類學注釋，確定序列的分類信息，獲得門、綱、目、科和屬各水平下的分類單元。利用QIIME軟件生成不同分類水平上的物種豐度表，再利用R語言工具繪制成樣品各分類學水平下的群落結構圖。

1.3.3 多樣性分析 Alpha多樣性反映的是單個樣品物種豐度及物種多樣性，有Chao1、Ace、Shannon、Simpson等多種衡量指標。基于OTU的結果，利用Mothur軟件計算豐富度指數(shù)（Chao、Ace）和Alpha多樣性指數(shù)（Shannon、Simpson）來進行生物多樣性分析。使用QIIME軟件進行Beta多樣性分析，比較不同樣品在物種多樣性方面存在的相似程度?；诟鳂悠沸蛄虚g的進化信息的差異來計算樣品間的差異，利用距離算法（binary、bray、weighted、unweighted）得到樣品間的距離矩陣，通過R語言工具繪制樣品熱圖（Heatmap）和采用非加權配對平均法（UPGMA）對樣品進行層次聚類繪制系統(tǒng)發(fā)生樹。

2 結果與分析

2.1 光葉珙桐根際土壤細菌和真菌OTU水平分析

經(jīng)數(shù)據(jù)前處理，3個樣品（永善三江口、彝良小草壩和鎮(zhèn)雄塘房）總DNA的16S rDNA V3+V4區(qū)經(jīng)測序分析得到67 712～68 829條有效序列數(shù)據(jù)（表1），序列長度集中在390～450 bp之間，平均長度分別為419、415、421 bp，從序列長度來看，與16S rDNAV3+V4區(qū)序列長度大致吻合，進行聚類后分為2 293個OTUs（序列相似性＞97%）。3個樣品總DNA的ITS區(qū)經(jīng)測序分析得到69 510～69 717條有效序列數(shù)據(jù)（表1），序列長度集中在160～340 bp之間，平均長度分別為230、218、244 bp，進行聚類后分為880個OTUs。從3個土樣中細菌和真菌的rarefaction曲線（圖1）可以看出，隨著測序數(shù)量的增加，它們所能代表構建的OTU數(shù)目增加，曲線的斜率逐漸上升，后趨于平緩，同時樣品所得到的序列的覆蓋度均高于90%，說明本研究測序數(shù)據(jù)量合理，可反映樣品中微生物群落的種類和結構。

2.2 光葉珙桐根際土壤微生物物種組成分析

本研究檢測到的細菌和真菌種類繁多，其中很多物種的含量非常少，為使視圖效果最佳，柱狀圖只顯示豐度水平排名前10的真菌和細菌優(yōu)勢菌，并將其他物種合并為Others在圖中顯示，Unknown代表未得到分類學注釋的物種，以相對豐度為縱坐標繪制柱形圖。圖2A和圖2B分別展示了樣品中相對豐度排名前10個綱的細菌和真菌，表明光葉珙桐根際土壤中細菌主要分布于以下10個綱：α-變形菌綱Alphaproteobacteria、Subgroup_6、δ-變 形 菌 綱 Deltaproteobacteria、γ-變形菌綱Gammaproteobacteria、酸桿菌綱Acidobacteria、β-變形菌綱 Betaproteobacteria、索利氏菌綱Solibacteres、Subgroup_2、放線菌綱Actinobacteria、芽單胞菌綱Gemmatimonadetes（圖2A），其中α-變形菌綱Alphaproteobacteria最豐富，在3個土樣中相對豐度分別為19.37%、18.48%、17.02%。永善三江口、鎮(zhèn)雄塘房土樣中主要細菌種類的相對豐度比較相近，與彝良小草壩土樣差異較大，尤其是在Subgroup_2和算桿菌綱Acidobacteria的相對豐度分別高達14.93%和10.05%，而在永善三江口土樣（0.36%、2.50%）和鎮(zhèn)雄塘房土樣（0.04%、0.14%）中含量較低；真菌主要分布于以下10個綱：糞殼菌綱Sordariomycetes、傘菌綱Agaricomycetes、古菌綱 Archaeorhizomycetes、Mortierellomycotina_cls_Incertae_sedis、盤菌綱Pezizomycetes、散囊菌綱Eurotiomycetes、座囊菌綱Dothideomycetes、球囊菌綱Glomeromycetes、錘舌菌綱Leotiomycetes、銀耳綱Tremellomycetes（圖2B），不同地區(qū)的土樣中真菌菌群分布有明顯差異，在永善三江口土樣中相對豐度最高的3個種類是盤菌綱Pezizomycetes（17.01%）、糞殼菌綱Sordariomycetes（15.04%）、Mortierellomycotina_cls_Incertae_sedis（7.64%），在彝良小草壩土樣中相對豐度最高的3個種類是古菌綱Archaeorhizomycetes（29.79%）、傘菌綱Agaricomycetes（20.90%）、糞 殼 菌 綱 Sordariomycetes（4.87%）， 在 鎮(zhèn)雄塘房土樣中相對豐度最高的3個種類是糞殼 菌 綱 Sordariomycetes（20.77%）、 傘 菌 綱Agaricomycetes（10.53%）、Mortierellomycotina_cls_Incertae_sedis（8.35%）， 其 中 盤 菌 綱Pezizomycetes在永善三江口土樣中相對豐度為17.01%，而在彝良小草壩和鎮(zhèn)雄塘房土樣中含量很低或未檢測出，傘菌綱Archaeorhizomycetes在彝良小草壩土樣中相對豐度為29.79%，而在永善三江口、鎮(zhèn)雄塘房土樣中含量很低甚至未檢測出。相對豐度大于1%的種類在3個樣品中都有分布的只有傘菌綱Agaricomycetes、Mortierellomycotina_cls_Incertae_sedis和糞殼菌綱Sordariomycetes，其中傘菌綱Sordariomycetes在3個樣品中都具有較高的相對豐度，分別為15.04%、4.87%、20.77%。

表1 光葉珙桐根際土壤樣品測序的tags信息及OTU統(tǒng)計分析Table 1 Tags information and OTU statistical analysis of rhizosphere soil sample sequencing of D. involucrate

圖1 光葉珙桐根際土壤樣品細菌和真菌測序的稀釋曲線Fig. 1 Dilution curves of bacterial and fungi sequencing of rhizosphere soil samples of D. involucrate

圖2 光葉珙桐根際土壤樣品細菌群落和真菌群落分布Fig. 2 Distribution of bacterial community and fungi community in rhizosphere soil samples of D. involucrate

2.3 光葉珙桐根際土壤微生物多樣性分析

2.3.1 Alpha多樣性分析 從Alpha多樣性指數(shù)結果（表2）可以看出，3個樣品都具有較高的細菌和真菌群落多樣性。3個土樣中細菌群落多樣性Chao1指數(shù)平均值和Shannon指數(shù)平均值表現(xiàn)為永善三江口＞鎮(zhèn)雄塘房＞彝良小草壩，可見永善三江口光葉珙桐根際土壤細菌多樣性最豐富；真菌群落多樣性Chao1指數(shù)值表現(xiàn)為鎮(zhèn)雄塘房＞彝良小草壩＞永善三江口，Shannon指數(shù)值表現(xiàn)為鎮(zhèn)雄塘房＞永善三江口＞彝良小草壩?？梢?，光葉珙桐根際土壤真菌豐富度方面鎮(zhèn)雄塘房＞彝良小草壩＞永善三江口，多樣性方面鎮(zhèn)雄塘房＞永善三江口＞彝良小草壩。

表2 光葉珙桐根際土壤樣品根際土壤群落多樣性指數(shù)Table 2 Community diversity index of rhizosphere soil samples of D. involucrate

2.3.2 Beta多樣性分析 圖3為基于距離算法（Unweighted）得到3個土樣間的距離矩陣，通過R語言工具繪制樣品熱圖Heatmap，根據(jù)顏色變化可以看出樣品間的差異性，結果表明永善三江口和鎮(zhèn)雄塘房地區(qū)光葉珙桐根際土壤中細菌菌群（圖3A）和真菌菌群（圖3B）的相似性均較高。

圖3 光葉珙桐根際土壤樣品加權UniFrac的heatmap圖Fig. 3 Heatmap diagram of Weighted Unirac of rhizosphere soil samples of D.involucrate

根據(jù)每個樣品OTU的組成情況，基于Beta多樣性分析得到的距離矩陣，通過R語言工具采用非加權配對平均法（UPGMA）對樣品進行層次聚類，可以判斷樣品間物種組成的相似性，圖4是3個土樣中細菌和真菌OTU豐度的聚類分析結果，樣品越靠近，枝長越短，說明兩個樣品的物種組成越相似，結果顯示永善三江口和鎮(zhèn)雄塘房地區(qū)OTU水平相似，與彝良小草壩地區(qū)根際細菌和真菌菌群的相似性較低。

圖4 基于OTU的光葉珙桐根際土壤樣品樣本樹狀聚類結果Fig. 4 Clustering diagram of rhizosphere soil samples of D. involucrate based on the OUT

3 討論

珙桐根際土壤微生物多樣性研究報道較少，楊敬天等［15］對珙桐土壤微生物數(shù)量及其與土壤因子的關系做過報道，但未能反映微生物種類信息、豐度和群落結構等情況。本研究采用Illumina HiSeq高通量測序技術相較于傳統(tǒng)分子生物學方法，該技術覆蓋了整個土壤微生物群落的信息，能準確反映微生物群落的種類、豐度和分布。秦越等［16］研究認為，連作會導致馬鈴薯根際土壤從“細菌型”向“真菌型”土壤轉化。通過對永善三江口、彝良小草壩和鎮(zhèn)雄塘房3個地區(qū)光葉珙桐根際土壤微生物多樣性分析結果表明，永善三江口的光葉珙桐根際土壤細菌多樣性最豐富，鎮(zhèn)雄塘房的光葉珙桐根際土壤真菌多樣性最豐富。羅柏青等［17］對滇東北珙桐群落調查的結果顯示，永善三江口和彝良小草壩珙桐徑級分布呈現(xiàn)出穩(wěn)定增長的金字塔型，種群處于穩(wěn)定階段。夏嘉禹等［18］對鎮(zhèn)雄縣珙桐群落結構調查結果顯示呈現(xiàn)倒金字塔形更新較差，處于衰退階段。該珙桐群落結構的調查結果與本研究3個地區(qū)光葉珙桐根際土壤微生物多樣性結果相吻合，珙桐種群穩(wěn)定的地區(qū)（永善三江口和彝良小草壩）根際土壤細菌多樣性較豐富，種群處于衰退階段的地區(qū)（鎮(zhèn)雄）光葉珙桐根際土壤真菌多樣性較豐富。光葉珙桐根際土壤細菌和真菌多樣性的這種變化與光葉珙桐自然分布區(qū)種群結構變化是否也存在聯(lián)系還有待進一步研究。

4 結論

永善三江口、彝良縣小草壩和鎮(zhèn)雄塘房3個地區(qū)光葉珙桐根際土壤中細菌和真菌物種組成豐富，但不同地區(qū)光葉珙桐根際土壤中微生物種群結構有一定差異。在綱分類水平上，細菌中α-變形菌綱Alphaproteobacteria是最優(yōu)勢綱，永善三江口、彝良縣小草壩和鎮(zhèn)雄塘房3個地區(qū)的相對豐度分別為19.37%、18.48%、17.02%，真菌中只有糞殼菌綱Sordariomycetes在3個樣品中都具有較高的相對豐度，分別為15.04%、4.87%、20.77%，其他豐度較高的種類都存在明顯差異。通過微生物種群結構分析結果顯示，永善三江口與鎮(zhèn)雄塘房的光葉珙桐根際土壤微生物種群結構的相似性較高，其中永善三江口光葉珙桐根際土壤細菌多樣性最豐富，鎮(zhèn)雄塘房光葉珙桐根際土壤真菌多樣性最豐富。