基于高通量測序的磷石膏上芒草內部細菌多樣性分析

周煬　李劍峰　郭欣楊　張淑卿　倪明秀　張煜冕　郭金梅

周煬，李劍峰，郭欣楊，等.基于高通量測序的磷石膏上芒草內部細菌多樣性分析[J].南方農業，2023，17（11）：25-29.

摘 要 為了探究磷石膏堆渣上野生芒草（Miscanthus Sinensis）植株及種子內部微生物組成情況，給今后磷石膏堆渣無害化處理、利用提供一定的依據，以貴州省開陽縣露天磷石膏堆場上生長芒草植株及種子為試驗材料，通過高通量測序技術分析芒草植株及種子內部微生物多樣性。結果共獲得123 363條有效序列，通過聚類獲得274個OTUs，屬于17個門42個綱148個屬；在門水平上，芒草植株、種子內部的優勢細菌門為藍菌門（Cyanobacteria_Chlfaxianoroplast），分別占59.03%、79.29%；在屬水平上，優勢菌屬為鞘氨醇單胞菌屬細菌（Sphingomonas）、鏈霉菌屬細菌（Streptomyces），分別占7%、2%。

關鍵詞 磷石膏；芒草；高通量測序；細菌多樣性；環境修復

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

磷石膏是濕法磷酸生產過程中的副產物，主要成分為二水硫酸鈣，除此之外還有部分重金屬元素、可溶性磷酸鹽及氟化物等雜質[1]。經計算每生產1 t濕法磷酸，將排放出約5 t的磷石膏廢渣[2]。目前，國內對于磷石膏利用率不足40%，大部分是通過堆存方式處理，這種處理方式不僅占據了大量的土地資源，而且其潛在的污染物泄露和垮堆風險也會對當地居民的身體和生態環境造成威脅[3]。磷石膏的低成本處理與有效利用始終是我國磷化工業的一大難題。

芒草（Miscanthus Sinensis）為禾本科（Gramineae）芒屬植物[4]，通常分布在熱帶、亞熱帶及溫帶亞洲地區。芒草在我國分布較廣，其對于環境的適應能力非常強，在各種惡劣環境下都能夠生存[5]。目前已有將芒草應用于重金屬土壤修復的相關研究報道，并發現芒草能夠耐受一種甚至多種高濃度重金屬元素的脅迫[6]。貴州作為磷石膏生產排放的集中區域，且磷石膏的產出量在逐年增長，對于磷石膏的有效處理已成為當今貴州乃至全國的熱點問題。目前，有關利用芒草植株修復磷石膏堆渣污染的相關研究報道還較為少見，因此具有較大的研究意義。

1? 材料與方法

1.1? 材料

2021年7月，在貴州省貴陽市開陽縣龍景灣貴州磷化集團開陽露天堆場（東經106.852 3°，北緯27.151 3°）采集磷石膏上生長的芒草植株和種子樣品，采集時芒草呈深綠色，群落單一，株高在0.4～1.2 m。

1.2? 方法

1.2.1? 樣品處理

將植株根系挖出，用無菌水清洗表面后放入陰涼通風處自然風干，待植株風干后與種子一同裝入無菌自封袋內放入4 ℃冰箱暫存。

1.2.2? DNA提取及測序

根據試劑盒要求提取樣本的DNA，利用1%的瓊脂糖凝膠電泳測定DNA純度。采用16SrDNA對于V3-V4引物進行擴增，引物序列為（上游引物338F：ACTCCTACGGGAGGCAGCAG，下游引物806R：GGACTACHVGGGTWTCTAAT）PCR反應參數：a. 1×（3 min at 95 ℃）b.循環數×30 s at 95 ℃；30 s at 55 ℃；45 s at 72 ℃）c.10 min at 72 ℃，10 ℃ until halted by user。根據各樣品的測序量，將其配置成對應的配比，送至生工（上海）生物工程有限公司，利用Illumina公司的Miseq PE300，對其進行高通量測序。

2? 結果與分析

2.1? 細菌多樣性分析

2.1.1? 測序數據處理結果

如表1所示，樣品經過高通量測序共獲得124 097條原始序列，經對原始數據初步處理，獲得了123 363條有效的序列。供試樣本有效序列條數在51 383～71 980條，其中phps（芒草種子）有效序列條數大于php（芒草植株）。通過聚類一共獲得274個OTUS，通常將相似性值≥97%的序列定義為同一個OTUS。為了比較直觀地分析和統計，對于測序樣本的共有OTUS和特有OTUS進行分析比較，結果顯示，php與phps的共有OTUS有138個，php特有OTUS為66個，phps特有OTUS為70個，即php與phps的微生物群落組成較為相似。

