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

    健康與罹患青枯病的番茄土壤細菌群落特征比較①

    2021-03-23 07:47:56沈宗專操一凡王東升劉紅軍沈其榮
    土壤 2021年1期
    關鍵詞:土傳感病青枯病

    沈宗專,黃 炎,操一凡,王東升, 2,劉紅軍,李 榮*,沈其榮

    健康與罹患青枯病的番茄土壤細菌群落特征比較①

    沈宗專1,黃 炎1,操一凡1,王東升1, 2,劉紅軍1,李 榮1*,沈其榮1

    (1 南京農業(yè)大學資源與環(huán)境科學學院,江蘇省固體有機廢棄物資源化高技術研究重點實驗室,江蘇省有機固體廢棄物資源化協(xié)同創(chuàng)新中心,國家有機類肥料工程技術研究中心,作物免疫學重點實驗室,南京 210095;2 南京市蔬菜科學研究所,南京 210042)

    應用實時熒光定量PCR及MiSeq高通量測序技術,全面地研究了連作番茄田塊中健康與感染青枯病植株周圍土體及根際土壤細菌群落結構和組成。結果表明:健康番茄土體土壤的pH及全碳含量顯著高于感病番茄土體土壤;土體及根際土壤的細菌群落結構和組成明顯不同于感病番茄土體及根際土壤細菌群落。與感病番茄根際相比,健康番茄根際細菌的數量顯著升高而青枯菌數量顯著降低;細菌群落的Shannon多樣性指數顯著增高;擬桿菌門及其所含的噬幾丁質菌屬、金桿菌屬、動桿菌屬、黃桿菌屬及的相對豐度顯著增高而變形菌門及其所含的青枯菌屬的相對豐度顯著降低。綜上,抑制土傳青枯病發(fā)生的番茄根際土壤細菌群落特征明顯,其生物量及多樣性高,土著有益菌群數量多而病原菌數量少,為番茄土傳青枯病的生物防控提供了指導方向與理論依據。

    番茄青枯??;抑病型土壤;高通量測序;細菌群落多樣性;細菌群落組成

    土壤生態(tài)系統(tǒng)是連接大氣圈、水圈、巖石圈和生物圈的紐帶,與地上生態(tài)系統(tǒng)之間緊密聯(lián)系、相互依存,共同決定著陸地生態(tài)系統(tǒng)的特征、過程和功能[1]。土壤微生物作為土壤生態(tài)系統(tǒng)的主要組成部分,不僅在元素的生物地球化學循環(huán)中起著重要作用,也在作物健康生長和可持續(xù)生產中起著關鍵作用[2]。盡管植物與微生物的互作關系是土傳病害發(fā)生的重要決定因子[3],但目前對于根際微生物區(qū)系組成與生態(tài)系統(tǒng)功能之間的關聯(lián)仍有待深入研究。

    近年來農用化學品的大量投入是作物持續(xù)高產的重要支撐之一,但同時也給土壤及環(huán)境的健康帶來了挑戰(zhàn)[4]。高強度集約化農業(yè)生產一定程度上破壞了土壤微生物群落結構的穩(wěn)態(tài)平衡,削弱了土壤生態(tài)系統(tǒng)的功能。尤以抑制土傳病害發(fā)生的能力(抑病能力)下降為甚,導致作物根腐病、枯萎病、全蝕病、立枯病等土傳病害頻繁爆發(fā)[5]。番茄是全世界栽培最為普遍的果菜之一,設施番茄產業(yè)已成為農民增收的重要產業(yè)之一。然而,因青枯菌()侵染導致的番茄青枯病嚴重威脅著世界番茄種植產業(yè)的發(fā)展[6]。因此,土傳青枯病的有效防控仍是關系到世界番茄產業(yè)可持續(xù)發(fā)展的重大理論及技術問題。

    土壤微生物群落結構及組成對于番茄抵御青枯菌的侵染有著重要影響[7],因此,研究同一連作田塊中健康與罹患青枯病的番茄土壤微生物群落結構及組成差異對于番茄土傳青枯病的防控極其重要。然而,盡管已有大量研究表明調控土壤微生物區(qū)系可有效防控番茄青枯病的發(fā)生[6, 8-9],但對于同一連作條件下健康與患病番茄土壤微生物群落結構與組成差異的研究仍相對缺乏。楊尚東等[10]利用傳統(tǒng)的分子生態(tài)學手段(PCR-DGGE)初步比較了番茄青枯病罹病植株和健康植株根際土壤細菌群落結構及多樣性特征,但該研究手段不能揭示群落的組成差異;Li等[5]利用454焦磷酸測序技術研究比較了健康與患病番茄根際土壤及根內細菌群落結構和組成特征,但未揭示與根際細菌裝配密切相關的土體土壤細菌群落差異。

    因此,為全面揭示連作條件下健康與患病番茄的土體及根際細菌群落結構和組成差異及探究土體細菌群落的差異與根際細菌群落組成差異之間的關系,本研究擬采集設施大棚內長期連作番茄田塊中感染土傳青枯病植株的土體及根際土壤和健康植株的土體及根際土壤,利用熒光定量PCR 和 MiSeq 高通量測序技術比較設施大棚內健康番茄與感染土傳青枯病番茄的土體及根際土壤中細菌群落結構和組成差異,為番茄青枯病土壤的微生物區(qū)系調控提供理論依據和方向支撐。

    1 材料與方法

    1.1 供試材料

    1.1.1 供試種子 供試番茄種子為世紀紅冠(Mill., CV. Shijihongguan),由南京市蔬菜花卉科學研究所提供。

    1.1.2 試驗點基本信息 試驗點設在江蘇省南京市蔬菜花卉科學研究所(31°43′N,118°46′E)。試驗區(qū)域為典型亞熱帶季風氣候,年平均溫度為15.4 °C,年均降雨量達1 106 mm。試驗設施土壤的基本性狀如下:供試土壤為黃崗土,pH 6.33、有機質23.56 g/kg、全氮1.89 g/kg、有效磷138 mg/kg、速效鉀464 mg/kg。

    1.2 試驗設計

    1.2.1 番茄育苗 將供試番茄種子用2% 次氯酸鈉表面消毒3 min,再用無菌蒸餾水洗滌3次后在合適的溫度和水分環(huán)境下發(fā)芽。種子發(fā)芽后,播入育苗基質中于育苗車間進行培育。當幼苗長出3 ~ 4片真葉時,選取長勢一致的幼苗進行田間移植。

    1.2.2 田間試驗 供試設施大棚自2013年2月開始于每年春夏及秋冬兩季連續(xù)種植番茄。本試驗于2015年2—6月在此設施大棚開展,此前已連續(xù)種植4季番茄,部分蕃茄植株已呈現(xiàn)出感染青枯病的癥狀。試驗設4個小區(qū),每個小區(qū)面積為10 m2(長×寬=2 m×5 m)。種植前施用7 kg商品雞糞有機肥(N + P2O5+ K2O> 5%)及0.5 kg復合肥(N-P2O5-K2O,15- 15-15)作為基肥,視長勢追施尿素及氯化鉀。所用商品雞糞有機肥、復合肥、尿素及氯化鉀均為市售常見肥料?;适┤? d后,每個小區(qū)移栽32棵番茄種苗,種植密度為0.6 m × 0.6 m(株距×行距)。移栽定植番茄種苗20 d后,每周監(jiān)測番茄青枯病發(fā)生狀況。為盡量避免土壤過于濕潤所導致的潛在的細菌過度交叉污染,除定植移苗后一次性澆透以外,日常水分管理采取“少量多次”的策略,利用滴灌設施向番茄莖基部補充水分。

    1.2.3 番茄土壤樣品采集及DNA提取 于番茄收獲期每小區(qū)分別隨機選定5株青枯病感病或健康(未明顯感病)植株,用鐵鏟將距番茄莖部10 cm以內且0 ~ 20 cm深的土壤和根系一起挖取,將5株樣品混勻后帶回實驗室進行后續(xù)處理。輕抖根系,將附著根上的土壤抖落,連同土體土壤利用四分法保留約200 g土壤作為感病(DB)及健康(HB)土體土壤樣品。將去除輕松附著根系上的番茄根置入預先裝有150 ml無菌水的500 ml三角燒瓶中,室溫下以170 r/min震蕩30 min,震蕩后4 000 ×離心5 min,離心后沉淀即為根際土壤,將收集到的根際土壤混勻保存于–70 °C備用。稱取0.4 g土壤樣品,用土壤DNA提取試劑盒(Mo Bio Laboratories, Inc., Carlsbad, CA, USA)按操作說明提取土壤DNA并用核酸定量儀(NanoDrop 2000, USA)測定提取的土壤DNA核酸濃度。

    1.3 測定方法

    1.3.1 土壤理化性質的測定 土壤理化性質測定參照《土壤農化分析》[10]。利用玻璃電極酸度計測定土壤pH((水)︰(土)=2.5︰1);利用電導率儀測定土壤電導率(EC)((水)︰(土) =5︰1);利用元素分析儀(Vario EL, Germany)通過干燒法測定土壤的全碳(TOC)和全氮(TON)含量;0.01 mol/L氯化鈣溶液浸提,流動分析儀(Auto Analyzer 3, Germany)測定土壤懸液的NH4+-N和NO– 3-N含量;醋酸銨浸提,火焰分光光度計測定土壤有效鉀(AK)含量;碳酸氫鈉浸提,鉬銻抗比色法測定土壤有效磷(AP)含量。

    1.3.2 土壤細菌和青枯菌數量的測定 土壤細菌和青枯菌的數量采用實時熒光定量PCR法測定[12-13]。細菌擴增采用引物Eub338/Eub518 (5′- ACTCCTA-CGGGAGGCAGCAG- 3′/5′- ATTACCGCG GCTGC-TGG - 3′),青枯菌擴增采用F/R (5′- GAACG-CCAACGGTGCGAACT - 3′/5′- GGCGGCCTTCA-GGGAGGTC-3′)。將含大腸桿菌16 S全長及青枯菌基因的質粒進行梯度稀釋,采用ABI 7500熒光定量PCR儀,按照標準程序制成標準曲線。擴增條件為:10 μl SYBR?Premix Ex Taq,0.4 μl上游引物和下游引物,0.4 μl ROX Reference Dye,2 μl 模板DNA和 6.8 μl無菌水。PCR擴增程序為:95 °C預變性30 s;95 °C預變性5 s,60 °C延伸34 s,循環(huán)40次。所有樣品進行3次重復,用無菌超純水代替DNA模板設置為陰性對照。根據各樣品值計算每克土壤所含的拷貝數,并換算成每克干土形成的拷貝數,取對數值,以log (copies/g 干土) 表示。

    1.3.3 細菌測序文庫構建 參照Caporaso等[14]發(fā)表的方法構建細菌測序文庫。細菌16S rRNA擴增采用引物520F (5′- AYTGGGYDTAAAGNG - 3′)/802R (5′- TACNVGGGTATCTAATCC- 3′)。細菌16S rRNA擴增體系為(25 μl):5.0 μl 5 × Reaction Buffer、5.0 μl 5 × GC high Enhancer、2.0 μl dNTP (2.5 mmol/L)、1.0 μl上下游引物(10 μmol/L)、0.25 μl TaKaRa Polymerase (5 U/μl)、2 μl模板DNA及9.0 μl 無菌水。PCR擴增條件為:98 °C預變性5 min;98 °C變性30 s,50 °C退火30 s,72 °C延伸30 s,25個循環(huán);最終72 °C延伸5 min。對PCR產物進行割膠回收后利用BioTek酶標儀對純化產物進行定量。將純化的擴增子等摩爾濃度合并,然后用NEB Next? UltraTM DNA Library Prep Kit for Illumina(New England Biolabs, UK)建庫,使用Agilent 2100 Bioanalyzer Instruments、Pico green和熒光分光光度計對文庫進行質量和濃度的檢測。多樣品測序文庫均一化至10 nmol/L后等體積混合,逐步稀釋定量至4 ~ 5 pmol/L后進行上機測序。所有PCR擴增、文庫準備及上機測序均在上海派森諾生物科技股份有限公司完成。

    1.3.4 高通量數據分析 利用QIIME軟件對下機后的原始序列根據相應的barcode信息進行各樣品序列分配并去除接頭和引物序列;去除低質量序列后,利用Usearch軟件對正反向序列進行拼接;拼接后的序列參照UPARSE的標準流程分析,生成操作分類單元(OTU)表格并挑選出各OTU代表性序列[15];利用RDP classifier將代表性序列與RDP 16S rRNA數據庫比對獲取各OTU分類組成信息[16];最后,對去除所有樣品中相對豐度均小于0.01% 的OTU的數據進行后續(xù)多樣性分析。在Mothur軟件[17]內基于Bray-Curtis參數計算各樣品主成分分析(principal co-ordinates analysis,PCoA)比較各處理細菌群落結構差異;各樣品取平后,基于97% 相似度水平利用Mothur計算細菌群落豐富度Chao值及多樣性Shannon值;利用韋恩圖探究樣品之間組成上的差異性;各物種的相對豐度用該物種序列數除以該樣品總序列數的百分比表示。

    1.4 數據統(tǒng)計分析

    生物數據統(tǒng)計分析在IBM SPSS 20.0(SPSS Inc., Chicago, USA)及R語言中完成。采用多因素方差分析、多重比較及獨立樣本檢驗分析比較處理間的顯著性差異。

    2 結果與分析

    2.1 健康與感病番茄土壤基本理化性狀

    即使在同一設施大棚小區(qū)內,健康與感病番茄植株土體土壤的部分基本理化性質差異明顯(表1)。健康與感病番茄土壤中電導率、NH4+-N、NO– 3-N、有效磷和全氮含量無顯著差異,但健康植株土壤中pH和全碳含量顯著高于感病植株土壤,速效鉀的含量顯著低于感病植株土壤。

    表1 健康與感病番茄土體土壤基本理化性質

    注:表中同列數據小寫字母不同表示健康和感病番茄土體土壤間理化性質差異達0.05顯著水平。

    2.2 健康與感病番茄土壤細菌及青枯菌數量

    熒光定量 PCR 結果表明,根際土壤中細菌和青枯菌的數量顯著高于土體土壤(0.05),根際青枯菌的數量比土體青枯菌的數量高出 3 個數量級。健康番茄植株土體及根際土壤中細菌的數量分別顯著高于感病番茄植株土體及根際土壤,而健康番茄植株土體及根際土壤青枯菌數量顯著低于感病番茄植株土體及根際土壤(圖 1)。

    2.3 健康與感病番茄土壤細菌群落測序數據總體概況

    高通量原始序列經過基本質控過濾并去除古菌及低豐度序列后,16 個土壤樣品共有 856 374 條優(yōu)質16S rRNA序列,其中最少的樣品有 38 449 條序列,而最多的樣品有 74 009 條序列。在 97% 相似度水平下,856 374 條序列共劃分為 3 452 個 OTU,其中最少的樣品有 1 168 個OTU,而最多的樣品有 3 079 個 OTU。與 RDP 數據庫比對后,856 374 條序列可歸類于 38 個細菌門,其中酸桿菌門、放線菌門、擬桿菌門、厚壁菌門及變形菌門為相對豐度最高的前 5 個門,占序列總數的 86.7%。

    2.4 健康與感病番茄土壤細菌群落alpha多樣性

    如圖2 所示,根際土壤中細菌群落的 alpha 多樣性指標,如 Chao 及 Shannon 指數均顯著低于土體土壤(檢驗,<0.05)。盡管健康番茄土體土壤豐富度 Chao 值顯著高于感病土壤,但健康及感病番茄根際土壤中 Chao 值無顯著差異。與之相反,盡管健康與感病番茄土體土壤多樣性 Shannon 指數無顯著差異,但健康番茄根際土壤中 Shanon 指數顯著高于感病番茄根際土壤。

    2.5 健康與感病番茄土壤細菌群落beta多樣性

    健康及感病番茄根際土壤細菌群落主坐標成分在一軸上明顯區(qū)分于其土體土壤(圖 3)。健康番茄根際土壤細菌群落主成分在二軸上明顯區(qū)分于感病番茄根際土壤,盡管健康與感病番茄土體土壤細菌主成分差異不如根際土壤差異明顯,但健康番茄土體土壤細菌群落主成分仍明顯不同于感病番茄土體土壤。進一步的 PerMANOVA 分析,也表明土壤類型(土體與根際)及土壤健康(健康與感病)均是細菌群落差異的顯著性驅動因素(Adonis 檢驗,F(xiàn)土體vs根際= 40.8,土體vs根際= 0.001;健康vs感病= 6.7,健康vs感病= 0.009)。

    2.6 健康與感病番茄土壤細菌群落組成差異

    如圖 4 所示,與土體土壤相比,番茄根際土壤在細菌門水平上的組成類似,但門的相對豐度有所差異,如變形菌門的相對豐度顯著增高,而酸桿菌門、放線菌門、芽單胞菌門及疣微菌門的相對豐度顯著降低(檢驗,<0.05)。盡管健康與感病番茄土體土壤細菌主要門的相對豐度無顯著差異,但健康番茄根際土壤中擬桿菌門及疣微菌門的相對豐度顯著高于感病番茄根際土壤,而厚壁菌門、芽單胞菌門及變形菌門的相對豐度顯著低于感病番茄根際土壤(0.05)。

    在屬水平上健康與感病番茄土體及根際土壤的組成較為類似。在土體土壤中共鑒定出 691 個屬,其中有 657 個屬(95.1%)共同存在于健康及感病番茄土體土壤中;在根際土壤中共鑒定出 615 個屬,其中有 476 個屬(77.5%)共同存在于健康及感病番茄根際土壤中。健康及感病番茄土壤中共有屬相對豐度的檢驗結果表明:在健康與感病番茄土體土壤共有的 657 個屬中,健康番茄土體土壤中有 83 個屬的相對豐度顯著高于感病番茄土體土壤,而有 60 個屬的相對豐度顯著低于感病番茄土體土壤;在健康與感病番茄根際土壤共有的 476 個屬中,健康番茄根際土壤中有 32 個屬的相對豐度顯著高于感病番茄根際土壤,而有 60 個屬的相對豐度顯著低于感病番茄根際土壤(圖 5)。

    健康番茄土體土壤中相對豐度顯著增加的 83 個屬及顯著降低的 60 個屬均歸類于酸桿菌門、放線菌門、擬桿菌門、厚壁菌門、變形菌門、疣微菌門及部分低豐度門;而健康番茄根際土壤中相對豐度顯著降低的 60 個屬歸類于酸桿菌門、放線菌門、芽單胞菌門、厚壁菌門、變形菌門、疣微菌門及部分低豐度門,相對豐度顯著增加的 32 個屬歸類于擬桿菌門、厚壁菌門、變形菌門、疣微菌門及部分低豐度門(圖 6)。其中,健康番茄根際顯著增加和降低的屬分別多屬于擬桿菌門和變形菌門。

    如圖 7 所示,在健康番茄根際土壤中相對豐度顯著降低的變形菌門中,青枯病致病菌所在的青枯菌屬的相對豐度僅為 27.4%,顯著低于感病番茄根際土壤中青枯菌屬的相對豐度(54.6%),而土體土壤之間無顯著差異。在健康番茄根際土壤中相對豐度顯著增加的擬桿菌門中,噬幾丁質菌屬、金桿菌屬、動桿菌屬、黃桿菌屬及的相對豐度大于 1%,為優(yōu)勢菌屬,而土體土壤之間也無顯著差異。

    3 討論

    根際土壤中細菌群落的多樣性及穩(wěn)定性決定著土壤和植物的健康以及生態(tài)系統(tǒng)的可持續(xù)性[18]。抑病型土壤指的是病原菌不能定殖,或能定殖但危害很小或沒有危害,或者能定殖并一時造成危害隨后即使在病原菌存在的情況下發(fā)病也很輕的土壤;而容易感病的土壤即為導病型土壤[19]。土著微生物群落或者特定種群具有抑制病原菌侵染植物的能力是抑病型土壤的關鍵因素[20-21]。本研究通過比較連作番茄田塊中健康番茄土壤(抑病型土壤)及感病番茄土壤(導病型土壤)細菌組成差異,利用高通量測序技術深入揭示了抑制土傳青枯病發(fā)生的土壤細菌群落結構及組成特征。

    與Li等[5]的研究結果類似,本研究中健康番茄土體及根際土壤細菌群落明顯不同于感病番茄土體及根際土壤,表明抑制土傳青枯病發(fā)生的土壤細菌群落確與感病土壤細菌群落不同。Garbeva等[22]研究表明,根際土壤細菌多樣性高能抑制土傳病害的發(fā)生,與之類似,本研究中健康番茄根際土壤細菌群落的Shannon多樣性指數更高,證實了細菌多樣性與生態(tài)系統(tǒng)生產力之間存在的正相關關系。

    Li 等[5]利用454 焦磷酸測序技術研究番茄根際土壤細菌群落組成時發(fā)現(xiàn)番茄根際土壤細菌群落中變形菌門的相對豐度最高。本研究中根際細菌群落中變形菌門的相對豐度也占優(yōu)勢地位,可能是因為其相對較快的生長速率[23]。本研究結果表明,抑制番茄土傳青枯病發(fā)生的根際土壤中擬桿菌門的相對豐度顯著高于感病土壤。擬桿菌門通常棲居于土體土壤中,而且是微生態(tài)系統(tǒng)中相對穩(wěn)定的組分[24]。盡管擬桿菌門的相對豐度很少受土壤養(yǎng)分狀況的影響,但近年來越來越多的研究表明擬桿菌門可能與土壤抑病能力相關[25]。本研究中擬桿菌門所包含的噬幾丁質菌屬、金桿菌屬、動桿菌屬、黃桿菌屬及均為抑制土傳青枯病發(fā)生的根際土壤中的優(yōu)勢菌屬(相對豐度>1%),且與抑病能力密切相關。Liu 等[26]研究表明噬幾丁質菌屬與番茄青枯病之間呈負相關關系,而Inderbitzin 等[27]研究發(fā)現(xiàn)西蘭花和幾丁質混合物中接種噬幾丁質菌屬可提升土壤抑制輪枝菌侵染誘發(fā)的黃萎病的能力。金桿菌屬已被成功分離并應用于辣椒疫病的生物防治[28]。Kwak 等[25]研究表明,黃桿菌屬在番茄根際土壤抵御青枯病侵染的過程中發(fā)揮重要作用。盡管尚未有直接證據表明其具有生防功能,但近年來的研究表明動桿菌屬及在健康土壤中的相對豐度較高,可能也與抑病相關[29]。此外,本研究中健康番茄根際土壤中無論是青枯菌屬的相對豐度還是絕對數量雖顯著低于感病根際土壤,但仍屬于較高水平,說明健康根際的整體細菌群落對于提升土壤抑病能力也具有一定的作用。

    土壤微生物的裝配受植物基因型和土壤狀況等共同影響[30]。盡管感病與健康樣品在同一設施大棚內采集,但可能因土壤的異質性、施肥的不均勻、灌溉不均勻以及與作物互作關系不一致等原因,導致健康與感病土壤樣品的部分理化性狀有顯著差異。本研究中,抑制土傳青枯病發(fā)生的土壤擁有較高的pH及全碳含量。Deltour等[31]研究發(fā)現(xiàn)高土壤pH能增強土壤抑鐮刀菌枯萎病能力。土壤pH不僅能通過作用于病原菌或其他微生物直接影響土壤抑病型能力,還能通過影響土壤速效養(yǎng)分間接影響土壤抑病型能力[32]。此外,土壤全碳含量高通常預示著土壤有機質含量較高,而有機質含量較高的土壤中通常微生物生物量及活力也較高,也能在抑病過程中發(fā)揮重要作用[22]。Satoh和Toyota[33]研究證實提升土壤中有機質含量可降低土壤中青枯菌的數量,進而提升土壤抑制青枯病的能力??偠灾嵘寥纏H及全碳(有機質)含量,不僅可直接降低土壤中青枯菌,還可通過調控土體土壤微生物組成,進而影響根際微生物群落裝配,并最終決定群落整體的抑病能力。

    微生物、土壤及植物之間的互作關系非常復雜,影響了植物地下部微生物組成及生態(tài)系統(tǒng)功能,進而決定了地上部植物的生長與健康。本研究借助MiSeq高通量測序技術利用全面的分析比較了同一連作田塊中健康與患病番茄土體及根際土壤細菌群落的組成特征,結果表明,相較于感病土壤,健康土壤中pH及有機質含量較高,其棲居的細菌群落組成特異,所種植的番茄從土壤中招募微生物裝配形成的根際細菌群落多樣性高,且富集了噬幾丁質菌屬、金桿菌屬、動桿菌屬、黃桿菌屬及等土著有益微生物,不僅有效降低了根際土壤中青枯病的相對豐度,還提升了群落整體抑病能力。但番茄如何招募這些有益微生物在根際進行裝配將是未來研究的熱點與難點。

    [1] 賀紀正, 李晶, 鄭袁明.土壤生態(tài)系統(tǒng)微生物多樣性-穩(wěn)定性關系的思考[J].生物多樣性, 2013, 21(4): 411–420

    [2] Berendsen R L, Pieterse C M J, Bakker P A H M. The rhizosphere microbiome and plant health[J]. Trends in Plant Science, 2012, 17(8): 478–486.

    [3] Hassani M A, Durán P, Hacquard S. Microbial interactions within the plant holobiont[J]. Microbiome, 2018, 6: 58.

    [4] Khan S, Hanjra M A, Mu J X. Water management and crop production for food security in China: a review[J]. Agricultural Water Management, 2009, 96(3): 349–360.

    [5] Li J G, Ren G D, Jia Z J, et al. Composition and activity of rhizosphere microbial communities associated with healthy and diseased greenhouse tomatoes[J]. Plant and Soil, 2014, 380(1/2): 337–347.

    [6] 韋中. 生物有機肥防控土傳番茄青枯病的效果及其機制研究[D]. 南京: 南京農業(yè)大學, 2012.

    [7] Raaijmakers J M, Paulitz T C, Steinberg C, et al. The rhizosphere: A playground and battlefield for soilborne pathogens and beneficial microorganisms[J]. Plant and Soil, 2009, 321(1/2): 341–361.

    [8] 張鵬, 韋中, 朱震, 等. 生物有機肥對連作番茄和辣椒根際土壤微生物區(qū)系及茄科雷爾氏菌的影響[J]. 南京農業(yè)大學學報, 2013, 36(4): 77–82.

    [9] Wei Z, Yang X M, Yin S X, et al. Efficacy of Bacillus-fortified organic fertiliser in controlling bacterial wilt of tomato in the field[J]. Applied Soil Ecology, 2011, 48(2): 152–159.

    [10] 鮑士旦. 土壤農化分析(第3版)[M]. 北京: 中國農業(yè)出版社, 2000.

    [11] 楊尚東, 趙久成, 郭伊娟, 等. 番茄青枯病罹病植株和健康植株根際土壤細菌群落結構的初步分析[J]. 中國蔬菜, 2014(8): 25–29.

    [12] Fierer N, Jackson J A, Vilgalys R, et al. Assessment of soil microbial community structure by use of taxon-specific quantitative PCR assays[J]. Applied and Environmental Microbiology, 2005, 71(7): 4117–4120.

    [13] Sch?nfeld J, Heuer H, van Elsas J D, et al. Specific and sensitive detection of Ralstonia solanacearum in soil on the basis of PCR amplification of fliC fragments[J]. Applied and Environmental Microbiology, 2003, 69(12): 7248–7256.

    [14] Caporaso J G, Lauber C L, Walters W A, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample[J]. PNAS, 2011, 108(Supplement_1): 4516–4522.

    [15] Edgar R C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads[J]. Nature Methods, 2013, 10(10): 996–998.

    [16] Wang Q, Garrity G M, Tiedje J M, et al. Na?ve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy[J]. Applied and Environmental Microbiology, 2007, 73(16): 5261–5267.

    [17] Schloss P D, Westcott S L, Ryabin T, et al. Introducing mothur: Open-source, platform-independent, community- supported software for describing and comparing microbial communities[J]. Applied and Environmental Microbiology, 2009, 75(23): 7537–7541.

    [18] Trivedi P, He Z L, van Nostrand J D, et al. Huanglongbing alters the structure and functional diversity of microbial communities associated with citrus rhizosphere[J]. The ISME Journal, 2012, 6(2): 363–383.

    [19] 張瑞福, 沈其榮.抑病型土壤的微生物區(qū)系特征及調控[J]. 南京農業(yè)大學學報, 2012, 35(5): 125–132.

    [20] Weller D M, Raaijmakers J M, Gardener B B M, et al. Microbial populations responsible for specific soil suppressiveness to plant pathogens[J]. Annual Review of Phytopathology, 2002, 40(1): 309–348.

    [21] Mendes R, Kruijt M, de Bruijn I, et al. Deciphering the rhizosphere microbiome for disease-suppressive bacteria[J]. Science, 2011, 332(6033): 1097–1100.

    [22] Garbeva P, van Veen J A, van Elsas J D. Microbial diversity in soil: Selection of microbial populations by plant and soil type and implications for disease suppressiveness[J]. Annual Review of Phytopathology, 2004, 42(1): 243–270.

    [23] Fierer N, Bradford M A, Jackson R B. Toward an ecological classification of soil bacteria[J]. Ecology, 2007, 88(6): 1354–1364.

    [24] Yergeau E, Schoondermark-Stolk S A, Brodie E L, et al. Environmental microarray analyses of Antarctic soil microbial communities[J]. The ISME Journal, 2009, 3(3): 340–351.

    [25] Kwak M J, Kong H G, Choi K, et al. Rhizosphere microbiome structure alters to enable wilt resistance in tomato[J]. Nature Biotechnology, 2018, 36(11): 1100– 1109.

    [26] Liu H J, Xiong W, Zhang R F, et al. Continuous application of different organic additives can suppress tomato disease by inducing the healthy rhizospheric microbiota through alterations to the bulk soil microflora[J]. Plant and Soil, 2018, 423(1/2): 229–240.

    [27] Inderbitzin P, Ward J, Barbella A, et al. Soil microbiomes associated with verticillium wilt-suppressive broccoli and chitin amendments are enriched with potential biocontrol agents[J]. Phytopathology, 2018, 108(1): 31–43.

    [28] Kim H S, Sang M K, Jung H W, et al. Identification and characterization of Chryseobacterium wanjuense strain KJ9C8 as a biocontrol agent of Phytophthora blight of pepper[J]. Crop Protection, 2012, 32: 129–137.

    [29] Fu L, Penton C R, Ruan Y Z, et al. Inducing the rhizosphere microbiome by biofertilizer application to suppress banana Fusarium wilt disease[J]. Soil Biology and Biochemistry, 2017, 104: 39–48.

    [30] Lundberg D S, Lebeis S L, Paredes S H, et al. Defining the core Arabidopsis thaliana root microbiome[J]. Nature, 2012, 488(7409): 86–90.

    [31] Deltour P, C Fran?a S, Liparini Pereira O, et al. Disease suppressiveness to Fusarium wilt of banana in an agroforestry system: Influence of soil characteristics and plant community[J]. Agriculture, Ecosystems & Environ-ment, 2017, 239: 173–181.

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

    [33] Satoh K, Toyota K. Comparison of disease suppressiveness of different soils with or without repeated application of organic matters toward bacterial wilt of tomato caused by ralstonia solanacearum[J]. Microbes and Environments, 2004, 19(4): 310–314.

    Comparison of Bacterial Communities in Bulk and Rhizosphere Soils of Healthy and Diseased Tomato Infected by Bacterial Wilt

    SHEN Zongzhuan1, HUANG Yan1, CAO Yifan1, WANG Dongsheng1,2, LIU Hongjun1, LI Rong1*, SHEN Qirong1

    (1 College of Resources and Environmental Sciences, Jiangsu Provincial Key Laboratory for Organic Solid Waste Utilization, Jiangsu Collaborative Innovation Center for Solid Organic Waste Resource Utilization, National Engineering Research Center for Organic-based Fertilizers, Key Laboratory of Plant Immunity, Nanjing Agricultural University, Nanjing 210095, China;2 Nanjing Institute of Vegetable Science, Nanjing 210042, China)

    In this study, real time PCR and MiSeq sequencing technology was used to deeply compare the differences of bacterial community structures and compositions in bulk and rhizosphere soils of healthy and diseased tomato-planting fields under mono-cropping system. The results showed that healthy bulk soil displayed higher pH value and total carbon content. The bacterial structures and compositions in healthy bulk and rhizosphere soils were obviously different from those of diseased ones. Compared to diseased rhizosphere soil, healthy rhizosphere soil exhibited a significant higher number of total bacteria and a significant lower number of. Also, healthy rhizosphere soil showed a significant higher index of Shannon diversity. Furthermore, the relative abundance of Bacteroides and,,,andincluded in Bacteroides was significantly enriched while the relative abundance of Proteobacteria andincluded in Proteobacteria was significantly depleted in healthy rhizosphere soil compared to diseased one. In summary, soil suppressive to tomato bacterial wilt disease harbored a unique bacterial community, which showed high bacterial population, diversity and relative abundance of beneficial indigenous microorganisms and low relative abundance of pathogen. This research could provide a guideline and theoretical principle for biological control of tomato bacterial wilt disease.

    Tomato wilt disease; Disease-suppressive soil; High-throughput sequencing; Bacterial community diversity; Bacterial community composition

    S154

    A

    10.13758/j.cnki.tr.2021.01.002

    沈宗專, 黃炎, 操一凡, 等. 健康與罹患青枯病的番茄土壤細菌群落特征比較. 土壤, 2021, 53(1): 5–12.

    國家自然科學基金項目(41977044)、江蘇省自然科學基金項目(BK20160710)和江蘇省政策引導類計劃(國際科技合作/港澳臺科技合作)項目(BZ2020026)資助。

    (lirong@njau.edu.cn)

    沈宗專(1987—),男,江蘇寶應人,博士,副教授,主要從事土壤微生物與生物肥料研究。E-mail: victorshenzz@163.com

    猜你喜歡
    土傳感病青枯病
    土壤pH和主要養(yǎng)分含量與山核桃干腐病的相關性研究
    同種負密度制約效應與園林樹木土傳病害的有效防控
    北京園林(2020年4期)2020-01-18 05:16:50
    5個歐亞種葡萄品種感染霜霉病后4種酶活性的變化
    土壤真菌多樣性對土傳病害影響的研究進展
    察右后旗馬鈴薯土傳病害防治試驗總結
    土傳疾病的生物解決方案
    營銷界(2015年23期)2015-02-28 22:06:15
    淺談茄子青枯病防治技術
    葡萄白粉病與幾種酶活性的關系
    三種堆肥對番茄生長及青枯病防治效果的影響
    不同抗性甜瓜接種蔓枯病菌后若干生理指標的變化
    中國瓜菜(2012年1期)2012-04-29 00:44:03
    巨乳人妻的诱惑在线观看| 欧美色欧美亚洲另类二区| 亚洲欧美一区二区三区黑人| 日韩欧美免费精品| 叶爱在线成人免费视频播放| 偷拍熟女少妇极品色| 99国产综合亚洲精品| a级毛片a级免费在线| 国产三级黄色录像| 久久久国产精品麻豆| 亚洲欧美日韩卡通动漫| 久久天躁狠狠躁夜夜2o2o| 18禁美女被吸乳视频| 夜夜爽天天搞| 亚洲国产色片| 精品久久蜜臀av无| 99久久无色码亚洲精品果冻| 亚洲人成电影免费在线| 亚洲国产色片| av片东京热男人的天堂| 亚洲天堂国产精品一区在线| 国产人伦9x9x在线观看| 亚洲av成人精品一区久久| 亚洲国产精品合色在线| 黑人巨大精品欧美一区二区mp4| 日本成人三级电影网站| 亚洲在线自拍视频| 国产探花在线观看一区二区| 亚洲人成网站高清观看| 日本三级黄在线观看| 99热只有精品国产| 午夜福利欧美成人| 午夜a级毛片| 国产99白浆流出| 精品久久久久久久人妻蜜臀av| 国产精品美女特级片免费视频播放器 | 99国产精品99久久久久| 国产成人欧美在线观看| 久久久久久久精品吃奶| 一进一出好大好爽视频| 神马国产精品三级电影在线观看| 欧美乱码精品一区二区三区| 免费在线观看视频国产中文字幕亚洲| 一个人看的www免费观看视频| 在线十欧美十亚洲十日本专区| 国产成人欧美在线观看| 亚洲第一欧美日韩一区二区三区| 美女高潮的动态| 欧美中文日本在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 免费观看精品视频网站| 啦啦啦韩国在线观看视频| 国产亚洲精品一区二区www| 成年版毛片免费区| 又紧又爽又黄一区二区| 身体一侧抽搐| 伊人久久大香线蕉亚洲五| 国产精品,欧美在线| 麻豆av在线久日| 美女扒开内裤让男人捅视频| 国产乱人伦免费视频| 国产一区在线观看成人免费| 免费人成视频x8x8入口观看| 丰满人妻熟妇乱又伦精品不卡| 手机成人av网站| 在线观看一区二区三区| 午夜免费激情av| 久久热在线av| 三级男女做爰猛烈吃奶摸视频| 日韩 欧美 亚洲 中文字幕| 少妇裸体淫交视频免费看高清| 成年免费大片在线观看| 亚洲国产精品合色在线| 欧美黄色片欧美黄色片| 午夜激情福利司机影院| 国产亚洲精品综合一区在线观看| 日韩中文字幕欧美一区二区| 亚洲av美国av| 在线看三级毛片| 免费一级毛片在线播放高清视频| 男人舔奶头视频| 视频区欧美日本亚洲| 久99久视频精品免费| 精品久久久久久久人妻蜜臀av| 91字幕亚洲| 国产午夜精品久久久久久| 高清毛片免费观看视频网站| 免费在线观看成人毛片| 精品久久久久久久人妻蜜臀av| 亚洲电影在线观看av| 久久亚洲精品不卡| 亚洲18禁久久av| or卡值多少钱| 久久天躁狠狠躁夜夜2o2o| 久久久国产成人免费| 日韩大尺度精品在线看网址| 一本综合久久免费| 亚洲精品在线美女| 色在线成人网| 日本精品一区二区三区蜜桃| 日韩av在线大香蕉| 成人av在线播放网站| bbb黄色大片| 亚洲中文日韩欧美视频| av国产免费在线观看| 免费人成视频x8x8入口观看| 高潮久久久久久久久久久不卡| 日韩欧美 国产精品| 国产乱人视频| av女优亚洲男人天堂 | 夜夜夜夜夜久久久久| 精品久久蜜臀av无| 国产成人欧美在线观看| 夜夜躁狠狠躁天天躁| 麻豆一二三区av精品| 欧美极品一区二区三区四区| 波多野结衣巨乳人妻| 成在线人永久免费视频| 美女午夜性视频免费| 久久精品aⅴ一区二区三区四区| 我要搜黄色片| 无人区码免费观看不卡| 国产精品久久电影中文字幕| 亚洲成av人片在线播放无| 一二三四在线观看免费中文在| 搡老岳熟女国产| 香蕉av资源在线| 日本精品一区二区三区蜜桃| 成人特级av手机在线观看| 99久久无色码亚洲精品果冻| 亚洲人成伊人成综合网2020| 天堂影院成人在线观看| 国产精品国产高清国产av| 手机成人av网站| 国产三级在线视频| 国产精品,欧美在线| 18禁美女被吸乳视频| 俺也久久电影网| 日本黄色片子视频| 看片在线看免费视频| 精品国产乱子伦一区二区三区| 又爽又黄无遮挡网站| 成人18禁在线播放| 国产在线精品亚洲第一网站| 精品欧美国产一区二区三| 久9热在线精品视频| 我要搜黄色片| tocl精华| 国产精品永久免费网站| 国产毛片a区久久久久| 人人妻,人人澡人人爽秒播| 亚洲美女黄片视频| 淫妇啪啪啪对白视频| 真实男女啪啪啪动态图| 老司机在亚洲福利影院| 午夜a级毛片| 久久久久久国产a免费观看| 免费观看人在逋| 搡老妇女老女人老熟妇| 可以在线观看的亚洲视频| 亚洲欧洲精品一区二区精品久久久| 青草久久国产| 日韩三级视频一区二区三区| 制服丝袜大香蕉在线| 亚洲一区高清亚洲精品| 一本精品99久久精品77| 亚洲色图av天堂| 国产成人系列免费观看| 少妇裸体淫交视频免费看高清| 国产一区二区激情短视频| 欧美成人性av电影在线观看| 欧美中文综合在线视频| 日本熟妇午夜| 12—13女人毛片做爰片一| 亚洲专区字幕在线| 国产精品亚洲一级av第二区| 国产私拍福利视频在线观看| 999久久久精品免费观看国产| 国产精品免费一区二区三区在线| 久久精品91无色码中文字幕| 日韩欧美国产在线观看| 小说图片视频综合网站| 欧美在线黄色| 亚洲av熟女| 精品日产1卡2卡| 亚洲va日本ⅴa欧美va伊人久久| 国产又色又爽无遮挡免费看| 制服丝袜大香蕉在线| 欧美成人免费av一区二区三区| 欧美激情久久久久久爽电影| 国产一区二区在线观看日韩 | 亚洲电影在线观看av| www日本在线高清视频| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩亚洲国产一区二区在线观看| 久久精品国产综合久久久| 久久久久久人人人人人| 在线永久观看黄色视频| 日本 欧美在线| 欧美丝袜亚洲另类 | 变态另类成人亚洲欧美熟女| 午夜免费成人在线视频| 成人午夜高清在线视频| 在线播放国产精品三级| 麻豆久久精品国产亚洲av| 色综合站精品国产| x7x7x7水蜜桃| 人妻丰满熟妇av一区二区三区| 身体一侧抽搐| 国产精品日韩av在线免费观看| 亚洲成a人片在线一区二区| 欧美日韩国产亚洲二区| 最新中文字幕久久久久 | 亚洲在线自拍视频| 熟女电影av网| 99精品欧美一区二区三区四区| 免费看十八禁软件| xxx96com| 国产极品精品免费视频能看的| 黄片小视频在线播放| 亚洲精华国产精华精| 亚洲男人的天堂狠狠| 国产毛片a区久久久久| 欧美在线黄色| 国模一区二区三区四区视频 | 国产精品久久久人人做人人爽| 欧美成人一区二区免费高清观看 | 真实男女啪啪啪动态图| 一级毛片女人18水好多| 亚洲国产精品合色在线| 婷婷六月久久综合丁香| 91久久精品国产一区二区成人 | 国产成人福利小说| 午夜福利在线观看吧| 国产 一区 欧美 日韩| 成年女人看的毛片在线观看| 日本精品一区二区三区蜜桃| 香蕉国产在线看| 国产亚洲欧美98| 啦啦啦韩国在线观看视频| 亚洲专区中文字幕在线| 视频区欧美日本亚洲| 日韩欧美免费精品| 狂野欧美激情性xxxx| 国产精品一区二区三区四区免费观看 | 婷婷六月久久综合丁香| 午夜福利18| 成熟少妇高潮喷水视频| 欧美xxxx黑人xx丫x性爽| 我要搜黄色片| 淫秽高清视频在线观看| 久久久久久久久久黄片| av欧美777| 最好的美女福利视频网| 久久久久久九九精品二区国产| 国产视频一区二区在线看| 欧美精品啪啪一区二区三区| 亚洲av片天天在线观看| 啦啦啦韩国在线观看视频| 真实男女啪啪啪动态图| 久久伊人香网站| 99久国产av精品| 窝窝影院91人妻| 精品电影一区二区在线| 国产精品精品国产色婷婷| 国产成人影院久久av| 久久午夜综合久久蜜桃| 在线a可以看的网站| 天堂动漫精品| 亚洲欧美日韩高清专用| 国产熟女xx| 国产爱豆传媒在线观看| 日本免费a在线| 桃色一区二区三区在线观看| 成人亚洲精品av一区二区| 久久久久国产精品人妻aⅴ院| 午夜福利18| 亚洲九九香蕉| 久久精品综合一区二区三区| 日韩欧美三级三区| 久久天堂一区二区三区四区| 黄色女人牲交| 99久久精品热视频| 久久精品国产清高在天天线| 免费av毛片视频| 欧美午夜高清在线| 色综合亚洲欧美另类图片| 日本 av在线| 亚洲片人在线观看| 国产高清激情床上av| 久久精品aⅴ一区二区三区四区| 国产一区二区在线观看日韩 | 日本一本二区三区精品| 亚洲 国产 在线| 人人妻人人澡欧美一区二区| 亚洲美女黄片视频| 亚洲国产精品sss在线观看| 亚洲成人中文字幕在线播放| 啦啦啦免费观看视频1| 最新中文字幕久久久久 | 九色成人免费人妻av| 19禁男女啪啪无遮挡网站| 久久精品综合一区二区三区| 一个人免费在线观看的高清视频| or卡值多少钱| 男女之事视频高清在线观看| 免费人成视频x8x8入口观看| 久久天躁狠狠躁夜夜2o2o| 欧美在线一区亚洲| 91久久精品国产一区二区成人 | 亚洲专区国产一区二区| 国产精品1区2区在线观看.| 女同久久另类99精品国产91| 日韩中文字幕欧美一区二区| 精品国产亚洲在线| 一区二区三区高清视频在线| 国产一区二区三区视频了| 国内少妇人妻偷人精品xxx网站 | 国产精品 国内视频| 黄片大片在线免费观看| 精品欧美国产一区二区三| 最近最新免费中文字幕在线| 丝袜人妻中文字幕| 亚洲午夜理论影院| 国产av麻豆久久久久久久| 亚洲av中文字字幕乱码综合| 午夜日韩欧美国产| 国产午夜福利久久久久久| 夜夜看夜夜爽夜夜摸| 51午夜福利影视在线观看| 国产亚洲精品综合一区在线观看| 国产午夜精品久久久久久| 中文在线观看免费www的网站| 欧美在线黄色| 欧美日韩亚洲国产一区二区在线观看| 1024香蕉在线观看| 色综合欧美亚洲国产小说| 免费看a级黄色片| 黄色片一级片一级黄色片| 免费在线观看影片大全网站| www.www免费av| 色老头精品视频在线观看| 观看美女的网站| 日韩欧美在线乱码| 成在线人永久免费视频| 国产午夜精品论理片| 综合色av麻豆| 国产亚洲精品综合一区在线观看| 天天添夜夜摸| 真人做人爱边吃奶动态| 哪里可以看免费的av片| 欧美最黄视频在线播放免费| 长腿黑丝高跟| 九色国产91popny在线| 亚洲激情在线av| 欧美中文综合在线视频| 久久久成人免费电影| 国产精品 欧美亚洲| 日韩 欧美 亚洲 中文字幕| 老司机午夜福利在线观看视频| 国产男靠女视频免费网站| 三级毛片av免费| 国产精品综合久久久久久久免费| 国产真实乱freesex| 亚洲成人久久性| 又大又爽又粗| 免费av不卡在线播放| 嫩草影院入口| 九九在线视频观看精品| 桃红色精品国产亚洲av| 国产一级毛片七仙女欲春2| 日本一二三区视频观看| 999久久久精品免费观看国产| 禁无遮挡网站| 我的老师免费观看完整版| 久久性视频一级片| 麻豆国产av国片精品| 1000部很黄的大片| 母亲3免费完整高清在线观看| 日韩精品中文字幕看吧| 欧美黑人巨大hd| 国产精品一区二区精品视频观看| 国模一区二区三区四区视频 | 99国产精品一区二区蜜桃av| 国产精品99久久99久久久不卡| 一个人观看的视频www高清免费观看 | 久久精品人妻少妇| 99热6这里只有精品| 伊人久久大香线蕉亚洲五| 无人区码免费观看不卡| 日韩欧美国产一区二区入口| 欧美黄色淫秽网站| 91av网一区二区| 色av中文字幕| 精品国产美女av久久久久小说| 高清在线国产一区| 亚洲av美国av| 免费一级毛片在线播放高清视频| 亚洲性夜色夜夜综合| 久久久国产成人免费| 九九在线视频观看精品| 99在线视频只有这里精品首页| 午夜视频精品福利| 欧美色欧美亚洲另类二区| 色在线成人网| 99热这里只有精品一区 | 国产97色在线日韩免费| 亚洲国产看品久久| 欧洲精品卡2卡3卡4卡5卡区| 免费av毛片视频| 久久热在线av| 成年女人永久免费观看视频| 精品久久久久久久人妻蜜臀av| 综合色av麻豆| 亚洲一区二区三区不卡视频| 久久久久久久精品吃奶| 熟女电影av网| 精品一区二区三区视频在线 | 国产亚洲欧美在线一区二区| 香蕉久久夜色| 亚洲精品中文字幕一二三四区| 曰老女人黄片| 亚洲五月婷婷丁香| 久久久久久九九精品二区国产| 久久这里只有精品19| 又粗又爽又猛毛片免费看| 两人在一起打扑克的视频| 成人三级黄色视频| 国产精品永久免费网站| 美女高潮喷水抽搐中文字幕| 精品久久久久久久人妻蜜臀av| av在线天堂中文字幕| 亚洲最大成人中文| 丝袜人妻中文字幕| 日本免费一区二区三区高清不卡| 97超级碰碰碰精品色视频在线观看| 久久国产乱子伦精品免费另类| 色视频www国产| 啦啦啦免费观看视频1| svipshipincom国产片| 国产精品一及| 脱女人内裤的视频| 极品教师在线免费播放| 国产精品国产高清国产av| 亚洲精华国产精华精| 午夜两性在线视频| 亚洲国产色片| 伊人久久大香线蕉亚洲五| 亚洲欧美激情综合另类| avwww免费| 亚洲一区高清亚洲精品| 婷婷精品国产亚洲av在线| aaaaa片日本免费| 免费在线观看影片大全网站| 国产亚洲av嫩草精品影院| 国产午夜福利久久久久久| 国产淫片久久久久久久久 | 国产激情欧美一区二区| 俄罗斯特黄特色一大片| 白带黄色成豆腐渣| 国产亚洲av嫩草精品影院| 国产高潮美女av| 亚洲成人久久爱视频| 女人高潮潮喷娇喘18禁视频| 蜜桃久久精品国产亚洲av| 久久久久久久久久黄片| 国产1区2区3区精品| 黄色视频,在线免费观看| 日韩国内少妇激情av| 夜夜爽天天搞| 悠悠久久av| 日本黄色片子视频| 色综合亚洲欧美另类图片| 看免费av毛片| 日本免费a在线| 日本成人三级电影网站| bbb黄色大片| 日韩欧美精品v在线| 欧美中文综合在线视频| 国产av不卡久久| 欧美日韩国产亚洲二区| 最近最新中文字幕大全免费视频| 亚洲熟女毛片儿| 黄色片一级片一级黄色片| 亚洲成人中文字幕在线播放| av黄色大香蕉| 欧美一区二区国产精品久久精品| 9191精品国产免费久久| a级毛片在线看网站| av在线蜜桃| 麻豆成人午夜福利视频| 亚洲成人精品中文字幕电影| 一本精品99久久精品77| 久久久国产精品麻豆| 不卡av一区二区三区| 一本久久中文字幕| 国产精品久久久久久亚洲av鲁大| 亚洲精品乱码久久久v下载方式 | 精品一区二区三区视频在线 | а√天堂www在线а√下载| 欧美日韩综合久久久久久 | 在线观看66精品国产| 麻豆成人午夜福利视频| 精华霜和精华液先用哪个| 国产高清三级在线| 舔av片在线| 国产乱人伦免费视频| 别揉我奶头~嗯~啊~动态视频| 日本五十路高清| 国产精品 欧美亚洲| 天堂√8在线中文| 一个人看视频在线观看www免费 | 最新美女视频免费是黄的| 九九热线精品视视频播放| 手机成人av网站| 别揉我奶头~嗯~啊~动态视频| 成人三级黄色视频| 午夜久久久久精精品| 麻豆国产97在线/欧美| 欧美精品啪啪一区二区三区| 精品欧美国产一区二区三| 日韩欧美精品v在线| 久久国产精品人妻蜜桃| 偷拍熟女少妇极品色| 国产精品亚洲一级av第二区| 日本撒尿小便嘘嘘汇集6| 国产精品av视频在线免费观看| 99精品在免费线老司机午夜| 国产亚洲欧美在线一区二区| 91在线观看av| а√天堂www在线а√下载| 国产午夜福利久久久久久| 国产麻豆成人av免费视频| 亚洲成人精品中文字幕电影| 级片在线观看| 哪里可以看免费的av片| 欧美乱色亚洲激情| 午夜精品在线福利| 精品久久久久久,| 男女那种视频在线观看| 叶爱在线成人免费视频播放| 亚洲 欧美 日韩 在线 免费| 国产精品亚洲av一区麻豆| 久久天堂一区二区三区四区| 韩国av一区二区三区四区| 麻豆av在线久日| 久久精品人妻少妇| 国内久久婷婷六月综合欲色啪| 欧美日韩中文字幕国产精品一区二区三区| 国语自产精品视频在线第100页| 国产成人aa在线观看| 久久精品aⅴ一区二区三区四区| 麻豆一二三区av精品| 最近最新免费中文字幕在线| 国产精品久久久久久亚洲av鲁大| 天堂影院成人在线观看| 偷拍熟女少妇极品色| 日本黄色视频三级网站网址| 美女高潮的动态| 91麻豆精品激情在线观看国产| 天堂av国产一区二区熟女人妻| 成人午夜高清在线视频| 999久久久国产精品视频| 亚洲欧美精品综合久久99| 亚洲精品美女久久久久99蜜臀| 岛国在线观看网站| 麻豆成人午夜福利视频| 亚洲片人在线观看| 国产成人精品久久二区二区免费| 丝袜人妻中文字幕| 久久国产乱子伦精品免费另类| xxx96com| 最近最新中文字幕大全免费视频| 日本免费一区二区三区高清不卡| 国内久久婷婷六月综合欲色啪| 国产精品自产拍在线观看55亚洲| 欧美三级亚洲精品| 老司机午夜福利在线观看视频| 99在线人妻在线中文字幕| 成年版毛片免费区| 亚洲人成伊人成综合网2020| 日韩欧美三级三区| 国内精品一区二区在线观看| 操出白浆在线播放| 久久久成人免费电影| 亚洲一区二区三区不卡视频| 好男人在线观看高清免费视频| 国内揄拍国产精品人妻在线| 9191精品国产免费久久| 国产亚洲精品久久久com| www日本黄色视频网| 99国产精品99久久久久| 成人无遮挡网站| 色在线成人网| 99热这里只有是精品50| 成人午夜高清在线视频| 男人舔奶头视频| 两个人看的免费小视频| 最近最新中文字幕大全免费视频| 大型黄色视频在线免费观看| 亚洲天堂国产精品一区在线| 日本免费a在线| 精品电影一区二区在线| 亚洲第一电影网av| 免费大片18禁| 日本黄色视频三级网站网址| 欧美极品一区二区三区四区| 亚洲av日韩精品久久久久久密| 亚洲熟女毛片儿| 成人三级做爰电影| 日韩中文字幕欧美一区二区| 色视频www国产| 国产成人福利小说| 国产视频一区二区在线看|