耐鹽堿溶磷菌Y2R2的分離鑒定及溶磷特性

楊海霞 劉希旻 潘奕臣 趙香林 黃海東

（天津農(nóng)學(xué)院 農(nóng)學(xué)與資源環(huán)境學(xué)院，天津 300384）

土壤鹽漬化是影響農(nóng)業(yè)產(chǎn)量的主要因素之一，全球鹽漬土壤面積已達(dá)9.5×108hm2，其中由于灌溉、施肥不當(dāng)及城市融雪等人類(lèi)活動(dòng)造成的次生鹽漬化日益嚴(yán)重［1-3］。鹽漬土及堿性環(huán)境下，大量磷元素轉(zhuǎn)化為難溶性的磷酸鹽而被土壤固定，從而出現(xiàn)土壤鹽堿化和“遺傳學(xué)缺磷”并存的現(xiàn)象，新施入磷肥的當(dāng)季作物利用率低，造成土壤的潛在磷庫(kù)雖很大，但有效磷含量卻很低，植物很難直接吸收利用［4-5］。前人研究發(fā)現(xiàn)，自然界中存在大量具有溶磷能力的微生物，這類(lèi)溶磷微生物能分泌有機(jī)酸、磷酸酶或以其他方式將被土壤固定的難溶性磷轉(zhuǎn)化為植物可利用的有效磷［6］。將溶磷菌制成菌劑在農(nóng)田中施用，能在不增加土壤磷庫(kù)的情況下，提供更多的磷元素供作物生長(zhǎng)，溶磷菌還能改善土壤的鹽堿條件，增加土壤微生物群落多樣性，通過(guò)在作物根際的繁殖而抑制其他病原菌的生長(zhǎng)，從而減輕作物病害［7］。目前篩選到溶磷能力較強(qiáng)的細(xì)菌包括芽孢桿菌屬（Bacillus）、節(jié)桿菌屬（Arthrobacter）、曲霉屬（Penicillium）和假單胞菌屬（Pseudomonas）的一些菌株［8］，這些溶磷菌大多從耕作土壤或植物根際中篩選得到［9-10］，對(duì)鹽堿土壤的適應(yīng)能力較差。鹽堿土壤中有一些優(yōu)勢(shì)的種群，如鹽單胞菌屬（Halomonas）。這種革蘭氏陰性無(wú)芽孢的細(xì)菌大多可以耐受12%以上的鹽度［11］，但目前鹽堿土壤中鹽單胞菌屬溶磷特性的研究還未見(jiàn)報(bào)道。

天津?yàn)I海新區(qū)屬于海積平原地貌，土壤長(zhǎng)期受到海水和風(fēng)暴潮的浸蝕，地下水位高，排水不暢，形成高鹽堿度的濱海鹽漬土［12-13］，是耐鹽堿微生物篩選的理想采樣地。本實(shí)驗(yàn)從天津?yàn)I海新區(qū)鹽堿土中分離到一株耐鹽堿溶磷菌Y2R2，對(duì)該菌株進(jìn)行了分類(lèi)鑒定，并研究了其溶磷特性，旨為今后制作耐鹽堿溶磷菌劑提供菌種資源及濱海鹽堿地改良和解決鹽堿土壤固磷強(qiáng)度高的問(wèn)題提供支持。

1 材料與方法

1.1 材料

1.1.1 試劑 土樣采自于天津?yàn)I海新區(qū)鹽堿地土壤，北緯39°4'46''，東經(jīng)117°41'37''。細(xì)菌基因組DNA提取試劑盒，Taq酶，質(zhì)粒提取試劑盒，普通DNA產(chǎn)物純化試劑盒，pGM-T克隆試劑盒均等購(gòu)自北京天根公司，鉬酸銨、丙酮酸鈉、抗壞血酸等均為分析純。

1.1.2 培養(yǎng)基 R2A培養(yǎng)基（g/L）：酵母膏0.5，蛋白胨0.5，酸水解酪素0.5，磷酸氫二鉀0.3，葡萄糖0.5，可溶性淀粉0.5，丙酮酸鈉0.3，硫酸鎂0.05，氯化鈉20，瓊脂15，pH 9.0±0.2。無(wú)機(jī)磷培養(yǎng)基（g/L）：葡萄糖10，硫酸銨0.5，酵母浸粉0.5，氯化鉀0.3，硫酸鎂0.3，硫酸亞鐵0.03，硫酸錳0.03，磷酸三鈣5，瓊脂15，氯化鈉20，pH 9.0±0.2。

