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

    連作木薯對(duì)土壤理化性質(zhì)及細(xì)菌群落演替的影響

    2024-08-22 00:00:00彭曉輝黎亮武歐桂寧黃苑航彭曉雪陽太億甘李申章佑韋茂貴
    關(guān)鍵詞:根際土壤木薯

    摘要:【目的】研究連作木薯對(duì)非根際和根際土壤理化性質(zhì)變化及細(xì)菌群落演替的影響,為系統(tǒng)闡明木薯連作障礙形成機(jī)制及其緩解措施研發(fā)提供理論依據(jù)?!痉椒ā恳阅臼砥贩N華南205為試驗(yàn)材料,在撂荒地上連續(xù)種植3年,每年木薯塊根成熟期(植后250d)采集土壤樣品。結(jié)合土壤農(nóng)化分析、高通量測序技術(shù)和生物信息學(xué)手段,研究連作木薯對(duì)非根際和根際土壤理化性質(zhì)及細(xì)菌群演替的影響?!窘Y(jié)果】與種植第1年木薯相比,連作2年和3年的木薯非根際和根際土壤pH及有機(jī)質(zhì)、堿解氮和速效磷含量降低,脲酶活性提高,除速效磷外,其他指標(biāo)的變化均達(dá)顯著水平(Plt;0.05,下同)。連作木薯會(huì)引起土壤細(xì)菌群落多樣性及豐富度降低,其中根際土壤細(xì)菌群落多樣性及豐富度下降速度比非根際土壤慢。連作木薯非根際與根際土壤細(xì)菌群落的分類和系統(tǒng)發(fā)育Beta多樣性存在顯著差異。連作木薯非根際和根際土壤的優(yōu)勢菌門均為放線菌門(Actinobacteriota)、變形菌門(Proteobacteria)、綠彎菌門(Chloroflexi)、厚壁菌門(Firmicutes)、酸桿菌門(Acidobacteriota),主要優(yōu)勢菌屬均為芽孢桿菌屬(Bacillus)、norank f Roseiflexaceae、norank_f_norank_o_Vicinamibacterales、norank_f_norank_o_Gaiellales。在操作分類單元(OTU)水平,不同連作年限木薯的非根際和根際土壤細(xì)菌群落組成存在顯著或極顯著(Plt;0.01)差異。土壤的細(xì)菌群落結(jié)構(gòu)與理化性質(zhì)相互影響,對(duì)細(xì)菌群落影響較大的環(huán)境因子排序?yàn)閴A解氮含量gt;有機(jī)質(zhì)含量gt;pHgt;速效磷含量gt;脲酶活性。

    【結(jié)論】連作會(huì)導(dǎo)致木薯非根際和根際土壤理化性質(zhì)、細(xì)菌群落多樣性和豐富度發(fā)生變化,土壤理化性質(zhì)和細(xì)菌群落結(jié)構(gòu)均在連作第3年出現(xiàn)明顯劣化。木薯土壤理化性質(zhì)變化與土壤細(xì)菌群落結(jié)構(gòu)差異存在一定相關(guān)性,有益菌豐度降低可能是木薯連作造成土壤細(xì)菌群落劣化演替的重要原因之一。

    關(guān)鍵詞:木薯;連作;非根際土壤;根際土壤;細(xì)菌群落

    中圖分類號(hào):S533

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):2095-1191(2024)04-0942-12

    Effects of continuous cropping of cassava on soil physicochemi-cal properties and bacterial community succession

    PENG Xiao-hui1, LI Liang-wu1, OU Gui-ning1, HUANG Yuan-hang1, PENG Xiao-xue1,YANG Tai-yi1, GAN Li1, SHEN Zhang-you2*, WEI Mao-gui1,3*

    (1College of Agriculture, Guangxi University, Nanning, Guangxi 530004, China; 2Cash Crops Research Institute, Guangxi Academy of Agricultural Sciences,Nanning, Guangxi 530007, China; 3Guangxi Key Laboratory of Agro-environment and Agro-product Safety, Nanning, Guangxi 530004, China)

    Abstract:[Objective]The aim of this study was to investigate how continuous cropping of cassava affected the physicochemical properties of non-rhizosphere and rhizosphere soils, as well as the succession of bacterial communities, and thus providing a theoretical basis for understanding the formation of cassava continuous cropping obstacles. 【Method】The cassava variety South China 205 was continuously cultivated in an abandoned field for three years. The soil were sampled at cassava root tuber maturity stage (250 d after planting) in each year. Combining soil agrochemical analysis, high-throughput sequencing technology and bioinformatics methods, effects of continuous cassava on non-rhizosphere and rhizosphere soil physicochemical properties and bacterial group succession were studied. 【Result] The results showed that pH, organic matter, alkali-hydrolyzed nitrogen and available phosphorus contents of non-rhizosphere and rhizosphere soil of cassava were decreased and urease activity was increased in two and three years of continuous cropping compared with the first year of planting, and except available phosphorus content, the changes of other indexes reached significant level (Plt;0.05, the same below). Continuous cropping of cassava led to a decrease in the diversity and richness of soil bacterial communities. However, the decline rate in rhizosphere soil bacterial communities was slower than that in non-rhizosphere soils. There were significant differences in the classification and phylogenetic Beta diversity of non-rhizosphere and rhizo- sphere soil bacterial communities of continuous cropping cassava. The dominant bacterial phyla in non-rhizosphere and rhi- zosphere soils were Actinobacteriota, Proteobacteria, Chloroflexi, Firmicutes and Acidobacteriota. The main dominant bac- terial genera were Bacillus, norank f_Roseiflexaceae, norank f_norank_o_Vicinamibacterales and norank _f_norank_o_ Gaiellales. At the operational taxonomic unit (OTU) level, significant or extremely significant(Plt;0.01) differences were found in the composition of soil bacterial communities in non-rhizosphere and rhizosphere soils of cassava at different con- tinuous cropping years. The structure of soil bacterial communities was mutually affected by soil physicochemical proper- ties, the environmental factors affected bacterial community ranked as alkaline hydrolysis nitrogen content gt; organic mat- ter content gt; pH gt; available phosphorus content gt; urease activity. 【Conclusion ]Continuous cropping can lead to changes in non-rhizosphere and rhizosphere soil physicochemical properties, bacterial community diversity and richness, and soil physicochemical properties and bacterial community structure of cassava greatly deteriorate in the third year of continuous cropping. Soil physicochemical properties are correlated to structure of bacterial communities at certain extent. Notably, the decrease in the abundance of beneficial bacteria might be an important reason for the deterioration of soil bacterial community succession caused by cassava continuous cropping.

    Key words: cassava; continuous cropping; non-rhizosphere soil; rhizosphere soil; bacterial community

    Foundation items: National Natural Science Foundation of China (31960389,31860347)

    0 引言

    【研究意義】隨著人口的不斷增加和城鎮(zhèn)化進(jìn)程的加快,我國耕地資源日益減少。土地資源高度集約化規(guī)模管理、提高土地和肥料利用率已成為現(xiàn)代農(nóng)業(yè)發(fā)展的趨勢。現(xiàn)代農(nóng)業(yè)逐漸形成了復(fù)種指數(shù)高、作物規(guī)模化種植且種類單一的特點(diǎn)。受土地資源、氣候條件、種植習(xí)慣及經(jīng)濟(jì)效益等因素影響,同種或同科作物連作不可避免,許多作物存在著不同程度的連作障礙,已成為影響作物產(chǎn)量提高的關(guān)鍵問題之一(翁佩瑩和鄭紅艷,2020)。木薯是我國重要的熱帶作物之一,屬高光效C3作物,其理論產(chǎn)量最高可達(dá)120t/ha,大田種植高產(chǎn)記錄也高達(dá)80t/haa但實(shí)際生產(chǎn)中因木薯具有耐旱耐貧瘠特性,我國木薯大多連年種植在南方相對(duì)貧瘠的邊際土地上,廣種薄收,連作現(xiàn)象普遍,部分地塊連作年限長達(dá)15年以上,塊根平均產(chǎn)量僅16~17 t/ha(梁海波等,2017)。連作障礙已被證實(shí)是我國木薯單產(chǎn)顯著下降的重要因素之一(王戰(zhàn),2016;梁海波等,2017)。此外,為保證產(chǎn)量,化肥和農(nóng)藥使用量隨木薯連作年限的增加而逐年增加,進(jìn)而造成土壤酸化、重金屬污染和水體污染等農(nóng)業(yè)面源污染問題。因此,研究木薯連作障礙及其緩解措施,對(duì)提高木薯產(chǎn)量、減少化肥施用和保護(hù)生態(tài)環(huán)境均具有重要意義?!厩叭搜芯窟M(jìn)展】連作障礙的形成是土壤、植物和微生物綜合作用的結(jié)果,其形成機(jī)理相對(duì)復(fù)雜。根據(jù)前人研究結(jié)果,引起連作障礙的主要原因可歸納為作物化感物質(zhì)的自毒作用、土壤理化性質(zhì)劣化及土壤生物學(xué)環(huán) 境惡化(侯乾等,2019)。酚酸類、黃酮類和生物堿類 物質(zhì)是作物的主要自毒物質(zhì)種類(陳?;鄣?,2022)。 作物的化感物質(zhì)主要通過影響植株的細(xì)胞膜透性、 光合作用、呼吸作用、有絲分裂、酶活性、蛋白質(zhì)合成 等制約植物生長,進(jìn)而造成作物產(chǎn)量及品質(zhì)降低(楊 浩娜等,2022)。連作還會(huì)改變土壤理化性質(zhì),如連 作木薯使土壤容重變大,三相比例劣化,有機(jī)質(zhì)、堿 解氮、速效鉀和速效磷含量下降,根際土壤總酚酸含 量增加,土壤酸化嚴(yán)重(周貴靖,2017);連作大豆能 提高土壤中速效氮、速效磷和速效鉀含量(劉株秀 等,2019);連作靈芝時(shí)土壤pH、有機(jī)質(zhì)、速效鉀、速 效磷和堿解氮含量均下降(盧孟召等,2022)。土壤 理化性質(zhì)的改變會(huì)影響土壤微生物的豐富度(Hon- toria et al.,2019)。長年連作木薯的土壤中細(xì)菌種群 數(shù)量減少,真菌種群數(shù)量增加,根際微生物群落向真 菌型轉(zhuǎn)化(劉珊廷,2020);蘋果連作后土壤中的細(xì)菌 多樣性減少,真菌多樣性增加(王曉琪等,2020);連 作大豆和玉米均會(huì)顯著降低土壤中真菌的豐富度和 多樣性(宋秀麗等,2022)。此外,土壤微生物包含有 益菌群和有害菌群,有益菌群對(duì)作物生長發(fā)育起促 進(jìn)作用,而有害菌群會(huì)影響作物正常的生長發(fā)育(孫 子欣和蔡柏巖,2022)。由此可見,不同作物連作會(huì) 對(duì)土壤理化性質(zhì)及土壤微生物群落產(chǎn)生不同影響。 【本研究切入點(diǎn)】目前,關(guān)于木薯連作障礙的研究主 要集中在土壤理化性質(zhì)方面(周貴靖,2017),也有少量連作對(duì)土壤微生物影響的報(bào)道(劉珊廷,2020),但 木薯連作對(duì)土壤理化性質(zhì)和微生物群落演替及其相 關(guān)性影響的研究鮮有報(bào)道。【擬解決的關(guān)鍵問題】以 木薯品種華南205為試驗(yàn)材料,在撂荒地上連續(xù)種 植3年,利用高通量測序技術(shù)和生物信息學(xué)手段,研 究連作木薯對(duì)非根際和根際土壤理化性質(zhì)變化及細(xì) 菌群落演替的影響,為系統(tǒng)闡明木薯連作障礙形成 機(jī)制及其緩解措施研發(fā)提供理論依據(jù)。

    1材料與方法

    1.1試驗(yàn)材料

    供試材料為廣西主栽木薯品種華南205。

    1.2試驗(yàn)方法

    試驗(yàn)于2019年4月-2021年12月在廣西南寧 進(jìn)行。試驗(yàn)地撂荒多年,其0~20cm土層pH6.24, 含有機(jī)質(zhì)17.23 g/kg、堿解氮60.38 mg/kg、速效磷 45.09mg/kg、速效鉀42.17mg/kg。常規(guī)耕作地塊,犁地后使用臥式旋耕耙兩輪。田塊劃分成12個(gè)面 積為70m2(14m×5m)的試驗(yàn)小區(qū)網(wǎng)格,第1年隨機(jī) 挑選3個(gè)小區(qū)作為重復(fù),此后2年仍定位原小區(qū)作 為當(dāng)年重復(fù)。木薯分別于2019年4月16日、2020年4月10日和2021年4月11日種植。種莖長度約15cm, 株行距1m×1m,每小區(qū)共種植70株木薯。試驗(yàn)期間無人工灌水,木薯田間管理方式一致。試驗(yàn)期間 降水量和月平均氣溫如圖1所示。

    于木薯塊根成熟期(植后250d)采集土壤樣品。根際土壤:每小區(qū)隨機(jī)選擇3株木薯,挖出完整木薯 塊根和纖維根,采用抖根法去除大塊土壤,用毛刷采 集附著在根系上的土壤。非根際土壤:各小區(qū)按S 形隨機(jī)選取5個(gè)點(diǎn),分別用取土器采集0~20cm土 樣,混勻后作為該小區(qū)的非根際土壤樣品。土壤樣 品置于塑料薄膜上,剔除植物殘?bào)w、石礫等雜物,混 勻后以四分法分為3份。第1份立刻用液氮速凍,置于冰盒中帶回實(shí)驗(yàn)室于-80°C保存,用于土壤微生 物群落多樣性測定;第2份做好標(biāo)記,帶回實(shí)驗(yàn)室于 冰箱4°℃保存,用于土壤酶活性測定;第3份置于標(biāo) 記好的自封袋內(nèi),常溫保存帶回實(shí)驗(yàn)室,風(fēng)干后磨碎 過篩,用于土壤理化性質(zhì)測定。

    1.3土壤理化性質(zhì)及酶活性測定

    土壤pH采用酸度計(jì)法測定,有機(jī)質(zhì)含量采用重 鉻酸鉀容量法測定,堿解氮含量采用堿解擴(kuò)散法測 定,速效磷含量使用鉬酸比色法測定,速效鉀含量采 用NH2OAc浸提一火焰光度法測定。土壤脲酶活性 采用苯酚一次氯酸鈉比色法測定。以上指標(biāo)測定均 參照覃鋒燕等(2022)的方法。

    1.4土壤微生物高通量測序及生物信息學(xué)分析

    參照覃鋒燕等(2022)的方法,略有改動(dòng):使用E.Z.N.A.?Soil DNA Kit(美國Omega Bio-Tek公司) 提取土壤基因組DNA,用特異引物對(duì)16SrRNA序列的V3~V4可變區(qū)進(jìn)行PCR擴(kuò)增。PCR產(chǎn)物經(jīng)2% 瓊脂糖凝膠回收,用AxyPrep DNA Gel Extraction Kit(美國Axygen公司)進(jìn)行純化,2%瓊脂糖凝膠 電泳檢測,并用QuantusTMFluorometer(美國Promega 公司)進(jìn)行定量檢測。使用NEXTflex?Rapid DNA-Seq Kit(美國Bioo Scientific公司)構(gòu)建文庫,委托上海美吉生物醫(yī)藥科技有限公司利用Illumina MiSeqPE300平臺(tái)進(jìn)行土壤樣品高通量測序。原始數(shù)據(jù)上傳至NCBISRA數(shù)據(jù)庫。分別在 Trimmomatic(version 0.32)和FLASH(version 1.2.11)上進(jìn)行原始序列質(zhì)控和拼接,在97%相似度下使用UPARSE Version7.1將序列進(jìn)行操作分類單元(OTU)聚類?;赟ilva數(shù)據(jù)庫,設(shè)置70%的比對(duì)閾值,在RDP Classifier(http://rdp.cme.msu.edu/)上對(duì)篩選后的非嵌合序列進(jìn)行物種分類注釋。細(xì)菌OTU豐富度及Alpha多樣性分析分別使用Past3(http://fold.uio.no/ohammer/past)和Mothur(ver-sion v.1.30.2)進(jìn)行計(jì)算。細(xì)菌Beta多樣性分析包含主成分分析、置換多因素方差分析和冗余分析,均使用R(version 3.3.1)vegan包完成分析并制作相應(yīng)圖片。用stats包進(jìn)行細(xì)菌在屬水平上的組間差異顯著性檢驗(yàn)并制圖。

    1.5統(tǒng)計(jì)分析

    使用Excel 2010和SPSS 25.0進(jìn)行數(shù)據(jù)整理和統(tǒng)計(jì)分析;使用SIMCA14.0進(jìn)行主成分分析及制圖。

    2結(jié)果與分析

    2.1連作木薯對(duì)非根際和根際土壤理化性質(zhì)的影響

    由表1可知,隨連作年限增加,木薯非根際土壤的pH逐年下降,根際土壤的pH先下降后上升,年際間差異顯著(Plt;0.05,下同);非根際土壤和根際土壤的有機(jī)質(zhì)、堿解氮和速效磷含量均隨連作年限的增加逐年下降,而速效鉀含量在年際間均無顯著差異(Pgt;0.05,下同)。隨連作年限增加,非根際土壤的脲酶活性先升高后降低,根際土壤的脲酶活性逐年升高,連作2年和3年處理間無顯著差異,但均顯著高于連作1年處理。

    由圖2可看出,第一主成分(PC1)解釋了總方差 的59.10%,按年份將根際土壤與非根際土壤大致分 成3個(gè)模塊,說明連作年限是引起木薯根際土壤和 非根際土壤養(yǎng)分含量及脲酶活性差異的主要原因, 連作會(huì)引起木薯土壤理化性質(zhì)劣化。

    2.2連作木薯對(duì)非根際和根際土壤細(xì)菌群落多樣性及結(jié)構(gòu)的影響

    2.2.1Alpha多樣性分析表2為不同處理細(xì)菌群

    落多樣性指數(shù)(Shannon指數(shù)和Simpson指數(shù))和群落豐富度指數(shù)(Sobs指數(shù)、ACE指數(shù)和Chaol指數(shù))。各處理的OTU文庫覆蓋率均超過95.00%。隨連作年限增加,木薯非根際土壤的細(xì)菌Shannon指數(shù)呈 下降趨勢,連作1年與連作3年間差異顯著,而根際土壤的細(xì)菌Shannon指數(shù)無顯著變化;非根際和根 際土壤的細(xì)菌Simpson指數(shù)在年際間均無顯著差 異。隨連作年限增加,非根際和根際土壤的細(xì)菌 Sobs指數(shù)均呈下降趨勢,ACE指數(shù)先上升后下降; 非根際土壤細(xì)菌Chao1指數(shù)隨連作年限增加而逐年 下降,根際土壤細(xì)菌Chao1指數(shù)則先上升后下降。由此可知,連作木薯根際與非根際土壤細(xì)菌群落的 Alpha多樣性逐年降低,但根際土壤細(xì)菌群落的Alpha多樣性降低速度較非根際土壤慢。

    2.2.2Beta多樣性分析

    主成分分析結(jié)果(圖3)表明,PC1將非根際土壤與根際土壤連作1年和2年的 樣本與連作3年的樣本分為2個(gè)群體,分別解釋總方 差的17.18%和17.87%,故推測連作可引起細(xì)菌群落 Beta多樣性變化,連作3年差異達(dá)顯著水平。由表3可知,基于細(xì)菌群落半度量距離(Bray-Curtis)和未加權(quán)Unifrac距離,使用置換多因素方差分析方法分 別計(jì)算連作年限與土壤類型對(duì)細(xì)菌群落變化的貢 獻(xiàn),再次驗(yàn)證了連作年限是影響土壤細(xì)菌群落分類 和系統(tǒng)發(fā)育Beta多樣性差異的主要因素,而土壤類型對(duì)Beta多樣性差異的影響不顯著。

    2.2.3不同處理細(xì)菌群落組成分析

    從圖4可看 出,不同連作年限木薯的非根際和根際土壤細(xì)菌門 水平的群落結(jié)構(gòu)組成基本相同,主要有放線菌門(Actinobacteriota)、綠彎菌門(Chloroflexi)、變形菌 門(Proteobacteria)、酸桿菌門(Acidobacteriota)、厚壁 菌門(Firmicutes)、芽單胞菌門(Gemmatimonadota)、 Myxococcota、Methylomirabilota、擬桿菌門(Bacte- roidota)、浮霉菌門(Planctomycetota)、Latescibacte- rota、藍(lán)藻門(Cyanobacteria)和疣微菌門(Verrucomi- crobiota)等。其中,占主導(dǎo)地位的菌門為放線菌門 (占比19.18%~32.74%)、綠彎菌門(占比13.09%~21.56%)、變形菌門(占比14.08%~20.73%)、酸桿菌門(占比10.66%~22.83%)和厚壁菌門(占比1.16%~8.72%),占總相對(duì)豐度的77.14%~85.34%。從圖5可看出,在屬分類水平中,芽孢桿菌屬(Bacillus)、norank f_ Roseiflexaceae、 norank f_norank_o_Vicina-mibacterales 、norank_f_norank_o_Gaiellales為土壤細(xì)菌的主要優(yōu)勢菌群。其中,根際土壤中芽孢桿菌屬的相對(duì)豐度隨著連作年限的增加呈先增加后下降變化趨勢,非根際土壤中則是呈逐年下降趨勢。norank f_Roseiflexaceae的相對(duì)豐度在根際土壤中呈先升高后降低趨勢,在非根際土壤中呈連續(xù)下降趨勢。非根際土壤中norank f_norank_o_Vicina-mibacterales的相對(duì)豐度隨著連作年限增加呈下降趨勢,根際土壤中則先升高后降低。norankfnorank o Gaiellales的相對(duì)豐度在非根際土壤與根際土壤中均先升高后降低。綜上所述,不同連作年限木薯非根際和根際土壤細(xì)菌在門分類水平和屬分類水平的群落組成基本相同,但其細(xì)菌群落相對(duì)豐度隨連作年限增加而呈不同變化趨勢。

    2.2.4細(xì)菌群落組成差異分析

    從所有樣品的細(xì)菌群落中選取相對(duì)豐度均存在顯著或極顯著(Plt;0.01,下同)差異且排名前10位的 OTUs進(jìn)行細(xì)菌群落組成差異分析,結(jié)果(圖 6)表明,連作木薯非根際土壤細(xì)菌群落的前 10 位 OTUs 被注釋為RB41屬(OTU1199、OTU886)、Bacillus_sp._X12014種(OTU2085)、阿氏芽孢桿菌(OTU2064)、C0119目(OTU2693)、地桿菌(OTU2281)、Methyloligellaceae科(OTU1246)、Roseiflexaceae科(OTU646)、羅庫菌目(OTU2219)和Gaiellales目(OTU1045),分別屬于酸桿菌門、擬桿菌門、綠彎菌門、放線菌門和Methy-lomirabilota(圖 6-A)。其中 OTU1199、OTU2693、OTU646、OTU886 和OTU2219的相對(duì)豐度均在連作第年升高,且與前2年相比存在顯著或極顯著差異,OTU2085、OTU2064、OTU2281和OTU1246 的相對(duì)豐度均在連作第3年降低,且與前2年相比存在顯著或極顯著差異。

    從圖6-B可看出,連作木薯根際土壤細(xì)菌群落的前10位OTUs被注釋為SC-I-84科(OTU922)、Vic-inamibacteraceae科(OTU2006)、Bacillus_sp._X12014種(OTU2085)、阿氏芽孢桿菌(OTU2064)、芽單胞菌科(OTU1258)、C0119目(OTU2693)、地桿菌(OTU2281)、TK10綱(OTU2632)、Roseiflexaceae科(OTU1206)和RB41屬(OTU646),分別屬于變形菌門、酸桿菌門、擬桿菌門、芽單胞菌門和綠彎菌門。其中,OTU922、OTU2006、OTU2693、OTU1206和OTU886的相對(duì)豐度均在連作第3年升高,且與前2年相比存在顯著差異,OTU2085、OTU2064、OTU2281和OTU2632的相對(duì)豐度均在連作第3年降低,且與前2年相比存在顯著或極顯著差異??梢姡B作木薯非根際和根際土壤細(xì)菌群落中相對(duì)豐度均存在顯著或極顯著差異且排名前10位的0TUs大多在連作第3年出現(xiàn)顯著或極顯著變化。

    2.3連作木薯對(duì)土壤養(yǎng)分與細(xì)菌群落相關(guān)關(guān)系的影響

    在非根際土中(圖7-A),芽孢桿菌屬(Bacillus)與堿解氮和速效磷含量,norank_f_norank_o_Gaiel-lales、norank f JG30-KF-CM45、Gaiella與脲酶活性,分別呈顯著或極顯著正相關(guān)。norank_f_Gemmati-monadaceae與速效磷含量,norank f_norank_onorank c KD4-96與有機(jī)質(zhì)含量,norank f noranko_C0119和norank f SC-I-84與pH、堿解氮和速效磷含量,呈顯著或極顯著負(fù)相關(guān);在根際土中(圖7-B),norank f_norank o C0119與脲酶活性,norank fnorank_o_norank_c_KD4-96、Gaiella 和norank_fnorank o norank c TK10與速效鉀含量,norank fnorank o Vicinamibacterales與堿解氮含量,呈顯著或極顯著正相關(guān);RB41與有機(jī)質(zhì)含量,norankf_Gemmatimonadaceae與有機(jī)質(zhì)和速效磷含量,norank_fnorank_o_norank_c_KD4-96與pH,norank f_norank_o_C0119、norank_f_SC-I-84、Haliangium 和 norank_f_norank_o_Rokubacteriales 與有機(jī)質(zhì)、堿解 氮和速效磷含量,呈顯著或極顯著負(fù)相關(guān)??梢姡?壤細(xì)菌群落組成同時(shí)受土壤理化性質(zhì)和酶活性的影 響,在根際與非根際土環(huán)境中的影響程度存在差異。

    從圖8可看出,細(xì)菌群落豐度值最高的前5個(gè)細(xì) 菌門分別為變形菌門、酸桿菌門、綠彎菌門、厚壁菌 門和放線菌門,第一軸(RDA1)和第二軸(RDA2)的 解釋變量分別為40.31%和9.99%,2個(gè)變量在坐標(biāo) 對(duì)角線上把根際與非根際土壤樣本按年際分成3個(gè) 群體,其中第3年樣本與第1和第2年樣本距離較遠(yuǎn),說明木薯連作第3年開始對(duì)土壤理化性質(zhì)及微生物群落產(chǎn)生明顯影響。土壤pH、堿解氮含量、有機(jī)質(zhì)含量、速效磷含量和脲酶活性是對(duì)細(xì)菌群落結(jié)構(gòu)影響較大的環(huán)境因子,其中,脲酶活性與其他土壤理化性質(zhì)指標(biāo)均呈負(fù)相關(guān)。

    進(jìn)一步選擇門分類水平相對(duì)豐度大于0.10%的土壤細(xì)菌群落與理化性質(zhì)進(jìn)行相關(guān)分析,結(jié)果(表4)表明,有9個(gè)細(xì)菌門與土壤理化性質(zhì)有顯著或極顯著相關(guān)性。其中綠彎菌門與土壤pH呈極顯著負(fù)相關(guān),與有機(jī)質(zhì)、堿解氮和速效磷含量呈顯著負(fù)相關(guān);厚壁菌門與有機(jī)質(zhì)含量呈顯著正相關(guān),與堿解氮含量呈極顯著正相關(guān);芽單胞菌門與有機(jī)質(zhì)和堿解氮含量呈顯著負(fù)相關(guān);Methylomirabilota與堿解氮含量呈顯著負(fù)相關(guān);浮霉菌門與 pH.有機(jī)質(zhì)和堿解氮含量呈極顯著正相關(guān),與速效磷含量呈顯著正相關(guān),與脲酶活性呈極顯著負(fù)相關(guān);Unclassified_k_norankd_Bacteria與堿解氮含量呈顯著負(fù)相關(guān);疣微菌門和藍(lán)藻門與 pH、有機(jī)質(zhì)和堿解氮含量呈極顯著正相關(guān),與速效磷含量呈顯著正相關(guān),與脲酶活性呈極顯著或顯著負(fù)相關(guān)。

    綜上所述,土壤的細(xì)菌群落分布與土壤環(huán)境因 子相互影響,對(duì)細(xì)菌群落影響較大的環(huán)境因子排序 為堿解氮含量gt;有機(jī)質(zhì)含量gt;pHgt;速效磷含量gt;脲酶 活性。連作木薯會(huì)導(dǎo)致土壤pH及堿解氮、有機(jī)質(zhì)和 速效磷含量降低,進(jìn)而影響微生物群落演替,且連作 障礙在第3年開始顯現(xiàn)。

    3討論

    土壤理化性質(zhì)是評(píng)價(jià)土壤優(yōu)劣的重要指標(biāo)。本 研究中,與頭茬木薯相比,連作木薯的非根際土壤和根際土壤pH及有機(jī)質(zhì)、堿解氮和速效磷含量均降低,其中,第2年非根際和根際土壤有機(jī)質(zhì)含量急劇下降。前人研究表明,當(dāng)自然生態(tài)系統(tǒng)受到人為干擾轉(zhuǎn)換為農(nóng)業(yè)生態(tài)系統(tǒng)后,由于大量地面裸露,導(dǎo)致土壤有機(jī)碳庫在溫帶和熱帶地區(qū)分別降低60%和75%以上(陳虎等,2012;李娟等,2022)。本研究中,試驗(yàn)地位于亞熱帶地區(qū),在種植木薯前為撂荒地,近似于自然生態(tài)系統(tǒng)。第1年種植木薯時(shí),因前期除草效果較差,生育前期(4-6月)大量雜草覆蓋地表,減少了地面裸露引起的土壤有機(jī)質(zhì)消耗;塊根膨大及成熟期(9—12月)遭遇長期干旱,導(dǎo)致木薯生長緩慢,塊根平均產(chǎn)量低,木薯生長對(duì)土壤養(yǎng)分的消耗較小,且前期雜草、間苗秸稈及后期落葉還田對(duì)土壤有機(jī)質(zhì)碳庫有一定補(bǔ)充;收獲期塊根及莖稈生物量低,碳移除量較少,故第1年種植木薯后,土壤有機(jī)質(zhì)還保留在10~20g/kg的水平(張童瑤等,2020;張葉葉等,2021)。而種植木薯的第2年,地塊雜草相對(duì)較少,加上前期防控措施得當(dāng),自3月整地后,雜草更少,一直呈裸露狀態(tài),至7月后才開始封行,長期的土地裸露導(dǎo)致有機(jī)質(zhì)下降(陳虎等,2012);塊根膨大和成熟期雨量分布相對(duì)均勻,木薯整個(gè)生育期長勢較好,塊根平均產(chǎn)量提高,木薯生長對(duì)土壤養(yǎng)分消耗巨大,木薯和土壤微生物對(duì)土壤養(yǎng)分的競爭劇烈,進(jìn)而促進(jìn)微生物礦化分解有機(jī)質(zhì)釋放養(yǎng)分(渠晨晨等,2022);此外,第2年雜草較少,還田葉片也遠(yuǎn)少于移除的莖稈和塊根,因此土壤有機(jī)質(zhì)碳庫一直處于虧損狀態(tài)。另外,張葉葉等(2021)的研究表明,秸稈還田所引起的碳激發(fā)效應(yīng)也可能造成土壤有機(jī)質(zhì)下降;碳激發(fā)效應(yīng)與土壤有機(jī)碳含量存在一定關(guān)系,當(dāng)土壤有機(jī)碳含量在10~20g/kg時(shí)激發(fā)強(qiáng)度最高,高于或低于該閾值時(shí)激發(fā)變?nèi)?;長期單一作物秸稈還田不利于構(gòu)建多樣的微生物環(huán)境,易造成同一營養(yǎng)型微生物富集,對(duì)土壤有機(jī)碳造成長期損耗。本研究中,第1年收獲木薯后的土壤有機(jī)質(zhì)含量剛好處于最高激發(fā)強(qiáng)度區(qū)間,且第2年木薯長勢較好,相對(duì)強(qiáng)大的根系所產(chǎn)生的根系分泌物也較多,加上連作還田的落葉,可能導(dǎo)致同一營養(yǎng)型微生物富集,進(jìn)而對(duì)土壤有機(jī)碳造成損耗。此外,連作第2年土壤速效鉀含量較前、后2年偏高,可能是因?yàn)榈?年6月降雨量較大(為第1年和第3年6月降雨量的2倍),此時(shí)為木薯塊根形成初期,根系活動(dòng)旺盛,受 澇害脅迫會(huì)分泌大量有機(jī)酸和有機(jī)碳,在有機(jī)酸的 酸化作用下土壤速效鉀含量緩慢提高,但因?yàn)樽饔?時(shí)間短且作用面積小,速效鉀含量變化未達(dá)顯著水 平,該結(jié)果與Henry等(2007)、何冰等(2015)的研究 結(jié)果一致。綜上所述,土壤理化性質(zhì)在連作木薯第 3年明顯劣化,建議木薯連作不超過3年。該結(jié)果 也與梁海波等(2017)基于多年生產(chǎn)實(shí)際的預(yù)測結(jié)果 一致。

    在健康的土壤系統(tǒng)中,植物會(huì)分泌可供微生物 利用的營養(yǎng)物質(zhì),同時(shí)有益微生物也會(huì)以各種方式 協(xié)助植物吸收土壤養(yǎng)分,如改變土壤養(yǎng)分有效性、形 成從枝菌根等,且有益微生物能同時(shí)協(xié)助植物抵御 病蟲害和適應(yīng)逆境(孫子欣和蔡柏巖,2022)。而連 作引起的土壤理化性質(zhì)劣化和土壤微生物群落從 “有益菌”向“有害菌”改變是連作障礙形成的重要原 因(戴傳超和張偉,2022)。在連作馬鈴薯(謝奎忠, 2021)、三七(劉雨艷等,2022)、切花菊(王坦,2022) 等作物中均發(fā)現(xiàn)根際分泌物會(huì)改變微生物群落多樣 性。本研究結(jié)果表明,連作木薯根際與非根際土壤 細(xì)菌群落的Alpha多樣性會(huì)逐年降低,但根際土壤細(xì)菌群落的Alpha多樣性降低速度比非根際土壤慢。這可能是因?yàn)檫B作導(dǎo)致土壤養(yǎng)分下降,增強(qiáng)細(xì)菌群落之間的拮抗作用,進(jìn)而導(dǎo)致群落穩(wěn)定性下降;木薯根際分泌物中豐富的糖類和氨基酸能提供豐富的碳源和氮源,緩解根際土壤細(xì)菌群落對(duì)養(yǎng)分的競爭,故而保持較高的豐富度,從而導(dǎo)致非根際與根際土壤細(xì)菌群落結(jié)構(gòu)的多樣性和豐富度差異。

    本研究中,不同連作年限木薯非根際與根際土壤中的優(yōu)勢菌門均包含放線菌門、變形菌門、綠彎菌門、擬桿菌門和酸桿菌門,與劉詩蓉等(2022)對(duì)半夏連作根際土細(xì)菌群落結(jié)構(gòu)的研究結(jié)果相似?;诓煌B作年限間有顯著差異且最豐富的前 10個(gè)OTUs對(duì)比研究發(fā)現(xiàn),非根際與根際土壤中共有的優(yōu)勢菌群為 Bacillus_sp._X12014 種、阿氏芽孢桿菌、地桿菌、C0119 目、RB41屬和 Roseiflexaceae 科。Bacil-lus_sp._X12014種與阿氏芽孢桿菌均屬于芽孢桿菌屬,該菌屬中有許多能降解土壤有毒物質(zhì)的益生菌,能增強(qiáng)植物對(duì)外界環(huán)境的適應(yīng)和抵抗能力,其中,阿氏芽孢桿菌具有促進(jìn)植物吸收鋅離子、防治植物病害、促進(jìn)植物生長等作用(熊乙等,2018)。本研究發(fā)現(xiàn),非根際與根際土壤的 Bacillus_sp._X12014種與阿氏芽孢桿菌的相對(duì)豐度均在連作第3 年大幅下降,導(dǎo)致土壤降解有毒物質(zhì)的能力下降,這可能是引起木薯連作障礙的重要原因之一。地桿菌的相對(duì)豐度也在連作第3年顯著下降,其在厭氧條件下的亞硝酸鹽和硝酸鹽轉(zhuǎn)化為 NH的過程中扮演重要角色(田麗麗和王鑫,2020),這可能是引起土壤堿解氮含量隨連作年限增加而下降的原因之一。C0119目與土壤CO,的相互作用可促進(jìn)植物生長(Miao et al.,2021),而RB41被認(rèn)為是土壤中參與碳循環(huán)的關(guān)鍵菌屬(Stone et al.,2021),兩者的相對(duì)豐度均在連作木薯第3年升高,但由于其相對(duì)豐度較?。╨t;1.00%),有可能是由于相對(duì)豐度較高的芽孢桿菌屬大幅下降引起的。此外,木薯根際土壤中特有的優(yōu)勢菌群有SC-I-84科、Vicinamibacteraceae科、芽單胞菌科和TK10綱。其中,SC-I-84被證實(shí)為作物根際土壤特有優(yōu)勢菌群或根系內(nèi)生菌(龐師嬋等,2020;肖健等,2022),其相對(duì)豐度隨木薯連作年限逐年上升,這可能與某些木薯根系分泌物的逐年積累相關(guān);而Vici-namibacteraceae科屬酸桿菌門,其相對(duì)豐度隨連作年限呈先降后升趨勢,可能是由于連作第2年6月降雨過多使耕地受漬,甲烷菌大量繁殖,進(jìn)而導(dǎo)致第2年Vicinamibacteraceae科相對(duì)豐度下降。綜上所述,木薯非根際與根際土壤中相對(duì)豐度較高的細(xì)菌OTUs均在連作第3年出現(xiàn)顯著的相對(duì)豐度變化,其中Bacillus_sp._X12014、阿氏芽孢桿菌和地桿菌等有益菌豐度大幅下降,這可能是木薯連作障礙在連續(xù)種植第3年后開始顯現(xiàn)的重要原因之一,也直接印證了梁海波等(2017)的模型預(yù)測結(jié)果。此外,本研究還發(fā)現(xiàn),在連作木薯非根際與根際土壤中,pH、有機(jī)質(zhì)含量、堿解氮含量、速效磷含量、速效鉀含量和脲酶活性對(duì)土壤細(xì)菌群落結(jié)構(gòu)有顯著影響,但不 同環(huán)境因子對(duì)土壤細(xì)菌群落的影響也不同,表明在 不同生境中,土壤細(xì)菌對(duì)環(huán)境因子的響應(yīng)存在差異。

    4結(jié)論

    連作會(huì)導(dǎo)致木薯非根際和根際土壤理化性質(zhì)、 細(xì)菌群落多樣性和豐富度發(fā)生變化,土壤理化性質(zhì) 和細(xì)菌群落結(jié)構(gòu)均在連作第3年出現(xiàn)明顯劣化。土 壤理化性質(zhì)變化與根際土壤細(xì)菌群落結(jié)構(gòu)的差異存 在一定相關(guān)性,有益菌豐度降低可能是木薯連作造 成土壤細(xì)菌群落劣化演替的重要原因之一。

    參考文獻(xiàn)(References):

    陳?;郏昴死?,姜明國,王一兵.2022.作物重茬連作障礙中 自毒物質(zhì)的研究進(jìn)展[J].中國農(nóng)業(yè)科技導(dǎo)報(bào),24(10): 125-132. [Chen F H, Shen N K, Jiang M G, Wang Y B. 2022. Research progress of autotoxic secretions in crops replant successive cropping obstacles [J]. Journal of Agri- cultural Science and Technology, 24 (10) : 125-132.] doi: 10.13304/j.nykjdb.2021.0889.

    陳虎,王嘉學(xué),胡燦燦,馬照倫.2012.滇東喀斯特地區(qū)石漠化 裸露紅土的表層有機(jī)質(zhì)含量分析[J].中國巖溶,31(4): 423-425. [Chen H, Wang J X, Hu C C, Ma Z L. 2012. Organic matter content of the bare lateritic surface-layer in east Yunnan karst rocky deserts [J]. Carsologica Sinica,31(4):423-425.] doi:10.3969/j.issn.1001-4810.2012.04.011.

    戴傳超,張偉.2022.有益微生物緩解花生連作障礙機(jī)理研究 進(jìn)展[J].微生物學(xué)雜志,42(1):1-8.[DaiCC,Zhang W.

    2022. Research progress on the beneficial microbe-mediated alleviation of peanut continuous cropping obstacle[J]. Journal of Microbiology,42(1):1-8.] doi:10.3969/j.issn.1005-7021.2022.01.001.

    何冰,薛剛,張小全,許曉敬,姚健,楊鐵釗.2015.有機(jī)酸對(duì)土 壤鉀素活化過程的化學(xué)分析[J].土壤,47(1):74-79.[He B, Xue G, Zhang X Q, Xu X J, Yao J, Yang T Z. 2015. Analysis on chemical mechanism of potassium release pro- cess from soil as influenced by organic acids[J]. Soils,47(1):74-79.] doi:10.13758/j.cnki.tr.2015.01.012.

    侯乾,王萬興,李廣存,熊興耀.2019.馬鈴薯連作障礙研究進(jìn) 展[J].作物雜志,(6):1-7.[Hou Q,Wang WX,Li G C, Xiong X Y. 2019. Advances in the research on potato con- tinuous cropping obstacles [J]. Crops, (6) : 1-7.] doi: 10.16035/j.issn.1001-7283.2019.06.001.

    李娟,魏甲彬,楊寧.2022.土地利用變化對(duì)土壤有機(jī)質(zhì)影響 的研究進(jìn)展及展望[J].湖南生態(tài)科學(xué)學(xué)報(bào),9(3):106-111.李軍,魏建波,楊。土地利用變化對(duì)土壤有機(jī)質(zhì)影響的研究進(jìn)展與展望[J].湖南生態(tài)學(xué)報(bào),9 (3) : 106-111?!筪oi:10.3969/j . ISSN . 2095-7300 . 2022 . 03 . 014。

    梁海波,黃潔,魏云霞,陸昆典,劉翠娟.2017.基于農(nóng)戶尺度 的木薯產(chǎn)量差及生產(chǎn)限制因素分析[J].云南農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)),32(6):975-984.[Liang H B,HuangJ,Wei Y X, Lu K D,Liu C J. 2017. Analysis of yield gap andlimiting factors for cassava on the farmland[J]. Journal ofYunnan Agricultural University,32(6) : 975-984.] doi: 10.16211/j.issn.1004-390X(n).2017.06.003.

    劉珊廷.2020.木薯連作與輪作對(duì)土壤理化性狀及微生物群落和產(chǎn)量的影響[D].南寧:廣西大學(xué).[LiuST.2020.Causes and mechanism on crop continuous monocultureproblems and its control strategy [D]. Nanning: GuangxiUniversity.] doi:10.27034/d.cnki.ggxiu.2020.000417.

    劉詩蓉,王紅蘭,孫輝,楊萍,蔣舜媛,蔣桂華.2022.半夏連作對(duì)根際土壤微生物群落的影響研究[J].中草藥,53(4):1148-1155. [Liu S R, Wang H L, Sun H, Yang P, Jiang SY, Jiang G H. 2022. Effects of continuous cropping ofPinellia ternata on rhizospheric microbial community [J].Chinese Traditional and Herbal Drugs, 53 (4) : 1148-1155.]doi: 10.7501/j.issn.0253-2670.2022.04.023.

    劉雨艷,保麗美,魏云林,字富庭,譚勇.2022.連作三七根際土壤自毒物質(zhì)與微生物互作研究進(jìn)展[J].中藥材,45(1):242-247.[Liu Y Y,Bao L M,Wei Y L,Zi F T,Tan Y.2022. Research progress on interaction between autotoxicsubstances and microorganisms in rhizosphere soil ofPanax notoginseng in continuous cropping[J]. Journal ofChinese Medicinal Materials , 45 (1) : 242-247. ] doi : 10.13863/j.issn1001-4454.2022.01.043.

    劉株秀,劉俊杰,徐艷霞,張武,米剛,姚欽,王光華.2019.不 同大豆連作年限對(duì)黑土細(xì)菌群落結(jié)構(gòu)的影響[J].生態(tài)學(xué)報(bào),39(12):4337-4346.[Liu ZX,Liu JJ,XuYX,Zhang W, Mi G, Yao Q, Wang G H. 2019. Effects of continuous cropping years of soybean on the bacterial community structure inblack soil [J]. Acta Ecologica Sinica, 39 (12) :

    4337-4346.] doi:10.5846/stxb201801270212.

    盧孟召,劉梅,陳光,王雪峰.2022.靈芝連作對(duì)土壤理化性質(zhì) 及線蟲群落的影響[J].吉林農(nóng)業(yè)大學(xué)學(xué)報(bào),44(5):586- 594. [Lu M Z, Liu M, Chen G, Wang X F. 2022. Effect of biochar-based fertilizer on soil bacteria and fungi quantity and community structure in acidified tea garden [J]. Jour- nal of Jilin Agricultural University, 44 (5) : 586-594.] doi: 10.13327/j.jjlau.2022.1487.

    龐師嬋,郭霜,任奎瑜,王帥帥,楊尚東.2020.番茄/茄子嫁接 對(duì)其根際土壤生物學(xué)性狀及細(xì)菌群落結(jié)構(gòu)的影響[J].園藝學(xué)報(bào),47(2):253-263.[Pang S C,Guo S,Ren KY, Wang S S, Yang S D. 2020. Impact of grafting on soil microbial properties and bacterial community structure in tomato rhizosphere [J]. Acta Horticulturae Sinica, 47 (2) :253-263.] doi: 10.16420/j.issn.0513-353x.2019-0444.

    覃鋒燕,楊慰賢,彭曉輝,黎亮武,陽太億,周佳,韋茂貴,申章佑.2022.粉壟耕作木薯根際與非根際土壤的細(xì)菌群落結(jié)構(gòu)多樣性差異[J].西南農(nóng)業(yè)學(xué)報(bào),35(4):729-739.[QinF Y,Yang W X,Peng X H,Li L W,Yang T Y,Zhou J,WeiM G, Shen Z Y. 2022. Difference in the diversity of bacte-rial community structure in rhizosphere and non-rhizospheresoil of cassava in Fenlong tillage[J]. Southwest China Jour-nal of Agricultural Sciences, 35 (4) : 729-739.] doi: 10.16213/j.cnki.scjas.2022.4.001.

    渠晨晨,任穩(wěn)燕,李秀秀,蔡鵬,陳雯莉,黃巧云.2022.重新認(rèn)識(shí)土壤有機(jī)質(zhì)[J].科學(xué)通報(bào),67(10):913-923.[QuCC, Ren W Y, Li X X, Cai P, Chen W L, Huang Q Y. 2022. Revisit soil organic matter[J]. Chinese Science Bulletin,67(10):913-923.] doi:10.1360/TB-2021-0704.

    宋秀麗,吳舒婷,李錦輝,楊錫朗,周芮,況丹妮,李警督,于海 濤,楊德光.2022.土壤微生物群落對(duì)連作種植制度的響 應(yīng)[J].玉米科學(xué),30(1):172-181.[Song XL,Wu S T,Li J H,Yang XL,Zhou R,Kuang D N,Li J D,Yu H T,Yang D G. 2022. Response of soil microbial community to con- tinuous cropping system[J]. Journal of Maize Sciences,30(1) : 172-181. ] doi : 10.13597/j.cnki.maize.science.20220124.

    孫子欣,蔡柏巖.2022.連作對(duì)土壤微生物菌群影響及修復(fù)研 究進(jìn)展[J].作物雜志,(6):7-13.[SunZX,Cai B Y.2022.Research advances on the effects of continuous cropping on soil microbial community and restoration techniques [J]. Crops, (6):7-13.] doi:10.16035/j.issn.1001-7283.2022.06.002.

    田麗麗,王鑫.2020.地桿菌屬微生物在厭氧污水處理中的應(yīng) 用潛力[J].給水排水,56(S2):201-206.[Tian LL,Wang X. 2020. Potential applications of Geobacterin anaerobic wastewater treatment[J]. Water amp; Wastewater Engineering, 56 (S2) : 201-206. ] doi : 10.13789/j.cnki. wwe1964.2020.S2.030

    王坦.2022.連作下切花菊根際微生物群落組成和功能的動(dòng) 態(tài)變化[D].武漢:華中農(nóng)業(yè)大學(xué).[Wang T.2022.Dynamic changes of rhizosphere microbial community composition and function in cut chrysanthemum under continuous crop- ping[D]. Wuhan: Huazhong Agricultural University.] doi: 10.27158/d.cnki.ghznu.2022.000706.

    王曉琪,姜偉濤,姚媛媛,尹承苗,陳學(xué)森,毛志泉.2020.蘋果 連作障礙土壤微生物的研究進(jìn)展[J].園藝學(xué)報(bào),47(11): 2223-2237. [Wang X Q, Jiang W T, Yao Y Y, Yin C M, Chen X S, Mao Z Q. 2020. Research advance of apple replant disease based on soil microorganism[J]. Acta Hor- ticulturae Sinica,47(11):2223-2237.] doi: 10.16420/j.issn.0513-353x.2019-0941.

    王戰(zhàn).2016.木薯磷營養(yǎng)與連作障礙的關(guān)系研究[D].南寧:廣西大學(xué).[Wang Z.2016.Study on the relationship be- tween phosphorus nutrition and continuous cropping obstacle in cassava [D]. Nanning:Guangxi University.]

    翁佩瑩,鄭紅艷.2020.作物連作障礙的成因與機(jī)制及其消減策略[J].亞熱帶植物科學(xué),49(2):157-162.[Weng PY,Zheng H Y. 2020. Causes and mechanism on crop conti-nuous monoculture problems and its control strategy [J].Subtropical Plant Science, 49 (2) : 157-162.] doi: 10.3969/j.issn.2095-1191.2021.06.005.

    肖健,梁闐,楊尚東,譚宏偉.2022.甘蔗植株莖部內(nèi)生細(xì)菌群落組成對(duì)不同配方緩釋肥的響應(yīng)[J].熱帶作物學(xué)報(bào),44(3):614-627. [Xiao J, Liang T, Yang S D, Tan H W. 2022.

    Response of endophytic bacterial community compositions in stems of sugarcanes under different slow-release fertilizers applications[J]. Chinese Journal of Tropical Crops,44

    (3):614-627.] doi:10.3969/j.issn.1000-2561.2023.03.019.

    謝奎忠.2021.連作馬鈴薯根系分泌物介導(dǎo)的枯萎病發(fā)生機(jī) 制及根際互作[D].蘭州:甘肅農(nóng)業(yè)大學(xué).[XieKZ.2021.

    Mechanism of Fusarium wilt and rhizosphere interaction mediated by root exudates in continuous cropping potato [D]. Lanzhou: Gansu Agricultural University.] doi: 10.27025/d.cnki.ggsnu.2021.000041.

    熊乙,歐翔,賈蓉,玉柱,許慶方.2018.阿氏芽孢桿菌應(yīng)用研 究進(jìn)展[J].生物技術(shù),28(3):302-306.[Xiong Y,Ou X, Jia R, Yu Z, Xu Q F. 2018. Research progress in applica- tion of Bacillus aryabhattai [J]. Biotechnology, 28 (3) :302-306.] doi:10.16519/j.cnki.1004-311x.2018.03.0052.

    楊浩娜,周成言,鄔臘梅,王立峰,李祖任.2022.植物化感物 質(zhì)的作用機(jī)理研究進(jìn)展[J].湖南農(nóng)業(yè)科學(xué),(3):108-112.[Yang H N,Zhou C Y,Wu L M,Wang L F,Li Z R. 2022. Research progress of plant allelochemicals mechanism[J].Hunan Agricultural Sciences, (3): 108-112.] doi:10.16498/ j.cnki.hnnykx.2022.003.025.

    張童瑤,胡月明,任向?qū)?,陳飛香,馮雪珂.2020.棄耕行為對(duì) 亞熱帶農(nóng)田土壤有機(jī)質(zhì)時(shí)空變化作用機(jī)理研究[J].農(nóng)業(yè)資源與環(huán)境學(xué)報(bào),37(6):805-817.[ZhangTY,Hu Y M, Ren X N,Chen F X,F(xiàn)eng X K. 2020. Study on the spatio- temporal variation of soil organic matter induced by aban- doned tillage behavior[J]. Journal of Agricultural Resour- ces and Environment, 37 (6) : 805-817.] doi: 10.13254/j.jare.2020.0480.

    張葉葉,莫非,韓娟,溫曉霞,廖允成.2021.秸稈還田下土壤 有機(jī)質(zhì)激發(fā)效應(yīng)研究進(jìn)展[J].土壤學(xué)報(bào),58(6):1381- 1392. [Zhang Y Y, Mo F, Han J, Wen X X, Liao Y C.2021. Research progress on the native soil carbon priming after straw addition [J]. Acta Pedologica Sinica, 58 (6) :1381-1392.] doi:10.11766/trxb202006260259.

    周貴靖.2017.木薯根際土壤理化性狀與連作障礙關(guān)系研究[D].南寧:廣西大學(xué).[Zhou G J.2017.The relationship between the continous cropping obstacle and physical and chemical properties of cassava rhizosphere soil abstract[D]. Nanning:Guangxi University.]

    Henry A, Doucette W, Norton J, Bugbee B. 2007. Changes in crested wheatgrass root exudation caused by flood,drought, and nutrient stress[J]. Journal of Environmental Quality,36(3):904-912. doi: 10.2134/jeq2006.0425sc.

    Hontoria C, Garcia-Gonzalez I, Quemada M, Roldan A, Alguacil M M. 2019. The cover crop determines the AMF community composition in soil and in roots of maize after a ten-year continuous crop rotation [J]. Science of the Total Environment, 660: 913-922. doi: 10.1016/j.scitotenv.2019.01.095.

    Miao S J, Zhao H F, Pan S Q, Zhang N, Qiao Y F. 2021. Elevated CO2 weakens the shift in bacterial community structure in response to 8-year soybean straw return in the same experiment [J]. International Journal of Phytoreme-diation, 23 (5): 505-510. doi: 10.1080/15226514.2020.1825332.

    Stone B W, Li J H,Koch B J, Blazewicz S J,Dijkstra P, Hayer M, Hofmockel K S, Liu X A, Mau R L, Morrissey E M, Pett-Ridge J, Schwartz E, Hungate B A. 2021. Nutrients cause consolidation of soil carbon flux to small proportion of bacterial community[J]. Nature Communications, 12

    (1):3381. doi:10.1038/s41467-021-23676-x.

    (責(zé)任編輯王暉)

    猜你喜歡
    根際土壤木薯
    刮木薯
    四方聯(lián)動(dòng)將剛果(布)打造成木薯全產(chǎn)業(yè)鏈提升的案例國
    柬埔寨拜靈木薯喜獲大豐收,市場價(jià)格保持穩(wěn)定
    挖木薯
    不同種植模式對(duì)烤煙根際土壤養(yǎng)分、酶活性及煙草生長的影響
    馬纓丹入侵對(duì)根際和非根際土壤理化性質(zhì)的影響
    小蓬竹根際土壤可培養(yǎng)細(xì)菌多樣性
    華重樓根際土芽胞桿菌多樣性研究
    廣藿香根際土壤微生物總DNA提取方法的優(yōu)化
    尼日利亞木薯農(nóng)民合作聯(lián)盟簽署協(xié)議促進(jìn)木薯生產(chǎn)
    亚洲成a人片在线一区二区| 国产亚洲精品久久久com| 久久6这里有精品| 18禁美女被吸乳视频| 午夜亚洲福利在线播放| 给我免费播放毛片高清在线观看| 国产单亲对白刺激| 国产亚洲欧美98| 亚洲久久久久久中文字幕| 午夜久久久久精精品| xxx96com| 成人永久免费在线观看视频| 欧美区成人在线视频| 在线观看av片永久免费下载| 午夜a级毛片| 他把我摸到了高潮在线观看| 久久久精品欧美日韩精品| 欧美成人a在线观看| 禁无遮挡网站| 国产97色在线日韩免费| 国产又黄又爽又无遮挡在线| 99精品在免费线老司机午夜| 国产探花极品一区二区| 欧美日韩黄片免| 两个人的视频大全免费| 欧美av亚洲av综合av国产av| 岛国视频午夜一区免费看| 亚洲国产欧美人成| 日韩高清综合在线| 最新在线观看一区二区三区| 国产成人啪精品午夜网站| 国产精品一区二区免费欧美| 国产免费男女视频| 国内精品美女久久久久久| 亚洲专区中文字幕在线| 久久久成人免费电影| 长腿黑丝高跟| 午夜亚洲福利在线播放| 欧美成狂野欧美在线观看| 一进一出抽搐动态| av黄色大香蕉| 91在线观看av| 999久久久精品免费观看国产| 日韩欧美国产一区二区入口| 久久精品国产自在天天线| 五月玫瑰六月丁香| 搡老妇女老女人老熟妇| 久久午夜亚洲精品久久| 校园春色视频在线观看| 国产探花在线观看一区二区| 老司机在亚洲福利影院| 老熟妇仑乱视频hdxx| 91字幕亚洲| 男人舔奶头视频| 在线观看一区二区三区| 久久精品综合一区二区三区| 一夜夜www| 舔av片在线| 亚洲五月天丁香| 少妇人妻精品综合一区二区 | 1000部很黄的大片| 日韩成人在线观看一区二区三区| 一二三四社区在线视频社区8| av视频在线观看入口| 午夜精品在线福利| 亚洲国产精品久久男人天堂| 久久精品国产亚洲av涩爱 | 亚洲第一电影网av| 久久国产精品人妻蜜桃| 好看av亚洲va欧美ⅴa在| 又紧又爽又黄一区二区| 国产午夜精品久久久久久一区二区三区 | 精品熟女少妇八av免费久了| 久久久久精品国产欧美久久久| 啦啦啦观看免费观看视频高清| 午夜久久久久精精品| 午夜免费男女啪啪视频观看 | 亚洲第一电影网av| 免费人成在线观看视频色| 欧美日韩瑟瑟在线播放| 亚洲精品在线美女| 他把我摸到了高潮在线观看| 一个人观看的视频www高清免费观看| 亚洲精品粉嫩美女一区| 嫁个100分男人电影在线观看| 成人鲁丝片一二三区免费| 99热这里只有是精品50| 欧美精品啪啪一区二区三区| 嫩草影视91久久| 黄色女人牲交| 一个人观看的视频www高清免费观看| 成年人黄色毛片网站| 亚洲人成网站在线播| 国产精品电影一区二区三区| 亚洲专区国产一区二区| 亚洲精品日韩av片在线观看 | 波多野结衣高清无吗| 久久久久久久久久黄片| 制服丝袜大香蕉在线| 欧美色视频一区免费| 色综合欧美亚洲国产小说| 国产 一区 欧美 日韩| 久久精品国产亚洲av涩爱 | 亚洲av电影在线进入| 老鸭窝网址在线观看| 亚洲人与动物交配视频| 国内久久婷婷六月综合欲色啪| 日韩欧美在线乱码| av福利片在线观看| 3wmmmm亚洲av在线观看| 精品国产三级普通话版| 久久久久久国产a免费观看| 男女视频在线观看网站免费| 欧美色欧美亚洲另类二区| 久久久久久久亚洲中文字幕 | 看片在线看免费视频| 一级毛片女人18水好多| 国产伦一二天堂av在线观看| 观看免费一级毛片| 一本综合久久免费| 久久久成人免费电影| tocl精华| 国产伦精品一区二区三区四那| 国产三级中文精品| 免费观看人在逋| 欧美一级毛片孕妇| 精品国产亚洲在线| 日本精品一区二区三区蜜桃| 欧美性猛交黑人性爽| 日本 欧美在线| 亚洲第一欧美日韩一区二区三区| 亚洲无线观看免费| 国产精品影院久久| 亚洲国产精品999在线| 精华霜和精华液先用哪个| 床上黄色一级片| 亚洲五月天丁香| e午夜精品久久久久久久| 日韩大尺度精品在线看网址| 久久久久久九九精品二区国产| 国产精品香港三级国产av潘金莲| 久久久久久久精品吃奶| 国产视频内射| 午夜免费男女啪啪视频观看 | 亚洲第一欧美日韩一区二区三区| 99热这里只有是精品50| av专区在线播放| 亚洲国产精品合色在线| 最好的美女福利视频网| 精品乱码久久久久久99久播| 国模一区二区三区四区视频| 国产成人啪精品午夜网站| 黄色丝袜av网址大全| 少妇熟女aⅴ在线视频| 精品人妻偷拍中文字幕| 国产精品爽爽va在线观看网站| 一本一本综合久久| 国产97色在线日韩免费| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利高清视频| 亚洲av成人av| 一边摸一边抽搐一进一小说| 日本在线视频免费播放| 国产精品久久久久久久电影 | 一个人观看的视频www高清免费观看| 天天添夜夜摸| 色尼玛亚洲综合影院| 日本黄色片子视频| 色视频www国产| 一本精品99久久精品77| 很黄的视频免费| 麻豆国产av国片精品| 亚洲欧美激情综合另类| 欧美中文日本在线观看视频| 狂野欧美白嫩少妇大欣赏| 国产精品久久视频播放| 人妻久久中文字幕网| 波多野结衣巨乳人妻| 成人av一区二区三区在线看| 在线免费观看的www视频| 夜夜爽天天搞| 成人18禁在线播放| 变态另类成人亚洲欧美熟女| 九九久久精品国产亚洲av麻豆| 亚洲av一区综合| 熟女人妻精品中文字幕| 午夜福利成人在线免费观看| 欧美日韩精品网址| 香蕉av资源在线| 亚洲精品影视一区二区三区av| 亚洲成a人片在线一区二区| 在线播放无遮挡| 久久久久久九九精品二区国产| 亚洲,欧美精品.| 日本a在线网址| 精品无人区乱码1区二区| 精品一区二区三区av网在线观看| 欧美三级亚洲精品| 国产精品久久久人人做人人爽| 99国产综合亚洲精品| 18禁黄网站禁片午夜丰满| 亚洲在线自拍视频| 熟女电影av网| 国产91精品成人一区二区三区| www国产在线视频色| 精品久久久久久久久久免费视频| 日本成人三级电影网站| 欧美午夜高清在线| 精品免费久久久久久久清纯| 欧美丝袜亚洲另类 | 两个人看的免费小视频| 白带黄色成豆腐渣| 美女黄网站色视频| 午夜福利视频1000在线观看| 在线播放国产精品三级| 成人av一区二区三区在线看| 88av欧美| 亚洲最大成人手机在线| 国产成人av激情在线播放| 三级男女做爰猛烈吃奶摸视频| 亚洲精品成人久久久久久| 国产精品精品国产色婷婷| www.www免费av| 国产精品永久免费网站| 精品福利观看| 精品国产三级普通话版| 18禁裸乳无遮挡免费网站照片| 日本在线视频免费播放| 日本成人三级电影网站| 18禁国产床啪视频网站| 桃红色精品国产亚洲av| 国产高清三级在线| 国产视频一区二区在线看| 国产精品久久久人人做人人爽| 国产精品久久久久久人妻精品电影| 亚洲成人久久爱视频| 麻豆一二三区av精品| 国产高潮美女av| 国产精品免费一区二区三区在线| 亚洲人成电影免费在线| 成人一区二区视频在线观看| 99久久久亚洲精品蜜臀av| 老汉色av国产亚洲站长工具| 久99久视频精品免费| 日日夜夜操网爽| 国产乱人视频| 日本a在线网址| 国产精品久久久久久久久免 | 97碰自拍视频| 成年人黄色毛片网站| 国产毛片a区久久久久| 19禁男女啪啪无遮挡网站| 丰满人妻熟妇乱又伦精品不卡| 男女午夜视频在线观看| 给我免费播放毛片高清在线观看| 麻豆成人av在线观看| 99国产综合亚洲精品| 欧美一级a爱片免费观看看| 一区福利在线观看| 欧美黄色淫秽网站| 成人国产一区最新在线观看| tocl精华| 精品国产超薄肉色丝袜足j| 日本a在线网址| 看片在线看免费视频| 免费电影在线观看免费观看| 国产精品久久电影中文字幕| 3wmmmm亚洲av在线观看| 国内精品久久久久久久电影| 12—13女人毛片做爰片一| 免费人成视频x8x8入口观看| 成人高潮视频无遮挡免费网站| 欧美性猛交黑人性爽| 国产精品一区二区免费欧美| 一区福利在线观看| 午夜福利成人在线免费观看| 2021天堂中文幕一二区在线观| 久久九九热精品免费| 亚洲av五月六月丁香网| 精品免费久久久久久久清纯| 久久精品国产综合久久久| 男女床上黄色一级片免费看| 国产高清三级在线| 国产精品,欧美在线| 99久久综合精品五月天人人| 午夜影院日韩av| 99在线视频只有这里精品首页| 亚洲色图av天堂| 香蕉丝袜av| 麻豆一二三区av精品| 美女高潮喷水抽搐中文字幕| 国产av一区在线观看免费| 在线视频色国产色| 首页视频小说图片口味搜索| 午夜免费成人在线视频| 精品无人区乱码1区二区| 两个人的视频大全免费| 国产淫片久久久久久久久 | 亚洲精品亚洲一区二区| 国产麻豆成人av免费视频| 熟女电影av网| 黄片大片在线免费观看| 麻豆成人av在线观看| 成人午夜高清在线视频| 亚洲国产精品合色在线| 久久久久久久精品吃奶| 国产伦一二天堂av在线观看| 超碰av人人做人人爽久久 | 亚洲av免费在线观看| 色尼玛亚洲综合影院| 五月伊人婷婷丁香| www国产在线视频色| 日本五十路高清| 麻豆成人午夜福利视频| 一区二区三区高清视频在线| 国产精品嫩草影院av在线观看 | 亚洲人成网站高清观看| 高清日韩中文字幕在线| 国产 一区 欧美 日韩| 婷婷六月久久综合丁香| 亚洲人成网站在线播| 国产成人系列免费观看| 91麻豆av在线| 男人舔女人下体高潮全视频| 法律面前人人平等表现在哪些方面| 日本与韩国留学比较| 亚洲人成电影免费在线| 又粗又爽又猛毛片免费看| 一级作爱视频免费观看| 在线观看av片永久免费下载| 国产麻豆成人av免费视频| 三级国产精品欧美在线观看| 午夜精品在线福利| 成人av一区二区三区在线看| 久99久视频精品免费| 国产免费一级a男人的天堂| 久久性视频一级片| 成年女人看的毛片在线观看| 男插女下体视频免费在线播放| 久久国产精品影院| 嫁个100分男人电影在线观看| 波野结衣二区三区在线 | 91在线观看av| 我要搜黄色片| 亚洲最大成人中文| 两个人视频免费观看高清| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清有码在线观看视频| 亚洲色图av天堂| 丝袜美腿在线中文| 黄色成人免费大全| 久久久久九九精品影院| 日韩中文字幕欧美一区二区| 久久久久九九精品影院| 亚洲一区二区三区色噜噜| 俺也久久电影网| 精品福利观看| 十八禁人妻一区二区| 亚洲人成网站在线播放欧美日韩| 嫩草影院入口| 麻豆成人av在线观看| 欧美成人a在线观看| 日本三级黄在线观看| www日本黄色视频网| 女人被狂操c到高潮| 色在线成人网| 人妻久久中文字幕网| 国产一区二区在线观看日韩 | 黄色成人免费大全| 一本久久中文字幕| 欧美最新免费一区二区三区 | 在线观看免费午夜福利视频| 日韩欧美国产一区二区入口| 欧美中文日本在线观看视频| 欧美bdsm另类| 日韩欧美在线二视频| 搡女人真爽免费视频火全软件 | 3wmmmm亚洲av在线观看| 久久伊人香网站| 国产精品久久久久久精品电影| 亚洲不卡免费看| 日韩欧美一区二区三区在线观看| 国产在线精品亚洲第一网站| 免费高清视频大片| 中文字幕av在线有码专区| 国产亚洲欧美在线一区二区| 级片在线观看| 91av网一区二区| 99久国产av精品| 神马国产精品三级电影在线观看| 美女cb高潮喷水在线观看| 欧美日韩精品网址| 男人的好看免费观看在线视频| 国产午夜福利久久久久久| 给我免费播放毛片高清在线观看| 国产乱人视频| 女人被狂操c到高潮| 97人妻精品一区二区三区麻豆| 国产单亲对白刺激| 男女做爰动态图高潮gif福利片| av黄色大香蕉| 国产av不卡久久| 99精品在免费线老司机午夜| 亚洲无线在线观看| 一个人看视频在线观看www免费 | ponron亚洲| 在线免费观看不下载黄p国产 | 国产爱豆传媒在线观看| 在线免费观看不下载黄p国产 | 欧美+亚洲+日韩+国产| 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 日本黄色片子视频| 热99在线观看视频| 国产成人av激情在线播放| 国产欧美日韩精品亚洲av| 免费人成在线观看视频色| 成人特级av手机在线观看| 中文字幕av在线有码专区| 俄罗斯特黄特色一大片| 757午夜福利合集在线观看| 亚洲精品影视一区二区三区av| 国产亚洲精品久久久久久毛片| 一本久久中文字幕| 精品人妻偷拍中文字幕| 久久久国产成人精品二区| 精华霜和精华液先用哪个| 国产精品久久久久久久久免 | 国产精品久久久人人做人人爽| 九色国产91popny在线| 一级a爱片免费观看的视频| 麻豆国产av国片精品| 两个人看的免费小视频| 久久人妻av系列| 精品人妻1区二区| 性色avwww在线观看| 啦啦啦观看免费观看视频高清| 精华霜和精华液先用哪个| 精品99又大又爽又粗少妇毛片 | 99视频精品全部免费 在线| 久久性视频一级片| 此物有八面人人有两片| 蜜桃亚洲精品一区二区三区| 国产真人三级小视频在线观看| 亚洲不卡免费看| 久久婷婷人人爽人人干人人爱| 日本在线视频免费播放| 国产av不卡久久| 在线播放无遮挡| 叶爱在线成人免费视频播放| 精品福利观看| 久久久国产成人免费| av天堂在线播放| 久久草成人影院| 他把我摸到了高潮在线观看| 国产欧美日韩精品一区二区| 床上黄色一级片| 亚洲中文日韩欧美视频| 51午夜福利影视在线观看| 日本成人三级电影网站| 国产精品美女特级片免费视频播放器| 色在线成人网| av黄色大香蕉| 色综合站精品国产| 少妇高潮的动态图| 真人一进一出gif抽搐免费| 色哟哟哟哟哟哟| 国产精品美女特级片免费视频播放器| 丁香欧美五月| 中文字幕人妻熟人妻熟丝袜美 | 国产精品99久久久久久久久| 国产精品 国内视频| 老汉色∧v一级毛片| 美女 人体艺术 gogo| 亚洲国产精品999在线| 精品人妻一区二区三区麻豆 | 久久国产精品人妻蜜桃| 琪琪午夜伦伦电影理论片6080| 内射极品少妇av片p| 一边摸一边抽搐一进一小说| 久久精品人妻少妇| 蜜桃亚洲精品一区二区三区| 久久久久久久久久黄片| 韩国av一区二区三区四区| 69av精品久久久久久| 全区人妻精品视频| 51国产日韩欧美| 日韩 欧美 亚洲 中文字幕| 两个人看的免费小视频| 极品教师在线免费播放| 国产老妇女一区| 黄色片一级片一级黄色片| 亚洲国产精品sss在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美一区二区国产精品久久精品| 日韩精品中文字幕看吧| 日本一二三区视频观看| 免费看a级黄色片| 国产精品永久免费网站| 欧美日韩中文字幕国产精品一区二区三区| 国产精品国产高清国产av| 日韩亚洲欧美综合| 成人三级黄色视频| 国产成人a区在线观看| 极品教师在线免费播放| 亚洲人成网站在线播放欧美日韩| 老司机午夜十八禁免费视频| 精品国产三级普通话版| 国产激情偷乱视频一区二区| 日日干狠狠操夜夜爽| 国产99白浆流出| 欧美黄色淫秽网站| 99热只有精品国产| 欧美一级a爱片免费观看看| 免费大片18禁| 成年女人看的毛片在线观看| 成人鲁丝片一二三区免费| 美女黄网站色视频| 一本综合久久免费| 毛片女人毛片| 亚洲国产精品999在线| 国产精品一区二区免费欧美| 久久久国产精品麻豆| 少妇高潮的动态图| 亚洲成人久久爱视频| av欧美777| 免费在线观看日本一区| 黄色视频,在线免费观看| 在线观看免费视频日本深夜| xxxwww97欧美| 香蕉av资源在线| netflix在线观看网站| 在线观看午夜福利视频| 国产男靠女视频免费网站| 国产综合懂色| 极品教师在线免费播放| 国产精华一区二区三区| av天堂在线播放| 亚洲天堂国产精品一区在线| 国产精品女同一区二区软件 | 老司机在亚洲福利影院| 国产麻豆成人av免费视频| 国产伦精品一区二区三区视频9 | 欧美黑人欧美精品刺激| 久久精品亚洲精品国产色婷小说| 欧美乱妇无乱码| 亚洲成人精品中文字幕电影| 亚洲在线观看片| 免费观看人在逋| 亚洲一区二区三区色噜噜| 国产v大片淫在线免费观看| 长腿黑丝高跟| 老司机在亚洲福利影院| 99热6这里只有精品| 国产免费av片在线观看野外av| 真人做人爱边吃奶动态| 欧美成人a在线观看| 男人舔奶头视频| 日韩欧美国产在线观看| 制服丝袜大香蕉在线| 我的老师免费观看完整版| 校园春色视频在线观看| 久久香蕉国产精品| 在线a可以看的网站| 精品无人区乱码1区二区| 高清日韩中文字幕在线| 熟妇人妻久久中文字幕3abv| 久久久国产成人免费| 日本 av在线| 免费av观看视频| 亚洲五月婷婷丁香| 亚洲成人中文字幕在线播放| 黄色丝袜av网址大全| 亚洲熟妇中文字幕五十中出| 日韩欧美在线乱码| 日本一二三区视频观看| 亚洲精品日韩av片在线观看 | 人妻夜夜爽99麻豆av| 成熟少妇高潮喷水视频| 九九热线精品视视频播放| 亚洲成人免费电影在线观看| 99久久综合精品五月天人人| 欧美一区二区国产精品久久精品| 国产视频一区二区在线看| 成年女人毛片免费观看观看9| 一二三四社区在线视频社区8| 亚洲色图av天堂| e午夜精品久久久久久久| 欧美黄色片欧美黄色片| 亚洲va日本ⅴa欧美va伊人久久| 国产真实乱freesex| 日韩欧美一区二区三区在线观看| 久久精品国产综合久久久| 欧美成狂野欧美在线观看| 在线观看一区二区三区| 91在线观看av| 99在线视频只有这里精品首页| 男女做爰动态图高潮gif福利片| 91九色精品人成在线观看| 给我免费播放毛片高清在线观看| 男女下面进入的视频免费午夜| 9191精品国产免费久久| 国产精华一区二区三区| 一区福利在线观看| 亚洲无线观看免费| 男人和女人高潮做爰伦理| 热99在线观看视频| 国产伦一二天堂av在线观看| 成人av在线播放网站| 国产欧美日韩一区二区三| 成人永久免费在线观看视频| 在线观看免费午夜福利视频| 宅男免费午夜| 99久久九九国产精品国产免费| 日本撒尿小便嘘嘘汇集6|