2.1.2? α多樣性

樣品稀釋曲線傾斜度在不斷降低逐漸趨于平緩（見圖1），同時樣品中測序覆蓋度均為1（見表2），說明測序數據量比較合理，基本上能夠反映出樣品細菌群落的組成情況[7]。通過計算Ace指數、Chao指數、Shannon指數等三個指標結果表明，php細菌多樣性大于phps。

Rank Abundance曲線圖能夠直觀明確地展示各樣本中OTUS的均勻度及豐富度，如圖2所示，當OTU數量小于100時，php與phps的物種相對豐度較高，但曲線下降速度比較快，說明php與phps的分布均勻程度較低。當OTU數量超過100的情況下，物種的相對豐度較低，但是曲線比較平坦，表明物種分布較為均衡。圖中php與phps樣本間的曲線寬度、平滑程度都比較一致，表明其物種豐度及均勻度都無顯著差異。

2.1.3? 樣本主成分

通過主成分（PCA）分析可知，主成分1（PC1）和主成分2（PC2）貢獻率分別為58.23%和40.5%，合計為98.73%。如圖3所示，兩樣本點的距離十分接近，且位于第一象限有明顯的聚集，表明php與phps的細菌相似性較高。

2.2? 細菌群落結構差異分析

2.2.1? 門水平下陳化磷石膏上生長芒草植株及種子細菌群落結構

在門水平下，一共鑒定出17個門類細菌，如表3所示，豐度低于1%的細菌門合并為Others，其中，芒草植株（php）及種子（phps）共有的優勢細菌門為藍菌門（Cyanobacteria_Chloroplast）。如圖4所示，藍菌門細菌在芒草植株、種子內部的相對豐度分別為59%、79.3%，即藍菌門的相對豐度phps大于php。原因可能是：1）種子相對于植株而言，營養物質積累較多，較多的營養物質有利于藍菌門的生長，所以種子內部的藍菌門豐度較大；2）在細菌功能方面，藍菌門的菌類具有較強的抗紫外線輻射能力，在較強紫外線的輻射下能夠形成一種代謝響應機制，使其適應紫外輻射的脅迫[8]。

2.2.2? 屬水平下陳化磷石膏上生長芒草植株及種子細菌相對豐度

在屬水平下，共檢測出148個屬的細菌，在芒草植株中鏈型植物（Streptophyta 59%）占比最大，其次是未知菌屬細菌（unclassified_Bacteria 39.5%）；鞘氨醇單胞菌屬（Sphingomonas 0.07%）為植株菌群內部的重要細菌群落，在芒草種子（phps）中占比最大的同樣為鏈型植物（Streptophyta 79.2%），其次為未知菌屬（unclassified_Bacteria 19.8%）。如表4、圖5所示，在植株到種子過程中，其內部鏈型植物、未知菌豐度逐漸增大，但鞘氨醇單胞菌屬、鏈霉菌屬等細菌豐度相對減少。鏈型植物占比較高的原因可能為芒草植株及種子內部葉綠體信息沒有被排除，導致測序結果中鏈型植物數據占了絕大部分。芒草植株之所以能夠在磷石膏上生長，原因可能是鞘氨醇單胞菌屬、鏈霉菌屬細菌對植物的生長也起著一定的促進作用。鞘氨醇單胞菌屬細菌常被研究者用于修復重金屬污染土壤[9]。

2.3? 16s功能預測分析

如圖6所示，在php與phps中，其功能性強的前五類均為：信號轉導組氨酸激酶、糖基轉移酶、由CheY樣受體域和翼螺旋DNA結合結構域組成的響應調節劑、藍藻中保守的未表征蛋白質、未表征的低復雜度蛋白質，且種子內部功能豐度值高于植株。其中植株與種子中功能性最強的均為信號轉導組氨酸激酶。組氨酸激酶可以作為細菌細胞中的重要信號感知分子，細菌可以利用組氨酸激酶來感受胞外環境的變化，作為一種自身具有激酶、磷酸轉移酶和磷酸酶活性的多功能酶，組氨酸激酶能夠作為信號分子的細胞受體，且具有跨細胞膜信號轉導的功能。這種轉導信號在細菌、酵母菌和高等植物中都有發現[10-11]。

