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

    炭化蘋果枝對蘋果根區(qū)土壤細(xì)菌和真菌多樣性的影響

    2016-12-01 08:10:22李燕歌周春然寧留芳楊洪強(qiáng)
    中國農(nóng)業(yè)科學(xué) 2016年17期
    關(guān)鍵詞:根區(qū)炭化群落

    曹 輝,李燕歌,周春然,寧留芳,楊洪強(qiáng)

    (山東農(nóng)業(yè)大學(xué)園藝科學(xué)與工程學(xué)院/作物生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室,山東泰安 271018)

    炭化蘋果枝對蘋果根區(qū)土壤細(xì)菌和真菌多樣性的影響

    曹 輝,李燕歌,周春然,寧留芳,楊洪強(qiáng)

    (山東農(nóng)業(yè)大學(xué)園藝科學(xué)與工程學(xué)院/作物生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室,山東泰安 271018)

    【目的】根區(qū)土壤微生物是影響根系環(huán)境的重要因素,炭化蘋果枝是廢棄果樹枝條低氧高溫?zé)峤猱a(chǎn)物,研究施用炭化蘋果枝對蘋果根區(qū)土壤細(xì)菌和真菌的群落結(jié)構(gòu)及其多樣性的影響,為炭化蘋果枝的合理應(yīng)用以及改善果園土壤生物學(xué)性狀提供理論依據(jù)?!痉椒ā吭诖杭荆瑢㈤L勢一致的2年生‘富士’蘋果幼樹(砧木為平邑甜茶)移栽到含有不同質(zhì)量比(0—4%)炭化蘋果枝的盆栽土壤中,于移栽120 d后采集土樣,提取基因組DNA,通過PCR擴(kuò)增建立文庫,利用Miseq平臺Illumina第二代高通量測序技術(shù)并結(jié)合相關(guān)生物信息學(xué)分析土壤細(xì)菌16S rRNA基因V3+V4區(qū)域和真菌ITS1區(qū)域的豐富度和多樣性指數(shù)以及群落結(jié)構(gòu)。【結(jié)果】從15個(gè)蘋果根區(qū)土壤樣本中獲得16 656個(gè)細(xì)菌分類操作單元(OTU)和435個(gè)真菌OTU,其中,變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)和酸桿菌門(Acidobacteria)是優(yōu)勢細(xì)菌,其相對豐度共占70.68%—72.80%;擔(dān)子菌門(Basidiomycota)、子囊菌門(Ascomycota)和接合菌門(Zygomycota)是優(yōu)勢真菌,其相對豐度共占68.00%—75.14%。群落物種豐富度指數(shù)(Chao指數(shù)和Ace指數(shù))分析顯示,1%炭化蘋果枝增加了細(xì)菌的群落豐富度,其Chao指數(shù)提高了15.42%,Ace指數(shù)提高3.89%;0.5%炭化蘋果枝增加了真菌的群落豐富度,其Chao指數(shù)提高了2.80%,Ace指數(shù)提高了3.61%。群落多樣性指數(shù)(Shannon指數(shù)和Simpson指數(shù))分析表明,0.5%—4%炭化蘋果枝降低了土壤細(xì)菌的多樣性,增加了土壤真菌的多樣性;其中,細(xì)菌群落多樣性Shannon指數(shù)在炭化蘋果枝用量為1%時(shí)最低,真菌Shannon指數(shù)在炭化蘋果枝用量為0.5%時(shí)最高。1%、2%和4%的炭化蘋果枝均不同程度地降低了根區(qū)土壤變形菌門、酸桿菌門和擔(dān)子菌門的相對豐度,提高了擬桿菌門和接合菌門的相對豐度;真菌擔(dān)子菌門中的銹革孔菌科(Hymenochaetaceae)的相對豐度最大(數(shù)值在31.99%—46.74%),施用1%、2%和4%的炭化蘋果枝可使其降低。0.5%—4%的炭化蘋果枝能增加不同處理所獨(dú)有的細(xì)菌OTU數(shù)目,改變細(xì)菌類群組成,其中獨(dú)有的細(xì)菌OTU數(shù)目是共有OTU數(shù)目的1—3倍,但對真菌類群組成沒有明顯影響?!窘Y(jié)論】施用0.5%—4%(w/w)炭化蘋果枝明顯改變蘋果根區(qū)土壤細(xì)菌和真菌的豐富度和多樣性,增加各用量下所特有的細(xì)菌物種,其中1%的炭化蘋果枝明顯提高根區(qū)土壤細(xì)菌的豐富度。

    炭化蘋果枝;根區(qū)土壤;Illumina高通量測序;細(xì)菌;真菌

    0 引言

    【研究意義】微生物是構(gòu)成土壤肥力的重要因素,但由于氮素化肥的大量施用,導(dǎo)致果園土壤氮素含量偏高,有機(jī)質(zhì)(有機(jī)碳)含量偏低,致使土壤微生物結(jié)構(gòu)惡化,土壤肥力普遍下降,嚴(yán)重威脅果樹生產(chǎn)的可持續(xù)發(fā)展,因此,迫切需要通過增碳減氮等措施改善微生物群落結(jié)構(gòu),提高土壤質(zhì)量?!厩叭搜芯窟M(jìn)展】有機(jī)物料經(jīng)過高溫炭化后施入土壤,不僅促進(jìn)廢棄物的資源化利用,還具有給土壤增碳和促進(jìn)植物生長的作用。閆麗娟等[1]報(bào)道炭化蘋果枝施入土壤可促進(jìn)平邑甜茶葉片的光合作用、根系發(fā)育和植株生長。有機(jī)物料高溫炭化產(chǎn)物是生物炭,生物炭作為一種土壤改良劑在農(nóng)林和環(huán)境等領(lǐng)域具有重要應(yīng)用價(jià)值[2]。大量研究發(fā)現(xiàn),向土壤中添加生物炭,可改善土壤結(jié)構(gòu)[3],干擾土壤硝酸鹽代謝[4],提高土壤pH、電導(dǎo)率(EC)、陽離子交換量[5]、養(yǎng)分持留能力、保水能力等[6],并因此能夠?yàn)橥寥牢⑸飾⑻峁┝己玫摹拔h(huán)境”[7]。土壤微生物不僅是土壤的重要組成部分,更是土壤養(yǎng)分循環(huán)的主要推動者,其在土壤有機(jī)物質(zhì)分解、養(yǎng)分釋放、能量轉(zhuǎn)移等生物地化循環(huán)中起著重要作用[8]。生物炭具有眾多孔隙和巨大比表面積,還有碳、磷和鉬等多種營養(yǎng)元素,這不僅可有效吸附微生物,還能為微生物提供所需營養(yǎng)物質(zhì)和優(yōu)越的棲息地[9]。LEHMANN等[10]認(rèn)為生物炭通過為土壤細(xì)菌提供更多碳源及改善土壤細(xì)菌生存環(huán)境條件而促進(jìn)了細(xì)菌的生長和繁衍。土壤微生物特別是細(xì)菌對自然或人為活動引起的環(huán)境變化十分敏感,曾有人把細(xì)菌群落作為評估土壤環(huán)境變化的指示性標(biāo)志[11]。SOLAIMAN等[12]發(fā)現(xiàn)腐生菌因其豐富的菌絲和高活性胞外酶而便于在生物炭中生長,而生物炭對腐生菌也起保護(hù)作用,它能促進(jìn)真菌與根系共生以及真菌菌根的生長;WARNOCK等[13]研究發(fā)現(xiàn)生物炭和菌根真菌互生能影響微生物菌群,適量菌根真菌和生物炭添加還能改善植物養(yǎng)分吸收及緩解鹽脅迫[14]。物種多樣性是分析土壤微生物狀況的重要參數(shù),常采用群落豐富度指數(shù)(Chao和Ace)和群落多樣性指數(shù)(Shannon和Simpson)等反映[15]。在土壤微生物多樣性研究中,前人主要采用瓊脂培養(yǎng)基培養(yǎng)方法、Biolog平板法、磷脂脂肪酸和核酸分析法等來分析微生物,但對不可培養(yǎng)的微生物無能為力,試驗(yàn)中還會受到主觀因素、操作繁瑣、數(shù)據(jù)量小等難題的困擾[8]。而高通量測序可通過檢測土壤微生物細(xì)胞內(nèi)特定遺傳物質(zhì)(原核微生物16s rDNA/rRNA、真核微生物18s rDNA/rRNA)的堿基序列,經(jīng)過對這些序列的對比分析,可探究并揭示土壤微生物物種和群落結(jié)構(gòu)的多樣性[16]。【本研究切入點(diǎn)】近年來高通量測序技術(shù)的不斷發(fā)展,其無需培養(yǎng)、能夠客觀還原菌群結(jié)構(gòu)和數(shù)據(jù)量大等優(yōu)點(diǎn)為深入分析微生物種類及其多樣性提供了可能[17]。此外,盡管生物炭可改善土壤微生物群落結(jié)構(gòu),但是其效果常因生物炭的制備材料等不同而存在很大差異[18-19];炭化蘋果枝是一種不同于以稻殼、秸稈等為材料制備的生物炭,其對土壤細(xì)菌和真菌多樣性的影響目前還不清楚。【擬解決的關(guān)鍵問題】本研究以冬季修剪下來的蘋果枝為原料,通過高溫低氧處理后施入土壤,采用高通量測序技術(shù),探討炭化蘋果枝對蘋果根區(qū)土壤細(xì)菌、真菌多樣性的影響,揭示根區(qū)土壤微生物菌落結(jié)構(gòu)在施用炭化蘋果枝后的變化,以期為炭化蘋果枝的合理應(yīng)用以及改善果園土壤生物學(xué)性狀提供理論依據(jù)。

    1 材料與方法

    試驗(yàn)于2014年3—11月在山東農(nóng)業(yè)大學(xué)園藝試驗(yàn)站和山東省高校果樹生物學(xué)重點(diǎn)實(shí)驗(yàn)室進(jìn)行。

    1.1 試驗(yàn)材料

    供試植株為長勢一致的2年生盆栽紅富士蘋果幼樹,砧木為平邑甜茶;盆栽土壤由壤土、兔糞、沙按照3∶1∶1(v/v/v)比例均勻混合而成,土壤基本理化指標(biāo):容重1.45 g·cm-3、pH 7.1、速效磷32.9 mg·kg-1、堿解氮26.6 mg·kg-1、速效鉀102.5 mg·kg-1、有機(jī)質(zhì)12.41 g·kg-1。炭化蘋果枝是將冬季修剪下來的蘋果枝條粉碎后用高溫炭化爐在700℃下低氧活化2 h制得;采用JSM-5610LV型掃描電鏡,先將炭化蘋果枝放入鍍膜機(jī)內(nèi)鍍膜,之后將其放入載物臺并用導(dǎo)電膠將其固定后裝入機(jī)箱進(jìn)行照片采集,其超微結(jié)構(gòu)見圖1、理化性質(zhì)見表1。

    1.2 試驗(yàn)設(shè)計(jì)

    2014年4月8日將盆栽蘋果盆土倒出,分別與質(zhì)量比(w/w)0(CK)、0.5%、1.0%、2.0%、4.0%炭化蘋果枝混勻后裝盆、栽樹,使樹根在盆中均勻分布。每個(gè)炭化蘋果枝用量為一個(gè)處理,每處理重復(fù)3次。處理120 d后開始取樣,每盆采集3點(diǎn)混合土樣,用直徑2 cm不銹鋼土鉆鉆取0—20 cm深度土壤,去除根系、雜草、土壤動物和石塊等雜質(zhì),混勻后分成2份,一份于-20℃冷凍保存用于土壤微生物分析,一份放進(jìn)4℃冰箱備測。

    圖1 炭化蘋果枝橫、縱徑電鏡掃描Fig. 1 SEM figure of carbonized apple branches

    表1 炭化蘋果枝的化學(xué)特性Table 1 Chemical characteristics of carbonized apple branches used in this study

    1.3 試驗(yàn)方法

    1.3.1 土壤微生物基因組DNA的提取 采用E.Z.N.A.?Soil DNA Kit(OMEGA,美國)的試劑盒方法,稱取0.5 g于-20℃保存的土壤樣品,按試劑盒的試驗(yàn)步驟進(jìn)行土壤微生物總DNA的提取,DNA樣品于-20℃冰箱保存待用。

    1.3.2 試驗(yàn)流程 提取樣品總DNA后,根據(jù)設(shè)計(jì)得到細(xì)菌V3+V4(F:5′-ARACTYCTACGGRAGGCW G-3′;R:5′-GACTACNVGGGTATCTAATCC-3′)和真菌ITS1(F:5′-AACCTGCGGAAGGATCATT-3′;R:5′-GARCCAAGAGATCCRTTG-3′)合成的引物,合并引物接頭,進(jìn)行PCR擴(kuò)增并對其產(chǎn)物進(jìn)行純化、定量和均一化形成測序文庫,建好的文庫先進(jìn)行文庫質(zhì)檢,質(zhì)檢合格的文庫用Illumina MiSeq PE300進(jìn)行測序,由北京百邁客生物科技有限公司完成測序。

    1.3.3 信息分析流程 對原始測序序列進(jìn)行質(zhì)量評估,拼接雙端Reads并過濾得到優(yōu)化序列,用于后續(xù)生物信息學(xué)的分析。將優(yōu)化序列進(jìn)行物種注釋并聚類得到OTU(Operational Taxonomic Units,操作分類單元),根據(jù)OTU的序列組成得到其物種分類?;贠TU聚類分析結(jié)果,對樣品的多樣性及組間差異性進(jìn)行分析,并結(jié)合物種分類信息,在各個(gè)分類水平上進(jìn)行分類學(xué)分析(圖2)。

    圖2 微生物多樣性生物信息分析流程圖Fig. 2 Bioinformatics analysis of microbial diversity process

    1.4 生物信息分析方法

    1.4.1 數(shù)據(jù)過濾及質(zhì)量評估 測序得到的原始序列(Sequenced reads)里含有帶接頭的、低質(zhì)量的reads,為了保證信息分析質(zhì)量,需要對原始測序序列進(jìn)行過濾,得到高質(zhì)量的reads,用于后續(xù)信息分析。數(shù)據(jù)過濾的主要步驟[20]如下:(1)去除帶接頭的reads;(2)若一條reads上N(未能確定出具體的堿基類型)的比例大于5%,則過濾掉該雙端reads;(3)去除低質(zhì)量reads。

    1.4.2 序列優(yōu)化 利用Mothur軟件(version 1.34.4,http://www.mothur.org/),進(jìn)行序列優(yōu)化,主要步驟[21]如下:(1)對序列進(jìn)行篩選,去除模糊堿基數(shù)大于0、單堿基高重復(fù)區(qū)大于8、重疊區(qū)錯(cuò)配數(shù)大于0及長度大于97.5%的序列(細(xì)菌)、重疊區(qū)長小于20 bp的序列(真菌);(2)進(jìn)行去冗余處理。

    1.4.3 OTU聚類分析 OTU即操作分類單元,是在系統(tǒng)發(fā)生學(xué)研究或群體遺傳學(xué)研究中,為了便于進(jìn)行分析,人為給某一個(gè)分類單元(品系、種、屬、分組等)設(shè)置的同一標(biāo)志。一般情況下,如果序列之間的相似性高于97%就可以把它定義為一個(gè)OTU,每個(gè)OTU對應(yīng)于一種不同的序列[22]??筛鶕?jù)不同的相似度水平,對所有序列進(jìn)行OTU劃分,通常在97%的相似水平下對OTU進(jìn)行生物信息統(tǒng)計(jì)分析。利用Mothur軟件以平均鄰近聚類算法(Average Neighbor Clustering Algorithm)在0.03(或97%的相似度)水平下進(jìn)行OTU的聚類,并統(tǒng)計(jì)獲得OTU的個(gè)數(shù)。

    1.4.4 豐富度指數(shù)和多樣性指數(shù)計(jì)算 選取Chao和Ace指數(shù)反映微生物豐富度,用Simpson和Shannon指數(shù)反映微生物多樣性[23-24]。

    Chao指數(shù):是用Chao1算法估計(jì)樣品中所含的OTU數(shù)目的指數(shù),Chao1在生態(tài)學(xué)中常用來估計(jì)物種總數(shù)。計(jì)算公式:

    式中,Schao1=估計(jì)的OTU數(shù);Sobs=實(shí)際觀測到的OTU數(shù);n1=只含有一條序列OTU數(shù)目;n2=只含有兩條序列的OTU數(shù)目。

    Ace指數(shù):用來估計(jì)群落中含有OTU數(shù)目的指數(shù),由Chao提出,是生態(tài)學(xué)中估計(jì)物種總數(shù)的常用指數(shù)之一,與Chao1的算法不同。計(jì)算公式:

    式中,ni=含有i條序列的OTU數(shù);Srare=低于閾值的OTU數(shù);Sabund=超過閾值的OTU數(shù);abund=含有豐富OTU數(shù)的閾值。

    Simpson指數(shù):用來估計(jì)樣品中微生物多樣性指數(shù)之一,在生態(tài)學(xué)中常用來定量描述一個(gè)區(qū)域的多樣性。Simpson指數(shù)越大,說明群落多樣性越低。

    式中,Sobs=實(shí)際觀測到的OTU數(shù);ni=含有i條序列的OTU數(shù);N=所有的序列數(shù)。

    Shannon指數(shù):用來估算樣品中微生物多樣性指數(shù)之一,它與Simpson多樣性指數(shù)常用于反映alpha多樣性指數(shù),Shannon值越大,說明群落多樣性越高。式中,Sobs=實(shí)際觀測到的OTU數(shù);ni=含有i條序列的OTU數(shù);N=所有的序列數(shù)。

    1.5 數(shù)據(jù)分析

    隨機(jī)選取相似度在97%條件下的OTU生成稀釋曲線,并利用軟件Mothur計(jì)算豐富度指數(shù)Chao和ACE,多樣性指數(shù)Simpson和Shannon?;赗DP和UNITE分類學(xué)數(shù)據(jù)庫對OTU進(jìn)行物種注釋,并用Excel和SPSS進(jìn)行數(shù)據(jù)處理,利用Excel和R語言工具對樣品物種組成及相對豐度統(tǒng)計(jì)結(jié)果繪制柱狀圖和Veen圖。

    2 結(jié)果

    2.1 土壤樣品測序深度評估

    Illumina第二代高通量測序結(jié)果顯示,15個(gè)土壤細(xì)菌樣品測序、過濾后共獲得482 Mb數(shù)據(jù)量,雙端reads拼接后共產(chǎn)生892 071條Raw Tags,序列優(yōu)化后共得到78 237條Clean Tags,聚類一共獲得16 656個(gè)OTU;15個(gè)土壤真菌樣品測序、過濾后共獲得1Gb數(shù)據(jù)量,雙端reads拼接后共產(chǎn)生2 049 567條Raw Tags,序列優(yōu)化后共得到1 103 018條Clean Tags,聚類一共獲得435個(gè)OTU。

    參照文獻(xiàn)[25]隨機(jī)抽取測序序列,將抽到的序列數(shù)與它們所能代表OTU的數(shù)目構(gòu)建曲線,在97%相似性水平下聚類OTU并制作各樣品的稀釋曲線圖(圖3),由圖可知,細(xì)菌和真菌曲線逐漸趨向平坦,說明測序數(shù)量合理。

    2.2 炭化蘋果枝對根區(qū)土壤細(xì)菌和真菌群落豐富度和多樣性的影響

    Ace和Chao指數(shù)可反映群落物種豐富度,由表2可知,土壤細(xì)菌Ace和Chao指數(shù)在施用1%(w/w)炭化蘋果枝后增加,增加幅度分別為15.42%和3.89%;在施用0.5%、2%和4%后降低,其中,Ace指數(shù)分別降低52.12%、10.87%和44.67%,Chao指數(shù)分別降低48.69%、18.54%和41.55%。真菌的Ace和Chao指數(shù)在施用0.5%炭化蘋果枝后分別增加,增加幅度分別為2.80%和3.61%;在施用1%、2%和4%炭化蘋果枝后降低。其中,Ace指數(shù)分別降低0.36%、0.75%和0.26%,Chao指數(shù)分別降低0.21%、1.80%和1.95%。顯示土壤細(xì)菌和真菌群落物種豐富度的變化因炭化蘋果枝的用量而有較大差異,施用1%(w/w)炭化蘋果枝可增加土壤細(xì)菌群落豐富度、0.5%的炭化蘋果枝則增加土壤真菌群落豐富度,施用2%和4%炭化蘋果枝均降低細(xì)菌和真菌的群落豐富度。

    表2數(shù)據(jù)顯示,施用0.5%、1%、2%和4%炭化蘋果枝后,土壤細(xì)菌Shannon指數(shù)均有所下降,下降幅度分別為2.90%、4.01%、0.14%和3.18%,在1%的用量時(shí)細(xì)菌Shannon指數(shù)最低;土壤真菌Shannon指數(shù)在施用0.5%、1%、2%和4%炭化蘋果枝后均有所提高,提高幅度分別為14.69%、3.13%、10.31%和7.81%,在0.5%的用量下真菌Shannon指數(shù)最高;Simpson指數(shù)的變化與之相反??梢?,向土壤中添加0.5%—4%的炭化蘋果枝降低土壤細(xì)菌多樣性而提高真菌多樣性。

    圖3 土壤樣品細(xì)菌(A)、真菌(B)多樣品稀釋曲線Fig. 3 Bacterial (A) and fungal (B) rarefaction curves of different soil samples

    表2 炭化蘋果枝處理下根區(qū)土壤樣本群落豐富度和多樣性指數(shù)Table 2 OTUs’ abundance and diversity of root-zone soil samples with carbonized apple branches applied

    2.3 炭化蘋果枝對根區(qū)土壤細(xì)菌和真菌相對豐度的影響

    圖4為根區(qū)土壤細(xì)菌和真菌在不同分類水平上的物種柱狀圖。由圖4-A可知,蘋果根區(qū)土壤樣品中的細(xì)菌包含變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)、酸桿菌門(Acidobacteria)、放線菌門(Actinobacteria)、厚壁菌門(Firmicutes)等菌群,各處理間細(xì)菌群落組成相似;除了未確定種屬的其他菌類,變形菌門、擬桿菌門、酸桿菌門的相對豐度最高,分別為39.59%—44.03%、13.38%—18.50%和11.85%—15.39%,屬于優(yōu)勢細(xì)菌,它們的相對豐度合計(jì)為70.68%—72.80%。

    圖4-A還顯示,施用0.5%、1%、2%和4%炭化蘋果枝能不同程度的降低根區(qū)土壤變形菌門和酸桿菌門的相對豐度,提高了擬桿菌門的相對豐度,其中,0.5%、1%、2%和4%炭化蘋果枝使變形菌門的相對豐度分別比對照降低10.08%、5.56%、2.43%和5.93%,使擬桿菌門的相對豐度分別比對照提高了21.08%、26.01%、28.03%和38.27%,使酸桿菌門的相對豐度分別比對照降低了3.25%、12.61%、23.00%、20.60%。0.5%和1%炭化蘋果枝使放線菌門相對豐度分別提高17.96%和1.48%;2%和4%炭化蘋果枝則使放線菌門相對豐度分別降低24.55%和19.60%。

    圖4 炭化蘋果枝處理下土壤樣品細(xì)菌門水平(A)、真菌門(B)和科(C)水平相對豐度變化Fig. 4 Changes in relative abundance of soil bacterial species at phylum level (A) and fungal species at phylum (B) and family (C) levels with carbonized apple branches applied

    圖4-B顯示,蘋果根區(qū)土壤樣品中的真菌包含擔(dān)子菌門(Basidiomycota)、子囊菌門(Ascomycota)、接合菌門(Zygomycota)、壺菌門(Chytridiomycota)、球囊菌門(Glomeromycota)等,除了未確定種屬的其他菌類,擔(dān)子菌門(Basidiomycota)、子囊菌門(Ascomycota)、接合菌門(Zygomycota)相對豐度最高,分別為46.09%—56.15%、8.57%—15.11%和7.03%—12.75%,屬于優(yōu)勢細(xì)菌,它們的相對豐度合計(jì)為68.00%—75.14%。

    圖4-B還顯示,施用1%、2%和4%的炭化蘋果枝降低了擔(dān)子菌門的相對豐度(分別比對照降低15.09%、15.62%和3.15%),施用0.5%和2%炭化蘋果枝降低了子囊菌門的相對豐度(比對照降低6.64%和0.22%),而施用1%和4%炭化蘋果枝使子囊菌門的相對豐度分別提高了59.15%和64.60%。0.5%、1%、2%和4%炭化蘋果枝增加了接合菌門的相對豐度,分別增加了34.28%、43.39%、81.37%和1.42%。

    圖4-C顯示,根區(qū)土壤真菌在科水平上,除了未確定科類外,擔(dān)子菌門中的銹革孔菌科(Hymenochaetaceae)的相對豐度在各處理中均最大(數(shù)值在31.99%—46.74%)。施入不同用量炭化蘋果枝可引發(fā)土壤真菌在不同分類水平上相對豐度的差異性變化,其中施用0.5%炭化蘋果枝提高了銹革孔菌科相對豐度(比對照提高9.44%),而1%、2%和4%炭化蘋果枝使銹革孔菌科相對豐度比對照分別降低了9.72%、25.10%和19.01%。

    可見,變形菌門、擬桿菌門和酸桿菌門是根區(qū)土壤的優(yōu)勢細(xì)菌,擔(dān)子菌門、子囊菌門和接合菌門是根區(qū)土壤的優(yōu)勢真菌,銹革孔菌科是相對豐度最大的土壤擔(dān)子菌;施入1%、2%和4%的炭化蘋果枝均可不同程度地降低根區(qū)土壤變形菌門、酸桿菌門和擔(dān)子菌門以及擔(dān)子菌門中的銹革孔菌科的相對豐度,提高了擬桿菌門和接合菌門的相對豐度。

    2.4 炭化蘋果枝處理后根區(qū)土壤細(xì)菌、真菌類群分析

    在97%的相似度下,得到了每個(gè)樣品的OTU個(gè)數(shù),如表3所示,0、0.5%、1%、2%和4%炭化蘋果枝處理后分別獲得6 374、4 093、5 378、10 598和4 476個(gè)細(xì)菌OTU,以及422、421、419、423和426個(gè)真菌OTU。

    表3 蘋果根區(qū)土壤樣本測序獲得的細(xì)菌、真菌OTU序列讀數(shù)Table 3 Reads of observed soil bacterial and fungal OTUs

    Venn圖能夠反映組間或樣品之間共有和特有OTU數(shù)目,能夠直觀地表現(xiàn)出組間或樣品間OTU的重疊情況[26]。結(jié)合OTU所代表的物種,可以找出不同環(huán)境中的核心微生物。從圖5-A中可以看出,不同處理之間共有的細(xì)菌OTU數(shù)目為1 192個(gè),代表的物種分別屬于α-變形菌綱(Alphaproteobacteria)、β-變形菌綱(Betaproteobacteria)、酸桿菌綱-Gp6(Acidobacteria_Gp6)、γ-變形菌綱(Gammaproteobacteria)、δ-變形菌綱(Deltaproteobacteria)、放線菌綱(Actinobacteria)、擬桿菌綱(Bacteroidetes)等。未施用炭化蘋果枝(CK)的根區(qū)土壤不含有其獨(dú)有的細(xì)菌OTU,而施用不同量的炭化蘋果枝后,各用量下所獨(dú)有的細(xì)菌OTU數(shù)目是共有OTU數(shù)目的1—3倍,說明施用炭化蘋果枝能引起土壤細(xì)菌特異性變化,在根區(qū)土壤能夠增加該處理所特有細(xì)菌物種,但增加的物種,因用量而不同,其中0.5%的用量增加了柔膜菌綱(Mollicutes)、Spartobacteria等特有物種,1%的用量增加了迷蹤菌綱(Elusimicrobia)、裝甲菌綱(Armatimonadia)等特有物種,2%的用量增加了疣微菌綱(Verrucomicrobiae)、Chthonomonadetes等特有物種,4%的用量增加了Ignavibacteria等特有物種。

    從圖5-B中可以看出,各處理之間共有的真菌OTU數(shù)目為391個(gè),主要共有物種分別屬于傘菌綱(Agaricomycetes)、座囊菌綱(Dothideomycetes)、銀耳綱(Tremellomycetes)、子囊菌綱(Ascomycota)、糞殼菌綱(Sordariomycetes)、散囊菌綱(Eurotiomycetes)、壺菌綱(Chytridiomycetes)、球囊菌綱(Glomeromycetes)、圓盤菌綱(Orbiliomycetes)等。不同處理之間沒有各自獨(dú)有的真菌OTU,說明施用不同用量的炭化蘋果枝對土壤真菌類群組成沒有明顯影響。

    圖5 炭化蘋果枝處理后根區(qū)土壤細(xì)菌(A)、真菌(B)Venn圖Fig. 5 Venn diagrams of soil bacterial (A) and fungal (B) community with carbonized apple branches applied

    可見,土壤施入0.5%—4%的炭化蘋果枝能增加蘋果根區(qū)土壤中獨(dú)有的細(xì)菌OTU數(shù)目,改變細(xì)菌類群組成,但對真菌類群組成沒有明顯影響。

    3 討論

    生物炭施入土壤通過改變土壤理化性質(zhì)而直接或間接影響土壤微生物。本研究發(fā)現(xiàn)施入1%炭化蘋果枝提高了細(xì)菌豐富度指數(shù),施入0.5%、2%和4%時(shí)則降低了該指數(shù)(表2),說明不同用量炭化蘋果枝對土壤中細(xì)菌豐富度的影響不同,較低(0.5%)和較高(2%和4%)用量都不利于細(xì)菌物種數(shù)量的增加。生物炭的孔徑、揮發(fā)物及灰分的組成與含量等都會影響細(xì)菌在土壤中的繁衍和生存狀態(tài)[13],PIETIK?INEN等[27]認(rèn)為細(xì)菌能夠被吸附到生物炭的表面,使它們不易受土壤淋洗的影響,從而增加了土壤中細(xì)菌的數(shù)量;但由于生物炭的來源和熱解條件的不同,它們會因土壤類型和生物炭用量等因素的差異而對微生物種群和數(shù)量產(chǎn)生不同的影響[18-19]。不同用量的生物炭施入土壤,土壤的結(jié)構(gòu)和理化性狀也會發(fā)生不同改變,這種改變并不與用量呈現(xiàn)簡單的線性關(guān)系[6,28],因此,細(xì)菌物種豐富度也不會因生物炭用量的逐漸增加呈簡單的線性變化。比如,炭化蘋果枝比土壤輕,施入土壤后會使土壤容重下降,這樣土壤氧氣含量會因炭化蘋果枝用量的增加而逐漸提高,土壤細(xì)菌既有厭氧類群,也有好氧類群,它們各自的豐富度會因土壤容重的減小而互有升降;炭化蘋果枝pH為堿性,隨著施用量的增加,土壤微區(qū)pH會隨之增加,這樣對pH要求不同的細(xì)菌的豐富度也會互有升降。因此,施入不同量的炭化蘋果枝后整個(gè)細(xì)菌類群豐富度的變化會很復(fù)雜,要揭示細(xì)菌豐富度與炭化蘋果枝用量之間更細(xì)致的關(guān)系,還有待將土壤細(xì)菌分類由門具體到科、屬、種甚至小種的層面。真菌物種豐富度隨炭化蘋果枝用量的變化也存在類似的問題,不過,有報(bào)道指出在土壤pH較高時(shí),真菌豐度會降低[29],高用量的炭化蘋果枝會提高土壤微區(qū)pH,因此當(dāng)用量提高到1%以上時(shí)會降低土壤真菌物種的豐富度。

    本研究發(fā)現(xiàn)施用1%炭化蘋果枝增加了細(xì)菌物種豐富度,但降低了細(xì)菌物種多樣性;施用1%、2%和4%炭化蘋果枝降低了真菌物種豐富度,卻增加真菌物種多樣性,這主要是因?yàn)槲⑸锏亩鄻有杂晌锓N豐富度和均勻度共同決定,它與物種豐富度的關(guān)系會因豐富度和均勻度二者所占比重而不同[30]。施用炭化蘋果枝一方面促進(jìn)某些微生物生長,同時(shí)也會抑制某類微生物生長,使微生物群落結(jié)構(gòu)發(fā)生變化,導(dǎo)致物種均勻度改變,并進(jìn)而影響微生物多樣性與物種豐富度的關(guān)系。炭化蘋果枝像其他生物炭一樣可吸附微生物,改變微生物生長環(huán)境等[31],這必然會引起微生物數(shù)量差異或是群落結(jié)構(gòu)的變化,影響群落的物種均勻度和豐富度,并進(jìn)而導(dǎo)致多樣性指數(shù)的變化(表2、圖4)。

    基于RDP分類學(xué)數(shù)據(jù)庫和UNITE分類學(xué)數(shù)據(jù)庫的分析顯示,在蘋果根區(qū)土壤中,變形菌門、擬桿菌門和酸桿菌門等是主要細(xì)菌類群,擔(dān)子菌門、子囊菌門、接合菌門等是主要真菌類群(圖4),這與前人對土壤微生物群落的研究結(jié)果一致[32]。0.5%—4%炭化蘋果枝能不同程度的降低變形菌門和酸桿菌門的相對豐度,提高擬桿菌門的相對豐度。變形菌門是細(xì)菌中最大的一類,所試樣品中的變形菌主要有α-變形菌、β-變形菌、γ-變形菌、δ-變形菌等,多樣性很高,數(shù)量也很大,屬于蘋果根區(qū)土壤中的優(yōu)勢細(xì)菌。酸桿菌門大多都是嗜酸菌,其豐富度與土壤pH呈顯著的負(fù)相關(guān),在土壤pH較低的環(huán)境中酸桿菌豐度最高[33],因此,炭化蘋果枝可通過提高土壤微區(qū)pH而降低酸桿菌門的相對豐度。擬桿菌門的細(xì)菌與DNA、蛋白質(zhì)和脂類等有機(jī)物質(zhì)的轉(zhuǎn)換緊密相關(guān),這些有機(jī)物質(zhì)的吸收和利用是碳循環(huán)的重要組成部分[34],炭化蘋果枝作為一種碳源,施入土壤有利于提高土壤擬桿菌門的相對豐度。不同用量炭化蘋果枝能引起真菌群落差異性變化(圖4-B),其中,0.5%—4%炭化蘋果枝能不同程度提高接合菌門的相對豐度,這種不一致的變化可能是由于炭化蘋果枝本身有效碳的缺乏抑制了某些真菌的定殖,而表面吸附的溶解性有機(jī)碳可以選擇性地被接合菌門利用[19]。大多數(shù)的病原菌來自于真菌,本研究顯示1%—4%炭化蘋果枝能降低真菌擔(dān)子菌門銹革孔菌科的相對豐度,而銹革孔菌科包含了很多與林木病有關(guān)的真菌[35],其相對豐度的降低會在一定程度上減少植株受病菌侵害的幾率。

    土壤細(xì)菌和真菌相對豐度的變化也可改變土壤微生物類群組成。本研究發(fā)現(xiàn)添加0.5%—4%炭化蘋果枝能增加蘋果根區(qū)土壤中獨(dú)有的細(xì)菌OTU數(shù)目,增加細(xì)菌特有物種,進(jìn)而改變了細(xì)菌類群組成,但對真菌類群組成沒有明顯影響(圖5),這主要在于炭化蘋果枝會像其他生物炭一樣能影響土壤理化性質(zhì)和生物學(xué)特性,能夠?yàn)槲⑸飫?chuàng)造特定的微環(huán)境[19],不同類群的微生物(比如真菌和細(xì)菌)對最適環(huán)境的要求不會完全一樣,炭化蘋果枝對不同菌群生存和繁衍的影響也不可能完全一樣。微生物對外界環(huán)境非常敏感,土壤環(huán)境的微小變動就會引起微生物多樣性的變化[11],土壤施入不同量炭化蘋果枝會使土壤環(huán)境發(fā)生不同程度的改變,這必然會引起微生物多樣性的差異性變化(表2)。通常認(rèn)為熟化程度高、肥力好的土壤微生物數(shù)量以及細(xì)菌所占的比例均較高,干旱和貧瘠土壤中的真菌和放線菌所占比例較高[36],本研究發(fā)現(xiàn),施用0.5%—4%的炭化蘋果枝使各用量下所特有的細(xì)菌物種增多(圖5),2%和4%的炭化蘋果枝降低真菌群落豐富度(表2)以及放線菌門的相對豐度(圖4-A)等,這些變化應(yīng)該有利于土壤肥力狀況的改善,具體有待進(jìn)一步研究。

    本研究還發(fā)現(xiàn)蘋果根區(qū)土壤細(xì)菌和真菌中還有很多未分類或是未確定種屬的物種,它們也受到炭化蘋果枝的影響,但其功能和特點(diǎn)尚不清楚,還有待通過深度測序或利用其他先進(jìn)手段對微生物進(jìn)行更細(xì)致地分類研究,并結(jié)合土壤的相關(guān)指標(biāo)進(jìn)行深入探索。

    4 結(jié)論

    高通量測序技術(shù)在蘋果根區(qū)土壤微生物研究中具有極大優(yōu)勢和可行性,施用0.5%—4%(w/w)炭化蘋果枝明顯改變蘋果根區(qū)土壤細(xì)菌和真菌的豐富度和多樣性,增加各用量下所特有細(xì)菌物種。炭化蘋果枝對土壤微生物的影響因施用劑量而呈現(xiàn)不同的變化趨勢,其中,1%的炭化蘋果枝明顯提高根區(qū)土壤細(xì)菌的豐富度,2%和4%的炭化蘋果枝降低真菌群落豐富度。施用適量炭化蘋果枝可改變微生物群落結(jié)構(gòu),將有利于改善土壤生物學(xué)性狀。

    [1] 閆麗娟, 楊洪強(qiáng), 蘇倩, 門秀巾, 張瑋瑋. 施用炭化蘋果枝粉末對平邑甜茶生長及根系構(gòu)型的影響. 園藝學(xué)報(bào), 2014, 41(7): 1436-1442. YAN L J, YANG H Q, SU Q, MEN X J, Zhang W W. Effects of carbonized powder of apple branch on the growth and root architecture of Malus hupehensis. Acta Horticulturae Sinica, 2014, 41(7): 1436-1442. (in Chinese)

    [2] RAJAPAKSHA A U, CHEN S S, TSANG D C, ZHANG M, VITHANAGE M, MANDAL S, GAO B, BOLAN N S, OK Y S. Engineered/designer biochar for contaminant removal/immobilization from soil and water: Potential and implication of biochar modification. Chemosphere, 2016, 148: 276-291.

    [3] LEHMANN J, RILLIG M C,THIES J, MASIELLO C A, HOCKADAY W C, CROWLEY D. Biochar effects on soil biota-A review. Soil Biology and Biochemistry, 2011, 43(9): 1812-1836.

    [4] 閆麗娟, 楊洪強(qiáng), 蘇倩, 門秀巾, 張瑋瑋. 炭化秸稈對蘋果根系一氧化氮生成及根區(qū)土壤硝酸鹽代謝的影響. 中國農(nóng)業(yè)科學(xué), 2014, 47(19): 3850-3856. YAN L J, YANG H Q, SU Q, MEN X J, ZHANG W W. Effects of carbonized straw on the nitric oxide formation and nitrate metabolism in apple roots and its root zone soil. Scientia Agricultura Sinica, 2014, 47(19): 3850-3856. (in Chinese)

    [5] JANUS A, PELFRêNE A, HEYMANS S, DEBOFFE C, DOUAY F, WATERLOT C. Elaboration, characteristics and advantages of biochars for the management of contaminated soils with a specific overview on Miscanthus biochars. Journal of Environmental Management, 2015, 162: 275-289.

    [6] MOLNáR M, VASZITA E, FARKAS é, UJACZKI é, FEKETEKERTéSZ I, TOLNER M, KLEBERCZ O, KIRCHKESZNER C, GRUIZ K, UZINGER N, FEIGL V. Acidic sandy soil improvement with biochar: A microcosm study. Science of the Total Environment, 2016, 563/564: 855-865.

    [7] LIANG B Q, LEHMANN J, SOHI S P, THIES J E, O’NEILL B, TRUJILLO L, GAUNT J, SOLOMON D, GROSSMAN J, NEVES E G, LUIZ?O F J. Black carbon affects the cycling of non-black carbon in soil. Organic Geochemistry, 2010, 41(2): 206-213.

    [8] 周桔, 雷霆. 土壤微生物多樣性影響因素及研究方法的現(xiàn)狀與展望. 生物多樣性, 2007, 15(3): 306-311. ZHOU J, LEI T. Review and prospects on methodology and affecting factors of soil microbial diversity. Biodiversity Science, 2007, 15(3): 306-311. (in Chinese)

    [9] FARRELL M, KUHN T K, MACDONALD L M, MADDERN T M, MURPHY D V, HALL P A, SINGH B P, BAUMANN K, KRULL E S, BALDOCK J A. Microbial utilisation of biochar-derived carbon. Science of Total Environment, 2013, 465: 288-297.

    [10] LEHMANN J, RILLIG M C, THIES J, MASIELLO C A, HOCKADAY W C, CROWLEY D. Biochar effects on soil biota-A review. Soil Biology and Biochemistry, 2011, 43(9): 1812-1836.

    [11] SHARMA S K, RAMESH A, SHARMA M P, JOSHI O P, GOVAERTS B, STEENWERTH K L, KARLEN D L. Microbial community structure and diversity as indicators for evaluating soil quality//Lichtfouse E. ed. Biodiversity, Biofuels, Agroforestry and Conservation Agriculture. Springer Netherlands, 2010, 5: 317-358.

    [12] SOLAIMAN Z M, BLACKWELL P, ABBOTT L K, STORER P. Direct and residual effect of biochar application on mycorrhizal root colonization, growth and nutrition of wheat. Australian Journal of Soil Research, 2010, 48(7): 546-554.

    [13] WARNOCK D D, LEHMANN J, KUYPER T W, RILLIG M C. Mycorrhizal responses to biochar in soil-concepts and mechanisms. Plant Soil, 2007, 300: 9-20.

    [14] HAMMER E C, FORSTREUTER M, RILLIG M C, KOHLEr J. Biochar increases arbuscular mycorrhizal plant growth enhancement and ameliorates salinity stress. Applied Soil Ecology, 2015, 96: 114-121.

    [15] 唐玉姝, 魏朝富, 顏廷梅, 楊林章, 慈恩. 土壤質(zhì)量生物學(xué)指標(biāo)研究進(jìn)展. 土壤, 2007(2): 157-163. TANG Y S, WEI C F, YAN T M, YANG L Z, CI E. Biological indicator of soil quality: A review. Soils, 2007(2): 157-163. (in Chinese)

    [16] CUMMINGS P J, AHMED R, DUROCHER J A, JESSEN A, VARDI T, OBOM K M. Pyrosequencing for microbial Identification and characterization. Journal of Visualized Experiments, 2013, 78: e50405.

    [17] 段曌, 肖煒, 王永霞, 賴泳紅, 崔曉龍. 454測序技術(shù)在微生物生態(tài)學(xué)研究中的應(yīng)用. 微生物學(xué)雜志, 2011, 31(5): 76-81. DUAN Z, XIAO W, WANG Y X, LAI Y H, CUI X L. Application of 454 sequencing technique in microbial ecology. Journal of Microbiology, 2011, 31(5): 76-81. (in Chinese)

    [18] KUPPUSAMY S, THAVAMANI P, MEGHARAJ M, VENKATESWARLU K, NAIDU R. Agronomic and remedial benefits and risks of applying biochar to soil: Current knowledge and future research directions. Environment International, 2016, 87: 1-12.

    [19] 丁艷麗, 劉杰, 王瑩瑩. 生物炭對農(nóng)田土壤微生物生態(tài)的影響研究進(jìn)展. 應(yīng)用生態(tài)學(xué)報(bào), 2013, 24(11): 3311-3317. DING Y L, LIU J, WANG Y Y. Effects of biochar on microbial ecology in agriculture soil: A review. Chinese Journal of Applied Ecology, 2013, 24(11): 3311-3317. (in Chinese)

    [20] EDGAR R C, HAAS B J, CLEMENTE J C, QUINCE C, KNIGHT R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics, 2011, 27(16): 2194-2200.

    [21] KOZICH J J, WESTCOTT S L, BAXTER N T, HIGHLANDER S K, SCHLOSS P D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Applied Environmental Microbiology, 2013, 79(17): 5112-5120.

    [22] BLAXTER M, MANN J, CHAPMAN T, THOMAS F, WHITTON C, FLOYD R, ABEBE E. Defining operational taxonomic units using DNA barcode data. Philosophical transactions of the Royal Society of London Series B Biological sciences, 2005, 360(1462): 1935-1943.

    [23] PITTA D W, PARMAR N, PATEL A K, INDUGU N, KUMAR S, PRAJAPATHI K B, PATEL A B, REDDY B, JOSHI C. Bacterial diversity dynamics associated with different diets and different primer pairs in the rumen of Kankrej cattle. PLoS ONE, 2014, 9(11):e111710.

    [24] SUN R B, ZHANG X X, GUO X S, WANG D Z, CHU H Y. Bacterial diversity in soils subjected to long-term chemical fertilization can be more stably maintained with the addition of livestock manure than wheat straw. Soil Biology and Biochemistry, 2015, 88: 9-18.

    [25] AMATO K R, YEOMAN C J, KENT A, RIGHINI N, CARBONERO F, ESTRADA A, GASKINS H R, STUMPF R M, YILDIRIM S, TORRALBA M, GILLIS M, WILSON B A, NELSON K E, WHITE B A, LEIGH S R. Habitat degradation impacts black howler monkey (Alouatta pigra) gastrointestinal microbiomes. The ISME Journal, 2013, 7: 1344-1353.

    [26] CHEN H, BOUTROS P C. Venn Diagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics, 2011, 12(1): 35.

    [27] PIETIK?INEN J, KIIKKIL? O, FRITZE H. Charcoal as a habitat for microbes and its effect on the microbial community of the underlying humus. Oikos, 2000, 89(2): 231-242.

    [28] CHEN X W, WONG J T, NG C W, WONG M H. Feasibility of biochar application on a landfill final cover: A review on balancing ecology and shallow slope stability. Environmental Science and Pollution Research, 2016, 23(8): 7111-7125.

    [29] ROUSK J, B??TH E, BROOKES P C, LAUBER C L, LOZUPONE C, CAPORASO J G, KNIGHT R, FIERER N. Soil bacterial and fungal communities across a pH gradient in an arable soil. The ISME Journal, 2010, 4(10): 1340-1351.

    [30] 劉根林. 轉(zhuǎn)基因大豆對根際根際土壤微生物群落的影響及其多樣性指數(shù)的度量[D]. 南京: 南京大學(xué), 2012. LIU G L. Effects of transgenic soybeans on the rhizospheric soil microbial communities and their diversity indices measurement [D]. Nanjing: Nanjing University, 2012. (in Chinese)

    [31] 張又弛, 李會丹. 生物炭對土壤中微生物群落結(jié)構(gòu)及其生物地球化學(xué)功能的影響. 生態(tài)環(huán)境學(xué)報(bào), 2015, 24(5): 898-905. ZHANG Y C, LI H D. Influence of biochar on the community structure and biogeochemical functions of microorganisms in soils. Ecology and Environmental Sciences, 2015, 24(5): 898-905. (in Chinese)

    [32] ROESCH L F, FULTHORPE R R, RIVA A, CASELLA G, HADWIN A K, KENT A D, DAROUB S H, CAMARGO F A, FARMERIE W G, TRIPLETT E W. Pyrosequencing enumerates and con-trasts soil microbial diversity. The ISME Journal, 2007, 1(4): 283-290.

    [33] JONES R T, ROBESON M S, LAUBER C L, HAMADY M, KNIGHT R, FIERER N. A comprehensive survey of soil acidobacterial diversity using pyrosequencing and clone library analyses. International Society for Microbial Ecology Journal, 2009, 3(4): 442-453.

    [34] MICHAND L, LO GINDICE A, TROUSSELLIER M, SMEDILE F, BRUNI V, BLANCHETON J P. Phylogenetic characterization of the heterotrophic bacterial communities inhabiting a marine recirculating aquaculture system. Journal of Applied Microbiology, 2009, 107(6): 1935-1946.

    [35] 戴玉成. 中國木本植物病原木材腐朽菌研究. 菌物學(xué)報(bào), 2012, 31(4): 493-509. DAI Y C. Pathogenic wood-decaying fungi on woody plants in China. Mycosystema, 2012, 31(4): 493-509. (in Chinese)

    [36] 周麗霞, 丁明懋. 土壤微生物學(xué)特性對土壤健康的指示作用. 生物多樣性, 2007, 15(2): 162-171. ZHOU L X, DING M M. Soil microbial characteristics as bioindicators of soil health. Biodiversity Science, 2007, 15(2): 162-171. (in Chinese)

    (責(zé)任編輯 趙伶俐)

    Effect of Carbonized Apple Branches on Bacterial and Fungal Diversities in Apple Root-Zone Soil

    CAO Hui, LI Yan-ge, ZHOU Chun-ran, NING Liu-fang, YANG Hong-qiang
    (College of Horticulture Science and Engineering, Shandong Agricultural University/State Key Laboratory of Crop Biology, Tai’an 271018, Shandong)

    【Objective】The microorganisms in soil of root-zone are important factors affecting root environment, carbonized apple branches are the low oxygen pyrolysis products of the abandoned fruit trees. The purpose of this study was to understand the structure of soil bacteria and fungi in apple root-zone and the response of their diversities to carbonized apple branches, and to provide a theoreticalbasis for the reasonable application of carbonized apple branches and the improvement of soil biological characters in orchard.【Method】 In the spring, the 2-year-old ‘Fuji’ apple trees (rootstock for Malus hupehensis Rehd) in similar growth were transplanted to the potting soil, the soil was mixed with different mass ratios (0-4%) of carbonized apple branches beforehand. Soil samples were collected after 120 days of transplanting, genomic DNA was extracted, and PCR amplification was made to establish libraries. In this study, the 16S rRNA genes V3+V4 regions of soil bacteria and fungal ITS1 regions were sequenced by Illumina high-throughput sequencing technology on Miseq platform, and related biological analysis was conducted to explore the changes of soil bacterial and fungal abundances, diversities and structures.【Result】A total of 16 656 bacterial operational taxonomic units (OTUs) and 435 fungal OTUs were obtained from 15 apple root-zone soil samples, among them, Proteobacteria, Bacteroidetes and Acidobacteria were the dominant bacteria which the total relative abundance was 70.68%-72.80%, and Basidiomycota, Ascomycota and Zygomycota were dominant fungi which the total relative abundance was 68.00%-75.14%. The richness indices of Chao and Ace showed that 1% (w/w) carbonized apple branches increased the abundance of bacteria by 15.42% and 3.89% compared with the control, respectively. 0.5% (w/w) carbonized apple branches increased the richness of fungi by 2.80% and 3.61%, respectively. Simpson and Shannon diversity index analysis showed that 0.5%-4% (w/w) carbonized apple branches reduced the diversity of soil bacteria, increased the diversity of soil fungi; bacteria Shannon diversity index was the lowest at 1% carbonized apple branches and fungi Shannon diversity index was the highest at 0.5% carbonized apple branches. Application of 1%, 2% and 4% carbonized apple branches to the soil, the relative abundance of Proteobacteria, Acidobacteria and Basidiomycota all decreased and the relative abundance of Bacteroidetes and Zygomycota all increased in apple root-zone soil at different degrees. The relative abundance of Hymenochaetaceae (31.99%-46.74%) was the highest in Basidiomycota. 1%, 2% and 4% carbonized apple branches all made it decrease. Carbonized apple branches at 0.5%-4% dosage could increase the number of bacteria unique OTUs significantly and affect the bacteria groups, wherein the number of bacterial unique OTU was one to three times of the common OTU numbers, but they had no significant effect on the fungal groups.【Conclusion】The 0.5%-4% (w/w) carbonized apple branch applied to the soil of apple root-zone, the abundance and diversity of soil bacteria and fungi had obviously changed, and the specific bacterial species under each dosage increased too. The soil bacterial richness increased significantly in root-zone applied 1% carbonized apple branches.

    carbonized apple branches; root-zone soil; Illumina high-throughput sequencing; soil bacteria; soil fungi

    2016-01-31;接受日期:2016-04-19

    國家“十二五”科技支撐計(jì)劃(2014BAD16B02)、國家自然科學(xué)基金(31372016)

    聯(lián)系方式:曹輝,E-mail:caohui716@163.com。通信作者楊洪強(qiáng),E-mail:hqyang@sdau.edu.cn, labft@sdau.edu.cn

    猜你喜歡
    根區(qū)炭化群落
    熱風(fēng)管道加溫下日光溫室根區(qū)溫度場的CFD模擬
    桉樹人工幼齡林根區(qū)和非根區(qū)土壤屬性特征分析
    大學(xué)生牙齦炎齦上菌斑的微生物群落
    合成微生物群落在發(fā)酵食品中的應(yīng)用研究
    水稻秸稈制備生物制活性碳
    市政污泥炭化時(shí)間與溫度研究
    安徽建筑(2018年4期)2019-01-29 02:16:14
    LED補(bǔ)光和根區(qū)加溫對日光溫室起壟內(nèi)嵌式基質(zhì)栽培甜椒生長及產(chǎn)量的影響*
    春季和夏季巢湖浮游生物群落組成及其動態(tài)分析
    炭化米糠經(jīng)臭氧活化制備活性炭及其去除Cr(VI)離子
    新型炭材料(2015年3期)2015-01-01 08:20:20
    樹盤施肥區(qū)域大小對 15N吸收利用及桃幼樹生長的影響
    极品教师在线免费播放| 欧美一区二区精品小视频在线| 欧美性长视频在线观看| 亚洲国产精品合色在线| 岛国在线观看网站| 欧美激情久久久久久爽电影| 日本熟妇午夜| 亚洲一区中文字幕在线| 超碰成人久久| 欧美日韩国产亚洲二区| 国内精品久久久久精免费| 法律面前人人平等表现在哪些方面| 欧美中文综合在线视频| 蜜桃久久精品国产亚洲av| 国产激情久久老熟女| 色尼玛亚洲综合影院| 亚洲男人的天堂狠狠| 法律面前人人平等表现在哪些方面| svipshipincom国产片| 丰满人妻熟妇乱又伦精品不卡| 女人被狂操c到高潮| 少妇的丰满在线观看| 欧美成人免费av一区二区三区| 又爽又黄无遮挡网站| 久久久国产成人精品二区| xxxwww97欧美| 欧美中文日本在线观看视频| 最新在线观看一区二区三区| 成年免费大片在线观看| 久久久国产成人精品二区| 欧美中文日本在线观看视频| 99久久综合精品五月天人人| 又爽又黄无遮挡网站| 久久久久国产一级毛片高清牌| 免费在线观看成人毛片| 97超级碰碰碰精品色视频在线观看| 夜夜夜夜夜久久久久| 精品欧美一区二区三区在线| 国产片内射在线| 最近最新中文字幕大全电影3| 欧美一级a爱片免费观看看 | 国产av一区二区精品久久| 黄色女人牲交| 黄片小视频在线播放| 波多野结衣高清作品| 亚洲熟妇熟女久久| 亚洲av日韩精品久久久久久密| 欧美性猛交╳xxx乱大交人| 两个人视频免费观看高清| 久久精品91无色码中文字幕| 制服诱惑二区| 99热6这里只有精品| 成人国语在线视频| 亚洲精品国产一区二区精华液| 日韩欧美国产在线观看| 窝窝影院91人妻| 国产亚洲精品一区二区www| 熟妇人妻久久中文字幕3abv| 老司机在亚洲福利影院| 日本成人三级电影网站| 精品国产超薄肉色丝袜足j| 亚洲精品中文字幕在线视频| 97碰自拍视频| 亚洲真实伦在线观看| 久久久久久九九精品二区国产 | 九九热线精品视视频播放| 亚洲国产欧美人成| 久久精品影院6| 九九热线精品视视频播放| 听说在线观看完整版免费高清| 亚洲欧美一区二区三区黑人| 亚洲国产精品成人综合色| 色精品久久人妻99蜜桃| 青草久久国产| 法律面前人人平等表现在哪些方面| 18禁黄网站禁片午夜丰满| 狂野欧美白嫩少妇大欣赏| 久久久国产成人精品二区| 国产片内射在线| 可以在线观看毛片的网站| 99久久99久久久精品蜜桃| 少妇熟女aⅴ在线视频| 国产亚洲精品久久久久久毛片| 动漫黄色视频在线观看| 人妻夜夜爽99麻豆av| 好男人电影高清在线观看| 嫩草影视91久久| 一a级毛片在线观看| 免费电影在线观看免费观看| 久久久久久大精品| 18禁美女被吸乳视频| 俺也久久电影网| 国产成人aa在线观看| 精品日产1卡2卡| 制服人妻中文乱码| 麻豆一二三区av精品| 老司机深夜福利视频在线观看| 又紧又爽又黄一区二区| 又粗又爽又猛毛片免费看| 欧美av亚洲av综合av国产av| 亚洲av美国av| 亚洲无线在线观看| 69av精品久久久久久| 欧美乱妇无乱码| 国产成人欧美在线观看| 亚洲男人的天堂狠狠| 欧美最黄视频在线播放免费| 日本熟妇午夜| 亚洲一区二区三区色噜噜| 波多野结衣高清无吗| 两个人免费观看高清视频| 国产免费av片在线观看野外av| 女生性感内裤真人,穿戴方法视频| 香蕉丝袜av| 午夜福利视频1000在线观看| av福利片在线| 国产v大片淫在线免费观看| 最近最新中文字幕大全电影3| 亚洲国产欧美网| 777久久人妻少妇嫩草av网站| 国产高清视频在线播放一区| 国产av在哪里看| 久久亚洲精品不卡| 色在线成人网| 欧美黑人欧美精品刺激| 视频区欧美日本亚洲| 久久天躁狠狠躁夜夜2o2o| 非洲黑人性xxxx精品又粗又长| 欧美日韩亚洲国产一区二区在线观看| 最近最新中文字幕大全电影3| 在线国产一区二区在线| 国产精品1区2区在线观看.| 我的老师免费观看完整版| 亚洲成人精品中文字幕电影| 国产精品,欧美在线| 琪琪午夜伦伦电影理论片6080| 色精品久久人妻99蜜桃| 无人区码免费观看不卡| 久久伊人香网站| 在线国产一区二区在线| 国产精品亚洲一级av第二区| 国产av在哪里看| 亚洲,欧美精品.| 国内毛片毛片毛片毛片毛片| www.www免费av| 少妇的丰满在线观看| 脱女人内裤的视频| 久热爱精品视频在线9| 亚洲精品一区av在线观看| 18禁国产床啪视频网站| 天天添夜夜摸| 久久性视频一级片| 欧美一级毛片孕妇| 高潮久久久久久久久久久不卡| 日韩国内少妇激情av| 色哟哟哟哟哟哟| 中文字幕人成人乱码亚洲影| АⅤ资源中文在线天堂| 人成视频在线观看免费观看| 久久久国产成人免费| 99热这里只有精品一区 | 真人一进一出gif抽搐免费| 最新在线观看一区二区三区| 两个人看的免费小视频| 国产伦一二天堂av在线观看| 亚洲av成人不卡在线观看播放网| 亚洲,欧美精品.| 日本一二三区视频观看| 精品午夜福利视频在线观看一区| 国产成人av激情在线播放| 亚洲国产精品合色在线| 亚洲精品在线观看二区| 99精品欧美一区二区三区四区| 欧美zozozo另类| 很黄的视频免费| 国内精品久久久久精免费| 成年版毛片免费区| 国产精品影院久久| 天堂√8在线中文| av超薄肉色丝袜交足视频| 美女 人体艺术 gogo| 国产高清videossex| 搞女人的毛片| 亚洲人成77777在线视频| 亚洲黑人精品在线| 中出人妻视频一区二区| 日日夜夜操网爽| 日本成人三级电影网站| av国产免费在线观看| 欧美性猛交黑人性爽| 禁无遮挡网站| 在线播放国产精品三级| 成人永久免费在线观看视频| 亚洲电影在线观看av| 99re在线观看精品视频| 一个人免费在线观看电影 | 色精品久久人妻99蜜桃| 国产精品亚洲av一区麻豆| 中文字幕久久专区| 午夜两性在线视频| 午夜a级毛片| 免费看美女性在线毛片视频| 最近在线观看免费完整版| 蜜桃久久精品国产亚洲av| 欧美日韩国产亚洲二区| 国产爱豆传媒在线观看 | 久久午夜综合久久蜜桃| 欧美日韩瑟瑟在线播放| 日本免费一区二区三区高清不卡| 香蕉av资源在线| 99久久无色码亚洲精品果冻| 在线看三级毛片| 窝窝影院91人妻| 亚洲电影在线观看av| 日韩欧美三级三区| 最近最新免费中文字幕在线| 精品午夜福利视频在线观看一区| netflix在线观看网站| 精品国产亚洲在线| 欧美黑人精品巨大| 亚洲av第一区精品v没综合| 久久久精品国产亚洲av高清涩受| 精品第一国产精品| 一级黄色大片毛片| 草草在线视频免费看| 亚洲av第一区精品v没综合| 精品欧美一区二区三区在线| 国产熟女xx| 两个人的视频大全免费| 亚洲av美国av| 99riav亚洲国产免费| 亚洲国产中文字幕在线视频| 欧美绝顶高潮抽搐喷水| av天堂在线播放| 久久久久久久午夜电影| 国产精华一区二区三区| xxx96com| 又紧又爽又黄一区二区| 久久精品国产亚洲av香蕉五月| 国产成人精品无人区| 欧美黄色淫秽网站| 国产午夜精品论理片| 亚洲精品在线观看二区| 亚洲片人在线观看| 亚洲aⅴ乱码一区二区在线播放 | 动漫黄色视频在线观看| 国产精品香港三级国产av潘金莲| 国产亚洲欧美在线一区二区| 欧美成狂野欧美在线观看| 精品国内亚洲2022精品成人| 大型黄色视频在线免费观看| 亚洲国产欧洲综合997久久,| 久久性视频一级片| 日韩欧美国产在线观看| 曰老女人黄片| 国内揄拍国产精品人妻在线| 69av精品久久久久久| 成熟少妇高潮喷水视频| 9191精品国产免费久久| 又粗又爽又猛毛片免费看| 亚洲狠狠婷婷综合久久图片| 国产成人一区二区三区免费视频网站| 国产精品野战在线观看| 亚洲成人久久性| 久久精品影院6| 中文字幕熟女人妻在线| 别揉我奶头~嗯~啊~动态视频| 一级毛片女人18水好多| 女同久久另类99精品国产91| 一进一出好大好爽视频| 俺也久久电影网| 久久热在线av| 久久婷婷成人综合色麻豆| 一二三四社区在线视频社区8| 精品熟女少妇八av免费久了| 亚洲精品av麻豆狂野| 日本免费一区二区三区高清不卡| 19禁男女啪啪无遮挡网站| bbb黄色大片| 99国产精品99久久久久| 神马国产精品三级电影在线观看 | 怎么达到女性高潮| 亚洲成人久久性| 欧美一区二区国产精品久久精品 | 人人妻人人看人人澡| 午夜免费成人在线视频| 免费看美女性在线毛片视频| 成人三级做爰电影| 国产黄片美女视频| 99久久精品热视频| 国模一区二区三区四区视频 | 91大片在线观看| 精品少妇一区二区三区视频日本电影| 男女下面进入的视频免费午夜| 一级片免费观看大全| 精品午夜福利视频在线观看一区| 老熟妇乱子伦视频在线观看| www日本黄色视频网| 一个人免费在线观看电影 | 欧美一级毛片孕妇| 亚洲色图av天堂| 亚洲精品久久成人aⅴ小说| 一本久久中文字幕| 国产99久久九九免费精品| 精品电影一区二区在线| 国产一区二区三区在线臀色熟女| 夜夜夜夜夜久久久久| 日日爽夜夜爽网站| 亚洲专区中文字幕在线| 九九热线精品视视频播放| 夜夜夜夜夜久久久久| 超碰成人久久| 一区福利在线观看| 99riav亚洲国产免费| 黑人欧美特级aaaaaa片| 麻豆久久精品国产亚洲av| 中文字幕熟女人妻在线| 很黄的视频免费| 一本综合久久免费| 在线视频色国产色| 99国产精品一区二区蜜桃av| 国产精品 欧美亚洲| 夜夜看夜夜爽夜夜摸| 日本精品一区二区三区蜜桃| 欧美色视频一区免费| 波多野结衣巨乳人妻| 国产av在哪里看| av视频在线观看入口| 久久伊人香网站| 欧美三级亚洲精品| 日韩三级视频一区二区三区| 国产麻豆成人av免费视频| 黄色毛片三级朝国网站| 97碰自拍视频| 国产黄色小视频在线观看| 成人亚洲精品av一区二区| 全区人妻精品视频| АⅤ资源中文在线天堂| 精品无人区乱码1区二区| 亚洲专区字幕在线| 啦啦啦免费观看视频1| 久久精品91无色码中文字幕| 久久国产精品人妻蜜桃| 又粗又爽又猛毛片免费看| 在线观看美女被高潮喷水网站 | 99riav亚洲国产免费| 成人一区二区视频在线观看| 99国产极品粉嫩在线观看| 免费搜索国产男女视频| 精品第一国产精品| 久久久久久人人人人人| 日本a在线网址| 中文字幕最新亚洲高清| 日韩欧美 国产精品| 一个人免费在线观看电影 | 黄色丝袜av网址大全| 免费在线观看影片大全网站| 免费看美女性在线毛片视频| 久久久久精品国产欧美久久久| 悠悠久久av| 国产黄色小视频在线观看| 久久久久精品国产欧美久久久| 成人av一区二区三区在线看| 国内毛片毛片毛片毛片毛片| 日韩欧美国产在线观看| 天天躁夜夜躁狠狠躁躁| 久久精品国产亚洲av高清一级| 亚洲第一电影网av| 黄色 视频免费看| 欧美日韩亚洲综合一区二区三区_| 欧美日韩亚洲国产一区二区在线观看| 中文亚洲av片在线观看爽| 琪琪午夜伦伦电影理论片6080| 狂野欧美激情性xxxx| 成人特级黄色片久久久久久久| 亚洲中文av在线| 女人被狂操c到高潮| 一区二区三区激情视频| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕精品亚洲无线码一区| 又黄又爽又免费观看的视频| 亚洲18禁久久av| 最新在线观看一区二区三区| 亚洲av五月六月丁香网| 丰满人妻一区二区三区视频av | 国产伦一二天堂av在线观看| 国产免费av片在线观看野外av| 精品免费久久久久久久清纯| 成年免费大片在线观看| 国产乱人伦免费视频| 波多野结衣高清作品| 国产区一区二久久| 亚洲精品一卡2卡三卡4卡5卡| 免费高清视频大片| 成人亚洲精品av一区二区| 日韩欧美 国产精品| 91av网站免费观看| 欧美另类亚洲清纯唯美| 色哟哟哟哟哟哟| 丝袜美腿诱惑在线| 精品久久久久久久久久免费视频| 久久欧美精品欧美久久欧美| 香蕉丝袜av| 91大片在线观看| 国产男靠女视频免费网站| 国产精品久久久久久人妻精品电影| 欧美最黄视频在线播放免费| 成人三级黄色视频| 亚洲欧美日韩高清专用| 久久久久免费精品人妻一区二区| 在线永久观看黄色视频| 亚洲人与动物交配视频| 婷婷精品国产亚洲av| 国产三级黄色录像| √禁漫天堂资源中文www| 999久久久国产精品视频| 精品免费久久久久久久清纯| 精品国内亚洲2022精品成人| 色噜噜av男人的天堂激情| 国产真实乱freesex| 两性午夜刺激爽爽歪歪视频在线观看 | 妹子高潮喷水视频| 最近视频中文字幕2019在线8| 日韩欧美国产在线观看| 波多野结衣高清无吗| 18禁裸乳无遮挡免费网站照片| 国产熟女xx| 国产成年人精品一区二区| 欧美成人免费av一区二区三区| 999久久久精品免费观看国产| 精品免费久久久久久久清纯| 国产探花在线观看一区二区| 欧美在线黄色| 在线观看免费午夜福利视频| 亚洲欧美精品综合一区二区三区| 91成年电影在线观看| 国产高清videossex| 色在线成人网| 国产欧美日韩一区二区精品| 色老头精品视频在线观看| 90打野战视频偷拍视频| 窝窝影院91人妻| av国产免费在线观看| 在线观看日韩欧美| 在线a可以看的网站| 国产一区二区三区在线臀色熟女| 一夜夜www| 久久 成人 亚洲| 欧美黄色片欧美黄色片| 亚洲av电影不卡..在线观看| 男人舔女人的私密视频| www.www免费av| 超碰成人久久| 九九热线精品视视频播放| 国产又黄又爽又无遮挡在线| 精品久久久久久,| 欧美日韩一级在线毛片| 一级黄色大片毛片| 一夜夜www| 亚洲av成人一区二区三| 亚洲av美国av| 欧美国产日韩亚洲一区| 国产在线观看jvid| 精品无人区乱码1区二区| 欧美三级亚洲精品| 亚洲av第一区精品v没综合| 午夜日韩欧美国产| 中文字幕人妻丝袜一区二区| 成人三级做爰电影| 婷婷精品国产亚洲av在线| 真人一进一出gif抽搐免费| 久热爱精品视频在线9| 国产黄色小视频在线观看| 99国产精品99久久久久| 亚洲精品色激情综合| 国产片内射在线| 在线观看一区二区三区| 久久久久久亚洲精品国产蜜桃av| 老熟妇仑乱视频hdxx| 男插女下体视频免费在线播放| 白带黄色成豆腐渣| 亚洲av第一区精品v没综合| 国产亚洲精品一区二区www| 亚洲精品国产精品久久久不卡| 精品第一国产精品| 久久99热这里只有精品18| 成熟少妇高潮喷水视频| 婷婷精品国产亚洲av在线| 久99久视频精品免费| 在线观看www视频免费| 一边摸一边抽搐一进一小说| 欧美激情久久久久久爽电影| 一个人免费在线观看电影 | 色综合婷婷激情| 黄色视频不卡| 欧美国产日韩亚洲一区| 国产又色又爽无遮挡免费看| 久久精品国产99精品国产亚洲性色| 97碰自拍视频| 国产成人av教育| 久久精品人妻少妇| 欧美色视频一区免费| 国产精品久久久久久久电影 | 国产黄片美女视频| 丝袜美腿诱惑在线| 国产av麻豆久久久久久久| 老汉色av国产亚洲站长工具| 两个人看的免费小视频| 黄色女人牲交| 可以免费在线观看a视频的电影网站| 精品人妻1区二区| 很黄的视频免费| 我的老师免费观看完整版| 久久午夜亚洲精品久久| 中文亚洲av片在线观看爽| 国产精品久久久久久人妻精品电影| 亚洲欧美日韩东京热| 亚洲性夜色夜夜综合| 成人一区二区视频在线观看| 两个人免费观看高清视频| 一卡2卡三卡四卡精品乱码亚洲| 人人妻人人澡欧美一区二区| 精品无人区乱码1区二区| 亚洲专区中文字幕在线| 别揉我奶头~嗯~啊~动态视频| 欧美日韩黄片免| 国产亚洲av嫩草精品影院| 欧美成人午夜精品| 最近最新中文字幕大全免费视频| 国产伦一二天堂av在线观看| 精品国产乱码久久久久久男人| netflix在线观看网站| 淫秽高清视频在线观看| 看片在线看免费视频| 国产黄a三级三级三级人| 精品一区二区三区av网在线观看| 在线免费观看的www视频| 一本精品99久久精品77| 亚洲成人国产一区在线观看| 日韩欧美国产在线观看| 午夜免费观看网址| 久久热在线av| 国产高清视频在线观看网站| 日本一本二区三区精品| 丰满人妻一区二区三区视频av | 日韩免费av在线播放| 日韩有码中文字幕| 成人av一区二区三区在线看| 精品欧美国产一区二区三| 亚洲欧美精品综合一区二区三区| 最近最新免费中文字幕在线| 成人国产综合亚洲| 两个人视频免费观看高清| 激情在线观看视频在线高清| 老鸭窝网址在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成人久久爱视频| 午夜福利高清视频| 无人区码免费观看不卡| 久久久久亚洲av毛片大全| 热99re8久久精品国产| av天堂在线播放| 免费在线观看亚洲国产| 亚洲美女黄片视频| 女人被狂操c到高潮| 成人三级黄色视频| 特大巨黑吊av在线直播| 国产不卡一卡二| 黄色丝袜av网址大全| 在线免费观看的www视频| 国产成人av激情在线播放| 天堂√8在线中文| 国产精品电影一区二区三区| 国产三级黄色录像| 欧美乱码精品一区二区三区| 欧美性猛交╳xxx乱大交人| 熟妇人妻久久中文字幕3abv| 在线观看美女被高潮喷水网站 | 看片在线看免费视频| 亚洲狠狠婷婷综合久久图片| 搡老熟女国产l中国老女人| 搡老妇女老女人老熟妇| 欧美不卡视频在线免费观看 | 国产精品一区二区三区四区免费观看 | 国产成人精品无人区| 欧美不卡视频在线免费观看 | 精品一区二区三区av网在线观看| 老汉色av国产亚洲站长工具| 99精品在免费线老司机午夜| 欧美性猛交黑人性爽| 国产午夜精品论理片| 欧美zozozo另类| 亚洲中文字幕日韩| 精品熟女少妇八av免费久了| 免费看a级黄色片| 老司机福利观看| 精品国产美女av久久久久小说| 久久国产精品人妻蜜桃| 九色成人免费人妻av| 韩国av一区二区三区四区| 欧美日韩乱码在线| 老熟妇乱子伦视频在线观看| 午夜免费激情av| 又粗又爽又猛毛片免费看| 可以在线观看的亚洲视频| 日本免费一区二区三区高清不卡| 三级毛片av免费| 18禁黄网站禁片免费观看直播| www.精华液| 一进一出好大好爽视频| 一级毛片女人18水好多| 淫秽高清视频在线观看| 一本大道久久a久久精品| 夜夜爽天天搞| 人妻夜夜爽99麻豆av|