• <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)
    日本五十路高清| 国产真人三级小视频在线观看| 天堂中文最新版在线下载| av天堂久久9| 亚洲精品在线观看二区| 一夜夜www| 美女福利国产在线| kizo精华| 国产单亲对白刺激| 午夜精品国产一区二区电影| 亚洲精品中文字幕一二三四区 | 久久99热这里只频精品6学生| 日韩视频在线欧美| 亚洲国产av新网站| 国产精品久久久久久精品古装| 纵有疾风起免费观看全集完整版| 黄色毛片三级朝国网站| 一区福利在线观看| 18禁美女被吸乳视频| 我的亚洲天堂| 午夜福利一区二区在线看| 久久久久久亚洲精品国产蜜桃av| 999久久久精品免费观看国产| 免费女性裸体啪啪无遮挡网站| 亚洲色图综合在线观看| 久久天堂一区二区三区四区| 久久久久视频综合| 日韩欧美国产一区二区入口| tocl精华| 黑人巨大精品欧美一区二区mp4| 狠狠婷婷综合久久久久久88av| 免费一级毛片在线播放高清视频 | 99精品在免费线老司机午夜| 国产不卡av网站在线观看| 成人三级做爰电影| 老熟女久久久| 久久久国产成人免费| 黄片播放在线免费| 亚洲第一青青草原| 欧美大码av| 在线永久观看黄色视频| 精品午夜福利视频在线观看一区 | 人妻一区二区av| 亚洲精品久久午夜乱码| 国产成人精品在线电影| 两性午夜刺激爽爽歪歪视频在线观看 | 久久午夜亚洲精品久久| 嫩草影视91久久| av网站在线播放免费| 亚洲,欧美精品.| 在线观看免费高清a一片| 肉色欧美久久久久久久蜜桃| 国产精品一区二区精品视频观看| 中文欧美无线码| 天天躁夜夜躁狠狠躁躁| 日韩中文字幕视频在线看片| 一进一出好大好爽视频| 欧美日韩av久久| 大型黄色视频在线免费观看| 热99国产精品久久久久久7| 两个人看的免费小视频| 99热网站在线观看| 免费人妻精品一区二区三区视频| 亚洲熟女精品中文字幕| 欧美变态另类bdsm刘玥| 在线观看舔阴道视频| 国产精品成人在线| 亚洲精品在线观看二区| 正在播放国产对白刺激| 99精国产麻豆久久婷婷| 亚洲三区欧美一区| 十分钟在线观看高清视频www| 另类亚洲欧美激情| 国产成人精品久久二区二区免费| 一本大道久久a久久精品| 91九色精品人成在线观看| 丰满迷人的少妇在线观看| 一本—道久久a久久精品蜜桃钙片| 久久精品亚洲av国产电影网| 中文字幕人妻丝袜制服| 精品免费久久久久久久清纯 | 国产在视频线精品| 欧美人与性动交α欧美精品济南到| 老司机深夜福利视频在线观看| 男女边摸边吃奶| 色在线成人网| 精品国内亚洲2022精品成人 | 热re99久久精品国产66热6| 成年女人毛片免费观看观看9 | 美女福利国产在线| 中国美女看黄片| 男女午夜视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲av第一区精品v没综合| 久久久久久人人人人人| 欧美国产精品一级二级三级| 99国产综合亚洲精品| 国产无遮挡羞羞视频在线观看| av又黄又爽大尺度在线免费看| 久久久久精品人妻al黑| 777久久人妻少妇嫩草av网站| 一区二区av电影网| 热99久久久久精品小说推荐| av电影中文网址| 久久人妻熟女aⅴ| 在线十欧美十亚洲十日本专区| 久久精品成人免费网站| 亚洲一区中文字幕在线| 色94色欧美一区二区| 国产视频一区二区在线看| 亚洲成国产人片在线观看| 捣出白浆h1v1| 久久国产精品影院| 亚洲性夜色夜夜综合| 纵有疾风起免费观看全集完整版| 天天躁夜夜躁狠狠躁躁| 欧美大码av| 高清黄色对白视频在线免费看| 日韩制服丝袜自拍偷拍| 9色porny在线观看| 亚洲欧美一区二区三区黑人| 曰老女人黄片| 十八禁人妻一区二区| 久久精品熟女亚洲av麻豆精品| 亚洲午夜理论影院| 波多野结衣一区麻豆| 国产免费av片在线观看野外av| 99国产精品一区二区三区| 极品人妻少妇av视频| 成年人免费黄色播放视频| 国产亚洲av高清不卡| 高清av免费在线| 丝袜美腿诱惑在线| 国产欧美日韩一区二区三区在线| 久久人妻av系列| 色精品久久人妻99蜜桃| 欧美av亚洲av综合av国产av| 99国产精品99久久久久| 美女主播在线视频| 91精品三级在线观看| 亚洲专区国产一区二区| 亚洲第一欧美日韩一区二区三区 | 色94色欧美一区二区| 12—13女人毛片做爰片一| 欧美久久黑人一区二区| 999久久久国产精品视频| 天堂俺去俺来也www色官网| 亚洲国产中文字幕在线视频| 亚洲av第一区精品v没综合| 91精品三级在线观看| 性高湖久久久久久久久免费观看| 69av精品久久久久久 | 高清在线国产一区| 亚洲午夜精品一区,二区,三区| videos熟女内射| 亚洲,欧美精品.| 亚洲欧美激情在线| 亚洲第一欧美日韩一区二区三区 | 王馨瑶露胸无遮挡在线观看| 成年人午夜在线观看视频| 国产亚洲欧美在线一区二区| 无人区码免费观看不卡 | 免费黄频网站在线观看国产| 国产成人精品无人区| 最新的欧美精品一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 高潮久久久久久久久久久不卡| 飞空精品影院首页| 亚洲男人天堂网一区| 最新的欧美精品一区二区| 多毛熟女@视频| a级毛片黄视频| 午夜福利视频在线观看免费| 99国产综合亚洲精品| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区国产精品乱码| 18禁观看日本| 色94色欧美一区二区| 久久久精品国产亚洲av高清涩受| 9色porny在线观看| 男人舔女人的私密视频| 久久久久精品人妻al黑| 一区在线观看完整版| 人妻久久中文字幕网| 在线观看舔阴道视频| 免费观看av网站的网址| 久久人妻福利社区极品人妻图片| 免费观看人在逋| 日韩人妻精品一区2区三区| e午夜精品久久久久久久| 久久人妻av系列| 欧美久久黑人一区二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲七黄色美女视频| 日韩大片免费观看网站| av网站在线播放免费| 久久久国产精品麻豆| 新久久久久国产一级毛片| 免费日韩欧美在线观看| 亚洲欧洲精品一区二区精品久久久| 大香蕉久久成人网| 在线天堂中文资源库| av天堂在线播放| 亚洲午夜精品一区,二区,三区| av天堂久久9| 免费在线观看黄色视频的| 国产单亲对白刺激| 久久中文字幕人妻熟女| av欧美777| 亚洲精品在线美女| 大片免费播放器 马上看| 久久久国产一区二区| 一区二区三区精品91| 亚洲精品国产精品久久久不卡| 色在线成人网| 亚洲熟女毛片儿| 色综合婷婷激情| 欧美另类亚洲清纯唯美| 久久精品91无色码中文字幕| 亚洲精品美女久久av网站| 精品亚洲成国产av| 757午夜福利合集在线观看| 精品第一国产精品| 日本wwww免费看| 黄频高清免费视频| 精品视频人人做人人爽| 黄色 视频免费看| 一区二区三区乱码不卡18| 美女高潮到喷水免费观看| 国产熟女午夜一区二区三区| 国产有黄有色有爽视频| 久久久欧美国产精品| 国产一区二区三区视频了| 无限看片的www在线观看| 男女高潮啪啪啪动态图| 久久香蕉激情| 欧美精品一区二区大全| 99国产精品免费福利视频| 色在线成人网| 不卡一级毛片| 免费日韩欧美在线观看| 亚洲 欧美一区二区三区| 精品久久蜜臀av无| 久久午夜亚洲精品久久| 最近最新中文字幕大全免费视频| 精品少妇黑人巨大在线播放| 国产野战对白在线观看| 亚洲熟女毛片儿| 久久精品亚洲熟妇少妇任你| 精品国产乱子伦一区二区三区| 精品欧美一区二区三区在线| 波多野结衣一区麻豆| 欧美精品一区二区大全| 国产一区二区三区在线臀色熟女 | 91成年电影在线观看| 国产欧美日韩一区二区三| 国产男靠女视频免费网站| 999久久久精品免费观看国产| 久久精品亚洲熟妇少妇任你| 啦啦啦在线免费观看视频4| 色尼玛亚洲综合影院| 日韩一卡2卡3卡4卡2021年| 国产精品秋霞免费鲁丝片| 国产成人av激情在线播放| 悠悠久久av| 精品视频人人做人人爽| 国产在线免费精品| 天天操日日干夜夜撸| av福利片在线| 午夜成年电影在线免费观看| 亚洲中文日韩欧美视频| 久久青草综合色| 精品第一国产精品| 男女高潮啪啪啪动态图| 人妻久久中文字幕网| 欧美黄色淫秽网站| 18禁黄网站禁片午夜丰满| 日本av免费视频播放| 这个男人来自地球电影免费观看| 高清黄色对白视频在线免费看| 国产精品自产拍在线观看55亚洲 | 自线自在国产av| 亚洲精品美女久久久久99蜜臀| 免费在线观看完整版高清| 电影成人av| aaaaa片日本免费| 久久热在线av| 国产老妇伦熟女老妇高清| 久久久久久久大尺度免费视频| 亚洲欧美色中文字幕在线| 在线天堂中文资源库| e午夜精品久久久久久久| 9热在线视频观看99| 手机成人av网站| 国产精品国产高清国产av | 久久久久精品国产欧美久久久| 国产男女内射视频| 久久久国产一区二区| 欧美一级毛片孕妇| 成人国语在线视频| 夜夜夜夜夜久久久久| 国产精品秋霞免费鲁丝片| 国产极品粉嫩免费观看在线| 国产欧美日韩一区二区精品| 日本av免费视频播放| 久久性视频一级片| a在线观看视频网站| 超色免费av| 又大又爽又粗| 国产成人影院久久av| 国产黄色免费在线视频| 最近最新免费中文字幕在线| 欧美精品人与动牲交sv欧美| 纵有疾风起免费观看全集完整版| 国产成人精品无人区| 丝袜美腿诱惑在线| 亚洲精品久久成人aⅴ小说| 精品国产亚洲在线| www日本在线高清视频| 色老头精品视频在线观看| 成年版毛片免费区| 一级黄色大片毛片| 建设人人有责人人尽责人人享有的| 欧美成人午夜精品| 在线观看www视频免费| 蜜桃国产av成人99| 老司机靠b影院| 男人舔女人的私密视频| 国产在线一区二区三区精| videos熟女内射| 天天影视国产精品| 国产精品香港三级国产av潘金莲| 日本欧美视频一区| 在线观看免费视频网站a站| 另类亚洲欧美激情| 美女高潮喷水抽搐中文字幕| 十八禁网站网址无遮挡| av有码第一页| 精品国产一区二区三区久久久樱花| 国产精品av久久久久免费| 亚洲av第一区精品v没综合| 亚洲黑人精品在线| 久久精品91无色码中文字幕| 国产精品偷伦视频观看了| 男人操女人黄网站| 18禁黄网站禁片午夜丰满| 国产成人免费观看mmmm| 人妻久久中文字幕网| 国产人伦9x9x在线观看| 高潮久久久久久久久久久不卡| 美女福利国产在线| 不卡av一区二区三区| 九色亚洲精品在线播放| 一区二区三区激情视频| 久久久久久久国产电影| 97在线人人人人妻| 男女无遮挡免费网站观看| 国产一卡二卡三卡精品| 19禁男女啪啪无遮挡网站| 久久中文字幕人妻熟女| 啪啪无遮挡十八禁网站| 国产黄频视频在线观看| 免费在线观看黄色视频的| 国产精品亚洲av一区麻豆| 久久久精品94久久精品| 午夜免费鲁丝| 香蕉久久夜色| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲av高清不卡| 男人舔女人的私密视频| 99re在线观看精品视频| 国产伦理片在线播放av一区| 国产免费现黄频在线看| 69av精品久久久久久 | 老司机深夜福利视频在线观看| 一二三四在线观看免费中文在| 999精品在线视频| 午夜福利视频在线观看免费| 久久天堂一区二区三区四区| 欧美日韩中文字幕国产精品一区二区三区 | 久久天躁狠狠躁夜夜2o2o| 在线观看免费午夜福利视频| 欧美一级毛片孕妇| 我的亚洲天堂| 91大片在线观看| 国产成人精品久久二区二区免费| 激情在线观看视频在线高清 | 亚洲av成人一区二区三| 九色亚洲精品在线播放| 久久久久久久国产电影| 国产精品香港三级国产av潘金莲| 精品人妻熟女毛片av久久网站| 亚洲人成77777在线视频| 欧美精品人与动牲交sv欧美| 精品熟女少妇八av免费久了| 人妻久久中文字幕网| 久久99热这里只频精品6学生| 老司机深夜福利视频在线观看| 激情视频va一区二区三区| 日韩成人在线观看一区二区三区| 精品乱码久久久久久99久播| 肉色欧美久久久久久久蜜桃| 高潮久久久久久久久久久不卡| 男女无遮挡免费网站观看| 最近最新免费中文字幕在线| 两个人看的免费小视频| 国产高清视频在线播放一区| 国产精品一区二区在线观看99| 日韩大片免费观看网站| 激情在线观看视频在线高清 | 麻豆成人av在线观看| 丝袜喷水一区| 国产一区二区激情短视频| 亚洲一区二区三区欧美精品| e午夜精品久久久久久久| 最新在线观看一区二区三区| 中文字幕色久视频| 国产成人免费观看mmmm| 大码成人一级视频| 成年女人毛片免费观看观看9 | 男女下面插进去视频免费观看| 欧美亚洲 丝袜 人妻 在线| 亚洲av成人一区二区三| 精品视频人人做人人爽| 电影成人av| 黑人欧美特级aaaaaa片| 国产精品98久久久久久宅男小说| 精品一品国产午夜福利视频| 纵有疾风起免费观看全集完整版| 国产激情久久老熟女| 国产一区二区三区视频了| av在线播放免费不卡| 免费高清在线观看日韩| 嫁个100分男人电影在线观看| 国产97色在线日韩免费| 国产在线精品亚洲第一网站| 高清欧美精品videossex| 在线 av 中文字幕| 国产高清视频在线播放一区| 国产欧美日韩一区二区三| 黄色怎么调成土黄色| 国产成人免费无遮挡视频| 一二三四在线观看免费中文在| 在线观看免费视频日本深夜| 日韩免费av在线播放| 久久久精品免费免费高清| 久久中文字幕一级| 国产国语露脸激情在线看| 国产高清国产精品国产三级| 国产老妇伦熟女老妇高清| 精品亚洲成a人片在线观看| 久久精品人人爽人人爽视色| 波多野结衣一区麻豆| 最近最新中文字幕大全电影3 | 大片免费播放器 马上看| 老熟妇仑乱视频hdxx| 天堂动漫精品| 青草久久国产| 久久国产精品人妻蜜桃| 一二三四在线观看免费中文在| 国产不卡一卡二| 最近最新中文字幕大全免费视频| 一本色道久久久久久精品综合| 丝袜美足系列| 久久精品国产a三级三级三级| 国产三级黄色录像| 99热国产这里只有精品6| 99精品在免费线老司机午夜| 老汉色∧v一级毛片| 人妻一区二区av| 午夜久久久在线观看| 欧美av亚洲av综合av国产av| 少妇的丰满在线观看| 欧美在线黄色| 国产成人影院久久av| aaaaa片日本免费| 丝瓜视频免费看黄片| 又大又爽又粗| 亚洲三区欧美一区| 一区二区三区精品91| 亚洲成人免费电影在线观看| 麻豆av在线久日| 亚洲专区国产一区二区| 国产精品久久久久久人妻精品电影 | 最近最新中文字幕大全免费视频| 亚洲午夜理论影院| 国产三级黄色录像| 亚洲精品av麻豆狂野| 三上悠亚av全集在线观看| videosex国产| 十分钟在线观看高清视频www| 亚洲专区字幕在线| 99re在线观看精品视频| 久久久久久久久久久久大奶| 国产高清videossex| 精品一区二区三区视频在线观看免费 | 美女国产高潮福利片在线看| tube8黄色片| 母亲3免费完整高清在线观看| 高清在线国产一区| 国产精品久久久久久精品古装| 亚洲精品粉嫩美女一区| 老司机午夜十八禁免费视频| 免费在线观看视频国产中文字幕亚洲| 国产在线一区二区三区精| 中文字幕另类日韩欧美亚洲嫩草| 精品久久蜜臀av无| 欧美日韩av久久| 免费在线观看完整版高清| 中文字幕人妻熟女乱码| 黑人操中国人逼视频| 在线观看一区二区三区激情| av网站免费在线观看视频| 久久影院123| 久久人人97超碰香蕉20202| 国产有黄有色有爽视频| 欧美另类亚洲清纯唯美| 亚洲中文字幕日韩| 国产三级黄色录像| 少妇的丰满在线观看| 一二三四在线观看免费中文在| 天堂动漫精品| 久久久久久亚洲精品国产蜜桃av| 老熟妇乱子伦视频在线观看| 国产伦理片在线播放av一区| 女警被强在线播放| 欧美 亚洲 国产 日韩一| 国产精品国产av在线观看| 成人黄色视频免费在线看| 一区二区三区乱码不卡18| 午夜福利一区二区在线看| 最黄视频免费看| 欧美国产精品va在线观看不卡| 宅男免费午夜| 国产一区二区在线观看av| 9热在线视频观看99| 日本黄色日本黄色录像| 成年女人毛片免费观看观看9 | 午夜精品久久久久久毛片777| 无限看片的www在线观看| 宅男免费午夜| 操美女的视频在线观看| 国产精品免费视频内射| 成人免费观看视频高清| 高清av免费在线| 欧美日韩亚洲高清精品| 国产有黄有色有爽视频| 成人18禁高潮啪啪吃奶动态图| 国产激情久久老熟女| 在线亚洲精品国产二区图片欧美| 久久精品国产亚洲av高清一级| 黄片播放在线免费| 欧美精品高潮呻吟av久久| 一区二区三区精品91| 成人永久免费在线观看视频 | 亚洲精品一卡2卡三卡4卡5卡| 久久精品91无色码中文字幕| 999精品在线视频| 国产欧美日韩一区二区精品| 动漫黄色视频在线观看| 日韩大码丰满熟妇| 国产精品免费视频内射| 欧美亚洲 丝袜 人妻 在线| 少妇粗大呻吟视频| 国产极品粉嫩免费观看在线| 女人高潮潮喷娇喘18禁视频| 色综合欧美亚洲国产小说| 久久99热这里只频精品6学生| 日韩欧美一区二区三区在线观看 | 亚洲黑人精品在线| 黑人巨大精品欧美一区二区mp4| 日本黄色视频三级网站网址 | 成年人免费黄色播放视频| 精品国内亚洲2022精品成人 | 久久99热这里只频精品6学生| 一区福利在线观看| 亚洲免费av在线视频| 老司机靠b影院| 亚洲av国产av综合av卡| 亚洲自偷自拍图片 自拍| 成年人免费黄色播放视频| www.熟女人妻精品国产| 精品亚洲成国产av| 美女主播在线视频| 久久国产精品影院| 国产伦人伦偷精品视频| 美女福利国产在线| 天堂动漫精品| 正在播放国产对白刺激| 国产真人三级小视频在线观看| 国产麻豆69| 国产高清视频在线播放一区| 十八禁人妻一区二区| 老司机在亚洲福利影院| 久久人妻福利社区极品人妻图片| 高清欧美精品videossex| 色精品久久人妻99蜜桃| 妹子高潮喷水视频| 18禁黄网站禁片午夜丰满| videos熟女内射| 大码成人一级视频| 一区二区日韩欧美中文字幕| 色尼玛亚洲综合影院| 宅男免费午夜| 精品国产乱子伦一区二区三区| 日本黄色日本黄色录像| 免费一级毛片在线播放高清视频 | 肉色欧美久久久久久久蜜桃| 在线观看舔阴道视频| 国产精品一区二区在线观看99| 啦啦啦在线免费观看视频4| 亚洲久久久国产精品| 国产日韩一区二区三区精品不卡| 久久久精品免费免费高清| 一进一出好大好爽视频|