3? 討論與結論

本研究以磷石膏上芒草植株和種子為研究對象，利用高通量測序技術對微生物多樣性進行細致研究，結果顯示：在門水平上，其共同的優勢菌門是藍菌門。在屬水平上，鞘氨醇單胞菌屬細菌（Sphingomonas）、鏈霉菌屬細菌（Streptomyces）為優勢菌屬，且php中兩種優勢細菌相對豐度高于phps。這兩種細菌在植物生長及抗逆性方面都起著較為重要的作用。黃媛林等對一株紅樹林鏈霉菌所產抑菌活性化合物的分離及其生物合成基因簇的研究表明，鏈霉菌的次生代謝產物具有多種抗性的功效，且多被用作生物防治劑[12]。何曉紅等對鞘氨醇單胞菌的各項功能的綜述研究表明，鞘氨醇單胞菌能夠耐受惡劣環境及貧瘠營養環境，且能夠促進植物生長與抵抗多種植物病原菌[13]。

目前，鏈霉菌屬與鞘氨醇單胞菌屬已廣泛用于重金屬污染的土壤植物與植物促生菌的聯合修復研究[14]，微生物環境修復仍然是目前研究的熱點問題之一。然而現階段對于磷石膏的處理利用及其內部元素對于環境污染的系統研究還比較缺乏，這些都將成為未來磷石膏堆渣污染研究的重要課題。芒草植株能夠在磷石膏堆渣上面生存的原因可能是由于鏈霉菌與鞘氨醇單胞菌的存在，這兩株細菌具備開發成為磷石膏堆渣治理與利用的植物促生菌的潛在價值，值得進一步開發與研究。

參考文獻：

[1]? 李東旭，焦嘉偉，廖大龍.我國磷石膏在新型建筑材料行業的發展及思考[J].新型建筑材料，2021，48（11）：1-4，11.

[2]? 龍謹，楊鳳龍，吳欽雪，等.磷石膏對土壤的性能影響研究綜述[J].廣州化工，2021，49（19）：35-37.

[3]? 彭曉曦，吳怡，黃祥，等.磷石膏快速干法穩定化的礦物特征研究[J].無機鹽工業，2022，54（4）：24-28.

[4]? 趙銳名，成小英.鉛脅迫下產鐵載體細菌對芒草幼苗的促生特性研究[J].安徽農業科學，2019，47（16）：93-97.

[5]? 陳許兵，朱廣龍，任志強，等.能源植物芒草耐鹽性研究進展[J].安徽農業科學，2020，48（13）：12-15，20.

[6]? 吳道銘，陳曉陽，曾曙才.芒屬植物重金屬耐性及其在礦山廢棄地植被恢復中的應用潛力[J].應用生態學報，2017，28（4）：1397-1406.

[7]? 張守梅，郭玉秋，張娜娜，等.基于高通量測序分析玉米浸泡副產物中細菌群落結構和多樣性[J].山東農業科學，2021，53（12）：143-148.

[8]? 劉榮，于海峰，逄圣慧，等.藍細菌中紫外吸收物質研究進展[J].食品工業科技，2012，33（11）：388-390.

[9]? BAO C. Improvement of cadmium uptake and accumulation in Sedum alfredii by endophytic bacteria Sphingomonas SaMR12： Effects on plant growth and root exudates[J]. Chemosphere， 2014， 117： 367-373.

[10]? 陳一.擬南芥跨膜組氨酸激酶ATHK1研究進展[J].安徽農業科學，2014，42（7）：1923-1924，2036.

[11]? 高雅，金華燕，付鳴佳，等.蛹蟲草組氨酸激酶基因（CmHK）的序列分析和受藍光影響的表達[J].應用與環境生物學報，2014，20（1）：50-55.

[12]? 黃媛林，潘信利，陸璐，等.一株紅樹林鏈霉菌所產抑菌活性化合物的分離及其生物合成基因簇的研究[J].中國抗生素雜志，2021，46（12）：1090-1101.

[13]? 胡杰，何曉紅，李大平，等.鞘氨醇單胞菌研究進展[J].應用與環境生物學報，2007（3）：431-437.

[14]? 陳兆進，林立安，李英軍，等.鎘脅迫對芒草根際細菌群落結構、共發生網絡和功能的影響[J].環境科學，2021，42（8）：3997-4004.

（責任編輯：丁志祥）