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

    輪作對香草蘭病原菌及根際土壤微生物群落的影響

    2025-01-12 00:00:00邢詒彰洪珊楊勁明趙青云蘇凡莊輝發(fā)王輝
    熱帶作物學報 2025年1期
    關鍵詞:宋體方正香草

    摘""要:香草蘭是重要的天然香料經濟作物,由尖孢鐮刀菌引起的土傳病害嚴重限制了其可持續(xù)發(fā)展,而作物輪作是防控土傳病害的有效途徑。本研究采用實時熒光定量PCR結合高通量測序方法,在溫室盆栽條件下,研究輪作胡椒/斑蘭葉/糯米香茶對香草蘭病原菌及根際土壤微生物群落的影響。結果表明:與撂荒和連作香草蘭處理相比,輪作斑蘭葉和糯米香茶顯著降低了土壤中尖孢鐮刀菌(Fusarium"oxysporum)的數量,并顯著提高了根際土壤細菌群落的豐富度和多樣性。撂荒、連作和輪作顯著誘導形成3種不同的群落結構,但輪作斑蘭葉和糯米香茶結構相似。輪作后門水平上顯著增加了優(yōu)勢細菌酸桿菌門(Acidobacteria)、綠彎菌門(Chloroflexi)及Rokubacteria門,以及優(yōu)勢真菌子囊菌門(Ascomycota)的相對豐度,且屬水平上增加了優(yōu)勢細菌類諾卡氏屬(Nocardioides)的相對豐度,即輪作胡椒后該屬的相對豐度為64.26%,輪作斑蘭葉后相對豐度達到98.54%,輪作糯米香茶后則為63.07%。隨機森林結果表明,輪作斑蘭葉能特異激發(fā)Terriglobus屬、克雷伯桿菌屬(Kribbella)、粘球菌屬(Myxococcus)細菌核心類群,以及枝頂孢霉屬(Acremonium)、帚枝霉屬(Sarocladium)真菌核心類群,且輪作斑蘭葉后真菌群落物種間網絡互作更緊密。輪作后顯著提高土壤pH,而土壤pH、細菌和真菌群落結構是顯著抑制病原菌的最重要指示因子。綜合表明,在長期連作的香草蘭土壤中輪作短期且經濟價值高的香料作物斑蘭葉,對緩解香草蘭土傳病害的發(fā)生效果較好。

    關鍵詞:輪作;香草蘭;斑蘭葉;根際微生物;微生物群落中圖分類號:Q939.96;S573.9""""""文獻標志碼:A

    Crop"Rotation"Effects"on"Pathogen"Dynamics"and"Rhizosphere"Microbial"Assemblages"of"Vanilla

    XING"Yizhang1,"HONG"Shan2,3,4*,"YANG"Jinming5,"ZHAO"Qingyun1,"SU"Fan1,"ZHUANG"Huifa1,"WANG"Hui1

    1."Spice"and"Beverage"of"Institute,"Chinese"Academy"of"Tropical"Agricultural"Sciences"/"Key"Laboratory"of"Genetic"Resources"Utilization"of"Spice"and"Beverage"Crops,"Ministry"of"Agriculture"and"Rural"Affairs"/"Hainan"Provincial"Key"Laboratory"of"Genetic"Improvement"and"Quality"Regulation"for"Tropical"Spice"and"Beverage"Crops,"Wanning,"Hainan"571533,"China;"2."Institute"of"Genetics"and"Developmental"Biology,"Chinese"Academy"of"Sciences"/"State"Key"Laboratory"of"Plant"Genomics,"Beijing"100101,"China;"3."Hainan"Seed"Industry"Laboratory,"Sanya,"Hainan"572025,"China;"4."Sanya"Institute,"Hainan"Academy"of"Agricultural"Sciences"/"Institute"of"Vegetables,"Hainan"Academy"of"Agricultural"Sciences,"Sanya,"Hainan"572025,"China;"5."School"of"Tropical"Agriculture"and"Forestry,"Hainan"University,"Haikou,"Hainan"570228,"China

    Abstract:"Vanilla,"a"vital"spice"crop,"faces"significant"challenges"from"soil-borne"diseases"caused"by"Fusarium"oxysporum,"affecting"its"sustainable"cultivation."This"study"utilized"quantitative"polymerase"chain"reaction"(qPCR)"and"high-throughput"sequencing"to"assess"the"impact"of"rotating"black"pepper,"pandan,"and"sweet"rice"tea"on"pathogen"levels"and"the"rhizosphere"soil"microbial"communities"of"vanilla"plants"cultivated"in"a"pot"environment."Incorporating"pandan"and"sweet"rice"tea"into"the"crop"rotation"significantly"reduced"the"prevalence"of"F."oxysporum"and"enhanced"both"the"abundance"and"diversity"of"the"rhizosphere"soil"bacteria."Three"distinct"microbial"community"structures"associated"with"fallow,"monoculture,"and"crop"rotation"conditions"were"identified."Notably,"the"rotations"involving"pandan"and"sweet"rice"tea"showed"no"significant"differences"in"the"bacterial"and"fungal"community"compositions."Crop"rotation"notably"increased"the"relative"abundance"of"key"phyla"such"as"Acidobacteria,"Chloroflexi,"Rokubacteria"and"Ascomycota."Moreover,"the"relative"abundance"of"the"dominant"bacterial"genus"Nocardioides"increased"at"the"genus"level,"with"a"relative"abundance"of"64.26%"after"pepper"rotation,"98.54%"after"pandan"rotation,"and"63.07%"after"sweet"rice"tea"rotation."The"core"microbiome,"featuring"species"such"as"Terriglobus,"Kribbella,"Myxococcus,"Acremonium"and"Sarocladium,"showed"a"particularly"strong"response"to"the"presence"of"pandan,"suggesting"a"closer"network"interaction"among"fungal"species"post-rotation."Additionally,"crop"rotation"was"found"to"significantly"raise"the"soil"pH,"which,"along"with"the"altered"bacterial"and"fungal"community"structures,"emerged"as"critical"factors"in"disease"suppression."Collectively,"our"results"suggest"that"integrating"the"economically"valuable"spice"crop"pandan"into"the"rotation"schedule"in"vanilla"monoculture"systems"can"significantly"reduce"the"incidence"of"soil-borne"diseases,"offering"a"sustainable"cultivation"strategy"for"vanilla.

    Keywords:"crop"rotation;"vanilla;"pandan;"rhizosphere"microorganisms;"microbial"community

    DOI:"10.3969/j.issn.1000-2561.2025.01.020

    "

    香草蘭(Vanilla"planifolia),隸屬于蘭科香草蘭屬,是一種多年生的熱帶藤本攀緣植物。以其獨特的香氣和風味被譽為天然食品香料之王[1],享有極高的聲譽。但因其易受土傳病原真菌尖孢鐮刀菌(Fusarium"oxysporum"f."sp."vanilla)侵染而引發(fā)嚴重的連作障礙,造成產量下降和品質降低。前期研究表明,長期連作后加劇香草蘭枯萎病的發(fā)生可歸因于土壤微生物群落組成和結構的改變,即有益微生物的減少和真菌病原菌的積累[2]。

    土壤微生物群落是增強作物抗病和抗逆能力的關鍵因素[3]。土壤微生物群落的多樣性與作物的健康狀況緊密相關[4]。研究指出,土壤中功能微生物的多樣性可以提高作物對病害的抵抗力,并通過促進營養(yǎng)吸收和產生抗菌化合物來增強作物的抗逆性[5]。土壤微生物尤其是真菌,對土壤結構的形成和維持起著重要作用[6]。它們通過分泌膠結物質促使形成土壤微團聚體,進而影響土壤的水分保持能力和通氣性,從而促進作物健康生長[7]。此外,農業(yè)管理措施包括耕作方式和作物輪作,對土壤微生物群落的組成和功能有顯著的正向調控作用[8-9]。

    作物輪作通過調控土壤微生物群落的豐富度和多樣性以保障土壤和作物健康,是環(huán)境友好且經濟可行的防控土傳病害的綠色措施[10]。基于此,團隊前期開發(fā)了胡椒(Piper"nigrum"L.-香草蘭和咖啡(Coffea"L.-香草蘭輪作模式,發(fā)現輪作胡椒能顯著增加真菌木霉屬(Trichoderma)和青霉屬(Penicillium)的有益類群。通過重塑真菌群落結構和組成,從而實現高效抗病[11]。但胡椒也是多年生香料作物,如何在有限的耕地上結合時下農業(yè)生產短平快的經濟種植效益需求,開發(fā)出更多既能緩解香草蘭連作障礙又能收獲經濟效益高的輪作作物,對香草蘭產業(yè)的可持續(xù)發(fā)展具有重要的作用。

    根際核心微生物是一類關鍵的指示性物種,它們通過與植物的互動或在微生物之間的相互作用,對微生物群落的結構進行調控[12]。團隊前期研究發(fā)現,在香草蘭長期連作系統(tǒng)中,與抑制香草蘭枯萎病較相關的是真菌群落的變化,且是以抑病核心微生物被孢霉菌屬(Mortierella)為主[13]。此外,研究團隊還發(fā)現了土壤理化性質,尤其是pH的變化對微生物群落結構和根際核心微生物的豐度有顯著影響[14]。鑒于此,將土壤理化特性與微生物群落分析相結合,對于深入理解輪作對香草蘭土傳枯萎病的抑制機制至關重要[15]。

    以往有關輪作提高香草蘭抗病機制的研究大多聚焦于多年生香料作物,且側重從土壤真菌微生物群落角度解析香草蘭抑病機制,而缺乏解析土壤細菌微生物群落及將土壤微生物和土壤理化因子聯(lián)合深度挖掘其相關機制的研究。據此,本研究在團隊前期開發(fā)胡椒的基礎上,結合新形勢下香草蘭的產業(yè)需求,篩選出熱帶地區(qū)極具特色,生產上種植周期短、經濟效益價值高的香料作物——斑蘭葉(Pandanus"amaryllifolius)和糯米香茶(Strobilanthes"tonkinensis"Lindau)構建胡椒-香草蘭/斑蘭葉-香草蘭/糯米香茶-香草蘭輪作新模式。通過熒光定量PCR和高通量測序的方法,評估這3種輪作模式對香草蘭病原菌及根際土壤微生物群落的影響,以期為香草蘭土傳病害的防控提供新視角和理論依據。"""

    1""材料與方法

    1.1""材料

    供試連作香草蘭品種為墨西哥香草蘭Vanilla"planifolia"Andrews);輪作作物選取熱帶地區(qū)特色的香料作物——胡椒/斑蘭葉/糯米香茶,均由中國熱帶農業(yè)科學院香料飲料研究所提供。

    1.2""方法

    1.2.1""試驗設計""盆栽試驗于20224月至20234月在中國熱帶農業(yè)科學院香料飲料研究所溫室大棚中開展。采集香草蘭枯萎病發(fā)病嚴重且連續(xù)種植10"a以上的香草蘭土壤。以連種香草蘭(X)為連作對照,同時設置撂荒負對照(CK),以種植胡椒(H)、斑蘭葉(B)及糯米香茶(C)為輪作模式作物。采用隨機區(qū)組設計,每個處理設立3次生物學重復,每個重復種植6盆,每盆1株。每盆填充15"kg土壤,花盆規(guī)格為32"cm×"25"cm,確保每盆植物之間有足夠的空間,并且相互獨立,避免潛在的交叉影響。

    每個花盆施加90"g普通有機肥料,肥料一次性與采集土壤混合,確保肥料在每盆土壤中分布均勻。斑蘭葉作物遵循農業(yè)生產操作,于盆栽4個月時進行1次收割,采收后繼續(xù)留苗種植。化學肥料施用方案則參照傳統(tǒng)農事操作,于盆栽12個月后采集根際土壤樣品。

    1.2.2""根際土壤樣品采集""每個處理隨機選取9盆(9個重復),將香草蘭/胡椒/斑蘭葉/糯米香茶植株從土壤中輕輕拔出,保存至4"℃冰盒中迅速帶回實驗室,輕抖植物根系,抖落下的與植物根系緊密結合的土壤視為根際土壤[16]。

    1.2.3""土壤樣品的理化性質測定""土壤理化性質測定參照《土壤農化分析》[17]中的相關方法。使用玻璃電極酸度計測定土壤懸濁液pH;采用重鉻酸鉀氧化結合外加熱方法測定土壤有機質含量;利用堿解擴散法測定土壤堿解氮含量;采用鉬銻抗比色法分析土壤速效磷含量;利用火焰光度計測定土壤速效鉀含量。

    1.2.4""樣品總DNA的提取""準確稱取0.4"g各處理的根際土壤。使用強力土壤DNA提取試劑盒(MoBio"Laboratories,"Carlsbad,"CA,"USA)提取DNA。參照試劑盒說明書進行試驗操作。提取后將所得DNA樣本存放于-70"冰箱中,備用。

    1.2.5""香草蘭根際土壤病原菌總豐度測定""采用實時熒光定量PCR技術測定香草蘭病原菌的總豐度。病原菌特異性擴增引物為AFP308RITS1F[13]。采用含有尖孢鐮刀菌(Fusarium"oxysporumITS區(qū)域序列的質粒,通過10倍稀釋法建立標準曲線。利用ABI7500實時熒光定量PCR儀,按照PCR模板程序對標準曲線和樣本進行測定。PCR擴增體系:總體積20"μL,包含10"μL"SYBR?"Premix"Ex"Taq?"2倍濃度,TaKaRa"Bio"Inc.,"Japan),每種引物0.4"μL(濃度達到10"μmol/L),0.4"μL"ROX"Reference"Dye"50倍濃度),2"μL"DNA模板,6.8"μL無菌水。通過分析溶解曲線和擴增效率來評估PCR的擴增效果。每個樣本進行3次獨立測定,所得數據進行對數轉換,以每克干土對數拷貝數呈現結果。

    1.2.6""香草蘭根際土壤樣品擴增子建庫測序""土壤細菌的16S"rRNA基因V4區(qū)域的擴增采用特異性引物520F802R[18]。對于土壤真菌,擴增ITS區(qū)域的ITS1片段,使用的引物為ITS5FITS1R[19]。擴增過程中樣本會加上Barcodes/LinkersAdapters以備測序。PCR擴增體系:25"μL"總混合體系,5倍反應緩沖液5"μL,5GC緩沖液5"μL,10"μmol/L的引物1"μL,模板DNA"2"μL,100"mmol/L"dNTP"5"μL,以及無酶水8.75"μL。PCR擴增參數設置:98"預變性,持續(xù)2"min98"變性,持續(xù)15"s;細菌擴增的退火溫度為55",持續(xù)30"s;真菌擴增的退火溫度為50",同樣持續(xù)30"s;72"延伸,30"s28~30個循環(huán)。

    PCR擴增后的產物通過QIAquick"PCR"Purification"Kit(德國QIAGEN公司)進行純化處理。純化后的DNA樣本利用Qubit?2.0"Fluorometer(美國Invitrogen公司)測定其濃度。隨后,采用NEB"Next?"UltraTM"DNA"Library"Prep"Kit"for"Illumina(英國New"England"Biolabs公司)對等濃度混合的測序樣本構建測序文庫。

    構建的文庫質量通過Agilent"2100"Bioanalyzer"Instruments(美國Agilent"Technologies"Co."Ltd)和KAPA"Library"Quantification"Kits(美國Kapa"Biosystems公司)進行檢測和驗證。從文庫構建到測序均由上海派森諾生物科技股份有限公司完成。

    1.2.7""香草蘭根際土壤樣品數據生物信息學分析""原始測序數據首先通過Trimmomatic軟件[20]v0.33)進行質量篩選,去除低質量的序列;利用Cutadapt[21]v1.9.1)識別并剪除引物序列;利用USEARCH[22]v10)對成對的序列進行合并,并使用UCHIME[23]v8.1)去除嵌合體,獲得純凈的高質量序列。

    使用QIIME2[24]v2020.6)中的DADA2[25]算法對經過質量控制的數據進行去噪處理。以測序所得序列總數的0.005%作為閾值,用于篩選和過濾掉異常序列變體。使用RDP"classifier工具對擴增的序列變體(ASV)進行序列比對分析,細菌的序列與RDP"Bacterial"16S"rRNA數據庫[26]進行匹配,真菌的序列與UNITE"Fungal"ITS數據庫[27]進行比對,確保準確識別其分類。"

    在完成測序數據的質量控制后,共鑒定出327"794個細菌的序列變體(ASVs)和12"633真菌的ASVs。Chao1ACE指數用于評估微生物群落的物種豐富度,而ShannonInvsimpson指數則用于衡量群落的多樣性。采用基于Bray-Curtis距離矩陣的主坐標分析(PCoA)和層次聚類分析測定微生物群落的Beta多樣性,并通過多元置換方差分析(PERMANOVA)來檢驗不同群落間的差異是否具有統(tǒng)計學意義。每個ASV的相對豐度以其在樣本中序列數占總序列數的百分比來表示。

    利用randomForest包分析輪作作物間差異特異富集的核心微生物。網絡分析部分使用ggCluster Net包制圖,篩選出相關系數(rgt;0.8Plt;0.05

    的數據點構建微生物間的相關性網絡圖。以上分析均在R軟件中完成。

    1.3""數據處理

    采用Duncan’s新復極差法統(tǒng)計分析不同處理間的病原菌豐度、微生物多樣性指數、微生物門和屬水平的組成,以及土壤理化因子之間的差異顯著性;利用R語言中的vegan包進行線性擬合,探究病原菌與微生物多樣性、群落結構及組成之間的關系;采用一般線性模型(逐步回歸分析)的方法評估細菌和真菌群落及其多樣性和核心物種在抑制病原菌方面的潛力。

    2""結果與分析

    2.1""輪作對香草蘭枯萎病病原菌數量的影響

    ITS測序結果發(fā)現,與撂荒(CK)和香草蘭連作(X)處理相比,輪作胡椒(H)、斑蘭葉(B)和糯米香茶(C)后根際土壤中的鐮刀菌屬(Fusarium)相對豐度均無顯著差異,但輪作斑蘭葉后呈降低趨勢,且輪作斑蘭葉顯著低于輪作糯米香茶處理(Plt;0.051A)。qPCR熒光定量結果則表明,與撂荒和香草蘭連作處理相比,輪作斑蘭葉和糯米香茶后則能顯著降低土壤中病原菌尖孢鐮刀菌(F."oxysporum)的拷貝數(Plt;0.05,1B)。以上結果表明,輪作斑蘭葉能有效降低土壤中病原菌的數量。

    "

    "

    2.2""輪作對根際土壤細菌和真菌群落多樣性的影響

    由表1可知,與香草蘭連作處理相比,輪作斑蘭葉和糯米香茶均能顯著提高土壤細菌群落的豐富度和多樣性(Plt;0.05),但二者之間無顯著差異;輪作胡椒能顯著提高土壤細菌群落的豐富度。對真菌群落而言,與撂荒和香草蘭連作處理相比,輪作斑蘭葉能顯著降低根際土壤微生物細菌群落的豐富度和多樣性(Plt;0.05);而與香草蘭連作處理相比,輪作胡椒顯著降低土壤真菌群落的豐富度。

    2.3""輪作對根際土壤細菌和真菌群落結構和組成的影響

    土壤細菌方面,通過PCoA分析發(fā)現,根際土壤細菌群落變異形成了撂荒、香草蘭連作及輪

    "

    作(胡椒/斑蘭葉/糯米香茶)3種不同典型栽培模式的結構(圖2A)。而基于置換多元方差分析(PERMANOVA)表明,撂荒、連作和輪作分組間的群落結構呈顯著差異(R2=0.34,"Plt;0.01)。進一步層次聚類分析也表明,撂荒處理和連作香草蘭處理樣本間重復性良好,均能獨立聚為一簇,胡椒輪作也能聚為一簇,但輪作斑蘭葉和輪作糯米香茶群落結構極為相似,無明顯分區(qū)(圖2B)。細菌高通量測序結果表明,與連作香草蘭相比,輪作在優(yōu)勢門水平上顯著增加酸酐菌門(Acidoba cteria)、綠彎菌門(Chloroflexi)及Rokubacteria門類群的相對豐度(Plt;0.05,圖3A)。屬水平分析則發(fā)現,與連作香草蘭相比,輪作胡椒/斑蘭葉/糯米香茶處理則不同程度上顯著增加類諾卡氏屬No cardioides)微生物,其增幅分別達到64.26%、98.54%63.07%(圖2B)。

    土壤真菌方面,通過PCoA分析表明,根際真菌土壤微生物群落也形成了撂荒、香草蘭連作及輪作(胡椒/斑蘭葉/糯米香茶)3種不同的結構(圖2C)。且經過置換多元方差檢驗也證實3不同的群落結構呈顯著差異(R2=0.32,"Plt;0.01)。層次聚類結果表明,撂荒處理和香草蘭連作處理單獨聚為一類,輪作胡椒/斑蘭葉/糯米香茶的群落結構相似,錯綜交叉聚集,無明顯區(qū)分(圖2D)。基于ITS高通量測序則發(fā)現,與連作香草蘭相比,輪作在優(yōu)勢門水平上顯著增加子囊菌門Ascomycota)類群的相對豐度(Plt;0.05,圖3B)。屬水平分析則發(fā)現,與連作香草蘭相比,輪作處理能特異顯著增加Ascobolus屬微生物的相對豐度,輪作胡椒、斑蘭葉和糯米香茶處理的Ascobolus屬相對豐度分別為4.69%、4.85%、1.75%,而連作香草蘭處理的相對豐度僅為0.01%(圖2D)。

    2.4""微生物群落多樣性、結構及組成與香草蘭枯萎病病原菌的線性關系

    線性回歸分析結果表明,土壤細菌微生物群落豐富度(Chao指數)、多樣性(Shannon指數)及群落結構(PCoA1)均與香草蘭枯萎病病原菌F."oxysporum拷貝數存在顯著負相關關系(Plt;"0.05,圖4A~4C);而土壤細菌中僅有門水平上的放線菌門(Actinobacteria)和綠彎菌門(Chloroflex)類群的相對豐度與病原菌F."oxysporum拷貝數呈顯著負相關(Plt;0.05,圖4D~4K);土壤真菌方面,真菌群落結構(PCoA1)、群落豐富度(Chao指數)及真菌組成子囊菌門(Ascomycota)類群相對豐度均與病原菌F."oxysporum拷貝數顯著負相關(Plt;0.05,圖4L~4O)。

    2.5""輪作對根際土壤核心微生物特異富集的影響

    為了識別撂荒、連作和輪作樣本間的差異關鍵核心物種,采用隨機森林分析并根據物種的貢獻度權重高低進行排序,排名前10ASV物種如圖5所示。在細菌中,與連作香草蘭處理相比,輪作斑蘭葉和糯米香茶能特異富集ASV122517、ASV8386ASV275757物種(Plt;0.05),且均與病原菌F."oxysporum顯著負相關(Plt;0.05,5A~圖5B)。而在真菌中,與連作香草蘭相比,輪作斑蘭葉能特異顯著富集ASV5998ASV7292物種(Plt;0.05),且與病原菌F."oxysporum呈顯著負相關(Plt;0.05);而輪作糯米香茶后則能顯著富集ASV5434Plt;"0.05),與病原菌F."oxysporum同樣呈顯著負相關關系(Plt;0.05,圖5C~圖5D)。

    2.6""輪作不同香料作物后根際土壤微生物網絡比較

    撂荒、連作及輪作處理的根際細菌和真菌相關性網絡具有明顯差異(圖6)。每個圓點代表1ASV物種,并被注釋到門水平;橘黃色代表物種間呈正相關,藍色代表物種間呈負相關。每個處理的網絡拓撲性質如表2所示。結果表明,在細菌中,與連作香草蘭負相關連接數(43.55%相比,輪作胡椒處理節(jié)點的負相關連接為38.89%而輪作斑蘭葉和糯米香茶的則分別為45.45%50.00%。連作與輪作之間的節(jié)點平均度差異不顯著,但輪作斑蘭葉和糯米香茶處理的平均聚類系數高于連作處理,特別是輪作斑蘭葉后平均聚類系數達0.50,而連作香草蘭的僅為0.19(圖6A)。

    在真菌中(表2),與連作香草蘭相比,輪作胡椒/斑蘭葉/糯米香茶均不同程度地減少了節(jié)點負相關的連接數,增加了平均度和平均聚類系數。與連作香草蘭負相關連接數(9.89%)相比,輪作胡椒處理的負相關連接數為2.68%,輪作斑蘭葉的僅為1.23%,而輪作糯米香茶的為3.47%,表明輪作后根際土壤真菌群落結構更穩(wěn)定。與連作香草蘭節(jié)點平均度相比,輪作胡椒/斑蘭葉/糯米香茶的節(jié)點平均度均不同程度的增加,特別是輪作斑蘭葉后的真菌網絡節(jié)點平均度達到5.74,增加了2倍,菌群連接更緊密。此外,真菌類群在微生物網絡中的節(jié)點平均度普遍高于細菌類群,特別是子囊菌門(Ascomycota)類群(圖6B),表明子囊菌門在微生物互作中具有更重要的生態(tài)位。

    綜合分析表明,在長期連作的香草蘭土壤上輪作胡椒/斑蘭葉/糯米香茶均會不同程度地增加細菌和真菌群落物種間互惠共生的正向作用,且輪作斑蘭葉效果更好。"

    2.7""抑制香草蘭病原菌數量降低的土壤理化因子和微生物指示因子預測

    通過測定輪作不同香料作物后的土壤理化性質,結果表明,與連作香草蘭處理相比,輪作胡椒/斑蘭葉/糯米香茶均能顯著提高土壤pH,以及土壤有機質、堿解氮、速效磷和速效鉀的含量(Plt;0.05),但輪作斑蘭葉和糯米香茶后的pH和速效鉀含量無顯著差異(表3)。

    將有效的土壤理化因子與尖孢鐮刀菌構建一般線性模型,用于評估土壤理化因子抑制病原菌增殖的可能性大小及預測其相對重要性。結果表明,5個土壤理化因子解釋了總模型變量的87.33%,其中,pH是影響病原菌濃度的一個顯著且關鍵的微生物生態(tài)指標(Plt;0.05),它對模型總變量的解釋度達到26.72%(表4)。

    將對病原菌具有顯著負相關關系的微生物指"

    示因子與病原菌再次構建線性模型,用于預測微生物因子抑制病原菌濃度的相對重要性。結果表明,細菌群落、真菌群落及子囊菌門(Ascomycota)類群相對豐度是抑制病原菌的顯著指示因子(Plt;0.05),分別解釋了17.01%、8.05%11.22%的模型總變量(表4)。

    3""討論

    連作障礙是指在同一片土地上連續(xù)種植相同作物,導致土壤養(yǎng)分失衡、病原菌累積和微生物群落失衡,進而限制作物生長,降低產量和品質的現象。香草蘭是一種熱帶多年生藤本植物,尤其容易受到連作障礙的影響。本團隊前期研究表明,通過輪作種植多年生作物胡椒可以有效改善土壤真菌群落結構及其組成,增強香草蘭的抗病性[11]。進一步篩選熱帶特色、種植周期短且經濟效益高的香料作物斑蘭葉和糯米香茶進行輪作盆栽試驗,結果顯示,輪作斑蘭葉能顯著降低土壤中鐮刀菌屬的相對豐度,特別是減少了尖孢鐮刀菌的數量。這一發(fā)現與先前的輪作研究相符,證實了輪作可以有效控制病原菌的增長[11]。此外,類似的通過輪作減少作物土傳病害的研究亦廣泛開展。如辣椒-香蕉輪作[12]、茄子-香蕉輪作[28]和菠蘿-香蕉輪作[29]模式均能有效降低根際土壤中尖孢鐮刀菌的數量,增強香蕉的抗病性。輪作改變了植物宿主,打斷了病原菌的營養(yǎng)循環(huán)[3],這可能是提升香草蘭抗病能力的關鍵因素。

    土壤微生物的多樣性和豐富度在植物應對環(huán)境壓力時發(fā)揮著關鍵作用,它們協(xié)助植物抵御外來侵害并促進植物恢復[3]。本研究發(fā)現,輪作周期短且經濟價值高的斑蘭葉和糯米香茶顯著提高了土壤細菌群落的豐富度和多樣性,而降低了土壤真菌群落的豐富度和多樣性。前期研究也發(fā)現,輪作多年生胡椒和咖啡能提高土壤真菌群落的豐富度和多樣性[11]。表明不同作物的輪作模式對細菌和真菌群落的豐富度和多樣性的影響有差異。特定的輪作作物通過激發(fā)特定的有益微生物類群的繁殖和生長,進而影響土壤微生物群落結構。而這些微生物與病原菌爭奪植物根部的生態(tài)位資源,限制了病原菌的營養(yǎng)來源,從而能有效控制其增長[30]。

    已有研究表明,長期作物輪作后微生物群落結構和組成與連作栽培模式間有顯著差異[10]。本研究中,無論細菌還是真菌群落結構,撂荒、連作和輪作均能顯著區(qū)分,形成3種完全不同的群落結構,但輪作胡椒/斑蘭葉/糯米香茶之間的群落結構很相似,無法區(qū)分。這一發(fā)現與前期的輪作研究結果[11]一致。推測可能是輪作的時長不夠,還不足以富集到更多特異的微生物。由于輪作效應,微生物群落的變化將影響特定微生物的響應。進一步研究發(fā)現,輪作能顯著富集細菌門水平上的酸桿菌門(Acidobacteria)和綠彎菌門Chloroflexi)微生物,屬水平則顯著富集類諾卡氏屬(Nocar dioides)微生物,且與病原菌呈顯著負相關;而真菌則顯著富集門水平上的子囊菌門(Ascomy cota)和屬水平上的Ascobolus屬微生物。雖然目前尚無足夠的證據表明這些微生物具有直接的抗病功能。但值得注意的是,這些微生物都是土壤中的優(yōu)勢菌門,可能是通過間接影響其他微生物的定殖而發(fā)揮抗病功能[31]。

    為鑒定識別出各輪作作物特異富集且對病原菌有潛在拮抗功能的核心微生物,采用專業(yè)尋找生物標志物的隨機森林分析,研究發(fā)現,輪作斑蘭葉后能特異富集Terriglobus屬、克雷伯桿菌屬(Kribbella)和粘球菌屬(Myxococcus)核心細菌類群,以及枝頂孢霉屬(Acremonium)和帚枝霉屬(Sarocladium)真菌核心類群,且均與病原菌呈顯著負相關關系。值得注意的是,前期的胡椒輪作能顯著富集有益的拮抗真菌類群木霉屬(Trichoderma)和青霉屬(Penicillium),表明不同的輪作作物根系代謝物有差異,富集的特定有益微生物類群也有差異[32]。已有研究表明,克雷伯桿菌屬和粘球菌屬的微生物在土壤中廣泛分布,而枝頂孢霉屬及帚枝霉屬均屬于叢梗孢科Hypomycetes)真菌微生物,均能產生多樣的次級代謝產物,而這些產物具有抗菌、抗生素和其他活性物質,是生物防治方面重要的微生物資源[33-35]。因此推測:這些具備生物防治特性的微生物種群豐度的變化可能影響了輪作過程中的病原菌數量,進而增強了香草蘭對病害的抵抗力。

    基于數學算法構建的微生物共發(fā)生關系網絡,可用于研究微生物群落的多樣性、復雜性及穩(wěn)定性,近年來已被廣泛應用[36-38]?;诖耍瑢α袒?、連作、輪作胡椒、輪作斑蘭葉及輪作糯米香茶5個處理的根際土壤細菌和真菌分別構建了微生物網絡,通過綜合考量網絡特征系數,如正相關連接邊數、平均度、平均聚類系數等,發(fā)現輪作斑蘭葉后的細菌和真菌物種間互惠共生的正相關作用更強,可能通過ASV物種間更緊密的互作關系提高了根際微生物多樣性[39],從而提升了菌群協(xié)同抵抗病原菌入侵根際能力。

    前期研究得出,在高發(fā)枯萎病的香蕉園輪作辣椒后,通過提高土壤pH可重塑微生物組成,富集有益的假單胞菌屬(Pseudomonas)微生物,增強了香蕉的抗病性[14]。本研究也發(fā)現,輪作斑蘭葉和糯米香茶能顯著提升根際土壤pH,而土壤pH是限制病原菌的最顯著因子,可解釋模型變量的26.72%。進一步利用一般線性模型分析揭示了細菌群落、真菌群落以及子囊菌門的相對豐度在抑制病原菌方面起著顯著的作用。這與TAO[40]的研究結果相符合,即根際細菌群落的變化對于抑制香蕉土傳枯萎病至關重要。

    4""結論

    在土傳病害頻發(fā)的香草蘭園輪作斑蘭葉可顯著降低鐮刀菌屬的相對豐度和尖孢鐮刀菌的拷貝數,增加細菌群落的豐富度和多樣性。通過特異富集細菌Terriglobus屬、克雷伯桿菌屬(Krib bella)和粘球菌屬(Myxococcus)核心類群,以及真菌枝頂孢霉屬(Acremonium)及帚枝霉屬(Sarocladium)核心類群可重塑微生物群落組成。輪作斑蘭葉作物后可顯著提高土壤pH,土壤pH提升的同時影響了細菌和真菌的群落結構,從而抑制了病原菌的定殖,推測這些變化是輪作提高香草蘭抗病能力的重要原因。

    參考文獻

    [1]"MINOO"D,"JAYAKUMAR"V,"VEENA"S,"VIMALA"J,"BASHA"A,"SAJI"K,"BABU"K"N,"PETER"K."Genetic"variations"and"interrelationships"in"Vanilla"planifolia"and"few"related"species"as"expressed"by"RAPD"polymorphism[J]."Genetic"Resources"and"Crop"Evolution,"2008,"55(3):"459-470.

    [2]"XIONG"W,"ZHAO"Q"Y,"ZHAO"J,"XUN"W,"LI"R,"ZHANG"R,"WU"H,"SHEN"Q."Different"continuous"cropping"spans"significantly"affect"microbial"community"membership"and"structure"in"a"vanilla-grown"soil"as"revealed"by"deep"pyrosequencing[J]."Microbial"Ecology,"2015,"70(1):"209-218.

    [3]"TRIVEDI"P,"LEACH"J"E,"TRINGE"S"G,"SA"T"M,"SINGH"B"K."Plant-microbiome"interactions:"from"community"assembly"to"plant"health[J]."Nature"Reviews"Microbiology,"2020,"18(11):"607-621.

    [4]"PEDRINHO"A,"MENDES"L,"DE"ARAUJO"PEREIRA"A,"VAISHNV"A,"KARPOUZAS,"SINGH"B"K."Soil"microbial"diversity"plays"an"important"role"in"resisting"and"restoring"degraded"ecosystems[J]."Plant"and"Soil,"2024,"500(1):"325-"349.

    [5]"ROMERO"F,"HILFIKER"S,"EDLINGER"A,"HELD"A,"HARTMAN"K,"LABOUYRIE"M,"VAN"DER"M."Soil"microbial"biodiversity"promotes"crop"productivity"and"agro-"ecosystem"functioning"in"experimental"microcosms[J]."Science"of"The"Total"Environment,"2023,"885:"163683.

    [6]"YARZáBAL"R,"áLVAREZ"P,"GUNDE-CIMERMAN"N,"CIANCAS"J,"GUTIéRREZ-CEPEDA"A,"OCA?A"A,"BATISTA-GARCíA"R."Exploring"extremophilic"fungi"in"soil"mycobiome"for"sustainable"agriculture"amid"global"change[J]."Nature"Communications,"2024,"15(1):"6951.

    [7]"HARTMANN"M,"SIX"J."Soil"structure"and"microbiome"functions"in"agroecosystems[J]."Nature"Reviews"Earth"amp;"Environment,"2023,"4(1):"4-18.

    [8]"LONGEPIERRE"M,"WIDMER"F,"KELLER"T,"WEISSKOPF"P,"COLOMBI"T,"SIX"J,"HARTMANN"M."Limited"resilience"of"the"soil"microbiome"to"mechanical"compaction"within"four"growing"seasons"of"agricultural"management[J]."ISME"Communications,"2021,"1(1):"44.

    [9]"YANG"T"K,"SIDDIQUE"H"M,"LIU"K."Cropping"systems"in"agriculture"and"their"impact"on"soil"health:"a"review[J]."Global"Ecology"and"Conservation,"2020,"23:"e01118.

    [10]"ZHOU"Y,"YANG"Z,"LIU"J,"LI"X,"WANG"X,"DAI"C,"ZHANG"T,"CARRION"V"J,"WEI"Z,"CAO"F."Crop"rotation"and"native"microbiome"inoculation"restore"soil"capacity"to"suppress"a"root"disease[J]."Nature"Communications,"2023,"14(1):"8126.

    [11]"XIONG"W,"ZHAO"Q"Y,"XUE"C,"XUN"W"B,"ZHAO"J,"WU"H"S,"LI"R,"SHEN"Q"R."Comparison"of"fungal"community"in"black"pepper-vanilla"and"vanilla"monoculture"systems"associated"with"vanilla"Fusarium"wilt"disease"[J]."Frontiers"in"Microbiology,"2016,"7:"117.

    [12]"HONG"S,"YUAN"X"F,"YANG"J"M,"YANG"Y,"JV"H"L,"LI"R,"JIA"Z"J,"RUAN"Y"Z."Selection"of"rhizosphere"communities"of"diverse"rotation"crops"reveals"unique"core"microbiome"associated"with"reduced"banana"Fusarium"wilt"disease[J]."New"Phytologist,"2023,"238(5):"2194-2209.

    [13]"XIONG"W,"LI"R,"REN"Y,"LIU"C,"ZHAO"Q"Y,"WU"H"S,"JOUSSET"A,"SHEN"Q"R."Distinct"roles"for"soil"fungal"and"bacterial"communities"associated"with"the"suppression"of"vanilla"Fusarium"wilt"disease[J]."Soil"Biology"and"Biochemistry,"2017,"107:"198-207.

    [14]"HONG"S,"JV"H"L,"LU"M,"WANG"B"B,"ZHAO"Y,"RUAN"Y"Z."Significant"decline"in"banana"Fusarium"wilt"disease"is"associated"with"soil"microbiome"reconstruction"under"chilli"pepper-banana"rotation[J]."European"Journal"of"Soil"Biology,"2020,"97:"103154.

    [15]"CHEN"Y,"BONKOWSKI"M,"SHEN"Y,"GRIFFITHS"B"S,"JIANG"Y,"WANG"X,"SUN"B."Root"ethylene"mediates"rhizosphere"microbial"community"reconstruction"when"chemi cally"detecting"cyanide"produced"by"neighbouring"plants[J]."Microbiome,"2020,"8:"4.

    [16]"鮑士旦."土壤農化分析[M]."3"."北京:"中國農業(yè)出版社,"2000.BAO"S"D."Soil"analysis"in"agricultural"chemistry[M]."3rd"ed."Beijing:"China"Agricultural"Press,"2000."(in"Chinese)

    [17]"CAPORASO"J"G,"LAUBER"C"L,"WALTERS"W"A,"BERG"L"D,"LOZUPONE"C"A,"TURNBAUGH"P"J,"FIERER"N,"KNIGHT"R."Global"patterns"of"16S"rRNA"diversity"at"a"depth"of"millions"of"sequences"per"sample[J]."Proceedings"of"the"National"Academy"of"Sciences"of"the"United"States"of"America,"2011,"108:"4516-4522.

    [18]"SCHOCH"C"L,"SEIFERT"K"A,"HUHNDORF"S,"ROBERT"V,"SPOUGE"J"L,"LEVESQUE"C"A,"CHEN"W."Nuclear"ribosomal"internal"transcribed"spacer"(ITS)"region"as"a"universal"DNA"barcode"marker"for"fungi[J]."Proceedings"of"the"National"Academy"of"Sciences"of"the"United"States"of"America,"2012,"109(16):"6241-6246.

    [19]"BOLGER"A"M,"LOHSE"M,"USADEL"B."Trimmomatic:"a"flexible"trimmer"for"Illumina"sequence"data[J]."Bioinformatics,"2014,"30(15):"2114-2120.

    [20]"MARTIN"M."Cutadapt"removes"adapter"sequences"from"high-throughput"sequencing"reads[J]."EMBnet"Journal,"2011,"17(1):"10-12.

    [21]"EDGAR"R"C."UPARSE:"highly"accurate"OTU"sequences"from"microbial"amplicon"reads[J]."Nature"Methods,"2013,"10(10):"996-999

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

    [23]"BOLYEN"E,"RIDEOUT"J"R,"DILLON"M"R,"BOKULICH"N"A,"ABNET"C"C,"GHALITH"G"A,"ALEXANDER"H,"ALM"E"J,"ARUMUGAM"M,"ASNICAR"F."Reproducible,"interactive,"scalable"and"extensible"microbiome"data"science"using"QIIME"2[J]."Nature"Biotechnology,"2019,"37(8):"852-857.

    [24]"CALLAHAN"B"J,"MCMURDIE"P"J,"ROSEN"M"J,"HAN"A"W,"JOHNSON"A"J"A,"HOLMES"S"P."DADA2:"high-resolu tion"sample"inference"from"Illumina"amplicon"data[J]."Nature"Methods,"2016,"13(7):"581-583.

    [25]"WANG"Q,"GARRITY"G"M,"TIEDJE"J"M,"COLE"J"R."Naive"bayesian"classifier"for"rapid"assignment"of"rRNA"sequences"into"the"new"bacterial"taxonomy[J]."Applied"and"Environmental"Microbiology,"2007,"73(16):"5261-5267.

    [26]"K?LJALG"U,"NILSSON"R"H,"ABARENKOV"K,"TEDERSOO"L,"TAYLOR"A"F"S,"BAHRAM"M,"BATES"S"T,"BRUNS"T"D,"BENGTSSON-PALME"J,"CALLAGHAN"T"M."Towards"a"unified"paradigm"for"sequence-based"identification"of"fungi[J]."Molecular"Ecology,"2013,"22(21):"5271-"5277.

    [27]"趙娜,"李榮,"辛侃,"趙艷,"阮云澤,"符常明."茄科蔬菜輪作對高發(fā)枯萎病蕉園土壤可培養(yǎng)微生物的影響[J]."熱帶作物學報,"2014,"35(8):"1469-1474.ZHAO"N,"LI"R,"XIN"K,"ZHAO"Y,"RUAN"Y"Z,"FU"C"M."Effects"of"different"Solanaceae"crop"rotations"on"the"soil"culturable"microbes"in"an"orchard"with"serious"Fusarium"wilt"disease[J]."Chinese"Journal"of"Tropical"Crops,"2014,"35(8):"1469-1474."(in"Chinese)"

    [28]"YUAN"X"F,"WANG"B"B,"HONG"S,"XIONG"W,"SHEN"Z"Z,"RUAN"Y"Z,"LI"R,"SHEN"Q"R,"DINI-ANDREOTE"F."Promoting"soil"microbial-mediated"suppressiveness"against"Fusarium"wilt"disease"by"the"enrichment"of"specific"fungal"taxa"via"crop"rotation[J]."Biology"and"Fertility"of"Soils,"2021,"57(8):"1137-1153.

    [29]"CHA"J"Y,"HAN"S,"HONG"H"J,"CHO"H,"KIM"D,"KWON"Y,"KWON"S"K,"CRüSEMANN"M,"LEE"Y"B,"KIM"J"F."Microbial"and"biochemical"basis"of"a"Fusarium"wilt-suppressive"soil[J]."The"ISME"Journal,"2016,"10(1):"119-129.

    [30]"WEN"T,"DING"Z,"THOMASHOW"L"S,"HALE"L,"YANG"S,"XIE"P,"LIU"X,"WANG"H,"SHEN"Q,"YUAN"J."Deciphering"the"mechanism"of"fungal"pathogen-induced"disease-suppres sive"soil[J]."New"Phytologist,"2023,"238(6):"2634-2650.

    [31]"HUBER"K"J,"PESTER"M,"EICHORST"S"A,"NAVARRETE"A"A,"FOESEL"B"U."Acidobacteria:"towards"unraveling"the"secrets"of"a"widespread,"though"enigmatic,"phylum[J]."Frontiers"in"Microbiology,"2022,"13:"960602.

    [32]"WEN"T,"XIE"P"H,"LIU"H"W,"LIU"T,"ZHAO"M"L,"YANG"S"D,"NIU"G"Q,"HALE"L,"SINGH"B"K,"KOWALCHUK"G"A,"SHEN"Q"R,"YUAN"J."Tapping"the"rhizosphere"metabolites"for"the"prebiotic"control"of"soil-borne"bacterial"wilt"disease[J]."Nature"Communications,"2023,"14:"4497.

    [33]"VIRUéS-SEGOVIA"J"R,"REYES"F,"RUíZ"S,"MARTíN"J,"FERNáNDEZ-PASTOR"I,"JUSTICIA"C,"DE"M,"DíAZ"C,"MACKENZIE"T"A,"GENILLOUD"O."Kribbellichelins"A"and"B,"two"new"antibiotics"from"Kribbella"sp."CA-293567"with"activity"against"several"human"pathogens[J]."Molecules"(Basel,"Switzerland),"2022,"27(19):"6355.

    [34]"SAGGU"S"K,"NATH"A,"KUMAR"S."Myxobacteria:"biology"and"bioactive"secondary"metabolites[J]."Research"in"Microbiology,"2023,"174(7):"104079.

    [35]"PERDOMO"H,"SUTTON"D"A,"GARCíA"D,"FOTHERGILL"A"W,"CANO"J,"GENé"J,"SUMMERBELL"R"C,"RINALDI"M"G,"GUARRO"J."Spectrum"of"clinically"relevant"Acremonium"species"in"the"United"States[J]."Journal"of"Clinical"Microbiology,"2011,"49(1):"243-256.

    [36]"XIONG"C,"ZHU"Y"G,"WANG"J"T,"SINGH"B,"HAN"L"L,"SHEN"J"P,"LI"P"P,"WANG"G"B,"WU"C"F,"GE"A"H."Host"selection"shapes"crop"microbiome"assembly"and"network"complexity[J]."New"Phytologist,"2021,"229(2):"1091-1104.

    [37]"DE"VRIES"F"T,"GRIFFITHS"R"I,"BAILEY"M,"CRAIG"H,"GIRLANDA"M,"GWEON"H"S,"HALLIN"S,"KAISERMANN"A,"KEITH"A"M,"KRETZSCHMAR"M."Soil"bacterial"networks"are"less"stable"under"drought"than"fungal"networks[J]."Nature"Communications,"2018,"9:"3033.

    [38]"GUSEVA"K,"DARCY"S,"SIMON"E,"ALTEIO"L"V,"MONTESINOS-NAVARRO"A,"KAISER"C."From"diversity"to"complexity:"microbial"networks"in"soils[J]."Soil"Biology"and"Biochemistry,"2022,"169:"108604.

    [39]"WU"H,"GAO"T"H,"HU"A,"WANG"J"J."Network"complexity"and"stability"of"microbes"enhanced"by"microplastic"diversity[J]."Environmental"Science"amp;"Technology,"2024,"58:"4334-4345.

    [40]"TAO"C"Y,"LI"R,"XIONG"W,"SHEN"Z"Z,"LIU"S"S,"WANG"B"B,"RUAN"Y"Z,"GEISEN"S"F,"SHEN"Q"R,"KOWALCHUK"G"A."Bio-organic"fertilizers"stimulate"indigenous"soil"Pseudomonas"populations"to"enhance"plant"disease"suppression[J]."Microbiome,"2020,"8(1):"137.

    "

    猜你喜歡
    宋體方正香草
    不同濃度微量元素處理對斑蘭葉組培過程中揮發(fā)性物質及角鯊烯相關酶活性的影響
    乙烯利處理菠蘿不同成熟期果實中AcMYB103基因表達規(guī)律研究
    淺析射陽大米的由來
    蘑菇頭和小香草
    蘑菇頭和小香草
    蘑菇頭和小香草
    方正
    香草
    天池小小說(2009年9期)2009-09-22 08:04:06
    福布斯中國大陸50富豪排行榜
    全國新書目(2001年5期)2001-06-20 18:27:40
    武漢迅通集郵收藏品有限公司——售品目錄
    收藏(2001年3期)2001-06-14 05:12:36
    亚洲精品国产av蜜桃| 亚洲欧美一区二区三区黑人 | 中国国产av一级| 欧美人与性动交α欧美软件| 999久久久国产精品视频| 欧美+日韩+精品| 成人毛片a级毛片在线播放| 色婷婷av一区二区三区视频| 99九九在线精品视频| 婷婷色av中文字幕| 免费大片黄手机在线观看| 国产精品成人在线| 黄色怎么调成土黄色| 亚洲精品成人av观看孕妇| 亚洲国产看品久久| 母亲3免费完整高清在线观看 | 天天躁夜夜躁狠狠躁躁| 亚洲av电影在线进入| 69精品国产乱码久久久| 中文字幕av电影在线播放| 国产成人aa在线观看| 欧美成人午夜精品| 韩国高清视频一区二区三区| 亚洲国产精品999| 国产熟女欧美一区二区| 午夜av观看不卡| av免费在线看不卡| 看免费av毛片| 亚洲国产精品成人久久小说| 韩国av在线不卡| 久久精品国产a三级三级三级| 新久久久久国产一级毛片| 欧美+日韩+精品| 18在线观看网站| 大码成人一级视频| 国产av国产精品国产| 一级a爱视频在线免费观看| 久久久久久人妻| 亚洲人成网站在线观看播放| 性高湖久久久久久久久免费观看| freevideosex欧美| 爱豆传媒免费全集在线观看| 综合色丁香网| 欧美精品亚洲一区二区| 乱人伦中国视频| 日本欧美视频一区| 黄色配什么色好看| 日产精品乱码卡一卡2卡三| 国产一区亚洲一区在线观看| 色视频在线一区二区三区| 亚洲欧美清纯卡通| 久久久精品区二区三区| 国产精品不卡视频一区二区| 欧美精品高潮呻吟av久久| h视频一区二区三区| 久久毛片免费看一区二区三区| 99国产精品免费福利视频| 国产不卡av网站在线观看| 欧美黄色片欧美黄色片| 国产精品二区激情视频| 99国产综合亚洲精品| 国产福利在线免费观看视频| 七月丁香在线播放| 女性被躁到高潮视频| 亚洲三级黄色毛片| 亚洲av国产av综合av卡| 久久久久人妻精品一区果冻| 丝袜在线中文字幕| 亚洲成人av在线免费| 亚洲av成人精品一二三区| 国产免费视频播放在线视频| 波多野结衣av一区二区av| 精品国产乱码久久久久久男人| 最近最新中文字幕免费大全7| 一级毛片我不卡| 新久久久久国产一级毛片| 精品国产乱码久久久久久男人| 日韩三级伦理在线观看| 一级黄片播放器| 宅男免费午夜| 国产精品偷伦视频观看了| 中国国产av一级| 免费看av在线观看网站| 免费高清在线观看日韩| 久久久国产一区二区| 亚洲国产色片| 热99久久久久精品小说推荐| 最新中文字幕久久久久| 一区二区三区乱码不卡18| 国产亚洲av片在线观看秒播厂| 精品国产乱码久久久久久男人| 国产精品人妻久久久影院| 欧美国产精品va在线观看不卡| 男人添女人高潮全过程视频| 只有这里有精品99| 亚洲 欧美一区二区三区| 三上悠亚av全集在线观看| 国产免费现黄频在线看| 亚洲综合色网址| 日韩 亚洲 欧美在线| 国产精品国产av在线观看| 国产极品天堂在线| 国产成人av激情在线播放| 国产精品蜜桃在线观看| 成年美女黄网站色视频大全免费| 久久ye,这里只有精品| 麻豆精品久久久久久蜜桃| 在线精品无人区一区二区三| 菩萨蛮人人尽说江南好唐韦庄| 最近手机中文字幕大全| 国产又色又爽无遮挡免| 搡女人真爽免费视频火全软件| 精品国产国语对白av| 免费黄频网站在线观看国产| 日韩av在线免费看完整版不卡| 激情五月婷婷亚洲| 国产成人欧美| 在线看a的网站| 男女国产视频网站| 久久久久久人人人人人| 你懂的网址亚洲精品在线观看| 一级片'在线观看视频| 97人妻天天添夜夜摸| 国产精品国产av在线观看| 波多野结衣一区麻豆| 男男h啪啪无遮挡| 韩国精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 国产片内射在线| 丝袜在线中文字幕| 国产精品av久久久久免费| 亚洲综合色网址| 久久人人爽人人片av| 亚洲,欧美,日韩| 欧美日韩一级在线毛片| 在线 av 中文字幕| 国产一区二区激情短视频 | 欧美成人精品欧美一级黄| 亚洲经典国产精华液单| 国产精品亚洲av一区麻豆 | 久久婷婷青草| 日本色播在线视频| 女性被躁到高潮视频| 熟女av电影| 日韩精品有码人妻一区| 下体分泌物呈黄色| 国产一区二区激情短视频 | 欧美激情 高清一区二区三区| 日本色播在线视频| 男人爽女人下面视频在线观看| 宅男免费午夜| 极品人妻少妇av视频| www日本在线高清视频| 国产不卡av网站在线观看| 99香蕉大伊视频| 亚洲国产最新在线播放| 久久亚洲国产成人精品v| 久久精品国产自在天天线| 一级毛片我不卡| 午夜福利一区二区在线看| 欧美日韩视频精品一区| 夜夜骑夜夜射夜夜干| 久热这里只有精品99| 两性夫妻黄色片| 王馨瑶露胸无遮挡在线观看| 熟女av电影| 少妇人妻久久综合中文| 欧美日韩av久久| 日本黄色日本黄色录像| 一边亲一边摸免费视频| 自线自在国产av| 满18在线观看网站| 国产成人免费无遮挡视频| 国产xxxxx性猛交| 欧美日韩成人在线一区二区| 中文字幕精品免费在线观看视频| 91精品三级在线观看| 色婷婷av一区二区三区视频| 在线观看www视频免费| 国产成人精品在线电影| 国产男女内射视频| 纯流量卡能插随身wifi吗| 肉色欧美久久久久久久蜜桃| av国产精品久久久久影院| 一个人免费看片子| 18在线观看网站| 精品少妇内射三级| 一区二区日韩欧美中文字幕| 99热国产这里只有精品6| 国产精品熟女久久久久浪| 1024香蕉在线观看| 国产麻豆69| 美女国产高潮福利片在线看| 久久国产精品大桥未久av| av一本久久久久| 丝瓜视频免费看黄片| 色网站视频免费| 一级毛片电影观看| 免费黄网站久久成人精品| 在线看a的网站| 日韩成人av中文字幕在线观看| 精品亚洲成a人片在线观看| 一级,二级,三级黄色视频| 91精品伊人久久大香线蕉| 日韩视频在线欧美| 久久久精品国产亚洲av高清涩受| 国产亚洲一区二区精品| 国产女主播在线喷水免费视频网站| www.熟女人妻精品国产| 国产爽快片一区二区三区| 精品福利永久在线观看| 咕卡用的链子| 国产在线视频一区二区| 男女高潮啪啪啪动态图| 波多野结衣av一区二区av| 青草久久国产| 欧美日韩精品网址| 肉色欧美久久久久久久蜜桃| 男男h啪啪无遮挡| 午夜福利一区二区在线看| 久久久精品国产亚洲av高清涩受| 国产精品久久久久久精品古装| 免费观看av网站的网址| 久久 成人 亚洲| 男女无遮挡免费网站观看| 一区在线观看完整版| 国产老妇伦熟女老妇高清| 看免费成人av毛片| 男女高潮啪啪啪动态图| 欧美在线黄色| 精品一区二区免费观看| 欧美 日韩 精品 国产| 国产福利在线免费观看视频| 美国免费a级毛片| 亚洲国产精品999| av又黄又爽大尺度在线免费看| 亚洲成人一二三区av| 亚洲第一区二区三区不卡| 黄片播放在线免费| 国产成人午夜福利电影在线观看| 亚洲综合色网址| 国产爽快片一区二区三区| 麻豆av在线久日| 一级毛片我不卡| 91精品三级在线观看| 免费久久久久久久精品成人欧美视频| 欧美成人午夜免费资源| 婷婷色综合大香蕉| 99热全是精品| 交换朋友夫妻互换小说| 高清欧美精品videossex| 日本黄色日本黄色录像| 亚洲av电影在线观看一区二区三区| 91成人精品电影| 一区在线观看完整版| 九色亚洲精品在线播放| 亚洲av日韩在线播放| 久久 成人 亚洲| 国产不卡av网站在线观看| 中国三级夫妇交换| 九九爱精品视频在线观看| 最新的欧美精品一区二区| 日韩免费高清中文字幕av| 亚洲久久久国产精品| 中文精品一卡2卡3卡4更新| 久久久久久久国产电影| av免费观看日本| 国产精品99久久99久久久不卡 | 国产成人免费观看mmmm| 晚上一个人看的免费电影| 久久久久久久久久久久大奶| 亚洲精品日本国产第一区| 久久女婷五月综合色啪小说| 久久青草综合色| 午夜福利乱码中文字幕| 精品午夜福利在线看| 精品99又大又爽又粗少妇毛片| 视频在线观看一区二区三区| 女性被躁到高潮视频| 老司机影院毛片| 久久综合国产亚洲精品| 宅男免费午夜| 中文字幕最新亚洲高清| 日韩欧美一区视频在线观看| 在线观看三级黄色| 免费黄网站久久成人精品| 欧美精品人与动牲交sv欧美| 国产男人的电影天堂91| 五月伊人婷婷丁香| 亚洲婷婷狠狠爱综合网| 精品亚洲成国产av| 国产高清不卡午夜福利| 色94色欧美一区二区| 成人国语在线视频| 人妻人人澡人人爽人人| 欧美在线黄色| 91午夜精品亚洲一区二区三区| 精品午夜福利在线看| 一级毛片我不卡| 天堂俺去俺来也www色官网| 男人添女人高潮全过程视频| 亚洲一码二码三码区别大吗| 最新中文字幕久久久久| 国产免费视频播放在线视频| 九九爱精品视频在线观看| 性高湖久久久久久久久免费观看| 久久精品国产自在天天线| 久久人人爽av亚洲精品天堂| 国产无遮挡羞羞视频在线观看| 免费大片黄手机在线观看| 大码成人一级视频| 在线观看一区二区三区激情| 久久久久视频综合| 亚洲国产毛片av蜜桃av| 一本色道久久久久久精品综合| 黑人巨大精品欧美一区二区蜜桃| 日本色播在线视频| 丝袜在线中文字幕| 欧美日韩精品成人综合77777| 一级黄片播放器| 美女大奶头黄色视频| 国产精品嫩草影院av在线观看| 亚洲人成77777在线视频| 久久 成人 亚洲| 精品人妻熟女毛片av久久网站| 欧美精品av麻豆av| 久久亚洲国产成人精品v| 一级毛片电影观看| 国产又色又爽无遮挡免| 中国国产av一级| 欧美变态另类bdsm刘玥| 国产精品 欧美亚洲| 国产成人精品一,二区| 日本vs欧美在线观看视频| 熟妇人妻不卡中文字幕| 亚洲美女视频黄频| 亚洲少妇的诱惑av| 亚洲av男天堂| 免费黄频网站在线观看国产| 亚洲综合色惰| av在线老鸭窝| 国产精品人妻久久久影院| 久久久国产欧美日韩av| 久久久久视频综合| 美女中出高潮动态图| 免费不卡的大黄色大毛片视频在线观看| 亚洲少妇的诱惑av| 一本—道久久a久久精品蜜桃钙片| 欧美成人午夜免费资源| 日本vs欧美在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 精品一区二区免费观看| 久久久a久久爽久久v久久| 亚洲伊人久久精品综合| 精品亚洲成国产av| 在线观看三级黄色| 精品一品国产午夜福利视频| 日韩三级伦理在线观看| 国产xxxxx性猛交| 国产精品人妻久久久影院| 国产一区有黄有色的免费视频| 日韩大片免费观看网站| 国产日韩欧美在线精品| 高清av免费在线| 久久狼人影院| 亚洲欧洲日产国产| kizo精华| 日本av手机在线免费观看| 亚洲男人天堂网一区| 久久久久久久久免费视频了| 捣出白浆h1v1| 久久 成人 亚洲| 女性生殖器流出的白浆| 嫩草影院入口| 国产av码专区亚洲av| 亚洲欧美清纯卡通| 纵有疾风起免费观看全集完整版| 最近的中文字幕免费完整| 久久久久久久久久久久大奶| 欧美成人午夜精品| 免费观看a级毛片全部| 欧美日韩亚洲高清精品| 日韩制服骚丝袜av| 在线观看人妻少妇| 亚洲精品乱久久久久久| h视频一区二区三区| 国产精品欧美亚洲77777| h视频一区二区三区| 欧美国产精品一级二级三级| 日韩在线高清观看一区二区三区| 国产视频首页在线观看| 久久精品国产亚洲av涩爱| 久久97久久精品| 女性生殖器流出的白浆| 亚洲精品国产av成人精品| 一级毛片黄色毛片免费观看视频| 国产老妇伦熟女老妇高清| 制服丝袜香蕉在线| 一二三四在线观看免费中文在| 亚洲一区二区三区欧美精品| 国产精品久久久久久精品古装| av免费观看日本| 天天躁夜夜躁狠狠躁躁| 精品国产一区二区久久| 国产乱人偷精品视频| 七月丁香在线播放| 一区福利在线观看| 久久午夜综合久久蜜桃| 九草在线视频观看| av电影中文网址| 亚洲av日韩在线播放| 国产精品99久久99久久久不卡 | 久久久久久伊人网av| 少妇精品久久久久久久| 色哟哟·www| 国产成人午夜福利电影在线观看| 精品少妇久久久久久888优播| 国产一区二区激情短视频 | 国产成人精品无人区| 水蜜桃什么品种好| 欧美精品一区二区免费开放| 纯流量卡能插随身wifi吗| 菩萨蛮人人尽说江南好唐韦庄| 久久久a久久爽久久v久久| 一级片'在线观看视频| 久久毛片免费看一区二区三区| 国产一区二区 视频在线| 亚洲国产日韩一区二区| 色播在线永久视频| 免费观看在线日韩| 国产精品99久久99久久久不卡 | 久久久久久久精品精品| 久久精品人人爽人人爽视色| 超碰成人久久| 亚洲av中文av极速乱| 一区二区三区激情视频| 看非洲黑人一级黄片| av在线app专区| 久久久久国产一级毛片高清牌| 精品视频人人做人人爽| 青青草视频在线视频观看| 伦理电影大哥的女人| 亚洲综合精品二区| 麻豆av在线久日| 亚洲精品中文字幕在线视频| 菩萨蛮人人尽说江南好唐韦庄| 巨乳人妻的诱惑在线观看| 有码 亚洲区| 亚洲第一青青草原| 精品卡一卡二卡四卡免费| 亚洲欧美一区二区三区黑人 | 一区二区av电影网| 国产免费又黄又爽又色| 美女xxoo啪啪120秒动态图| 九草在线视频观看| 亚洲av免费高清在线观看| 国产亚洲av片在线观看秒播厂| 亚洲精品乱久久久久久| 18+在线观看网站| 久久青草综合色| 人体艺术视频欧美日本| 欧美日韩av久久| 香蕉国产在线看| 精品亚洲乱码少妇综合久久| 有码 亚洲区| 欧美人与善性xxx| 国产免费一区二区三区四区乱码| 中文天堂在线官网| 国产一级毛片在线| 亚洲av.av天堂| 日韩一区二区视频免费看| 成人漫画全彩无遮挡| 捣出白浆h1v1| 99久久精品国产国产毛片| 精品人妻偷拍中文字幕| 久久久久久久大尺度免费视频| 搡女人真爽免费视频火全软件| 欧美在线黄色| 日韩成人av中文字幕在线观看| 老司机影院毛片| 人妻少妇偷人精品九色| 成人黄色视频免费在线看| 国产野战对白在线观看| av一本久久久久| www.av在线官网国产| 一本大道久久a久久精品| 美女大奶头黄色视频| 香蕉精品网在线| 女人精品久久久久毛片| 亚洲精品一二三| 国产1区2区3区精品| 国产免费福利视频在线观看| 建设人人有责人人尽责人人享有的| 人人妻人人澡人人看| 欧美97在线视频| 精品第一国产精品| 熟女电影av网| 黄网站色视频无遮挡免费观看| 精品一区二区三区四区五区乱码 | 在线免费观看不下载黄p国产| 欧美日韩精品成人综合77777| 亚洲欧美色中文字幕在线| 国产日韩一区二区三区精品不卡| 波野结衣二区三区在线| 久久久久国产精品人妻一区二区| 看非洲黑人一级黄片| 日本欧美视频一区| 精品少妇久久久久久888优播| 中文字幕另类日韩欧美亚洲嫩草| 建设人人有责人人尽责人人享有的| 亚洲国产av影院在线观看| 日韩精品免费视频一区二区三区| 国精品久久久久久国模美| 亚洲欧美精品自产自拍| 婷婷成人精品国产| 777米奇影视久久| 亚洲在久久综合| 我要看黄色一级片免费的| 最新的欧美精品一区二区| 乱人伦中国视频| 精品一区在线观看国产| 国产男女超爽视频在线观看| 国产成人aa在线观看| 视频区图区小说| 久久精品熟女亚洲av麻豆精品| 久久99热这里只频精品6学生| 午夜福利在线观看免费完整高清在| 麻豆精品久久久久久蜜桃| 国产成人91sexporn| 久久久欧美国产精品| 亚洲婷婷狠狠爱综合网| 18在线观看网站| av在线app专区| 成人漫画全彩无遮挡| 香蕉丝袜av| 精品国产超薄肉色丝袜足j| 国产成人免费观看mmmm| 欧美黄色片欧美黄色片| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产av影院在线观看| 欧美变态另类bdsm刘玥| 99九九在线精品视频| 欧美日韩成人在线一区二区| av在线老鸭窝| 日本爱情动作片www.在线观看| 亚洲av日韩在线播放| 久久久久久久大尺度免费视频| 熟女少妇亚洲综合色aaa.| 五月开心婷婷网| 一二三四在线观看免费中文在| 精品一区二区三区四区五区乱码 | 国产精品一二三区在线看| 国产97色在线日韩免费| 91精品三级在线观看| 在线 av 中文字幕| 久久久久精品久久久久真实原创| 国产av一区二区精品久久| 99久久综合免费| 免费高清在线观看视频在线观看| av在线老鸭窝| 啦啦啦中文免费视频观看日本| 永久免费av网站大全| 中文字幕亚洲精品专区| 久久鲁丝午夜福利片| 在线天堂中文资源库| 久久 成人 亚洲| 叶爱在线成人免费视频播放| 9热在线视频观看99| 日产精品乱码卡一卡2卡三| 亚洲av综合色区一区| 观看av在线不卡| 日日啪夜夜爽| 久久久久国产一级毛片高清牌| 日韩欧美一区视频在线观看| 久久人人爽av亚洲精品天堂| 精品酒店卫生间| 精品卡一卡二卡四卡免费| 亚洲欧美一区二区三区黑人 | 亚洲,欧美,日韩| 在线观看免费高清a一片| 黄色一级大片看看| 国产免费一区二区三区四区乱码| 最近中文字幕2019免费版| 五月天丁香电影| 老鸭窝网址在线观看| 国产日韩欧美视频二区| 黄片无遮挡物在线观看| 妹子高潮喷水视频| 观看av在线不卡| 国产在线一区二区三区精| 99久久中文字幕三级久久日本| 国产淫语在线视频| 美女xxoo啪啪120秒动态图| av有码第一页| 超碰成人久久| 2022亚洲国产成人精品| 晚上一个人看的免费电影| 国产97色在线日韩免费| 日韩欧美一区视频在线观看| 久久久久久久精品精品| 人体艺术视频欧美日本| 狂野欧美激情性bbbbbb| 欧美人与性动交α欧美精品济南到 | 亚洲av综合色区一区| 菩萨蛮人人尽说江南好唐韦庄| 青春草视频在线免费观看| 久久久久久久久免费视频了| 夫妻性生交免费视频一级片| 欧美国产精品va在线观看不卡| 一区二区三区精品91| 久久精品国产a三级三级三级| 国产亚洲欧美精品永久| 久久亚洲国产成人精品v| 精品少妇内射三级| 交换朋友夫妻互换小说|