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

    連作煙田健株與感染根腐病煙株根際土壤細(xì)菌群落多樣性研究

    2022-03-11 21:30:10敖金成李永梅李博
    江蘇農(nóng)業(yè)科學(xué) 2022年4期
    關(guān)鍵詞:根際土壤高通量測(cè)序

    敖金成 李永梅 李博

    摘要:為探明連作與感染根腐病煙株根際土壤細(xì)菌群落結(jié)構(gòu)變化,采用Illumina高通量技術(shù)對(duì)細(xì)菌16S rDNA的V3-V4區(qū)域進(jìn)行測(cè)序,對(duì)比研究連作煙田健株與感染根腐病煙株根際土壤細(xì)菌群落結(jié)構(gòu)。結(jié)果表明,變形菌門、放線菌門、浮霉菌門、酸桿菌門、芽單胞菌門、綠彎菌門是連作煙田健株和病株根際土壤主要的細(xì)菌優(yōu)勢(shì)種群,相對(duì)豐度累計(jì)總和為89.1%~93.9%,隨連作年限的延長(zhǎng)各優(yōu)勢(shì)種群發(fā)生趨向性變化;隨連作年限的延長(zhǎng),健株根際土壤細(xì)菌群落豐富度和多樣性呈先增后降趨勢(shì),病株根際土壤細(xì)菌群落豐富度和多樣性呈降低趨勢(shì),短期連作(2~4 年)煙田健株和病株根際土壤細(xì)菌群落組成較為相似,與連作8年和撂荒2年以上(CK)土壤樣本的細(xì)菌群落組成差異較大;冗余分析結(jié)果表明,土壤化學(xué)因子與細(xì)菌群落分布關(guān)聯(lián)緊密,其中pH值是影響土壤細(xì)菌群落分布的核心因子,其貢獻(xiàn)度為15.77%。綜上,連作和染病降低了土壤細(xì)菌群落豐富度和多樣性,優(yōu)勢(shì)菌群趨向性演化,有益菌群相對(duì)豐度降低。在生產(chǎn)實(shí)踐中,提高煙田土壤pH值,降低連作年限并及時(shí)清除煙株病殘?bào)w有利于維持土壤微生態(tài)環(huán)境的穩(wěn)定性。

    關(guān)鍵詞:高通量測(cè)序;連作煙田;細(xì)菌群落;冗余分析;根際土壤;多樣性分析

    中圖分類號(hào): S154.36;S154.37文獻(xiàn)標(biāo)志碼: A

    文章編號(hào):1002-1302(2022)04-0198-07

    收稿日期:2021-05-26

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(編號(hào):32060445、41807524);云南省自然科學(xué)基金(編號(hào):202001AU070114);廣西中煙工業(yè)有限責(zé)任公司科技項(xiàng)目(編號(hào):GXZYCX2019b004)。

    作者簡(jiǎn)介:敖金成(1984—),男,云南曲靖人,博士研究生,高級(jí)工程師,主要從事植物營(yíng)養(yǎng)與病害控制技術(shù)研究。E-mail:89693180@qq.com。

    通信作者:李永梅,博士,教授,博士生導(dǎo)師,主要從事植物營(yíng)養(yǎng)與病害控制技術(shù)研究。E-mail:youngmaylee@126.com。

    烤煙(Nicotiana tobacum L.)是云南重要的經(jīng)濟(jì)作物,是典型的忌連作作物。受土地資源短缺和其他經(jīng)濟(jì)作物種植面積的擠壓,云南烤煙連作現(xiàn)象十分嚴(yán)重,多個(gè)核心煙區(qū)烤煙連作年限多在10年以上。連作引發(fā)許多土壤問題[1-4],加上近年來煙葉秸稈就地粉碎還田帶來的弊端[5],致使連作煙田土傳病害危害、煙葉產(chǎn)量和品質(zhì)下降等問題持續(xù)加重。煙草病害一直是嚴(yán)重制約世界煙草生產(chǎn)和科研發(fā)展的重要難題[6]。在諸多的土傳病害中,煙草鐮刀菌根腐病是由鐮刀菌屬(Fusarium)病原菌引起的一種典型的土傳病害,常與黑脛病、青枯病、黑根腐病復(fù)合發(fā)生,使得該病具有隱秘性、易混淆、難防治、危害重等特點(diǎn),近年來全國(guó)各煙區(qū)發(fā)病率呈上升趨勢(shì)。探尋有效防控?zé)煵莞〉姆椒?,是促進(jìn)煙草產(chǎn)業(yè)可持續(xù)發(fā)展的關(guān)鍵。

    土壤細(xì)菌是土壤養(yǎng)分循環(huán)重要的驅(qū)動(dòng)者[7],具有重要的生態(tài)功能[8-9]。連作可以改變根際細(xì)菌的群落結(jié)構(gòu),降低土壤細(xì)菌物種數(shù)和群落α多樣性[10]。長(zhǎng)期連作嚴(yán)重影響了土壤理化性質(zhì),進(jìn)而顯著降低了土壤細(xì)菌群落和組成多樣性[11],并且土壤優(yōu)勢(shì)菌門的相對(duì)豐度和多樣性存在時(shí)空異質(zhì)性[12],反映了土壤微生態(tài)的復(fù)雜性,易受環(huán)境因子[13-14]、作物類型[15]、施肥[16]等因素的調(diào)控,因而了解土壤細(xì)菌群落結(jié)構(gòu)及多樣性對(duì)維持土壤生態(tài)功能具有重要意義。史普酉等研究指出,不同發(fā)病程度煙株根際細(xì)菌群落組成和結(jié)構(gòu)發(fā)生趨向性變化,土壤生物多樣性降低[17]。煙草鐮刀菌根腐病發(fā)病煙株根際真菌群落多樣性顯著降低,病原菌相對(duì)豐度顯著高于健康煙株[18]。鄭元仙等研究也指出,煙株根際土壤真菌群落結(jié)構(gòu)改變及物種多樣性降低是根腐病發(fā)生的重要特征[19]。然而,關(guān)于土壤細(xì)菌群落結(jié)構(gòu)及多樣性對(duì)連作條件下煙株感染根腐病的響應(yīng)特征鮮見報(bào)道。鑒于此,本研究以不同連作年限煙田健株和感染根腐病煙株根際土壤為研究對(duì)象,探索連作條件下健株和病株根際土壤細(xì)菌群落結(jié)構(gòu)及多樣性特征,以期為連作煙田煙草根腐病的科學(xué)防治提供參考依據(jù)。

    1 材料與方法

    1.1 連作土壤采集區(qū)概況

    試驗(yàn)點(diǎn)位于云南省紅河州瀘西縣(103°52′26″E,24°40′03″N),地處滇南低山丘陵烤煙區(qū),地勢(shì)平緩,多數(shù)栽煙區(qū)域海拔在1 100~1 600 m之間,熱量條件優(yōu)越,降水充沛,年均溫為16.0~18.0℃,烤煙大田生育期降水量在1 000.0 mm左右,是云南省典型的清甜香風(fēng)格煙葉核心產(chǎn)區(qū),植煙歷史悠久。2010—2020年間,試驗(yàn)地區(qū)域年均溫為15.5 ℃,年均降水量為810.0 mm,年均日照時(shí)數(shù)為2 042.0 h,海拔為1 803.0 m。

    1.2 樣品采集

    試驗(yàn)于2020年7月烤煙成熟期,選取2、4、8 年連作植煙地塊,隨機(jī)取健康(以下簡(jiǎn)稱為“健株”)和感染根腐病煙株(以下簡(jiǎn)稱為“病株”)各4株,去除地表雜物,然后將煙株整株挖起,去除根系外圍土壤,采用抖根法采集根際土[20],輕輕抖落并收集須根2 mm范圍內(nèi)的土壤。對(duì)照(CK)為同區(qū)域撂荒2年以上未種植烤煙的土壤。每株烤煙根際土為1次重復(fù),各4次重復(fù)。健株分別編號(hào)為HT2、HT4、HT8,病株分別編號(hào)為ST2、ST4、ST8。試驗(yàn)地塊均為同一農(nóng)戶自己栽種,植煙歷史可考證。供試烤煙品種為云煙87。

    1.3 檢測(cè)項(xiàng)目及方法

    1.3.1 土壤化學(xué)性狀 取樣后將樣品帶回實(shí)驗(yàn)室,土壤經(jīng)風(fēng)干、磨細(xì)、過篩后備用。土壤有機(jī)質(zhì)、堿解氮、速效磷、速效鉀含量和pH值的測(cè)定參照文獻(xiàn)[21]中的方法。

    1.3.2 土壤細(xì)菌高通量測(cè)序 用于高通量測(cè)序的樣本取樣后立即放入冰盒帶回實(shí)驗(yàn)室,于-80 ℃保存?zhèn)溆??;贗llumina高通量測(cè)序技術(shù)檢測(cè)分析連作條件下煙田健株和感染根腐病煙株根際土壤細(xì)菌群落結(jié)構(gòu)的豐度和多樣性。

    1.3.2.1 基因組DNA的提取 樣本由廣州基迪奧生物科技有限公司平臺(tái)(https://www. Omic-smart.com)檢測(cè)。采用HiPure Soil DNA Mini Kits(廣州美基,中國(guó),cat#3412)土壤DNA提取試劑盒進(jìn)行土壤微生物DNA提取。首先在2 mL 珠管中加入0.25~0.50 g土壤樣品和0.6 mL 緩沖溶液,利用渦旋儀(米歐mix-28+,廣州圍谷潤(rùn)儀器有限公司)充分勻漿裂解,然后70 ℃水浴10 min,并離心收集管壁上的液滴,然后加入200 μL 聚苯乙烯緩沖液 和150 μL 吸收溶液,渦旋混勻20 s,靜置5 min后離心5 min,取上清液至2 mL離心管,加入等體積5′-二磷酸鳥苷二鈉緩沖液,混勻,獲得DNA柱。然后用1%瓊脂糖(Biowest Agarose,西班牙)凝膠電泳(DYY-6C,北京市六一儀器廠)檢測(cè)基因組DNA的完整性及是否發(fā)生蛋白等污染,DNA樣品于 -80 ℃ 保存待用。

    1.3.2.2 目的基因擴(kuò)增及Illumina測(cè)序 采用具有條形碼標(biāo)記的引物341F:C C T A C G G G N G G T C W G C A G、806R:G G A C T A C H V G G G T A C T A A T對(duì)細(xì)菌16S rDNA 的V3-V4區(qū)進(jìn)行MiSeq擴(kuò)增子測(cè)序。第1輪擴(kuò)增體系包括:10×Buffer KOD 5 μL,2 mmol/L dNTPs 5.0 μL,25 mmol/L MgSO4 3.0 μL,引物341F(10 μmol/L) 1.5 μL,引物806 R (10 μmol/L) 1.5 μL,KOD酶 1.0 μL,模板X μL(100 ng),用蒸餾水定容至 50 μL。然后利用AMPure XP Beads 進(jìn)行PCR純化,純化后用Qubit 3.0定量。進(jìn)入第二輪擴(kuò)增,擴(kuò)增體系包括:10×Buffer KOD 5 μL,2 mmol/L dNTPs 5.0 μL,25 mmol/L MgSO4 3.0 μL,索引引物(10 μmol/L)1 μL,PCR通用引物(10 μmol/L) 1 μL,KOD酶1 μL,模板 X μL(100 ng),再用蒸餾水定容至 50 μL。使用AMPure XP Beads 對(duì)第二輪擴(kuò)增產(chǎn)物進(jìn)行純化,用ABI StepOnePlus Real-Time PCR System(Life Technologies,美國(guó))進(jìn)行定量,根據(jù)Novaseq 6000的樣品要求規(guī)范化操作后,采用PE250模式上機(jī)測(cè)試。

    1.4 數(shù)據(jù)分析

    采用Microsoft Excel 2016進(jìn)行數(shù)據(jù)處理,用IBM Statistics SPSS 19.0進(jìn)行方差分析。利用軟件平臺(tái)Usearch(version 7.0)對(duì)相似度在97%條件下的操作分類單元(OUT)進(jìn)行質(zhì)控拼接和Tag聚類去嵌合體,獲得OTU的豐度和OTU代表序列,生成稀釋曲線,由圖1可知,各樣本的稀釋曲線基本趨于平坦,說明本次測(cè)序數(shù)據(jù)量合理,能夠比較真實(shí)地反映土壤樣本的細(xì)菌群落。利用軟件Mothur(version v.1.30.1)計(jì)算反映群落豐富度的Sobs指數(shù)、Chao1指數(shù)、Ace指數(shù)和反映群落多樣性的Shannon-Wiener指數(shù)、Simpson指數(shù)。然后利用R語(yǔ)言vegan包進(jìn)行群落柱狀圖的統(tǒng)計(jì)和作圖,以及主坐標(biāo)分析(principal coordinate analysis,簡(jiǎn)稱PCoA)、非度量多維尺度分析(nonmetric multidimensional scaling,簡(jiǎn)稱NMDS)、非加權(quán)組平均法(unweighted pair-group method with arithmetic means,簡(jiǎn)稱UPGMA)分析,利用冗余分析(redundancy analysis,簡(jiǎn)稱RDA)評(píng)價(jià)環(huán)境因子與細(xì)菌群落分布的關(guān)聯(lián)性。

    2 結(jié)果與分析

    2.1 連作煙田健株和感染根腐病煙株根際土壤的化學(xué)性質(zhì)

    由表1可知,不同連作年限煙田健株和病株根際土壤樣本的有機(jī)質(zhì)、堿解氮、速效磷、速效鉀含量及pH值存在顯著(P<0.05)或極顯著(P<0.01)差異。與CK相比,健株和病株根際土壤樣本的pH值隨連作年限增加極顯著(P<0.01)降低,有機(jī)質(zhì)和堿解氮含量呈先升后降的趨勢(shì),速效磷和速效鉀含量呈極顯著(P<0.01)增加的趨勢(shì)。與CK相比,HT2、HT4和HT8樣本的pH值分別下降0.8、0.9、2.5,ST2、ST4和ST8樣本的pH值分別下降0.8、0.7、2.1。除連作8年煙田健株和病株根際土壤樣本的有機(jī)質(zhì)含量均低于CK外,其他處理顯著高于CK或與CK一致。總體看,連作2年和連作4年煙田健株和病株根際土壤樣本的土壤養(yǎng)分含量整體一致,而連作8年煙田健株和病株根際土壤樣本的土壤pH值和有機(jī)質(zhì)、堿解氮含量整體較低。健株根際土壤pH值整體低于染病煙株,但速效磷含量均明顯高于染病煙株,短期連作(2~4年)煙田健株根際土壤樣本的速效鉀含量明顯低于病株,連作8年的煙田則以健株稍高。

    2.2 連作煙田健株和感染根腐病煙株根際土壤細(xì)菌群落組成分析

    在門水平上,不同連作年限煙田健株和病株根際土壤細(xì)菌優(yōu)勢(shì)菌群相對(duì)豐度存在一定差異。由物種分布堆疊圖(圖2)和相對(duì)豐度值(表2)可以看出,變形菌門(Proteobactera)、放線菌門(Actinobacteria)、浮霉菌門(Planctomycetes)、酸桿菌門(Acidobacteria)、芽單胞菌門(Gemmatimonadetes)、綠彎菌門(Chloroflexi)是不同連作年限煙田健株和病株根際土壤的6個(gè)優(yōu)勢(shì)菌群,累積相對(duì)豐度總和為89.1%~93.9%。隨連作年限的延長(zhǎng),HT2、HT4和HT8樣本變形菌門和髕骨細(xì)菌門的相對(duì)豐度呈增加趨勢(shì),且均大于CK,放線菌門、芽單胞菌門的相對(duì)豐度呈先降后升趨勢(shì),浮霉菌門、酸桿菌門、擬桿菌門、己科河菌門的相對(duì)豐度呈先升后降趨勢(shì),綠彎菌門的相對(duì)豐度呈降低趨勢(shì)。ST2、ST4和ST8樣本的變形菌門、放線菌門、髕骨細(xì)菌門的相對(duì)豐度呈增加趨勢(shì),浮霉菌門、酸桿菌門、芽單胞菌門、綠彎菌門、擬桿菌門和己科河菌門的相對(duì)豐度呈降低趨勢(shì),僅疣微菌門的相對(duì)豐度呈先升后降趨勢(shì)。

    2.3 連作煙田健株和感染根腐病煙株根際土壤細(xì)菌群落的α多樣性分析

    Sobs指數(shù)是細(xì)菌群落豐富度的實(shí)際觀測(cè)值,Ace指數(shù)和Chao1指數(shù)用于反映細(xì)菌群落豐富度,Shannon-Wiener指數(shù)用于反映群落多樣性。由表3可知,連作煙田健株和病株根際土壤細(xì)菌群落豐富度和多樣性存在差異。HT2、HT4和HT8樣本的Sobs指數(shù)、Chao1指數(shù)和Ace指數(shù)及Shannon-Wiener指數(shù)均表現(xiàn)為隨連作年限的增加呈先增后降的趨勢(shì),說明隨連作年限的延長(zhǎng),連作煙田健株根際土壤細(xì)菌群落豐富度和多樣性先增后降,其中HT8樣本細(xì)菌群落豐富度和多樣性均明顯小于HT2和HT4樣本。ST2、ST4和ST8樣本的Sobs指數(shù)、Chao1指數(shù)和Ace指數(shù)及Shannon-Wiener指數(shù)均表現(xiàn)為隨連作年限的增加呈降低趨勢(shì),說明隨連作年限的延長(zhǎng),煙田病株根際土壤細(xì)菌群落的豐富度和多樣性降低。健株和病株根際土壤樣本間物種覆蓋度均無(wú)顯著性差異,均大于0.99,說明樣本測(cè)序飽和度均較高。

    2.4 連作煙田健株和感染根腐病煙株根際土壤細(xì)菌群落的β多樣性分析

    主成分分析(principal component analysis,簡(jiǎn)稱PCA)和非度量多維尺度分析是基于OTU列表的物種豐度信息反映樣本間群落組成距離關(guān)系的方法。PCA用于研究樣本間的組成距離關(guān)系,NMDS分析可以更好地反映生態(tài)學(xué)數(shù)據(jù)的非線性結(jié)構(gòu),脅強(qiáng)系數(shù)(stress)越?。ㄐ∮?.1較好),說明模型越可靠。從圖3-a可以看出,第1主成分(PC1)和第2主成分(PC2)對(duì)樣品差異的貢獻(xiàn)值為82.18%,且PC1對(duì)樣品差異的貢獻(xiàn)值最高。NMDS分析結(jié)果(圖3-b)顯示,stress為0.050,小于0.1,說明模型可靠性高。圖3-a和圖3-b均表現(xiàn)出短期連作(2~4年)煙田的健株和病株樣本距離較近,群落組成較相似,而與CK、HT8和ST8樣本的距離較遠(yuǎn),說明群落結(jié)構(gòu)差異較大。進(jìn)一步進(jìn)行Anosim檢驗(yàn),結(jié)果表明,相同連作年限煙田健株和病株樣本細(xì)菌群落組成無(wú)顯著差異,但均與CK樣本有顯著差異(P<0.05)。3種連作年限的煙田健株(HT2、HT4和HT8)樣本、病株(ST2、ST4和ST8)樣本與CK樣本之間,及HT2、HT4和HT8樣本之間,ST2、ST4和ST8樣本之間的細(xì)菌組成均有極顯著差異(P<0.01)。

    2.5 連作煙田健株和感染根腐病煙株根際土壤理化性質(zhì)對(duì)細(xì)菌群落的影響

    為了解土壤理化因子對(duì)連作煙田健株和病株根際土壤微生物群落分布的影響,在屬水平上進(jìn)行冗余分析。從圖4-a可以看出,在屬水平上,RDA的前2個(gè)軸總共解釋了95.60%的群落變化,第1個(gè)軸和第2個(gè)軸分別解釋了85.15%、10.45%的群落變化,說明觀察樣本聚類特征與分組一致。從圖 4-b 和表4可以看出,在門水平上,土壤pH值與變形菌門的相對(duì)豐度呈極顯著負(fù)相關(guān)關(guān)系(r=-0.505,P<0.01),與酸桿菌門(r=0.578,P<0.01)和綠彎菌門(r=0.581,P<0.01)的相對(duì)豐度呈極顯著正相關(guān)關(guān)系,與浮霉菌門(r=0.402,P<0.05)和裝甲菌門(r=0.429,P<0.05)的相對(duì)豐度呈顯著正相關(guān)關(guān)系,與藍(lán)藻細(xì)菌門的相對(duì)豐度呈顯著負(fù)相關(guān)關(guān)系(r=-0.417,P<0.05);土壤有機(jī)質(zhì)含量與擬桿菌門的相對(duì)豐度呈顯著正相關(guān)關(guān)系(r=0.379,P<0.05);土壤堿解氮含量與細(xì)菌群落的相對(duì)豐度無(wú)顯著相關(guān)性;速效磷含量與酸桿菌門(r=-0.395,P<0.05)和綠彎菌門(r=-0.407,P<0.05)的相對(duì)豐度呈顯著負(fù)相關(guān)關(guān)系;土壤速效鉀含量?jī)H與硝化螺旋菌門的相對(duì)豐度呈極顯著正相關(guān)關(guān)系(r=0.479,P<0.01)。從圖4-c可以看出,在屬水平上,根際土壤有機(jī)質(zhì)、速效磷、速效鉀含量及pH值對(duì)細(xì)菌群落的分布影響以pH值的貢獻(xiàn)度(15.77%)最大,堿解氮和速效鉀含量的貢獻(xiàn)度為負(fù)值,表示對(duì)根際土壤細(xì)菌物種分布無(wú)意義。

    3 討論與結(jié)論

    本研究結(jié)果表明,連作煙田健株與感染根腐病煙株根際土壤細(xì)菌群落的相對(duì)豐度及多樣性存在差異。變形菌門、放線菌門、浮霉菌門、酸桿菌門、芽單胞菌門、綠彎菌門為幾種連作煙田土壤健株和病株根際土壤的優(yōu)勢(shì)菌群。研究表明,變形菌門是最為普遍的細(xì)菌門[22],同時(shí)包含了大量的動(dòng)植物致病菌[23],其豐度的增加可能引起土壤環(huán)境抗逆性降低,植株發(fā)病率增加。本試驗(yàn)中,隨連作年限的增加,變形菌門的相對(duì)豐度呈明顯增加的趨勢(shì),且感染根腐病煙株根際土壤變形菌門的相對(duì)豐度明顯高于健株。試驗(yàn)中,除了連作2、4年的健株和連作2年的病株放線菌菌門相對(duì)豐度較低外,其余樣本均較高,且表現(xiàn)出病株大于健株。說明隨連作年限的延長(zhǎng),土壤變形菌門逐漸演化成優(yōu)勢(shì)菌群。絕大多數(shù)放線菌為腐生菌,但少數(shù)寄生性放線菌則能引起某些動(dòng)植物的病害[24],其含量的高低可作為土壤健康狀況的評(píng)價(jià)指標(biāo)[25]。隨連作年限的增加,總體上連作煙田健株和病株根煙株根際土壤酸桿菌門、綠彎菌門的相對(duì)豐度呈明顯降低趨勢(shì),且整體表現(xiàn)為病株下降幅度大于健株。有研究認(rèn)為,土壤細(xì)菌是環(huán)境變化的敏感指標(biāo),其組成和多樣性與氣溫、降水以及土壤pH值、有機(jī)碳和全氮等含量密切相關(guān)[26],土壤微生物群落結(jié)構(gòu)差異可指示土壤環(huán)境的變化[27]。岳思君等研究認(rèn)為,引起硒砂瓜連作障礙的發(fā)生主要因素是隨著連作時(shí)間的增加,放線菌門、綠彎菌門、酸桿菌門、厚壁菌門等有益微生物豐度下降造成的[28]。該研究說明隨連作年限的延長(zhǎng),土壤細(xì)菌群落結(jié)構(gòu)趨向性改變,土壤細(xì)菌性病原菌相對(duì)豐度可能增加,而細(xì)菌有益菌群落相對(duì)豐度降低。

    有研究表明,土壤中微生物多樣性越高、結(jié)構(gòu)越復(fù)雜,系統(tǒng)的穩(wěn)定性也越高[29]。試驗(yàn)中,烤煙短期連作(2~4年)煙田健株和病株根際土壤細(xì)菌群落豐度和多樣性較高,群落組成較為相近,與撂荒土壤(CK)的距離較遠(yuǎn),而連作8年煙田健株(ST8)和病株(ST8)根際土壤細(xì)菌群落相對(duì)豐度和多樣性均顯著減小。說明隨連作年限的延長(zhǎng),土壤微生態(tài)穩(wěn)定性降低,進(jìn)而影響土壤系統(tǒng)抗性。土傳病害的發(fā)生與土壤環(huán)境密切相關(guān)[30-31],土壤微生物與植物健康的相互關(guān)系研究備受關(guān)注。土壤微生物尤其是根際微生物的組成、結(jié)構(gòu)、多樣性和生態(tài)網(wǎng)絡(luò)關(guān)系與植物土傳病害的發(fā)生密切相關(guān)[32-35]。汪洋等綜述了農(nóng)田管理措施對(duì)土壤微生物的影響,認(rèn)為提高有機(jī)肥和磷肥的施用比例有利于土壤中微生物豐富度的提高和微生物量碳的積累[36]。納小凡等研究認(rèn)為,因?yàn)橥寥纏H值是影響微生物群落結(jié)構(gòu)組成最重要的因素之一[10],本試驗(yàn)中根際土壤pH值對(duì)健株和病株根際土壤細(xì)菌群落分布影響最大,說明pH值是影響細(xì)菌群落分布及豐富度的核心環(huán)境因子。本研究的供試土壤隨連作年限的延長(zhǎng),pH值明顯降低。因而針對(duì)連作植煙地塊,生產(chǎn)中可通過撒施生石灰,增施有機(jī)肥、生物質(zhì)炭、過磷酸鈣等方式以提高連作土壤pH值,進(jìn)而改善連作土壤微生態(tài)環(huán)境。

    綜上所述,本研究利用Illumina高通量測(cè)序技術(shù)探明了連作煙田健株和感染根腐病煙株根際土壤細(xì)菌群落結(jié)構(gòu)和多樣性的變化,為連作煙田烤煙根腐病的防治及土壤改良提供了借鑒。然而本研究只關(guān)注了連作煙田健株和病株根際細(xì)菌群落的豐富度、多樣性及其細(xì)菌群落與連作土壤理化因子的關(guān)聯(lián)性,未探明根腐病發(fā)病與相關(guān)細(xì)菌菌群的相關(guān)性,以及真菌群落豐富度、多樣性的變化情況,在今后的研究中將加以完善。

    參考文獻(xiàn):

    [1]孫敬國(guó),王昌軍,孫光偉,等. 連作年限對(duì)植煙根際土壤化感物質(zhì)積累的影響——以湖北黃棕壤煙田為例[J]. 土壤,2021,53(1):148-153.

    [2]程亞東,白羽祥,史普酉,等. 連作植煙土壤具有顯著積累特征的兩種酚酸的化感效應(yīng)評(píng)價(jià)[J]. 核農(nóng)學(xué)報(bào),2020,34(10):2307-2315.

    [3]岳冰冰,李 鑫,張會(huì)慧,等. 連作對(duì)黑龍江烤煙土壤生物功能多樣性的影響[J]. 土壤,2013,45(1):116-119.

    [4]何 川,劉國(guó)順,李祖良,等. 連作植煙土壤有機(jī)碳和酶活性的影響及其與土傳病害的關(guān)系[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報(bào),2011,45(6):701-705.

    [5]閆 寧,郭東鋒,姚忠達(dá),等. 煙桿還田對(duì)煙草生長(zhǎng)、產(chǎn)量、質(zhì)量及病毒病發(fā)生的影響[J]. 江西農(nóng)業(yè)學(xué)報(bào),2016,28(7):68-72.

    [6]高福宏,陳 靜,陶 瓊,等. 2011年昆明市煙草主要病害發(fā)生特點(diǎn)[J]. 植物保護(hù),2013,35(2):158-164.

    [7]Gabriel C,Paulo J P L T,Sur H P,et al. Root microbiota drive direct integration of phosphate stress and immunity[J]. Nature,2017,543(7646):513-518.

    [8]Li J Y,Zhang Q C,Li Yong,et al. Effects of long-term mowing on the fractions and chemical composition of soil organic matter in a semiarid grassland[J]. Biogeosciences,2017,14(10):2685-2696.

    [9]Bulgarelli D,Schlaeppi K,Spaepen S,et al. Structure and functions of the bacterial microbiota of plants[J]. Annual Review of Plant Biology,2013,64(1):807-838.

    [10]納小凡,鄭國(guó)旗,刑正操,等. 連作對(duì)再植枸杞根際細(xì)菌群落多樣性和群落結(jié)構(gòu)的影響[J]. 土壤學(xué)報(bào),2017,54(5):1280-1292.

    [11]寇智瑞,周鑫斌. 不同連作年限黃壤煙田土壤細(xì)菌群落的差異[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2020,26(3):511-521.

    [12]敖金成,李 博,閻 凱,等. 連作對(duì)云南典型煙區(qū)植煙土壤細(xì)菌群落多樣性的影響[J]. 農(nóng)業(yè)資源與環(huán)境學(xué)報(bào),2022,39(1):46-54.

    [13]Sheng R,Meng D L,Wu M N,et al. Effect of agricultural land use change on community composition of bacteria and ammonia oxidizers[J]. Journal of Soil and Sediments,2013,13(7):1246-1256.

    [14]Van Horn D J,Van Horn M L,Barrett J E,et al. Factors controlling soil microbial biomass and bacterial diversity and community composition in a cold desert ecosystem:Role of geographic scale[J]. PLoS One,2013,8(6):e66103.

    [15]Marques J M,da Silva T F,Vollu R E,et al. Plant age and genotype affect the bacterial community composition in the tuber rhizosphere of filed-grown sweet potato plants[J]. FENS Microbiology Ecology,2014,88(2):424-433.

    [16]顧美英,劉洪亮,李志強(qiáng),等. 新疆連作棉田施用生物炭對(duì)土壤養(yǎng)分及微生物群落多樣性的影響[J]. 中國(guó)農(nóng)業(yè)科學(xué),2014,47(20):4128-4138.

    [17]史普酉,楊成翠,賈 孟,等. 不同黑脛病發(fā)病程度下植煙根際土壤酶活性及細(xì)菌群落結(jié)構(gòu)差異比較[J]. 中國(guó)土壤與肥料,2020(1):179-187.

    [18]李朋發(fā),楊 龍,李桂龍,等. 基于FUNGuild的鐮刀菌根腐病發(fā)病煙株根際真菌群落研究[J]. 中國(guó)煙草學(xué)報(bào),2019,25(2):63-68.

    [19]鄭元仙,楊 敏,王繼明,等. 烤煙根腐病對(duì)煙株根際土壤真菌群落結(jié)構(gòu)的影響[J]. 中國(guó)煙草科學(xué),2021,42(2):50-55.

    [20]Epelde L,Becerril J M,Hernandez-Allica J,et al. Functional diversity as indicator of the recovery of soil health from Thlaspi caerulescens growth and metal phytoextraction[J]. Applied Soil Ecology,2008,39(3):299-310.

    [21]鮑士旦. 土壤農(nóng)化分析[M]. 北京:中國(guó)農(nóng)業(yè)出版社,2000.

    [22]許艷蕊,方志軍,盧曉平,等. 基于高通量測(cè)序技術(shù)分析使它隆對(duì)玉米土壤細(xì)菌多樣性的影響[J]. 微生物學(xué)報(bào),2017,57(7):985-993.

    [23]Wagg C,Dudenhffer J H,Widmer F,et al. Linking diversity,synchrony and stability in soil microbial communities[J]. Functional Ecology,2018,32(5):1280-1292.

    [24]李文均,職曉陽(yáng),唐蜀昆. 我國(guó)放線菌系統(tǒng)學(xué)研究歷史、現(xiàn)狀及未來發(fā)展趨勢(shì)[J]. 微生物學(xué)通報(bào),2013,40(10):1860-1873.

    [25]張亞培,張嫣蕓,李梅梅,等. 變形菌多重耐藥性與整合子相關(guān)性分析[J]. 中國(guó)抗生素雜質(zhì),2015,40(2):141-144.

    [26]Van Veen J A,Ladd J N,Amato M. Turnover of carbon and nitrogen through the microbial biomass in a sandy loam and a clay soil incubated wtih[14C(U)]glucose and[15N](NH4)2SO4 under different moisture regimes[J]. Soil Biology and Biochemistry,1985,17(6):747-756.

    [27]李 飛,劉振恒,賈甜華,等. 高寒濕地和草甸退化及恢復(fù)對(duì)土壤微生物碳代謝功能多樣性的影響[J]. 生態(tài)學(xué)報(bào),2018,38(17):6006-6015.

    [28]岳思君,馮翠娥,楊彥研,等. 不同連作年限硒砂瓜土壤細(xì)菌群落結(jié)構(gòu)特征[J]. 干旱地區(qū)農(nóng)業(yè)研究,2020,38(1):230-236.

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

    [30]王 弋,楊煥文,趙正雄,等. 不同抗性烤煙品種根際微生物數(shù)量及多樣性差異研究[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2012,18(2):451-458.

    [31]Benizri E,Piutti S,Verger S,et al. Replant diseases:bacterial community structure and diversity in peach rhizosphere as determined by metabolic and genetic fingerprinting[J]. Soil Biology & Biochemistry,2005,537(9):1738-1746.

    [32]Niu L,Chan J,Xiao Y,et al. Insight into the effects of different cropping systems on soil bacterial community and tobacco bacterial wilt rate[J]. J Basic Microbiol,2017,57(1):3-11.

    [33]Latz E,Eisenhauer N,Rall B C,et al. Unravelling linkages between plant community composition and the pathogen-suppressive potential of soils[J]. Scientific Reports,2016,6(1):23584.

    [34]Niu J,Rang Z,Zhang C,et al. The succession pattern of soil microbial communities and its relationship with tobacco bacterial wilt[J]. BMC Microbiology,2016,16(1):1-10.

    [35]李 赟,劉 迪,范如芹,等. 土壤改良劑的研究進(jìn)展[J]. 江蘇農(nóng)業(yè)科學(xué),2020,48(10):63-69.

    [36]汪 洋,楊殿林,王麗麗,等. 農(nóng)田管理措施對(duì)土壤有機(jī)碳周轉(zhuǎn)及微生物量的影響[J]. 農(nóng)業(yè)資源與環(huán)境學(xué)報(bào),2020,37(3):340-352.

    3761501908218

    猜你喜歡
    根際土壤高通量測(cè)序
    環(huán)狀RNA在疾病發(fā)生中的作用
    川明參輪作對(duì)煙地土壤微生物群落結(jié)構(gòu)的影響
    不同種植模式對(duì)烤煙根際土壤養(yǎng)分、酶活性及煙草生長(zhǎng)的影響
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    人參根際真菌群落多樣性及組成的變化
    LncRNAs作為miRNA的靶模擬物調(diào)節(jié)miRNA
    馬纓丹入侵對(duì)根際和非根際土壤理化性質(zhì)的影響
    小蓬竹根際土壤可培養(yǎng)細(xì)菌多樣性
    華重樓根際土芽胞桿菌多樣性研究
    廣藿香根際土壤微生物總DNA提取方法的優(yōu)化
    亚洲中文av在线| 极品教师在线免费播放| 一区二区三区激情视频| 黄色丝袜av网址大全| 啪啪无遮挡十八禁网站| 色在线成人网| 欧美 亚洲 国产 日韩一| 精品久久蜜臀av无| 亚洲国产精品一区二区三区在线| 99re在线观看精品视频| 757午夜福利合集在线观看| 免费不卡黄色视频| 国产精品 国内视频| 精品一品国产午夜福利视频| 妹子高潮喷水视频| 丝袜美足系列| 亚洲男人的天堂狠狠| 桃红色精品国产亚洲av| 人人妻人人爽人人添夜夜欢视频| 男人舔女人的私密视频| 岛国在线观看网站| 91在线观看av| 老司机福利观看| 夜夜爽天天搞| 天堂√8在线中文| 亚洲精品一区av在线观看| 成年版毛片免费区| 亚洲avbb在线观看| 成年人黄色毛片网站| 久久国产精品男人的天堂亚洲| 黄色a级毛片大全视频| 高清毛片免费观看视频网站 | 国产亚洲av高清不卡| 操出白浆在线播放| 女生性感内裤真人,穿戴方法视频| 欧美乱色亚洲激情| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品综合一区在线观看 | a级毛片在线看网站| 免费av毛片视频| 97碰自拍视频| 亚洲中文字幕日韩| 老司机深夜福利视频在线观看| svipshipincom国产片| 日本一区二区免费在线视频| 91av网站免费观看| 深夜精品福利| 老熟妇仑乱视频hdxx| 777久久人妻少妇嫩草av网站| 丁香六月欧美| aaaaa片日本免费| 亚洲精华国产精华精| 久久久久久大精品| 国产高清激情床上av| 成人av一区二区三区在线看| 99国产精品一区二区三区| 国产av一区在线观看免费| 老司机午夜十八禁免费视频| 人成视频在线观看免费观看| 亚洲av成人av| 午夜福利欧美成人| 久久人人97超碰香蕉20202| 亚洲自偷自拍图片 自拍| xxxhd国产人妻xxx| 大码成人一级视频| 国产精华一区二区三区| 国产免费现黄频在线看| 国产主播在线观看一区二区| 欧美日韩亚洲高清精品| 免费日韩欧美在线观看| 精品久久久久久成人av| 视频区图区小说| 国产精品 国内视频| 国产成人精品在线电影| 啦啦啦免费观看视频1| 亚洲欧美日韩另类电影网站| 久久 成人 亚洲| 欧美日韩亚洲国产一区二区在线观看| 69精品国产乱码久久久| 久久香蕉激情| 如日韩欧美国产精品一区二区三区| 淫秽高清视频在线观看| 80岁老熟妇乱子伦牲交| 每晚都被弄得嗷嗷叫到高潮| 美国免费a级毛片| 免费在线观看日本一区| 日本三级黄在线观看| 免费在线观看黄色视频的| 日韩大码丰满熟妇| 久久精品国产亚洲av高清一级| 90打野战视频偷拍视频| 1024香蕉在线观看| 午夜a级毛片| 乱人伦中国视频| 每晚都被弄得嗷嗷叫到高潮| 黄色成人免费大全| 色综合婷婷激情| 脱女人内裤的视频| 正在播放国产对白刺激| 精品久久久久久电影网| 波多野结衣一区麻豆| 国产精品久久视频播放| 狠狠狠狠99中文字幕| 久久久国产精品麻豆| 99热只有精品国产| 国产成人啪精品午夜网站| 亚洲九九香蕉| 女生性感内裤真人,穿戴方法视频| 久久中文看片网| 欧美黑人精品巨大| 国产av在哪里看| www.自偷自拍.com| 身体一侧抽搐| 精品国产亚洲在线| 少妇 在线观看| 亚洲午夜精品一区,二区,三区| 国产亚洲精品一区二区www| 国产aⅴ精品一区二区三区波| 最新在线观看一区二区三区| 国产麻豆69| 国产av又大| 久久精品91无色码中文字幕| 国产伦人伦偷精品视频| www.999成人在线观看| 丝袜在线中文字幕| 黑人巨大精品欧美一区二区mp4| 欧美精品啪啪一区二区三区| 精品高清国产在线一区| 可以在线观看毛片的网站| 高清在线国产一区| 法律面前人人平等表现在哪些方面| e午夜精品久久久久久久| 久99久视频精品免费| 欧美色视频一区免费| tocl精华| 日韩大码丰满熟妇| 视频在线观看一区二区三区| 久久亚洲真实| 久久99一区二区三区| 久久久精品欧美日韩精品| 88av欧美| 天堂√8在线中文| 亚洲av成人av| 亚洲av成人av| 十八禁网站免费在线| 黄片小视频在线播放| 久久中文看片网| 亚洲av成人av| 在线免费观看的www视频| 亚洲第一青青草原| 国产一卡二卡三卡精品| 日韩欧美一区二区三区在线观看| 久久人人精品亚洲av| 日韩成人在线观看一区二区三区| 成人国产一区最新在线观看| 波多野结衣高清无吗| 999精品在线视频| 精品国产超薄肉色丝袜足j| av免费在线观看网站| 亚洲五月色婷婷综合| 老司机靠b影院| 亚洲精品在线美女| 国产成人精品久久二区二区免费| 自拍欧美九色日韩亚洲蝌蚪91| av超薄肉色丝袜交足视频| 男女之事视频高清在线观看| 一夜夜www| 亚洲黑人精品在线| 精品福利永久在线观看| 国产黄色免费在线视频| 亚洲一区二区三区欧美精品| 麻豆av在线久日| 黄片大片在线免费观看| 99在线视频只有这里精品首页| 精品久久蜜臀av无| 欧洲精品卡2卡3卡4卡5卡区| 99久久综合精品五月天人人| 成熟少妇高潮喷水视频| 亚洲精品国产一区二区精华液| 国产精品 欧美亚洲| 亚洲精品国产一区二区精华液| 日本撒尿小便嘘嘘汇集6| 很黄的视频免费| 亚洲成国产人片在线观看| 成熟少妇高潮喷水视频| 国产成人啪精品午夜网站| 亚洲国产中文字幕在线视频| 美女扒开内裤让男人捅视频| 国产成人啪精品午夜网站| 正在播放国产对白刺激| 国产日韩一区二区三区精品不卡| 9色porny在线观看| 夜夜躁狠狠躁天天躁| 国产高清国产精品国产三级| 男女高潮啪啪啪动态图| 91老司机精品| 黑人欧美特级aaaaaa片| 怎么达到女性高潮| 亚洲欧美一区二区三区久久| 美国免费a级毛片| 欧美人与性动交α欧美软件| 不卡av一区二区三区| 精品少妇一区二区三区视频日本电影| 精品国产亚洲在线| 18禁观看日本| 国产三级在线视频| 日韩 欧美 亚洲 中文字幕| 国产欧美日韩一区二区三| 亚洲成av片中文字幕在线观看| 国产成年人精品一区二区 | 欧美中文日本在线观看视频| 日韩欧美在线二视频| 欧美久久黑人一区二区| 欧美最黄视频在线播放免费 | 欧美亚洲日本最大视频资源| 好男人电影高清在线观看| 18禁观看日本| 国产在线观看jvid| 80岁老熟妇乱子伦牲交| 成年人免费黄色播放视频| 免费观看人在逋| 母亲3免费完整高清在线观看| 欧美精品一区二区免费开放| 国产精品久久久久久人妻精品电影| 一二三四社区在线视频社区8| 美女高潮到喷水免费观看| 亚洲av第一区精品v没综合| 又大又爽又粗| 国产国语露脸激情在线看| 夜夜夜夜夜久久久久| 国产激情久久老熟女| 人人妻人人添人人爽欧美一区卜| 一区二区三区国产精品乱码| 在线国产一区二区在线| 国产精品一区二区精品视频观看| 精品国产超薄肉色丝袜足j| 亚洲精品在线观看二区| 久久亚洲精品不卡| 免费高清视频大片| 91av网站免费观看| 精品福利永久在线观看| 欧美色视频一区免费| 最新在线观看一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕高清在线视频| 亚洲一区中文字幕在线| 亚洲国产欧美日韩在线播放| 日韩精品青青久久久久久| 欧美在线黄色| 女生性感内裤真人,穿戴方法视频| 日本a在线网址| 青草久久国产| 成人国语在线视频| 国产人伦9x9x在线观看| www日本在线高清视频| 国产精品98久久久久久宅男小说| 日本黄色视频三级网站网址| 亚洲精品在线观看二区| 婷婷六月久久综合丁香| 国产欧美日韩精品亚洲av| 国产伦人伦偷精品视频| 一边摸一边抽搐一进一小说| 国产熟女午夜一区二区三区| 亚洲成国产人片在线观看| 国产成人啪精品午夜网站| 日日摸夜夜添夜夜添小说| 欧美日韩精品网址| 久久人人97超碰香蕉20202| 亚洲第一av免费看| 久久久水蜜桃国产精品网| 国产精品久久久久久人妻精品电影| 欧美日韩视频精品一区| 亚洲成国产人片在线观看| 久久久久久久久免费视频了| 成人永久免费在线观看视频| 成人av一区二区三区在线看| 男人舔女人下体高潮全视频| 19禁男女啪啪无遮挡网站| 中亚洲国语对白在线视频| x7x7x7水蜜桃| 国产亚洲精品综合一区在线观看 | 国产99久久九九免费精品| 亚洲色图av天堂| 精品人妻1区二区| 深夜精品福利| 国产成人精品久久二区二区91| 久久香蕉激情| 免费久久久久久久精品成人欧美视频| 亚洲av片天天在线观看| 欧美成人午夜精品| 搡老岳熟女国产| 国产精品爽爽va在线观看网站 | √禁漫天堂资源中文www| 午夜福利影视在线免费观看| 亚洲专区国产一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av成人av| 久久精品91无色码中文字幕| 极品人妻少妇av视频| 丰满人妻熟妇乱又伦精品不卡| 丝袜美足系列| 激情视频va一区二区三区| 久久狼人影院| 国产精品自产拍在线观看55亚洲| 嫁个100分男人电影在线观看| 欧美成人免费av一区二区三区| 在线观看免费高清a一片| 777久久人妻少妇嫩草av网站| www.自偷自拍.com| 国产精品二区激情视频| 亚洲黑人精品在线| www国产在线视频色| 中文字幕高清在线视频| 亚洲精品中文字幕在线视频| 国产av一区二区精品久久| 中文字幕另类日韩欧美亚洲嫩草| 国产欧美日韩一区二区精品| 正在播放国产对白刺激| 亚洲va日本ⅴa欧美va伊人久久| 国产欧美日韩一区二区三区在线| 欧美午夜高清在线| 日本a在线网址| 免费在线观看亚洲国产| 久久青草综合色| 国产亚洲av高清不卡| 精品人妻在线不人妻| 国产精品av久久久久免费| 国产欧美日韩综合在线一区二区| 窝窝影院91人妻| 亚洲三区欧美一区| 一区福利在线观看| 亚洲av第一区精品v没综合| 视频区图区小说| 国产精品电影一区二区三区| 丰满饥渴人妻一区二区三| 久久狼人影院| 免费高清视频大片| 日韩高清综合在线| 日韩欧美在线二视频| 激情在线观看视频在线高清| 无遮挡黄片免费观看| 精品久久久久久久毛片微露脸| 麻豆国产av国片精品| 最近最新中文字幕大全电影3 | 午夜福利欧美成人| 久久精品国产亚洲av高清一级| 首页视频小说图片口味搜索| 在线观看免费高清a一片| 亚洲精品中文字幕在线视频| 久久中文看片网| 一级,二级,三级黄色视频| 国产激情久久老熟女| 老熟妇乱子伦视频在线观看| 亚洲自拍偷在线| 91av网站免费观看| 欧美黄色淫秽网站| 精品卡一卡二卡四卡免费| 嫩草影视91久久| 一个人免费在线观看的高清视频| 别揉我奶头~嗯~啊~动态视频| 人人妻人人添人人爽欧美一区卜| 午夜两性在线视频| 在线免费观看的www视频| 国产成人影院久久av| 伦理电影免费视频| 久热爱精品视频在线9| 国产不卡一卡二| 亚洲黑人精品在线| 国产亚洲欧美精品永久| 欧美精品啪啪一区二区三区| 丁香欧美五月| 久久精品国产亚洲av高清一级| 免费不卡黄色视频| av电影中文网址| 91成人精品电影| 大陆偷拍与自拍| 亚洲成a人片在线一区二区| 欧美乱码精品一区二区三区| 国产精品影院久久| 国产免费现黄频在线看| 精品一区二区三区四区五区乱码| 99久久99久久久精品蜜桃| 精品人妻1区二区| 欧美激情高清一区二区三区| 中文字幕色久视频| 午夜激情av网站| 精品人妻1区二区| 亚洲国产看品久久| 国产有黄有色有爽视频| 视频在线观看一区二区三区| 国产精品一区二区在线不卡| 亚洲熟女毛片儿| 深夜精品福利| av福利片在线| 午夜日韩欧美国产| 免费看a级黄色片| 天堂中文最新版在线下载| 精品日产1卡2卡| 少妇 在线观看| 国内毛片毛片毛片毛片毛片| 欧美成人免费av一区二区三区| 欧美不卡视频在线免费观看 | 国产xxxxx性猛交| 免费搜索国产男女视频| 亚洲精品国产色婷婷电影| 午夜两性在线视频| 天堂动漫精品| 视频区图区小说| 无限看片的www在线观看| 日韩精品免费视频一区二区三区| 日韩人妻精品一区2区三区| 天堂影院成人在线观看| av免费在线观看网站| 精品乱码久久久久久99久播| 免费在线观看黄色视频的| 91在线观看av| 国产真人三级小视频在线观看| 电影成人av| 久久人人爽av亚洲精品天堂| 每晚都被弄得嗷嗷叫到高潮| 亚洲av电影在线进入| www国产在线视频色| 欧美精品亚洲一区二区| 亚洲男人天堂网一区| av福利片在线| 中文字幕精品免费在线观看视频| av在线播放免费不卡| 十八禁人妻一区二区| 窝窝影院91人妻| 无遮挡黄片免费观看| 满18在线观看网站| 精品久久久久久久久久免费视频 | 午夜福利一区二区在线看| 欧美乱妇无乱码| 精品卡一卡二卡四卡免费| 日韩 欧美 亚洲 中文字幕| 一区二区三区精品91| 美女福利国产在线| 国产在线观看jvid| 18禁裸乳无遮挡免费网站照片 | 久久久久久免费高清国产稀缺| 国产精品九九99| 三级毛片av免费| 国产精品乱码一区二三区的特点 | 亚洲美女黄片视频| 欧美在线一区亚洲| 黄色a级毛片大全视频| 人人妻,人人澡人人爽秒播| 日韩av在线大香蕉| 亚洲精品久久午夜乱码| 欧美日本亚洲视频在线播放| 一个人观看的视频www高清免费观看 | 亚洲精华国产精华精| 中文字幕av电影在线播放| 国产欧美日韩一区二区精品| 久久影院123| 999久久久国产精品视频| 99精国产麻豆久久婷婷| 日韩欧美在线二视频| 老司机靠b影院| 日本a在线网址| 午夜免费鲁丝| 国产真人三级小视频在线观看| 国产精华一区二区三区| 欧美日韩国产mv在线观看视频| 午夜福利一区二区在线看| bbb黄色大片| 777久久人妻少妇嫩草av网站| 国产精品日韩av在线免费观看 | 在线观看免费午夜福利视频| 亚洲国产精品合色在线| 亚洲精品美女久久av网站| 99久久综合精品五月天人人| 妹子高潮喷水视频| 亚洲国产精品一区二区三区在线| 久久国产精品人妻蜜桃| 国产野战对白在线观看| 这个男人来自地球电影免费观看| 亚洲午夜理论影院| 国产精品国产av在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲人成网站在线播放欧美日韩| 国产成人精品在线电影| 亚洲精品国产色婷婷电影| 亚洲少妇的诱惑av| 高清毛片免费观看视频网站 | 男人舔女人的私密视频| 免费搜索国产男女视频| 真人做人爱边吃奶动态| 欧美午夜高清在线| 亚洲欧美激情综合另类| 侵犯人妻中文字幕一二三四区| 国产精品日韩av在线免费观看 | 婷婷丁香在线五月| 身体一侧抽搐| 人成视频在线观看免费观看| 91精品国产国语对白视频| 午夜福利一区二区在线看| 一进一出抽搐动态| 在线观看免费视频日本深夜| 国内毛片毛片毛片毛片毛片| 我的亚洲天堂| 一区二区日韩欧美中文字幕| 日本vs欧美在线观看视频| 亚洲精品在线美女| 69av精品久久久久久| 久久久久久久久中文| 久久中文看片网| 国产高清videossex| 1024香蕉在线观看| 久久午夜综合久久蜜桃| 巨乳人妻的诱惑在线观看| 日本一区二区免费在线视频| 麻豆成人av在线观看| 99国产综合亚洲精品| xxxhd国产人妻xxx| 国产欧美日韩一区二区三| 久久精品国产99精品国产亚洲性色 | 国内毛片毛片毛片毛片毛片| 久久久久久久久免费视频了| 一区福利在线观看| 国产97色在线日韩免费| ponron亚洲| 一区二区三区精品91| 欧美精品啪啪一区二区三区| 色综合站精品国产| 精品国产乱子伦一区二区三区| 国产高清激情床上av| 亚洲欧美激情综合另类| 久久婷婷成人综合色麻豆| 一边摸一边抽搐一进一小说| 欧美性长视频在线观看| 90打野战视频偷拍视频| 91老司机精品| 国产蜜桃级精品一区二区三区| 久久久国产成人精品二区 | 黄片大片在线免费观看| 老熟妇乱子伦视频在线观看| 国内毛片毛片毛片毛片毛片| 久久国产亚洲av麻豆专区| 亚洲精品久久成人aⅴ小说| 十分钟在线观看高清视频www| 日本免费一区二区三区高清不卡 | 夜夜爽天天搞| 老司机亚洲免费影院| 国产精品爽爽va在线观看网站 | 精品一区二区三区视频在线观看免费 | 国产成人精品在线电影| 黄色怎么调成土黄色| 亚洲精品一二三| 日韩欧美在线二视频| 国产又爽黄色视频| 国产97色在线日韩免费| 国产深夜福利视频在线观看| 中文欧美无线码| 亚洲欧美日韩另类电影网站| 国产成人精品久久二区二区91| 国产成人精品久久二区二区免费| 人人妻人人爽人人添夜夜欢视频| 在线播放国产精品三级| av欧美777| 成人国语在线视频| 日本五十路高清| 美女 人体艺术 gogo| 午夜视频精品福利| 午夜精品在线福利| 黄网站色视频无遮挡免费观看| 1024视频免费在线观看| 老熟妇乱子伦视频在线观看| 51午夜福利影视在线观看| 成人永久免费在线观看视频| 日韩视频一区二区在线观看| 精品高清国产在线一区| 欧美激情极品国产一区二区三区| 亚洲 欧美一区二区三区| 深夜精品福利| 91老司机精品| 人人妻人人爽人人添夜夜欢视频| 级片在线观看| 欧美午夜高清在线| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产亚洲av香蕉五月| 国产亚洲欧美精品永久| 一a级毛片在线观看| 久久久国产一区二区| 美女午夜性视频免费| 两人在一起打扑克的视频| 国产亚洲精品第一综合不卡| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产成人精品二区 | 久久久国产成人免费| 激情视频va一区二区三区| www.www免费av| 一级黄色大片毛片| 麻豆av在线久日| 黄色成人免费大全| 亚洲五月天丁香| 国产视频一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 美女国产高潮福利片在线看| 亚洲av熟女| a级毛片黄视频| 日韩av在线大香蕉| 亚洲av五月六月丁香网| 最近最新中文字幕大全免费视频| 久久精品亚洲熟妇少妇任你| 51午夜福利影视在线观看| 国产成+人综合+亚洲专区| 叶爱在线成人免费视频播放| 91在线观看av| 国产av精品麻豆| 亚洲午夜理论影院| 久久中文字幕人妻熟女| 欧美乱码精品一区二区三区| 欧美日韩亚洲综合一区二区三区_|