1.2 方法

1.2.1 耐鹽堿溶磷菌的篩選 取10 g鹽堿地土樣，加入90 mL無(wú)菌水渦旋振蕩10 min，靜置5 min后將上層懸液進(jìn)行倍比稀釋?zhuān)?00 μL涂布于pH 9.0、鹽度4%的R2A培養(yǎng)基平板上，30℃培養(yǎng)72 h，篩選耐鹽堿細(xì)菌，分離純化的耐鹽堿菌再轉(zhuǎn)接到無(wú)機(jī)磷固體培養(yǎng)基上，篩選有溶磷圈的菌落，分離純化后保藏于-80℃冰箱中。

1.2.2 菌株的形態(tài)學(xué)及生理生化特征測(cè)定 使用JEM1230透射電子顯微鏡進(jìn)行菌株的形態(tài)學(xué)觀察；接觸酶、脲酶、七葉苷水解、生長(zhǎng)的溫度、pH和NaCl范圍等生理生化指標(biāo)的測(cè)定參考文獻(xiàn)［14］進(jìn)行；使用BioMérieux公司的API ZYM試劑條測(cè)定19種酶的活性。

1.2.3 菌株的16S rDNA基因序列測(cè)定及系統(tǒng)進(jìn)化分析 利用細(xì)菌DNA提取試劑盒提取試驗(yàn)菌株的基因組DNA，采用細(xì)菌16S rDNA擴(kuò)增通用引物27F和1492R進(jìn)行PCR擴(kuò)增，程序?yàn)?5℃預(yù)變性5 min；94℃變性45 s，59℃退火45 s，72℃延伸1 min，共30個(gè)循環(huán)；72℃延伸10 min。PCR產(chǎn)物進(jìn)行瓊脂糖凝膠電泳檢測(cè)后，用DNA產(chǎn)物純化試劑盒純化PCR產(chǎn)物，與pGM-T載體連接后，轉(zhuǎn)化至TOP10感受態(tài)細(xì)胞中，在含 IPTG和X-Gal的LB平板上37℃過(guò)夜培養(yǎng)，挑取白色單菌落，驗(yàn)證并送上海生工進(jìn)行測(cè)序。將得到的16S rDNA序列與GenBank中的核酸數(shù)據(jù)進(jìn)行BLAST比對(duì)，得到與其同源性較高的對(duì)照菌株及序列，采用軟件Clustal X1.83進(jìn)行多序列比對(duì)分析［15］，得到序列之間的相似值；用軟件Mega 6計(jì)算出序列的系統(tǒng)進(jìn)化距離，采用鄰位相連法構(gòu)建進(jìn)化樹(shù)［16］，Bootstrap值為1 000。

1.2.4 菌株Y2R2的全基因組測(cè)序及分析 提取菌株Y2R2的基因組DNA，在北京百邁客公司的Nanopore和Illumina平臺(tái)進(jìn)行文庫(kù)構(gòu)建和基因組測(cè)序，測(cè)得4 221 Mb的原始數(shù)據(jù)，通多過(guò)濾接頭、短片段及低質(zhì)量數(shù)據(jù)后，共得到4 146 Mb的Clean Data，測(cè)序深度為847 X，拼接得到0 gap的菌株Y2R2基因組完成圖。參考文獻(xiàn)［17］的方法進(jìn)行平均核苷酸一致性（Average Nucleotide Identity，ANI）和DNA-DNA雜交（DNA-DNA hybridization，DDH）分析。

1.2.5 全細(xì)胞脂肪酸分析 按照Sasser［18］的方法進(jìn)行樣品的提取、皂化、甲基化和萃取，使用Agient 7890氣相色譜儀檢測(cè)脂肪酸組分，參考MIDI Sherlock細(xì)菌系統(tǒng)數(shù)據(jù)庫(kù)進(jìn)行分析鑒定。

