• <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)系分析
    美女国产视频在线观看| 国产午夜精品久久久久久一区二区三区| 国产成人aa在线观看| 久久久久久久久大av| 日本黄大片高清| 国产v大片淫在线免费观看| 高清不卡的av网站| 我的女老师完整版在线观看| 国产免费福利视频在线观看| 久久久精品免费免费高清| 久久人人爽人人片av| 亚洲国产欧美人成| 一区二区三区精品91| 人妻一区二区av| 性高湖久久久久久久久免费观看| av在线播放精品| 一级毛片黄色毛片免费观看视频| 男女国产视频网站| 成人黄色视频免费在线看| 亚洲欧美成人精品一区二区| 色视频www国产| 日本一二三区视频观看| 成人影院久久| 国产伦精品一区二区三区视频9| 午夜福利网站1000一区二区三区| 国产精品99久久久久久久久| 亚洲精品一区蜜桃| 97精品久久久久久久久久精品| 国产精品久久久久久av不卡| 久久久国产一区二区| 亚洲精品久久久久久婷婷小说| 国产探花极品一区二区| 美女中出高潮动态图| 国产91av在线免费观看| 乱系列少妇在线播放| 久久 成人 亚洲| 一本色道久久久久久精品综合| 国产黄色免费在线视频| 夫妻性生交免费视频一级片| 日韩在线高清观看一区二区三区| 国产一级毛片在线| 最近手机中文字幕大全| 亚洲精品国产色婷婷电影| 国产精品秋霞免费鲁丝片| 日韩av在线免费看完整版不卡| a级毛片免费高清观看在线播放| 蜜臀久久99精品久久宅男| 性高湖久久久久久久久免费观看| 男女无遮挡免费网站观看| 中文字幕免费在线视频6| 亚洲欧美中文字幕日韩二区| 综合色丁香网| 一本—道久久a久久精品蜜桃钙片| 亚洲精品日韩在线中文字幕| 久久鲁丝午夜福利片| 成人午夜精彩视频在线观看| 91精品国产九色| 亚洲国产精品成人久久小说| 在线观看免费高清a一片| 精品一区在线观看国产| 久久久久精品久久久久真实原创| 午夜福利在线观看免费完整高清在| 欧美人与善性xxx| 国产色爽女视频免费观看| 免费看日本二区| 国产精品一区www在线观看| 久久青草综合色| 国产视频首页在线观看| 99久久精品热视频| 久久精品人妻少妇| 在线看a的网站| 黄色日韩在线| 91精品一卡2卡3卡4卡| 国产乱人偷精品视频| 久久精品国产亚洲网站| 嫩草影院入口| 亚洲精品国产成人久久av| 精品一区二区免费观看| 欧美日韩视频精品一区| 欧美精品一区二区大全| 伦理电影大哥的女人| 国产精品人妻久久久久久| 高清在线视频一区二区三区| 麻豆乱淫一区二区| 高清欧美精品videossex| 男女国产视频网站| 日本vs欧美在线观看视频 | 亚洲不卡免费看| 老司机影院成人| 少妇裸体淫交视频免费看高清| 精品久久久久久久久亚洲| 男女边摸边吃奶| 永久网站在线| 青春草视频在线免费观看| 极品教师在线视频| 国产黄片视频在线免费观看| 婷婷色综合大香蕉| 国产免费福利视频在线观看| 国产精品人妻久久久久久| 纯流量卡能插随身wifi吗| 一级黄片播放器| 国产成人免费无遮挡视频| 久久久久久久久大av| 黄片无遮挡物在线观看| 亚洲熟女精品中文字幕| 自拍偷自拍亚洲精品老妇| 伦理电影大哥的女人| 久久这里有精品视频免费| 国产久久久一区二区三区| 狂野欧美白嫩少妇大欣赏| 好男人视频免费观看在线| 91久久精品国产一区二区三区| 成人国产麻豆网| 精品人妻偷拍中文字幕| 久久鲁丝午夜福利片| 秋霞伦理黄片| 中文字幕人妻熟人妻熟丝袜美| 欧美丝袜亚洲另类| 日本欧美国产在线视频| 日韩成人av中文字幕在线观看| 成人18禁高潮啪啪吃奶动态图 | 一本—道久久a久久精品蜜桃钙片| 国产精品福利在线免费观看| 亚洲国产高清在线一区二区三| 亚洲无线观看免费| 成人综合一区亚洲| 亚洲人成网站在线播| 麻豆乱淫一区二区| 亚洲国产精品成人久久小说| 午夜福利影视在线免费观看| 综合色丁香网| 国产一区亚洲一区在线观看| 亚洲欧美日韩东京热| 久久久久久人妻| 秋霞在线观看毛片| 久久国内精品自在自线图片| 少妇裸体淫交视频免费看高清| 亚洲精品国产色婷婷电影| 国产国拍精品亚洲av在线观看| av播播在线观看一区| 永久网站在线| 亚洲最大成人中文| 亚洲久久久国产精品| kizo精华| 网址你懂的国产日韩在线| 免费看光身美女| 色5月婷婷丁香| 亚洲丝袜综合中文字幕| 国产日韩欧美在线精品| 精品亚洲成国产av| 中文字幕免费在线视频6| 国产人妻一区二区三区在| 激情 狠狠 欧美| 欧美精品国产亚洲| 亚洲av福利一区| 高清日韩中文字幕在线| 国内揄拍国产精品人妻在线| 观看av在线不卡| 亚洲一区二区三区欧美精品| 久久久久久久久久人人人人人人| 亚洲欧美日韩东京热| 午夜免费观看性视频| 国产乱人偷精品视频| 蜜桃久久精品国产亚洲av| 亚洲av在线观看美女高潮| 亚洲av综合色区一区| 少妇人妻精品综合一区二区| 热99国产精品久久久久久7| 国产成人a∨麻豆精品| 中文精品一卡2卡3卡4更新| 五月伊人婷婷丁香| 国产精品国产三级专区第一集| 亚洲av日韩在线播放| 成年美女黄网站色视频大全免费 | 少妇丰满av| 黄片无遮挡物在线观看| 成人黄色视频免费在线看| av免费观看日本| 免费人妻精品一区二区三区视频| 久久99精品国语久久久| 国产乱人视频| 国产熟女欧美一区二区| 麻豆成人av视频| 中文字幕人妻熟人妻熟丝袜美| 免费黄频网站在线观看国产| 成人综合一区亚洲| 80岁老熟妇乱子伦牲交| 丝袜喷水一区| 成人美女网站在线观看视频| 国产成人精品一,二区| 少妇精品久久久久久久| 联通29元200g的流量卡| 国产成人午夜福利电影在线观看| 中文资源天堂在线| 91狼人影院| 亚洲精品乱码久久久久久按摩| 日本wwww免费看| 少妇被粗大猛烈的视频| 在线观看一区二区三区| 国产成人精品婷婷| 18禁在线播放成人免费| 久久久久久久精品精品| 久久 成人 亚洲| 亚洲欧美清纯卡通| 国产片特级美女逼逼视频| 麻豆乱淫一区二区| 久久久久久久大尺度免费视频| 蜜桃久久精品国产亚洲av| 美女中出高潮动态图| 中国美白少妇内射xxxbb| 国产欧美亚洲国产| 九九在线视频观看精品| 男的添女的下面高潮视频| 日韩av不卡免费在线播放| 3wmmmm亚洲av在线观看| xxx大片免费视频| 91久久精品国产一区二区三区| 天天躁日日操中文字幕| 成年人午夜在线观看视频| 建设人人有责人人尽责人人享有的 | 丰满人妻一区二区三区视频av| 啦啦啦在线观看免费高清www| 精品一区在线观看国产| 一级二级三级毛片免费看| av国产精品久久久久影院| 一本一本综合久久| 亚洲成人一二三区av| 内地一区二区视频在线| 高清av免费在线| 日日撸夜夜添| 国产精品一区二区性色av| 中文乱码字字幕精品一区二区三区| 热re99久久精品国产66热6| 亚洲av不卡在线观看| 我的女老师完整版在线观看| 七月丁香在线播放| 97在线视频观看| 久久6这里有精品| 夜夜骑夜夜射夜夜干| a级毛色黄片| 在线观看人妻少妇| 日韩成人伦理影院| 黄色日韩在线| 午夜福利高清视频| 亚洲真实伦在线观看| 99热网站在线观看| 多毛熟女@视频| 亚洲,一卡二卡三卡| 亚洲va在线va天堂va国产| 久久久久精品性色| 男人爽女人下面视频在线观看| 麻豆精品久久久久久蜜桃| 亚洲av在线观看美女高潮| 久久这里有精品视频免费| 久久亚洲国产成人精品v| 女人久久www免费人成看片| 岛国毛片在线播放| 久久久久久久久久久免费av| 永久网站在线| 中文字幕av成人在线电影| 国产一区二区在线观看日韩| 久久久久久久久大av| 777米奇影视久久| 成人18禁高潮啪啪吃奶动态图 | 一级毛片我不卡| 夜夜骑夜夜射夜夜干| 国产美女午夜福利| 久久久久久人妻| 免费看不卡的av| h日本视频在线播放| 男女啪啪激烈高潮av片| 国产精品秋霞免费鲁丝片| 高清av免费在线| av卡一久久| 亚洲欧美一区二区三区国产| 在线亚洲精品国产二区图片欧美 | 男女下面进入的视频免费午夜| 亚洲内射少妇av| 国产免费视频播放在线视频| 日本欧美国产在线视频| 极品教师在线视频| 97精品久久久久久久久久精品| 成人二区视频| 一级毛片aaaaaa免费看小| 色综合色国产| 国产永久视频网站| 日本vs欧美在线观看视频 | 久久精品久久久久久噜噜老黄| 中文精品一卡2卡3卡4更新| 国产成人免费无遮挡视频| av女优亚洲男人天堂| 成年av动漫网址| 日韩欧美一区视频在线观看 | 免费大片黄手机在线观看| 男女免费视频国产| 久久精品人妻少妇| av国产免费在线观看| 亚洲av福利一区| 亚洲婷婷狠狠爱综合网| 久久久久国产精品人妻一区二区| videos熟女内射| 国产在线视频一区二区| 欧美精品人与动牲交sv欧美| 国产精品av视频在线免费观看| 又粗又硬又长又爽又黄的视频| 91精品国产国语对白视频| 久久精品国产亚洲网站| av在线老鸭窝| 性高湖久久久久久久久免费观看| 黄片无遮挡物在线观看| 七月丁香在线播放| av在线播放精品| 亚洲精品第二区| 男人添女人高潮全过程视频| 色综合色国产| 成人毛片60女人毛片免费| 欧美日韩精品成人综合77777| 成人影院久久| 国产成人aa在线观看| 成人美女网站在线观看视频| 人人妻人人澡人人爽人人夜夜| 免费观看a级毛片全部| 日韩强制内射视频| 久久国内精品自在自线图片| videos熟女内射| 免费大片黄手机在线观看| 日日啪夜夜撸| 色哟哟·www| 中文乱码字字幕精品一区二区三区| 高清在线视频一区二区三区| 老司机影院毛片| 男女下面进入的视频免费午夜| 黄色怎么调成土黄色| 色婷婷久久久亚洲欧美| 最近最新中文字幕大全电影3| .国产精品久久| av免费观看日本| 一级毛片久久久久久久久女| 韩国av在线不卡| 少妇精品久久久久久久| 波野结衣二区三区在线| 亚洲怡红院男人天堂| 国产精品一区二区在线不卡| 久久久久网色| 热re99久久精品国产66热6| 日韩一区二区视频免费看| 欧美xxⅹ黑人| 美女主播在线视频| 亚洲av男天堂| 国产精品久久久久久av不卡| 色视频www国产| 欧美一级a爱片免费观看看| 国产精品一区二区在线观看99| 男人舔奶头视频| 高清日韩中文字幕在线| 天天躁日日操中文字幕| 亚洲国产成人一精品久久久| 欧美激情国产日韩精品一区| 久久久久久九九精品二区国产| 久久久久视频综合| 毛片女人毛片| 亚洲三级黄色毛片| 插逼视频在线观看| 99热国产这里只有精品6| 激情 狠狠 欧美| 久久久久精品性色| 观看免费一级毛片| 纵有疾风起免费观看全集完整版| 亚洲精品色激情综合| 久久国内精品自在自线图片| 亚洲精品乱码久久久v下载方式| 亚洲中文av在线| 少妇人妻 视频| 一区二区三区免费毛片| 丝瓜视频免费看黄片| 一个人看的www免费观看视频| 亚洲国产精品专区欧美| 王馨瑶露胸无遮挡在线观看| 乱码一卡2卡4卡精品| 午夜福利在线观看免费完整高清在| 久久久久精品性色| 又黄又爽又刺激的免费视频.| 日本免费在线观看一区| 色吧在线观看| 亚洲欧美日韩卡通动漫| 纯流量卡能插随身wifi吗| 黑丝袜美女国产一区| 久久久久久伊人网av| 99久久精品国产国产毛片| 一级毛片aaaaaa免费看小| 亚洲精品乱码久久久久久按摩| 搡老乐熟女国产| 欧美丝袜亚洲另类| 欧美激情国产日韩精品一区| 3wmmmm亚洲av在线观看| 超碰av人人做人人爽久久| 啦啦啦啦在线视频资源| 18+在线观看网站| 男的添女的下面高潮视频| 色哟哟·www| 建设人人有责人人尽责人人享有的 | 2021少妇久久久久久久久久久| 久久久久久人妻| 成人高潮视频无遮挡免费网站| 一区在线观看完整版| 在线观看免费高清a一片| 99热这里只有精品一区| 18+在线观看网站| 91aial.com中文字幕在线观看| 91在线精品国自产拍蜜月| 狠狠精品人妻久久久久久综合| 国产精品一区二区三区四区免费观看| av国产免费在线观看| 免费观看无遮挡的男女| 中文字幕亚洲精品专区| 在线天堂最新版资源| 成人特级av手机在线观看| 妹子高潮喷水视频| 一级毛片aaaaaa免费看小| 制服丝袜香蕉在线| 七月丁香在线播放| 国产亚洲精品久久久com| a级一级毛片免费在线观看| 中文欧美无线码| 久久久久久久精品精品| 久久久久久久久久久丰满| 国产爱豆传媒在线观看| 一级二级三级毛片免费看| 香蕉精品网在线| 午夜福利在线观看免费完整高清在| 少妇的逼好多水| 在线免费十八禁| 久久婷婷青草| 亚洲精品乱码久久久久久按摩| 免费大片18禁| 国产乱人偷精品视频| 身体一侧抽搐| 国产精品欧美亚洲77777| 亚洲av综合色区一区| av天堂中文字幕网| 国产精品一区二区在线不卡| 欧美 日韩 精品 国产| 夜夜看夜夜爽夜夜摸| 欧美精品国产亚洲| 性色av一级| 99精国产麻豆久久婷婷| 亚洲精品国产成人久久av| 中文字幕免费在线视频6| 三级国产精品欧美在线观看| 久久ye,这里只有精品| 美女福利国产在线 | 黄色欧美视频在线观看| 国产精品福利在线免费观看| 欧美成人一区二区免费高清观看| av视频免费观看在线观看| 欧美bdsm另类| 亚洲av免费高清在线观看| 高清视频免费观看一区二区| 久久久久久久久大av| 亚洲人与动物交配视频| 大香蕉久久网| 亚洲内射少妇av| 91精品国产九色| 国产精品久久久久久久电影| a级毛片免费高清观看在线播放| 国产精品久久久久久精品古装| 最近中文字幕2019免费版| 亚州av有码| 热re99久久精品国产66热6| 欧美日韩一区二区视频在线观看视频在线| 精品一区在线观看国产| 老师上课跳d突然被开到最大视频| 亚洲av国产av综合av卡| 国产黄片视频在线免费观看| 国产爽快片一区二区三区| 日韩av不卡免费在线播放| 少妇猛男粗大的猛烈进出视频| 成人综合一区亚洲| 国产真实伦视频高清在线观看| 毛片女人毛片| 亚洲欧美日韩东京热| 日韩av免费高清视频| 午夜精品国产一区二区电影| 三级国产精品欧美在线观看| 国产av国产精品国产| 日韩欧美一区视频在线观看 | 久久99精品国语久久久| 亚洲欧美清纯卡通| 一个人看视频在线观看www免费| 日韩中文字幕视频在线看片 | 久久久久久九九精品二区国产| 大片免费播放器 马上看| 18+在线观看网站| 日韩大片免费观看网站| 欧美丝袜亚洲另类| 亚洲国产高清在线一区二区三| 综合色丁香网| 国产欧美日韩精品一区二区| av国产免费在线观看| 夜夜骑夜夜射夜夜干| 熟女av电影| 777米奇影视久久| 久久精品国产亚洲av天美| 欧美精品一区二区大全| 国产成人a∨麻豆精品| 亚洲欧美成人精品一区二区| 爱豆传媒免费全集在线观看| 蜜桃亚洲精品一区二区三区| 日韩人妻高清精品专区| 精品久久久精品久久久| 免费大片黄手机在线观看| 激情 狠狠 欧美| 精品久久久噜噜| 又大又黄又爽视频免费| 三级国产精品片| 国产女主播在线喷水免费视频网站| 免费观看无遮挡的男女| 亚洲人成网站在线观看播放| 日本午夜av视频| av在线播放精品| 日韩 亚洲 欧美在线| 亚洲人成网站在线观看播放| 日韩av不卡免费在线播放| 九九久久精品国产亚洲av麻豆| 亚洲精品乱久久久久久| 久久久精品94久久精品| 在线 av 中文字幕| 人妻 亚洲 视频| 街头女战士在线观看网站| 在线 av 中文字幕| 老司机影院毛片| 啦啦啦中文免费视频观看日本| a级一级毛片免费在线观看| 国产精品秋霞免费鲁丝片| 免费看日本二区| 久久久午夜欧美精品| 国产亚洲精品久久久com| 夫妻性生交免费视频一级片| 综合色丁香网| 亚洲精品一二三| 欧美精品国产亚洲| 日日啪夜夜爽| 3wmmmm亚洲av在线观看| 日韩一区二区三区影片| 国产爱豆传媒在线观看| 欧美日韩国产mv在线观看视频 | 亚洲国产欧美在线一区| 国产在线一区二区三区精| av国产免费在线观看| 青春草国产在线视频| 中文天堂在线官网| 国产精品久久久久成人av| 一个人看的www免费观看视频| 日本爱情动作片www.在线观看| 日韩成人伦理影院| 人人妻人人添人人爽欧美一区卜 | 亚洲激情五月婷婷啪啪| 大香蕉97超碰在线| 网址你懂的国产日韩在线| 久久99热这里只频精品6学生| 18+在线观看网站| 精品久久国产蜜桃| 在线观看一区二区三区激情| 久久久久久伊人网av| 中文字幕制服av| 国产一区二区在线观看日韩| 亚洲av综合色区一区| av在线观看视频网站免费| 亚洲国产色片| 边亲边吃奶的免费视频| 国内精品宾馆在线| 国产av精品麻豆| 99久久综合免费| 国产成人freesex在线| 超碰av人人做人人爽久久| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 免费看不卡的av| 国产亚洲最大av| 亚洲无线观看免费| 18禁在线播放成人免费| 久久毛片免费看一区二区三区| 国产精品欧美亚洲77777| h视频一区二区三区| 亚洲欧洲日产国产| 人妻 亚洲 视频| 一级毛片电影观看| 精品视频人人做人人爽| 2018国产大陆天天弄谢| 国产淫语在线视频| 又粗又硬又长又爽又黄的视频| 国产成人精品一,二区| 亚洲av福利一区| 黄色欧美视频在线观看| 大片电影免费在线观看免费| 在线观看人妻少妇| 色吧在线观看| 国产精品麻豆人妻色哟哟久久| 日韩中文字幕视频在线看片 | 亚洲伊人久久精品综合| 成年av动漫网址| 亚洲国产av新网站| 久久精品国产a三级三级三级| 视频区图区小说| 久久精品国产鲁丝片午夜精品| 久久99热这里只有精品18| 亚洲av在线观看美女高潮| 一级毛片 在线播放| 少妇精品久久久久久久| 一级片'在线观看视频| 80岁老熟妇乱子伦牲交| 噜噜噜噜噜久久久久久91| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久精品精品| 国产欧美亚洲国产|