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

    不同烤煙套作模式對烤煙根際土壤細菌群落特征及烤煙產(chǎn)質(zhì)量的影響

    2017-05-30 10:48:04韋俊楊煥文徐照麗鄧小鵬王戈白羽祥呂世保譚小兵
    南方農(nóng)業(yè)學報 2017年4期
    關(guān)鍵詞:產(chǎn)質(zhì)量高通量測序套作

    韋俊 楊煥文 徐照麗 鄧小鵬 王戈 白羽祥 呂世?!∽T小兵

    摘要:【目的】探討云南不同烤煙套作模式對烤煙根際土壤細菌群落特征及烤煙產(chǎn)質(zhì)量的影響,為優(yōu)化當?shù)乜緹煼N植模式提供科學依據(jù)?!痉椒ā坎捎脝我蛩仉S機區(qū)組試驗,設(shè)烤煙套種白菜(T1)、烤煙套種結(jié)球甘藍(T2)和烤煙套種豌豆(T3)3個處理,以烤煙單作為對照(CK),采用16S rRNA測序技術(shù)分析不同套作模式烤煙根際土壤細菌群落結(jié)構(gòu)及多樣性差異,同時進行土壤優(yōu)勢細菌與土壤理化性狀的相關(guān)性分析及不同處理的烤煙產(chǎn)質(zhì)量分析?!窘Y(jié)果】16S rRNA測序結(jié)果表明,套作處理整體可提高烤煙根際土壤細菌多樣性指數(shù),Shannon指數(shù)排序為T3>T2>T1>CK,Simpson指數(shù)排序為T3=T2>T1>CK,Chao1和ACE指數(shù)排序為T3>T2>CK>T1。各處理烤煙根際土壤細菌主要分布在10個菌門,其中變形菌門(Proteobacteria)的豐度最高,占各樣本總有效序列的50%以上。CK處理根際土壤中副球菌屬(Paracocccus)、假單胞菌屬(Pseudomonas)、馬賽菌屬(Massilia)、Psychrobacter屬、短波單胞菌屬(Brevundimonas)和產(chǎn)卟啉桿菌屬(Porphyrobacter)豐度較套作處理顯著提高(P<0.05,下同),而鞘氨醇單胞菌屬(Sphingomonas)和Haliangium屬豐度顯著降低;烤煙套作處理顯著降低了烤煙根際土壤pH,同時提高了烤煙根際土壤的主要養(yǎng)分含量。相關(guān)性分析結(jié)果表明,CK處理土壤中豐度較高的假單胞菌屬(Pseudomonas)與土壤總有機碳含量呈顯著負相關(guān),豐度較低的鞘氨醇單胞菌屬(Sphingomonas)和Haliangium屬與土壤總有機碳含量呈顯著正相關(guān)。不同套作處理可不同程度地提高烤煙產(chǎn)質(zhì)量,產(chǎn)值、產(chǎn)量、均價和上等煙比例排序均為T3>T2>T1>CK,以T3處理與優(yōu)質(zhì)煙葉標準最接近,化學成分最協(xié)調(diào)?!窘Y(jié)論】烤煙套種白菜、結(jié)球甘藍和豌豆較烤煙單作可明顯提高烤煙根際細菌多樣性,不同程度地降低土壤中潛在致病菌豐度,提高烤煙產(chǎn)質(zhì)量,其中以套種豌豆的效果最佳。

    關(guān)鍵詞: 烤煙;套作;高通量測序;細菌;多樣性;產(chǎn)質(zhì)量

    中圖分類號: S572;S154.3 文獻標志碼:A 文章編號:2095-1191(2017)04-0601-08

    Abstract:【Objective】Effects of different intercropping patterns on rhizosphere soil bacteria colony characteristics and flue-cured tobacco yield and quality were studied in order to provide reference for improving local tobacco planting pattern. 【Method】Single-factor randomized block experiment was adopted and tobacco-Chinese cabbage intercropping(T1), toba-

    cco-common head cabbage intercropping(T2), tobacco-pea intercropping(T3) were set, and tobacco single planting was taken as control(CK). Rhizosphere bacteria colony characteristics and diversity index were analyzed by 16S rRNA sequen-

    cing. Correlation between dominant bacteria and soil physical and chemical properties, yield and quality of flue-cured tobacco in different treatments were also studied. 【Result】16S rRNA sequencing showed that intercropping could increase diversity indexes of rhizosphere soil bacteria. Shannon index ranking was as follows:T3>T2>T1>CK, Simpson index ranking was: T3=T2>T1>CK, Chaol index and ACE index ranking were the same: T3>T2>CK>T1. Rhizosphere soil bacteria in different treatments fell into ten phyla, and abundance of Proteobacteria was the highest and consisted of over 50% total effective sequences. In CK, abundance of Paracocccus,Pseudomonas,Massilia,Psychrobacter,Brevundimonas and Porphyrobacter in rhizosphere soil significantly increased compared with those in intercropping treatments(P<0.05, the same below), but abundance of Sphingomonas and Haliangium significantly decreased. Intercropping treatments significantly decreased pH value but increased major nutrient contents of rhizosphere soil. Correlation analysis showed that:in CK, there was significant negative correlation between high-abundance Pseudomonas and total soil organic carbon content, and significant positive correlation between low-abundance Sphingomonas and Haliangium and soil total organic carbon content. Different intercropping treatments could increase tobacco quality to various extents, and the ranks of output value, yield, average price and proportion of first-class tobacco were the same: T3>T2>T1>CK. The tobacco quality in T3 was the closest to high-quality tobacco standard with the most balanced chemical components. 【Conclusion】Tobacco intercropping Chinese cabbage, common head cabbage and pea can obviously increase diversity of rhizosphere soil bacteria, lower abundance of potential pathogenic bacteria in soil to certain extent, and improve tobacco quality. Tobacco-pea intercropping has the optimal effects.

    Key words: flue-cured tobacco; intercropping; high throughput sequencing; bacteria; diversity; yield and quality

    0 引言

    【研究意義】微生物作為土壤的重要組成成分與土壤營養(yǎng)轉(zhuǎn)化、土壤肥力及土壤質(zhì)量密切相關(guān)(陽顯斌等,2016)。土壤微生物群落結(jié)構(gòu)和多樣性在土壤有機質(zhì)分解、養(yǎng)分循環(huán)利用及維持生態(tài)系統(tǒng)的穩(wěn)定性方 面發(fā)揮著重要作用(賈志紅等,2004)。烤煙是云南重要的經(jīng)濟作物之一,但隨著種煙面積的逐年增加,近年來糧煙爭地矛盾日益凸顯,加之煙區(qū)連作現(xiàn)象普遍,給植煙土壤的健康及煙葉可持續(xù)生產(chǎn)帶來了巨大威脅。間套種栽培模式在我國有著悠久歷史(吳存浩,1996),其中適宜的套作模式對緩解用地矛盾、增加土地產(chǎn)出、降低病蟲害、維持土壤生態(tài)系統(tǒng)穩(wěn)定性及促進土地可持續(xù)利用均具有重要作用。煙/白菜、煙/結(jié)球甘藍和煙/豌豆是云南煙區(qū)烤煙生產(chǎn)中常見的套作搭配模式。因此,研究這3種套作模式對煙草根際土壤微生物數(shù)量和多樣性的影響,進一步探明不同套作模式烤煙根際微生物結(jié)構(gòu)差異及其與土壤理化性質(zhì)的相關(guān)性,可為今后優(yōu)化烤煙種植模式、促進烤煙產(chǎn)質(zhì)量提高提供參考。【前人研究進展】近年來,套作模式的研究正在從宏觀的作物產(chǎn)質(zhì)量和病蟲害防治向微觀的土壤微生態(tài)和土壤成分互作方面延伸。套作植物間不僅競爭著地上、地下部的生長空間,還通過根系向土壤釋放次生代謝物質(zhì)來影響土壤微生物群落結(jié)構(gòu),從而間接影響作物的生長發(fā)育和養(yǎng)分吸收。套作可以提高土壤微生物群落多樣性,微生物群落多樣性又與作物生長發(fā)育息息相關(guān)(Zuo,2004)。李秋紅等(2010)研究了黃瓜套作毛蔥和蒜對其根際土壤細菌群落結(jié)構(gòu)多樣性的影響,結(jié)果表明,套作改變了黃瓜根際土壤細菌群落結(jié)構(gòu)的多樣性和優(yōu)勢菌群結(jié)構(gòu),并且可提高黃瓜產(chǎn)量。徐強等(2013)研究發(fā)現(xiàn),線辣椒套種玉米可顯著提高線辣椒根際土壤微生物群落的種間相遇幾率和Shannon、McIntosh多樣性指數(shù),且線辣椒的根際土壤微生物多樣性指數(shù)分別與其生物學產(chǎn)量呈顯著或極顯著正相關(guān)。涂勇等(2015)研究指出,烤煙—大豆帶狀套作可顯著增加作物根際土壤的細菌、放線菌數(shù)量,減少真菌數(shù)量,同時可提高與土壤氮素代謝相關(guān)的氨化細菌、硝酸細菌、亞硝酸細菌、自生固氮菌等功能性微生物數(shù)量,能降低烤煙主要土傳病害的發(fā)生。陽顯斌等(2016)研究指出,煙/蒜套作可提高烤煙根際土壤細菌和放線菌數(shù)量,降低土壤真菌數(shù)量,在烤煙生長發(fā)育的中后期顯著提高烤煙根際土壤解鉀菌和解磷菌含量,從而提高根際可吸收養(yǎng)分,有效改善土壤根際微生態(tài)。張黎杰等(2016)研究表明,西瓜—草菇—辣椒輪作模式可增加土壤中細菌和放線菌數(shù)量,同時降低有害病原真菌數(shù)量,改善土壤理化性質(zhì)和恢復土壤微生物種群數(shù)量,有利于消除土壤連作障礙的影響?!颈狙芯壳腥朦c】目前,關(guān)于烤煙套作白菜、豌豆和結(jié)球甘藍對烤煙根際細菌微生物群落的影響及其與土壤理化性質(zhì)相關(guān)性的研究尚無報道?!緮M解決的關(guān)鍵問題】基于Illumina MiSeq測序平臺(Luo et al.,2012),通過構(gòu)建16S rDNA基因文庫(Youssef et al.,2009;Caporaso et al.,2011;Hess et al.,2011),從細菌豐度差異和群落多樣性方面分析烤煙套作白菜、豌豆和結(jié)球甘藍的烤煙根際微生物群落優(yōu)劣及其與土壤理化性質(zhì)的關(guān)系,旨在揭示這3種烤煙套作模式對烤煙根際土壤質(zhì)量變化影響的內(nèi)部生態(tài)學機理,為優(yōu)化烤煙種植模式提供科學依據(jù)。

    1 材料與方法

    1. 1 試驗材料

    試驗用烤煙品種為當?shù)刂髟云贩NK326,采用漂浮育苗方式育苗,烤煙種子由玉溪中煙種子有限責任公司提供;白菜、結(jié)球甘藍和豌豆種子由中農(nóng)種子中心提供。

    1. 2 試驗方法

    大田試驗于2015年在云南省楚雄州姚安縣進行。采用單因素隨機區(qū)組試驗,設(shè)烤煙套種白菜(T1)、烤煙套種結(jié)球甘藍(T2)和烤煙套種豌豆(T3)3個處理,以烤煙單作為對照(CK),每處理重復3次。煙苗于4月20日統(tǒng)一移栽,以純氮120 kg/ha為施肥標準,種植密度16500株/ha,參考優(yōu)質(zhì)烤煙栽培技術(shù)進行田間管理。套作作物白菜、結(jié)球甘藍、豌豆于8月5日(烤煙生長的中后期)在煙墑中下部移栽,10月20日開始采收,12月20日采收完畢,田間管理按作物優(yōu)產(chǎn)田栽培技術(shù)統(tǒng)一進行。

    1. 3 測定項目及方法

    1. 3. 1 樣品采集

    1. 3. 1. 1 煙葉樣品 采收各處理上部成熟煙葉進行烘烤,烤后分等級稱重,選取上部葉片(等級均為上桔二)粉碎過篩常溫保存,用于煙葉化學成分分析。

    1. 3. 1. 2 土壤樣品 烤煙和套作作物收獲后用螺旋取土鉆垂直鉆取烤煙根際0~20 cm土壤,將土壤分為兩份,一份避光風干至恒重后用于土壤理化性質(zhì)分析,另一份置于-80 ℃冰箱存儲,用于細菌16S rRNA測序。

    1. 3. 2 土壤細菌16S rRNA測序 應(yīng)用CTAB法提取土壤樣本細菌基因組DNA,用2%瓊脂糖凝膠電泳檢測DNA濃度與純度,將樣品放入離心管中,無菌水稀釋樣品至樣品基因組DNA濃度為1 ng/μl。以提取的DNA為模板,根據(jù)細菌測序區(qū)域的選擇515F(16S V4區(qū)引物)鑒定細菌多樣性。使用帶Barcode的特異性引物(New England Biolabs公司的Phusion High-Fidelity PCR Master Mix with GC Buffer)和高效保真酶進行PCR。應(yīng)用瓊脂糖凝膠電泳檢測PCR產(chǎn)物DNA的完整性、純度和片段大小,依據(jù)PCR產(chǎn)物濃度等量混合樣品,充分混勻后使用瓊脂糖凝膠電泳檢測PCR產(chǎn)物,使用Qiagen公司提供的膠回收試劑盒對目的條帶產(chǎn)物進行回收。

    1. 3. 3 土壤細菌基因組DNA信息分析

    1. 3. 3. 1 測序數(shù)據(jù)處理 首先從Illumina MiSeq測序平臺下機數(shù)據(jù)中分離出樣品數(shù)據(jù),使用FLASH(Mago and Salzberg,2011)拼接截去Barcode和引物序列后的樣品Reads,得到Raw Tags,Raw Tags經(jīng)過濾處理(Bokulich et al.,2012)后得到Clean Tags。然后按照Tags質(zhì)量控制流程(Caporaso et al.,2010)通過數(shù)據(jù)庫Gold Database與UCHIME Algorithm(Edgar et al.,2011)比對后將Clean Tags進行檢測和去除嵌合體序列(Haas et al.,2011)處理,最終得到Effective Tags。

    1. 3. 3. 2 OTU聚類和物種注釋 利用UPARSE對各處理土壤細菌樣品的Effective Tags進行聚類分析處理,以97%的一致性將序列聚類成為OTUs,選取OTUs的代表性序列進行物種注釋,用Mothur方法結(jié)合SILVA(Altschul et al.,1990)SSUrRNA數(shù)據(jù)庫(Kljalg et al.,2013)進行物種注釋分析(設(shè)定閾值為0.8~1.0),獲得分類學信息并分別在界、門、綱、目、科、屬、種各個分類水平統(tǒng)計樣本的群落組成。使用PyNAST(Yilmaz et al.,2014)結(jié)合Green Gene數(shù)據(jù)庫中的Core Set數(shù)據(jù)信息進行快速多序列比對,得到所有OTUs代表序列的系統(tǒng)發(fā)生關(guān)系。最后以數(shù)據(jù)量最少的樣本為標準進行均一化處理。

    1. 3. 3. 3 樣品多樣性分析 使用Qiime 1.7.0計算各處理樣本的Observed-species、Goods-coverage、Shannon、Simpson、Chao1和ACE指數(shù),其中Goods-coverage指數(shù)用于分析測序深度,Shannon和Simpson指數(shù)用于分析菌群多樣性,Chao和ACE指數(shù)用于分析菌群豐富度。

    1. 3. 4 土壤理化性質(zhì)測定 參照Liu(1996)測定土壤pH及有機質(zhì)、全氮、全磷、全鉀、速效氮、速效磷和速效鉀含量。

    1. 3. 5 煙葉化學成分測定 依照YC/T 159162- 2002采用連續(xù)流動法測定煙葉的總糖、還原糖、總氮、煙堿和氯含量,依照YC/T 217-2007應(yīng)用連續(xù)流動法測定煙葉鉀含量。

    1. 4 統(tǒng)計分析

    應(yīng)用Exce1 2007和SPSS 19.0進行數(shù)據(jù)處理和統(tǒng)計分析;利用UPARSE和Mothur方法結(jié)合SILVASSUrRNA數(shù)據(jù)庫進行物種注釋分析;使用PyNAST結(jié)合GreenGene數(shù)據(jù)庫中的Core Set數(shù)據(jù)信息進行快速多序列比對。

    2 結(jié)果與分析

    2. 1 不同套作模式對烤煙根際土壤細菌α多樣性指數(shù)的影響

    由表1可知,烤煙根區(qū)的細菌OTUs數(shù)為2759.33~ 3246.00個,其中T3處理最高,與其他處理間差異顯著(P<0.05,下同),T1、T2和CK處理間OTUs數(shù)差異未達顯著水平(P>0.05,下同)。各處理樣本覆蓋率較高,為96.2%~97.1%,說明構(gòu)建的文庫包含樣品中大部分細菌種類,能較好地反映土壤細菌群落結(jié)構(gòu)。不同套作模式對土壤α多樣性指數(shù)的影響存在差異,Shannon指數(shù)排序為T3>T2>T1>CK,其中T3與T2處理差異不顯著,但顯著高于T1和CK處理;Simpson指數(shù)大小排序為T3=T2>T1>CK,3個套作處理差異不顯著,但均顯著高于CK,說明烤煙套作處理較烤煙單作可提高其根際細菌的多樣性指數(shù),其中套種豌豆的提高幅度最大。Chao1和ACE指數(shù)排序均為T3>T2>CK>T1,說明套種結(jié)球甘藍和豌豆有利于提高烤煙根際土壤細菌豐富度指數(shù),套種白菜則降低了烤煙根際土壤細菌豐富度指數(shù),但套種結(jié)球甘藍和白菜與烤煙單作間細菌豐富度指數(shù)差異不顯著。

    2. 2 不同套作模式對烤煙根際土壤細菌豐度的影響

    由圖1可看出,各處理烤煙根際土壤細菌主要分布在10個菌門,包括變形菌門(Proteobacteria)、擬桿菌門(Bacteroidetes)、放線菌門(Actinobacteria)、酸桿菌門(Acidobacteria)、厚壁菌門(Firmicutes)、芽單胞菌門(Gemmatimonadetes)、藍藻門(Cyanobacteria)、Saccharibacteria門、疣微門(Verrucomicrobia)和浮霉菌門(Planctomycetes);其中變形菌門的豐度最高,占各樣本總有效序列的50%以上。

    對樣本豐度前32的菌屬進行方差分析,結(jié)果(表2)發(fā)現(xiàn),不同套作處理對烤煙根際土壤細菌豐度影響差異明顯。CK的副球菌屬(Paracocccus)、馬賽菌屬(Massilia)、Psychrobacter屬、短波單胞菌屬(Brevundimonas)和產(chǎn)卟啉桿菌屬(Porphyrobacter)豐度顯著高于3個套作處理;假單胞菌屬(Pseudomonas)豐度最高,顯著高于T3處理;鞘氨醇單胞菌屬(Sphingomonas)和Haliangium屬顯著低于3個套作處理;鞘氨醇桿菌屬(Sphingobacterium)豐度與其他處理差異不顯著;芽孢桿菌屬(Bacillus)豐度低于3個套作處理,與T1處理差異顯著;新鞘氨醇桿菌屬(Novosphingobium)和溶桿菌屬(Lysobacter)豐度均低于T1處理,但與T2和T3處理差異不顯著。

    2. 3 不同套作處理對烤煙根際土壤理化性質(zhì)的影響

    由表3可知,烤煙套作處理顯著降低了烤煙根際土壤pH,說明套種3種作物有效改善了土壤pH,使其更接近于種植烤煙的適宜pH范圍(5.5~6.5),其中以T3處理的改善效果最好;同時,烤煙套種處理較CK可提高烤煙根區(qū)土壤的總有機碳、全氮、全磷、全鉀、速效氮、速效磷和速效鉀含量,其中3個套種處理的全氮、速效氮和速效鉀含量均顯著高于CK。

    2. 4 不同套作處理烤煙根際土壤細菌與土壤理化性質(zhì)的相關(guān)性分析結(jié)果

    將土壤理化性質(zhì)與各樣本豐度前32的菌屬進行相關(guān)性分析,結(jié)果(表4)發(fā)現(xiàn),土壤pH與寡養(yǎng)單胞菌屬(Stenotrophomonas)呈極顯著正相關(guān)(P<0.01,下同);總有機碳含量與鞘氨醇單胞菌屬(Sphingomonas)和Haliangium屬呈顯著正相關(guān),與假單胞菌屬(Pseudomonas)呈顯著負相關(guān);全氮含量與馬賽菌屬(Massilia)、Psychrobacter和Dyadobacter屬呈顯著負相關(guān);全鉀含量與Steroidobacter屬和亞硝化螺菌屬(Nitrosospira)呈顯著正相關(guān),與溶桿菌屬(Lysobacter)呈極顯著正相關(guān);速效氮含量與Haliangium屬呈顯著正相關(guān);速效鉀含量與鞘氨醇單胞菌屬(Sphingomonas)呈顯著正相關(guān),與馬賽菌屬(Massilia)、金黃桿菌屬(Chryseobacterium)和產(chǎn)卟啉桿菌屬(Porphyrobacter)呈顯著負相關(guān)。

    2. 5 不同套作模式對烤煙產(chǎn)質(zhì)量的影響

    2. 5. 1 經(jīng)濟性狀 由表5可知,不同套作模式較單作可不同程度地提高烤煙產(chǎn)質(zhì)量,產(chǎn)值、產(chǎn)量、均價和上等煙比例排序均為T3>T2>T1>CK,其中3個套作處理的產(chǎn)值和均價均顯著高于CK,T3處理的產(chǎn)量和上等煙比例顯著高于CK。

    2. 5. 2 化學成分 中高海拔煙區(qū)優(yōu)質(zhì)煙葉以總氮含量1.5%~3.5%、煙堿1.5%~3.5%、總糖含量20.0%~28.0%、還原糖含量18.0%~25.0%、鉀含量在2.0%以上、氯含量在1.0%以下、糖堿比10左右、氮堿比1左右、鉀氯比4~10為宜(劉國順,2003;胡榮海,2007)。由表6可知,T1和T3處理的總糖含量、CK和T3處理的還原糖含量符合優(yōu)質(zhì)煙葉標準,說明T3處理煙葉糖含量較適宜;各處理總氮、煙堿和氯含量均處于適宜范圍內(nèi),鉀含量均偏低,其中,T3處理的總氮、煙堿和鉀含量最高,氯含量最低。各處理糖堿比均高于優(yōu)質(zhì)煙葉標準;氮堿比均小于1;各處理兩糖差值均高于優(yōu)質(zhì)煙葉標準(2%~5%);除T處理外,鉀氯比均在適宜范圍內(nèi)。綜合分析各項化學成分含量及相關(guān)比值,以T3處理與優(yōu)質(zhì)煙葉標準最接近,化學成分最協(xié)調(diào)。

    3 討論

    3. 1 套作對烤煙根際土壤細菌群落及土壤理化性質(zhì)的影響

    套作處理可顯著影響土壤微生物的數(shù)量及多樣性。馮俊喜等(2011)研究發(fā)現(xiàn),小麥/烤煙套作處理較烤煙單作可有效提高各取樣時期的土壤微生物多樣性指數(shù)。徐強等(2013)研究表明,線辣椒與玉米套作可提高線辣椒根際土壤微生物多樣性指數(shù),從而增強土壤微生物對單一碳源的利用能力。本研究發(fā)現(xiàn),烤煙套種白菜、結(jié)球甘藍和碗豆后,烤煙根際土壤細菌多樣性指數(shù)均表現(xiàn)為T3>T2>T1>CK,且T3處理的Shannon指數(shù)及各套作處理的Simpson指數(shù)顯著高于CK。說明烤煙套種處理較烤煙單作可不同程度地提高烤煙根際細菌多樣性指數(shù),與馮俊喜等(2011)和徐強等(2013)的研究結(jié)果相似。不同的套作搭配品種對微生物多樣性的影響效果存在差異,本研究發(fā)現(xiàn)套種豌豆對烤煙根際細菌多樣性指數(shù)的提高幅度最大,其原因可能是烤煙與豌豆套作,烤煙吸收土壤中的氮元素從而降低豌豆可吸收養(yǎng)分,氮脅迫條件下促使豌豆形成根瘤,合成氮素,固氮作用有效提高了土壤養(yǎng)分含量,進一步優(yōu)化土壤質(zhì)量,從而提高烤煙根際細菌多樣性。

    烤煙栽培模式不同程度地影響著土壤微生物數(shù)量,進而影響烤煙的養(yǎng)分吸收和病蟲害的發(fā)生。如煙/蒜套作可提高根際土壤細菌含量,降低真菌數(shù)量,在烤煙生長的中后期顯著提高根際土壤解鉀菌和解磷菌含量(陽顯斌等,2016),而烤煙單作使益生菌豐度降低,潛在的致病菌增加(石秋環(huán)等,2009)。本研究發(fā)現(xiàn),相對于套作處理,烤煙單作的芽孢桿菌屬豐度最低,假單胞菌豐度最高。芽孢桿菌可通過同化作用產(chǎn)生抑菌物質(zhì),進而有效抑制有害病原物的生長,具有很強的抗逆性能和抗菌作用(Elliott et al.,2001),而假單胞菌是能夠寄生于植物的潛在病原菌(Fahlgren et al.,2011)。因此,烤煙套種3種作物不同程度地提高了益生菌的數(shù)量,降低了潛在病原菌的數(shù)量。然而,同種菌屬中不同菌種在功能上并非完全一致,且同種菌種含量不同也可能導致作用效果不同,至于是何種細菌在一定豐度范圍內(nèi)影響作物生長有待進一步研究。

    間套作模式可顯著影響土壤的養(yǎng)分含量。如禾本科與豆科作物間作可提高豆科作物固氮能力,有效增加土壤的氮素營養(yǎng)(Cong et al.,2015);花生/玉米(Inal et al.,2007)、青蔥/黃瓜(Xiao et al.,2013)間作可顯著提高土壤中的速效鉀含量與植株中的鉀含量;煙/蒜套作可提高成熟期烤煙下部葉和中部葉根際土壤的有效磷含量(唐彪等,2015)。本研究發(fā)現(xiàn),烤煙套種白菜、結(jié)球甘藍和豌豆均可降低烤煙根際土壤pH,使土壤酸堿度趨近適宜范圍,同時提高烤煙根際的全氮、全磷、全鉀、速效氮、速效磷、速效鉀及總有機碳的含量。經(jīng)相關(guān)分析可知,土壤pH和主要養(yǎng)分含量與不同菌屬豐度的相關(guān)性差異明顯,其中潛在致病菌假單胞菌屬與總有機碳含量呈顯著負相關(guān),而單作處理的假單胞菌屬豐度均高于套作處理,可能是套作處理尤其是套種豌豆處理較烤煙單作的總養(yǎng)分含量較高的原因之一。

    3. 2 套作對烤煙產(chǎn)質(zhì)量的影響

    套作處理可影響烤煙的產(chǎn)質(zhì)量。付利波等(2005)研究表明,煙田套種白菜和馬鈴薯可使煙葉產(chǎn)值比烤煙單作處理提高13.73%~33.20%。時安東等(2011)研究表明,烤煙套種紅薯、花生可使煙葉產(chǎn)量、產(chǎn)值連續(xù)6年顯著高于烤煙單作,同時煙葉化學成分含量適中,比例更協(xié)調(diào)。本研究發(fā)現(xiàn),不同套作模式較單作可不同程度地提高烤煙的產(chǎn)質(zhì)量,產(chǎn)值、產(chǎn)量、均價和上等煙比例排序均為T3>T2>T1>CK,其中3個套作處理的產(chǎn)值和均價均顯著高于CK,T3處理的產(chǎn)量和上等煙比例顯著高于CK;同時綜合分析各項化學成分含量及相關(guān)比值,以T3處理與優(yōu)質(zhì)煙葉標準最接近,化學成分最協(xié)調(diào)。說明套作處理可提高煙葉品質(zhì),尤其以套作豌豆對烤煙品質(zhì)的改善作用最明顯。

    4 結(jié)論

    本研究結(jié)果表明,烤煙套種白菜、結(jié)球甘藍和豌豆較烤煙單作可明顯提高烤煙根際細菌多樣性,不同程度地降低土壤中潛在致病菌豐度,提高烤煙產(chǎn)質(zhì)量,其中以套種豌豆的效果最佳。

    參考文獻:

    馮俊喜,王樹聲,石屹,張繼光,鄭林林,李湘?zhèn)?,胡慶輝,王浩浩. 2011. 山東煙區(qū)不同種植模式土壤微生物群落特征研究[J]. 中國煙草科學,32(2):38-42.

    Feng J X,Wang S S,Shi Y,Zhang J G,Zheng L L,Li X W,Hu Q H,Wang H H. 2011. Characteristics of soil microbial community under different cropping systems in flue-cured tobacco growing areas of Shandong province[J]. Chinese Journal of Tobacco Science,32(2):38-42.

    付利波,王毅,楊躍,瞿興,左天龍,耿明建,蘇帆,洪麗芳. 2005. 利用煙田套作調(diào)控高肥力土壤烤煙生產(chǎn)[J]. 植物營養(yǎng)與肥料學報,11(1):128-132.

    Fu L B,Wang Y,Yang Y,Qu X,Zuo T L,Geng M J,Su F,Hong L F. 2005. Regulating the growth of flue-cured tobacco by interplanting in fertility soil[J]. Plant Nutrition and Fertilizer Science,11(1):128-132.

    賈志紅,楊珍平,張永清,苗果園. 2004. 麥田土壤微生物三大種群數(shù)量的研究[J]. 麥類作物學報,24(3):53-56.

    Jia Z H,Yang Z P,Zhang Y Q,Miao G Y. 2004. Study on quantity of three main colony of soil microbe in wheat farmland[J]. Journal of Triticeae Crops,24(3):53-56.

    胡榮海. 2007. 云南煙草栽培學[M]. 北京:科學出版社.

    Hu R H. 2007. Yunnan Tobacco Cultivation[M]. Beijing:Science Press.

    李秋紅,吳鳳芝,景芳毅. 2010. 套作對黃瓜根際土壤細菌群落結(jié)構(gòu)的多樣性影響[J]. 土壤通報,41(1):47-52.

    Li Q H,Wu F Z,Jing F Y. 2010. T-RFLP analysis of effects of intercropping on in cucumber rhizosphere soil bacterial community structure diversity[J]. Chinese Journal of Soil Science,41(1):47-52.

    劉國順. 2003. 煙草栽培學[M]. 北京:中國農(nóng)業(yè)出版社.

    Liu G S. 2003. Tobacco Cultivation[M]. Beijing:China Agriculture Press.

    時安東,李建偉,袁玲. 2011. 輪間作系統(tǒng)對烤煙產(chǎn)量、品質(zhì)和土壤養(yǎng)分的影響[J]. 植物營養(yǎng)與肥料學報,17(2):411-418.

    Shi A D,Li J W,Yuan L. 2011. Effects of rotation and intercropping systems on yield,quality of flue-cured tobacco and soil nutrients[J]. Plant Nutrition and Fertilizer Science,17(2):411-418.

    石秋環(huán),焦楓,耿偉,李群平,時向東. 2009. 烤煙連作土壤環(huán)境中的障礙因子研究綜述[J]. 中國煙草學報,15(6):81-84.

    Shi Q H,Jiao F,Geng W,Li Q P,Shi X D. 2009. An overview on research into factors hindering continuous cropping in flue-cured tobacco[J]. Acta Tabacaria Sinica,15(6):81-84.

    唐彪,張錫洲,陽顯斌. 2015. 煙蒜輪作與套作對烤煙產(chǎn)量及根際土壤磷組分的影響[J]. 應(yīng)用生態(tài)學報,26(7):1977-1984.

    Tang B,Zhang X Z,Yang X B. 2015. Effects of tobacco garlic crop rotation and intercropping on tobacco yield and rhizosphere soil phosphorus fractions[J]. Chinese Journal of A-

    pplied Ecology,26(7):1977-1984.

    涂勇,楊文鈺,劉衛(wèi)國,雍太文,江連強,王小春. 2015. 大豆與烤煙不同套作年限對根際土壤微生物數(shù)量的影響[J]. 作物學報,41(5):733-742.

    Tu Y,Yang W Y,Liu W G,Yong T W,Jiang L Q,Wang X C. 2015. Effects of relay strip intercropping years between flue-cured tobacco and soybean on rhizospheric microbes quantities[J]. Acta Agronomica Sinica,41(5):733-742.

    吳存浩. 1996. 中國農(nóng)業(yè)史[M]. 北京:警官教育出版社.

    Wu C H. 1996. Agricultural History of China[M]. Beijing:Police Education Press.

    徐強,劉艷君,陶鴻. 2013. 間套作玉米對線辣椒根際土壤微生物生態(tài)特征的影響[J]. 中國生態(tài)農(nóng)業(yè)學報,21(9):1078-1087.

    Xu Q,Liu Y J,Tao H. 2013. Effects of relay intercropping maize on rhizosphere soil microbial ecological characteristics in capsicum fields[J]. Chinese Journal of Eco-Agriculture,21(9):1078-1087.

    陽顯斌,李廷軒,張錫洲,陳代榮,梁永江,張長華. 2016. 煙蒜輪作與套作對土壤微生物類群數(shù)量的影響[J]. 土壤,48(4):698-704.

    Yang X B,Li T X,Zhang X Z,Chen D R,Liang Y J,Zhang C H. 2016. Effects of tobacco garlic crop rotation and tobacco garlic crop intercropping on soil microbial groups in tobacco fields[J]. Soil,48(4):698-704.

    張黎杰,周玲玲,姜若勇,田福發(fā),劉金兵. 2016. 大棚西瓜—草菇—辣椒輪作對解除設(shè)施西瓜連作障礙的效果[J]. 江西農(nóng)業(yè)學報,28(4):11-14.

    Zhang L J,Zhou L L,Jiang R Y,Tian F F,Liu J B. 2016. Effect of watre melon-straw mushroom-pepper rotation on relie-

    ving continuous cropping obstacle of watermelon in greenhouse[J]. Acta Agriculturae Jiangxi,28(4):11-14.

    Altschul S F,Gish W,Miller W,Myers E W,Lipman D J. 1990. Basic local alignment search tool[J]. Journal of Molecular Biology,215(3):403-410.

    Bokulich N A,Subramanian S,F(xiàn)aith J J,Gevers D,Gordon J I,Knight R,Mill D A,Caporaso J G. 2012. Quality-filtering vastly improves diversity estimates from lllumina amplicon sequencing[J]. Nature Methods,10(1):57-59.

    Caporaso J G,Lauber C L,Walters W A,Berg-Lyons D,Lozupone C A,Turnbaugh P J,F(xiàn)ierer N,Knight R. 2011. Global pa-

    tterns of 16S rRNA diversity at a depth of millions of sequences per sample[J]. Proceedings of the National Academy of Sciences of the United States of America,108(S1):4516-4522.

    Caporaso J G,Kuczynski J,Stombaugh J,Bittinger K,Bushman F D,Costello E K, Fierer N,Pe■a A G,Goodrich J K,Gordon J I,Huttley G A,Kelley S T,Knights D,Koenig J E,Ley R E,Lozupone C A,McDonald D,Muegge B D,Pirrung M, Reeder J,Sevinsky J R,Turnbaugh P J,Walters W A,Widmann J,Yatsunenko T,Zaneveld J,Knight R. 2010. QIIME allows analysis of high-throughput community sequencing data[J]. Nature Methods,7(5):335-336.

    Cong W F,Hoffland E,Li L,Six J,Sun J H,Bao X G,Zhang F S,van der Werf W. 2015. Intercropping enhances soil carbon and nitrogen[J]. Global Change Biology,21(4):1715-1726.

    Edgar R C,Haas B J,Clemente J C,Quince C,Knight R. 2011. UCHIME improves sensitivity and speed of chimera detection[J]. Bioinformatics,27(16):2194-2200.

    Elliott M L,Des J E A,Batson W E Jr,Caceres J,Brannen P M,Howell C R,Benson D M,Conway K E,Rothrock C S,Schneider R W,Ownley B H,Canaday C H,Keinath A P,Huber D M,Sumner D R,Motsenbocker C E,Thaxton P M,Cubeta M A,Adams P D,Backman P A,F(xiàn)ajardo J,Newman M A,Pereira R M. 2001. Viability and stability of biological control agents on cotton and snap bean seeds[J]. Pest Mana-

    gement Science,57(8):695-706.

    Fahlgren C,Bratbak G,Sandaa R A,Thyrhaug R,Zweifel U L. 2011. Diversity of airborne bacteria in samples collected using different devices for aerosol collection[J]. Aerobiologia,27(2):107-120.

    Haas B J,Gevers D,Earl A M,F(xiàn)eldgarden M,Ward D V,Giannoukos G,Ciulla D,Tabbaa D,Highlander S K,Sodergren E,Methé B,DeSantis T Z,the Human Microbiome Consortium,Petrosino J F,Knight R,Birren B W. 2011. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons[J]. Genome Research,21(3):494-504.

    Hess M,Sczyrba A,Egan R,Kim T W,Chokhawala H,Schroth G,Luo S,Clark D S,Chen F,Zhang T,Mackie R I,Pennacchio L A,Tringe S G,Visel A,Woyke T,Wang Z,Rubin E M. 2011. Metagenomic discovery of biomass-degrading genes and genomes from cow rumen[J]. Science,331(6016):463-467.

    Inal A, Gunes A, Zhang F, Cakmak I. 2007. Peanut/maize intercropping induced changes in rhizosphere and nutrient concentrations in shoots[J]. Plant Physiology & Bioche-

    mistry,45(5):350-356.

    K■ljalg U,Nilsson R H,Abarenkov K,Tedersoo L,Taylor A F,Bahram M,Bates S T,Bruns T D,Bengtsson-Palme J,Callaghan T M,Douglas B,Drenkhan T,Eberhardt U,Due■as M, Grebenc T,Griffith G W,Hartmann M,Kirk P M,Kohout P, Larsson E,Lindahl B D,Lücking R,Martín M P,Matheny P B, Nguyen N H,Niskanen T,Oja J,Peay K G,Peintner U, Peterson M,P■ldmaa K,Saag L,Saar I,Schü■ler A,Scott J A,Senés C,Smith M E,Suija A,Taylor D L,Telleria M T,Weiss M,Larsson K H. 2013. Towards a unified paradigm for sequence-based identification of fungi[J]. Molecular Ecology,22(21):5271-5277.

    Liu G S. 1996. Standard Method of Observation and Analysis of the Chinese Ecosystem Research Network-Analysis of the Soil Physicochemical Characteristics and Delineation of the Soil Profiles[M]. Beijing:Standard Press.

    Luo C,Tsementzi D,Kyrpides N,Read T,Konstantinidis K T. 2012. Direct comparisons of lllumina vs. Roche 454 sequen-

    cing technologies on the same microbial community DNA sample[J]. Plos One,7(2):e30087.

    Mago■ T,Salzberg S L. 2011. FLASH:Fast length adjustment of short reads to improve genome assemblies[J]. Bioinformatis,27(21):2957-2963.

    Xiao X,Cheng Z,Meng H,Liu L,Li H,Dong Y. 2013. Intercropping of green garlic(Allium sativum L.)induces nutrient concentration changes in the soil and plants in continuously cropped cucumber(Cucumis sativus L.)in a plastic tunnel[J]. Plos One,8(4):e62173.

    Yilmaz P,Parfrey L W,Yarza P,Gerken J,Pruesse E,Quast C,Schweer T, Peplies J, Ludwig W, Gl■ckner F O. 2014. The SILVA and “All-species Living Tree Project(LTP)” taxonomic frameworks[J]. Nucleic Acids Research,42(Database issue):643-648.

    Youssef N,Sheik C S,Krumholz L R,Najar F Z,Roe B A,Elshahed M S. 2009. Comparison of species richness estimates obtained using nearly complete fragments and simulated pyrosequencing-generated fragments in 16S rRNA gene-based environmental surveys[J]. Applied & Environmental Microbiology,75(16):5227-5236.

    Zuo Y M. 2004. The mechanisms of root exudates of maize in improvement of iron nutrition of peanut in peanut/maize intercropping system by 14C tracer technique[J]. Acta Agriculturae Nucleatae Sinica,18(1):43-46.

    (責任編輯 王 暉)

    猜你喜歡
    產(chǎn)質(zhì)量高通量測序套作
    桑園全年免耕周年套作技術(shù)進行專家測產(chǎn)
    遠離套作之嫌,力創(chuàng)真我之文—以2019年浙江高考考場作文為例
    川明參輪作對煙地土壤微生物群落結(jié)構(gòu)的影響
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    人參根際真菌群落多樣性及組成的變化
    LncRNAs作為miRNA的靶模擬物調(diào)節(jié)miRNA
    磷素配比及一體肥對烤煙產(chǎn)質(zhì)量的影響研究
    南陽煙區(qū)密集烘烤不同裝煙方式對烤后煙葉產(chǎn)質(zhì)量的影響
    套作條件下不同品種小麥產(chǎn)量及其構(gòu)成分析
    移栽期及覆蓋物對煙草生長發(fā)育及產(chǎn)質(zhì)量的影響
    好男人在线观看高清免费视频| 色尼玛亚洲综合影院| 久久久久久国产a免费观看| 亚洲美女黄片视频| 两性午夜刺激爽爽歪歪视频在线观看| 真人做人爱边吃奶动态| 人妻丰满熟妇av一区二区三区| 日韩欧美国产在线观看| 日韩高清综合在线| 久久中文看片网| 动漫黄色视频在线观看| 熟妇人妻久久中文字幕3abv| 怎么达到女性高潮| 欧美丝袜亚洲另类 | 国产精品久久久久久久电影| 婷婷亚洲欧美| 欧美日韩国产亚洲二区| 观看免费一级毛片| 欧美中文日本在线观看视频| 1024手机看黄色片| 男人的好看免费观看在线视频| 日本黄大片高清| 老女人水多毛片| 性色avwww在线观看| 久久精品国产清高在天天线| 亚洲最大成人中文| 欧美色视频一区免费| 男女做爰动态图高潮gif福利片| 亚洲美女黄片视频| 女同久久另类99精品国产91| 国产熟女xx| netflix在线观看网站| 高潮久久久久久久久久久不卡| 国产久久久一区二区三区| 亚洲avbb在线观看| 亚洲欧美精品综合久久99| 18禁裸乳无遮挡免费网站照片| 日韩欧美一区二区三区在线观看| 国产毛片a区久久久久| 人妻丰满熟妇av一区二区三区| 一级作爱视频免费观看| 国产精品女同一区二区软件 | 身体一侧抽搐| 国产成+人综合+亚洲专区| 国产大屁股一区二区在线视频| 国产成人av教育| 精品久久久久久久末码| 亚洲无线在线观看| 少妇被粗大猛烈的视频| 精品一区二区三区人妻视频| 精品国产亚洲在线| 亚洲电影在线观看av| 九九在线视频观看精品| 亚洲第一区二区三区不卡| 日韩欧美 国产精品| 日韩国内少妇激情av| 国产av麻豆久久久久久久| 美女被艹到高潮喷水动态| 啦啦啦韩国在线观看视频| 国产av麻豆久久久久久久| 成年女人毛片免费观看观看9| 99久久精品国产亚洲精品| 午夜激情福利司机影院| 给我免费播放毛片高清在线观看| 久久天躁狠狠躁夜夜2o2o| 97热精品久久久久久| 听说在线观看完整版免费高清| 波多野结衣高清作品| 国产成人福利小说| 亚洲国产欧美人成| 一个人免费在线观看的高清视频| 午夜福利免费观看在线| 亚洲人成网站在线播放欧美日韩| 国产伦精品一区二区三区四那| 人妻丰满熟妇av一区二区三区| 午夜久久久久精精品| 人妻久久中文字幕网| 成年免费大片在线观看| 国产麻豆成人av免费视频| 免费在线观看亚洲国产| 成人鲁丝片一二三区免费| 欧美中文日本在线观看视频| 日本在线视频免费播放| 国产单亲对白刺激| 午夜精品久久久久久毛片777| 高清毛片免费观看视频网站| 欧美xxxx黑人xx丫x性爽| 国产精华一区二区三区| 真人一进一出gif抽搐免费| 亚洲精品影视一区二区三区av| 一个人免费在线观看电影| 国产精品久久久久久人妻精品电影| 男女视频在线观看网站免费| 亚洲精品久久国产高清桃花| 国模一区二区三区四区视频| 青草久久国产| 亚洲国产色片| 免费在线观看成人毛片| 97超级碰碰碰精品色视频在线观看| 婷婷丁香在线五月| 岛国在线免费视频观看| 男女做爰动态图高潮gif福利片| 久久久久亚洲av毛片大全| 国产午夜精品论理片| 国产v大片淫在线免费观看| 国模一区二区三区四区视频| 欧美日韩中文字幕国产精品一区二区三区| 村上凉子中文字幕在线| 日韩人妻高清精品专区| 国内少妇人妻偷人精品xxx网站| 男人舔奶头视频| 亚洲乱码一区二区免费版| 大型黄色视频在线免费观看| 国产亚洲精品久久久久久毛片| 欧美三级亚洲精品| 久久精品91蜜桃| 激情在线观看视频在线高清| 特级一级黄色大片| 老女人水多毛片| 看免费av毛片| 国产成人啪精品午夜网站| 两人在一起打扑克的视频| 最近在线观看免费完整版| 午夜福利在线观看吧| 日本免费a在线| 亚洲av电影在线进入| 亚洲精品粉嫩美女一区| 黄片小视频在线播放| 丁香欧美五月| 亚洲 国产 在线| 国内精品美女久久久久久| 九九在线视频观看精品| 亚洲专区国产一区二区| h日本视频在线播放| 老熟妇乱子伦视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产成+人综合+亚洲专区| 两性午夜刺激爽爽歪歪视频在线观看| 婷婷亚洲欧美| 亚洲av五月六月丁香网| 全区人妻精品视频| 久久热精品热| 日本与韩国留学比较| 搡老岳熟女国产| 亚洲片人在线观看| 亚洲,欧美精品.| 久久精品综合一区二区三区| 嫩草影院新地址| a在线观看视频网站| 一级黄片播放器| 亚洲欧美精品综合久久99| 99热这里只有是精品50| 欧美精品啪啪一区二区三区| 成人毛片a级毛片在线播放| 国产私拍福利视频在线观看| 一个人免费在线观看电影| 一区二区三区四区激情视频 | 精品人妻偷拍中文字幕| 国产精品久久久久久亚洲av鲁大| 久久精品国产亚洲av香蕉五月| 97超视频在线观看视频| 90打野战视频偷拍视频| 免费人成视频x8x8入口观看| 久久久久九九精品影院| 国产淫片久久久久久久久 | 国产精品永久免费网站| 中文字幕久久专区| 日本熟妇午夜| 成年人黄色毛片网站| 在线观看午夜福利视频| 18禁黄网站禁片午夜丰满| 成人亚洲精品av一区二区| 99久久精品一区二区三区| 一本精品99久久精品77| 婷婷丁香在线五月| 能在线免费观看的黄片| 18美女黄网站色大片免费观看| 99久久无色码亚洲精品果冻| 久久久久久国产a免费观看| 亚洲av.av天堂| 久久精品国产亚洲av香蕉五月| 波多野结衣高清无吗| 亚洲综合色惰| 天堂动漫精品| 我要搜黄色片| 午夜精品在线福利| 欧美3d第一页| 亚洲在线观看片| 国产毛片a区久久久久| 精品一区二区三区人妻视频| 少妇被粗大猛烈的视频| 精品久久国产蜜桃| 久久午夜福利片| 美女 人体艺术 gogo| 尤物成人国产欧美一区二区三区| 怎么达到女性高潮| 搞女人的毛片| 日本黄色视频三级网站网址| 婷婷精品国产亚洲av在线| 国产伦精品一区二区三区视频9| 丰满的人妻完整版| 五月玫瑰六月丁香| 欧美日韩中文字幕国产精品一区二区三区| 我要看日韩黄色一级片| 日日干狠狠操夜夜爽| 在线观看午夜福利视频| 别揉我奶头~嗯~啊~动态视频| 在线看三级毛片| 国产 一区 欧美 日韩| 看黄色毛片网站| 在线观看66精品国产| 午夜免费成人在线视频| 久久中文看片网| 啪啪无遮挡十八禁网站| 欧美色欧美亚洲另类二区| 国产高清三级在线| 伦理电影大哥的女人| 亚洲精品在线观看二区| 亚洲激情在线av| 欧美一区二区国产精品久久精品| 久久久久久久精品吃奶| 51国产日韩欧美| 国产在视频线在精品| 91九色精品人成在线观看| 国产探花极品一区二区| 亚洲成人免费电影在线观看| 亚洲精品成人久久久久久| 亚洲无线在线观看| 亚洲成av人片免费观看| 亚洲第一欧美日韩一区二区三区| 国产aⅴ精品一区二区三区波| 亚洲无线观看免费| 1000部很黄的大片| 热99在线观看视频| 国产伦人伦偷精品视频| 欧美激情在线99| 欧美成狂野欧美在线观看| 久久久久性生活片| 女人被狂操c到高潮| 99久久成人亚洲精品观看| 亚洲国产日韩欧美精品在线观看| 国语自产精品视频在线第100页| av在线蜜桃| 精品一区二区免费观看| а√天堂www在线а√下载| 一级黄片播放器| 亚洲在线观看片| 99视频精品全部免费 在线| 精品久久久久久久久亚洲 | 在线天堂最新版资源| 免费看a级黄色片| 男女床上黄色一级片免费看| 在线观看午夜福利视频| 99热这里只有精品一区| 99久久99久久久精品蜜桃| 动漫黄色视频在线观看| 热99在线观看视频| 99国产综合亚洲精品| 真人一进一出gif抽搐免费| 久久久久精品国产欧美久久久| 成人高潮视频无遮挡免费网站| 中文字幕av成人在线电影| 99热只有精品国产| 色精品久久人妻99蜜桃| 国产视频内射| 中文字幕久久专区| 日本黄色视频三级网站网址| 国产精品,欧美在线| 最近中文字幕高清免费大全6 | 精品国产亚洲在线| 精品一区二区三区人妻视频| 精品人妻视频免费看| 日本与韩国留学比较| 亚洲国产精品999在线| 少妇被粗大猛烈的视频| 久久精品久久久久久噜噜老黄 | 日本免费a在线| 男人的好看免费观看在线视频| 亚洲男人的天堂狠狠| 伊人久久精品亚洲午夜| av天堂在线播放| 亚洲国产欧美人成| 嫁个100分男人电影在线观看| 亚洲av电影不卡..在线观看| 熟妇人妻久久中文字幕3abv| 精品人妻偷拍中文字幕| 91麻豆精品激情在线观看国产| 我要看日韩黄色一级片| 啪啪无遮挡十八禁网站| 国产精华一区二区三区| 又爽又黄无遮挡网站| 成人一区二区视频在线观看| 黄色配什么色好看| 欧美又色又爽又黄视频| 国产精品一及| 嫩草影院新地址| bbb黄色大片| 亚洲天堂国产精品一区在线| 欧美乱色亚洲激情| 悠悠久久av| 九九久久精品国产亚洲av麻豆| 亚洲性夜色夜夜综合| 久久午夜亚洲精品久久| 91麻豆av在线| 午夜免费男女啪啪视频观看 | 午夜亚洲福利在线播放| 简卡轻食公司| 最近最新中文字幕大全电影3| 亚洲国产精品合色在线| 2021天堂中文幕一二区在线观| 国产精品一区二区三区四区久久| 天天躁日日操中文字幕| 欧美3d第一页| 精品久久久久久久人妻蜜臀av| 国产精品1区2区在线观看.| av在线观看视频网站免费| 99在线人妻在线中文字幕| 俺也久久电影网| 直男gayav资源| 亚洲精品一卡2卡三卡4卡5卡| 国产不卡一卡二| 性色av乱码一区二区三区2| 桃色一区二区三区在线观看| 国产激情偷乱视频一区二区| 日韩欧美 国产精品| 男女床上黄色一级片免费看| 1024手机看黄色片| 国产 一区 欧美 日韩| 两性午夜刺激爽爽歪歪视频在线观看| ponron亚洲| 美女高潮喷水抽搐中文字幕| 成人无遮挡网站| 亚洲不卡免费看| 欧美最新免费一区二区三区 | 国产高清激情床上av| 欧美在线一区亚洲| 日本三级黄在线观看| 免费无遮挡裸体视频| 国内揄拍国产精品人妻在线| 国产精品女同一区二区软件 | 成年人黄色毛片网站| 国产综合懂色| 亚洲国产高清在线一区二区三| 一进一出好大好爽视频| 黄色视频,在线免费观看| 久久午夜福利片| 老司机午夜十八禁免费视频| 久久精品国产亚洲av天美| 免费一级毛片在线播放高清视频| 久久国产精品影院| 日韩精品青青久久久久久| 国产精品三级大全| 亚洲精华国产精华精| 欧美高清性xxxxhd video| 天堂av国产一区二区熟女人妻| 亚洲第一区二区三区不卡| 欧美黄色片欧美黄色片| 亚洲欧美清纯卡通| 久久久久久久亚洲中文字幕 | 亚洲七黄色美女视频| 99久久精品国产亚洲精品| 91九色精品人成在线观看| 亚洲真实伦在线观看| 深爱激情五月婷婷| 精品午夜福利视频在线观看一区| 欧美激情久久久久久爽电影| 3wmmmm亚洲av在线观看| av视频在线观看入口| 亚洲第一欧美日韩一区二区三区| 9191精品国产免费久久| 国产高潮美女av| 欧美精品国产亚洲| 午夜福利视频1000在线观看| 国产久久久一区二区三区| 少妇高潮的动态图| 男女下面进入的视频免费午夜| 69av精品久久久久久| 亚洲久久久久久中文字幕| 国产精品国产高清国产av| avwww免费| 国产人妻一区二区三区在| 人妻制服诱惑在线中文字幕| 国产精品av视频在线免费观看| 一区二区三区免费毛片| 国产精品久久电影中文字幕| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧洲综合997久久,| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 午夜a级毛片| 久久久国产成人免费| 99久久99久久久精品蜜桃| av欧美777| 精品欧美国产一区二区三| 欧美日本视频| 亚洲精品成人久久久久久| 一边摸一边抽搐一进一小说| 十八禁国产超污无遮挡网站| 在线国产一区二区在线| 一区福利在线观看| 最近在线观看免费完整版| 搡老熟女国产l中国老女人| 久久婷婷人人爽人人干人人爱| 亚洲最大成人手机在线| 欧美日本视频| 男女下面进入的视频免费午夜| 最后的刺客免费高清国语| 国产在线男女| 国产主播在线观看一区二区| 麻豆成人午夜福利视频| 成人精品一区二区免费| 午夜免费成人在线视频| 极品教师在线免费播放| 精品久久久久久成人av| 全区人妻精品视频| 国产一区二区在线av高清观看| 色综合欧美亚洲国产小说| 国产一级毛片七仙女欲春2| www.999成人在线观看| 国产淫片久久久久久久久 | 在线观看舔阴道视频| av欧美777| 成年人黄色毛片网站| 国产亚洲精品av在线| 欧美bdsm另类| 国产亚洲精品久久久com| 亚洲中文日韩欧美视频| 97人妻精品一区二区三区麻豆| 十八禁人妻一区二区| 色综合婷婷激情| 久久亚洲精品不卡| bbb黄色大片| 亚洲精品一卡2卡三卡4卡5卡| 男女床上黄色一级片免费看| 国产成人av教育| 变态另类成人亚洲欧美熟女| 精品午夜福利在线看| 欧美黄色淫秽网站| 亚洲成人免费电影在线观看| 国产一区二区三区在线臀色熟女| 亚洲人成网站在线播| 久久久久久九九精品二区国产| 亚洲七黄色美女视频| 国产成人a区在线观看| 少妇人妻一区二区三区视频| 亚洲无线观看免费| 欧美+日韩+精品| 男女下面进入的视频免费午夜| 久久6这里有精品| 国产精品一区二区免费欧美| 在线天堂最新版资源| 丰满乱子伦码专区| 黄色丝袜av网址大全| 免费看a级黄色片| 在现免费观看毛片| 日韩精品中文字幕看吧| 高清日韩中文字幕在线| 麻豆av噜噜一区二区三区| 久久香蕉精品热| 国产黄色小视频在线观看| 亚洲在线观看片| 欧美bdsm另类| 无人区码免费观看不卡| 亚洲人成网站在线播放欧美日韩| 18禁在线播放成人免费| 精品久久久久久久久久免费视频| www日本黄色视频网| 国产美女午夜福利| 熟女电影av网| 午夜亚洲福利在线播放| 99热这里只有是精品50| 男人舔女人下体高潮全视频| 热99re8久久精品国产| 亚洲人成网站在线播放欧美日韩| 日韩欧美免费精品| 亚洲三级黄色毛片| 精品国产三级普通话版| 一进一出抽搐动态| 看十八女毛片水多多多| 精品久久国产蜜桃| 国产成人啪精品午夜网站| 不卡一级毛片| 性色av乱码一区二区三区2| 在线免费观看不下载黄p国产 | 亚洲精华国产精华精| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩无卡精品| 久久午夜亚洲精品久久| 日本熟妇午夜| 国产黄a三级三级三级人| 中文字幕免费在线视频6| 久久热精品热| 在线观看午夜福利视频| 色5月婷婷丁香| www.999成人在线观看| 校园春色视频在线观看| 嫩草影视91久久| 一边摸一边抽搐一进一小说| 亚洲无线观看免费| 五月玫瑰六月丁香| 亚洲avbb在线观看| 国产三级中文精品| 国产亚洲av嫩草精品影院| 少妇高潮的动态图| 亚洲av电影在线进入| 九色国产91popny在线| 成人三级黄色视频| 免费在线观看亚洲国产| 亚洲精品亚洲一区二区| 18禁黄网站禁片午夜丰满| 身体一侧抽搐| 91久久精品国产一区二区成人| 中文字幕av在线有码专区| 欧美黑人欧美精品刺激| 日韩欧美国产在线观看| 日韩中字成人| 两人在一起打扑克的视频| 欧美一区二区精品小视频在线| 亚洲欧美清纯卡通| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av日韩精品久久久久久密| 九色国产91popny在线| 日本 av在线| 老司机午夜福利在线观看视频| 自拍偷自拍亚洲精品老妇| 亚洲人成电影免费在线| 亚洲一区二区三区色噜噜| 最近最新免费中文字幕在线| 亚洲欧美日韩东京热| 国产麻豆成人av免费视频| 999久久久精品免费观看国产| 男女视频在线观看网站免费| 91狼人影院| 变态另类丝袜制服| 观看免费一级毛片| 久久精品国产自在天天线| 国产视频内射| 神马国产精品三级电影在线观看| 别揉我奶头~嗯~啊~动态视频| 久久精品91蜜桃| 18禁裸乳无遮挡免费网站照片| 欧美绝顶高潮抽搐喷水| 精品午夜福利在线看| 久久人妻av系列| 久久精品国产清高在天天线| 午夜福利成人在线免费观看| 深夜a级毛片| 99热这里只有是精品在线观看 | 露出奶头的视频| 国产精品电影一区二区三区| 一二三四社区在线视频社区8| 国产又黄又爽又无遮挡在线| 久久精品国产清高在天天线| 亚洲美女搞黄在线观看 | 国产v大片淫在线免费观看| 亚洲欧美日韩无卡精品| 成人鲁丝片一二三区免费| 日韩免费av在线播放| 俺也久久电影网| 少妇高潮的动态图| 18禁在线播放成人免费| 美女被艹到高潮喷水动态| 十八禁网站免费在线| 别揉我奶头 嗯啊视频| 亚洲国产精品sss在线观看| 成年免费大片在线观看| 国产高清视频在线播放一区| 在线观看午夜福利视频| 亚洲av五月六月丁香网| 美女xxoo啪啪120秒动态图 | 最新中文字幕久久久久| 亚洲一区高清亚洲精品| 精品人妻偷拍中文字幕| 精品人妻熟女av久视频| 亚洲精品亚洲一区二区| 日韩亚洲欧美综合| 久久热精品热| 狂野欧美白嫩少妇大欣赏| 成人国产一区最新在线观看| 国产成人av教育| 久久精品人妻少妇| 99精品久久久久人妻精品| 成人av一区二区三区在线看| 国产淫片久久久久久久久 | 他把我摸到了高潮在线观看| 亚洲av熟女| 如何舔出高潮| 国产爱豆传媒在线观看| 在线免费观看不下载黄p国产 | av在线蜜桃| 欧美日韩福利视频一区二区| 人妻夜夜爽99麻豆av| 给我免费播放毛片高清在线观看| 亚洲 欧美 日韩 在线 免费| 国产亚洲av嫩草精品影院| 成年女人看的毛片在线观看| 久久九九热精品免费| 国产一区二区三区视频了| 最近在线观看免费完整版| 91久久精品电影网| 18禁黄网站禁片午夜丰满| 国产精品久久久久久亚洲av鲁大| 国内精品久久久久久久电影| 亚洲精品一区av在线观看| 成人鲁丝片一二三区免费| 亚洲国产精品sss在线观看| 精品午夜福利视频在线观看一区| 午夜福利欧美成人| 国产久久久一区二区三区| 一进一出好大好爽视频| 成人美女网站在线观看视频| 99在线人妻在线中文字幕| 波多野结衣巨乳人妻| 深爱激情五月婷婷| 日本熟妇午夜|