1.2.6 可溶性磷含量的測(cè)定 取1.0 mL發(fā)酵液，于4℃，8 000 r/min離心5 min，取上清液稀釋至合適倍數(shù)，用鉬藍(lán)比色法測(cè)定可溶性磷的含量［19］，計(jì)算公式如下：Y=m×15.91A，式中Y為發(fā)酵液中可溶磷含量（mg/L），m為稀釋倍數(shù)，A為660 nm處吸光度值。

1.2.7 數(shù)據(jù)處理 可溶磷含量的測(cè)定，每個(gè)處理3個(gè)重復(fù)，采用Excel 2016進(jìn)行數(shù)據(jù)分析和作圖。

2 結(jié)果

2.1 耐鹽堿溶磷菌的篩選

在天津?yàn)I海新區(qū)鹽堿地土壤中分離得到多株耐鹽堿細(xì)菌，將這些菌種轉(zhuǎn)接到無(wú)機(jī)磷固體培養(yǎng)基上，篩選得到1株耐鹽堿溶磷菌（圖1），命名為Y2R2，進(jìn)一步試驗(yàn)表明該菌株可以耐受的最高pH為14，NaCl濃度為15%，具有較強(qiáng)的耐鹽堿性能。

圖1 菌株Y2R2產(chǎn)生的溶磷圈

2.2 菌株Y2R2的分類(lèi)鑒定

2.2.1 菌株Y2R2的形態(tài)特征 在R2A培養(yǎng)基平板上30℃培養(yǎng)2 d，菌株Y2R2形成邊緣整齊的乳白色菌落，質(zhì)地均勻、邊緣整齊，直徑1.5-2.3 mm。菌株Y2R2的細(xì)胞大小為0.4-0.7 μm × 1.7-2.6 μm，桿狀，無(wú)鞭毛（圖2），革蘭氏陰性。

圖2 菌株Y2R2的細(xì)胞形態(tài)（5000×）

2.2.2 菌株Y2R2的16S rDNA序列測(cè)定及系統(tǒng)進(jìn)化分析 擴(kuò)增得到菌株Y2R2的16S rDNA序列，長(zhǎng)度為1 460 bp，在GenBank/EMBL/DDBJ中的序列登記號(hào)為MK660018。將該序列在GenBank數(shù)據(jù)庫(kù)中進(jìn)行BLAST比對(duì)，得到與其同源性高的序列信息，進(jìn)而進(jìn)行Clustal X1.83比對(duì)，并利用軟件MEGA 6的鄰位相連法構(gòu)建系統(tǒng)進(jìn)化樹(shù)（圖3），結(jié)果表明菌株Y2R2屬于鹽單胞菌屬（Halomonas），與Halomonas huangheensisBJGMM-B45劃分在同一個(gè)簇內(nèi)，與該菌株的同源性也最高，為98.1%。

2.2.3 菌株Y2R2的全基因組測(cè)序及分析 菌株Y2R2的基因組長(zhǎng)度為4 895 267 bp，在GenBank/EMBL/DDBJ中的登記號(hào)為CP038437，基因組數(shù)據(jù)分析表明菌株Y2R2的G+C含量為57.02 mol%，基因組中非編碼RNA中有rRNA 15個(gè)，tRNA 65個(gè)；基因組中編碼基因?yàn)? 286個(gè)，總長(zhǎng)度4 327 563 bp，編碼基因的平均長(zhǎng)度為1 010 bp；有33個(gè)磷酸酶相關(guān)的編碼基因，其中有2個(gè)堿性磷酸酶編碼基因。利用ANI分析，可以從全基因組水平評(píng)估菌株Y2R2與同屬中相似菌株的親緣關(guān)系，將菌株Y2R2與16S rDNA序列最接近菌株H. huangheensisBJGMM-B45的基因組數(shù)據(jù)（CP013106）進(jìn)行ANI分析，按照BLAST（ANIb），MUMmer（ANIm）和OrthoANIu三種算法，ANI數(shù)值分別為74.9%，84.7%和75.1%，均小于95%的新菌判定閾值。按照基因組間距離計(jì)算（Genome-to-genome distance calculator，GGDC）方法，進(jìn)一步將菌種Y2R2與H.huangheensisBJGMM-B45的基因組進(jìn)行digital DDH分析，數(shù)值為22.2%。

