• <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
    91大片在线观看| 久久精品国产亚洲av高清一级| 最近最新免费中文字幕在线| 狂野欧美激情性xxxx| 丰满迷人的少妇在线观看| 日韩国内少妇激情av| 国产99白浆流出| 亚洲精品在线观看二区| 人人妻人人爽人人添夜夜欢视频| 久久精品亚洲av国产电影网| 亚洲avbb在线观看| 在线观看舔阴道视频| 一区二区三区激情视频| 久久国产乱子伦精品免费另类| 深夜精品福利| 久久精品亚洲av国产电影网| 欧美大码av| 免费日韩欧美在线观看| 在线观看一区二区三区激情| 天天添夜夜摸| aaaaa片日本免费| xxx96com| 国产精品98久久久久久宅男小说| 亚洲国产精品999在线| 91精品三级在线观看| 麻豆av在线久日| 色综合欧美亚洲国产小说| 午夜日韩欧美国产| 男女午夜视频在线观看| 欧美日韩乱码在线| 亚洲成a人片在线一区二区| xxx96com| 日韩中文字幕欧美一区二区| 亚洲激情在线av| 免费观看人在逋| 午夜老司机福利片| 搡老熟女国产l中国老女人| 精品国产超薄肉色丝袜足j| avwww免费| 国产亚洲精品久久久久5区| 99久久国产精品久久久| 黄片播放在线免费| 少妇粗大呻吟视频| 色综合站精品国产| 18禁黄网站禁片午夜丰满| 亚洲欧美日韩高清在线视频| 亚洲人成网站在线播放欧美日韩| 日韩欧美一区二区三区在线观看| 啪啪无遮挡十八禁网站| av网站在线播放免费| 午夜福利影视在线免费观看| 啦啦啦免费观看视频1| 午夜福利免费观看在线| 性欧美人与动物交配| 国产精品野战在线观看 | 亚洲精品美女久久久久99蜜臀| cao死你这个sao货| 91麻豆av在线| 欧美人与性动交α欧美软件| 一区福利在线观看| 中亚洲国语对白在线视频| 久久草成人影院| 啦啦啦在线免费观看视频4| 欧美日韩黄片免| 午夜福利,免费看| 18禁裸乳无遮挡免费网站照片 | 男女下面插进去视频免费观看| 亚洲成国产人片在线观看| 亚洲国产欧美网| 好男人电影高清在线观看| 一二三四社区在线视频社区8| 一区二区三区精品91| 成人影院久久| 超碰97精品在线观看| 天堂动漫精品| 中文字幕人妻熟女乱码| 国产激情欧美一区二区| 天天躁夜夜躁狠狠躁躁| 国内毛片毛片毛片毛片毛片| 夜夜夜夜夜久久久久| 日日摸夜夜添夜夜添小说| 丰满饥渴人妻一区二区三| 1024香蕉在线观看| 波多野结衣一区麻豆| xxxhd国产人妻xxx| 97碰自拍视频| 国产一区二区三区在线臀色熟女 | 男人舔女人下体高潮全视频| 欧美日韩乱码在线| 人妻久久中文字幕网| 美国免费a级毛片| 国产野战对白在线观看| 国产有黄有色有爽视频| 国产成人啪精品午夜网站| 手机成人av网站| 亚洲精品久久成人aⅴ小说| 中亚洲国语对白在线视频| 亚洲avbb在线观看| 国产精品免费视频内射| 十八禁网站免费在线| 最近最新免费中文字幕在线| www.精华液| 9热在线视频观看99| 亚洲精品久久成人aⅴ小说| 老鸭窝网址在线观看| 午夜成年电影在线免费观看| 黑人巨大精品欧美一区二区mp4| 一区福利在线观看| 丝袜人妻中文字幕| 丁香六月欧美| 亚洲国产精品999在线| av有码第一页| 老司机午夜福利在线观看视频| 香蕉久久夜色| 日韩 欧美 亚洲 中文字幕| 日韩av在线大香蕉| 国产av一区在线观看免费| 777久久人妻少妇嫩草av网站| 80岁老熟妇乱子伦牲交| 在线观看日韩欧美| 久久精品国产综合久久久| 色在线成人网| 久久久久久久精品吃奶| 国产人伦9x9x在线观看| 在线观看www视频免费| 国产一区二区三区综合在线观看| 亚洲一区高清亚洲精品| 国产成+人综合+亚洲专区| 涩涩av久久男人的天堂| 日韩人妻精品一区2区三区| 成人影院久久| av中文乱码字幕在线| 国产成人精品在线电影| 国产97色在线日韩免费| 欧美精品亚洲一区二区| 99在线人妻在线中文字幕| 99国产精品一区二区蜜桃av| 亚洲精品一二三| 国产精品偷伦视频观看了| xxxhd国产人妻xxx| 免费少妇av软件| 母亲3免费完整高清在线观看| 美女扒开内裤让男人捅视频| 欧美日韩亚洲高清精品| 51午夜福利影视在线观看| 色尼玛亚洲综合影院| 久久国产精品影院| 99在线人妻在线中文字幕| 国产精品美女特级片免费视频播放器 | 少妇裸体淫交视频免费看高清 | www.自偷自拍.com| 久久亚洲真实| 亚洲欧美日韩高清在线视频| 免费看a级黄色片| 国内毛片毛片毛片毛片毛片| 婷婷丁香在线五月| 丝袜在线中文字幕| 亚洲一码二码三码区别大吗| 国产激情欧美一区二区| 亚洲自拍偷在线| 久久精品国产亚洲av高清一级| 国产av一区在线观看免费| 夜夜看夜夜爽夜夜摸 | 一级a爱视频在线免费观看| 88av欧美| 91av网站免费观看| 婷婷丁香在线五月| 欧美一区二区精品小视频在线| 午夜福利免费观看在线| 伦理电影免费视频| 国产av在哪里看| 国内毛片毛片毛片毛片毛片| 日本三级黄在线观看| 亚洲人成网站在线播放欧美日韩| 91成年电影在线观看| 在线av久久热| 999久久久精品免费观看国产| 国产一区二区三区综合在线观看| 亚洲男人的天堂狠狠| 亚洲欧美日韩另类电影网站| 一个人免费在线观看的高清视频| 精品午夜福利视频在线观看一区| 国产99白浆流出| 五月开心婷婷网| 狂野欧美激情性xxxx| 黑人巨大精品欧美一区二区mp4| 日韩欧美三级三区| 精品国产一区二区久久| 日韩精品中文字幕看吧| 亚洲av美国av| 99久久综合精品五月天人人| 夜夜看夜夜爽夜夜摸 | 久久久久久久久免费视频了| 亚洲自拍偷在线| 一个人观看的视频www高清免费观看 | 欧美丝袜亚洲另类 | 国产精品久久久av美女十八| 国产熟女午夜一区二区三区| av免费在线观看网站| 亚洲av美国av| 色婷婷av一区二区三区视频| 国产亚洲av高清不卡| 国产精品日韩av在线免费观看 | 老熟妇乱子伦视频在线观看| 国产一区二区在线av高清观看| 桃色一区二区三区在线观看| 9热在线视频观看99| 亚洲精品美女久久av网站| 欧美丝袜亚洲另类 | 欧美亚洲日本最大视频资源| 人成视频在线观看免费观看| 看黄色毛片网站| 在线视频色国产色| 成人手机av| netflix在线观看网站| 两个人免费观看高清视频| 99在线视频只有这里精品首页| 在线观看舔阴道视频| 久久国产精品男人的天堂亚洲| 久久精品亚洲av国产电影网| 香蕉久久夜色| 88av欧美| 日本撒尿小便嘘嘘汇集6| 免费日韩欧美在线观看| 91成年电影在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 激情视频va一区二区三区| 国内毛片毛片毛片毛片毛片| 夜夜躁狠狠躁天天躁| 神马国产精品三级电影在线观看 | 12—13女人毛片做爰片一| 国产精品自产拍在线观看55亚洲| 国产片内射在线| 国产精品影院久久| 午夜免费激情av| 久久国产亚洲av麻豆专区| 午夜福利欧美成人| 女警被强在线播放| 国产精品99久久99久久久不卡| 欧美日韩中文字幕国产精品一区二区三区 | 中文字幕人妻熟女乱码| 亚洲午夜精品一区,二区,三区| 亚洲熟妇熟女久久| 99久久综合精品五月天人人| 999久久久国产精品视频| 精品福利永久在线观看| 咕卡用的链子| 午夜福利,免费看| 91麻豆av在线| 老熟妇仑乱视频hdxx| 狠狠狠狠99中文字幕| 亚洲情色 制服丝袜| 亚洲自拍偷在线| 十八禁人妻一区二区| 欧美乱码精品一区二区三区| 亚洲精品粉嫩美女一区| 丰满人妻熟妇乱又伦精品不卡| 乱人伦中国视频| 十八禁人妻一区二区| 又黄又粗又硬又大视频| 国产黄a三级三级三级人| 中出人妻视频一区二区| 老熟妇仑乱视频hdxx| 最近最新中文字幕大全免费视频| 精品国内亚洲2022精品成人| 国产精品香港三级国产av潘金莲| 久9热在线精品视频| 国产精品国产高清国产av| 亚洲成av片中文字幕在线观看| 99久久综合精品五月天人人| 成人国产一区最新在线观看| 99riav亚洲国产免费| 亚洲国产欧美网| 午夜福利,免费看| 亚洲欧美精品综合久久99| 欧美日韩黄片免| 久久久久久久久久久久大奶| 色综合婷婷激情| 午夜福利一区二区在线看| 性少妇av在线| 日韩免费av在线播放| av网站免费在线观看视频| 亚洲av成人一区二区三| 国产精品一区二区三区四区久久 | av网站免费在线观看视频| 午夜91福利影院| 日本精品一区二区三区蜜桃| 亚洲精品国产一区二区精华液| 久久欧美精品欧美久久欧美| 久久午夜亚洲精品久久| 亚洲色图av天堂| 久久久国产精品麻豆| 国产aⅴ精品一区二区三区波| 男女午夜视频在线观看| 亚洲人成伊人成综合网2020| 欧美日韩黄片免| 色精品久久人妻99蜜桃| 色老头精品视频在线观看| 欧美激情极品国产一区二区三区| 国产精品一区二区免费欧美| 天天躁狠狠躁夜夜躁狠狠躁| 极品人妻少妇av视频| 热re99久久精品国产66热6| 久久亚洲真实| 超碰成人久久| 久久人妻熟女aⅴ| 真人做人爱边吃奶动态| 欧美日本中文国产一区发布| 黄片播放在线免费| 欧美激情 高清一区二区三区| 亚洲精品国产一区二区精华液| 涩涩av久久男人的天堂| 99久久精品国产亚洲精品| 少妇粗大呻吟视频| 欧美人与性动交α欧美精品济南到| 九色亚洲精品在线播放| 亚洲精品中文字幕一二三四区| 中文字幕av电影在线播放| 在线免费观看的www视频| 精品久久蜜臀av无| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品亚洲av一区麻豆| 自线自在国产av| 大香蕉久久成人网| 亚洲五月色婷婷综合| 狠狠狠狠99中文字幕| 成人精品一区二区免费| 99国产精品免费福利视频| 中文字幕最新亚洲高清| 80岁老熟妇乱子伦牲交| 精品久久久久久成人av| 美女大奶头视频| 女警被强在线播放| 好男人电影高清在线观看| 国产熟女午夜一区二区三区| 悠悠久久av| 亚洲人成电影观看| 天堂俺去俺来也www色官网| av网站免费在线观看视频| 99国产精品一区二区三区| 成在线人永久免费视频| 国产亚洲精品综合一区在线观看 | 免费日韩欧美在线观看| 久久中文字幕一级| 欧美人与性动交α欧美软件| 亚洲精品中文字幕一二三四区| 精品高清国产在线一区| 亚洲精品美女久久av网站| 国产人伦9x9x在线观看| 宅男免费午夜| 国产日韩一区二区三区精品不卡| 免费在线观看完整版高清| 中出人妻视频一区二区| 夜夜看夜夜爽夜夜摸 | 欧美日韩国产mv在线观看视频| 国产三级黄色录像| 在线观看免费高清a一片| 日韩三级视频一区二区三区| 午夜免费鲁丝| 亚洲精品国产精品久久久不卡| 91大片在线观看| av欧美777| av网站在线播放免费| 国产又爽黄色视频| www.精华液| 一区福利在线观看| 伦理电影免费视频| 人人妻人人爽人人添夜夜欢视频| 亚洲精品久久午夜乱码| 国产黄a三级三级三级人| 久久国产精品人妻蜜桃| 日本精品一区二区三区蜜桃| 亚洲欧洲精品一区二区精品久久久| 国产单亲对白刺激| 神马国产精品三级电影在线观看 | avwww免费| 亚洲avbb在线观看| av在线播放免费不卡| 久久久久久亚洲精品国产蜜桃av| 美女大奶头视频| 色婷婷av一区二区三区视频| 国产成人系列免费观看| 高清av免费在线| 久久狼人影院| 在线播放国产精品三级| 久久精品国产亚洲av香蕉五月| www国产在线视频色| 久久久久亚洲av毛片大全| 99香蕉大伊视频| 99久久久亚洲精品蜜臀av| 成年人黄色毛片网站| 在线观看www视频免费| 欧美乱妇无乱码| 国产亚洲精品久久久久久毛片| a级毛片在线看网站| 国产精品免费视频内射| 91国产中文字幕| 国产一区二区三区视频了| 这个男人来自地球电影免费观看| 国产99白浆流出| 亚洲五月婷婷丁香| 一进一出抽搐gif免费好疼 | 日韩精品免费视频一区二区三区| 亚洲精品国产色婷婷电影| 免费av毛片视频| svipshipincom国产片| av视频免费观看在线观看| 国内久久婷婷六月综合欲色啪| 国产不卡一卡二| 中出人妻视频一区二区| 精品一区二区三卡| 丰满人妻熟妇乱又伦精品不卡| 亚洲视频免费观看视频| av国产精品久久久久影院| 女性生殖器流出的白浆| 又黄又粗又硬又大视频| a级毛片在线看网站| 亚洲欧美激情在线| 欧美色视频一区免费| 国产99白浆流出| 久久性视频一级片| 十八禁网站免费在线| 午夜免费鲁丝| 淫秽高清视频在线观看| 精品久久久精品久久久| 久久久久久久精品吃奶| 黄色成人免费大全| 一边摸一边抽搐一进一小说| 黄网站色视频无遮挡免费观看| 成人免费观看视频高清| 日韩欧美三级三区| 女警被强在线播放| av中文乱码字幕在线| 久久久精品国产亚洲av高清涩受| 婷婷精品国产亚洲av在线| 人人妻人人澡人人看| 国产亚洲精品第一综合不卡| 亚洲精品国产一区二区精华液| 久久天堂一区二区三区四区| 欧美最黄视频在线播放免费 | 操美女的视频在线观看| 国产又色又爽无遮挡免费看| 亚洲avbb在线观看| 国产欧美日韩一区二区三| 正在播放国产对白刺激| 在线免费观看的www视频| 亚洲情色 制服丝袜| 亚洲精品国产精品久久久不卡| 午夜福利一区二区在线看| 亚洲欧美精品综合久久99| 亚洲国产精品999在线| 91麻豆av在线| 丰满的人妻完整版| 淫妇啪啪啪对白视频| 99在线人妻在线中文字幕| 亚洲精品av麻豆狂野| 成在线人永久免费视频| 欧美日韩精品网址| 欧美成狂野欧美在线观看| 韩国精品一区二区三区| 欧美成人免费av一区二区三区| 国产不卡一卡二| 久久伊人香网站| 国产一卡二卡三卡精品| 日日夜夜操网爽| 大陆偷拍与自拍| 亚洲av五月六月丁香网| 嫁个100分男人电影在线观看| 天堂俺去俺来也www色官网| 亚洲人成电影免费在线| 搡老岳熟女国产| 男女高潮啪啪啪动态图| 国产精品九九99| 香蕉久久夜色| 亚洲黑人精品在线| 丝袜在线中文字幕| 国产成人欧美在线观看| 亚洲精品久久成人aⅴ小说| 淫妇啪啪啪对白视频| 国产一区二区激情短视频| 精品高清国产在线一区| 人成视频在线观看免费观看| 麻豆久久精品国产亚洲av | 日本 av在线| 国产高清国产精品国产三级| 视频区欧美日本亚洲| 国产精品久久久久成人av| 久久亚洲真实| 999久久久精品免费观看国产| 欧美日韩福利视频一区二区| 成人三级黄色视频| 日本三级黄在线观看| 99热只有精品国产| 黄色怎么调成土黄色| av在线天堂中文字幕 | 久久人妻熟女aⅴ| 亚洲精品成人av观看孕妇| 一级作爱视频免费观看| 亚洲av成人不卡在线观看播放网| 国产区一区二久久| 99re在线观看精品视频| 成人18禁在线播放| 日本精品一区二区三区蜜桃| 麻豆一二三区av精品| 欧美激情高清一区二区三区| 变态另类成人亚洲欧美熟女 | 成人精品一区二区免费| 免费观看人在逋| 午夜两性在线视频| 精品国产一区二区久久| 男人舔女人下体高潮全视频| 国产日韩一区二区三区精品不卡| 淫秽高清视频在线观看| 国产精品乱码一区二三区的特点 | 免费观看精品视频网站| 一个人免费在线观看的高清视频| 超色免费av| 波多野结衣高清无吗| 亚洲欧美精品综合久久99| 一进一出好大好爽视频| 精品国产乱子伦一区二区三区| 久久精品人人爽人人爽视色| 成年版毛片免费区| 一区二区三区精品91| av天堂在线播放| 高清在线国产一区| 日本黄色日本黄色录像| 老司机深夜福利视频在线观看| 亚洲黑人精品在线| 成人av一区二区三区在线看| 亚洲成a人片在线一区二区| 国产成人系列免费观看| 国产aⅴ精品一区二区三区波| 美女高潮到喷水免费观看| 热99国产精品久久久久久7| 久久久久国内视频| 悠悠久久av| av网站在线播放免费| 长腿黑丝高跟| 日韩免费高清中文字幕av| 一区二区三区精品91| 91成人精品电影| 国产精品一区二区在线不卡| 99精国产麻豆久久婷婷| 中文亚洲av片在线观看爽| 精品国产一区二区久久| 少妇被粗大的猛进出69影院| 91精品国产国语对白视频| 一级片免费观看大全| 纯流量卡能插随身wifi吗| 一级片'在线观看视频| 亚洲精品粉嫩美女一区| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产国语对白av| 日韩高清综合在线| 国产深夜福利视频在线观看| 国产av一区在线观看免费| 琪琪午夜伦伦电影理论片6080| 免费日韩欧美在线观看| 久久香蕉激情| 欧美av亚洲av综合av国产av| 国产在线精品亚洲第一网站| av在线天堂中文字幕 | av超薄肉色丝袜交足视频| 欧美亚洲日本最大视频资源| 成人三级做爰电影| 亚洲精品美女久久av网站| 亚洲一区中文字幕在线| bbb黄色大片| 国产免费av片在线观看野外av| 50天的宝宝边吃奶边哭怎么回事| 在线看a的网站| 999久久久国产精品视频| 国内久久婷婷六月综合欲色啪| 琪琪午夜伦伦电影理论片6080| 亚洲午夜精品一区,二区,三区| 黄色片一级片一级黄色片| 久久天堂一区二区三区四区| 日韩免费高清中文字幕av| av国产精品久久久久影院| 老鸭窝网址在线观看| 日本五十路高清| 国产在线精品亚洲第一网站| 欧美日韩亚洲综合一区二区三区_| 在线观看www视频免费| 色综合欧美亚洲国产小说| 成人免费观看视频高清| 一边摸一边做爽爽视频免费| 中文字幕色久视频| av有码第一页| 亚洲人成网站在线播放欧美日韩| 亚洲成a人片在线一区二区| 午夜免费观看网址| 咕卡用的链子| 欧美色视频一区免费| 啦啦啦免费观看视频1| 免费在线观看日本一区| 日韩人妻精品一区2区三区| 国产三级在线视频| 国产日韩一区二区三区精品不卡| а√天堂www在线а√下载| 老司机靠b影院| 免费av中文字幕在线| 香蕉丝袜av| 三上悠亚av全集在线观看| 高清在线国产一区| 国产精品自产拍在线观看55亚洲| 女同久久另类99精品国产91| aaaaa片日本免费| 少妇被粗大的猛进出69影院| 国产免费男女视频| 男人操女人黄网站| 国产蜜桃级精品一区二区三区| 在线国产一区二区在线| 日韩大码丰满熟妇|