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

    西瓜連作根際土壤微生物群落演替特征

    2023-11-18 07:00:24郭晗玥王東升阮楊喬亦鑄張蕓滔李玲黃啟為郭世偉凌寧沈其榮
    關(guān)鍵詞:菌門根際群落

    郭晗玥,王東升,阮楊,喬亦鑄,張蕓滔,李玲,黃啟為,郭世偉,凌寧,沈其榮

    西瓜連作根際土壤微生物群落演替特征

    1南京農(nóng)業(yè)大學(xué)資源與環(huán)境科學(xué)學(xué)院/國(guó)家有機(jī)肥重點(diǎn)實(shí)驗(yàn)室/江蘇省固體有機(jī)廢物利用重點(diǎn)實(shí)驗(yàn)室,南京 210095;2南京市蔬菜科學(xué)研究所,南京210042

    【目的】研究連作對(duì)西瓜根際土壤中細(xì)菌、真菌群落構(gòu)建和潛在功能的影響,明確根際微生物對(duì)環(huán)境變化的適應(yīng)性,為西瓜連作障礙的生態(tài)防治和農(nóng)田土壤的健康保持提供理論基礎(chǔ)?!痉椒ā恳晕催B作(對(duì)照)、連作2茬和連作 6茬的西瓜根際土壤為研究對(duì)象,利用16S rRNA和ITS高通量測(cè)序技術(shù)解析西瓜連作對(duì)根際土壤細(xì)菌、真菌群落特征和演替的影響。【結(jié)果】隨著連作茬次的增加,西瓜根際土壤中細(xì)菌多樣性指數(shù)呈現(xiàn)先下降后上升趨勢(shì),真菌多樣性指數(shù)顯著下降。在細(xì)菌屬水平,多茬次的連作降低了西瓜根際土壤中鞘氨醇單胞菌屬()和溶桿菌屬()的相對(duì)豐度;在真菌屬水平,鐮刀菌屬()相對(duì)豐度隨連作茬次增加呈增加趨勢(shì)。與對(duì)照相比,連作的網(wǎng)絡(luò)復(fù)雜度更高,但網(wǎng)絡(luò)結(jié)構(gòu)穩(wěn)定性低。此外,相比于對(duì)照,連作6茬顯著減少有害物質(zhì)的生物降解途徑以及氨基酸的代謝途徑的相對(duì)豐度,并且在連作6茬后真菌群落中病原真菌相對(duì)豐度顯著增多。在群落演替過程中,隨機(jī)性過程主導(dǎo)連作下西瓜根際細(xì)菌群落構(gòu)建,而確定性過程主導(dǎo)連作下西瓜根際真菌群落構(gòu)建?!窘Y(jié)論】連作引起根際細(xì)菌和真菌的群落特征、功能組成和演替過程變化,細(xì)菌群落關(guān)鍵功能的減少、病原真菌相對(duì)豐度的增加以及微生物群落網(wǎng)絡(luò)穩(wěn)定性的減小可能是導(dǎo)致西瓜連作障礙發(fā)生的重要因素。

    西瓜;連作土壤;根際微生物群落;功能預(yù)測(cè);群落演替

    0 引言

    【研究意義】西瓜(Thunb. Matsum and Naka)是我國(guó)乃至世界主要的夏季水果之一[1],其具有較高的經(jīng)濟(jì)價(jià)值以及巨大的市場(chǎng)需求,促使農(nóng)戶利用現(xiàn)代農(nóng)業(yè)技術(shù)進(jìn)行集約化生產(chǎn)。然而西瓜連年單作現(xiàn)象十分普遍[2],導(dǎo)致西瓜連作障礙愈發(fā)嚴(yán)重且成為了限制西瓜產(chǎn)業(yè)持續(xù)穩(wěn)定發(fā)展的重要因素[3]。長(zhǎng)期連作也會(huì)引起土壤理化性質(zhì)改變、土壤病原菌富集以及自毒物質(zhì)的積累,最終導(dǎo)致作物生長(zhǎng)受限[4-5]。土壤微生物在保持土壤功能和根際生態(tài)系統(tǒng)可持續(xù)性方面起著至關(guān)重要的作用,且通常被認(rèn)為是土壤健康的敏感生物指標(biāo)[6-7]。研究連作對(duì)西瓜根際土壤中細(xì)菌、真菌群落結(jié)構(gòu)和多樣性的影響,對(duì)防治西瓜連作障礙和保持農(nóng)田土壤的健康具有重要意義。【前人研究進(jìn)展】有研究表明,對(duì)連作土壤進(jìn)行滅菌能夠顯著提高植物的生長(zhǎng)和品質(zhì)[8],這表明生物因素是植物連作土壤病害發(fā)生的主要原因。植物根際擁有大量微生物,它們參與土壤養(yǎng)分循環(huán),影響植物對(duì)養(yǎng)分的吸收與利用,對(duì)于植物生長(zhǎng)起著至關(guān)重要的作用[9-10]。植物根際同時(shí)也是土傳病原菌進(jìn)入植物的必經(jīng)之地,在入侵植物體之前,病原菌會(huì)先與土著微生物競(jìng)爭(zhēng)生態(tài)位,再?gòu)母H感染植物體[5,11]。過往研究發(fā)現(xiàn),長(zhǎng)期的單一作物栽培所引起的土壤性質(zhì)惡化以及植物根系分泌物在內(nèi)的化感作用會(huì)改變根際微生物群落特征[12-14]。有研究表明,連作會(huì)簡(jiǎn)化土壤微生物群落生態(tài)網(wǎng)絡(luò)并降低網(wǎng)絡(luò)穩(wěn)定性,這可能與土壤理化性質(zhì)惡化、土壤肥力下降密切相關(guān)[15]。此外,了解微生物群落組裝的生態(tài)過程有助于識(shí)別微生物群落組成如何響應(yīng)環(huán)境變化,其中確定性過程和隨機(jī)性過程基于生態(tài)位理論和中性理論,通常用于解釋微生物群落構(gòu)建[16]。確定性過程包括非生物(pH、鹽度、有機(jī)碳含量[17-18])和生物(捕食、共生和寄生等[19])因素,決定物種的存在與否和相對(duì)豐度,與生態(tài)選擇相關(guān);隨機(jī)過程包括不可預(yù)測(cè)的擾動(dòng)、概率性的散布和隨機(jī)的出生-死亡事件等,這些變化不是由環(huán)境決定的適應(yīng)性結(jié)果。確定性過程和隨機(jī)過程在構(gòu)建生態(tài)群落中同時(shí)起作用來調(diào)節(jié)群落構(gòu)建[20-21]。【本研究切入點(diǎn)】目前對(duì)西瓜長(zhǎng)期連作體系下根際土壤微生物群落特征及其構(gòu)建過程研究較少,這對(duì)于預(yù)測(cè)根際微生物在調(diào)節(jié)植物-微生物生態(tài)系統(tǒng)功能中的作用是至關(guān)重要的?!緮M解決的關(guān)鍵問題】本研究以不同連作茬次的西瓜為研究對(duì)象,采用16S rRNA和ITS高通量測(cè)序技術(shù),研究連作對(duì)西瓜根際土壤細(xì)菌和真菌群落組成的影響,同時(shí)結(jié)合PICRUSt2和FUNGuild分別分析連作土壤細(xì)菌群落功能和真菌營(yíng)養(yǎng)型變化,并探究連作對(duì)根際微生物分子生態(tài)網(wǎng)絡(luò)復(fù)雜性和穩(wěn)定性的影響,最后通過計(jì)算Modified Stochasticity Ratio(校正的隨機(jī)率)來解析根際微生物群落構(gòu)建機(jī)制??傊?,研究有助于揭示西瓜連作的障礙機(jī)制,從而為西瓜的種植和產(chǎn)業(yè)可持續(xù)發(fā)展提供理論支持。

    1 材料與方法

    1.1 研究地點(diǎn)概況

    在南京市蔬菜花卉科學(xué)研究所(31°43′N,118°47′E)建立了一個(gè)長(zhǎng)期試驗(yàn)大棚。試驗(yàn)地土壤為黃棕壤,年平均日照時(shí)數(shù)2 016.4 h,年均氣溫14.7℃。

    1.2 試驗(yàn)設(shè)計(jì)和土壤樣品采集

    本試驗(yàn)在同一地塊種植品種為8424西瓜,種植面積為28 m2,共種植32株,種植間距為0.45 m,其生長(zhǎng)期間日均溫度38℃,夜均溫度24℃。自2014年,每年進(jìn)行春(3月)、秋(8月)兩茬種植。本試驗(yàn)共設(shè)置3個(gè)連作時(shí)間點(diǎn)進(jìn)行采樣,分別于2014(西瓜未連作,對(duì)照組CK)、2015(連作2茬,CC2)、2017年(連作6茬,CC6)的5月采集西瓜根際土壤樣品。在每一個(gè)時(shí)間點(diǎn)隨機(jī)采取9株長(zhǎng)勢(shì)相似西瓜的根際土,其中每3株合并為一個(gè)混合樣本,即每個(gè)時(shí)間點(diǎn)有3個(gè)重復(fù)樣本。共獲得9個(gè)樣本(3個(gè)時(shí)間點(diǎn)×3次重復(fù))。對(duì)于根際土壤的采集,將根部完整取出,并輕微搖動(dòng)去除松散的可見土壤[22]。將根樣品切成約2 cm小段,裝入50 mL無菌無酶離心管中,加入適量無菌蒸餾水,置于搖床中以180 r/min轉(zhuǎn)速振蕩30 min,洗滌緊貼根部表面的土壤,之后將離心管中的根取出,并以8 000 r/min離心10 min得到根際土壤樣品[23],將收集的根際土壤保存在-80℃冰箱中進(jìn)行后續(xù)分析。

    1.3 DNA提取和Illumina MiSeq測(cè)序

    使用FastDNA?(MP Biomedicals, Cleveland, OH, USA)土壤DNA提取試劑盒從0.5 g濕重土壤中提取根際土壤微生物總DNA。之后使用Qubit?DNA濃度檢測(cè)試劑盒測(cè)定提取的DNA質(zhì)量。分別使用引物515F/907R和ITS1F/ITS2R引物擴(kuò)增細(xì)菌基因V4—V5區(qū)和真菌基因ITS1區(qū)。使用Illumina Miseq系統(tǒng)進(jìn)行測(cè)序。

    1.4 數(shù)據(jù)統(tǒng)計(jì)與分析

    使用MOTHUR軟件對(duì)下機(jī)數(shù)據(jù)進(jìn)行質(zhì)量控制[24]。具體而言,將短序列以及低質(zhì)量序列進(jìn)行去除,通過數(shù)據(jù)庫(kù)比對(duì),去除存在的嵌合體序列。將剩余的高質(zhì)量序列以97%的閾值生成可操作分類單元(operational taxonomic units,OTUs)并選取代表序列進(jìn)行分類學(xué)分析,再基于樣本最小序列數(shù)對(duì)所有樣本抽平后計(jì)算chao1、richness和shannon多樣性指數(shù)。基于 Bray-Curtis 距離用“R(v4.0.2)”中的“vegan”包進(jìn)行主坐標(biāo)分析(PCoA)。使用“pheatmap”和“ggplot2”包在R環(huán)境中揭示不同土壤樣品中細(xì)菌和真菌物種在門水平和屬水平上的差異。采用PICRUSt2軟件對(duì)細(xì)菌OTU豐度表進(jìn)行標(biāo)準(zhǔn)化后進(jìn)行功能基因預(yù)測(cè)分析,將OTU信息與KEGG(Kyoto Encyclo-pediaof Genesand Genomes)數(shù)據(jù)庫(kù)進(jìn)行比對(duì),獲得三級(jí)功能基因的豐度,參考KEGG的二級(jí)代謝通路的功能基因類別劃分歸類,獲得功能基因組成。為了確定真菌營(yíng)養(yǎng)模式分類群,使用FUNGuild進(jìn)行注釋[25],并保留probable和highly probable置信度的數(shù)據(jù),以及single mode和mixed mode的數(shù)據(jù)。

    使用隨機(jī)矩陣?yán)碚摚≧MT)自動(dòng)確定合適的網(wǎng)絡(luò)閾值[26],基于所有處理的細(xì)菌和真菌總OTU豐度矩陣進(jìn)行物種共發(fā)生網(wǎng)絡(luò)分析并使用Gephi0.9.2得到細(xì)菌總網(wǎng)絡(luò)和真菌總網(wǎng)絡(luò)。之后從總節(jié)點(diǎn)和邊文件中抽提出每個(gè)處理所擁有的節(jié)點(diǎn)和邊,得到每個(gè)處理的細(xì)菌子網(wǎng)絡(luò)和真菌子網(wǎng)絡(luò)。通過計(jì)算網(wǎng)絡(luò)各節(jié)點(diǎn)的模塊內(nèi)連通性()和模塊間連通性(),識(shí)別網(wǎng)絡(luò)中潛在關(guān)鍵物種[27]。通過計(jì)算和值的大小,將網(wǎng)絡(luò)關(guān)鍵物種劃分為三類:模塊中心(模塊內(nèi)高度連接的節(jié)點(diǎn),>2.5和<0.62)、網(wǎng)絡(luò)中心(整個(gè)網(wǎng)絡(luò)內(nèi)高度連接的節(jié)點(diǎn),>2.5和>0.62)和網(wǎng)絡(luò)連接節(jié)點(diǎn)(連接不同模塊的節(jié)點(diǎn),<2.5和>0.62)[26]。使用Gephi0.9.2繪制網(wǎng)絡(luò)圖。使用自然連通度指標(biāo)確定網(wǎng)絡(luò)結(jié)構(gòu)的穩(wěn)定性[28],通過去除網(wǎng)絡(luò)中的節(jié)點(diǎn)來模擬受到的外界擾動(dòng),并計(jì)算節(jié)點(diǎn)去除后網(wǎng)絡(luò)的穩(wěn)定性,以評(píng)估網(wǎng)絡(luò)在受到隨機(jī)攻擊后穩(wěn)定性下降的速度。最后,使用校正的隨機(jī)率(modified stochasticity ratio,MST)來量化確定性和隨機(jī)過程在群落組裝過程中的相對(duì)重要性,以50%作為更確定性(MST<50%)和更隨機(jī)(MST>50%)集合之間的分界點(diǎn)[16]。

    2 結(jié)果

    2.1 連作對(duì)西瓜根際細(xì)菌和真菌群落多樣性和組成的影響

    不同連作茬次下西瓜根際細(xì)菌群落的chao1、richness和shannon指數(shù)隨著連作茬次的延長(zhǎng)呈現(xiàn)先下降后上升的趨勢(shì)(圖1-A),而真菌群落的相應(yīng)指數(shù)顯著降低(圖1-B)。這表明,連作在一定程度上有利于細(xì)菌群落豐富度的增加,而降低了真菌群落豐富度?;贐ray-Curtis距離的主坐標(biāo)分析(PCoA)對(duì)對(duì)照、連作2茬、連作6茬的根際細(xì)菌和真菌群落組成特征進(jìn)行可視化,結(jié)果表明,不同連作茬次中根際細(xì)菌和真菌群落聚集成不同的簇,也說明了隨著連作茬次的增加,微生物群落組成發(fā)生顯著改變(圖2-A,Adonis2=0.71,=0.003;圖2-B,Adonis2=0.84,=0.004)。

    我們將平均相對(duì)豐度大于1%的細(xì)菌門和真菌門劃分為優(yōu)勢(shì)菌門,3種不同連作茬次的根際土壤中細(xì)菌優(yōu)勢(shì)門分別是變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)、厚壁菌門(Firmicutes)、放線菌門(Actinobacteria)(圖2-C)。在細(xì)菌群落中,變形菌門在3種不同連作茬次中土壤的比例均最高,占比為83.3%—47.1%,并且隨著連作茬次增加呈先減少后增加趨勢(shì)。相比于對(duì)照,擬桿菌門和厚壁菌門的相對(duì)豐度在連作2茬中增加了22.2%和11.2%,而在連作6茬中分別增加了5.1%和0.76%(圖2-C)。屬于優(yōu)勢(shì)真菌門的是子囊菌門(Ascomycota)、擔(dān)子菌門(Basidiomycota)、壺菌門(Chytridiomycota)(圖2-D)。真菌群落中,子囊菌門和擔(dān)子菌門不僅相對(duì)豐度最高,且隨連作茬次表現(xiàn)出不同的差異變化。在真菌群落中,子囊菌門在3個(gè)不同連作茬次土壤中的占比均最高,為32.16%—66.1%,并且其相對(duì)豐度隨連作茬次增加呈先減少后增加趨勢(shì)。而擔(dān)子菌門的相對(duì)豐度呈先增加后減少趨勢(shì)(圖2-D)。在屬水平,相比于對(duì)照,連作2茬和6茬細(xì)菌群落中鞘氨醇單胞菌屬()相對(duì)豐度分別減少了4.5%和0.6%,鞘氨醇菌屬()相對(duì)豐度分別減少了7.6%和7.2%,溶桿菌屬()相對(duì)豐度分別減少了2.5%和1.6%;連作6茬與對(duì)照相比,劍菌屬()和芽孢桿菌屬()相對(duì)豐度分別減少了1.1%和0.2%(圖2-E)。真菌群落中,相比于對(duì)照,連作2茬和6茬的鐮刀菌屬()相對(duì)豐度分別增加了4.1%和60.6%(圖2-F)。這些結(jié)果表明連作過程明顯改變了根際細(xì)菌、真菌群落結(jié)構(gòu),且群落組成相似度隨連作茬次的增加而降低。

    A:不同連作茬次下根際細(xì)菌群落chao1、richness和shannon指數(shù);B:不同連作茬次下根際真菌群落chao1、richness和shannon指數(shù)。不同的字母表示不同處理間細(xì)菌、真菌豐度的顯著差異(P<0.05)

    2.2 連作對(duì)西瓜根際細(xì)菌和真菌群落預(yù)測(cè)功能的影響

    采用PICRUSt2對(duì)各個(gè)連作茬次中根際細(xì)菌群落進(jìn)行功能預(yù)測(cè)。在三級(jí)功能水平的預(yù)測(cè)功能中,對(duì)照的西瓜根際細(xì)菌群落共得到297個(gè)代謝通路,相較于其他兩個(gè)時(shí)期,富集了50個(gè)代謝通路;連作2茬共得到295個(gè)代謝通路,相較于其他兩個(gè)時(shí)期,富集了41個(gè)代謝通路;連作6茬共得到298個(gè)代謝通路,相較于其他兩個(gè)時(shí)期,富集了26個(gè)代謝通路(圖3-a)。將所獲得相對(duì)豐度高的代謝通路進(jìn)行二級(jí)功能水平歸類,結(jié)果顯示,除其他次生代謝物的生物合成外,其余功能基因的相對(duì)豐度隨連作茬次變化均有不同的響應(yīng)(圖3-c)。相比于對(duì)照,連作2茬會(huì)顯著減少碳水化合物、萜類和多酮類化合物的代謝、細(xì)胞生長(zhǎng)和死亡等代謝途徑相對(duì)豐度,而增加了復(fù)制和修復(fù)、輔助因子和維生素的代謝、運(yùn)輸和分解代謝等途徑;連作6茬會(huì)顯著減少有害物質(zhì)的生物降解和代謝及其他氨基酸的代謝等途徑的相對(duì)豐度,反而增加了多糖的生物合成和代謝途徑(圖3-c)。真菌功能群使用FUNGuild進(jìn)行預(yù)測(cè),其中連作6茬會(huì)顯著富集病理型(pathotroph)真菌(圖3-b)。其中相比于對(duì)照,連作6茬會(huì)顯著富Plant pathogen(植物病原體)類型真菌,而Soil saprotroph(土壤腐生生物)會(huì)顯著減少(圖3-d)。

    A:不同連作茬次下根際細(xì)菌群落主坐標(biāo)分析(PCoA);B:不同連作茬次下根際真菌群落主坐標(biāo)分析(PCoA);C:不同連作茬次下根際細(xì)菌群落門水平組成;D:不同連作茬次下根際真菌群落門水平組成;E:不同連作茬次下根際細(xì)菌群落屬水平組成;F:不同連作茬次下根際真菌群落屬水平組成

    a:不同連作茬次中KEGG三級(jí)代謝通路相對(duì)豐度差異變化。括號(hào)中的數(shù)字分別代表著比其他兩個(gè)處理顯著富集的代謝通路數(shù)量和預(yù)測(cè)得到的代謝通路總數(shù);b:不同連作茬次中根際真菌營(yíng)養(yǎng)型相對(duì)豐度的差異分析,括號(hào)中的數(shù)字代表著比其他兩個(gè)處理顯著富集的真菌營(yíng)養(yǎng)型數(shù)量;c:餅圖表示將圖a中富集的三級(jí)代謝通路歸類于二級(jí)代謝通路;d:不同連作茬次中根際真菌營(yíng)養(yǎng)亞型相對(duì)豐度的差異分析

    2.3 連作對(duì)西瓜根際細(xì)菌和真菌群落共生網(wǎng)絡(luò)結(jié)構(gòu)和關(guān)鍵物種的影響

    基于OTU豐度矩陣進(jìn)行物種共發(fā)生網(wǎng)絡(luò)分析并得到細(xì)菌總網(wǎng)絡(luò)和子網(wǎng)絡(luò)以及真菌總網(wǎng)絡(luò)和子網(wǎng)絡(luò)(圖4-a)(見方法2.4)。結(jié)果顯示,連作2茬與連作6茬的根際細(xì)菌和真菌網(wǎng)絡(luò)比對(duì)照具有更多的網(wǎng)絡(luò)節(jié)點(diǎn)、平均度和種間相互作用,而模塊化指數(shù)隨著連作茬次增加而下降(表1)。變形菌門在3個(gè)連作茬次(對(duì)照組、連作2茬和連作6茬,分別占62.6%、63.4%和65.9%)的網(wǎng)絡(luò)互作關(guān)系中占主導(dǎo)地位,與其他物種相互聯(lián)系的數(shù)量最多(圖4-a)。這表明與對(duì)照相比,連作增加了網(wǎng)絡(luò)中變形菌門物種的比例。在真菌群落網(wǎng)絡(luò)中,子囊菌門在3個(gè)連作茬次中(對(duì)照組、連作2茬和連作6茬,分別占73.2%、70.7%和69.6%)的網(wǎng)絡(luò)互作關(guān)系中占主導(dǎo)地位,與其他物種相互聯(lián)系的數(shù)量最多(圖4-d)。這表明與對(duì)照相比,連作減少了網(wǎng)絡(luò)中子囊菌門物種的比例。此外,對(duì)照、連作2茬與連作6茬的真菌共生網(wǎng)絡(luò)中正相關(guān)連接數(shù)占比分別為41.33%、51.48%和55.30%(表1),表明在連作過程中根際真菌間趨于競(jìng)爭(zhēng)關(guān)系而不是合作關(guān)系。

    表1 不同連作茬次下西瓜根際細(xì)菌和真菌網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)

    網(wǎng)絡(luò)結(jié)構(gòu)穩(wěn)定性的結(jié)果顯示,在細(xì)菌群落中,雖然連作6茬的根際初始網(wǎng)絡(luò)自然連通性最高(0.0078),其次是連作2茬(0.0075)和對(duì)照網(wǎng)絡(luò)(0.0068),但是當(dāng)去除20%的網(wǎng)絡(luò)節(jié)點(diǎn)時(shí),連作6茬的自然連通性急劇降至0.0028;當(dāng)去除40%的網(wǎng)絡(luò)節(jié)點(diǎn)時(shí),連作2茬的自然連通性急劇降至0.0024(圖4-b)。隨著節(jié)點(diǎn)不斷被移除,對(duì)照網(wǎng)絡(luò)的連通度下降趨勢(shì)(Slope=-0.004)比連作2茬(Slope=-0.005)和連作6茬(Slope=-0.006)更加平穩(wěn)(圖4-b)。而真菌群落中,連作6茬的根際初始網(wǎng)絡(luò)自然連通性仍是最高(0.085),其次是對(duì)照組(0.023)和連作2茬(0.013)網(wǎng)絡(luò),但是當(dāng)去除40%的網(wǎng)絡(luò)節(jié)點(diǎn)時(shí),連作6茬的自然連通性急劇降至0.0071,而連作2茬網(wǎng)絡(luò)隨著節(jié)點(diǎn)不斷被移除,其自然連通性下降趨勢(shì)(Slope=-0.011)比連作2茬(Slope=-0.021)和連作6茬(Slope=-0.077)更加平穩(wěn)(圖4-e);對(duì)照網(wǎng)絡(luò)在去除20%的網(wǎng)絡(luò)節(jié)點(diǎn)后,其自然連通性穩(wěn)定在一定水平(圖4-e)。這些結(jié)果表明在移除網(wǎng)絡(luò)節(jié)點(diǎn)后,對(duì)照的根際微生物網(wǎng)絡(luò)結(jié)構(gòu)仍保持相對(duì)較高的穩(wěn)定性。

    對(duì)照的細(xì)菌群落共擁有9個(gè)關(guān)鍵物種(圖4-c),網(wǎng)絡(luò)模塊中心(module hubs)主要來源于變形菌門和厚壁菌門,其中主要來自于假單胞菌屬()、芽孢桿菌屬()和鞘氨醇單胞菌屬()等。對(duì)照的真菌群落共擁有1個(gè)關(guān)鍵物種(圖4-f),網(wǎng)絡(luò)模塊中心(module hubs)主要來源于子囊菌門的閉小囊菌屬()。

    2.4 連作對(duì)西瓜根際細(xì)菌和真菌群落構(gòu)建過程的影響

    使用校正的隨機(jī)率(MST)定量生態(tài)過程中的隨機(jī)性和確定性。結(jié)果表明,在細(xì)菌群落中,隨機(jī)性過程在對(duì)照和連作6茬的細(xì)菌群落內(nèi)占主導(dǎo),而在連作2茬的細(xì)菌群落中確定性過程占主導(dǎo)(圖5-A);在真菌群落中,確定性過程在連作2茬和連作6茬中占主導(dǎo)(圖5-B)。

    A:不同連作茬次下根際細(xì)菌群落的MST柱狀圖;B:不同連作茬次下根際真菌群落的MST柱狀圖。不同的字母表示不同處理間的顯著差異(P<0.05)

    3 討論

    3.1 西瓜連作對(duì)根際土壤細(xì)菌和真菌群落組成和功能的影響

    在本研究中,利用Illumina Miseq高通量技術(shù)分析了連作后西瓜根際土壤細(xì)菌和真菌群落結(jié)構(gòu)。我們發(fā)現(xiàn)隨著連作茬次的增加,細(xì)菌和真菌群落特征差異較大。這可能是由于在連作系統(tǒng)中,相同的田間管理模式、植物化感物質(zhì)的長(zhǎng)期累積以及同類型根系分泌物長(zhǎng)期釋放的生態(tài)效應(yīng),會(huì)使微生物區(qū)系發(fā)生顯著改變[13,29]。具體來說,相比于對(duì)照,西瓜根際細(xì)菌群落的richness和shannon指數(shù)在連作2茬和6茬中表現(xiàn)出先下降后上升的趨勢(shì),而chao1指數(shù)沒有顯著的變化,這與部分前人的研究結(jié)果一致[30-31]。多樣性指數(shù)出現(xiàn)下降的原因可能是在2茬的連作過程中肥料的添加會(huì)暫時(shí)的增加營(yíng)養(yǎng)物質(zhì),土壤中的部分物種對(duì)營(yíng)養(yǎng)物質(zhì)的添加反應(yīng)迅速,進(jìn)而通過快速繁殖占領(lǐng)寡營(yíng)養(yǎng)微生物的生態(tài)位,導(dǎo)致微生物多樣性立即下降[32]。而連作6茬過程中,隨著養(yǎng)分長(zhǎng)期大量的積累,包括肥料的連續(xù)性投入和植物分泌物和殘?bào)w遺留,土壤中生態(tài)位會(huì)在一定程度上擴(kuò)張,會(huì)有更多的有機(jī)物質(zhì)和養(yǎng)分易被微生物利用,從而促進(jìn)了細(xì)菌的生長(zhǎng),有利于根際土壤中細(xì)菌多樣性的提高。有研究表明,變形菌門和酸桿菌門在養(yǎng)分循環(huán)中起重要作用,并且也與抑制病害等功能相關(guān)[33-35]。其中鞘氨醇單胞菌屬(變形菌門)具有降解多環(huán)芳烴和降低纖維素的能力[36],溶桿菌屬(變形菌門)對(duì)真菌、革蘭氏陰性細(xì)菌、革蘭氏陽(yáng)性細(xì)菌和線蟲都有突出的拮抗作用[37]。而本研究結(jié)果中,相對(duì)于對(duì)照,連作后的變形菌門和酸桿菌門的相對(duì)豐度都呈現(xiàn)降低的趨勢(shì),而且鞘氨醇單胞菌屬和芽孢桿菌屬等的相對(duì)豐度降低,從而導(dǎo)致植物抵御病原菌感染的能力下降,并且這種降低也會(huì)增加發(fā)生土壤病害的風(fēng)險(xiǎn)。而連作后的真菌群落相比于對(duì)照,chao1、richness和shannon指數(shù)顯著降低,這與吳鳳芝等的研究結(jié)果相一致[38]。有研究指出,微生物多樣性的減少是發(fā)生土壤病害的重要原因[5],而在本研究中真菌群落的多樣性指數(shù)在連作后顯著降低,可能增加發(fā)生土壤病害的風(fēng)險(xiǎn)。在真菌群落中,有研究表明鐮刀菌屬(子囊菌門)為致病菌,會(huì)引起西瓜、黃瓜枯萎病的發(fā)生[39-40],而本研究結(jié)果顯示出鐮刀菌屬的相對(duì)豐度在連作6茬后顯著增加。在CHEN等[41]的研究中得出,真菌病原體的增加和有益真菌的減少可能是導(dǎo)致連作植物生長(zhǎng)和產(chǎn)量下降的重要因素。因此,連作會(huì)改變微生物群落多樣性和組成,破壞根際生態(tài)平衡,從而更加容易引起連作障礙。

    我們利用PICRUSt2對(duì)不同連作茬次下根際細(xì)菌群落進(jìn)行功能預(yù)測(cè),發(fā)現(xiàn)根際細(xì)菌群落積極參與了基本的代謝過程。其中碳水化合物代謝途徑與固氮、溶磷等作用有著密切的關(guān)聯(lián),而且該途徑的產(chǎn)物為其他物質(zhì)的合成提供了原料[42];有害物質(zhì)的生物降解和代謝途徑包括了苯甲酸降解、多環(huán)芳烴降解和二甲苯降解途徑,這些途徑與植物抗逆能力密切相關(guān);氨基酸代謝與脫氨和轉(zhuǎn)氨等作用息息相關(guān),并且主要參與了植物氮素循環(huán)[43]。而本研究的結(jié)果表明,連作會(huì)減少根際細(xì)菌對(duì)碳源的利用、有害物質(zhì)的降解及次生代謝物的合成代謝途徑。尤其是相比于對(duì)照,連作6茬后的根際細(xì)菌氨基酸代謝相對(duì)豐度顯著降低,這可能表明連作后根際環(huán)境不利于細(xì)菌的繁殖、代謝和合成物質(zhì),并且土壤養(yǎng)分循環(huán)可能受到一定程度的抑制。因此,在連作下,一些比較重要的代謝功能(如,碳水化合物代謝途徑,有害物質(zhì)的降解和次生代謝物的合成代謝途徑)的減少或者喪失會(huì)是形成連作障礙的主要原因之一[44]。此外,我們利用FUNGuild來確定真菌分類群的營(yíng)養(yǎng)模式。結(jié)果顯示,相比于對(duì)照,連作6茬會(huì)顯著增加病原真菌,這種真菌會(huì)通過損害宿主細(xì)胞而獲取營(yíng)養(yǎng)[45],并且營(yíng)養(yǎng)亞型中植物病原體(如鐮刀菌屬)的相對(duì)豐度也會(huì)顯著增加,這些也可能是引起連作障礙的主要原因之一。

    3.2 西瓜連作對(duì)根際土壤細(xì)菌和真菌群落網(wǎng)絡(luò)結(jié)構(gòu)的影響

    物種間相互作用是影響土壤病害的發(fā)生一種重要的生物因素,在群落構(gòu)建過程中發(fā)揮重要作用[46],與關(guān)鍵生態(tài)過程和宿主健康密切相關(guān)[47]。本研究結(jié)果表明,多茬次的連作使細(xì)菌和真菌群落產(chǎn)生更復(fù)雜的潛在相互作用,這可能是因?yàn)橹参餁報(bào)w量的增加,以及前茬作物產(chǎn)生的各種分泌物和遺傳效應(yīng)影響土壤微生物群落[48]。對(duì)照的網(wǎng)絡(luò)直徑和模塊化均大于連作2茬和連作6茬,并且有研究指出具有更高模塊化的微生物網(wǎng)絡(luò)結(jié)構(gòu)更能抵抗環(huán)境干擾[44,49]。因此,連作會(huì)導(dǎo)致根際微生物對(duì)環(huán)境變化的抵抗力減弱,容易引起連作障礙的發(fā)生。另外,本研究結(jié)果表明,在去除一定比例的網(wǎng)絡(luò)節(jié)點(diǎn)后,連作6茬的網(wǎng)絡(luò)穩(wěn)定性指數(shù)急劇下降,表明連作使西瓜根際網(wǎng)絡(luò)結(jié)構(gòu)更加脆弱。因此,連作時(shí)長(zhǎng)可能是影響根際網(wǎng)絡(luò)穩(wěn)定性的關(guān)鍵因素。

    網(wǎng)絡(luò)中潛在的關(guān)鍵物種在維持多物種共存和生態(tài)系統(tǒng)功能方面可能至關(guān)重要[50-51]。我們觀察到,在對(duì)照的西瓜根際微生物網(wǎng)絡(luò)中,關(guān)鍵類群有芽孢桿菌屬()、鞘氨醇菌屬()、劍菌屬()和鞘氨醇單胞菌()等菌屬。據(jù)報(bào)道芽孢桿菌屬具有分解土壤中有機(jī)物的能力,如纖維素和木質(zhì)素[52];鞘氨醇單胞菌屬的成員具有降解復(fù)雜有機(jī)物質(zhì)、拮抗植物病原真菌以及分泌胞外多糖等多種生態(tài)功能[53];鞘氨醇菌屬能降解多環(huán)芳烴或酚類有機(jī)物[54],適應(yīng)環(huán)境能力強(qiáng)[55];劍菌屬能夠作為有益菌接種至土壤中促進(jìn)植物的生長(zhǎng)[56],有助于維持植物激素平衡、控制根系發(fā)育、促進(jìn)營(yíng)養(yǎng)獲取以及預(yù)防宿主植物疾病[57]。而在連作6茬后這些關(guān)鍵類群的相對(duì)豐度相比于對(duì)照發(fā)生了明顯的減少,因此連作土壤病害的發(fā)生可能是由于關(guān)鍵類群的減少或者缺失。

    3.3 西瓜連作對(duì)根際土壤細(xì)菌和真菌群落組裝的影響

    為了更好地了解微生物群落組成如何響應(yīng)環(huán)境變化,我們探究了西瓜根際微生物群落組裝的生態(tài)過程。在本研究中,當(dāng)短期連作時(shí)(連作2茬),一些非生物環(huán)境因素(如pH、鹽度、有機(jī)碳含量)的選擇(環(huán)境過濾和物種間相互的拮抗與協(xié)同作用)會(huì)使細(xì)菌群落構(gòu)建過程從隨機(jī)性轉(zhuǎn)變?yōu)榇_定性[20];但在較長(zhǎng)連作茬次下,細(xì)菌種群更容易受到隨機(jī)性事件的影響,這可以用體型-可塑性假說來解釋[58],即較小的生物更可能受到擴(kuò)散限制而不是物種分類的影響,因?yàn)檩^小的生物(如,細(xì)菌類群)具有更大的代謝可塑性和更大的環(huán)境耐受性,因此可以生活在不同的棲息地[59];相對(duì)于真菌的生活史策略而言,細(xì)菌類群的生存策略更加傾向于r策略,且在長(zhǎng)期連作下,養(yǎng)分較為充足,細(xì)菌的生長(zhǎng)沒有受到環(huán)境的限制,趨向于隨機(jī)性過程。而在真菌群落中,隨機(jī)性過程在對(duì)照中更為突出,但連作后逐漸轉(zhuǎn)變?yōu)榇_定性過程。這可能是一些環(huán)境變量對(duì)根際微生物群落的組裝的潛在影響,如植物分泌物和殘?bào)w的增加、土壤理化性質(zhì)的變化等。環(huán)境因素的差異也可能隨著連作茬次的增加而增加,并在根際微生物群落構(gòu)建中發(fā)揮關(guān)鍵作用[60]。值得注意的是,當(dāng)?shù)丨h(huán)境的變化逐漸增加了真菌群落的確定性選擇的重要性,而細(xì)菌群落的隨機(jī)性卻占主導(dǎo),這可能是因?yàn)殚L(zhǎng)期連作引起環(huán)境變化是由真菌性病害引起的(即由尖孢鐮刀菌引起的西瓜枯萎?。?,植物可能會(huì)向環(huán)境釋放特定的次生物質(zhì)來招募更多的病原性真菌[12-14],真菌群落因此會(huì)受到更強(qiáng)的環(huán)境選擇。同時(shí)細(xì)菌和真菌群落的MST值也呈現(xiàn)先降低后增大的趨勢(shì),表明隨機(jī)性和確定性過程共同調(diào)控著連作環(huán)境下根際微生物群落的構(gòu)建。

    4 結(jié)論

    多茬次的連作顯著影響了西瓜根際細(xì)菌和真菌群落結(jié)構(gòu),包括有益菌屬的相對(duì)豐度減少和致病菌屬的相對(duì)豐度增加,也會(huì)導(dǎo)致微生物群落網(wǎng)絡(luò)的相互作用更加復(fù)雜,但群落穩(wěn)定性和對(duì)環(huán)境變化的抵抗性減弱。此外,細(xì)菌群落功能組成發(fā)生顯著變化,包括有害物質(zhì)的生物降解和代謝,氨基酸的代謝等功能基因相對(duì)豐度在長(zhǎng)期連作后降低;而真菌群落中病原真菌的相對(duì)豐度在連作6茬后顯著增多。多茬次的連作也會(huì)使微生物群落組裝的機(jī)制發(fā)生變化,隨機(jī)性和確定性過程共同調(diào)控著連作環(huán)境下根際微生物群落的構(gòu)建。因此,在實(shí)際生產(chǎn)中,應(yīng)注意通過輪作或施用微生物肥料來維持西瓜根際土壤微生態(tài)的穩(wěn)定,以緩解連作障礙。本研究結(jié)果為理解連作與微生物群落結(jié)構(gòu)差異的關(guān)系提供了新視角,并揭示了長(zhǎng)期連作下微生物群落構(gòu)建過程和潛在機(jī)制,同時(shí)為開發(fā)微生物肥料進(jìn)行土壤改良,進(jìn)而緩解連作障礙提供理論依據(jù)。

    [1] 田晴, 高丹美, 李慧, 劉守偉, 周新剛, 吳鳳芝. 小麥根系分泌物對(duì)西瓜連作土壤真菌群落結(jié)構(gòu)的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2020, 53(5): 1018-1028. doi:10.3864/j.issn.0578-1752.2020.05.013.

    TIAN Q, GAO D M, LI H, LIU S W, ZHOU X G, WU F Z. Effects of wheat root exudates on the structure of fungi community in continuous cropping watermelon soil. Scientia Agricultura Sinica, 2020, 53(5): 1018-1028. doi:10.3864/j.issn.0578-1752.2020.05.013. (in Chinese)

    [2] 耿士均, 劉刊, 商海燕, 權(quán)俊嬌, 陸小平, 王波. 園藝作物連作障礙的研究進(jìn)展. 北方園藝, 2012(7): 190-195.

    GENG S J, LIU K, SHANG H Y, QUAN J J, LU X P, WANG B. Research progress of continuous cropping obstacle in horticultural plants. Northern Horticulture, 2012(7): 190-195. (in Chinese)

    [3] MARTYN R. Fusarium wilt of watermelon: 120 years of research. Horticultural Reviews, 2014, 42: 349-442.

    [4] 滕凱, 陳前鋒, 周志成, 向青松, 張敏, 尹華群, 劉勇軍. 煙草連作障礙與土壤理化性質(zhì)及微生物多樣性特征的關(guān)聯(lián). 微生物學(xué)通報(bào), 2020, 47(9): 2848-2856.

    TENG K, CHEN Q F, ZHOU Z C, XIANG Q S, ZHANG M, YIN H Q, LIU Y J. Effect of soil physical and chemical properties and microbial community on continuous cropping obstacles in tobacco field. Microbiology China, 2020, 47(9): 2848-2856. (in Chinese)

    [5] YUAN J, WEN T, ZHANG H, ZHAO M L, PENTON C R, THOMASHOW L S, SHEN Q R. Predicting disease occurrence with high accuracy based on soil macroecological patterns of Fusarium wilt. The ISME Journal, 2020, 14(12): 2936-2950.

    [6] DESSAUX Y, GRANDCLéMENT C, FAURE D. Engineering the rhizosphere. Trends in Plant Science, 2016, 21(3): 266-278.

    [7] BORDENSTEIN S R, THEIS K R. Host biology in light of the microbiome: Ten principles of holobionts and hologenomes. PLoS Biology, 2015, 13(8): e1002226.

    [8] 張樹生, 楊興明, 茆澤圣, 黃啟為, 徐陽(yáng)春, 沈其榮. 連作土滅菌對(duì)黃瓜(Cucumis sativus)生長(zhǎng)和土壤微生物區(qū)系的影響. 生態(tài)學(xué)報(bào), 2007, 27(5): 1809-1817.

    ZHANG S S, YANG X M, MAO Z S, HUANG Q W, XU Y C, SHEN Q R. Effects of sterilization on growth of cucumber plants and soil microflora in a continuous mono-cropping soil. Acta Ecologica Sinica, 2007, 27(5): 1809-1817. (in Chinese)

    [9] BRADFORD M A, MCCULLEY R L, CROWTHER T W, OLDFIELD E E, WOOD S A, FIERER N. Cross-biome patterns in soil microbial respiration predictable from evolutionary theory on thermal adaptation. Nature Ecology & Evolution, 2019, 3(2): 223-231.

    [10] BAHRAM M, HILDEBRAND F, FORSLUND S K, ANDERSON J L, SOUDZILOVSKAIA N A, BODEGOM P M, BENGTSSON-PALME J, ANSLAN S, COELHO L P, HAREND H, HUERTA-CEPAS J, MEDEMA M H, MALTZ M R, MUNDRA S, OLSSON P A, PENT M, P?LME S, SUNAGAWA S, RYBERG M, TEDERSOO L, BORK P. Structure and function of the global topsoil microbiome. Nature, 2018, 560(7717): 233-237.

    [11] GAO Z Y, HAN M K, HU Y Y, LI Z Q, LIU C F, WANG X, TIAN Q, JIAO W J, HU J M, LIU L F, GUAN Z J, MA Z M. Effects of continuous cropping of sweet potato on the fungal community structure in rhizospheric soil. Frontiers in Microbiology, 2019, 10: 2269.

    [12] HU L F, ROBERT C A M, CADOT S, ZHANG X, YE M, LI B B, MANZO D, CHERVET N, STEINGER T, VAN DER HEIJDEN M G A, SCHLAEPPI K, ERB M. Root exudate metabolites drive plant-soil feedbacks on growth and defense by shaping the rhizosphere microbiota. Nature Communications, 2018, 9: 2738.

    [13] LI X G, DING C F, HUA K, ZHANG T L, ZHANG Y N, ZHAO L, YANG Y R, LIU J G, WANG X X. Soil sickness of peanuts is attributable to modifications in soil microbes induced by peanut root exudates rather than to direct allelopathy. Soil Biology and Biochemistry, 2014, 78: 149-159.

    [14] XU W H, LIU D, WU F Z, LIU S W. Root exudates of wheat are involved in suppression of Fusarium wilt in watermelon in watermelon-wheat companion cropping. European Journal of Plant Pathology, 2015, 141(1): 209-216.

    [15] 曾維愛, 楊昭玥, 黃洋, 谷亞冰, 陶界錳, 劉勇軍, 謝鵬飛, 蔡海林, 尹華群. 長(zhǎng)期連作農(nóng)田土壤細(xì)菌群落結(jié)構(gòu)和共現(xiàn)網(wǎng)絡(luò)拓?fù)湫再|(zhì)對(duì)土壤理化性質(zhì)的響應(yīng). 微生物學(xué)報(bào), 2022, 62(6): 2403-2416.

    ZENG W A, YANG Z Y, HUANG Y, GU Y B, TAO J M, LIU Y J, XIE P F, CAI H L, YIN H Q. Response of soil bacterial community structure and co-occurrence network topology properties to soil physicochemical properties in long-term continuous cropping farmland. Acta Microbiologica Sinica, 2022, 62(6): 2403-2416. (in Chinese)

    [16] NING D L, DENG Y, TIEDJE J M, ZHOU J Z. A general framework for quantitatively assessing ecological stochasticity. Proceedings of the National Academy of Sciences of the United States of America, 2019, 116(34): 16892-16898.

    [17] GRYTA A, FR?C M, OSZUST K. The application of the biolog EcoPlate approach in ecotoxicological evaluation of dairy sewage sludge. Applied Biochemistry and Biotechnology, 2014, 174(4): 1434-1443.

    [18] GUO J J, LIU W B, ZHU C, LUO G W, KONG Y L, LING N, WANG M, DAI J Y, SHEN Q R, GUO S W. Bacterial rather than fungal community composition is associated with microbial activities and nutrient-use efficiencies in a paddy soil with short-term organic amendments. Plant and Soil, 2018, 424(1): 335-349.

    [19] GUO H P, DONG P S, GAO F, HUANG L, WANG S P, WANG R Y, YAN M C, ZHANG D M. Sucrose addition directionally enhances bacterial community convergence and network stability of the shrimp culture system. NPJ Biofilms and Microbiomes, 2022, 8: 22.

    [20] CHEN Q L, HU H W, YAN Z Z, LI C Y, NGUYEN B A T, SUN A Q, ZHU Y G, HE J Z. Deterministic selection dominates microbial community assembly in termite mounds. Soil Biology and Biochemistry, 2021, 152: 108073.

    [21] BURNS A R, STEPHENS W Z, STAGAMAN K, WONG S, RAWLS J F, GUILLEMIN K, BOHANNAN B J. Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. The ISME Journal, 2016, 10(3): 655-664.

    [22] TOJU H, OKAYASU K, NOTAGUCHI M. Leaf-associated microbiomes of grafted tomato plants. Scientific Reports, 2019, 9: 1787.

    [23] LING N, SONG Y, RAZA W, HUANG Q W, GUO S W, SHEN Q R. The response of root-associated bacterial community to the grafting of watermelon. Plant and Soil, 2015, 391(1): 253-264.

    [24] SCHLOSS P D, WESTCOTT S L, RYABIN T, HALL J R, HARTMANN M, HOLLISTER E B, LESNIEWSKI R A, OAKLEY B B, PARKS D H, ROBINSON C J, SAHL J W, STRES B, THALLINGER G G, VAN HORN D J, WEBER C F. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology, 2009, 75(23): 7537-7541.

    [25] NGUYEN N H, SONG Z W, BATES S T, BRANCO S, TEDERSOO L, MENKE J, SCHILLING J S, KENNEDY P G. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecology, 2016, 20: 241-248.

    [26] DENG Y, JIANG Y H, YANG Y F, HE Z L, LUO F, ZHOU J Z. Molecular ecological network analyses. BMC Bioinformatics, 2012, 13: 113.

    [27] SHI S J, NUCCIO E E, SHI Z J, HE Z L, ZHOU J Z, FIRESTONE M K. The interconnected rhizosphere: high network complexity dominates rhizosphere assemblages. Ecology Letters, 2016, 19(8): 926-936.

    [28] PENG G S, WU J. Optimal network topology for structural robustness based on natural connectivity. Physica A: Statistical Mechanics and Its Applications, 2016, 443: 212-220.

    [29] LI X G, JOUSSET A, DE BOER W, CARRIóN V J, ZHANG T L, WANG X X, KURAMAE E E. Legacy of land use history determines reprogramming of plant physiology by soil microbiome. The ISME Journal, 2019, 13(3): 738-751.

    [30] ZHU S Y, WANG Y Z, XU X M, LIU T M, WU D Q, ZHENG X, TANG S W, DAI Q Z. Potential use of high-throughput sequencing of soil microbial communities for estimating the adverse effects of continuous cropping on ramie (L. Gaud). PLoS ONE, 2018, 13(5): e0197095.

    [31] 李晶晶, 續(xù)勇波. 連作年限對(duì)設(shè)施百合土壤微生物多樣性的影響. 土壤通報(bào), 2020, 51(2): 343-351.

    LI J J, XU Y B. Effects of continuous cropping years of lily on soil microbial diversities under greenhouse cultivation. Chinese Journal of Soil Science, 2020, 51(2): 343-351. (in Chinese)

    [32] VAN BRUGGEN A H C, SEMENOV A M. In search of biological indicators for soil health and disease suppression. Applied Soil Ecology, 2000, 15(1): 13-24.

    [33] ZHANG Y, CAO C Y, PENG M, XU X J, ZHANG P, YU Q J, SUN T. Diversity of nitrogen-fixing, ammonia-oxidizing, and denitrifying bacteria in biological soil crusts of a revegetation area in Horqin Sandy Land, Northeast China. Ecological Engineering, 2014, 71: 71-79.

    [34] MUNOZ R, TEELING H, AMANN R, ROSSELLó-MóRA R. Ancestry and adaptive radiation of Bacteroidetes as assessed by comparative genomics. Systematic and Applied Microbiology, 2020, 43(2): 126065.

    [35] KWAK M J, KONG H G, CHOI K, KWON S K, SONG J Y, LEE J, LEE P A, CHOI S Y, SEO M, LEE H J, JUNG E J, PARK H, ROY N, KIM H, LEE M M, RUBIN E M, LEE S W, KIM J F. Rhizosphere microbiome structure alters to enable wilt resistance in tomato. Nature Biotechnology, 2018, 36(11): 1100-1109.

    [36] 王濤, 藍(lán)慧, 田云, 盧向陽(yáng). 多環(huán)芳烴的微生物降解機(jī)制研究進(jìn)展. 化學(xué)與生物工程, 2016, 33(2): 8-14.

    WANG T, LAN H, TIAN Y, LU X Y. Research progress on microbial degradation mechanisms for polycyclic aromatic hydrocarbons. Chemistry & Bioengineering, 2016, 33(2): 8-14. (in Chinese)

    [37] 姜英華, 胡白石, 劉鳳權(quán). 植物土傳病原菌拮抗細(xì)菌的篩選與鑒定. 中國(guó)生物防治, 2005, 21(4): 260-264.

    JIANG Y H, HU B S, LIU F Q. Selection and identification of antagonistic bacteria against soil-borne plant pathogens. Chinese Journal of Biological Control, 2005, 21(4): 260-264. (in Chinese)

    [38] 吳鳳芝, 王學(xué)征. 設(shè)施黃瓜連作和輪作中土壤微生物群落多樣性的變化及其與產(chǎn)量品質(zhì)的關(guān)系. 中國(guó)農(nóng)業(yè)科學(xué), 2007, 40(10): 2274-2280. doi:10.3321/j.issn: 0578-1752.2007.10.021.

    WU F Z, WANG X Z. Effect of monocropping and rotation on soil microbial community diversity and cucumber yield, quality under protected cultivation. Scientia Agricultura Sinica, 2007, 40(10): 2274-2280. doi:10.3321/j.issn: 0578-1752.2007.10.021. (in Chinese)

    [39] PLIEGO C, RAMOS C, DE VICENTE A, CAZORLA F M. Screening for candidate bacterial biocontrol agents against soilborne fungal plant pathogens. Plant and Soil, 2011, 340(1): 505-520.

    [40] ?í?KOVá D, ?R?TKA P, KOLA?íK M, KUBáTOVá A, PA?OUTOVá S. Assessing the pathogenic effect of Fusarium, Geosmithia and Ophiostoma fungi from broad-leaved trees. Folia Microbiologica, 2005, 50(1): 59-62.

    [41] CHEN M N, LI X, YANG Q L, CHI X Y, PAN L J, CHEN N, YANG Z, WANG T, WANG M, YU S L. Soil eukaryotic microorganism succession as affected by continuous cropping of peanut: Pathogenic and beneficial fungi were selected. PLoS ONE, 2012, 7(7): e40659.

    [42] 徐偉, 葛陽(yáng)陽(yáng), 陳翠婷, 馬婷婷. 基于宏基因組技術(shù)分析傳統(tǒng)紅茶菌中菌群組成及其主要代謝通路. 食品工業(yè)科技, 2018, 39(5): 119-123, 129.

    XU W, GE Y Y, CHEN C T, MA T T. Microorganism composition and main metabolic pathways analysis of traditional Kombucha by metagenomic technology. Science and Technology of Food Industry, 2018, 39(5): 119-123, 129. (in Chinese)

    [43] 丁鈺珮, 杜宇佳, 高廣磊, 張英, 曹紅雨, 朱賓賓, 楊思遠(yuǎn), 張儆醒, 邱業(yè), 劉惠林. 呼倫貝爾沙地樟子松人工林土壤細(xì)菌群落結(jié)構(gòu)與功能預(yù)測(cè). 生態(tài)學(xué)報(bào), 2021, 41(10): 4131-4139.

    DING Y P, DU Y J, GAO G L, ZHANG Y, CAO H Y, ZHU B B, YANG S Y, ZHANG J X, QIU Y, LIU H L. Soil bacterial community structure and functional prediction ofvar. mongolica plantations in the Hulun Buir Sandy Land. Acta Ecologica Sinica, 2021, 41(10): 4131-4139. (in Chinese)

    [44] 李曉雪. 辣椒對(duì)土壤連作障礙的響應(yīng)及其成因研究[D]. 保定: 河北農(nóng)業(yè)大學(xué), 2021.

    LI X X. The mechanisms of continuous cropping obstacle formation of chilli pepper and the response of chilli pepper to continuously mono-cropped system[D]. Baoding: Hebei Agricultural University, 2021. (in Chinese)

    [45] TOJU H, KISHIDA O, KATAYAMA N, TAKAGI K. Networks depicting the fine-scale co-occurrences of fungi in soil horizons. PLoS ONE, 2016, 11(11): e0165987.

    [46] KONOPKA A, LINDEMANN S, FREDRICKSON J. Dynamics in microbial communities: Unraveling mechanisms to identify principles. The ISME Journal, 2015, 9(7): 1488-1495.

    [47] LAYEGHIFARD M, HWANG D M, GUTTMAN D S. Disentangling interactions in the microbiome: A network perspective. Trends in Microbiology, 2017, 25(3): 217-228.

    [48] NUCCIO E E, STARR E, KARAOZ U, BRODIE E L, ZHOU J Z, TRINGE S G, MALMSTROM R R, WOYKE T, BANFIELD J F, FIRESTONE M K, PETT-RIDGE J. Niche differentiation is spatially and temporally regulated in the rhizosphere. The ISME Journal, 2020, 14(4): 999-1014.

    [49] YUAN M M, GUO X, WU L W, ZHANG Y, XIAO N J, NING D L, SHI Z, ZHOU X S, WU L Y, YANG Y F, TIEDJE J M, ZHOU J Z. Climate warming enhances microbial network complexity and stability. Nature Climate Change, 2021, 11(4): 343-348.

    [50] HERREN C M, MCMAHON K D. Keystone taxa predict compositional change in microbial communities. Environmental Microbiology, 2018, 20(6): 2207-2217.

    [51] CHEN L J, JIANG Y J, LIANG C, LUO Y, XU Q S, HAN C, ZHAO Q G, SUN B. Competitive interaction with keystone taxa induced negative priming under biochar amendments. Microbiome, 2019, 7(1): 77.

    [52] SUN X L, XU Z H, XIE J Y, HESSELBERG-THOMSEN V, TAN T M, ZHENG D Y, STRUBE M L, DRAGO? A, SHEN Q R, ZHANG R F, KOVáCS á T.stimulates resident rhizospherefor plant health through metabolic interactions. The ISME Journal, 2022, 16(3): 774-787.

    [53] WHITE D C, SUTTON S D, RINGELBERG D B. The genus Sphingomonas: Physiology and ecology. Current Opinion in Biotechnology, 1996, 7(3): 301-306.

    [54] 閻潔, 余雪巍, 李鑒博, 顧海萍, 郭二輝. 一株菲降解細(xì)菌產(chǎn)生生物表面活性劑特性的研究. 生態(tài)環(huán)境學(xué)報(bào), 2021, 30(8): 1683-1694.

    YAN J, YU X W, LI J B, GU H P, GUO E H. Research on the characterization of surfactant produced by a phenanthrene-degrading strain. Ecology and Environmental Sciences, 2021, 30(8): 1683-1694. (in Chinese)

    [55] ZHOU X Y, LIANG Y, REN G F, ZHENG K W, WU Y, ZENG X Y, ZHONG Y, YU Z Q, PENG P A. Biotransformation of tris(2- chloroethyl) phosphate (TCEP) in sediment microcosms and the adaptation of microbial communities to TCEP. Environmental Science & Technology, 2020, 54(9): 5489-5497.

    [56] SCHMITZ L, YAN Z C, SCHNEIJDERBERG M, DE ROIJ M, PIJNENBURG R, ZHENG Q, FRANKEN C, DECHESNE A, TRINDADE L M, VAN VELZEN R, BISSELING T, GEURTS R, CHENG X. Synthetic bacterial community derived from a desert rhizosphere confers salt stress resilience to tomato in the presence of a soil microbiome. The ISME Journal, 2022, 16(8): 1907-1920.

    [57] XU J, ZHANG Y Z, ZHANG P F, TRIVEDI P, RIERA N, WANG Y Y, LIU X, FAN G Y, TANG J L, COLETTA-FILHO H D, CUBERO J, DENG X L, ANCONA V, LU Z J, ZHONG B L, ROPER M C, CAPOTE N, CATARA V, PIETERSEN G, VERNIèRE C, AL-SADI A M, LI L, YANG F, XU X, WANG J, YANG H M, JIN T, WANG N. The structure and function of the global citrus rhizosphere microbiome. Nature Communications, 2018, 9: 4894.

    [58] ISABWE A, YAO H F, ZHANG S X, JIANG Y J, BREED M F, SUN X. Spatial assortment of soil organisms supports the size-plasticity hypothesis. ISME Communications, 2022, 2: 102.

    [59] WU W X, LU H P, SASTRI A, YEH Y C, GONG G C, CHOU W C, HSIEH C H. Contrasting the relative importance of species sorting and dispersal limitation in shaping marine bacterial versus protist communities. The ISME Journal, 2018, 12(2): 485-494.

    [60] DELGADO-BAQUERIZO M, MAESTRE F T, REICH P B, JEFFRIES T C, GAITAN J J, ENCINAR D, BERDUGO M, CAMPBELL C D, SINGH B K. Microbial diversity drives multifunctionality in terrestrial ecosystems. Nature Communications, 2016, 7: 10541.

    Characteristics and Succession of Rhizosphere Soil Microbial Communities in Continuous Cropping Watermelon

    GUO HanYue1, WANG DongSheng2, RUAN Yang1, QIAO YiZhu1, ZHANG YunTao1, LI Ling1, HUANG QiWei1, GUO ShiWei1, LING Ning1, SHEN QiRong1

    1College of Resources and Environment Science, Nanjing Agricultural University/Key Laboratory of Organic-based Fertilizers of China/Jiangsu Provincial Key Laboratory for Solid Organic Waste Utilization, Nanjing 210095;2Nanjing Institute of Vegetable Science, Nanjing 210042

    【Objective】The aim of this study was to investigate the effects of continuous cropping on the construction and potential functions of bacterial and fungal communities in the rhizosphere soil of watermelon, and to clarify the adaptability of rhizosphere microorganisms to environmental changes, so as to provide a theoretical basis for ecological control of watermelon continuous cropping obstacles and healthy maintenance of farmland. 【Method】In this study, the rhizosphere soil of watermelon without continuous cropping (CK), continuous cropping for 2 times and continuous cropping for 6 times was used as the research object. 16S rRNA and ITS high-throughput sequencing were used to analyze the effects of continuous cropping on the bacterial and fungal communities of rhizosphere soil of watermelon. 【Result】With the increasing continuous cropping times, the bacterial diversity index in the rhizosphere soil of watermelon showed a trend of first decreasing and then increasing, while the fungal diversity index decreased significantly. At the bacterial genus level, multiple times of continuous cropping decreased the relative abundance ofandin watermelon rhizosphere soil; at the fungal genus level, the relative abundance ofincreased with the continuous cropping times. Compared with CK, the network complexity of continuous cropping was higher, but the stability of network structure was lower. In addition, compared with CK, the relative abundance of biodegradation pathways of harmful substances and metabolic pathways of amino acids could be significantly reduced after continuous cropping for 6 times; the relative abundance of pathogenic fungi significantly increased in the fungal community after continuous cropping for 6 times. During community succession, the stochastic processes dominated rhizosphere bacterial community construction in watermelon under continuous cropping, while the deterministic processes dominated rhizosphere fungal community construction in watermelon. 【Conclusion】Continuous cropping caused changes in community characteristics, functional composition and succession process of rhizosphere bacteria and fungi. The decrease of key functions of bacterial community, the increase of pathologic fungi and the decrease of stability of microbial community network might be the important factors leading to occurrence of watermelon continuous cropping obstacles.

    watermelon; continuous cropping soil; rhizosphere microbial community; function prediction; community succession

    2022-11-19;

    2023-01-09

    國(guó)家自然科學(xué)基金(31972506)

    郭晗玥,E-mail:2020103110@stu.njau.edu.cn。通信作者凌寧,E-mail:nling@njau.edu.cn

    (責(zé)任編輯 李云霞)

    猜你喜歡
    菌門根際群落
    大學(xué)生牙齦炎齦上菌斑的微生物群落
    根際微生物對(duì)植物與土壤交互調(diào)控的研究進(jìn)展
    合成微生物群落在發(fā)酵食品中的應(yīng)用研究
    野生樹鼩與人工飼養(yǎng)樹鼩消化道不同部位微生物組成的比較研究
    饑餓與重?cái)z食對(duì)河蟹腸道菌群結(jié)構(gòu)的影響
    昆蟲體內(nèi)微生物多樣性的影響因素研究進(jìn)展
    妊娠期糖尿病腸道菌群變化及臨床價(jià)值分析
    黃花蒿葉水提物對(duì)三七根際尖孢鐮刀菌生長(zhǎng)的抑制作用
    促植物生長(zhǎng)根際細(xì)菌HG28-5對(duì)黃瓜苗期生長(zhǎng)及根際土壤微生態(tài)的影響
    春季和夏季巢湖浮游生物群落組成及其動(dòng)態(tài)分析
    亚洲人成电影观看| 飞空精品影院首页| 飞空精品影院首页| 性高湖久久久久久久久免费观看| 亚洲国产欧美一区二区综合| 各种免费的搞黄视频| 国产成人精品久久二区二区免费| 妹子高潮喷水视频| 精品少妇黑人巨大在线播放| 丝袜人妻中文字幕| www.熟女人妻精品国产| 亚洲三区欧美一区| 日本精品一区二区三区蜜桃| 免费高清在线观看视频在线观看| 青青草视频在线视频观看| 亚洲五月色婷婷综合| 亚洲成人免费av在线播放| www.av在线官网国产| 男女边摸边吃奶| 大码成人一级视频| 老司机影院成人| tube8黄色片| 国产男女超爽视频在线观看| av片东京热男人的天堂| 夜夜骑夜夜射夜夜干| 97在线人人人人妻| 天天添夜夜摸| 国产97色在线日韩免费| 国产亚洲精品一区二区www | 少妇猛男粗大的猛烈进出视频| 侵犯人妻中文字幕一二三四区| 日韩视频在线欧美| 日本a在线网址| 成人国语在线视频| 又黄又粗又硬又大视频| 日韩视频在线欧美| 免费看十八禁软件| 国产精品一区二区在线不卡| 97精品久久久久久久久久精品| 久热这里只有精品99| 国产免费福利视频在线观看| 欧美一级毛片孕妇| 亚洲专区国产一区二区| 精品乱码久久久久久99久播| 色婷婷久久久亚洲欧美| 欧美在线黄色| 视频在线观看一区二区三区| 91老司机精品| 欧美+亚洲+日韩+国产| av有码第一页| 久久精品久久久久久噜噜老黄| 啦啦啦视频在线资源免费观看| 黄色 视频免费看| 亚洲第一青青草原| 日韩有码中文字幕| 9191精品国产免费久久| 十八禁高潮呻吟视频| 国产男女超爽视频在线观看| 中文精品一卡2卡3卡4更新| 无限看片的www在线观看| 丝袜喷水一区| 一级片免费观看大全| 亚洲精品乱久久久久久| 夜夜骑夜夜射夜夜干| 国产一区有黄有色的免费视频| 最近中文字幕2019免费版| 99久久精品国产亚洲精品| 在线观看舔阴道视频| av在线app专区| 成在线人永久免费视频| 国产男人的电影天堂91| 亚洲精品第二区| a级片在线免费高清观看视频| 黄片大片在线免费观看| 丁香六月欧美| 无遮挡黄片免费观看| 免费观看av网站的网址| av一本久久久久| 午夜两性在线视频| 国产精品九九99| 久久久精品国产亚洲av高清涩受| 无限看片的www在线观看| 国产精品麻豆人妻色哟哟久久| 国产伦人伦偷精品视频| av天堂久久9| 国产精品一二三区在线看| 男女床上黄色一级片免费看| 免费观看a级毛片全部| 窝窝影院91人妻| 久久人人97超碰香蕉20202| 亚洲人成77777在线视频| 极品人妻少妇av视频| 一区二区三区乱码不卡18| 波多野结衣av一区二区av| 狠狠狠狠99中文字幕| 亚洲国产中文字幕在线视频| 精品一品国产午夜福利视频| 免费观看av网站的网址| 日本wwww免费看| 麻豆av在线久日| 亚洲精品美女久久久久99蜜臀| 考比视频在线观看| 亚洲av电影在线进入| 18在线观看网站| 中文字幕色久视频| 亚洲国产欧美网| 男人爽女人下面视频在线观看| 人妻人人澡人人爽人人| 咕卡用的链子| 动漫黄色视频在线观看| 一边摸一边做爽爽视频免费| 国产精品熟女久久久久浪| 精品久久久精品久久久| 少妇猛男粗大的猛烈进出视频| 国产日韩一区二区三区精品不卡| 免费久久久久久久精品成人欧美视频| e午夜精品久久久久久久| 国产成人av教育| 国产成人免费观看mmmm| 欧美大码av| 美女国产高潮福利片在线看| 91字幕亚洲| 国产精品香港三级国产av潘金莲| 日本精品一区二区三区蜜桃| 免费观看av网站的网址| 久久久久国内视频| 两性夫妻黄色片| 欧美黄色淫秽网站| 亚洲精品国产色婷婷电影| 五月天丁香电影| 老熟女久久久| 亚洲男人天堂网一区| 90打野战视频偷拍视频| 狠狠婷婷综合久久久久久88av| 色精品久久人妻99蜜桃| 欧美人与性动交α欧美精品济南到| 少妇的丰满在线观看| 国产麻豆69| 欧美黄色淫秽网站| 日韩中文字幕视频在线看片| 最新在线观看一区二区三区| 色婷婷av一区二区三区视频| 菩萨蛮人人尽说江南好唐韦庄| 人人妻人人澡人人看| 亚洲av男天堂| svipshipincom国产片| 国产激情久久老熟女| 国产亚洲午夜精品一区二区久久| 黑人猛操日本美女一级片| av一本久久久久| 久久午夜综合久久蜜桃| 丝袜美足系列| 两性午夜刺激爽爽歪歪视频在线观看 | 纵有疾风起免费观看全集完整版| 五月开心婷婷网| 亚洲中文av在线| 女人久久www免费人成看片| 天天添夜夜摸| 久久午夜综合久久蜜桃| av一本久久久久| 亚洲欧洲精品一区二区精品久久久| 欧美另类一区| 正在播放国产对白刺激| 国产人伦9x9x在线观看| 久久影院123| 国产精品熟女久久久久浪| 精品亚洲成a人片在线观看| 在线天堂中文资源库| 51午夜福利影视在线观看| 欧美性长视频在线观看| 国产亚洲一区二区精品| 欧美在线黄色| 一级,二级,三级黄色视频| 法律面前人人平等表现在哪些方面 | 久久 成人 亚洲| www日本在线高清视频| 亚洲精品国产av蜜桃| 亚洲欧洲精品一区二区精品久久久| 亚洲美女黄色视频免费看| 香蕉国产在线看| 99精品久久久久人妻精品| 成人国产av品久久久| 国产在线视频一区二区| kizo精华| 18禁国产床啪视频网站| 久久女婷五月综合色啪小说| 亚洲人成电影观看| 一本色道久久久久久精品综合| 精品国产乱子伦一区二区三区 | 国产精品1区2区在线观看. | 在线观看www视频免费| 国产精品久久久久成人av| 一区二区三区四区激情视频| 国产成人免费无遮挡视频| 久久国产精品人妻蜜桃| 日日夜夜操网爽| 777久久人妻少妇嫩草av网站| 欧美精品av麻豆av| 久久亚洲国产成人精品v| 日韩欧美一区视频在线观看| 亚洲精品中文字幕一二三四区 | 国产精品亚洲av一区麻豆| 美国免费a级毛片| 91字幕亚洲| 美女高潮到喷水免费观看| 国产精品.久久久| 91精品国产国语对白视频| 美女国产高潮福利片在线看| 欧美+亚洲+日韩+国产| 国产伦理片在线播放av一区| 男女之事视频高清在线观看| 男女之事视频高清在线观看| 桃红色精品国产亚洲av| av福利片在线| 最新在线观看一区二区三区| 亚洲欧洲精品一区二区精品久久久| 国精品久久久久久国模美| 欧美人与性动交α欧美软件| 热re99久久精品国产66热6| 国产高清视频在线播放一区 | 天天添夜夜摸| 国产av又大| 高清av免费在线| 丝袜美足系列| 国产精品国产三级国产专区5o| 亚洲熟女毛片儿| 十八禁高潮呻吟视频| 国产欧美日韩精品亚洲av| 亚洲精品美女久久av网站| 久久精品人人爽人人爽视色| 亚洲精品久久久久久婷婷小说| 中文精品一卡2卡3卡4更新| 在线观看免费高清a一片| 美女国产高潮福利片在线看| 天天躁夜夜躁狠狠躁躁| 国产又色又爽无遮挡免| 日韩三级视频一区二区三区| 99国产综合亚洲精品| 国产精品一区二区在线观看99| 国产福利在线免费观看视频| 成人手机av| 少妇人妻久久综合中文| 国产精品免费大片| 看免费av毛片| 久久精品亚洲熟妇少妇任你| 性色av一级| 日本91视频免费播放| 精品乱码久久久久久99久播| 黑丝袜美女国产一区| 国产精品影院久久| 国产精品麻豆人妻色哟哟久久| 亚洲国产精品999| 中亚洲国语对白在线视频| 成人手机av| 一本色道久久久久久精品综合| 欧美乱码精品一区二区三区| 自线自在国产av| 亚洲精品粉嫩美女一区| 91av网站免费观看| 亚洲五月婷婷丁香| 久久久国产成人免费| 日韩一卡2卡3卡4卡2021年| 黄片大片在线免费观看| 麻豆国产av国片精品| a在线观看视频网站| 色婷婷av一区二区三区视频| 男女高潮啪啪啪动态图| 精品国产国语对白av| 韩国高清视频一区二区三区| 午夜福利视频在线观看免费| 十八禁网站网址无遮挡| 欧美精品一区二区大全| 三级毛片av免费| 91老司机精品| 亚洲中文字幕日韩| 精品国产一区二区三区四区第35| 母亲3免费完整高清在线观看| 中文字幕人妻熟女乱码| av网站在线播放免费| 脱女人内裤的视频| 9热在线视频观看99| 三上悠亚av全集在线观看| 日韩欧美一区视频在线观看| 成人免费观看视频高清| 亚洲激情五月婷婷啪啪| 日韩一卡2卡3卡4卡2021年| 成人18禁高潮啪啪吃奶动态图| 日韩免费高清中文字幕av| 性色av一级| 国产真人三级小视频在线观看| 亚洲精品粉嫩美女一区| 色精品久久人妻99蜜桃| 操美女的视频在线观看| 淫妇啪啪啪对白视频 | 精品少妇一区二区三区视频日本电影| 九色亚洲精品在线播放| 午夜影院在线不卡| 叶爱在线成人免费视频播放| 又大又爽又粗| 国产激情久久老熟女| 精品福利观看| 欧美黑人精品巨大| 亚洲国产欧美一区二区综合| 爱豆传媒免费全集在线观看| 日日夜夜操网爽| 国产视频一区二区在线看| 一区二区日韩欧美中文字幕| 亚洲国产欧美一区二区综合| 国产成人免费观看mmmm| 色视频在线一区二区三区| 丝袜脚勾引网站| 亚洲成人免费av在线播放| 黄色视频不卡| 汤姆久久久久久久影院中文字幕| 中文字幕色久视频| 丝袜美足系列| 美女视频免费永久观看网站| 久久毛片免费看一区二区三区| 俄罗斯特黄特色一大片| 99精品久久久久人妻精品| 亚洲av电影在线观看一区二区三区| 欧美人与性动交α欧美软件| 成人国产一区最新在线观看| 黄片小视频在线播放| 亚洲av欧美aⅴ国产| 国产精品亚洲av一区麻豆| 91精品伊人久久大香线蕉| 另类精品久久| kizo精华| a级毛片在线看网站| 少妇的丰满在线观看| 精品卡一卡二卡四卡免费| 黄色 视频免费看| 精品乱码久久久久久99久播| a级毛片在线看网站| 色综合欧美亚洲国产小说| 精品少妇久久久久久888优播| 日韩一区二区三区影片| 欧美黑人欧美精品刺激| 后天国语完整版免费观看| 美女福利国产在线| 久久精品国产a三级三级三级| 91成年电影在线观看| 久久99一区二区三区| 国产精品亚洲av一区麻豆| 欧美乱码精品一区二区三区| 欧美精品高潮呻吟av久久| 国产成人精品无人区| 中文精品一卡2卡3卡4更新| 黄色毛片三级朝国网站| 亚洲欧美日韩高清在线视频 | 黑人操中国人逼视频| 天天躁夜夜躁狠狠躁躁| 51午夜福利影视在线观看| 精品国产乱子伦一区二区三区 | 交换朋友夫妻互换小说| 91麻豆av在线| www日本在线高清视频| 亚洲国产毛片av蜜桃av| 91av网站免费观看| 99热网站在线观看| 国产在视频线精品| av天堂在线播放| 91av网站免费观看| 亚洲成国产人片在线观看| 老司机深夜福利视频在线观看 | 亚洲黑人精品在线| 我的亚洲天堂| 国产日韩欧美在线精品| 亚洲熟女精品中文字幕| 岛国在线观看网站| 国产伦理片在线播放av一区| 极品人妻少妇av视频| 丝袜脚勾引网站| 高清视频免费观看一区二区| 国产成人一区二区三区免费视频网站| 国产成+人综合+亚洲专区| 欧美老熟妇乱子伦牲交| 18禁黄网站禁片午夜丰满| av电影中文网址| 欧美+亚洲+日韩+国产| 欧美精品一区二区大全| 久久久久久久久免费视频了| 嫩草影视91久久| 69av精品久久久久久 | 国产亚洲欧美精品永久| 亚洲精品成人av观看孕妇| 午夜福利免费观看在线| av又黄又爽大尺度在线免费看| 国产精品免费视频内射| 久久久国产精品麻豆| 欧美另类一区| 亚洲成人免费av在线播放| 日韩人妻精品一区2区三区| 人人妻人人添人人爽欧美一区卜| 少妇精品久久久久久久| 亚洲精品日韩在线中文字幕| 大片免费播放器 马上看| 亚洲,欧美精品.| 最新在线观看一区二区三区| 老熟妇乱子伦视频在线观看 | 最近中文字幕2019免费版| 国产一区二区在线观看av| 一本久久精品| 精品视频人人做人人爽| 久久精品成人免费网站| 免费av中文字幕在线| 久久国产精品影院| 国产亚洲欧美在线一区二区| 亚洲精品国产一区二区精华液| 国产精品久久久久久精品古装| 欧美日韩亚洲综合一区二区三区_| 狂野欧美激情性xxxx| 岛国在线观看网站| 国产主播在线观看一区二区| 一级片免费观看大全| 中文字幕最新亚洲高清| 国产精品久久久久久精品古装| www.av在线官网国产| 性少妇av在线| 欧美激情久久久久久爽电影 | 色综合欧美亚洲国产小说| 亚洲欧美精品综合一区二区三区| 久久天堂一区二区三区四区| 久久九九热精品免费| 啦啦啦视频在线资源免费观看| 国产在线一区二区三区精| 天堂中文最新版在线下载| 午夜成年电影在线免费观看| 国产精品久久久久久精品古装| 俄罗斯特黄特色一大片| 中文字幕色久视频| 亚洲成人免费电影在线观看| 亚洲性夜色夜夜综合| 国产男女内射视频| 欧美 日韩 精品 国产| 久久人妻熟女aⅴ| 久久青草综合色| 一区福利在线观看| 桃红色精品国产亚洲av| 制服人妻中文乱码| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美精品自产自拍| 色视频在线一区二区三区| 国产日韩欧美亚洲二区| 人成视频在线观看免费观看| 欧美久久黑人一区二区| 纯流量卡能插随身wifi吗| 国产av又大| 国产精品欧美亚洲77777| 99国产极品粉嫩在线观看| 欧美日韩成人在线一区二区| 亚洲欧美清纯卡通| 亚洲视频免费观看视频| 亚洲avbb在线观看| 精品亚洲成a人片在线观看| 国产av又大| 久热爱精品视频在线9| 午夜福利免费观看在线| 丰满迷人的少妇在线观看| 波多野结衣av一区二区av| 久久久久久久久久久久大奶| 日日爽夜夜爽网站| 日韩视频在线欧美| 69精品国产乱码久久久| 成年人午夜在线观看视频| av有码第一页| 国产日韩一区二区三区精品不卡| 免费黄频网站在线观看国产| 在线观看免费日韩欧美大片| 精品久久久精品久久久| 国产不卡av网站在线观看| 欧美另类亚洲清纯唯美| 午夜福利视频在线观看免费| 国产黄频视频在线观看| 欧美乱码精品一区二区三区| kizo精华| 亚洲免费av在线视频| 国产精品国产三级国产专区5o| 黑丝袜美女国产一区| 伦理电影免费视频| 午夜免费成人在线视频| 高清欧美精品videossex| 亚洲av美国av| 久久久久国产一级毛片高清牌| 999久久久精品免费观看国产| 中文字幕另类日韩欧美亚洲嫩草| 欧美久久黑人一区二区| 亚洲精品国产区一区二| 美女高潮喷水抽搐中文字幕| 日本av免费视频播放| 91精品伊人久久大香线蕉| 亚洲伊人色综图| 狠狠婷婷综合久久久久久88av| 国产精品自产拍在线观看55亚洲 | 黄色片一级片一级黄色片| 久久热在线av| 亚洲五月色婷婷综合| 欧美黄色片欧美黄色片| 一本—道久久a久久精品蜜桃钙片| 亚洲精品在线美女| 麻豆av在线久日| 欧美性长视频在线观看| av有码第一页| 黄片播放在线免费| 国产淫语在线视频| 久久精品久久久久久噜噜老黄| 天堂俺去俺来也www色官网| 国产精品久久久久久人妻精品电影 | 国产成人a∨麻豆精品| 日本91视频免费播放| 别揉我奶头~嗯~啊~动态视频 | 99国产综合亚洲精品| 国产精品一区二区免费欧美 | 黄色片一级片一级黄色片| a 毛片基地| 国产av精品麻豆| 精品一区二区三区av网在线观看 | 亚洲av成人一区二区三| 亚洲国产精品999| 国产亚洲欧美精品永久| 91成人精品电影| 99热国产这里只有精品6| 99热全是精品| 欧美精品人与动牲交sv欧美| 无限看片的www在线观看| 性高湖久久久久久久久免费观看| 久久国产精品人妻蜜桃| 在线观看免费视频网站a站| 久久久久久亚洲精品国产蜜桃av| 亚洲精品国产区一区二| 1024视频免费在线观看| 老司机影院成人| 日韩欧美免费精品| 丰满少妇做爰视频| 美女高潮到喷水免费观看| 日本av手机在线免费观看| 少妇被粗大的猛进出69影院| 一本色道久久久久久精品综合| 男人舔女人的私密视频| 欧美日韩成人在线一区二区| 欧美激情久久久久久爽电影 | 人人澡人人妻人| 国产精品久久久av美女十八| 国产成人av教育| 999久久久精品免费观看国产| 亚洲黑人精品在线| 亚洲专区字幕在线| 国产亚洲欧美在线一区二区| 国产一区二区 视频在线| 亚洲精品久久成人aⅴ小说| 日韩大片免费观看网站| 伊人久久大香线蕉亚洲五| tocl精华| 午夜两性在线视频| 亚洲国产精品一区二区三区在线| 大片免费播放器 马上看| 久久久久精品人妻al黑| 久久午夜综合久久蜜桃| 亚洲精品久久成人aⅴ小说| 国产成人欧美在线观看 | 国产一区二区三区av在线| 亚洲精品第二区| 亚洲av欧美aⅴ国产| 不卡av一区二区三区| 黑丝袜美女国产一区| 国产免费视频播放在线视频| 高清视频免费观看一区二区| 亚洲av电影在线观看一区二区三区| 一二三四社区在线视频社区8| 日日摸夜夜添夜夜添小说| 日本精品一区二区三区蜜桃| 久久青草综合色| 一二三四社区在线视频社区8| 久久ye,这里只有精品| 国产淫语在线视频| 亚洲第一欧美日韩一区二区三区 | 窝窝影院91人妻| 精品少妇黑人巨大在线播放| 天堂8中文在线网| 亚洲成av片中文字幕在线观看| 人妻久久中文字幕网| 12—13女人毛片做爰片一| 亚洲性夜色夜夜综合| www.熟女人妻精品国产| 三上悠亚av全集在线观看| 亚洲欧美精品自产自拍| 日本av手机在线免费观看| 精品亚洲乱码少妇综合久久| 十八禁人妻一区二区| 我要看黄色一级片免费的| 女性生殖器流出的白浆| 50天的宝宝边吃奶边哭怎么回事| 国产精品久久久人人做人人爽| 高清av免费在线| 国产成人系列免费观看| 我的亚洲天堂| 午夜福利一区二区在线看| 777米奇影视久久| 亚洲精品一区蜜桃| 国产精品偷伦视频观看了| 两性午夜刺激爽爽歪歪视频在线观看 | 一级黄色大片毛片| 十八禁人妻一区二区| 啦啦啦啦在线视频资源| 制服诱惑二区| 国产一级毛片在线| 久9热在线精品视频| 成年av动漫网址| 爱豆传媒免费全集在线观看| 日韩大码丰满熟妇| 老司机靠b影院| 日本wwww免费看| 欧美精品高潮呻吟av久久| 国产精品免费大片|