圖3 菌株Y2R2的16S rDNA序列系統(tǒng)進(jìn)化樹(shù)

2.2.4 菌株Y2R2的全細(xì)胞脂肪酸分析 全細(xì)胞脂肪酸分析是細(xì)菌化學(xué)分類(lèi)鑒定的重要指標(biāo)，菌株Y2R2的細(xì)胞脂肪酸分析結(jié)果如圖4所示，經(jīng)MIDI Sherlock全自動(dòng)細(xì)菌鑒定系統(tǒng)分析，菌株Y2R2所含脂肪酸主要有C16：0（29.4%），C19：0 cyclo w8c（18.4%），C18：1 w7c/C18：1 w6c（17.8%），C16：1w7c/C16：1w6c（8.3%）；而與Y2R2同源性最高的菌株H. huangheensisBJGMM-B45的細(xì)胞脂肪酸組成主要為：C18：1w7c（32.4%），C16：0（24.6%），C19：0 cyclo w8c（12.6%），C16：1w7c/C16：1w6c（11.1%），C12：0 3-OH（7.8%），C12：0（5.2%）。二者在占比最高的脂肪酸組分上不同，分別為C16：0和C18：1w7c，在其他細(xì)胞脂肪酸比例和組成上也有明顯差別，說(shuō)明Y2R2與菌株H. huangheensisBJGMM-B45在化學(xué)分類(lèi)指標(biāo)上不同。

圖4 菌株Y2R2的全細(xì)胞脂肪酸分析

2.2.5 菌株Y2R2的生理生化特征 菌株Y2R2的生長(zhǎng)溫度范圍為10-45℃，生長(zhǎng)的NaCl范圍為0-15%，生長(zhǎng)的pH范圍是5-14。該菌株的其他生理生化指標(biāo)如表1所示，可以看出菌株Y2R2的脲酶為陽(yáng)性；酯酶（C14）為弱陽(yáng)性；酪蛋白水解、酪氨酸利用、胰凝乳蛋白酶、β-半乳糖苷酶、α-甘露糖苷酶、α-葡萄糖苷酶、α-半乳糖苷酶和β-葡萄糖苷酶為陰性；菌株Y2R2的這些生理生化指標(biāo)與H. huangheensis均不同。結(jié)合形態(tài)學(xué)、生理生化、基于16S rDNA系統(tǒng)進(jìn)化、全基因組數(shù)據(jù)和細(xì)胞脂肪酸分析，確定菌株Y2R2為鹽單胞菌屬的一個(gè)新種。

2.3 培養(yǎng)時(shí)間對(duì)菌株Y2R2溶磷量的影響

將菌種Y2R2接種于無(wú)機(jī)磷液體培養(yǎng)基中，每12 h取樣1次，檢測(cè)發(fā)酵液pH、細(xì)胞濃度和上清液中的可溶磷含量，結(jié)果如圖5所示，菌株Y2R2在發(fā)酵培養(yǎng)48 h達(dá)到最高菌量，細(xì)胞濃度為2.6×109CFU/mL，隨著培養(yǎng)時(shí)間的增加，Y2R2培養(yǎng)液的pH也隨之逐漸降低，36 h的pH為5.0，此后發(fā)酵液的pH穩(wěn)定在5.0左右，隨著發(fā)酵液pH的降低，溶磷量隨之增高，在36 h出現(xiàn)溶磷量降低的現(xiàn)象，是由于此時(shí)菌株Y2R2處于對(duì)數(shù)生長(zhǎng)期，快速繁殖的細(xì)胞需磷增加，將發(fā)酵液中的可溶磷重新攝入胞內(nèi)。在48 h培養(yǎng)液中的可溶磷含量達(dá)到最高167.9 mg/L，48 h后由于細(xì)胞代謝活力的降低，菌株Y2R2的溶磷量開(kāi)始下降。

表1 菌株Y2R2與H. huangheensis的生理生化特征比較

圖5 培養(yǎng)時(shí)間對(duì)菌株Y2R2溶磷量的影響

2.4 NaCl濃度對(duì)菌株Y2R2溶磷量的影響

