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

    4 種輪作模式下長期連作煙田土壤微生態(tài)的響應(yīng)特征

    2024-06-16 19:51:19周旭東韓天華申云鑫施竹鳳賀彪楊明英裴衛(wèi)華何永宏楊佩文
    關(guān)鍵詞:植煙土壤微生物群落土壤酶活性

    周旭東 韓天華 申云鑫 施竹鳳 賀彪 楊明英 裴衛(wèi)華 何永宏 楊佩文

    摘要:為探究烤煙與不同作物輪作對長期連作植煙土壤質(zhì)量調(diào)控的效果,以連作10年以上的植煙田塊為對照(CK),設(shè)置大麥-烤煙(YCDM)、大蒜-烤煙(YCDS)、油菜-烤煙(YCYC)和蠶豆-烤煙(YCCD)共4種不同作物與烤煙輪作,分析不同輪作處理下土壤的理化性質(zhì)、酶活性及微生物群落結(jié)構(gòu)。基于相關(guān)性分析和冗余分析,解析土壤理化性質(zhì)和酶活性與微生物群落結(jié)構(gòu)間的關(guān)系。結(jié)果表明,輪作處理的土壤容重較CK降低26.58%~30.29%,總孔隙度增加21.13%~48.26%;其中烤煙-蠶豆(YCCD)輪作模式下土壤的pH、水解氮(available nitrogen,AN)、有效磷(available phosphorus,AP)和速效鉀( available potassium,AK)含量較CK顯著提高11.84%、30.57%、6.42%和41.5 1%。輪作處理的土壤過氧化氫酶(catalase,CAT)、蔗糖酶(invertase,INV)、脲酶(urease,URE)和酸性磷酸酶(acid phosphatase,ACP)的活性較CK 分別提高16.81%~42.35%、38.09%~51.48%、7.69%~64.29%和5.82%~76.33%,蠶豆-烤煙(YCCD)輪作模式的提升效果最為顯著。高通量測序結(jié)果表明,輪作處理下細(xì)菌的OTUs數(shù)量顯著高于CK。α多樣性顯示,不同輪作處理間的細(xì)菌群落豐富度和多樣性差異顯著;真菌群落的豐富度差異顯著,而多樣性差異不顯著。β多樣性分析顯示,不同輪作處理之間土壤真菌群落差異較小,細(xì)菌群落差異較大。其中,蠶豆-烤煙(YCCD)輪作處理的細(xì)菌、真菌群落與CK差異較大。在門水平,烤煙-蠶豆(YCCD)輪作使酸桿菌門(Acidobacteria)和放線菌門(Actinobacteria)的相對豐度增加,子囊菌門(Ascomycota)的相對豐度減少。相關(guān)性分析和冗余分析表明,URE、AN、AK和pH是影響土壤微生物群落結(jié)構(gòu)的關(guān)鍵因子。綜上所述,烤煙與其他作物合理輪作可提高土壤速效養(yǎng)分含量和土壤酶活性,調(diào)節(jié)土壤微生物群落結(jié)構(gòu),消減烤煙連作障礙,最終達(dá)到烤煙穩(wěn)產(chǎn)、增產(chǎn)的目的。

    關(guān)鍵詞:植煙土壤;輪作;土壤理化性質(zhì);土壤酶活性;微生物群落

    doi:10.13304/j.nykjdb.2022.0743

    中圖分類號:S572 文獻(xiàn)標(biāo)志碼:A 文章編號:1008‐0864(2024)03‐0174‐14

    煙草(Nicotiana tabacum L.)是我國重要的經(jīng)濟(jì)作物之一[1]。據(jù)統(tǒng)計,云南是全國種植烤煙面積最大的省份,也是全球最大的優(yōu)質(zhì)烤煙產(chǎn)區(qū),云南煙葉占全國總量的45%左右[2]。隨著云南省農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和種植制度變革,烤煙種植業(yè)發(fā)展迅猛,隨之而來的烤煙連作現(xiàn)象也日益突出。同一作物或其近緣種長期連續(xù)種植會造成土壤養(yǎng)分異常積累或過度消耗,病原微生物迅速繁衍,土壤微生物種群結(jié)構(gòu)失衡,導(dǎo)致作物產(chǎn)量、品質(zhì)降低及病蟲害頻發(fā),嚴(yán)重影響烤煙產(chǎn)業(yè)的持續(xù)健康發(fā)展[3]。因此,如何合理有效地緩解云南烤煙連作障礙、均衡土壤地力水平以及改善土壤微生態(tài)環(huán)境,成為目前亟待解決的重要問題。

    輪作可以在一些方面抑制連作障礙的產(chǎn)生,對作物生長起到一定的積極作用。輪作通過用地養(yǎng)地相互配合的方法,在提供正常設(shè)施農(nóng)業(yè)生產(chǎn)的同時,利用不同作物對土壤養(yǎng)分吸收利用的差異,使用不同作物進(jìn)行輪換栽培,有利于平衡土壤中的營養(yǎng)元素,促進(jìn)土壤生態(tài)系統(tǒng)健康發(fā)展,以此達(dá)到土地資源轉(zhuǎn)換和保持土壤健康的目的[4]。近年來,輪作的研究重點正在從宏觀的土壤質(zhì)量、作物產(chǎn)量和病蟲害防治向微觀的土壤微生態(tài)和土壤成分互作方面延伸,不同輪作作物的根系與根際土壤微生物的相互作用可能比地上部分的相互作用更為重要。輪作可以提高土壤微生物群落多樣性,而微生物群落多樣性又與作物生長發(fā)育息息相關(guān)[5]。國內(nèi)外學(xué)者一直密切關(guān)注輪作對于作物生長、土壤質(zhì)量以及土壤微生態(tài)等方面的影響。朱維偉等[6]研究發(fā)現(xiàn),不同作物與黃瓜輪作有利于黃瓜幼苗生物量的積累,同時可以抑制土壤鹽漬化,減緩?fù)寥栏粻I養(yǎng)化,均衡土壤養(yǎng)分,降低黃瓜病害,為黃瓜生長構(gòu)建出相對優(yōu)良的生態(tài)環(huán)境。趙慶雷等[7]發(fā)現(xiàn),水旱輪作模式可以明顯優(yōu)化花生產(chǎn)量結(jié)構(gòu),顯著提高根瘤菌數(shù)量,對花生果腐病和蠐螬有較理想的防控效果。徐新雯等[8]研究發(fā)現(xiàn),烤煙-大蒜輪作可以提高煙田土壤細(xì)菌群落結(jié)構(gòu)的調(diào)節(jié)能力,使煙田的發(fā)病率由往年的50.0%降低到6.5%。陽顯斌等[9]發(fā)現(xiàn),煙蒜輪作和煙蒜套作能顯著提高根際土壤中細(xì)菌、放線菌、解磷菌和解鉀菌的數(shù)量,維持良好的土壤質(zhì)量狀態(tài)。蘇燕等[10]研究表明,輪作栽培可以提高馬鈴薯根際土壤細(xì)菌群落的豐富度指數(shù)及多樣性指數(shù)。姚小東等[11]研究表明,花生輪作種植可以改變土壤微生物群落結(jié)構(gòu),增加有益細(xì)菌比例,降低有害真菌比例。由此可見,合理的輪作可以有效均衡土壤養(yǎng)分,改善土壤環(huán)境,有助于烤煙良好循環(huán)種植。

    盡管對于烤煙合理輪作的研究已有相關(guān)報道,但針對云南區(qū)域內(nèi)烤煙輪作的研究較少。不同區(qū)域的種植環(huán)境因其受到光、溫、水、氣以及人為因素的差異,其土壤微環(huán)境也不同。因此,為了解決烤煙輪作作物選配不合理導(dǎo)致的土壤質(zhì)量退化問題,分別以油菜、大麥、蠶豆、大蒜4種作物與烤煙進(jìn)行輪作,解析烤煙不同輪作模式對土壤理化性狀、酶活性以及微生物多樣性等土壤質(zhì)量主控因子的影響及調(diào)控作用,為烤煙合理輪作提供理論依據(jù)和實踐支持。

    1 材料與方法

    1.1 試驗地點

    試驗地位于云南省麗江市石鼓鎮(zhèn)魯瓦村(26°50′N, 99°54′E),海拔1 859 m。地處金沙江河谷地帶,屬江邊干熱河谷氣候。其季風(fēng)特點明顯,四季分明,春季干旱少雨,夏季炎熱多雨,秋季涼爽干燥,冬季寒冷有雪。年均降水量868.2 mm,主要集中于4—8月;年均氣溫21.4 ℃;年日照時數(shù)2 321~2 554 h;年無霜期191~310 d。試驗地土壤為砂質(zhì)山地黃壤土,pH 6.12、有機質(zhì)43.70 g·kg-1、全氮2.36 g·kg-1、水解氮221.22 mg·kg-1、全磷1.20g·kg-1、有效磷100.41 mg·kg-1、全鉀16.10 g·kg-1、速效鉀267.38 mg·kg-1、電導(dǎo)率115.9 μS·cm-1。供試烤煙品種為‘云煙99。

    1.2 試驗設(shè)計

    選擇冬季空閑且連續(xù)種植烤煙10年以上的田塊,以烤煙連作為對照(CK),設(shè)置4個輪作處理,分別為烤煙-蠶豆輪作(YCCD)、烤煙-大麥輪作(YCDM)、烤煙-大蒜輪作(YCDS)和烤煙-油菜輪作(YCYC),共計5個處理,每個處理3次重復(fù)。每小區(qū)3壟(3 m×7 m),于周圍設(shè)立保護(hù)行。前茬作物油菜、大麥、大蒜和蠶豆均施用腐熟細(xì)碎的有機堆肥3 000 kg·hm-2和草木灰850 kg·hm-2,旋耕后平栽。大麥化肥施用量為N 285 kg·hm-2、P2O5 120 kg·hm-2、K2O5 285 kg·hm-2;油菜化肥施用量為N 150 kg·hm-2、P2O5 75 kg·hm-2、K2O 75 kg·hm-2;大蒜化肥施用量為N 225 kg·hm-2、P2O5 150 kg·hm-2、K2O 300 kg·hm-2;蠶豆化肥施用量為N 75 kg·hm-2、P2O5 75 kg·hm-2、K2O 120 kg·hm-2。其中,氮肥和鉀肥為硝酸鉀肥,磷肥為鈣鎂磷肥。大麥?zhǔn)斋@后將秸稈還田,其余作物收獲后將秸稈全部移出農(nóng)田,進(jìn)行機械滅茬。烤煙施用的有機肥為酒糟有機肥(N+P2O5+K2O=6.2%),化肥為烤煙專用肥(N∶P2O5∶k2O=15∶15∶15,質(zhì)量比)。

    2020年種植油菜、大麥、大蒜、蠶豆,收獲后于2021 年4 月15 日移栽煙苗,株行距為50 cm×50 cm,2021年8月15日進(jìn)行采收;試驗中各類作物采用常規(guī)田間管理。

    1.3 土壤樣品采集與分析

    1.3.1 土壤樣品的采集 每個小區(qū)采用“W”形取樣法隨機取5株煙,每株煙為1個樣點,共計75個取樣點。于2021年8月中旬(收獲期),首先去除煙株根際表面覆蓋物,用環(huán)刀采集壟上近根(0—10 cm)土壤,用于土壤容重、孔隙度的測定。采集土壤樣品時,用75%乙醇對土鉆進(jìn)行滅菌,然后用土鉆采集0—20 cm土層土壤,剔除石塊、草根、秸稈等雜物,采用四分法混合均勻,取1 kg裝入無菌密封袋中帶回實驗室。每個密封袋中的土壤樣品分為兩部分,一部分自然風(fēng)干,用于土壤理化性質(zhì)的測定;另一部分裝入50 mL離心管,超低溫凍干機凍干后保存于-80 ℃超低溫冰箱,用于細(xì)菌、真菌群落的高通量測序。

    1.3.2 土壤理化性質(zhì)的測定 土壤容重和孔隙度采用環(huán)刀法測定,土壤pH采用酸度計測定,土壤有機質(zhì)(soil organic matter,SOM)含量采用重鉻酸鉀容量法測定,電導(dǎo)率(electrical conductivity,EC)采用電導(dǎo)率儀測定,全氮(total nitrogen,TN)含量采用凱氏定氮法測定,堿解氮(available nitrogen,AN)含量采用堿解擴(kuò)散法測定,全磷(totalphosphorus,TP)采用H2SO4-HClO4 鉬銻抗比色法測定,速效磷(available phosphorus,AP)含量采用碳酸氫鈉提取法測定;全鉀(total potassium,TK)采用NaOH 熔融- 火焰光度法測定,速效鉀(available potassium,AK)含量采用采用 NH4OAc浸提結(jié)合火焰光度法測定[12]。

    1.3.3 土壤酶活性的測定 土壤蔗糖酶(invertase,INV)、酸性磷酸酶(acid phosphatase,ACP)、過氧化氫酶(catalase,CAT)和脲酶(urease,URE)活性采用試劑盒(蘇州瑞思公司生物科技有限公司)微板法測定。

    1.3.4 微生物多樣性分析 采用E.Z.N.A.?soilDNA kit 抽提樣品總DNA,1% 瓊脂糖凝膠電泳檢測DNA 質(zhì)量,NanoDrop 測定DNA 含量和純度。以此DNA 為模板,運用引物338F(5-ACTCCTACGGGAGGCAGCAG-3) 、806R(5-GGACTACHVGGGTWTCTAAT-3)和SSU0817F(5-TTAGCATGGAATAATRRAATA GGA-3)、81196R(5-TCTGGACCTGGTGAGTTT CC-3)分別對細(xì)菌16S rRNA 基因的V3~V4 可變區(qū)、真菌18S rRNA基因的V3~V4高變區(qū)片段進(jìn)行PCR擴(kuò)增引物序列為。PCR程序:95 ℃ 2 min;95 ℃ 30 s,55 ℃ 30 s,72 ℃ 30 s,25個循環(huán); 72 ℃ 5 min。每個樣本3次重復(fù),將同一樣本的PCR產(chǎn)物混合后用2% 瓊脂糖凝膠電泳檢測,使用AxyPrepDNA凝膠回收試劑盒(AXYGEN 公司)切膠回收PCR 產(chǎn)物,Tris-HCl 洗脫;2% 瓊脂糖電泳檢測。參照電泳初步定量結(jié)果,將PCR 產(chǎn)物用QuantiFluorTM-ST藍(lán)色熒光定量系統(tǒng)(Promega 公司)進(jìn)行定量檢測,按照每個樣本的測序量要求,進(jìn)行相應(yīng)比例的混合,送至上海美吉生物醫(yī)藥科技有限公司于Illumina Miseq PE300、PE250 平臺進(jìn)行測序分析。

    1.4 數(shù)據(jù)統(tǒng)計與分析

    采用Excel 2016和SPSS 22.0進(jìn)行數(shù)據(jù)分析,用OriginPro 2021繪圖。采用單因素方差分析評價不同輪作模式間土壤理化性質(zhì)、酶活性以及微生物多樣性的差異(P<0.05)。采用CANOCO4.5軟件(Microcomputer Power,Ithaca,NY,USA)分析環(huán)境因子對微生物群落分布的影響。

    2 結(jié)果與分析

    2.1 不同輪作模式對土壤理化性質(zhì)的影響

    2.1.1 不同輪作模式對土壤物理性質(zhì)的影響 不同輪作模式對土壤物理性狀造成的影響如圖1所示。4種輪作模式的土壤容重顯著低于CK,分別為CK的76.91%、83.02%、69.80%和78.60%;除了YCDS模式外,其余3種輪作模式的土壤總孔隙度均顯著高于CK,分別為CK 的148.25%、123.23%和121.10%。與CK相比,4種輪作模式均顯著增加了土壤毛孔隙度,其中,YCYC 模式的增幅最大,為CK的292.83%。YCYC模式下土壤的非毛孔隙度較CK顯著降低。

    2.1.2 不同輪作模式對土壤化學(xué)性質(zhì)的影響 不同輪作模式對土壤化學(xué)性質(zhì)的影響如表1所示。除YCDM 模式外,其余3 種模式下土壤的pH 較CK顯著上升,其中,YCCD模式的pH最高,為CK處理的131.97%。YCDS模式的土壤電導(dǎo)率顯著高于CK處理,為CK的132.37%;YCYC模式的土壤電導(dǎo)率顯著低于CK,為CK 的63.58%;其余輪作模式與CK差異不顯著。YCDM輪作模式土壤的有機質(zhì)含量顯著高于CK,為CK 的116.28%。4 種輪作模式土壤、全氮、全磷和速效鉀含量與CK 處理差異不顯著。YCCD 模式土壤堿解氮含量顯著高于CK;其余3 種輪作模下則顯著低于CK 處理。YCDM、YCYC 和YCCD 模式的土壤有效磷含量顯著高于CK。YCDM、YCDS和YCCD模式的土壤速效鉀含量顯著高于;YCYC模式土壤速效鉀含量則顯著低于CK。

    2.2 不同輪作模式對土壤酶活性的影響

    不同輪作模式對土壤酶活性的影響如圖2所示。YCDM、YCYC和YCCD模式土壤過氧化氫酶活性顯著增加,其中YCDM 模式最高,為CK 的203.31%;YCDM、YCDS和YCCD模式土壤蔗糖酶活性顯著增加,其中YCCD 模式最高,為CK 的258.86%;YCDS和YCCD模式土壤脲酶活性顯著增加,分別為CK的175.89%和165.25%,其他2種輪作模式土壤脲酶活性雖有升高,但與CK差異不顯著;YCYC和YCCD模式土壤酸性磷酸酶活性顯著上升,其中YCCD 模式最高,為CK 處理的433.66%。

    2.3 不同輪作模式對土壤微生物OTUs 數(shù)量的影響

    不同輪作模式土壤特有及共有細(xì)菌、真菌的OTUs數(shù)量如圖3所示。不同輪作模式下細(xì)菌共有OTUs數(shù)量為2 482個,YCCD模式特有的OTUs數(shù)量為224個;YCYC模式為77個;YCDM模式為103個;YCDS模式為185個;CK處理為49個。不同輪作模式下真菌共有OTUs數(shù)量為282個,YCCD模式特有的OTUs 數(shù)量為88 個;YCYC 模式為139個;YCDM 模式為69個;YCDS 模式為97 個;CK處理為134個。土壤細(xì)菌數(shù)量明顯大于真菌數(shù)量。在YCCD模式下,土壤特有細(xì)菌數(shù)量較CK顯著增多,而特有真菌數(shù)量明顯減少,表明YCCD模式能夠增加土壤細(xì)菌群落數(shù)量,降低真菌群落數(shù)量。

    2.4 不同輪作模式對土壤微生物α 多樣性的影響

    不同輪作模式土壤微生物α多樣性如表2所示。對于土壤細(xì)菌群落,YCDS模式下細(xì)菌的物種數(shù)目、Ace指數(shù)和Chao指數(shù)較CK顯著增加,而多樣性指數(shù)無顯著變化,表明YCDS模式可以增加土壤細(xì)菌豐富度,而對細(xì)菌多樣性則無明顯影響;YCCD模式下細(xì)菌的物種數(shù)目、Ace指數(shù)、Chao指數(shù)和Shannon 指數(shù)均較CK 相顯著增加,而Simpson指數(shù)顯著降低,表明YCCD模式不僅能夠增加土壤細(xì)菌豐富度,還能夠增加土壤細(xì)菌多樣性。對于土壤真菌群落,YCDM 和YCDS 模式下真菌的物種數(shù)目、Ace指數(shù)和Chao指數(shù)較CK顯著降低,而多樣性指數(shù)無明顯變化,表明YCDM 和YCDS模式可以顯著降低土壤真菌豐富度,而對真菌多樣性則無明顯影響。

    2.5 不同輪作模式對土壤微生物β 多樣性的影響

    由圖4A可知,第1坐標(biāo)軸和第2坐標(biāo)軸解釋了細(xì)菌總方差的34.22%;由圖4B可知,第1坐標(biāo)軸和第2坐標(biāo)軸解釋了真菌總方差的34.26%,表明不同輪作模式下的細(xì)菌、真菌群落結(jié)構(gòu)存在差異。不同輪作模式細(xì)菌群落之間距離較遠(yuǎn),且與CK距離較遠(yuǎn),而真菌群落距離較近;細(xì)菌群落的r 值為0.582 1,真菌群落的r 值為0.471 5。由此表明,不同輪作模式間細(xì)菌群落組成的差異性更為顯著。

    2.6 不同輪作模式對土壤微生物群落結(jié)構(gòu)的影響

    2.6.1 不同輪作模式對土壤細(xì)菌門、屬水平群落結(jié)構(gòu)的影響 由圖5A可知,在細(xì)菌門水平上,5個處理的優(yōu)勢菌門均為變形菌門(Proteobacteria)、酸桿菌門(Acidobacteria)、放線菌門(Actinobacteria)和綠彎菌門(Chloroflexi)。其中,變形菌門在YCDM、YCDS、YCYC和YCCD模式的相對豐度分別為31.61%、33.95%、32.72% 和28.41%,均高于CK(27.37%);酸桿菌門和放線菌門的相對豐度分別為17.93%~20.55% 和15.66%~25.77%;綠彎菌門在不同輪作模式處理間無明顯變化。4種輪作模式增加了土壤中變形菌門的相對豐度,降低了放線菌門的相對豐度,其中YCDS模式下土壤變形菌門相對豐度的增幅明顯高于其余3 種輪作模式。

    在屬水平下(圖5B),5個處理的優(yōu)勢菌屬均為Vicinamibacterales、嗜食單胞菌(Sphingomonas)、芽孢桿菌屬(Bacillus)和Gaiellales,其相對豐度分別為4.98%~6.55%、2.52%~3.21%、1.58%~2.93% 和1.73%~3.24%。與CK 相比,YCCD 模式使Vicinamibacterales、嗜食單胞菌(Sphingomonas)、芽孢桿菌屬(Bacillus)和Gaiellales 的相對豐度分別提高了0.1、0.2、1.5和0.2個百分點。

    2.6.2 不同輪作模式對土壤真菌門、屬水平群落結(jié)構(gòu)的影響 由圖6A可知,在門水平,不同輪作模式下均表現(xiàn)為子囊菌門(Ascomycota)、被孢霉門(Mortierellomycota)、擔(dān)子菌門(Basidiomycota)、unclassified_k__Fungi、油壺菌門(Olpidiomycota)和羅茲菌門/隱菌門(Rozellomycota)的相對豐度較高,其中,子囊菌門和被孢霉門相對豐度累計超過70%。與CK相比,子囊菌門在4種輪作模式下的相對豐度均有所降低,為32.10%~62.08%,其中,YCDS模式最低;被孢霉門在4種輪作模式下的相對豐度均有所增加,為13.38%~32.86%,其中,YCDS模式最高。

    在屬水平,排名前10的菌屬依次為:被孢霉菌屬(Mortierella)、unclassified_k__Fungi、鐮刀霉屬(Fusarium)、芽枝霉屬(Cladosporium)、赤霉菌屬(Gibberella)、Saitozyma、油壺菌屬(Olpidium)、錐蓋傘屬(Conocybe)、向基孢霉屬(Basipetospora)和尾柄孢殼屬(Cercophora)。YCCD 模式中芽枝霉屬和錐蓋傘屬相對豐度分別為11.13% 和9.36%,與CK差異較大。

    2.7 理化性質(zhì)對土壤微生物群落結(jié)構(gòu)的影響

    2.7.1 理化性質(zhì)與土壤微生物優(yōu)勢菌門間的相關(guān)性分析 理化性質(zhì)與土壤細(xì)菌優(yōu)勢菌門間的相關(guān)性分析(圖7A)表明,變形菌門的相對豐度與AK顯著負(fù)相關(guān);放線菌門的相對豐度與AN和AK顯著正相關(guān);擬桿菌門的相對豐度與AK顯著負(fù)相關(guān);酸桿菌門的相對豐度與AN和ACP顯著正相關(guān)。理化性質(zhì)與土壤真菌優(yōu)勢菌門間的相關(guān)性分析(圖7B)表明,子囊菌門的相對豐度與pH顯著正相關(guān),與TN 和TP 呈顯著負(fù)相關(guān);unclassified_k__Fungi的相對豐度與AN顯著負(fù)相關(guān)。

    2.7.2 土壤微生物優(yōu)勢菌門與理化性質(zhì)的冗余分析 土壤微生物優(yōu)勢菌門與理化性質(zhì)的冗余分析如圖8所示。對于細(xì)菌,第1軸(RDA1)解釋率為46.07%,第2軸(RDA2)的解釋率為35.42%,累積總解釋率為83.76%;對于真菌,第1軸(RDA1)解釋率為35.42%,第2軸(RDA2)解釋率為31.57%,累積總解釋率為66.99%,說明分析結(jié)果可靠。進(jìn)一步分析發(fā)現(xiàn),AN和AK是影響土壤細(xì)菌群落結(jié)構(gòu)的主要環(huán)境因子;AN、URE、pH、AK和TP是影響土壤真菌群落結(jié)構(gòu)的主要環(huán)境因子。

    3 討論

    3.1 不同輪作模式對土壤理化性質(zhì)和酶活性的影響

    作物長期連作會引起土壤質(zhì)量退化、土傳病害增加等問題[13]。合理的輪作制度不僅能夠提高土壤肥力,增加土壤有機質(zhì)和均衡土壤酸堿度,還能減少土傳病害的發(fā)生,為烤煙高質(zhì)量生產(chǎn)提供有力保障[14]。本研究發(fā)現(xiàn),輪作后的土壤容重較連作土壤明顯降低,而土壤總孔隙度則明顯升高,這可能是由于長期的連作,不重視養(yǎng)地,使得土壤耕層變淺、土壤板結(jié),最終導(dǎo)致土壤容重增加、土壤總孔隙度下降,而不同的輪作作物因其根系生長發(fā)育特性不同,促使土壤物理性狀發(fā)生變化??梢?,輪作能夠改善土壤松緊度,增加其通透性。

    輪作還能提升土壤酸堿度,4種輪作模式下土壤pH較長期連作土壤明顯升高,這可能是因為前茬作物種植時,施用了較多的鈣鎂磷肥、草木灰等堿性肥料,使土壤pH增加,而連作土壤長期施用煙草專用化肥(N∶P2O5∶K2O=15∶15∶15),氮素投入量過多,導(dǎo)致土壤酸化[15‐16]。輪作后土壤的有機質(zhì)、堿解氮、有效磷和速效鉀含量均較連作土壤有不同程度的升高,YCCD處理的效果最明顯,其次是YCDM處理。這是由于豆科植物具有固氮作用,可以優(yōu)化土壤養(yǎng)分結(jié)構(gòu),有效改善土壤質(zhì)量;在YCDM處理中,大麥稈還田是一項增加土壤有機質(zhì)、改良土壤結(jié)構(gòu)的重要舉措,顯著提高了土壤有機質(zhì)含量,這與前人研究結(jié)果一致[17‐18]。

    本研究中,YCDM、YCDS 和YCYC 模式土壤堿解氮含量較CK顯著降低;YCDS和YCCD模式土壤有效磷含量較CK顯著降低;YCCD模式土壤速效鉀含量較CK顯著降低,與前人研究結(jié)果一致[19‐20]。這可能是由于在土壤養(yǎng)分含量相對一致的情況下,不同輪作作物的生長特性及需肥量不同,導(dǎo)致不同處理間土壤養(yǎng)分含量出現(xiàn)差異;也可能是由于輪作提升了烤煙根系根長,提高了土壤養(yǎng)分與根系的接觸面積,增加了作物對土壤養(yǎng)分的攝入量,而不同輪作模式對烤煙根系發(fā)育的促進(jìn)作用不同,從而導(dǎo)致處理間土壤養(yǎng)分含量不同[21]。

    土壤中過氧化氫酶、蔗糖酶、脲酶和酸性磷酸酶的活性在一定程度上代表了土壤氧化還原反應(yīng)和碳、氮、磷的代謝強弱[22]。研究表明,不同作物輪作能夠增強土壤酶活性,而長期連作使得土壤酶活性受到抑制[23]。本研究表明,烤煙與不同作物輪作均提高了土壤中過氧化氫酶、蔗糖酶、脲酶和酸性磷酸酶活性,其中YCCD模式對土壤酶活性的提升效果最好。這可能是由于豆科植物與根瘤菌的共生關(guān)系有效地調(diào)節(jié)和改善了土壤微生物群落結(jié)構(gòu),增強了微生物活性,從而促進(jìn)了酶活性的提高[24‐25]。

    3.2 不同輪作模式對土壤微生物群落的影響

    土壤微生物在土壤質(zhì)量改善和促進(jìn)土壤養(yǎng)分循環(huán)等方面起重要作用[14]。輪作作物因其對水分、養(yǎng)分的吸收利用存在差異,能夠改變土壤微環(huán)境,進(jìn)而直接或間接影響土壤微生物群落多樣性[16]。本研究發(fā)現(xiàn),CK土壤細(xì)菌的OTUs數(shù)量和種類均較少,而土壤真菌的OTUs數(shù)量和種類則高于YCDM、YCDS和YCCD 輪作模式;烤煙輪作改變了土壤微生物群落豐富度,特別是增加了細(xì)菌群落豐富度。與CK相比,YCDS模式中變形菌門的相對豐度顯著提高,YCDS模式中放線菌門的相對豐度有所提高;4種輪作模式降低了子囊菌門的相對豐度,增加了擔(dān)子菌門的相對豐度,其中YCDS處理中子囊菌門相對豐度顯著低于CK,被孢霉門相對豐度顯著高于CK。這是由于烤煙在連作條件下根系分泌出的有毒、有害物質(zhì)積累,使根系受到傷害,影響了根際微生物的種群結(jié)構(gòu),使得土壤細(xì)菌種群數(shù)量降低、真菌種群數(shù)量上升,進(jìn)而破壞了根際微生物平衡,促使土壤從細(xì)菌型土壤轉(zhuǎn)變?yōu)椴唤】档恼婢屯寥繹26]。合理的輪作可以打破這一變化,提高土壤細(xì)菌群落的多樣性,隨著輪作年份的增加,輪作土壤中根際細(xì)菌群落含量高于連作土壤的趨勢越來越明顯[27]。

    烤煙連作降低了土壤細(xì)菌群落的Ace指數(shù)、Chao指數(shù)和Shannon指數(shù),而增加了真菌群落多樣性的Ace 指數(shù)、Chao 指數(shù)和Shannon 指數(shù)。本研究發(fā)現(xiàn),YCDS模式下細(xì)菌的物種數(shù)量、Ace指數(shù)和Chao指數(shù)均顯著高于CK,而多樣性指數(shù)則無顯著性差異,表明YCDS模式有利于增加土壤細(xì)菌豐富度,而對細(xì)菌多樣性無明顯影響。YCCD模式下細(xì)菌的物種數(shù)目、Ace指數(shù)、Chao指數(shù)和Shannon 指數(shù)均顯著高于CK,而Simpson 指數(shù)顯著低于CK,表明YCCD模式不僅能夠增加土壤細(xì)菌豐富度,還能夠增加土壤細(xì)菌多樣性。這是由于蠶豆為豆科植物,可以調(diào)節(jié)土壤C/N,增加土壤碳源和氮源,進(jìn)而影響土壤微生物活性,改變土壤微生物種群和結(jié)構(gòu),使得土壤中細(xì)菌Shannon指數(shù)升高[28]。YCDM和YCDS模式下真菌的物種數(shù)目、Ace 指數(shù)和Chao 指數(shù)均顯著低于CK,而多樣性指數(shù)無明顯變化,表明YCDM 和YCDS模式可以顯著降低土壤真菌群落豐富度,而對真菌多樣性無明顯影響。這可能是因為大麥稈還田能夠增加土壤中有機質(zhì)含量,改良土壤結(jié)構(gòu),使土壤疏松、孔隙度增加、容量減輕,促進(jìn)微生物活力和作物根系發(fā)育,降低病原菌數(shù)量[29];也可能是因為大蒜的根際分泌物——大蒜素能夠殺死土傳病害病原菌,而這些病原菌大多數(shù)為真菌,故大蒜輪作土壤中真菌群落數(shù)量減少[30]。因此,不同輪作模式下土壤細(xì)菌β 多樣性總體上大于真菌。相關(guān)性分析和冗余分析顯示,土壤中脲酶活性、堿解氮、速效鉀含量和pH是影響土壤微生物群落數(shù)量的關(guān)鍵因子,且與子囊菌門、放線菌門和酸桿菌門相對豐度呈顯著正相關(guān),表明土壤輪作后能夠明顯改變土壤養(yǎng)分含量及酶活性,進(jìn)而提高土壤中有益菌群的數(shù)量,改善土壤質(zhì)量,促進(jìn)烤煙健康生長。

    參考文獻(xiàn)

    [1] 芶久蘭,顧小鳳,張萌,等.不同烤煙種植模式對貴州土壤養(yǎng)

    分、酶活性及細(xì)菌群落結(jié)構(gòu)的影響[J]. 核農(nóng)學(xué)報,2022,

    36(7):475-484.

    GOU J L, GU X F, ZHANG M, et al .. Effects of different fluecured

    tobacco planting patterns on nutrients, enzyme activities

    and bacterial community structure in soil of Guizhou province [J]. J.

    Nucl. Agric. Sci., 2022, 36(7):475-484.

    [2] 魯韋坤,逄濤,余凌翔,等.市場經(jīng)濟(jì)對云南規(guī)劃煙區(qū)耕地資

    源的影響分析[J].中國農(nóng)學(xué)通報,2021,37(5):137-142.

    LU W K, PANG T, YU L X, et al .. The influence of market

    economy on cultivated land resources in planned tobacco area

    of Yunnan [J]. Chin. Agric. Sci. Bull., 2021, 37(5):137-142.

    [3] 陳華,趙文軍,王正旭,等.不同輪作模式下氮素調(diào)控對烤煙

    產(chǎn)質(zhì)量及氮肥利用的影響[J]. 河南農(nóng)業(yè)科學(xué),2021,50(9):

    87-95.

    CHEN H, ZHAO W J, WANG Z X, et al .. Effects of nitrogen

    management on yield,quality and nitrogen utilization of fluecured

    tobacco under different rotation patterns [J]. J. Henan

    Agric. Sci., 2021, 50(9):87-95.

    [4] CHEN S, QI G, LUO T, et al .. Continuous‐cropping tobacco

    caused variance of chemical properties and structure of

    bacterial network in soils [J]. Land Degrad. Dev., 2018, 29(11):

    4106-4120.

    [5] 陳浩,魏立本,王亞麒,等.烤煙不同種植施肥模式對土壤養(yǎng)

    分、酶活性及細(xì)菌多樣性的影響[J].南方農(nóng)業(yè)學(xué)報,2019,

    50(5):982-989.

    CHEN H, WEI L B, WANG Y Q, et al .. Effects of different

    planting and fertilizing modes on soil nutrient,enzyme activity

    and bacterial diversity of tobacco [J]. J. Southern Agric., 2019,

    50(5):982-989.

    [6] 朱維偉. 不同輪作模式對黃瓜幼苗生長及土壤環(huán)境的影

    響[D].哈爾濱:東北農(nóng)業(yè)大學(xué),2020.

    ZHU W W. Effects of Different rotation sysetems on growth

    and soil environment of cucumber seedling [D]. Harbin:

    Northeast Agricultural University, 2020.

    [7] 趙慶雷,信彩云,王瑜,等.不同輪作模式對花生病蟲害及產(chǎn)

    量的影響[J].植物保護(hù)學(xué)報,2018,45(6):1321-1327.

    ZHAO Q L, XIN C Y, WANG Y, et al .. Effects of different

    rotation patterns on peanut diseases, pests and yield [J]. J.

    Plant Prot., 2018, 45(6):1321-1327.

    [8] 徐新雯,林正全,拓陽陽,等.煙蒜輪作對煙株根際土壤細(xì)菌

    群落結(jié)構(gòu)的影響[J]. 西南農(nóng)業(yè)學(xué)報,2020,33(9):1917-1924,

    2137.

    XU X W, LIN Z Q, TUO Y Y, et al .. Effects of tobacco garlic

    crop rotation on bacterial community structure of tobacco

    rhizosphere soil [J]. Southwest China J. Agric. Sci., 2020, 33(9):

    1917-1924,2137.

    [9] 陽顯斌,李廷軒,張錫洲,等.煙蒜輪作與套作對土壤微生物

    類群數(shù)量的影響[J].土壤,2016,48(4):698-704.

    YANG X B, LI T X, ZHANG X Z, et al .. Effects of tobacco

    garlic crop rotation and tobacco garlic crop intercropping on

    soil microbial groups in tobacco fields [J]. Soils, 2016, 48(4):

    698-704.

    [10] 蘇燕,李婕,曹雪穎,等.水旱輪作模式下馬鈴薯根際土壤細(xì)

    菌群落多樣性分析[J]. 南方農(nóng)業(yè)學(xué)報,2020,51(10):2374-

    2382.

    SU Y, LI J, CAO X Y, et al .. Diversity analysis of bacterial

    community in potato rhizosphere soil under the mode of paddyupland

    rotation [J]. J. Southern Agric., 2020, 51(10): 2374-

    2382.

    [11] 姚小東,李孝剛,丁昌峰,等.連作和輪作模式下花生土壤微

    生物群落不同微域分布特征[J]. 土壤學(xué)報,2019,56(4):

    975-985.

    YAO X D, LI X G, DING C F, et al .. Microzone distribution

    characteristics of soil microbial community with peanut

    cropping system, monocropping or rotation [J]. Acta Pedol.

    Sin., 2019, 56(4):975-985.

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

    1-495.

    [13] 李潤根,曾慧蘭,李興杰,等.連作龍牙百合與鐵炮百合根際

    土壤真菌群落結(jié)構(gòu)的差異分析[J]. 生態(tài)科學(xué),2022,41(4):

    189-195.

    LI R G, ZENG H L, LI X J, et al .. Differences of rhizosphere

    soil fungi between Lilium brownii var and Lilium longifllorum

    in continuous cropping [J]. Acta Ecol. Sci., 2022, 41(4):

    189-195.

    [14] 李林蓉,馮建路,劉苗苗,等.作物種植模式對土壤微生物和

    農(nóng)田有害生物的影響[J].中國農(nóng)學(xué)通報,2021,37(29):99-106.

    LI L R, FENG J L, LIU M M, et al .. Effect of crop planting

    Patterns on soil microorganisms and crop pests in farm [J].

    Chin. Agric. Sci. Bull., 2021, 37(29):99-106.

    [15] 于春雷,高嵩,孫文松,等.連作對遼細(xì)辛土壤理化性質(zhì)和根

    際微生物群落特征的影響[J]. 江蘇農(nóng)業(yè)科學(xué),2022,50(14):

    250-258.

    [16] 劉會芳,韓宏偉,王強,等.不同蔬菜與番茄輪作對設(shè)施土壤

    微生物多樣性、酶活性及土壤理化性質(zhì)的影響[J].微生物

    學(xué)報,2021,61(1):167-182.

    LIU H F, HAN H W, WANG Q, et al .. Effect of vegetablestomato

    rotation on soil microbial diversity, enzyme activity and

    physicochemical properties of vegetables in greenhouse [J].

    Acta Microbiol. Sin., 2021, 61(1):167-182.

    [17] SHI G, FANG J, WEI S, et al .. Characteristics of rhizosphere

    fungal community in spring wheat under different rotation

    fallow modes [J/OL]. E3S Web Conf, 2021, 292:03093. [2023-

    07-12]. https://doi.org/10.1051/e3sconf/202129203093.

    [18] SU Y, ZI H, WEI X, et al .. Application of manure rather than

    plant-origin organic fertilizers alters the fungal community in

    continuous cropping tobacco soil [J/OL]. Front Microbiol.,

    2022:1158 [2022-08-23]. http://doi.org/10.3389/fmicb.2022.

    1158.

    [19] 謝濤.不同栽培模式下稻麥輪作土壤速效養(yǎng)分供應(yīng)特征[D].

    揚州:揚州大學(xué),2022.

    XIE T. Characteristics of soil available nutrient supply in ricewheat

    rotation under different cultivation modes [D].

    Yangzhou: Yangzhou University, 2022.

    [20] 張成君,師尚禮,康文娟,等.不同輪作模式土壤酶活性特征

    及與化學(xué)性質(zhì)的關(guān)系[J].中國草地學(xué)報,2020,42(5):92-102.

    ZHANG C J, SHI S L, KANG W J, et al .. Characteristics of soil

    enzyme activities and its relationship with chemical properties

    under different rotation patterns [J]. Chin. J. Grass., 2020,

    42(5):92-102.

    [21] 魏全全,芶久蘭,趙歡,等.黃壤區(qū)烤煙輪作與連作根系形態(tài)、

    產(chǎn)量及養(yǎng)分吸收的變化[J].西南農(nóng)業(yè)學(xué)報,2018,31(11):2294-

    2299.

    WEI Q Q, GOU J L, ZHAO H, et al .. Changes in root

    morphology, yield and nutrient uptake of flue-cured tobacco

    under rotation and continuous cropping in yellow soil [J].

    Southwest China J. Agric. Sci., 2018, 31(11):2294-2299.

    [22] 樊俊,譚軍,王瑞,等.秸稈還田和腐熟有機肥對植煙土壤養(yǎng)

    分、酶活性及微生物多樣性的影響[J]. 煙草科技,2019,

    52(2):12-18,61.

    FAN J, TAN J, WANG R, et al .. Effects of straw Returning and

    decomposed organic manure on soil nutrients, enzyme

    activities and microbial diversity for tobacco-planting [J]. Tob.

    Sci. Tech., 2019, 52(2):12-18,61.

    [23] 羅影,王立光,陳軍,等.不同種植模式對甘肅中部高寒區(qū)胡

    麻田土壤酶活性及土壤養(yǎng)分的影響[J]. 核農(nóng)學(xué)報,2017,

    31(6):1185-1191.

    LUO Y, WANG L G, CHEN J, et al .. Effects of different flex

    cropping modes on soil enzyme activities and soil nutrients in

    the cold area of middle part of Gansu [J]. J. Nucl. Agric. Sci.,

    2017, 31(6):1185-1191.

    [24] RAMONEDA J, LE ROUX J, STADELMANN S, et al .. Soil

    microbial community coalescence and fertilization interact to

    drive the functioning of the legume-rhizobium symbiosis [J]. J.

    Appl. Ecol., 2021, 58(11):2590-2602.

    [25] 索炎炎,張翔,司賢宗,等.AM真菌和根瘤菌對連作花生養(yǎng)分

    吸收及土壤微生物特性的影響[J/OL]. 中國土壤與肥料,

    2023(2):106-112.

    SUO Y Y, ZHANG X, SI X Z, et al .. Effects of Arbuscular

    mycorrhizal fungi and rhizobia on nutrient uptake and soil

    microbial characteristics of continuous cropping peanut [J/OL].

    Soil Fert. Sci. China, 2023(2):106-112.

    [26] 饒德安,劉潘洋,鄒路易,等.長期連作及強還原土壤滅菌處

    理對烤煙根際土壤真菌群落的影響[J]. 中國土壤與肥料,

    2022(4):47-56.

    RAO D A, LIU P Y, ZOU L Y, et al .. Effects of long-term

    continuous cropping and reductive soil disinfestation on fungal

    community in flue-cured tobacco rhizosphere [J]. Soil Fert. Sci.

    China, 2022(4):47-56.

    [27] 段玉琪,晉艷,陳澤斌,等.烤煙輪作與連作土壤細(xì)菌群落多

    樣性比較[J].中國煙草學(xué)報,2012,18(6):53-59.

    DUAN Y Q, JIN Y, CHEN Z B, et al . Comparison of bacteria

    diversity between tobacco plantation soils of rotational

    cropping and continuous cropping [J]. Acta Tab. Sin., 2012,

    18(6):53-59.

    [28] 蘇貝貝,張英,道日娜.4種豆科植物根際土壤真菌群落特征

    與土壤理化因子間相關(guān)性分析[J]. 草地學(xué)報,2021,29(12):

    2670-2677.

    SU B B, ZHAN Y, DAO R N. Correlation analysis between

    fungal community characteristics and soil physicochemical

    factors in the rhizosphere of four legumes [J]. Acta Agrestia

    Sin., 2021, 29(12):2670-2677.

    [29] 溫美娟,楊思存,王成寶,等.深松和秸稈還田對灌耕灰鈣土

    土壤細(xì)菌多樣性和群落結(jié)構(gòu)的影響[J].農(nóng)業(yè)資源與環(huán)境學(xué)

    報,2023,40(2):423-433.

    WEN M J, YANG S C, WANG C B, et al .. Effects of subsoiling

    and the return of straw on soil bacterial diversity and

    community structure in an irrigated sierozem farmland [J]. J.

    Agric. Resour. Environ., 2023, 40(2):423-433.

    [30] 葉嘉,張浩,郭海燕,等.蚓糞對大蒜中蒜氨酸和大蒜素含量

    及大蒜精油抗菌活性的影響[J]. 中國瓜菜,2021,34(10):

    98-103.

    YE J, ZHANG H, GUO H Y, et al .. Effects of vermicompost on

    alliin and allicin contents and antibacterial activity of garlic oil

    in garlic [J]. China Cucurbits Veget., 2021, 34(10):98-103.

    (責(zé)任編輯:張冬玲)

    猜你喜歡
    植煙土壤微生物群落土壤酶活性
    微生物群落變化與妊娠結(jié)局關(guān)系的研究進(jìn)展
    貴州省植煙土壤養(yǎng)分含量狀況分析
    土壤調(diào)理劑對植煙土壤微生物群落的影響
    海岸帶森林生態(tài)系統(tǒng)中土壤酶活性研究進(jìn)展
    衡陽煙區(qū)植煙土壤保育技術(shù)研究與應(yīng)用
    高碳基土壤修復(fù)肥對植煙土壤有效微量元素及煙葉品質(zhì)的影響
    重金屬污染對土壤微生物及土壤酶活性影響的研究進(jìn)展
    精河沙區(qū)土壤酶分布特征及其對土壤理化性狀的響應(yīng)
    刈割對不同植物群落特征、土壤微生物和土壤酶活性的影響
    間作玉米對馬鈴薯連作土壤微生物群落多樣性的調(diào)節(jié)作用
    免费观看a级毛片全部| 法律面前人人平等表现在哪些方面 | 亚洲精品国产av蜜桃| 两性夫妻黄色片| 黑人操中国人逼视频| 少妇裸体淫交视频免费看高清 | 人妻人人澡人人爽人人| 国产伦理片在线播放av一区| 久久精品亚洲熟妇少妇任你| 国产一区二区 视频在线| 黑人欧美特级aaaaaa片| 老司机午夜福利在线观看视频 | 人妻一区二区av| 精品人妻1区二区| 精品少妇一区二区三区视频日本电影| 成人国产av品久久久| 岛国在线观看网站| 成年人免费黄色播放视频| 在线十欧美十亚洲十日本专区| 国产亚洲av高清不卡| 午夜两性在线视频| 性色av一级| 日韩一区二区三区影片| 亚洲精品中文字幕在线视频| 亚洲全国av大片| 纯流量卡能插随身wifi吗| 91字幕亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 最新的欧美精品一区二区| 超色免费av| 如日韩欧美国产精品一区二区三区| 纵有疾风起免费观看全集完整版| 国产欧美日韩精品亚洲av| www.999成人在线观看| 午夜免费观看性视频| 搡老岳熟女国产| 天天躁夜夜躁狠狠躁躁| 丰满饥渴人妻一区二区三| a级毛片黄视频| 91九色精品人成在线观看| 黑人猛操日本美女一级片| 成人影院久久| 亚洲精品国产一区二区精华液| 久9热在线精品视频| 欧美亚洲 丝袜 人妻 在线| 男女免费视频国产| 少妇粗大呻吟视频| av国产精品久久久久影院| 国产一区二区在线观看av| 欧美xxⅹ黑人| 亚洲性夜色夜夜综合| 91精品三级在线观看| 国产亚洲精品久久久久5区| 伊人亚洲综合成人网| 黄片大片在线免费观看| 久久久国产成人免费| 国产成人免费观看mmmm| 免费久久久久久久精品成人欧美视频| 色视频在线一区二区三区| 亚洲av电影在线观看一区二区三区| 精品一区在线观看国产| svipshipincom国产片| 国产成人精品在线电影| 久久久精品免费免费高清| 亚洲少妇的诱惑av| 国产精品香港三级国产av潘金莲| 极品少妇高潮喷水抽搐| 午夜免费观看性视频| 人人妻人人澡人人爽人人夜夜| 中亚洲国语对白在线视频| 制服诱惑二区| 午夜激情久久久久久久| 两个人看的免费小视频| 热99国产精品久久久久久7| 精品人妻一区二区三区麻豆| 天天躁夜夜躁狠狠躁躁| 欧美av亚洲av综合av国产av| av有码第一页| 高清黄色对白视频在线免费看| 精品熟女少妇八av免费久了| 亚洲免费av在线视频| 黄频高清免费视频| 最近最新免费中文字幕在线| 女警被强在线播放| 国产片内射在线| 午夜福利视频精品| 国产精品1区2区在线观看. | 国内毛片毛片毛片毛片毛片| 精品福利永久在线观看| 免费不卡黄色视频| 久久香蕉激情| 亚洲黑人精品在线| 亚洲精品美女久久久久99蜜臀| 免费在线观看完整版高清| 91精品三级在线观看| 免费观看av网站的网址| 中国美女看黄片| 99热网站在线观看| 两个人看的免费小视频| 女性被躁到高潮视频| 纯流量卡能插随身wifi吗| 在线永久观看黄色视频| 天堂中文最新版在线下载| 性少妇av在线| 亚洲欧美一区二区三区久久| 一二三四在线观看免费中文在| 无遮挡黄片免费观看| 狠狠婷婷综合久久久久久88av| 午夜激情久久久久久久| 欧美亚洲日本最大视频资源| www.av在线官网国产| 麻豆av在线久日| 两个人看的免费小视频| 999久久久国产精品视频| 天天添夜夜摸| 免费在线观看黄色视频的| 亚洲成人手机| 人妻 亚洲 视频| 视频在线观看一区二区三区| 久久国产亚洲av麻豆专区| av在线app专区| 丝袜美腿诱惑在线| 制服人妻中文乱码| 手机成人av网站| 日韩欧美一区视频在线观看| 国产色视频综合| 精品国内亚洲2022精品成人 | 黄色片一级片一级黄色片| 国精品久久久久久国模美| 天天躁夜夜躁狠狠躁躁| 免费观看人在逋| 一二三四在线观看免费中文在| 久久亚洲精品不卡| 欧美日韩福利视频一区二区| 国产老妇伦熟女老妇高清| 午夜福利一区二区在线看| 免费人妻精品一区二区三区视频| 久久久久视频综合| 亚洲自偷自拍图片 自拍| 男女无遮挡免费网站观看| 超色免费av| 国产99久久九九免费精品| 日韩欧美国产一区二区入口| 精品国产一区二区三区久久久樱花| 女警被强在线播放| 亚洲午夜精品一区,二区,三区| 蜜桃在线观看..| 男人操女人黄网站| 午夜日韩欧美国产| 精品一品国产午夜福利视频| 99精品久久久久人妻精品| 日日爽夜夜爽网站| 99国产极品粉嫩在线观看| 久久久精品免费免费高清| 亚洲五月婷婷丁香| 久久精品国产综合久久久| 美女午夜性视频免费| 一区二区三区精品91| av国产精品久久久久影院| av天堂在线播放| 国产免费福利视频在线观看| 午夜免费鲁丝| 婷婷色av中文字幕| 在线看a的网站| 亚洲综合色网址| 激情视频va一区二区三区| 国产成人精品在线电影| 成人影院久久| 亚洲伊人色综图| 51午夜福利影视在线观看| 中文字幕色久视频| 亚洲情色 制服丝袜| 国产精品一区二区在线观看99| 欧美激情 高清一区二区三区| 色94色欧美一区二区| 国产色视频综合| 国产高清videossex| 韩国高清视频一区二区三区| 国产欧美日韩一区二区三 | 色婷婷av一区二区三区视频| 夫妻午夜视频| 久久久久国内视频| 欧美日韩福利视频一区二区| 99热网站在线观看| 国产精品亚洲av一区麻豆| 欧美少妇被猛烈插入视频| 免费日韩欧美在线观看| 在线观看免费视频网站a站| 老司机影院毛片| 777米奇影视久久| 9191精品国产免费久久| 亚洲专区中文字幕在线| 成人亚洲精品一区在线观看| 亚洲色图 男人天堂 中文字幕| 色综合欧美亚洲国产小说| 亚洲中文日韩欧美视频| 999久久久国产精品视频| 丝袜美腿诱惑在线| 91麻豆av在线| 亚洲精品成人av观看孕妇| 国产又色又爽无遮挡免| 老司机影院毛片| 欧美黄色淫秽网站| 欧美另类一区| 水蜜桃什么品种好| 精品一品国产午夜福利视频| 免费观看av网站的网址| 性高湖久久久久久久久免费观看| 免费观看a级毛片全部| 人人妻人人添人人爽欧美一区卜| 欧美亚洲日本最大视频资源| 欧美精品高潮呻吟av久久| 悠悠久久av| av在线app专区| 国产精品 国内视频| 9色porny在线观看| 韩国精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 国产亚洲午夜精品一区二区久久| 色老头精品视频在线观看| 亚洲精品国产区一区二| videosex国产| 18禁裸乳无遮挡动漫免费视频| 18在线观看网站| 久久青草综合色| 亚洲伊人色综图| 伊人亚洲综合成人网| 成人影院久久| 超碰97精品在线观看| 日本黄色日本黄色录像| 久久久久国内视频| 一级毛片电影观看| 精品久久久精品久久久| 亚洲avbb在线观看| 久久国产精品大桥未久av| 欧美精品人与动牲交sv欧美| 日韩一卡2卡3卡4卡2021年| 五月天丁香电影| 深夜精品福利| 超色免费av| 国产一区二区 视频在线| 热re99久久精品国产66热6| 久久久欧美国产精品| 免费观看人在逋| 精品久久久久久电影网| 欧美大码av| 亚洲欧美激情在线| 日本91视频免费播放| videosex国产| 一区二区av电影网| 色婷婷久久久亚洲欧美| 一区福利在线观看| 肉色欧美久久久久久久蜜桃| 狂野欧美激情性xxxx| 久久热在线av| 亚洲全国av大片| 欧美日本中文国产一区发布| 交换朋友夫妻互换小说| 日韩大片免费观看网站| 亚洲一区中文字幕在线| 最新的欧美精品一区二区| 国产男女内射视频| 亚洲五月色婷婷综合| 一区二区三区激情视频| 91成人精品电影| 欧美老熟妇乱子伦牲交| 视频区欧美日本亚洲| 欧美成狂野欧美在线观看| 欧美国产精品一级二级三级| 秋霞在线观看毛片| 午夜福利视频在线观看免费| 国产亚洲精品久久久久5区| 18禁裸乳无遮挡动漫免费视频| 老司机深夜福利视频在线观看 | 亚洲专区中文字幕在线| 黄色 视频免费看| 国产欧美亚洲国产| 99久久国产精品久久久| 大片电影免费在线观看免费| a 毛片基地| 精品亚洲成a人片在线观看| 50天的宝宝边吃奶边哭怎么回事| 一级片'在线观看视频| 无限看片的www在线观看| 欧美精品高潮呻吟av久久| 中亚洲国语对白在线视频| 国产精品免费大片| 一本大道久久a久久精品| 丰满人妻熟妇乱又伦精品不卡| 欧美另类一区| 激情视频va一区二区三区| 日韩熟女老妇一区二区性免费视频| 欧美日韩亚洲国产一区二区在线观看 | 精品人妻在线不人妻| 国产亚洲欧美在线一区二区| 菩萨蛮人人尽说江南好唐韦庄| 欧美激情久久久久久爽电影 | 午夜视频精品福利| 亚洲视频免费观看视频| 波多野结衣一区麻豆| 国产99久久九九免费精品| 日韩一区二区三区影片| 黄片播放在线免费| 一级毛片女人18水好多| 亚洲专区国产一区二区| 男女床上黄色一级片免费看| 欧美日韩福利视频一区二区| 久久精品亚洲熟妇少妇任你| 国产精品99久久99久久久不卡| 18禁观看日本| 日韩人妻精品一区2区三区| 欧美另类一区| 老汉色∧v一级毛片| 免费观看人在逋| 天堂8中文在线网| 日本91视频免费播放| 视频在线观看一区二区三区| 一级毛片精品| 少妇粗大呻吟视频| 另类亚洲欧美激情| 久久久精品国产亚洲av高清涩受| 亚洲精品第二区| 国产欧美日韩一区二区三 | 日本vs欧美在线观看视频| 狂野欧美激情性bbbbbb| 99精品欧美一区二区三区四区| 成年女人毛片免费观看观看9 | 成在线人永久免费视频| e午夜精品久久久久久久| 欧美激情久久久久久爽电影 | 欧美日韩黄片免| 欧美日韩精品网址| 午夜福利,免费看| 亚洲avbb在线观看| 亚洲国产欧美日韩在线播放| 国产有黄有色有爽视频| 一区二区三区四区激情视频| 男人添女人高潮全过程视频| 亚洲美女黄色视频免费看| 在线观看舔阴道视频| 国产免费av片在线观看野外av| 国产色视频综合| 在线观看舔阴道视频| a级片在线免费高清观看视频| 欧美成狂野欧美在线观看| 亚洲专区字幕在线| 国产一区二区三区综合在线观看| 中文精品一卡2卡3卡4更新| 飞空精品影院首页| 日韩熟女老妇一区二区性免费视频| 日韩一卡2卡3卡4卡2021年| 国产在视频线精品| 水蜜桃什么品种好| 久久久久精品人妻al黑| 一本色道久久久久久精品综合| 麻豆乱淫一区二区| 国产精品av久久久久免费| 91大片在线观看| 精品人妻1区二区| 免费看十八禁软件| av电影中文网址| 欧美激情 高清一区二区三区| e午夜精品久久久久久久| 欧美精品亚洲一区二区| 亚洲精品美女久久av网站| 男女国产视频网站| 18禁观看日本| 中文欧美无线码| 国产精品久久久久成人av| 国产不卡av网站在线观看| 午夜老司机福利片| 国产精品99久久99久久久不卡| 久久久精品免费免费高清| 99国产精品99久久久久| 免费在线观看影片大全网站| 高清黄色对白视频在线免费看| 日本av手机在线免费观看| 日本欧美视频一区| 国产成人影院久久av| 一级片免费观看大全| 日本av手机在线免费观看| 国产在线观看jvid| 亚洲激情五月婷婷啪啪| 又黄又粗又硬又大视频| 国产一区二区激情短视频 | a 毛片基地| 亚洲国产日韩一区二区| 老司机午夜福利在线观看视频 | 国产精品 欧美亚洲| 久久久久久久大尺度免费视频| 日本精品一区二区三区蜜桃| 高清av免费在线| 老司机靠b影院| 最近最新中文字幕大全免费视频| 亚洲精品国产一区二区精华液| 黑人操中国人逼视频| 亚洲成国产人片在线观看| 丰满少妇做爰视频| 18禁国产床啪视频网站| 夜夜夜夜夜久久久久| 女性被躁到高潮视频| 国产淫语在线视频| 两人在一起打扑克的视频| 婷婷成人精品国产| 国产免费福利视频在线观看| 啦啦啦在线免费观看视频4| 久久国产精品男人的天堂亚洲| 日韩人妻精品一区2区三区| 最黄视频免费看| 久久久久国产精品人妻一区二区| 国产av精品麻豆| 黑人操中国人逼视频| 啦啦啦 在线观看视频| 国产精品 欧美亚洲| 啦啦啦免费观看视频1| 最黄视频免费看| 老司机靠b影院| 中国美女看黄片| 成年人午夜在线观看视频| 国产精品自产拍在线观看55亚洲 | 色老头精品视频在线观看| av视频免费观看在线观看| 国产欧美日韩综合在线一区二区| 一二三四在线观看免费中文在| 男女免费视频国产| 欧美乱码精品一区二区三区| a级片在线免费高清观看视频| 亚洲免费av在线视频| 91字幕亚洲| 午夜福利在线观看吧| 午夜福利视频精品| 国产亚洲一区二区精品| 老鸭窝网址在线观看| 高清欧美精品videossex| 少妇被粗大的猛进出69影院| 亚洲欧美成人综合另类久久久| 男人爽女人下面视频在线观看| 国产精品二区激情视频| 老熟女久久久| 亚洲五月色婷婷综合| 人人妻人人爽人人添夜夜欢视频| 多毛熟女@视频| 亚洲全国av大片| 桃红色精品国产亚洲av| 国产亚洲精品一区二区www | av欧美777| 99国产精品一区二区蜜桃av | 中文字幕人妻丝袜一区二区| 国产精品国产av在线观看| 国产av又大| 日韩 亚洲 欧美在线| 国产精品自产拍在线观看55亚洲 | 国产主播在线观看一区二区| 精品国产一区二区久久| 欧美日韩av久久| 肉色欧美久久久久久久蜜桃| 亚洲精品av麻豆狂野| 自线自在国产av| 国产精品久久久久成人av| 搡老岳熟女国产| 精品国产一区二区久久| 少妇粗大呻吟视频| 18在线观看网站| 咕卡用的链子| 成年人免费黄色播放视频| 免费av中文字幕在线| 极品少妇高潮喷水抽搐| 99香蕉大伊视频| 亚洲专区国产一区二区| 国产精品免费大片| 午夜精品久久久久久毛片777| 狂野欧美激情性bbbbbb| 亚洲黑人精品在线| 手机成人av网站| 国产成人精品无人区| 日韩一卡2卡3卡4卡2021年| 久久久久久久国产电影| 亚洲人成电影免费在线| 老司机靠b影院| 男人操女人黄网站| 国产欧美日韩一区二区三区在线| 成人手机av| 另类精品久久| 中文欧美无线码| 天堂俺去俺来也www色官网| 最近最新免费中文字幕在线| 成年人午夜在线观看视频| 亚洲av日韩在线播放| www日本在线高清视频| 男女免费视频国产| 欧美人与性动交α欧美软件| 久久人妻福利社区极品人妻图片| 色播在线永久视频| videosex国产| 精品卡一卡二卡四卡免费| 一个人免费在线观看的高清视频 | 蜜桃国产av成人99| 69精品国产乱码久久久| 中文字幕制服av| 亚洲天堂av无毛| 午夜精品国产一区二区电影| 欧美精品人与动牲交sv欧美| 777久久人妻少妇嫩草av网站| 国产精品 国内视频| 日本av手机在线免费观看| 欧美在线黄色| 久久久国产一区二区| 两个人免费观看高清视频| 亚洲欧洲精品一区二区精品久久久| 国产一卡二卡三卡精品| 男女国产视频网站| 欧美老熟妇乱子伦牲交| 亚洲精品久久午夜乱码| av有码第一页| 王馨瑶露胸无遮挡在线观看| 黄色视频,在线免费观看| 99国产精品一区二区三区| 狠狠狠狠99中文字幕| 精品福利观看| 最近最新免费中文字幕在线| 天堂俺去俺来也www色官网| av福利片在线| 精品少妇一区二区三区视频日本电影| 久久性视频一级片| 日韩中文字幕欧美一区二区| 午夜精品久久久久久毛片777| 男人舔女人的私密视频| 日本撒尿小便嘘嘘汇集6| 亚洲第一av免费看| 国产成人啪精品午夜网站| 亚洲,欧美精品.| 成人国产av品久久久| 国产精品一区二区精品视频观看| 亚洲国产精品一区二区三区在线| 国产精品一区二区精品视频观看| 这个男人来自地球电影免费观看| 午夜日韩欧美国产| 色94色欧美一区二区| 亚洲欧美清纯卡通| 母亲3免费完整高清在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 两个人看的免费小视频| 99国产精品免费福利视频| 99久久综合免费| 国产成人精品无人区| 亚洲五月色婷婷综合| 欧美人与性动交α欧美精品济南到| 欧美 日韩 精品 国产| 一本一本久久a久久精品综合妖精| 青春草亚洲视频在线观看| 国产欧美日韩精品亚洲av| 99久久99久久久精品蜜桃| 久久人妻熟女aⅴ| 精品久久久精品久久久| 亚洲性夜色夜夜综合| 男人爽女人下面视频在线观看| 啦啦啦 在线观看视频| 中文字幕精品免费在线观看视频| 大型av网站在线播放| 老司机在亚洲福利影院| 欧美成狂野欧美在线观看| 成年av动漫网址| 狂野欧美激情性xxxx| 久久毛片免费看一区二区三区| 亚洲国产精品一区三区| 日韩三级视频一区二区三区| 亚洲中文日韩欧美视频| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 成人亚洲精品一区在线观看| 精品国产国语对白av| 欧美激情高清一区二区三区| 成年人午夜在线观看视频| 亚洲av欧美aⅴ国产| 国产日韩欧美视频二区| 成人黄色视频免费在线看| 老熟女久久久| 午夜激情久久久久久久| 国产又色又爽无遮挡免| 青青草视频在线视频观看| 国产精品一区二区在线观看99| 女性被躁到高潮视频| 狠狠狠狠99中文字幕| 777久久人妻少妇嫩草av网站| a级毛片在线看网站| 一二三四在线观看免费中文在| 成年人黄色毛片网站| 12—13女人毛片做爰片一| 久久这里只有精品19| 国产又爽黄色视频| 成人18禁高潮啪啪吃奶动态图| 久久天堂一区二区三区四区| 免费av中文字幕在线| 人妻一区二区av| 国产1区2区3区精品| 国内毛片毛片毛片毛片毛片| 亚洲一区中文字幕在线| 亚洲国产欧美一区二区综合| 一区二区三区激情视频| 波多野结衣一区麻豆| 精品亚洲乱码少妇综合久久| 欧美亚洲 丝袜 人妻 在线| 日本猛色少妇xxxxx猛交久久| 国产主播在线观看一区二区| av天堂久久9| 国产在线免费精品| 又黄又粗又硬又大视频| 狠狠狠狠99中文字幕| 亚洲欧美日韩另类电影网站| 欧美日韩成人在线一区二区| 日韩中文字幕欧美一区二区| 欧美日韩视频精品一区| 欧美日韩亚洲高清精品| 最近中文字幕2019免费版| av国产精品久久久久影院|