• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    玉米秸稈低溫降解復(fù)合菌系降解能力及微生物組成研究*

    2020-11-10 06:33:36青格爾于曉芳高聚林王志剛胡萬吉鬧干朝魯胡樹平孫繼穎屈佳偉
    關(guān)鍵詞:菌系菌源桿菌屬

    青格爾, 于曉芳**, 高聚林**, 王志剛, 胡萬吉, 鬧干朝魯, 王 振, 胡樹平, 孫繼穎, 屈佳偉

    玉米秸稈低溫降解復(fù)合菌系降解能力及微生物組成研究*

    青格爾1,2, 于曉芳1,2**, 高聚林1,2**, 王志剛1,2, 胡萬吉1, 鬧干朝魯2,3, 王 振2,3, 胡樹平2,4, 孫繼穎1,2, 屈佳偉1,2

    (1. 內(nèi)蒙古農(nóng)業(yè)大學(xué)農(nóng)學(xué)院 呼和浩特 010019; 2. 內(nèi)蒙古自治區(qū)作物栽培與遺傳改良重點實驗室 呼和浩特 010019; 3. 內(nèi)蒙古農(nóng)業(yè)大學(xué)園藝與植保學(xué)院 呼和浩特 010019; 4. 內(nèi)蒙古農(nóng)業(yè)大學(xué)職業(yè)技術(shù)學(xué)院 包頭 014100)

    根據(jù)還田秸稈配施尿素的生產(chǎn)實際, 對玉米秸稈低溫高效降解復(fù)合菌系GF-20進(jìn)行氮源培養(yǎng)基馴化, 探明其菌種組成和功能多樣性及其與菌源菌種結(jié)構(gòu)差異, 完善復(fù)合菌系篩選方法, 促進(jìn)其開發(fā)利用。本文以低溫高效降解復(fù)合菌系GF-20為研究對象, 在硫酸銨和尿素不同配比下連續(xù)繼代培養(yǎng)10代, 獲得不同氮源菌系(硫酸銨處理N1, 硫酸銨和尿素混合處理N2-N5, 尿素處理N6), 測定其玉米秸稈降解率, 評價復(fù)合菌系秸稈降解效率; 采用MiSeq高通量測序?qū)赐寥罉悠芳安煌聪吕^代培養(yǎng)的復(fù)合菌系菌種組成和功能多樣性進(jìn)行研究。結(jié)果顯示N2處理玉米秸稈降解率顯著高于其他處理; 菌源土壤的Alpha多樣性指數(shù)顯著高于繼代培養(yǎng)后的復(fù)合菌系, 不同處理間N2處理顯著高于其他氮處理; 菌源和復(fù)合菌間以及不同氮處理間菌種組成具有顯著差異, N2處理菌種組成多樣性較高, 菌群結(jié)構(gòu)更加豐富、均勻, 且碳水化合物的代謝通路豐度較高。菌源經(jīng)限制性繼代篩選后得到了參與玉米秸稈降解過程的功能菌, 能有效提高秸稈降解率, 其在硫酸銨和尿素氮源為0.16%+0.04%的條件下, 菌系的秸稈降解效率較高, 這為復(fù)合菌的生產(chǎn)實際開發(fā)利用提供了理論依據(jù)。

    菌源; 復(fù)合菌系; 氮源; 高通量測序; 菌種組成多樣性; 玉米秸稈降解

    秸稈還田是提高秸稈資源綜合利用率、平衡農(nóng)田生態(tài)系統(tǒng)和實現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展的有效措施[1]。由于玉米()秸稈木質(zhì)纖維素的結(jié)晶區(qū)域組成, 使得還田玉米秸稈難以降解[2]; 而北方高寒地區(qū)秋冬季低溫、干燥等氣候因素進(jìn)一步限制了還田玉米秸稈的有效腐解。因此, 在低溫條件下施加外源菌劑加速分解木質(zhì)纖維素對玉米秸稈資源的開發(fā)利用具有重要意義。根據(jù)秸稈酶解過程以及微生物間的協(xié)同關(guān)系, 直接從自然環(huán)境中篩選或組配秸稈高效降解復(fù)合菌群, 用以提高還田秸稈的降解效率是目前秸稈高效降解的有效途徑之一[3-4]。其中, 菌源的選擇、復(fù)合菌群的篩選馴化方法對于獲得秸稈高效降解微生物成敗至關(guān)重要。前人研究表明, 根據(jù)應(yīng)用目的從相似的生態(tài)環(huán)境中可篩選到更高效適宜的菌種[5]。Li等[6]和Feng等[7]以高溫期堆肥、腐爛作物秸稈和秸稈還田土壤為菌源, 篩選得到一系列中高溫復(fù)合菌系, 其秸稈降解率達(dá)30%~50%, 被應(yīng)用于堆肥腐熟或南方高溫區(qū)秸稈還田。

    本研究團隊為了獲得玉米秸稈低溫高效降解菌以用于北方低溫種植區(qū)玉米秸稈原位促腐還田, 從低溫(?7~8 ℃)生態(tài)環(huán)境中采集多年秸稈還田土壤作為玉米秸稈低溫高效降解復(fù)合菌系篩選的菌源材料,經(jīng)長期的碳源限制性繼代培養(yǎng)與低溫馴化培養(yǎng)獲得一組微生物復(fù)合菌系GF-20[8], 其實驗室秸稈降解率達(dá)30%。但菌源材料與復(fù)合菌系GF-20間的菌種組成相關(guān)性及不同尿素含量處理對菌種組成的差異性尚不明確。這無疑關(guān)系到該菌系秸稈高效降解機理的揭示, 以及其能否被繼續(xù)開發(fā)、應(yīng)用于生產(chǎn)實際。因此, 本研究以菌源樣品、篩選得到的復(fù)合菌系GF-20及其尿素氮源馴化后的系列菌系為材料, 通過對各復(fù)合菌系菌種多樣性、玉米秸稈分解能力和菌群結(jié)構(gòu)特性等指標(biāo)的測定及分析, 揭示復(fù)合菌系GF-20與菌源材料間, 以及不同尿素氮源含量馴化得到的復(fù)合菌系間的菌種異同, 通過與COG(蛋白質(zhì)直系同源簇)和KEGG(代謝通路)數(shù)據(jù)庫比對及菌株功能性分析從而確定復(fù)合菌系中的關(guān)鍵功能菌,以及菌系適宜的尿素氮源條件, 為該復(fù)合菌系的開發(fā)利用提供理論依據(jù)和技術(shù)支持, 也給其他高效降解菌系的篩選和利用提供借鑒。

    1 材料與方法

    1.1 試驗時間與地點

    試驗于2018—2019年在內(nèi)蒙古農(nóng)業(yè)大學(xué)玉米中心微生物實驗室(內(nèi)蒙古包頭市土默特右旗溝門鎮(zhèn))進(jìn)行。

    1.2 供試材料

    菌源土壤(編號為Q, 為本文復(fù)合菌系GF-20的篩選菌源)取自內(nèi)蒙古通遼市常年秸稈還田土壤(43.37°N、122.16°E, 海拔400 m, 年平均氣溫5.5 ℃)和吉林省長春市連續(xù)秸稈還田土壤(43.82°N、125.32°E, 海拔250 m, 年平均氣溫4.8 ℃)的混合樣品, 供試復(fù)合菌系為本實驗室篩選馴化的玉米秸稈降解復(fù)合菌系GF-20(編號為N1)[8]。玉米秸稈取自內(nèi)蒙古農(nóng)業(yè)大學(xué)蒙西綜合試驗站(內(nèi)蒙古包頭市土默特右旗溝門鎮(zhèn))試驗田收獲的玉米秸稈(C/N為48.98, 纖維素、半纖維素和木質(zhì)素含量分別為51.24%、32.72%和9.15%, 全氮6.67 g·kg-1, 全磷2.78 g·kg-1, 全鉀11.48 g·kg-1), 洗凈烘干后剪成1~2 cm短節(jié)。濾紙選用Whatman No.Ⅰ濾紙, 裁成1 cm×10 cm大小。

    1.3 培養(yǎng)條件

    以基礎(chǔ)培養(yǎng)基[(NH4)2SO42.0 g、K2HPO41.0 g、MgSO4·7H2O 0.05 g、CaCO32.0 g、NaCl 0.2 g、蒸餾水1 L]為繼代培養(yǎng)基, 復(fù)合菌以5%~10%(/)的接種量接入裝有40 mL基礎(chǔ)培養(yǎng)基和1 g玉米秸稈碳源的100 mL三角瓶中, 置于10 ℃恒溫條件下靜置培養(yǎng)[9]。

    1.4 試驗設(shè)計

    基礎(chǔ)培養(yǎng)基中氮源設(shè)硫酸銨+尿素, 含量分別為0.16%+0.04%(N2)、0.12%+0.08%(N3)、0.08%+0.12%(N4)、0.04%+0.16%(N5)、0+0.2%(N6)。對復(fù)合菌系GF-20進(jìn)行尿素馴化處理, 原始培養(yǎng)基為硫酸銨0.2%+尿素0 (N1), 在以上處理條件下連續(xù)繼代培養(yǎng)10代以上, 獲得不同菌系, 每處理均3次重復(fù)。

    1.5 主要測定指標(biāo)與方法

    1)復(fù)合菌系玉米秸稈降解率的測定。將所得復(fù)合菌系以5%的接種量接入玉米秸稈培養(yǎng)基[9]中, 以不接菌的為對照(CK), 在10 ℃恒溫條件下培養(yǎng)15 d和30 d后, 采用失重法測定玉米秸稈降解率, 3次重復(fù)。秸稈降解率(%)=(0-1)/0×100%。0表示接種前培養(yǎng)基中的秸稈重(g),1表示培養(yǎng)結(jié)束烘干后降解剩余物重(g)。

    2)MiSeq 16S rDNA擴增子測序。采用土壤基因組DNA提取試劑盒和細(xì)菌基因組DNA提取試劑盒(中國, 天根生化科技有限公司)提取菌源土壤和各復(fù)合菌系的基因組DNA, 利用1%瓊脂糖凝膠電泳檢測。獲得的基因組DNA由北京奧維森基因科技有限公司(Allwegene Tech. Co. Ltd., Beijing)完成MiSeq高通量測序, 對獲得的數(shù)據(jù)通過序列拼接、過濾和去嵌合體得到優(yōu)化序列, 進(jìn)行OTU(Operational Taxonomic Units)聚類及各分類水平注釋。

    3)功能注釋?;跍y序得到的序列, 與COG(蛋白質(zhì)直系同源簇)和KEGG(代謝通路)數(shù)據(jù)庫進(jìn)行比對[10]。

    4)數(shù)據(jù)統(tǒng)計分析。利用QIIME(v1.8.0)軟件以97%相似度水平將序列進(jìn)行聚類物種分類的OTU, 獲得物種分類水平的比例信息及不同分類水平的群落結(jié)構(gòu); 通過Uclust(Version 1.8)、Usearch(Version 8.1.1861)等數(shù)據(jù)庫進(jìn)行物種比對, 得到各水平的分類信息, 進(jìn)行樣本組成及樣本間群落結(jié)構(gòu)差異分析; 采用R(v3.4.4)軟件繪制各樣品在各分類水平的比較圖, 基于聚類結(jié)果, 進(jìn)行Alpha多樣性分析和Beta多樣性分析。

    2 結(jié)果與分析

    2.1 不同氮源下復(fù)合菌系玉米秸稈降解率

    不同氮源條件下復(fù)合菌系GF-20的玉米秸稈降解率如圖1所示。不同氮源處理降解率均顯著高于CK; 不同處理隨著尿素比例的增加, 復(fù)合菌系GF-20玉米秸稈降解率呈先上升后下降趨勢。培養(yǎng)至15 d和30 d, 均以N2條件下復(fù)合菌系的玉米秸稈降解率最高, 分別為22.64%和37.58%, 較CK高15.85%和21.43%, 且顯著高于其他處理。

    圖1 不同氮源處理復(fù)合菌系玉米秸稈降解率

    N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例為分別0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%、0+0.2%。CK為不加菌對照, 培養(yǎng)條件同N1。不同小寫字母表示不同處理間差異顯著。The nitrogen resources of N1, N2, N3, N4, N5 and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16%, and 0+0.2%, respectively. CK is non-inoculated microbial consortium (control), and the culture conditions are the same as N1. Different lowercase letters mean significant differences among different treatments.

    2.2 菌源樣品及不同氮源下復(fù)合菌系A(chǔ)lpha多樣性指數(shù)分析

    所有樣品及處理覆蓋深度指數(shù)均在0.99以上, 說明本次測序結(jié)果基本能代表樣品的真實情況。菌源土壤平均OUT數(shù)為1 705, 顯著高于篩選獲得的復(fù)合菌系, 處理N1-N6平均OUT數(shù)為58; 由Shannon、PD whole tree、Observed species和Chao指數(shù)可知, 菌源土壤顯著高于處理N1-N6; 不同氮源培養(yǎng)基獲得的復(fù)合菌系表現(xiàn)為混合氮源處理(N2-N5)高于單一氮源處理(N1, N6), 且N2處理顯著高于其他處理, 其Shannon、PD whole tree、Observed species和Chao多樣性指數(shù)分別為4.36、6.14、59.93和62.45(圖2)。

    圖2 菌源樣品及不同氮源處理復(fù)合菌系A(chǔ)lpha多樣性指數(shù)

    Q為菌源土壤。N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例為分別0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%、0+0.2%。Q is original microbe source. The nitrogen resources of N1, N2, N3, N4, N5 and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16%, and 0+0.2%, respectively.

    2.3 菌源樣品及不同氮源下復(fù)合菌系細(xì)菌群落結(jié)構(gòu)分析

    對所有樣品測序結(jié)果進(jìn)行物種比對, 得到各水平的分類信息。門分類水平, 變形菌門(Proteobacteria)在所有樣本中豐富度最高(圖3a)。菌源樣品中變形菌門豐度為36.98%, 其他優(yōu)勢門是放線菌門(Actinobacteria)和綠彎菌門(Chloroflexi), 豐度分別為24.06%和10.64%; 復(fù)合菌系GF-20中變形菌門豐度為85.52%, 擬桿菌門(Bacteroidetes)豐度為13.86%。綱分類水平, 在所有樣本中α變形菌(Alphaproteobacteria)β變形菌(Betaproteobacteria)和γ變形菌(Gammaproteobacteria)豐度較高。而黃桿菌綱(Flavobacteriia)在復(fù)合菌中豐富度較高, 在菌源中僅有0.32%(圖3b)。屬分類水平(圖4), 菌源主要由(5.00%)、(3.73%)、(3.13%)、節(jié)細(xì)菌屬(, 2.26%)、(2.24%)、壤霉菌屬(, 2.01%)組成, 復(fù)合菌系主要由纖維弧菌屬()固氮螺菌屬()黃桿菌屬()土地桿菌屬()假單孢菌屬()纖維菌屬()等組成, 且不同氮源處理豐度具有顯著差異。說明復(fù)合菌系及菌源樣品的菌群結(jié)構(gòu)差異較大, 長期限制性繼代培養(yǎng)及低溫馴化培養(yǎng)顯著影響菌系菌種組成, 大量菌屬被淘汰。

    不同氮源處理間, 在門分類水平, N1、N2和N6中擬桿菌門豐度較高, 為13.86%、37.01%和17.72%; 放線菌門在N2-N5中含量分別為1.18%、5.71%、6.29%和26.11%, 表現(xiàn)為隨著尿素含量的增加而增加。屬水平, N1-N6處理復(fù)合菌系的細(xì)菌菌群結(jié)構(gòu)差異顯著(圖3c, 圖4), 復(fù)合菌系GF-20在N1條件下的優(yōu)勢菌屬為纖維弧菌屬(74.54%), 其次為固氮螺菌屬(7.40%); N2的優(yōu)勢菌屬為纖維弧菌屬(22.24%)、黃桿菌屬(9.35%)、土地桿菌屬(8.60%)和norank_c__(9.03%); N3主要由纖維弧菌屬(14.3%)、假單胞菌屬(55.72%)和纖維菌屬(5.48%)組成; N4主要由假單胞菌屬(56.82%)、纖維菌屬(5.80%)、德沃斯氏菌屬(, 5.44%)、芽殖桿菌屬(, 4.57%)和纖維弧菌屬(4.44%)組成; N5主要由纖維菌屬(25.98%)、(28.11%)、食酸菌屬(, 10.97%)和噬氫菌屬(, 6.36%)組成; N6的優(yōu)勢菌屬為纖維弧菌屬(31.13%)、短波單胞菌屬(, 18.5%)、黃桿菌屬(12.8%)和無色桿菌屬(12.56%)等??梢娔蛩睾繉?fù)合菌系GF-20的細(xì)菌菌群結(jié)構(gòu)影響較大, 培養(yǎng)基中加入尿素后使優(yōu)勢菌種種類增加, 混合氮源效果尤其明顯。綜上所述可知, 菌源樣品和所篩選獲得的復(fù)合菌系, 以及不同氮源下的菌系間菌種組成差異顯著, 主要體現(xiàn)在放線菌門、擬桿菌門、變形菌門的纖維弧菌屬、假單胞菌屬、、黃桿菌屬、無色桿菌屬、纖維菌屬等參與秸稈降解微生物豐度的差異。

    圖3 菌源樣品及不同氮源處理復(fù)合菌系門(a)和綱(b)水平群落結(jié)果及屬水平(c)Heatmap圖分析

    將在所有樣本中豐度占比均小于0.02%的物種歸為Others。Q為菌源土壤。N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例分別為0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%和0+0.2%。Microbes with less than 0.02% abundance in all samples are classified as others. Q is original microbe source. The nitrogen resources of N1, N2, N3, N4, N5, and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16% and 0+0.2%, respectively.

    2.4 不同氮源下復(fù)合菌系菌群群落結(jié)構(gòu)Heatmap分析

    根據(jù)所有樣品在屬水平的物種注釋及豐度信息, 選取豐度top 20的屬, 繪制屬水平的Heatmap圖(圖3c)。在屬水平上, 不同氮處理復(fù)合菌系優(yōu)勢菌屬不同, 混合氮源處理N2-N5的優(yōu)勢屬數(shù)量多于單一氮源N1和N6處理, N1處理中的固氮螺菌屬等屬豐度在N6中顯著降低, 而豐度較低的如無色桿菌屬短波單胞菌屬噬氫菌屬根瘤菌屬()貪噬菌()等屬在混合氮源處理中成為優(yōu)勢屬。

    由圖3c和圖4可知, 復(fù)合菌系GF-20在N1條件下, 纖維弧菌屬的豐度為74.50%, 但在含有不同尿素比例氮源條件下繼代培養(yǎng)數(shù)量顯著變化, 尤其在混合氮源條件下隨著尿素含量的增加而逐漸減少至0.81%; 德沃斯氏菌屬在N1條件下豐度為0.99%, 但在N4處理中增加至5.61%; 黃桿菌屬、地桿菌屬、短波單胞菌屬、無色桿菌屬等在N1條件下豐度均很低, 但在含有尿素氮源的菌系中大量存在, 成為優(yōu)勢菌屬。說明限制性繼代培養(yǎng)是有選擇性的進(jìn)行篩選, 大量富集與秸稈纖維素降解相關(guān)的菌屬, 提高菌屬數(shù)量比, 從而促進(jìn)秸稈的降解。纖維菌屬在N1條件下豐度僅僅為0.03%, 但在混合氮源處理下含量逐漸增加, 在N5處理下含量達(dá)27.31%; 假單胞菌屬在N1條件下含量為0.04%, 但在混合氮源N3和N4處理下數(shù)量劇增至57.68%和59.25%; 無色桿菌屬和鞘脂單胞菌屬()在N1處理中為0.08%和0.06%, 但隨著尿素含量的增加, 含量逐漸增加, N6處理增加至12.84%和3.59%;和在N1處理中含量為0.70%和0.14%, 在混合氮源N2處理中增加至3.91%和2.35%;在混合氮源處理N5中含量最高, 為29.75%;在N1處理中含量為1.68%, 在處理N2中增加至5.42%, 隨著尿素氮源的增加含量逐漸降低, 在N6處理中減少至0; 固氮螺菌屬和金黃桿菌屬()在N1處理中的豐度分別為7.2%、3.78%和2.86%, 但隨著尿素氮源的增加, 其含量逐漸減少, 尿素處理N6的含量減少至0。由此可知, 硫酸銨和尿素含量對復(fù)合菌系GF-20繁殖有明顯的影響, 由Alpha多樣性指數(shù)及菌群組成(圖2, 圖4)可知, N2處理中富集了更多的參與秸稈降解的微生物, 菌種組成多樣性較高, 菌群結(jié)構(gòu)更加豐富。

    圖4 菌源及不同氮源處理復(fù)合菌系細(xì)菌群落結(jié)構(gòu)分布(屬水平)

    Q為菌源土壤。N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例分別為0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%、0+0.2%。Q is original microbe source. The nitrogen resources of N1, N2, N3, N4, N5 and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16%, 0+0.2%, respectively.

    2.5 不同氮源下復(fù)合菌系菌群結(jié)構(gòu)PCoA分析及Rank Abundance分析

    根據(jù)不同氮源處理菌種組成及進(jìn)化關(guān)系進(jìn)行PCoA分析, 結(jié)果表明樣品分布較為分散(圖5a), 說明不同氮源處理菌系間物種豐度及差異度均存在差異, 其中GF-20(N1)、N2及N6與N4和N5處理的距離較遠(yuǎn), 而N3、N4和N5的距離較近, 再次說明不同氮源處理, 顯著影響復(fù)合菌系菌種物種組成, 且相同氮源的組成菌群結(jié)構(gòu)更加相似。Rank-Abundance說明物種豐富度和均勻度。由圖5b可知, N2處理水平方向曲線寬度最寬且平滑, 結(jié)合OTU指數(shù)分析發(fā)現(xiàn), N2處理OUT數(shù)為60±0.88, 顯著高于其他處理, 較N1及N3-N6分別高13.12%、8.38%、5.84%、13.12%和47.15%(圖2), 說明N2處理物種組成更加豐富和均勻。

    2.6 不同氮源下復(fù)合菌系菌群Venn圖分析

    進(jìn)一步在屬水平進(jìn)行Venn圖分析, 統(tǒng)計樣本中所共有和獨有菌屬數(shù)目, 直觀體現(xiàn)各樣本的菌屬組成相似性及重疊情況。由圖6可知, 有32個屬為6個處理共有的, 被稱為核心微生物, 直接或間接參與玉米秸稈降解代謝, 主要由放線菌、α變形菌、變形菌、擬桿菌綱的纖維弧菌屬、假單胞菌屬、黃桿菌屬、纖維菌屬、德沃斯氏菌屬、無色桿菌屬等組成。

    2.7 不同氮源下復(fù)合菌系功能多樣性預(yù)測分析

    將不同氮源處理復(fù)合菌系測序獲得的序列與COG和KEGG數(shù)據(jù)庫進(jìn)行系統(tǒng)功能預(yù)測分析。通過COG分析對潛在功能進(jìn)行分類, 結(jié)果表明amino acid transport and metabolism (氨基酸運輸/代謝, E)、general function prediction only (一般功能預(yù)測, R)、signal transduction mechanisms (信號轉(zhuǎn)導(dǎo)機制, T)、function unknown (未知功能, S)的注釋豐度最多為主導(dǎo)功能(圖7)。將24個類別概括為4類: 信息存儲與處理(第1類)、細(xì)胞過程和信號傳導(dǎo)(第2類)、代謝(第3類)以及功能較差類別(第4類), 其中carbohydrate transport and metabolism(碳水化合物轉(zhuǎn)運和代謝)的相對豐度在N1和N2處理明顯高于其他處理, 說明其通過增加碳水化合物的運輸和代謝功能, 促進(jìn)秸稈的降解。

    圖5 不同氮源處理復(fù)合菌系屬水平PCoA (a)和RankAbundance (b)分析

    N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例分別為0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%和0+0.2%。The nitrogen resources of N1, N2, N3, N4, N5 and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16%, and 0+0.2%, respectively.

    圖6 不同氮源處理復(fù)合菌系屬水平菌種分布比較的Venn圖分析

    N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例分別為0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%和0+0.2%。The nitrogen resources of N1, N2, N3, N4, N5 and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16%, and 0+0.2%, respectively.

    通過KEGG數(shù)據(jù)庫比對(表1)發(fā)現(xiàn), metabolism (代謝)功能注釋到的結(jié)果最多, 占全功能分類的60%以上, 不同處理間N2處理明顯高于其他處理。為進(jìn)一步分析代謝途徑, 在KEGG第3分類水平進(jìn)行分析, 選取豐度大于1.00%的路徑進(jìn)行統(tǒng)計(圖8), 共獲得21條代謝路徑, ABC transporters(ABCT)和two-component system(TCSP)注釋最多,其中9條代謝途徑豐度呈N2大于其他處理, 其中包含fructose and mannose metabolism(果糖和甘露糖代謝)、starch and sucrose metabolism(淀粉和蔗糖代謝)和amino sugar and nucleotide sugar metabolism(氨基糖和核苷酸糖代謝)等相關(guān)碳水化合物的代謝路徑。再次說明復(fù)合菌系在氮源為硫酸銨和尿素含量為0.16+0.04% (N2)條件下, 增加碳水化合物代謝功能, 促進(jìn)秸稈物質(zhì)的降解。

    圖7 不同氮源處理復(fù)合菌系COG功能注釋相對豐度

    N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例分別為0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%和0+0.2%。The nitrogen resources of N1, N2, N3, N4, N5 and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16%, and 0+0.2%, respectively.

    表1 不同氮源處理復(fù)合菌系KEGG注釋結(jié)果的第2類水平統(tǒng)計

    續(xù)表1

    N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例分別為0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%和0+0.2%。RA%指相對豐度, log指注釋數(shù)量的log10。The nitrogen resources of N1, N2, N3, N4, N5 and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16%, 0+0.2%, and respectively. RA% refers to relative abundance, and log refers to log10 of the number of annotations.

    3 討論

    本研究中菌源材料取自內(nèi)蒙古通遼市常年秸稈還田和吉林長春市連續(xù)秸稈還田混合土壤樣品, 主要由放線菌(24.06%)、變形菌(36.98%)、酸桿菌(9.71%)、綠彎菌(10.64%)、擬桿菌(6.94%)、芽單胞菌(5.74%)等門組成, 富含降解纖維素、半纖維素、木質(zhì)素的菌屬, 如鞘氨醇單胞菌屬(, 0.99%)、芽孢桿菌屬(, 0.22%)、梭菌屬(, 0.1%)、纖維弧菌屬(, 0.01%)、(0.02%)、纖維菌屬(, 0.07%)、(0.04%)、(0.14%)等, 從菌源樣品中經(jīng)長期限制性繼代培養(yǎng)篩選出了玉米秸稈低溫高效降解菌系GF-20, 其在低溫10 ℃下的秸稈降解率達(dá)30%。復(fù)合菌系篩選馴化的最終目的是應(yīng)用到秸稈還田土壤中, 通常以生態(tài)條件與目標(biāo)菌系相近的常年秸稈還田土壤作為菌源篩選材料, 提高篩選效率及應(yīng)用效果。

    圖8 不同氮源處理復(fù)合菌系KEGG主要代謝途徑的相對豐度(>1.00%)

    N1、N2、N3、N4、N5和N6的氮源為硫酸銨+尿素, 比例分別為0.2%+0、0.16%+0.04%、0.12%+0.08%、0.08%+0.12%、0.04%+0.16%和0+0.2%。The nitrogen resources of N1, N2, N3, N4, N5 and N6 are ammonium sulfate+urea with ratios of 0.2%+0, 0.16%+0.04%, 0.12%+0.08%, 0.08%+0.12%, 0.04%+0.16%, and 0+0.2%, respectively.

    3.1 限制性繼代培養(yǎng)過程中保留菌屬分析

    限制性繼代培養(yǎng)篩選過程中保留了在低溫條件下直接或間接參與木質(zhì)纖維素降解代謝途徑的菌, 菌源樣品及復(fù)合菌系數(shù)量比高于0.1%的共有菌屬有9個(圖3, 圖4, 圖6), 其中纖維弧菌屬可分泌3~4種纖維素降解酶系, 具有木聚糖酶基因, 降解羧甲基纖維素、纖維素、半纖維素、幾丁質(zhì)等, 如混合纖維弧菌(subsp.)、J3-8、等[11-12], 其中混合纖維弧菌屬于產(chǎn)纖維素酶類的兼性厭氧微生物, 分泌的纖維素酶屬于典型的低溫纖維素酶[13]; 土地桿菌屬和黃桿菌屬大多源自于土壤, 可利用秸稈降解中間代謝產(chǎn)物促進(jìn)秸稈分解, 如sp.可分泌葡糖苷酶、糖苷水解酶, 可利用碳水化合物、醇類和糖苷類等[14-15],sp.sp.分泌過氧化氫酶、氧化酶、葡糖苷酶、半乳糖苷酶等, 利用丙酸、乳酸等秸稈降解中間產(chǎn)物, 消除底物反饋抑制作用[16-17]; 復(fù)合菌系中還含有鞘氨醇單胞菌屬根瘤菌氫噬菌屬柄桿菌屬()等含量低但起到關(guān)鍵作用的菌, 其中鞘氨醇單胞菌屬和鞘脂菌屬()是降解多環(huán)芳烴(PAHs)的重要功能微生物, 如sp. GD542參與石油降解途徑[18],分泌過氧化氫酶和氧化酶, 代謝產(chǎn)生丙酸、琥珀酸、L-丙氨酰胺等有機酸類, 在秸稈分解初期調(diào)節(jié)培養(yǎng)體系pH值[19]; 氫噬菌屬通常分泌水解酶, 如為兼氧細(xì)菌, 分泌氧化酶[20]; La等[21]報道sp.等菌分泌氧化酶、過氧化氫酶, 促進(jìn)秸稈物質(zhì)的分解; 德沃斯氏菌屬的sp.等分泌過氧化氫酶, 利用木糖、戊醛糖、纖維素等[22-23]; 柄桿菌屬中的能降解水楊苷[24],分泌糖苷酶、均可在4 ℃低溫條件下生長, 分泌水解酶類,以纖維二糖、木糖為碳源產(chǎn)酸, 分泌葡萄糖苷酶[25]; 根瘤菌屬大多存在于根瘤菌中, 既可合成纖維素也降解木質(zhì)素, 其中YS-1r具有木質(zhì)素降解能力, 能夠降解各種木質(zhì)素單體、二聚體以及柳枝草和紫花苜蓿的天然木質(zhì)素, 并且分泌木質(zhì)素過氧化物酶(LiP)[26-27]。由此可知, 復(fù)合菌系GF-20中保留了菌源土壤中的功能菌, 通過限制性培養(yǎng)使其大量擴繁, 大大提高豐度, 從而促進(jìn)玉米秸稈的降解。

    3.2 復(fù)合菌系菌種功能多樣性分析

    復(fù)合菌系的限制性繼代篩選過程既是淘汰與秸稈降解無關(guān)或不適應(yīng)培養(yǎng)基條件的菌屬過程, 也是富集木質(zhì)纖維素降解菌的過程。復(fù)合菌系中保留的菌屬均直接或間接參與秸稈降解代謝過程, 或是分泌纖維素、半纖維素、木質(zhì)素酶等直接降解木質(zhì)纖維素, 如纖維弧菌屬德沃斯氏菌屬根瘤菌屬鞘氨醇桿菌屬寡養(yǎng)單胞菌()等; 或是分泌葡萄糖苷酶、水解酶、過氧化氫酶等消耗中間代謝產(chǎn)物纖維二糖、木糖、醛糖等, 促進(jìn)秸稈的持續(xù)降解, 如氫噬菌屬土地桿菌屬黃桿菌屬、柄桿菌屬、、類芽孢桿菌屬()等; 有些菌產(chǎn)生氧化酶類, 調(diào)節(jié)發(fā)酵體系酸堿度維持中性水平, 協(xié)調(diào)降解秸稈, 如鞘氨醇桿菌屬黃桿菌屬等分泌氧化酶類, 利用乙酸、乳酸等中間產(chǎn)物, 調(diào)節(jié)發(fā)酵體系酸堿度維持中性水平, 協(xié)調(diào)降解秸稈; 有些菌可能與秸稈降解沒有直接關(guān)系, 是因為適應(yīng)篩選培養(yǎng)基而得以保留。

    通過COG和KEEG數(shù)據(jù)庫對復(fù)合菌系功能多樣性進(jìn)行預(yù)測比對發(fā)現(xiàn), ABCT和TSCP相對含量最多。就功能概況而言, ABCT可將氨基酸、脂質(zhì)、脂糖、肽和無機離子(如金屬離子)轉(zhuǎn)運至細(xì)胞外部, 從細(xì)胞中攝取寡糖、多糖、單糖(例如纖維二糖、纖維糊精、葡萄糖、半乳糖和甘露糖); 而TCSP感知細(xì)胞外可溶性糖的存在, 并調(diào)節(jié)大多數(shù)糖苷水解酶(GH)和相關(guān)的ABCT, 調(diào)節(jié)糖類的轉(zhuǎn)運, 以促進(jìn)復(fù)雜植物材料(例如半纖維素)的利用, 并將各種結(jié)構(gòu)的產(chǎn)物(例如木糖)轉(zhuǎn)運到細(xì)胞中[28]。代謝功能注釋相對豐度大多表現(xiàn)為N2處理含有更多的碳水化合物代謝功能, 包括果糖和甘露糖代謝、淀粉和蔗糖代謝以及氨基糖和核苷酸糖代謝等, 促進(jìn)秸稈降解中間產(chǎn)物的代謝, 從而提高玉米秸稈木質(zhì)纖維素的降解(圖1); 而部分代謝基因在N6處理中注釋豐度明顯高于其他處理, 說明在尿素氮源條件下, 復(fù)合菌系通過促進(jìn)自身代謝, 從而適應(yīng)培養(yǎng)環(huán)境(表1)。

    3.3 復(fù)合菌系間差異菌屬分析

    菌源土壤與其經(jīng)過長時間限制性繼代培養(yǎng)篩選馴化獲得的復(fù)合菌系細(xì)菌群落結(jié)構(gòu)具有顯著差異, 限制性繼代培養(yǎng)過程是微生物淘汰的過程, 具有隨機性。同樣, 將篩選獲得的復(fù)合菌系培養(yǎng)于不同條件下時, 細(xì)菌群落結(jié)構(gòu)也有差異, 說明菌系內(nèi)部關(guān)鍵菌株的數(shù)量比例隨著外界培養(yǎng)條件的變化而變化(圖5, 圖6)。根霉菌屬()在N2中含量較高,可直接利用纖維二糖、木聚糖等參與秸稈降解代謝途徑[29],而可利用纖維二糖、木糖等代謝中間產(chǎn)物[30]; 短波單胞菌屬和無色桿菌屬在N1和N2中豐度為0.05%、0.08%和0.59%、3.37%, 在N6中為18.90%和12.84%, 說明尿素為唯一氮源促進(jìn)其生長,等均可分泌纖維素酶, 參與秸稈物質(zhì)的降解, 短波單胞菌屬的sp.sp.等可分泌葡萄糖苷酶, 參與秸稈降解代謝[31],sp.sp.等可利用纖維二糖等代謝產(chǎn)物; 司美茹等[32]從被石油污染的土壤和焦化污泥中分離出無色桿菌屬的XL和SQ-1, 幾乎能徹底降解C12?C23和C27?C43的烷烴和多環(huán)芳香烴,能氧化木糖, 具有過氧化氫酶和氧化酶活性[32-33]; 纖維菌屬在混合氮源條件下豐度較高, 許多纖維菌屬的細(xì)菌應(yīng)用于秸稈生物降解和生物修復(fù)領(lǐng)域, 其中分泌的纖維素酶活性高, 包括β-葡萄糖苷酶(GH1, GH3)、內(nèi)切酶(GH6, GH9)、β-1,4-葡聚糖酶(GH48, GH5)和纖維二糖磷酸化酶(GH94)[33-34];在N2中的豐度顯著高于其他處理菌系,菌屬中的NS114T、Y620-1等具有纖維素降解能力[34-35]。

    根瘤菌屬等6個豐度高于0.1%的屬在N1和N2中存在, 但N6中未檢測到, 可能是由于以上菌屬不能以尿素為唯一氮源, 在繼代過程中逐漸被淘汰或含量顯著降低; 短波單胞菌屬等8個屬在N1處理中含量很低, 但在N6中含量顯著增加, 說明尿素為唯一氮源促進(jìn)以上菌種繁殖。另外,等9個屬在菌源土壤中未檢測到, 而在復(fù)合菌系中檢測出, 說明在繼代培養(yǎng)的條件下新增加了9個微生物屬, 可能來源于空氣, 隨繼代過程中被保留下, 或這9個微生物屬在繼代培養(yǎng)的條件下得到了大量繁殖, 數(shù)量增加, 在16S rDNA擴增測序中被檢測到。菌屬在N1中含量高于0.1%, 具有直接降解木質(zhì)素的功能, 固氮螺菌屬豐度為7.40%,與固氮螺菌屬具有協(xié)同降解作用[36]; 寡養(yǎng)單胞菌()屬于纖溶微生物, 多報道于石油降解菌種, 具有降解烷烴類物質(zhì)的作用[37];可分泌脲酶和糖苷酶, 常以尿素為氮源[38]; Lewin等[39]和Abt等[40]報道的纖維素降解復(fù)合菌系中有屬的; 在牛糞中富含降解纖維素和木質(zhì)素的菌屬[41]; 類芽孢桿菌屬中的等篩選自土壤, 具有較高的纖維素酶活性, 可直接降解纖維素類物質(zhì)[42];和在相關(guān)文獻(xiàn)中暫沒有報道具有秸稈或纖維素降解功能。

    復(fù)合菌系通過調(diào)節(jié)微生物組成, 進(jìn)一步調(diào)節(jié)其功能多樣性, 從而適應(yīng)不同尿素含量培養(yǎng)條件, 進(jìn)而促進(jìn)秸稈的降解。不同氮源處理間, N2處理復(fù)合菌系A(chǔ)lpha多樣性指數(shù)顯著高于其他處理, 菌種組成多樣性較高,根瘤菌等參與秸稈降解途徑的菌屬含量增加, 促進(jìn)其秸稈降解進(jìn)程, 導(dǎo)致N2處理秸稈降解效率較高。

    4 結(jié)論

    菌源樣品及復(fù)合菌系間的Alpha多樣指數(shù)、OTU豐度和差異度、物種結(jié)構(gòu)及含量均存在顯著差異, 不同氮源處理間, N2處理復(fù)合菌系的菌種組成多樣性較高, 菌群結(jié)構(gòu)更加豐富、均勻。在限制性繼代培養(yǎng)過程中保留含量高于0.1%菌屬有纖維弧菌根瘤菌土地桿菌噬氫菌黃桿菌德沃斯氏菌柄桿菌鞘氨醇單胞菌等, 均可分泌水解酶和氧化酶, 直接或間接參與秸稈降解進(jìn)程中; 在氮源硫酸銨和尿素含量0.16+0.04%(N2)配比條件下, 優(yōu)勢菌種為纖維弧菌、黃桿菌、土地桿菌, N2處理較其他處理避免了單一氮源造成的局限性, 多樣性指數(shù)增加, 菌群組成結(jié)構(gòu)優(yōu)化, 碳水化合物代謝功能注釋相對豐度較高, 提高了秸稈降解效率, 為進(jìn)一步開發(fā)秸稈還田促腐菌劑提供基礎(chǔ)。

    References

    [1] MARSCHNER P, KANDELER E, MARSCHNER B. Structure and function of the soil microbial community in a long-term fertilizer experiment[J]. Soil Biology and Biochemistry, 2003, 35(3): 453–461

    [2] PETERSSON L, KVIEN I, OKSMAN K. Structure and thermal properties of poly (lactic acid)/cellulose whiskers nanocomposite materials[J]. Composites Science and Technology, 2007, 67(11/12): 2535–2544

    [3] YANG H Y, WU H, WANG X F, et al. Selection and characteristics of a switchgrass-colonizing microbial community to produce extracellular cellulases and xylanases[J]. Bioresource Technology, 2011, 102(3): 3546–3550

    [4] WANG X J, YUAN X F, WANG H, et al. Characteristics and community diversity of a wheat straw-colonizing microbial community[J]. African Journal of Biotechnology, 2011, 10(40): 7853–7861

    [5] ULRICH A, KLIMKE G, WIRTH S. Diversity and activity of cellulose-decomposing bacteria, isolated from a sandy and a loamy soil after long-term manure application[J]. Microbial Ecology, 2008, 55(3): 512–522

    [6] LI P P, WANG X F, YUAN X F, et al. Screening of a composite microbial system and its characteristics of wheat straw degradation[J]. Agricultural Sciences in China, 2011, 10(10): 1586–1594

    [7] FENG Y J, YU Y L, WANG X, et al. Degradation of raw corn stover powder (RCSP) by an enriched microbial consortium and its community structure[J]. Bioresource Technology, 2011, 102(2): 742–747

    [8] QINGGEER, GAO J L, YU X F, et al. Screening of a microbial consortium with efficient corn stover degradation ability at low temperature[J]. Journal of Integrative Agriculture, 2016, 15(10): 2369–2379

    [9] YU X F, BORJIGIN Q, GAO J L, et al. Exploration of the key microbes and composition stability of microbial consortium GF-20 with efficiently decomposes corn stover at low temperatures[J]. Journal of Integrative Agriculture, 2019, 18(8): 1893–1904

    [10] PARKS D H, BEIKO R G. Identifying biologically relevant differences between metagenomic communities[J]. Bioinformatics, 2010, 26(6): 715–721

    [11] NELSON C E, GARDNER J G. In-frame deletions allow functional characterization of complex cellulose degradation phenotypes in[J]. Applied and Environmental Microbiology, 2015, 81(17): 5968–5975

    [12] MILLWARD-SADLER S J, DAVIDSON K, HAZLEWOOD G P, et al. Novel cellulose-binding domains, NodB homologues and conserved modular architecture in xylanases from the aerobic soil bacteriasubsp.and[J]. Biochemical Journal, 1995, 312(1): 39–48

    [13] WU Y R, HE J Z. Characterization of a xylanase-producingstrain J3-8 and its genome analysis[J]. Scientific Reports, 2015, 5: 10521

    [14] WANG W D, YAN L, CUI Z J, et al. Characterization of a microbialcapable of degrading lignocellulose[J]. Bioresource Technology, 2011, 102(19): 9321–9324

    [15] DA X Y, JIANG F, CHANG X L, et al.sp. nov., isolated from soil in Antarctica[J]. International Journal of Systematic and Evolutionary Microbiology, 2015, 65(11): 3841–3846

    [16] WON K H, KOOK M, YI T H.sp. nov., isolated from soil of a bamboo plantation[J]. Antonie van Leeuwenhoek, 2015, 107(2): 565–573

    [17] KANG J Y, CHUN J, JAHNG K Y.sp. nov., isolated from freshwater, and emended description of the genus[J]. International Journal of Systematic and Evolutionary Microbiology, 2013, 63(5): 1633–1638

    [18] 周麗沙, 李慧, 張穎, 等. 石油污染土壤鞘氨醇單胞菌遺傳多樣性16S rDNA-PCR-DGGE分析[J]. 土壤學(xué)報, 2011, 48(4): 804–812 ZHOU L S, LI H, ZHANG Y, et al. Analysis ofgenetic diversity in petroleum-contaminated soils by using PCR-DGGE technique[J]. Acta Pedologica Sinica, 2011, 48(4): 804–812

    [19] DUAN J, LIANG J D, DU W J, et al. Biodegradation of kraft lignin by a bacterial strainsp. HY-H[J]. Advanced Materials Research, 2014, 955/959: 548–553

    [20] XING W, LI J L, LI P, et al. Effects of residual organics in municipal wastewater on hydrogenotrophic denitrifying microbial communities[J]. Journal of Environmental Sciences, 2018, 65: 262–270

    [21] LA SCOLA B, MALLET M N, GRIMONT P A D, et al.sp. nov.,sp. nov. andsp. nov., isolated from hospital water supplies, and emendation of the genus(Das1996)[J]. International Journal of Systematic and Evolutionary Microbiology, 2003, 53(1): 15–20

    [22] XU L, ZHANG Y, READ N, et al.sp. nov., a novel species isolated from the roots ofin China[J]. Antonie van Leeuwenhoek, 2017, 110(11): 1475–1483

    [23] HASSAN Y I, LEPP D, ZHOU T. Genome assemblies of three soil-associatedspecies:,.,and.[J]. Genome Announcements, 2015, 3(3): e00514-15

    [24] JIN L, LEE H G, KIM H S, et al.sp. nov., a stalked bacterium isolated from a eutrophic reservoir[J]. International Journal of Systematic and Evolutionary Microbiology, 2013, 63(7): 2559–2564

    [25] SUN L N, YANG E D, WEI J C, et al.sp. nov., a stalked bacterium isolated from rhizosphere soil[J]. International Journal of Systematic and Evolutionary Microbiology, 2015, 65(12): 4374–4380

    [26] JACKSON C A, COUGER M B, PRABHAKARAN M, et al. Isolation and characterization ofsp. strain YS-1r that degrades lignin in plant biomass[J]. Journal of Applied Microbiology, 2017, 122(4): 940–952

    [27] LAUS M C, VAN BRUSSEL A A N, KIJNE J W. Role of cellulose fibrils and exopolysaccharides ofin attachment to and infection ofroot hairs[J]. Molecular Plant-Microbe Interactions, 2005, 18(6): 533–538

    [28] TER BEEK J, GUSKOV A, SLOTBOOM D J. Structural diversity of ABC transporters[J]. Journal of General Physiology, 2014, 143(4): 419–435

    [29] UEKI A, KODAMA Y, KAKU N, et al.gen. nov., sp. nov., a facultatively anaerobic, fermentative stalked bacterium in the classisolated from rice plant roots[J]. The Journal of General and Applied Microbiology, 2010, 56(3): 193–203

    [30] KODAMA Y, WATANABE K.sp. nov., a facultatively anaerobic, fermentative, prosthecate bacterium isolated from a cellulose-fed microbial fuel cell[J]. International Journal of Systematic and Evolutionary Microbiology, 2011, 61(8): 1781–1785

    [31] ABRAHAM W R, ESTRELA A B, NIKITIN D I, et al.sp. nov.,sp. nov. andsp. nov., prosthecate bacteria from aquatic habitats[J]. International Journal of Systematic and Evolutionary Microbiology, 2010, 60(8): 1837–1843

    [32] 司美茹, 趙云峰, 蘇濤. 高效降解烷烴的無色桿菌XL株的分離鑒定及其降解特性[J]. 土壤通報, 2011, 42(3): 562–567 SI M R, ZHAO Y F, SU T. Isolation and identification of a high-efficiency alkane-degradingXL strain and its degradation characteristics[J]. Chinese Journal of Soil Science, 2011, 42(3): 562–567

    [33] 唐玉斌, 孫常宇, 陳芳艷, 等. 一株[艸屈]高效降解菌的分離鑒定及其降解特性[J]. 微生物學(xué)通報, 2009, 36(4): 593–597 TANG Y B, SUN C Y, CHEN F Y, et al. Isolation and identification of a chrysene-degrading strain and its degradation characteristics[J]. Microbiology, 2009, 36(4): 593–597

    [34] WAKARCHUK W W, BROCHU D, FOOTE S, et al. Proteomic analysis of the secretome ofATCC 484 andATCC 482[J]. PLoS One, 2016, 11(3): e0151186

    [35] ZHUANG W P, ZHANG S Z, XIA X, et al. Draft genome sequence ofT26Tand comparative analysis of six[J]. Standards in Genomic Sciences, 2015, 10: 104

    [36] KUPRYASHINA M A, PETROV S V, PONOMAREVA E G, et al. Ligninolytic activity of bacteria of the generaand[J]. Microbiology, 2015, 84(6): 791–795

    [37] 肖建軍, 李亞龍, 楊琦. 寡養(yǎng)單胞菌降解石油污染土壤中的甲苯[J]. 環(huán)境工程, 2018, 36(4): 186–189 XIAO J J, LI Y L, YANG Q. Degradation of toluene in petroleum contaminated soil bysp. MJ-1[J]. Environmental Engineering, 2018, 36(4): 186–189

    [38] MADHAIYAN M, JIN T Y, ROY J J, et al.sp. nov., an endophytic N-fixing bacterium isolated from root tissue ofL.[J]. International Journal of Systematic and Evolutionary Microbiology, 2013, 63(7): 2477–2483

    [39] LEWIN G R, JOHNSON A L, SOTO R D M, et al. Cellulose-enriched microbial communities from leaf-cutter ant () refuse dumps vary in taxonomic composition and degradation ability[J]. PLoS One, 2016, 11(3): e0151840

    [40] ABT B, TESHIMA H, LUCAS S, et al. Complete genome sequence oftype strain (4M15T)[J]. Standards in Genomic Sciences, 2011, 4(1): 2–12

    [41] ZIEMER C J. Newly cultured bacteria with broad diversity isolated from eight-week continuous culture enrichments of cow feces on complex polysaccharides[J]. Applied and Environmental Microbiology, 2014, 80(2): 574–585

    [42] LIANG Y L, ZHANG Z, WU M, et al. Isolation, screening, and identification of cellulolytic bacteria from natural reserves in the subtropical region of China and optimization of cellulase production byME27–1[J]. Biomed Research International, 2014, 2014: 512497

    Straw degradation ability and composition of microbial consortium for corn straw decomposition at low temperature*

    Qinggeer1,2, YU Xiaofang1,2**, GAO Julin1,2**, WANG Zhigang1,2, HU Wanji1, Naoganchaolu2,3, WANG Zhen2,3, HU Shuping2,4, SUN Jiying1,2, QU Jiawei1,2

    (1. Agricultural College, Inner Mongolia Agricultural University, Hohhot 010019, China; 2. Key Laboratory of Crop Cultivation and Genetic Improvement in Inner Mongolia Autonomous Region, Hohhot 010019, China; 3. Horticulture and Plant Protection College, Inner Mongolia Agricultural University, Hohhot 010019, China; 4. Vocational and Technical College, Inner Mongolia Agricultural University, Baotou 014100, China)

    It is common practice to return field straw with urea to accelerate decomposition. To improve the microbial screening methodology and investigate the microbes responsible for decomposition of corn straw, nitrogen acclimatization of the microbial consortium GF-20 (GF-20) was performed. Compositional differences in the cultured microbial consortium and the microbe source were also evaluated. GF-20 was cultured until the 10thgeneration in variable nitrogen conditions [ammonium sulfate (N1), mixtures of ammonium sulfate and urea (N2-N5), and urea (N6)]. The corn straw decomposition ratio was determined to estimate the activity of the composite microbial system, and the composition diversity and function were analyzed by MiSeq high-throughput sequencing. The results showed that the N2 degradation rate was significantly higher than the other treatments. The bacterial source alpha diversity index (ADI) was significantly higher than the cultured microbial consortium, and the N2 ADI was significantly higher than the other nitrogen treatments. The bacterial composition also significantly differed between the source and consortium, as well as among the nitrogen treatments. The N2 treatment yielded the most diverse bacterial composition, with richer and more uniform flora structures and a higher carbohydrate metabolic pathway activity (which promotes corn straw degradation). Functional microbial strains involved in corn straw degradation were obtained after restrictive sub-generation of the microbial sources, which can accelerate corn straw degradation. The highest corn straw degradation efficiency of microbial consortium was observed with the 0.16% ammonium sulfate + 0.04% of urea nitrogen (N2) treatment. These findings provide a basis for developing microbial consortium used in commercial production decomposition.

    Original microbe source; Microbial consortium; Nitrogen resources; High throughput sequencing; Composition diversity; Corn straw degradation

    s: YU Xiaofang, E-mail: yuxiaofang75@163.com; GAO Julin, E-mail: nmgaojulin@163.com

    S145.6

    10.13930/j.cnki.cjea.200128

    Feb. 26, 2020;

    Jun. 21, 2020

    青格爾,于曉芳, 高聚林, 王志剛, 胡萬吉, 鬧干朝魯, 王振, 胡樹平, 孫繼穎, 屈佳偉. 玉米秸稈低溫降解復(fù)合菌系降解能力及微生物組成研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(中英文), 2020, 28(11): 1753-1765

    Qinggeer, YU X F, GAO J L, WANG Z G, HU W J, Naoganchaolu, WANG Z, HU S P, SUN J Y, QU J W. Straw degradation ability and composition of microbial consortium for corn straw decomposition at low temperature[J]. Chinese Journal of Eco-Agriculture, 2020, 28(11): 1753-1765

    * 國家自然科學(xué)基金項目(31760353)、內(nèi)蒙古自然基金項目(2020MS03086, 2018ZD02)、國家重點研發(fā)計劃項目(2017YFD0300804)、國家玉米產(chǎn)業(yè)技術(shù)體系項目(CARS-02-63)、農(nóng)業(yè)部華北黃土高原地區(qū)作物栽培科學(xué)觀測實驗站(25204120)和內(nèi)蒙古農(nóng)業(yè)大學(xué)高層次人才引進(jìn)科研啟動項目(NDYB2016-15)資助

    于曉芳, 主要從事作物生理生態(tài)研究, E-mail: yuxiaofang75@163.com; 高聚林, 主要從事作物生理生態(tài)研究, E-mail: nmgaojulin@163.com

    青格爾, 主要從事玉米生理生態(tài)研究。E-mail: qinggeer001@163.com

    2020-02-26

    2020-06-21

    * This study was supported by the National Natural Science Foundation of China (31760353), the Inner Mongolia Natural Sciences Foundation (2020MS03086, 2018ZD02), the National Key Research and Development Project of China (2017YFD0300804), the Earmarked Fund for China Agriculture Research System (CARS-02-63), the Crop Science Observation & Experiment Station in Loess Plateau of North China, Ministry of Agriculture, China (25204120) and the Advanced Talented Scholars of Inner Mongolia Agricultural University, China (NDYB2016-15).

    猜你喜歡
    菌系菌源桿菌屬
    鄢陵縣2023年玉米南方銹病發(fā)生與防控研究
    湖北植保(2024年4期)2024-01-01 00:00:00
    菌源同工酶分析揭示菌源DPP4同工酶是潛在的抗2型糖尿病靶點
    遺傳(2023年8期)2023-08-25 06:59:00
    人工組合復(fù)合菌系分解木質(zhì)纖維素特征研究
    潰瘍性結(jié)腸炎患者腸道菌群分布特征分析
    玉米秸稈低溫降解復(fù)合菌的篩選及其菌種組成
    利用高通量測序技術(shù)對水稻秸稈中、低溫降解菌系的比較分析
    養(yǎng)豬微生物發(fā)酵床芽胞桿菌空間生態(tài)位特性
    來安縣2018年小麥赤霉病發(fā)生特點與防效分析
    香辛料對哈爾濱風(fēng)干腸菌系及理化性質(zhì)的影響
    肉類研究(2017年12期)2018-01-22 17:48:10
    小麥赤霉病發(fā)生與菌源?氣象因子量關(guān)系分析
    别揉我奶头~嗯~啊~动态视频| 日本在线视频免费播放| 亚洲欧美日韩东京热| 男人舔奶头视频| 国产男人的电影天堂91| 男人和女人高潮做爰伦理| 久久久久精品国产欧美久久久| 国产伦一二天堂av在线观看| 91av网一区二区| 综合色av麻豆| 中文字幕精品亚洲无线码一区| 大又大粗又爽又黄少妇毛片口| 三级男女做爰猛烈吃奶摸视频| 亚洲美女黄片视频| 欧美性猛交╳xxx乱大交人| 久久午夜亚洲精品久久| 成人一区二区视频在线观看| 亚洲精品亚洲一区二区| 露出奶头的视频| 亚洲四区av| 中文在线观看免费www的网站| 色综合色国产| 国产真实伦视频高清在线观看| 97在线视频观看| 别揉我奶头~嗯~啊~动态视频| 欧美一级a爱片免费观看看| 欧美区成人在线视频| 婷婷六月久久综合丁香| 欧美成人一区二区免费高清观看| 老师上课跳d突然被开到最大视频| 少妇高潮的动态图| 国产蜜桃级精品一区二区三区| 欧美zozozo另类| 亚洲精品456在线播放app| 国产伦在线观看视频一区| 日产精品乱码卡一卡2卡三| 丰满乱子伦码专区| 国产视频内射| 久久久久国产网址| 色综合亚洲欧美另类图片| 男女视频在线观看网站免费| 国产高清三级在线| 天天躁夜夜躁狠狠久久av| 国产国拍精品亚洲av在线观看| 午夜精品在线福利| 小蜜桃在线观看免费完整版高清| 亚洲自拍偷在线| 午夜亚洲福利在线播放| 国产乱人偷精品视频| 黄色日韩在线| 国产一区亚洲一区在线观看| 国产精品美女特级片免费视频播放器| 真人做人爱边吃奶动态| 蜜臀久久99精品久久宅男| 国产一区二区三区在线臀色熟女| 中出人妻视频一区二区| 看十八女毛片水多多多| 男女做爰动态图高潮gif福利片| 91久久精品国产一区二区成人| 欧美极品一区二区三区四区| 插逼视频在线观看| 深夜精品福利| 简卡轻食公司| 亚洲人成网站在线观看播放| 人人妻人人澡欧美一区二区| 成人午夜高清在线视频| 99久久成人亚洲精品观看| 亚洲成人av在线免费| 亚洲综合色惰| 国产成人影院久久av| 国产综合懂色| 亚洲激情五月婷婷啪啪| 在线观看美女被高潮喷水网站| 女人被狂操c到高潮| 欧美+日韩+精品| av天堂在线播放| 精品久久久久久久久av| 国产69精品久久久久777片| 亚洲无线观看免费| 国内揄拍国产精品人妻在线| 免费观看的影片在线观看| 老女人水多毛片| 亚洲第一区二区三区不卡| 国产欧美日韩精品亚洲av| av女优亚洲男人天堂| 美女被艹到高潮喷水动态| 老熟妇仑乱视频hdxx| 欧美一区二区国产精品久久精品| 蜜桃亚洲精品一区二区三区| 欧美区成人在线视频| 日韩三级伦理在线观看| 97超级碰碰碰精品色视频在线观看| 国产成人一区二区在线| 国产精品一及| 乱系列少妇在线播放| 赤兔流量卡办理| 亚洲真实伦在线观看| 18禁在线播放成人免费| 国产av一区在线观看免费| 成年女人看的毛片在线观看| 露出奶头的视频| 国产69精品久久久久777片| 日韩国内少妇激情av| 少妇丰满av| videossex国产| av在线老鸭窝| 久久精品91蜜桃| 国产乱人视频| 看黄色毛片网站| 亚洲最大成人中文| 国产熟女欧美一区二区| ponron亚洲| 日本色播在线视频| 搡老熟女国产l中国老女人| 麻豆一二三区av精品| 黄色配什么色好看| 日韩人妻高清精品专区| 日韩精品有码人妻一区| 18禁在线播放成人免费| 成人av在线播放网站| 亚洲中文日韩欧美视频| 五月伊人婷婷丁香| 亚洲国产精品成人久久小说 | 国产在线精品亚洲第一网站| 97碰自拍视频| 中国美女看黄片| 国产在线精品亚洲第一网站| 欧美日韩在线观看h| 亚洲国产精品成人综合色| 精品久久久久久久久久久久久| 亚洲电影在线观看av| 国国产精品蜜臀av免费| 中出人妻视频一区二区| 国产精品一区二区性色av| 亚洲精品国产成人久久av| 老司机午夜福利在线观看视频| 精品福利观看| 久久久久久国产a免费观看| a级毛色黄片| 国产毛片a区久久久久| 99久久精品热视频| 欧美zozozo另类| 亚洲欧美日韩高清专用| 成年女人毛片免费观看观看9| 99精品在免费线老司机午夜| 一个人看的www免费观看视频| 国内少妇人妻偷人精品xxx网站| 美女免费视频网站| 久久精品综合一区二区三区| 禁无遮挡网站| 蜜桃久久精品国产亚洲av| 日韩,欧美,国产一区二区三区 | 在线a可以看的网站| 久久精品夜夜夜夜夜久久蜜豆| 最后的刺客免费高清国语| 人妻丰满熟妇av一区二区三区| 亚洲激情五月婷婷啪啪| 18禁裸乳无遮挡免费网站照片| 亚洲三级黄色毛片| 午夜a级毛片| 99久国产av精品国产电影| 午夜日韩欧美国产| 精品乱码久久久久久99久播| 啦啦啦啦在线视频资源| 免费在线观看影片大全网站| 亚洲aⅴ乱码一区二区在线播放| 久久人妻av系列| 97碰自拍视频| 国产成人a∨麻豆精品| 啦啦啦观看免费观看视频高清| 国产男靠女视频免费网站| 69av精品久久久久久| 老女人水多毛片| 久久久午夜欧美精品| 亚洲av熟女| av在线老鸭窝| 欧美bdsm另类| 欧美xxxx性猛交bbbb| 日本熟妇午夜| 欧美一区二区国产精品久久精品| 国产69精品久久久久777片| 久久人妻av系列| 国产精品精品国产色婷婷| 久久久久九九精品影院| 国产v大片淫在线免费观看| 搡老熟女国产l中国老女人| 99国产极品粉嫩在线观看| 丰满乱子伦码专区| 国产成人aa在线观看| 中文字幕久久专区| 免费观看在线日韩| 久久中文看片网| 久久国内精品自在自线图片| 中文在线观看免费www的网站| 精品不卡国产一区二区三区| 九九久久精品国产亚洲av麻豆| 日日摸夜夜添夜夜爱| 99久国产av精品| 99久久九九国产精品国产免费| 五月玫瑰六月丁香| 91精品国产九色| 亚洲第一区二区三区不卡| 一级毛片久久久久久久久女| 99精品在免费线老司机午夜| 男女下面进入的视频免费午夜| 国产人妻一区二区三区在| 少妇熟女aⅴ在线视频| 国产一区亚洲一区在线观看| 亚洲精品456在线播放app| 亚洲第一区二区三区不卡| 免费大片18禁| 国产激情偷乱视频一区二区| 99久久精品热视频| 国产v大片淫在线免费观看| 国内精品美女久久久久久| 69人妻影院| 亚洲欧美日韩高清在线视频| 搡女人真爽免费视频火全软件 | 免费av观看视频| 99热网站在线观看| 欧美成人免费av一区二区三区| 精品一区二区三区人妻视频| videossex国产| 中文资源天堂在线| 少妇的逼好多水| 国产精品久久电影中文字幕| 欧美日本亚洲视频在线播放| 好男人在线观看高清免费视频| 天天一区二区日本电影三级| 亚洲性夜色夜夜综合| 国产aⅴ精品一区二区三区波| 国产男靠女视频免费网站| 日本一二三区视频观看| 欧美绝顶高潮抽搐喷水| 日韩强制内射视频| 精品一区二区三区av网在线观看| www.色视频.com| 国产日本99.免费观看| 免费观看在线日韩| 午夜福利在线观看吧| 国产精品不卡视频一区二区| 国产一区亚洲一区在线观看| 又爽又黄a免费视频| 国产午夜精品久久久久久一区二区三区 | 国产成人aa在线观看| 亚洲精品一卡2卡三卡4卡5卡| 丰满乱子伦码专区| 中文字幕久久专区| 精品乱码久久久久久99久播| 久久久久久久久久久丰满| 亚洲精品在线观看二区| 日韩亚洲欧美综合| 人妻夜夜爽99麻豆av| а√天堂www在线а√下载| 亚洲三级黄色毛片| 免费不卡的大黄色大毛片视频在线观看 | 亚洲中文字幕日韩| 久久天躁狠狠躁夜夜2o2o| 国产伦一二天堂av在线观看| or卡值多少钱| 亚洲av不卡在线观看| 日本免费一区二区三区高清不卡| 久久国内精品自在自线图片| 老熟妇乱子伦视频在线观看| 男女之事视频高清在线观看| 亚洲性久久影院| 床上黄色一级片| 中文字幕人妻熟人妻熟丝袜美| 久久久午夜欧美精品| 亚洲aⅴ乱码一区二区在线播放| 亚洲天堂国产精品一区在线| 欧美人与善性xxx| 神马国产精品三级电影在线观看| 听说在线观看完整版免费高清| 干丝袜人妻中文字幕| 国产成人一区二区在线| 午夜a级毛片| 给我免费播放毛片高清在线观看| a级一级毛片免费在线观看| 久久人人精品亚洲av| 少妇猛男粗大的猛烈进出视频 | 干丝袜人妻中文字幕| 在线观看免费视频日本深夜| 美女内射精品一级片tv| 欧美不卡视频在线免费观看| av视频在线观看入口| 午夜免费男女啪啪视频观看 | 亚洲精华国产精华液的使用体验 | 在线免费十八禁| .国产精品久久| 日韩欧美免费精品| 日日摸夜夜添夜夜添av毛片| 欧美性猛交黑人性爽| av在线老鸭窝| 亚洲一级一片aⅴ在线观看| 99热这里只有精品一区| 黄片wwwwww| 麻豆乱淫一区二区| 亚洲精品日韩在线中文字幕 | 成人综合一区亚洲| 国产精品一区二区三区四区久久| 国产av一区在线观看免费| 日韩中字成人| 精品午夜福利视频在线观看一区| 搡老岳熟女国产| 秋霞在线观看毛片| 嫩草影视91久久| 在现免费观看毛片| 亚洲av五月六月丁香网| 天堂√8在线中文| 国产精品一二三区在线看| 日韩精品中文字幕看吧| 高清日韩中文字幕在线| 麻豆成人午夜福利视频| 18+在线观看网站| 看非洲黑人一级黄片| 日韩高清综合在线| 精品久久久久久久人妻蜜臀av| 亚洲乱码一区二区免费版| 亚洲欧美日韩卡通动漫| 观看免费一级毛片| 日韩精品青青久久久久久| 99国产精品一区二区蜜桃av| 乱码一卡2卡4卡精品| 又粗又爽又猛毛片免费看| 岛国在线免费视频观看| 十八禁国产超污无遮挡网站| 国产三级中文精品| 亚洲av熟女| 啦啦啦观看免费观看视频高清| 国产精品无大码| 国产精华一区二区三区| 日韩在线高清观看一区二区三区| 亚洲欧美成人综合另类久久久 | 一区二区三区免费毛片| 在线a可以看的网站| 夜夜爽天天搞| 特级一级黄色大片| 午夜亚洲福利在线播放| 少妇熟女aⅴ在线视频| 在线免费观看的www视频| 欧美xxxx性猛交bbbb| 国产高清视频在线播放一区| 97热精品久久久久久| 看片在线看免费视频| 欧美精品国产亚洲| 久久久久久久亚洲中文字幕| 国产欧美日韩精品亚洲av| 99热这里只有是精品50| 91av网一区二区| 久久久国产成人免费| 美女被艹到高潮喷水动态| 特大巨黑吊av在线直播| 欧美三级亚洲精品| 乱人视频在线观看| 亚洲国产精品国产精品| 人人妻人人澡人人爽人人夜夜 | 我要看日韩黄色一级片| 日韩成人av中文字幕在线观看 | 六月丁香七月| 禁无遮挡网站| 免费一级毛片在线播放高清视频| 亚洲国产精品sss在线观看| 美女cb高潮喷水在线观看| 欧美绝顶高潮抽搐喷水| 十八禁网站免费在线| 99久国产av精品| 91久久精品电影网| 丝袜美腿在线中文| 亚洲丝袜综合中文字幕| 久久亚洲精品不卡| 51国产日韩欧美| 欧美最新免费一区二区三区| 色综合站精品国产| 国产乱人偷精品视频| 国产精品福利在线免费观看| 久久九九热精品免费| 国产成人freesex在线 | 校园春色视频在线观看| 99久久精品热视频| 又黄又爽又刺激的免费视频.| 久久中文看片网| 最近最新中文字幕大全电影3| 97超视频在线观看视频| 国产色婷婷99| 亚洲熟妇熟女久久| 看黄色毛片网站| 我要看日韩黄色一级片| 亚洲最大成人中文| 中国美女看黄片| 全区人妻精品视频| 在线免费十八禁| 大型黄色视频在线免费观看| 日本黄色片子视频| 最新中文字幕久久久久| 中文在线观看免费www的网站| 亚洲一级一片aⅴ在线观看| 国产在线精品亚洲第一网站| 成年女人永久免费观看视频| 国产亚洲精品综合一区在线观看| 日本免费a在线| 色在线成人网| 美女高潮的动态| 18禁黄网站禁片免费观看直播| 成熟少妇高潮喷水视频| 一进一出好大好爽视频| 精品一区二区三区av网在线观看| 久久精品综合一区二区三区| 国产精品亚洲美女久久久| 69人妻影院| 免费无遮挡裸体视频| 中国美女看黄片| 亚洲电影在线观看av| 麻豆国产av国片精品| 国产乱人视频| 一级av片app| 午夜a级毛片| 级片在线观看| 精品久久久久久久久av| 中文资源天堂在线| 噜噜噜噜噜久久久久久91| 欧美激情国产日韩精品一区| 嫩草影院入口| 色哟哟哟哟哟哟| 51国产日韩欧美| 欧美日本亚洲视频在线播放| 亚洲av不卡在线观看| 久久99热这里只有精品18| 观看免费一级毛片| 欧美一区二区国产精品久久精品| 欧美绝顶高潮抽搐喷水| 一区二区三区免费毛片| 国产免费一级a男人的天堂| 三级毛片av免费| 无遮挡黄片免费观看| 午夜精品国产一区二区电影 | 女生性感内裤真人,穿戴方法视频| 日韩在线高清观看一区二区三区| 女人十人毛片免费观看3o分钟| av天堂中文字幕网| 亚洲国产精品成人久久小说 | 偷拍熟女少妇极品色| 欧美激情久久久久久爽电影| 午夜免费激情av| 老熟妇仑乱视频hdxx| 国产精品一二三区在线看| 黄片wwwwww| 老师上课跳d突然被开到最大视频| 99在线人妻在线中文字幕| 在线播放无遮挡| 成人午夜高清在线视频| 亚洲第一区二区三区不卡| 最近2019中文字幕mv第一页| 精品人妻视频免费看| 国产精品永久免费网站| 五月玫瑰六月丁香| 国产亚洲精品av在线| 丝袜美腿在线中文| 亚洲欧美精品自产自拍| 在线免费观看不下载黄p国产| 成人三级黄色视频| 国产又黄又爽又无遮挡在线| 欧美一区二区国产精品久久精品| 国内精品一区二区在线观看| 给我免费播放毛片高清在线观看| 男人舔女人下体高潮全视频| 久久精品夜色国产| 久久人人精品亚洲av| 日本三级黄在线观看| 国产亚洲91精品色在线| 欧美日韩国产亚洲二区| 18禁在线无遮挡免费观看视频 | 夜夜夜夜夜久久久久| 少妇丰满av| 成人美女网站在线观看视频| 久久综合国产亚洲精品| 女同久久另类99精品国产91| 狂野欧美白嫩少妇大欣赏| 97人妻精品一区二区三区麻豆| 久久精品国产鲁丝片午夜精品| 亚洲自偷自拍三级| 国产精品亚洲美女久久久| 久久午夜福利片| 一进一出抽搐gif免费好疼| 亚洲一区高清亚洲精品| 丰满人妻一区二区三区视频av| 亚洲中文日韩欧美视频| 一个人观看的视频www高清免费观看| 国产成人精品久久久久久| 99热6这里只有精品| 久久精品国产自在天天线| 99热这里只有是精品50| 一级毛片久久久久久久久女| 高清午夜精品一区二区三区 | 中文资源天堂在线| 成人无遮挡网站| 91久久精品电影网| 午夜福利成人在线免费观看| 国产欧美日韩精品亚洲av| 伦精品一区二区三区| 日韩高清综合在线| 久久久午夜欧美精品| 高清毛片免费看| 国产成人福利小说| 又黄又爽又免费观看的视频| 1024手机看黄色片| 日本色播在线视频| 人妻丰满熟妇av一区二区三区| 欧美一区二区亚洲| 成人特级黄色片久久久久久久| 99热6这里只有精品| 男女做爰动态图高潮gif福利片| 美女被艹到高潮喷水动态| av视频在线观看入口| 麻豆久久精品国产亚洲av| 亚洲精品456在线播放app| 亚洲一级一片aⅴ在线观看| 国产精品不卡视频一区二区| 久久久久国内视频| 国产av在哪里看| 中文字幕人妻熟人妻熟丝袜美| or卡值多少钱| 欧美区成人在线视频| 久久精品国产亚洲网站| 亚洲综合色惰| 中国美白少妇内射xxxbb| 午夜老司机福利剧场| 人人妻,人人澡人人爽秒播| 亚洲精品乱码久久久v下载方式| 国产一区二区亚洲精品在线观看| 国产毛片a区久久久久| 别揉我奶头~嗯~啊~动态视频| 3wmmmm亚洲av在线观看| 99热精品在线国产| 卡戴珊不雅视频在线播放| 国产爱豆传媒在线观看| 亚洲久久久久久中文字幕| 女生性感内裤真人,穿戴方法视频| 婷婷亚洲欧美| 人妻少妇偷人精品九色| 国产精品av视频在线免费观看| 成人精品一区二区免费| 97热精品久久久久久| 日产精品乱码卡一卡2卡三| 日韩欧美国产在线观看| 嫩草影院精品99| 中文亚洲av片在线观看爽| 欧美丝袜亚洲另类| 九九久久精品国产亚洲av麻豆| 看黄色毛片网站| 欧美成人一区二区免费高清观看| 亚洲真实伦在线观看| 久久人妻av系列| 日本免费a在线| 我的老师免费观看完整版| 午夜福利成人在线免费观看| 国国产精品蜜臀av免费| 亚洲成人av在线免费| 桃色一区二区三区在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲18禁久久av| 人妻制服诱惑在线中文字幕| 国产真实乱freesex| eeuss影院久久| 亚洲av中文字字幕乱码综合| 国产精品久久视频播放| 国产三级在线视频| 国产精品一区二区免费欧美| 嫩草影院新地址| 久久久久久九九精品二区国产| 日本在线视频免费播放| 国产精品精品国产色婷婷| 亚洲电影在线观看av| 特大巨黑吊av在线直播| 色综合站精品国产| 中文亚洲av片在线观看爽| 亚洲精品456在线播放app| 国产成人福利小说| 三级男女做爰猛烈吃奶摸视频| 最新在线观看一区二区三区| 亚洲自偷自拍三级| 欧美激情久久久久久爽电影| 少妇人妻精品综合一区二区 | 国产精品精品国产色婷婷| 永久网站在线| 亚洲久久久久久中文字幕| 噜噜噜噜噜久久久久久91| 国产午夜福利久久久久久| 校园春色视频在线观看| or卡值多少钱| 国产成人精品久久久久久| 欧美高清性xxxxhd video| 嫩草影院精品99| 色哟哟哟哟哟哟| 亚洲图色成人| 深爱激情五月婷婷| 两个人的视频大全免费| 国产av一区在线观看免费| 我要搜黄色片| 国产爱豆传媒在线观看| 午夜久久久久精精品| а√天堂www在线а√下载| 欧美+亚洲+日韩+国产| 又黄又爽又刺激的免费视频.| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品在线观看二区| 亚洲欧美日韩高清在线视频| 亚洲精品影视一区二区三区av| 欧美日韩在线观看h| 九九热线精品视视频播放| 看十八女毛片水多多多| 麻豆av噜噜一区二区三区| 色哟哟·www| 亚洲国产高清在线一区二区三| 色av中文字幕|