在NaCl濃度0-18%的培養(yǎng)條件下，檢測(cè)培養(yǎng)48 h的發(fā)酵液，可以看出菌株Y2R2可以耐受較高的鹽濃度（圖6），當(dāng)培養(yǎng)基中含有9%的NaCl時(shí)，細(xì)胞濃度為不含NaCl對(duì)照的91.3%，說(shuō)明該菌株具有良好的耐鹽性能。當(dāng)NaCl濃度在12%以下時(shí)，發(fā)酵終點(diǎn)的pH均在5.8以下，且具有較好的溶磷效果，在9%的NaCl濃度時(shí)可溶磷含量最高，為239.2 mg/L。當(dāng)鹽濃度達(dá)到15%，菌株Y2R2的生長(zhǎng)被抑制，可溶磷含量為9%鹽度時(shí)的13.4%。

2.5 初始pH對(duì)菌株Y2R2溶磷量的影響

為研究菌株Y2R2在鹽堿條件下的溶磷能力，在NaCl濃度9%，初始pH為7-14的無(wú)機(jī)磷液體培養(yǎng)基中進(jìn)行發(fā)酵培養(yǎng)，結(jié)果如圖7所示，可以看出菌株Y2R2的最適生長(zhǎng)pH為8，在pH 9時(shí)細(xì)胞濃度為中性條件下的100.8%，pH 13時(shí)細(xì)胞濃度為中性條件下的70.6%。

圖6 NaCl對(duì)菌株Y2R2溶磷量的影響

在初始pH 7-9的范圍內(nèi)，發(fā)酵終點(diǎn)pH為5.1-5.4，具有較高的溶磷量，最高為pH 8時(shí)的247.6 mg/L；在初始pH 13時(shí)的溶磷量為115.4 mg/L，說(shuō)明菌株Y2R2在鹽堿條件下具有較強(qiáng)的溶磷能力。

3 討論

圖7 初始pH對(duì)菌株Y2R2溶磷量的影響

本研究從天津?yàn)I海新區(qū)鹽堿土中分離得到菌株Y2R2，基于16S rDNA的分子分類(lèi)標(biāo)準(zhǔn)［20］及基因組、生理生化和全細(xì)胞脂肪酸的多相分類(lèi)鑒定結(jié)果，試驗(yàn)菌株與H. huangheensis等［21］鹽單胞菌屬的模式菌種都存在明顯不同，鑒定該菌株是鹽單胞菌屬的一個(gè)新種。鹽單胞菌屬是Vreeland等［11］根據(jù)從太陽(yáng)能海水制鹽廠分離到的菌株Halomonas elongata提出的，屬于變形菌門(mén)（Proteobacteria），γ-變形菌綱（Gammaproteobacteria），海洋螺菌目（Oceanospirillales），鹽單胞菌科（Halomonadaceae）的一個(gè)屬，革蘭氏陰性桿菌，能在0.1%-20%的NaCl濃度下生長(zhǎng)，常在鹽堿土、工業(yè)鹵水和深海等高鹽環(huán)境中分離得到［22-23］。王永妍等［24］研究表明在養(yǎng)殖海水中，鹽單胞菌屬在細(xì)菌種群中占比11.4%，是主要的優(yōu)勢(shì)菌屬，汪輝等［25］從連云港贛榆附近海泥中分離到鹽單胞菌GY1，任世英等［26］從黃海海域分離到鹽單胞菌YSR-3。鹽單胞菌屬的一些菌株在降解芳香族污染物，生物脫氮和多糖合成方面均有潛在的應(yīng)用價(jià)值，但在溶磷方面的應(yīng)用尚未見(jiàn)報(bào)道。本研究分離得到的鹽單胞菌Y2R2具有耐鹽堿和溶磷能力，為該菌屬提供了新的微生物資源。

