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

    有機(jī)肥對香蕉枯萎病及土壤主要理化性質(zhì)和微生物群落的影響

    2021-11-12 10:23朱志炎梁雪雁林鳳玲馬嘉榮田志宏李建雄
    關(guān)鍵詞:高通量測序枯萎病有機(jī)肥

    朱志炎 梁雪雁 林鳳玲 馬嘉榮 田志宏 李建雄

    摘要:【目的】研究香蕉園施用有機(jī)肥防治香蕉枯萎病對土壤理化性質(zhì)及其根際土壤微生物群落的影響?!痉椒ā肯憬队酌缫圃灾敛^(qū)大田土壤中,處理組植株施用有機(jī)肥,并將未施用有機(jī)肥的植株設(shè)置為對照組。移栽后6個月統(tǒng)計(jì)處理組與對照組植株香蕉枯萎病發(fā)病率;采集土壤樣本,測定根際土壤的土壤肥力;提取根際土壤DNA,采用高通量測序技術(shù),結(jié)合生物信息學(xué)分析,解析施用有機(jī)肥后香蕉根際土壤細(xì)菌和真菌群落結(jié)構(gòu)組成及多樣性的變化?!窘Y(jié)果】施用有機(jī)肥提高了土壤pH值(14.85%)、全氮(25%)和全磷(19.04%)的含量,降低土壤全鐵含量(2.62%),香蕉枯萎病發(fā)病率下降了75%。和對照相比,子囊菌門(Ascomycota)與壺菌門(Chytridiomycota)的相對豐度分別提高了43.84%和90.64%,變形菌門(Proteobacteria)的相對豐度則降低了18.4g%。施用有機(jī)肥料提高了青霉菌屬(Penicillium)、Gibellulopsis和籃狀菌屬(Talaromyces)等的相對豐度,比例分別為對照組的2 93倍、2 12倍和11.93倍。施用有機(jī)肥料后,香蕉根際土壤真菌群落Chaol指數(shù)、ACE指數(shù)與香農(nóng)(Shannon)指數(shù)得到提升,分別提升了39.81%、38.43%和86.85%?!窘Y(jié)論】施用有機(jī)肥料改善了土壤理化性質(zhì),改變了根際土壤微生物群落結(jié)構(gòu)和多樣性,降低了香蕉枯萎病發(fā)病率。

    關(guān)鍵詞:枯萎病;香蕉;高通量測序;細(xì)菌群落;真菌群落;有機(jī)肥

    中圖分類號:S 432.1

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1008-0384(2021) 07-0806-11

    Effects of Bio-organic Fertilizer on Physicochemical Properties and Microflora of

    Banana Field Infected by Fusarium Wilt Disease

    ZHU Zhiyan l.2,LIAN Xueyan 2, LIN Fengling 2, MA Jiarong 2, TIAN Zhihong 1, LI Jianxiong 2*

    (1.School of Life Sciences, Yangtze University Jingzhou. Hubei 434025, China;2. South China Botanical Garden.

    Chinese Academy of Sciences. Guanzhou, Guandong 510650, China)

    Abstract: 【Objective】 Effects of a bio-organic fertilizer on physicochemical properties and microbial community inthizosphere soil of a banana field infected by Fusarium wilt disease were studied. 【Methods】 In a Fusarium wilt infectedbanana field. an random block design experiment on the application of a bio-organic fertilizer. Biofert. was conducted. Sixmonths after banana seedlings were transplanted to the field, thizosphere soil samples from lots with and without Biofertapplication were collected to determine the nutrient contents by chemical analysis and the microbial composition and diversityby high-throughput sequencing and bioinformatics analysis. 【Results】Compared to control, the application of Biofertincreased pH by 14.85%. the total N by 25%. and total P by 19.04%, but decreased the total Fe by 2.62% in the soil. whilelowered the incidence of Fusarium wilt on the plants by 75%. In the thizosphere soil. the Biofert-treated lots showed therelative abundance of. other than an 18.4g% decrease on Proteobacteria. increases by 43.84% on Ascomvcota. 90.64% onChvtridiomvcota. 293% on Penicillium. 212% on Gibellulopsis, and 1 193% on Talaromvces. 39.81% on fungal Chaol index,38.43% on ACE. and 86.85% on Shannon index. 【Conclusion】 Biofert application not only improved the soil quality butalso significanty altered the structure and diversity of the microbial commuunity in thizosphere soil and contributed to thereduced incidence of banana wilt disease

    Key words: Fussarium wilt; Banana; high-throughput sequencing; bacterial community; fungal community; organic fertilizerBiofert

    0 引言

    【研究意義】香蕉枯萎?。˙anana FusariumWilt)是一種危害香蕉生長的土傳病害,嚴(yán)重限制了香蕉產(chǎn)量的提升[1.2]。香蕉枯萎病致病菌尖孢鐮刀菌(Fusarium oxysporumf sp.Cubense)分為4個生理小種,其中以4號生理小種(Fusarium oxysporum f.sp.Cubense Tropical Race 4)危害范圍最廣,毒性最強(qiáng)。致病菌存在于土壤中,可以伴隨著香蕉幼苗的移栽而跨區(qū)域傳播[3]。通常致病菌通過侵染香蕉根系,在香蕉根部增殖,最終在香蕉原球莖聚集,導(dǎo)致原球莖中部變褐。原球莖的損傷限制了營養(yǎng)從根部流入葉片,造成香蕉葉片枯萎、黃化。病情嚴(yán)重的香蕉整株倒伏,直至絕產(chǎn)[4]。目前尚無控制香蕉枯萎病病害發(fā)生的有效措施,因此,探討該病的有效防控途徑迫在眉睫?!厩叭搜芯窟M(jìn)展】香蕉枯萎病的防治措施主要包括選育抗性品種、化學(xué)防治、輪作、生物防治等方法[5]。生物防治作為研究綜合防治香蕉枯萎病的措施之一,具有兼防兼治、無污染、利于環(huán)境保護(hù)和人畜安全的優(yōu)點(diǎn),且符合發(fā)展有機(jī)農(nóng)業(yè)的要求,備受人們重視[6]。Li[2]、ZHU等[7]研究發(fā)現(xiàn)鏈霉菌H4的代謝物可以抑制鐮刀菌孢子的萌發(fā),促進(jìn)細(xì)胞膜的降解,導(dǎo)致細(xì)胞內(nèi)物質(zhì)的外泄,在盆栽試驗(yàn)中,具有一定防治香蕉枯萎病的效果。Cheng等[8]利用植物內(nèi)生真菌印度梨形孢定殖于香蕉根部,提高香蕉體內(nèi)抗氧化酶活性,抵御香蕉枯萎病病原菌TR4的侵染。生物防治在植物病蟲害防治研究方面起到積極作用,含有生防菌的有機(jī)肥料更是引起越來越多的關(guān)注。施用含有生防菌的有機(jī)肥具有改善土壤結(jié)構(gòu)、促進(jìn)植物生長等優(yōu)點(diǎn)[9]。研究發(fā)現(xiàn)有機(jī)肥料還有修復(fù)土壤的作用,能夠使遭到破壞的細(xì)菌群落得到恢復(fù)[10]。含有解淀粉芽孢桿菌NJN-6的香蕉專用生物有機(jī)肥的施用能夠降低香蕉枯萎病發(fā)病率,促進(jìn)香蕉植株生長,增加香蕉產(chǎn)量,減少土壤尖孢鐮刀菌和真菌數(shù)量,增加細(xì)菌尤其是芽孢桿菌數(shù)量,并能提高土壤pH值、酶活力和硝態(tài)氮含量[11]。Tao等[12]研究表明添加了芽孢桿菌的有機(jī)肥料不僅可以改變土壤的理化性質(zhì),而且能夠促進(jìn)土壤中假單孢菌的功能增強(qiáng)。作物輪作中施用有機(jī)肥料也能在一定程度上防控枯萎病的發(fā)生,如在辣椒和香蕉輪作過程中,輪作辣椒配施生物有機(jī)肥降低了香蕉枯萎病的發(fā)病率[13]。隨著在香蕉生產(chǎn)上施用化肥造成的土壤問題日益突出,通過降低化肥的施用量,增施有機(jī)肥料使土壤肥力及微生物菌群得到改善,成為生產(chǎn)上病區(qū)大田防治香蕉枯萎病的一種值得研究的手段?!颈疚难芯壳腥朦c(diǎn)】有機(jī)肥在病區(qū)大田土壤中防治香蕉枯萎病研究報(bào)道不多,大多數(shù)研究主要集中在抑制致病菌生長以及溫室盆栽苗的研究上,直接應(yīng)用于大田香蕉枯萎病防治的探索報(bào)道較少。有機(jī)肥在大田中是否可以降低香蕉枯萎病發(fā)病情況及施用有機(jī)肥料后對病區(qū)土壤中微生物群落多樣性的影響有待進(jìn)一步研究?!緮M解決的關(guān)鍵問題】本研究對采白廣東省惠州市龍門縣病區(qū)大田未施用有機(jī)肥的植株與施用有機(jī)肥的植株的發(fā)病率進(jìn)行統(tǒng)計(jì),對土壤理化性質(zhì)進(jìn)行分析,并基于高通量測序技術(shù)對根際土壤中真菌與細(xì)菌群落結(jié)構(gòu)進(jìn)行分析,以期揭示施用有機(jī)肥后根際土壤中真菌與細(xì)菌的變化情況,探尋香蕉枯萎病發(fā)生與根際土壤微生物群落變化的微生態(tài)機(jī)制,從而為指導(dǎo)香蕉枯萎病的防控提供基礎(chǔ)信息支持。

    1材料與方法

    1.1試驗(yàn)材料

    巴西蕉(Musa acuminateL AAA group)[14.15]由中國科學(xué)院華南植物園蘭花組培研究組提供,有機(jī)肥(黃腐酸鉀、氨基酸、腐殖酸、煙粉、發(fā)酵玉米芯、EDTA-B、EDTA-Zn、EDTA-Ca、神農(nóng)共生菌與EM菌)由湖北龍旗生物科技有限公司提供。

    1.2主要試劑和儀器

    土壤基因組提取試劑盒(E.Z.N.A.?Soil DNAKit),MP Bio medicals公司;測序建庫試劑盒(NEBNext?UltraTM DNA Library Prep Kit for Illumina@),Illumina公司。DNA定量儀和MiSeq測序儀,Illumina公司;NanoDrop 2000紫外可見分光光度計(jì)、高速冷凍離心機(jī)、小型臺式離心機(jī)與PCR儀,賽默飛世爾科技公司。

    1.3試驗(yàn)設(shè)計(jì)與樣品采集

    巴西蕉幼苗(感病品種)[16]生長至4葉幼苗,將長勢一致的幼苗移栽至病區(qū)大田,60株幼苗隨機(jī)分成處理組和對照組。處理組( BIO)植株用500 g有機(jī)肥在根附近采用穴施;對照組( CK)不做任何處理。巴西蕉生長6個月后,進(jìn)行樣品采集,隨機(jī)挑取6株,鏟除表層土壤后,將整株連根拔起,去除附著于根部的大塊土壤,抖落并收集須根2 mm范圍的土壤,進(jìn)行混合。處理組3次重復(fù),樣品分別為BIOI、BI02與BI03;對照組3次重復(fù)樣品分別為CK1、CK2與CK3。樣品采集后放入低溫保藏箱,迅速保存80℃冰箱備用。

    1.4土壤部分理化性質(zhì)分析及病區(qū)土壤病原菌的分離

    土壤全氮含量測定采用凱氏定氮法,有所修改;土壤中全磷含量測定采用酸溶法;土壤pH采用電位法;全鐵采用氫氟酸一高氯酸一硝酸消煮法;土壤含水量參考烘干法[17]。香蕉枯萎病菌(Fusariumoxysporum f.sp.Cubense Tropical Race 4)TR4的分離采用梯度稀釋法,分離培養(yǎng)基參考Diane等[18]。香蕉枯萎病病原菌TR4初步鑒定采用的引物序列為(TR4特異性擴(kuò)增片段大小為463 bp):TR4-F(5'_CACGTTTAAGGTGCCATGAGAG-3 '). TR4-R(5'-GCCAGGACTGCCTCGTGA-3,)[19]。

    1.5 ITS與16S文庫構(gòu)建及高通量測序

    提取樣品總DNA后,根據(jù)保守區(qū)設(shè)計(jì)得到引物,在引物末端加上測序接頭,進(jìn)行PCR擴(kuò)增并對其產(chǎn)物進(jìn)行純化、定量和均一化形成測序文庫,建好的文庫先進(jìn)行文庫質(zhì)檢,質(zhì)檢合格的文庫用IlluminaHiSeq 2500進(jìn)行測序。

    (1) PE reads拼接:使用FLASH vl.2.7軟件[20],通過overlap對每個樣品的reads進(jìn)行拼接,得到的拼接序列即原始Tags數(shù)據(jù)(RawTags);

    ( 2) Tags過濾:使用Trimmomatic v0.33軟件[21],對拼接得到的Raw Tags進(jìn)行過濾,得到高質(zhì)量的Tags數(shù)據(jù)(Clean Tags);

    (3)去除嵌合體:使用UCHIME v4.2軟件,鑒定并去除嵌合體序列,得到最終有效數(shù)據(jù)( EffectiveTags)。

    1.6生物信息學(xué)分析

    對97%相似度序列使用UPARSE軟件[22]進(jìn)行操作分類單元( operational taxonomic units,OTU)聚類。對每條序列進(jìn)行物種分類注釋。聚類后的OTU用于稀釋曲線的構(gòu)建、Alpha多樣性分析(ACE指數(shù)、Chaol指數(shù)、Shannon指數(shù)、Simpson指數(shù)以及Coverage指數(shù))、Beta多樣性分析、物種組成及差異分析等。使用Mothur軟件,對各個樣品的Alpha多樣性指數(shù)進(jìn)行評估。使用QIIME軟件進(jìn)行Beta多樣性(Betadiversity)分析?;贐eta多樣性分析得到的距離矩陣,使用R語言工具繪制的PCoA圖。采用LefSe( Line Discriminant Analysis( LDA) Effect Size)在不同組間尋找具有統(tǒng)計(jì)學(xué)差異的Biomarker。根據(jù)各個物種在各個樣品中的豐度以及變化情況,進(jìn)行Sparcc相關(guān)分析并篩選相關(guān)性大于0.1且P值小于0.05的數(shù)據(jù)構(gòu)建相關(guān)性網(wǎng)絡(luò)。使用PICRUSt2軟件采用將待預(yù)測的特征序列與軟件中已有的系統(tǒng)發(fā)育樹中進(jìn)行物種注釋,使用IMG微生物基因組數(shù)據(jù)進(jìn)行功能信息的輸出進(jìn)而推測樣本中的功能基因組成,從而分析不同樣本或分組之間在功能上的差異。上述分析均在百邁客生物云網(wǎng)站平臺( https://developer-platform.biocloud.net/)完成。

    1.7數(shù)據(jù)統(tǒng)計(jì)與處理

    土壤理化性質(zhì)數(shù)據(jù)采用Excel和SPSS軟件進(jìn)行數(shù)據(jù)整理與圖表制作,用單因素方差分析( ANOVA)進(jìn)行差異顯著性分析,顯著性水平設(shè)定為P<0.05。采用Usearch軟件進(jìn)行操作分類單元(OTU)聚類分析。對于97%相似水平下的OTU進(jìn)行生物信息統(tǒng)計(jì)分析?;贠TU聚類分析以及各個水平分類下的物種對比結(jié)果,利用QIIME2軟件對OTU進(jìn)行多樣性分析?;诜诸悓W(xué)信息,在各個分類水平上進(jìn)行群落結(jié)構(gòu)的統(tǒng)計(jì)分析并統(tǒng)計(jì)各樣品組成。細(xì)菌16S rRNA基因V3~V4區(qū)采用SILVA細(xì)菌數(shù)據(jù)庫,使用樸素貝葉斯分類器對特征序列進(jìn)行分類學(xué)注釋;真菌ITSI區(qū)以UNITE為參考數(shù)據(jù)庫使用樸素貝葉斯分類器對特征序列進(jìn)行分類學(xué)注釋。

    2結(jié)果與分析

    2.1有機(jī)肥對田間香蕉枯萎病發(fā)病率的影響

    在大田試驗(yàn)中,挑選一塊香蕉種植常年發(fā)病區(qū)域,進(jìn)行巴西蕉病區(qū)種植試驗(yàn)。通過施用有機(jī)肥可以有效降低香蕉枯萎病發(fā)病率(圖1)。施用有機(jī)肥后,有效降低香蕉枯萎病發(fā)病率,使發(fā)病率維持在20%(表1),香蕉枯萎病發(fā)病率下降75%。對土壤致病菌進(jìn)行分離,得到3株形態(tài)疑似致病菌尖孢鐮刀菌菌株。通過TR4特異性引物擴(kuò)增,1號分離株與2號分離株中擴(kuò)增出大小為463 bp的片段,初步得出菌株1、2為當(dāng)?shù)叵憬犊菸≈虏【闠R4(圖2)。

    2.2有機(jī)肥對土壤部分理化性質(zhì)分析

    采集香蕉根際土壤,進(jìn)行理化性質(zhì)分析。發(fā)現(xiàn)施用有機(jī)肥植株土壤pH為5.8,未施用有機(jī)肥植株土壤pH為5.05,表明施用有機(jī)肥能使植株土壤pH值增加0.75個pH值,差異達(dá)到顯著水平(P<0.05)。而各土壤樣品在含水量方面均無顯著差異,為23%--26%。未施用有機(jī)肥植株土壤含氮量為0.16%,施用有機(jī)肥后,土壤含氮量為0.20%。未施用有機(jī)肥的植株土壤含磷量為0.08%,施用有機(jī)肥后,土壤含磷量為0.10%。施用有機(jī)肥的植株土壤同未施用有機(jī)肥植株在含氮量與含磷量之間的差異達(dá)到顯著水平(P<0.05)。施用有機(jī)肥的植株含鐵量顯著低于未施用有機(jī)肥的植株含鐵量(表1)。

    2.3香蕉植株根際土壤細(xì)菌與真菌擴(kuò)增子測序與物種多樣性分析

    所有樣品獲得總共459741條16S rRNA基因序列與471 682條真菌ITS基因序列,平均每份土壤中真菌與細(xì)菌序列數(shù)分別為78613條與76623條?;?7%的序列一致性分析,這些序列分別得到1707個細(xì)菌OTU與577個真菌OTU。所有樣品的真菌與細(xì)菌稀釋曲線表明隨著測序條數(shù)的增加,曲線趨于平緩,樣品測序充分(圖3)。

    對細(xì)菌與真菌結(jié)構(gòu)的豐富度(ACE和Chaol)和群落結(jié)構(gòu)多樣性(Shannon指數(shù))進(jìn)行了分析(表2),研究表明BIO處理后,相對于CK處理,真菌微生物群落的豐富度與多樣性顯著增加。然而,細(xì)菌群落結(jié)構(gòu)的豐富度(ACE和Chaol)并未有顯著性差異。

    2.4香蕉植株根際土壤細(xì)菌與真菌的Beta多樣性分析

    基于Bray-Curtis算法的主坐標(biāo)分析法(PCoA)可以揭示土壤微生物群落結(jié)構(gòu)之間的差異(圖4)。在土壤真菌群落方面,BIO處理同CK處理間表現(xiàn)出明顯差異。在所有變量數(shù)據(jù)中,最主要的2種因素分別可以影響70.68%與72.82%的細(xì)菌與真菌群落結(jié)構(gòu)。其中,最重要的因素( PCl)分別影響39.92%與41.02%的細(xì)菌與真菌群落結(jié)構(gòu)。

    2.5香蕉植株根際土壤細(xì)菌與真菌的物種豐度以及顯著性差異分析

    用土壤樣品中豐度前10的真菌與細(xì)菌門類來分析土壤微生物群落種群組成的變化(圖5)。變形菌門(Proteobacteria)、藍(lán)細(xì)菌門(Cyanobacteria)、酸桿菌門(Acidobacteria)、放線菌門(Actinobacteria)、綠灣菌門(Chloroflexi)和擬桿菌門(Bacteroidetes)是所有樣品中豐度最高的6個細(xì)菌門類,占各樣品所有測序序列80%以上。施用有機(jī)肥料后,在這些豐度較高的細(xì)菌門類中,藍(lán)細(xì)菌門(Cyanobacteria)和擬桿菌門(Bacteroidetes)豐度提升,而變形菌門(Pro teobacteria)、酸桿菌門(Acidobacteria)與放線菌門(Actinobacteria)的豐度降低,綠灣菌門(Chloroflexi)差異不顯著。在樣品真菌門類中,子囊菌門(Ascomycota)和擔(dān)子菌門(Basidiomycota)是所有樣品中豐度最高的兩種真菌門類,子囊菌門(Ascomycota)在施用有機(jī)肥后,豐度顯著增加,而擔(dān)子菌門(Basidiomycota)明顯下降。其余8個真菌門類,施用有機(jī)肥后未表現(xiàn)出顯著差異。

    在微生物屬水平差異方面,施用有機(jī)肥料后,比較豐度前20個屬的差異,結(jié)果表明施用有機(jī)肥料的根際土壤微生物群落中Coleofasciculaceae與SC-I-84多于未施用有機(jī)肥料土壤。而土壤中Paraburkholderia和Elsterales豐度降低。在土壤微生物群落中真菌前20個屬中,施用有機(jī)肥料后籃狀菌屬(Talaromyces)、Xenomyrothecium、Plectosphaerella、 殼多孢屬(Stagonospora)和圓孢霉屬( Staphylotrichurm)豐度提升顯著,然而未有真菌屬水平表現(xiàn)出顯著降低。鐮刀菌屬作為所有樣品中最主要的屬,在BioFOC處理與FOC處理中分別占10.87%與8.26%,兩者差異未達(dá)到顯著水平。

    用LEfSe分析進(jìn)行組間真菌與細(xì)菌各分類水平的差異顯著性分析(圖6)。BIO根際土壤中,在芽單孢菌門(Gemmatimonadetes)菌門中表現(xiàn)出最大的豐度差異。而且BIO根際土壤中,例如溶桿菌屬(Lvsobacter)、壤紅桿菌屬(Solirubrobacter)、黃色土源菌屬(Flavisolibacter)和稻田土壤菌屬(Oryzihumus)等菌屬的豐度也提升了。在真菌方面,BIO根際土壤中子囊菌屬(Ascomiycota)、Gibellulopsis、籃狀菌屬(Talaromyces)與Condenascus變得更加豐富(圖6-B)。

    2.6香蕉植株根際土壤細(xì)菌與真菌的物種間相關(guān)分析

    用Sparcc相關(guān)性進(jìn)行BIO處理與CK處理根際土壤中微生物群落真菌與細(xì)菌之間的相關(guān)性分析(表3)。細(xì)菌微生物群落分析了50個屬之間的相關(guān)性,結(jié)果表明BIO處理根際土壤微生物中,50個屬之間的相互抑制關(guān)系占據(jù)主要狀態(tài);而在CK處理根際土壤微生物中,50個屬之間的互相促進(jìn)關(guān)系呈現(xiàn)主要趨勢。土壤真菌微生物群落也進(jìn)行了50個屬之間相互關(guān)系的分析,發(fā)現(xiàn)BIO處理中,土壤真菌微生物屬之間表現(xiàn)出更加緊密的聯(lián)系,在BIO處理土壤中,木霉屬(Trichodema)與鐮刀菌屬(Fusarium)呈現(xiàn)出負(fù)相關(guān)關(guān)系,而未有真菌屬體現(xiàn)出與鐮刀菌屬(Fusarium)呈現(xiàn)正相關(guān)關(guān)系(表3)。在CK處理土壤中,葡萄穗霉屬(Stachybotrvs)與鐮刀菌屬(Fusarium)呈現(xiàn)出正相關(guān)關(guān)系。而在BIO組土壤中未有真菌屬與鐮刀菌屬(Fusarium)呈現(xiàn)出負(fù)相關(guān)關(guān)系。

    2.7香蕉植株根際土壤細(xì)菌與真菌的物種豐度與香蕉枯萎病病情指數(shù)相關(guān)性分析

    香蕉枯萎病是由于致病菌尖孢鐮刀菌(Fusariumoxyspoum f.sp.Cubense)引起的。在屬水平,表4展示了豐度前20的細(xì)菌和真菌屬。在細(xì)菌屬水平上,黃色桿菌屬(Xanthobacteraceae)與香蕉枯萎病發(fā)病顯示出負(fù)相關(guān)(P<0.05)。在真菌屬水平上,鐮刀菌屬(Fusarium)同香蕉枯萎病發(fā)病表現(xiàn)出正相關(guān)(P<0.05),而青霉菌屬(Penicillium)同香蕉枯萎病發(fā)病表現(xiàn)出負(fù)相關(guān)(P<0.05)。

    2.8有機(jī)肥料對根際土壤微生物群落功能影響

    通過KEGG代謝途徑的差異分析,可以觀測不同分組的樣品之間微生物群落的功能基因在代謝途徑上的差異和變化,是研究群落樣本為適應(yīng)環(huán)境變化發(fā)生的代謝功能改變的有效手段。分析得到的結(jié)果如圖7。涉及能量代謝( Energy metabolism)、輔助因子與維生素代謝( Metabolism of cofactors andvitamins)、復(fù)制與修復(fù)(Replication and repair)、核苷酸代謝( Nucleotide metabolism)與免疫系統(tǒng)(Immunesystem)的通路,在有機(jī)肥(BIO)處理中都得到增強(qiáng)。3討論與結(jié)論

    香蕉枯萎病在全球呈現(xiàn)出不斷蔓延的態(tài)勢,嚴(yán)重影響了我國香蕉產(chǎn)業(yè)的健康發(fā)展。在植物病害防治方面,主要存在化學(xué)防治、物理防治與生物防治。由于化學(xué)防治涉及生態(tài)環(huán)境問題,目前研究主要集中采用生物防治手段[23]。Zuo等[24]研究發(fā)現(xiàn)韭菜浸出液可以抑制TR4生長,能夠降低香蕉枯萎病病情。對病原菌具有有效防控的有機(jī)肥料,可以作為重要的生物防治劑[25]。因此,對有機(jī)肥的作用機(jī)制的研究可以有效地幫助我們更好地分析與改進(jìn)有機(jī)肥在農(nóng)業(yè)生產(chǎn)上的生物防治功能[26]。

    本研究中,在大田試驗(yàn)下,有機(jī)肥料的應(yīng)用顯著提高香蕉根際土壤pH、增加土壤中全氮與全磷含量。這與前人研究有機(jī)肥料可以改善土壤結(jié)構(gòu)與改良土壤理化性質(zhì)相一致[27]。有機(jī)肥中的添加菌印度梨形孢(Piriformospora indica)具有多種多樣功能,可以提高土壤中氮與磷的含量,促進(jìn)植物生長[28]。施用有機(jī)肥后,提高香蕉根際土壤全氮與全磷含量,供給香蕉根部更多營養(yǎng)物質(zhì)。有機(jī)肥料施用后,香蕉枯萎病發(fā)病率顯著降低,土壤中全鐵含量減少也達(dá)到顯著水平。一些對病原菌具有抑制效果的細(xì)菌,能夠分泌嗜鐵素,通過結(jié)合土壤中的Fe,進(jìn)而降低土壤中Fe含量,抑制病原菌生長,降低病害發(fā)生比率[29]。香蕉枯萎病病原菌尖孢鐮刀菌產(chǎn)生鐮刀菌酸加速其侵染香蕉根部,并造成香蕉球莖的腐爛[30]。pH4.0~5.5的偏酸性土壤適宜尖孢鐮刀菌繁殖,而pH為中性或以上土壤均不利于該菌生長[31]。土壤pH為酸性的條件下有利于尖孢鐮刀菌產(chǎn)生鐮刀菌酸,鐮刀菌酸的產(chǎn)生與菌絲生長呈現(xiàn)出正相關(guān)[32]。保持偏中性的土壤pH值,對香蕉枯萎病病原菌的繁殖與鐮刀菌酸的代謝產(chǎn)生限制作用,從而預(yù)防香蕉枯萎病的發(fā)生值得進(jìn)一步研究。本研究的有機(jī)肥料具有在香蕉枯萎病病區(qū)大田中穩(wěn)定土壤pH值以及降低土壤中Fe含量的作用。研究結(jié)果與之前發(fā)表文章一致[33],表明本研究的有機(jī)肥料可以通過穩(wěn)定香蕉根際土壤pH、減少全鐵含量、增加土壤中全氮與全磷含量有效降低香蕉枯萎病的發(fā)病率。

    微生物群落組成與結(jié)構(gòu)的研究發(fā)現(xiàn),施用有機(jī)肥后,顯著提高土壤微生物群落中真菌的豐富度與多樣性。微生物群落的豐富性與多樣性的提高,可以促使微生物群落具有更加穩(wěn)定的結(jié)構(gòu)去應(yīng)對逆境環(huán)境。在微生物群落相關(guān)性網(wǎng)絡(luò)圖中,施用有機(jī)肥后的真菌之間也具有更加緊密的聯(lián)系[34]。在屬水平,表4展示了豐度前20的屬。涉及與香蕉枯萎病發(fā)病相關(guān)的細(xì)菌屬中,黃色桿菌(Xanthobacteraceae)與香蕉枯萎病發(fā)病顯示出負(fù)相關(guān)。在真菌屬水平上,香蕉枯萎病是由于致病菌尖孢鐮刀菌(Fusariumoxyspoum f sp.cubense)引起的,分析結(jié)果顯示鐮刀菌屬(Fusarium)同香蕉枯萎病發(fā)病表現(xiàn)出正相關(guān)(P<0.05),而真菌結(jié)構(gòu)中青霉菌屬(Penicillium)同香蕉枯萎病發(fā)病表現(xiàn)出負(fù)相關(guān)(P<0.05)。施用有機(jī)肥料后,Penicillium的相對豐度提升。青霉菌屬(Penicillium)中真菌可以分泌一種新型抗真菌蛋白,具有生物防治應(yīng)用前景[35]。根際土壤微生物相關(guān)性分析結(jié)果表明,有機(jī)肥料處理后,真菌結(jié)構(gòu)中體現(xiàn)出木霉菌屬(Trichodema)與鐮刀菌屬(Fusarium)呈現(xiàn)出負(fù)相關(guān)關(guān)系,前人研究木霉菌屬(Trichodema)對于多種病原菌具有良好的抑制作用[36]。綜合來看,通過微生物結(jié)構(gòu)的分析,為我們提供了進(jìn)行下一步研究的方向,青霉菌屬(Penicillium)與木霉菌屬(Trichodema)能夠可以成為今后生物防控香蕉枯萎病研究的候選菌屬。

    本研究表明,有機(jī)肥料可以改良香蕉根際土壤部分理化性質(zhì),施用有機(jī)肥料植株土壤能使其土壤pH值增加0.75,增加土壤全氮與全磷含量,降低全鐵含量,有效降低香蕉枯萎病發(fā)病率,改善土壤微生物結(jié)構(gòu)。通過以下改善微生物結(jié)構(gòu)途徑來降低香蕉枯萎病發(fā)病率:(1)增加土壤中真菌的豐富度與多樣性;(2)提升土壤中對致病菌尖孢鐮刀菌的抑制菌種豐度;(3)促進(jìn)土壤中與香蕉枯萎病發(fā)病顯著負(fù)相關(guān)的菌種生長。

    參考文獻(xiàn):

    [1] LIN Y H,LIN Y J,CHANG T D,et al Developinent of a TaqManProbe-Based Insulated Isothermal Polymerase Chain Reaction (iiPCR)Assav for Detection of Fusarium oxysporum f.sp cubense Race 4[J]PLoS One. 2016. 11(7):1-13

    [2]11 x J,LI K,ZHOU D B,et al Biological control of banana wiltdisease caused bv Fusarium oxyspoum f.sp Cubense usingStreptomycessp H4[J].Biological Control, 2021. 155(lO):1-9

    [3] LIN Y J,LIN H K,LIN Y H Construction of Raman spectroscopicfingerprints for the detection of Fusarium wilt of banana iiiTaiwan [J] PLoS One. 2020. 15(3):1-14

    [4] WARMAN N M,AITKEN E A B The Movement of Fusoriumoxysporum f sp cubense (Sub-Tropical Race 4)in SusceptibleCultivars of Banana [J]. Frontiers in Plant Science. 2018.9: 17481757

    [5] 李進(jìn),張立丹,劉芳,等堿性肥料對香蕉枯萎病發(fā)生及上壤微生物群落的影響[J]植物營養(yǎng)與肥料學(xué)報(bào),2016. 22 (2):429436

    LI J. ZHANG L D LIU F.et al Effects of alkaline fertilizer onoccurrence of banana wilt disease and soil microbial community[J]Journal of Plant Nutrition and Fertilizer. 2016. 22(2):429436.(in Chinese)

    [6] 王榮,劉吉青,周海霞,等生物有機(jī)肥與保水劑對沒施連作黃瓜生長和上壤肥力的影響[J]河南農(nóng)業(yè)科學(xué),2018. 47 (8):45-53

    WANG R. LIU J Q,ZHOU H X,et al Effects of bio-organic fertilizerand water-retaining agent on cucumber growvth and soil fertility undercontmuous cropping in greenhouse[J].Journal of Henan AgriculturalSciences. 2018. 47(8)45-53( in Chinese)

    [7] ZHU z Y,TIAN z H,LI J x A Streptomyces morookaensis strainpromotes plant growth and suppresses Fusarium wilt of banana [J]Tropical Plont Pathology. 2021, 46(2):175-185

    [8] CHENG c z,LI D,QI Q,et al The root endophytic fungusSerendipita indica improves resistance of Banana to Fusoriumoxysporum f sp cubense tropical race 4 [J]. European Journal ofPlant Pathology, 2020, 156(1):87-100

    [9] LUO G w,LI L, FRIMAN v P,et al Organic ainendments increasecrop yields by improving microbe-mediated soil functioning ofagroecosysteins:A meta-analysis[J] Soil Biology and Biochemistry.2018. 124:105-115

    [IO] CHENG H Y,ZHANG D Q,HUANG B,et al.Organic fertilizerimproves soil fertility and restores the bactellal community after l,3一dichloropropene fumnigation[J].The Science of Total Environment,2020.738(6):1-42

    [11]陶成圓.含解淀粉芽孢桿菌NIN一6的生物有機(jī)肥防控香蕉枯萎病研究[D]南京:南京農(nóng)業(yè)大學(xué),2016

    TAo c Y Effect of the combination of bio-orrganic,fertilizer withBacillus amvloliauefaciens NJN一6 0n the coIltrol of baIlana Fisariumwlh dlsease[D] Nanjing:Nanjing A91lcultural Uhlversity,2016(inChlnese)

    [12]TAo c Y,LI R,xIoNG w,et al Blo-organic femlizers stimulate1ndigenous soil Pseudomonas populanons to enllanc.plant diseas.suppresslon[J] Microbiome,2020,8(1):137

    [13]劇虹伶辣椒一香蕉輪作聯(lián)合生物有機(jī)肥減輕高發(fā)枯萎病蕉園連作障礙機(jī)制研究[D]海口:海南大學(xué),2017

    Ju H L ReSedrch on the mechanisms of mitigate coontinuous obStdcleof banana orchord with serious with disense bycombmed pepper—banana rotation togemer with application fo bio-organic fertilier[D]Halkou:Hainan Unlversity,2017(in Chinese)

    [14]吉福桑,楊振,徐亞,等鹽脅迫下巴西蕉葉片的轉(zhuǎn)錄組和蛋白質(zhì)組關(guān)聯(lián)分析[J]分子植物育種,2020,18(23):7671-7678

    JI F S,YANG Z,XU Y,et al Associatlon analysls on transcriptomlcsand pmteomlcs ofMusa paradisiaca banaIla leafunder salt stless[J]Molecular Plant Breeding,2020,18(23) 7671 7678(inChinese)

    [15]Guo L J,YANG L Y,LIANG c c,et al Dlfferentlal colonlzatlonPatterns of Bananas(Musa spp)by P11yslol091cal Race l and Race 4Isolates of Fussing arysparum f.sp cubense[J].Journal ofPhytopathology,2015,163(lO):807-817

    [16]FAN H Y,LEI z x,DoNG H H,et al Immune lesponses ln BrazilianbaIlaIla determiIlmg the. pathogenlc difkrences between the。physlologlcal mces l and 4 of Fusarium axysporum f sp cubense[J]Journal of Plant Pathology,2019,10l(2)225-234

    [17]劉光崧上壤理化分析與剖面描述[M].北京:中國標(biāo)準(zhǔn)出版社,1996:50-8l

    [18]MosTERT D,MoLINA A B,DANIELLs J,et al The distributionand 110st range of the banana Fusarium wilt fungus,F(xiàn)usariumoxysporum f.sp cubense,m Asla[J].PLoS One,2017,l2(7):l-lO

    [19]DITA M A,wAALwIJK c,BuDDENHAGEN l w,et al Amolecular diagnostic for tropical race 4 0f the banana Fusarium wiltpathogen[J].Plant,Pathology,2010,59(2):348-357

    [20]MAGOC T,sALzBERG s L FLAsH:fast length adjustment of shortreads t0 implove. genome assemblles[J] Bioinformatics,2011,27(21):2957-2963

    [21]BoLGER A M,LoHsE M,usADEL B Trimmuomatlc:a fiexibletrinmer for Illumina sequence data[J] Bioinformatics(Oxford,England),2014,30(15):2114-2120

    [22]EDGAR R c uPARsE:bighly accurate oTu s.quences flommicrobial amplicon reads[J].Nature Methods,2013,10(lO):996-998

    [23]cHEN D, LrLJ x, LI c Y, et al Isolatlon of Bacillusamyloliquefaciens S20 and lt8 applicatioIl ln coIltrol of.ggplantbactenal wilt[J].Journal of Environmental Management,2014,137(1):l20-127

    [24]zuo c w,LI c Y,LI B,et al The toxlc luechanlsln and bloactⅣecomponents of Chlnese leek root exudates acting against Fusariumoxysporum f sp cubense troplcal race 4[J].European Journal ofPlant Pathology,2015,l43(3):447-460

    [25]cHowDHuRY s P,HARTMANN A,GAo x w,et al.Biocontrolmechanlsm by root—assoclated Bacillus amyloliquefaciens FZB42一arevlew[J].Frontiers in Microbiology 2015,6(7):780-790

    [26]曹云soR 9微生物有機(jī)肥防治黃瓜上傳枯萎病的效應(yīng)與機(jī)制研究[D]南京:南京農(nóng)業(yè)大學(xué),2011

    CAO Y.Control of Fusarium wilt dlsease of cucumber by applicatlonof blo-organlc feltlllzer aIld lts workmg mechanlsm[D] NanjingNanjing Agriculrural Unlverslty,20l l(in Chinese)

    [27]KAI T,TAMAKI M Effect of organlc and chemlcal fernllzer8ppllcatlon on growth,yleld,aIld soll blochemlcal propemes ofLandrace Brassica napus L.leaf-and-stem vgetable and Landrace(norabona)[J] Journal of Agricultural Chemistry and Environment.2020.9(4):314-330

    [28]袁聽,覃金蘭,何珊,等印度梨形孢對植物生長調(diào)節(jié)劑在水稻上過量使用的緩解作用[J].福建農(nóng)業(yè)學(xué)報(bào),2020,35(4):398-405

    YUAN T,QIN J L.HE S,et al Mmgatmg ill—efkct of plant growthregulator ovemse on ric。plallts by Piriformospora indica [J].FujianJournal of Agricultural Sciences,2020,35(4):398 405(inChinese)

    [29]Gu s H,wEI z,sHAo z Y,et al.competltlon for iron drlvesphytopathogen control by nanlral lhlzosphere mlcroblomes[J]Nature Microbiology,2020,5(8):1002-10lO

    [30]LI c Y,zuo c w,DENG G M,et al contalnlnatlon ofbananas wlthbeauvellcm aIld fusaric acld produced by Fusarium Oxysporum F.SPcubense[J] PLos One,2013,8(7):l-ll

    [31]彭雙,王一明,葉旭紅,等上壤環(huán)境因素對致病性尖孢鐮刀菌生長的影響[J]上壤,2014. 46 (5):845850

    PENG S WANG Y M. YE X H et al Effects of Soil Habitat Factorson Growth of Fusarium oxysporum f sp. niveum and Fusariumoxysporum f sp cucumerinum [J]. Soils. 2014. 46 (5):845-850(in Chinese)

    [32]馬圍斌,林德佩,王葉筠,等培養(yǎng)條件對西瓜枯萎病菌鐮刀菌酸產(chǎn)生的影響[J]植物病理學(xué)報(bào),1996. 26 (2):92-96

    MA G B,LIN D P,WANG Y J,et al Effect of cultural conditions onproduction of fusaric acid of watermelon Fusorium wilt fungus [J]Acta Phytopathologica Sinica, 1996. 26(2):92-96( in Chinese)

    [33]薛超香蕉根際上壤微生物區(qū)系特征與上傳枯萎病防控研究[D]南京:南京農(nóng)業(yè)大學(xué),2015

    XUE C Manipulation of microbial community in banana thizosphererto suppress Fusarium wilt of banana[D]. Nanjing: NanjingAgricultural University 2015 (in Chinese)

    [34]柳凱,季倩茹,陳靜,等施用Streptomyces alfalfae XY25T對根腫病上壤性質(zhì)及微生物群落的影響[J]微生物學(xué)通報(bào),2020.47 (1):97-108

    LIU K,JI Q R,CHEN J,et al Effect of Streptomy'ces alfalfae XY25Ton soil properties and microflora iii clubroot-diseased soil [J]Microbiology China. 2020. 47(1):97-108 (in Chinese)

    [35] CHEN z T,AO J Q,YANG w c,et al Purification andcharacterization of a novel antifungal protein secreted bv Penicilliumchrvsogenum from an Arctic sediment [J] Applied Microbiology andBiotechnology. 2013. 97( 24): 10381-10390

    [36] NGUVO K J,GAO x Q Weapons hidden undemeath: Bio-controlagents and their potentials to activate plant induced svstemic resistancein controlling crop Fusarium diseases [J].Journal of Plant Diseasesand Protection. 2019. 126(3):177-190

    (責(zé)任編輯:黃愛萍)

    收稿日期:2020-11-19初稿;2021-03-16修改稿

    作者簡介:朱志炎(1989-),男,博士研究生,研究方向:植物病理(E-mail: 420905220@qq.com)

    *通信作者:李建雄(1969-),男,博士,研究員,研究方向:植物病理(E-inail: jxli@scbgaccn)

    基金項(xiàng)目:中國科學(xué)院STS區(qū)域重點(diǎn)項(xiàng)目(KFJ-STS-QYZX-044);廣東省農(nóng)業(yè)廳委托研究項(xiàng)日(Y434121002)

    猜你喜歡
    高通量測序枯萎病有機(jī)肥
    辣椒枯萎病生防木霉菌T21的分離鑒定及其生物學(xué)特性研究
    黃瓜穴盤苗期枯萎病抗性鑒定方法及枯萎病脅迫下的生理響應(yīng)
    4種農(nóng)家有機(jī)肥對當(dāng)歸增產(chǎn)效果試驗(yàn)初報(bào)
    覆蓋防草布的藍(lán)莓樹施用有機(jī)肥技術(shù)
    油頁巖灰渣有機(jī)肥在白樺地的試驗(yàn)研究
    川明參輪作對煙地土壤微生物群落結(jié)構(gòu)的影響
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    人參根際真菌群落多樣性及組成的變化
    LncRNAs作為miRNA的靶模擬物調(diào)節(jié)miRNA
    不同西瓜品種枯萎病抗性鑒定方法比較
    日本欧美视频一区| 亚洲精品自拍成人| 亚洲婷婷狠狠爱综合网| 国产精品三级大全| 日韩制服骚丝袜av| 日本-黄色视频高清免费观看| 日本黄色片子视频| 一边亲一边摸免费视频| 国产精品蜜桃在线观看| 国产成人精品久久久久久| 亚洲欧美一区二区三区国产| 一区在线观看完整版| 高清在线视频一区二区三区| 丰满迷人的少妇在线观看| 午夜视频国产福利| 青春草亚洲视频在线观看| 国产成人freesex在线| 亚洲高清免费不卡视频| 午夜免费男女啪啪视频观看| 亚洲成人手机| 国产一区有黄有色的免费视频| 久久久久网色| 我要看日韩黄色一级片| 国产一区二区三区av在线| 午夜视频国产福利| 在线免费观看不下载黄p国产| 久久综合国产亚洲精品| 亚洲精品乱久久久久久| 国产精品一区二区性色av| 国产精品蜜桃在线观看| 美女内射精品一级片tv| 亚洲色图综合在线观看| 久久影院123| 亚洲精品国产成人久久av| 国产亚洲午夜精品一区二区久久| 日韩三级伦理在线观看| 亚洲自偷自拍三级| av免费在线看不卡| 美女大奶头黄色视频| 国产免费又黄又爽又色| 高清欧美精品videossex| 日韩强制内射视频| 亚洲av福利一区| 99热网站在线观看| 欧美日韩在线观看h| 不卡视频在线观看欧美| 蜜桃在线观看..| 丝袜喷水一区| 精品亚洲成国产av| 丁香六月天网| 久久精品国产a三级三级三级| 色视频www国产| 成年人午夜在线观看视频| 男人和女人高潮做爰伦理| 久久97久久精品| 免费大片18禁| 美女xxoo啪啪120秒动态图| 亚洲精品国产成人久久av| 美女国产视频在线观看| 欧美激情国产日韩精品一区| 日韩 亚洲 欧美在线| 精品久久国产蜜桃| 亚洲av欧美aⅴ国产| 亚洲国产精品国产精品| 青青草视频在线视频观看| 国产一区二区三区av在线| 久久久久久久久久久免费av| 三上悠亚av全集在线观看 | 国产伦精品一区二区三区四那| 日本黄色片子视频| 免费观看性生交大片5| 97在线视频观看| 精品国产乱码久久久久久小说| 免费人成在线观看视频色| 51国产日韩欧美| 美女xxoo啪啪120秒动态图| av一本久久久久| 精品久久久久久久久av| 自拍偷自拍亚洲精品老妇| 欧美日韩综合久久久久久| 亚洲美女搞黄在线观看| 高清欧美精品videossex| 国产综合精华液| 在线观看av片永久免费下载| 看十八女毛片水多多多| 亚洲久久久国产精品| 国产色婷婷99| 麻豆精品久久久久久蜜桃| 男男h啪啪无遮挡| 美女国产视频在线观看| 自拍偷自拍亚洲精品老妇| 国产在线男女| 蜜桃久久精品国产亚洲av| 五月天丁香电影| 亚洲成人手机| 综合色丁香网| 国产精品国产三级国产av玫瑰| 日韩精品免费视频一区二区三区 | 久久99热这里只频精品6学生| av福利片在线观看| 国产乱来视频区| 菩萨蛮人人尽说江南好唐韦庄| 日韩视频在线欧美| 在线观看美女被高潮喷水网站| 一个人看视频在线观看www免费| 女人精品久久久久毛片| 亚洲精品中文字幕在线视频 | 极品少妇高潮喷水抽搐| 国产熟女午夜一区二区三区 | 亚洲真实伦在线观看| 又爽又黄a免费视频| 如日韩欧美国产精品一区二区三区 | 亚洲婷婷狠狠爱综合网| 免费高清在线观看视频在线观看| 日本欧美视频一区| 午夜福利影视在线免费观看| 亚洲欧美成人综合另类久久久| 日本黄大片高清| 午夜免费男女啪啪视频观看| 亚洲av成人精品一区久久| 观看免费一级毛片| 男人狂女人下面高潮的视频| 亚洲人与动物交配视频| 最后的刺客免费高清国语| 日韩精品有码人妻一区| h日本视频在线播放| 在现免费观看毛片| 一区二区三区四区激情视频| 国产高清国产精品国产三级| 各种免费的搞黄视频| 久久6这里有精品| 国产欧美亚洲国产| 一个人免费看片子| 亚洲av在线观看美女高潮| 九九爱精品视频在线观看| 中文字幕亚洲精品专区| 亚洲人成网站在线播| a 毛片基地| 亚洲,一卡二卡三卡| 国产精品人妻久久久影院| 日韩中字成人| 一区在线观看完整版| 这个男人来自地球电影免费观看 | 久久99一区二区三区| 国产在线免费精品| 免费看日本二区| 国产精品久久久久久精品电影小说| av福利片在线| 欧美丝袜亚洲另类| 视频区图区小说| 国产精品一二三区在线看| 精品久久久久久电影网| 国产欧美亚洲国产| 亚洲成人一二三区av| 91精品伊人久久大香线蕉| 免费av不卡在线播放| 一级毛片我不卡| 男的添女的下面高潮视频| 少妇精品久久久久久久| 国产成人a∨麻豆精品| 精品久久久久久久久亚洲| 精品人妻一区二区三区麻豆| 国产老妇伦熟女老妇高清| 国国产精品蜜臀av免费| 亚洲美女黄色视频免费看| 国产精品久久久久久久电影| 午夜免费观看性视频| 国产精品伦人一区二区| 亚洲精品国产av成人精品| 欧美日韩国产mv在线观看视频| 人人妻人人澡人人看| 日本免费在线观看一区| 久久精品熟女亚洲av麻豆精品| 亚洲av综合色区一区| 一区二区三区乱码不卡18| 91久久精品国产一区二区三区| 天天操日日干夜夜撸| 日韩在线高清观看一区二区三区| 狠狠精品人妻久久久久久综合| 日韩免费高清中文字幕av| 国产伦精品一区二区三区视频9| 永久网站在线| 日日摸夜夜添夜夜添av毛片| 啦啦啦啦在线视频资源| av国产久精品久网站免费入址| 日本猛色少妇xxxxx猛交久久| av播播在线观看一区| 免费av不卡在线播放| 亚洲人成网站在线观看播放| 久久久久久久国产电影| 五月天丁香电影| 国产精品伦人一区二区| 国产黄色免费在线视频| 久久毛片免费看一区二区三区| 亚洲av在线观看美女高潮| 中文在线观看免费www的网站| 国产一区有黄有色的免费视频| 在线观看一区二区三区激情| 欧美xxxx性猛交bbbb| a 毛片基地| 免费黄色在线免费观看| 卡戴珊不雅视频在线播放| 国产一区亚洲一区在线观看| 免费久久久久久久精品成人欧美视频 | 日本免费在线观看一区| 国产精品一区二区在线观看99| 午夜福利,免费看| 久久久久久久久久成人| 国产欧美亚洲国产| 熟女av电影| 国产午夜精品一二区理论片| 在线观看免费视频网站a站| av视频免费观看在线观看| 精品一品国产午夜福利视频| 免费看av在线观看网站| 国产 精品1| 好男人视频免费观看在线| 丰满迷人的少妇在线观看| 亚洲内射少妇av| 丰满人妻一区二区三区视频av| 美女主播在线视频| 久久99精品国语久久久| 黄色配什么色好看| 大香蕉久久网| 人妻制服诱惑在线中文字幕| 国产高清三级在线| 欧美日韩亚洲高清精品| 女性被躁到高潮视频| 六月丁香七月| 丰满迷人的少妇在线观看| 国产国拍精品亚洲av在线观看| 人人妻人人澡人人看| 国产一区二区三区av在线| 美女大奶头黄色视频| 久久久国产精品麻豆| 日韩成人伦理影院| 日韩视频在线欧美| 亚洲av在线观看美女高潮| 99热这里只有是精品在线观看| av女优亚洲男人天堂| 日日爽夜夜爽网站| 最近中文字幕高清免费大全6| 精品少妇内射三级| 久久久久精品性色| 91精品国产九色| 成人美女网站在线观看视频| 精品一品国产午夜福利视频| 伊人久久精品亚洲午夜| 国产有黄有色有爽视频| 精品人妻偷拍中文字幕| 22中文网久久字幕| 午夜视频国产福利| 亚洲国产日韩一区二区| 全区人妻精品视频| 成年美女黄网站色视频大全免费 | videossex国产| 精品久久久久久久久av| 六月丁香七月| av福利片在线观看| 美女中出高潮动态图| 新久久久久国产一级毛片| 午夜影院在线不卡| 国产淫片久久久久久久久| 国产精品99久久99久久久不卡 | 桃花免费在线播放| av在线观看视频网站免费| 成人黄色视频免费在线看| 交换朋友夫妻互换小说| 国产一区二区三区av在线| 丰满人妻一区二区三区视频av| 精华霜和精华液先用哪个| 欧美区成人在线视频| 国产淫语在线视频| 午夜影院在线不卡| 人妻制服诱惑在线中文字幕| 亚洲一区二区三区欧美精品| 男的添女的下面高潮视频| 欧美精品人与动牲交sv欧美| 亚洲精品aⅴ在线观看| 天堂中文最新版在线下载| 欧美xxⅹ黑人| 亚洲av欧美aⅴ国产| 国产女主播在线喷水免费视频网站| 国产免费福利视频在线观看| 永久网站在线| 99视频精品全部免费 在线| 伊人久久精品亚洲午夜| 亚洲国产精品一区二区三区在线| 汤姆久久久久久久影院中文字幕| 久久久久久伊人网av| 国产成人精品婷婷| 人人妻人人澡人人看| 国产一区二区在线观看日韩| 欧美日韩精品成人综合77777| 国产 一区精品| 五月伊人婷婷丁香| 免费播放大片免费观看视频在线观看| 美女主播在线视频| 婷婷色综合大香蕉| 3wmmmm亚洲av在线观看| 少妇被粗大猛烈的视频| 日本午夜av视频| 国内揄拍国产精品人妻在线| 午夜av观看不卡| 美女主播在线视频| 日本午夜av视频| 免费播放大片免费观看视频在线观看| www.色视频.com| 国产男女内射视频| 亚洲精品国产色婷婷电影| 成人影院久久| 成人漫画全彩无遮挡| 国产高清国产精品国产三级| 亚洲国产精品一区三区| 亚洲精品久久午夜乱码| 国产精品一区www在线观看| videos熟女内射| 国产毛片在线视频| 男女免费视频国产| 亚洲欧洲国产日韩| 日韩中文字幕视频在线看片| 午夜91福利影院| 国产乱来视频区| 久久久久久久亚洲中文字幕| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品一区蜜桃| 日韩视频在线欧美| 人人妻人人澡人人爽人人夜夜| 亚洲av欧美aⅴ国产| 婷婷色麻豆天堂久久| 亚洲美女视频黄频| 亚洲自偷自拍三级| 国产精品一区二区在线不卡| 日日爽夜夜爽网站| 熟女人妻精品中文字幕| 国产老妇伦熟女老妇高清| 久久久国产欧美日韩av| 久久精品国产亚洲av涩爱| 天天操日日干夜夜撸| av免费观看日本| 国产黄色免费在线视频| 日本午夜av视频| 色吧在线观看| 美女主播在线视频| 精品久久久久久久久亚洲| 亚洲国产av新网站| 丰满乱子伦码专区| 丝袜脚勾引网站| 免费黄频网站在线观看国产| 精品亚洲成a人片在线观看| 精品亚洲成a人片在线观看| 久久久久国产网址| 午夜日本视频在线| 国产成人91sexporn| 一级二级三级毛片免费看| 青春草国产在线视频| 欧美97在线视频| 高清视频免费观看一区二区| 亚洲性久久影院| 成年美女黄网站色视频大全免费 | 狠狠精品人妻久久久久久综合| 秋霞伦理黄片| 免费高清在线观看视频在线观看| 日韩中文字幕视频在线看片| 丰满乱子伦码专区| 一区二区三区乱码不卡18| 国产黄色视频一区二区在线观看| 蜜桃久久精品国产亚洲av| av专区在线播放| 99精国产麻豆久久婷婷| 久久99蜜桃精品久久| 精品久久国产蜜桃| 天美传媒精品一区二区| 美女福利国产在线| 男人爽女人下面视频在线观看| 下体分泌物呈黄色| 免费黄网站久久成人精品| 亚洲人成网站在线播| 丝瓜视频免费看黄片| 一本久久精品| 欧美 亚洲 国产 日韩一| 亚洲色图综合在线观看| 精品少妇内射三级| 51国产日韩欧美| 晚上一个人看的免费电影| 国产精品伦人一区二区| av在线老鸭窝| 亚洲国产毛片av蜜桃av| 国产伦理片在线播放av一区| 色94色欧美一区二区| 最近中文字幕2019免费版| 寂寞人妻少妇视频99o| 亚洲精品亚洲一区二区| 国产 一区精品| 国产伦精品一区二区三区视频9| 桃花免费在线播放| 国产成人免费观看mmmm| 丰满乱子伦码专区| 亚洲,一卡二卡三卡| 性色avwww在线观看| 国产精品久久久久久精品电影小说| 免费在线观看成人毛片| 亚洲国产成人一精品久久久| 久久精品国产鲁丝片午夜精品| 精品视频人人做人人爽| 亚洲,欧美,日韩| 插逼视频在线观看| 91午夜精品亚洲一区二区三区| 观看av在线不卡| 日韩中文字幕视频在线看片| 午夜日本视频在线| 精品国产乱码久久久久久小说| 大码成人一级视频| 国产精品女同一区二区软件| 18禁在线无遮挡免费观看视频| 亚洲四区av| 国产亚洲一区二区精品| 精品99又大又爽又粗少妇毛片| 男女无遮挡免费网站观看| 麻豆成人av视频| 国产女主播在线喷水免费视频网站| 日本wwww免费看| 青青草视频在线视频观看| 亚洲人与动物交配视频| 热re99久久精品国产66热6| 亚洲三级黄色毛片| 日韩成人伦理影院| 久久热精品热| 国产一区亚洲一区在线观看| 亚洲婷婷狠狠爱综合网| 亚洲美女搞黄在线观看| 日韩在线高清观看一区二区三区| 国产伦理片在线播放av一区| 在线天堂最新版资源| 看非洲黑人一级黄片| 9色porny在线观看| 丝袜脚勾引网站| 成年女人在线观看亚洲视频| 高清午夜精品一区二区三区| 91aial.com中文字幕在线观看| 18禁裸乳无遮挡动漫免费视频| 高清av免费在线| 少妇猛男粗大的猛烈进出视频| 麻豆精品久久久久久蜜桃| 老女人水多毛片| 特大巨黑吊av在线直播| 亚洲无线观看免费| a级毛片在线看网站| 中国三级夫妇交换| 色哟哟·www| 中文字幕精品免费在线观看视频 | 国产亚洲av片在线观看秒播厂| 日本91视频免费播放| 日本欧美视频一区| 精品国产露脸久久av麻豆| 日韩制服骚丝袜av| 成人毛片a级毛片在线播放| 少妇人妻久久综合中文| 2021少妇久久久久久久久久久| tube8黄色片| 日本午夜av视频| 美女内射精品一级片tv| 99热6这里只有精品| 亚洲国产日韩一区二区| 国产一区二区三区av在线| 免费久久久久久久精品成人欧美视频 | 免费黄色在线免费观看| 久久狼人影院| 美女脱内裤让男人舔精品视频| 亚洲av.av天堂| 校园人妻丝袜中文字幕| 高清毛片免费看| 天天操日日干夜夜撸| 国产日韩一区二区三区精品不卡 | 黄色日韩在线| 日韩人妻高清精品专区| 最近手机中文字幕大全| 一级毛片我不卡| 建设人人有责人人尽责人人享有的| 亚洲图色成人| 欧美最新免费一区二区三区| 少妇人妻 视频| 天美传媒精品一区二区| 国产精品伦人一区二区| 少妇的逼水好多| 国内揄拍国产精品人妻在线| 王馨瑶露胸无遮挡在线观看| 日韩一区二区视频免费看| 亚洲精品视频女| 免费观看a级毛片全部| 日韩av免费高清视频| 免费观看性生交大片5| 涩涩av久久男人的天堂| 91久久精品电影网| 伦理电影免费视频| 啦啦啦啦在线视频资源| 黑丝袜美女国产一区| 少妇被粗大猛烈的视频| 人妻人人澡人人爽人人| 一级毛片aaaaaa免费看小| 婷婷色麻豆天堂久久| av卡一久久| 国产欧美日韩精品一区二区| 全区人妻精品视频| 精品午夜福利在线看| 国产精品一区www在线观看| 伊人久久精品亚洲午夜| 亚洲伊人久久精品综合| 国产黄片视频在线免费观看| 新久久久久国产一级毛片| 欧美精品一区二区大全| 青青草视频在线视频观看| 亚洲国产成人一精品久久久| 少妇被粗大的猛进出69影院 | 深夜a级毛片| 色吧在线观看| 少妇猛男粗大的猛烈进出视频| 久久久国产欧美日韩av| a级毛片在线看网站| 最近中文字幕高清免费大全6| 精品人妻一区二区三区麻豆| 精品人妻偷拍中文字幕| 亚洲情色 制服丝袜| 国产无遮挡羞羞视频在线观看| 三级经典国产精品| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美中文字幕日韩二区| 熟女av电影| 人人妻人人澡人人看| 国产美女午夜福利| 亚洲国产精品999| 婷婷色综合大香蕉| 亚洲精品久久午夜乱码| 一本大道久久a久久精品| 丝瓜视频免费看黄片| 国产有黄有色有爽视频| 亚洲精品亚洲一区二区| av在线观看视频网站免费| 中文精品一卡2卡3卡4更新| 中文字幕精品免费在线观看视频 | 亚洲精品亚洲一区二区| 日韩三级伦理在线观看| 建设人人有责人人尽责人人享有的| 久久鲁丝午夜福利片| 亚洲精品自拍成人| 国精品久久久久久国模美| 在线观看免费视频网站a站| 哪个播放器可以免费观看大片| 色视频在线一区二区三区| 国产亚洲av片在线观看秒播厂| 美女大奶头黄色视频| 日日摸夜夜添夜夜爱| av在线老鸭窝| 国产黄色免费在线视频| 午夜91福利影院| 国产精品久久久久久久久免| 成人免费观看视频高清| 欧美精品国产亚洲| 六月丁香七月| 国产午夜精品一二区理论片| 日本wwww免费看| 免费黄色在线免费观看| 亚洲久久久国产精品| 日韩三级伦理在线观看| 国产日韩欧美视频二区| 少妇高潮的动态图| 五月玫瑰六月丁香| 免费看av在线观看网站| 亚洲在久久综合| 国产精品人妻久久久影院| 水蜜桃什么品种好| 国产精品.久久久| 五月玫瑰六月丁香| 免费看日本二区| 免费观看av网站的网址| 最近中文字幕高清免费大全6| 久久精品熟女亚洲av麻豆精品| 最新中文字幕久久久久| 亚洲欧美日韩另类电影网站| 日韩欧美一区视频在线观看 | 97在线视频观看| 中文在线观看免费www的网站| 国产免费一级a男人的天堂| 91在线精品国自产拍蜜月| 丁香六月天网| 黄片无遮挡物在线观看| 最近手机中文字幕大全| 久久精品国产亚洲av天美| av在线老鸭窝| 亚洲激情五月婷婷啪啪| 国产黄片视频在线免费观看| 国产精品一二三区在线看| 亚洲av综合色区一区| 精品一区二区三卡| 人体艺术视频欧美日本| 日韩中文字幕视频在线看片| 水蜜桃什么品种好| 国产伦在线观看视频一区| 免费观看a级毛片全部| av卡一久久| 乱码一卡2卡4卡精品| 人体艺术视频欧美日本| 黑人高潮一二区| 精品卡一卡二卡四卡免费| 中文字幕免费在线视频6| 国产一区二区三区综合在线观看 | 欧美日韩一区二区视频在线观看视频在线| 成人毛片60女人毛片免费| 国产精品免费大片| 国产日韩欧美视频二区| 国产在线视频一区二区| 少妇人妻一区二区三区视频| 国产免费视频播放在线视频| 精品亚洲成国产av| 在线看a的网站|