施用微生物菌肥是鹽堿土壤生物改良的重要措施，鄭曉梅等［27］將分離得到的黑曲霉制成菌劑，施加到實(shí)驗(yàn)室的鹽堿土壤中，56 d后土壤的pH值降低0.86-1.60，含鹽量也下降；楊美英等［28］將含腸桿菌屬NDW1和沙雷氏菌屬NDW3的外源溶磷菌液施用在鹽堿土中，盆栽試驗(yàn)表明水稻的凈光合速率提高，產(chǎn)量增加。目前微生物菌劑種類(lèi)繁多，在不同類(lèi)型土壤上的應(yīng)用效果相差較大，施入外來(lái)菌劑后，菌劑中的微生物在土壤中能否存活及繁殖數(shù)量受很多因素影響，如土壤環(huán)境變化的干擾、與土壤中其他土著微生物的生態(tài)位競(jìng)爭(zhēng)等［29］，外來(lái)菌株在鹽堿土壤中的定殖就更為困難。天津鹽堿地有4.9×105hm2，其中濱海新區(qū)鹽堿程度最重，鹽分主要以氯化物為主，土壤物理性狀差，植物生長(zhǎng)受到不同程度的影響［30］。本研究篩選到的菌株Y2R2是濱海新區(qū)鹽堿土中的土著微生物，能較好地適應(yīng)當(dāng)?shù)赝寥篮蜌夂颦h(huán)境，后續(xù)將研究其菌劑對(duì)土壤和植物生長(zhǎng)的影響，探索一種改良濱海鹽堿土壤，增加土壤養(yǎng)分，并促進(jìn)作物生長(zhǎng)的方法。

本研究發(fā)現(xiàn)，鹽單胞菌Y2R2可以在鹽堿條件下將難溶性無(wú)機(jī)磷化物轉(zhuǎn)變?yōu)橹参锟晌绽玫目扇苄杂行Я住＜?xì)菌對(duì)難溶性無(wú)機(jī)磷的解磷機(jī)制，普遍認(rèn)為是細(xì)胞向胞外分泌的各種有機(jī)酸降低了外界環(huán)境的pH，從而產(chǎn)生酸解效應(yīng)［31］。本實(shí)驗(yàn)中，在發(fā)酵液初始pH 10以下時(shí)，接入菌種培養(yǎng)48 h后，發(fā)酵液的pH均在5.7以下；在NaCl濃度12%以下時(shí)，發(fā)酵48 h后的pH均在5.6以下，說(shuō)明試驗(yàn)菌株在鹽堿條件下可產(chǎn)酸，從而造成難溶磷酸鹽的酸解。由于鹽單胞菌Y2R2的全基因組測(cè)序已經(jīng)完成，從基因組功能注釋和代謝分析可以找到與丙酮酸合成相關(guān)的基因5個(gè)，順烏頭酸酶編碼基因5個(gè)，檸檬酸合成酶編碼基因3個(gè)，蘋(píng)果酸脫氫酶編碼基因2個(gè)，葡萄糖醛酸合成所需的葡萄糖脫氫酶編碼基因2個(gè)。鑒于鹽單胞菌屬菌株產(chǎn)生有機(jī)酸的種類(lèi)和溶磷機(jī)制還未見(jiàn)報(bào)道，菌株Y2R2與溶磷相關(guān)的有機(jī)酸合成種類(lèi)和含量，還需要進(jìn)一步的試驗(yàn)測(cè)定，特別是鹽堿對(duì)這些有機(jī)酸合成的調(diào)控機(jī)制需要進(jìn)行深入研究。本研究篩選到的的菌株具有較高的耐鹽堿和溶磷能力，為濱海鹽堿地的土著微生物新種，能夠更好地適應(yīng)本地土壤和氣候環(huán)境，有望用于溶磷菌劑的制備，為濱海鹽堿地改良和解決鹽堿土壤固磷強(qiáng)度高的問(wèn)題提供了優(yōu)良菌種。

4 結(jié)論

菌株Y2R2為革蘭氏陰性桿菌，其16S rDNA序列與鹽單胞菌屬其他菌株的最高同源性為98.1%，基因組ANI和DDH值均小于新菌判定的閾值，結(jié)合生理生化和全細(xì)胞脂肪酸數(shù)據(jù)分析，鑒定菌株Y2R2是鹽單胞菌屬的一個(gè)新種。菌株Y2R2的生長(zhǎng)溫度范圍為10-45℃，具有一定的耐鹽堿能力，在0-15%的NaCl范圍，5-14的pH范圍內(nèi)均可生長(zhǎng)。菌株Y2R2在pH 8和NaCl濃度9%的無(wú)機(jī)磷培養(yǎng)基中，發(fā)酵48 h溶磷能力最高，達(dá)到247.6 mg/L。