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

    威百畝熏蒸后土壤微生物群落重建及功能恢復(fù)

    2021-06-29 00:15趙鵬宇燕平梅趙曉東白雪
    植物保護(hù) 2021年3期

    趙鵬宇 燕平梅 趙曉東 白雪

    摘要 :本研究旨在科學(xué)評價威百畝熏蒸對土壤微生物生態(tài)系統(tǒng)的影響,為環(huán)境友好型消毒劑的選擇提供理論線索。本研究以實驗室威百畝熏蒸的土壤為材料,通過高通量測序技術(shù),研究威百畝熏蒸對土壤細(xì)菌群落的影響及熏蒸后土壤微生物群落重建及功能修復(fù)的機(jī)制。結(jié)果表明:高通量測序總共測得1 062 241個高質(zhì)量序列,共比對出5 882個細(xì)菌分類單元(OTUs)。熏蒸處理后不同階段土壤細(xì)菌群落結(jié)構(gòu)在門水平上較為類似,但是在種水平差異顯著。威百畝熏蒸處理對微生物群落結(jié)構(gòu)和功能造成很大影響,使得細(xì)菌群落α多樣性、群落構(gòu)建的驅(qū)動因素、共發(fā)生網(wǎng)絡(luò)關(guān)系以及碳水化合物代謝、氨基酸代謝、能量代謝等代謝活性在熏蒸處理后的初、中期階段表現(xiàn)出上升或者下降的趨勢,并在處理后的中后期迅速恢復(fù)并趨于平穩(wěn)。以上結(jié)果表明,威百畝處理對細(xì)菌群落造成強烈的“生態(tài)擾動”,導(dǎo)致細(xì)菌群落表現(xiàn)出“抑制激活恢復(fù)”的階段性特征。微生物群落的穩(wěn)定性與群落的抵抗力和恢復(fù)力密切相關(guān)。本研究為科學(xué)評價威百畝土壤消毒對土壤微生態(tài)系統(tǒng)的影響提供理論支撐,為環(huán)境友好型消毒劑的選擇提供線索。

    關(guān)鍵詞 :土壤消毒; 威百畝; 土壤微生物群落; 功能重建

    中圖分類號:

    Q 939.96

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

    DOI: 10.16688/j.zwbh.2020461

    Reconstruction and functional recovery of soil microbial community

    after fumigation of metam-sodium

    ZHAO Pengyu*, YAN Pingmei, ZHAO Xiaodong, BAI Xue

    (Department of Biology, Taiyuan Normal University, Jinzhong 030619, China)

    Abstract

    The purpose of this study is to scientifically evaluate the impact of metam-sodium fumigation on the soil micro-ecosystem and provide theoretical clues for the selection of environment-friendly disinfectants. The influence of metam-sodium fumigation on soil bacterial community and the mechanism of reconstruction and functional restoration of soil microbial community after fumigation were studied in the laboratory through high-throughput sequencing technology using metam-sodium fumigated soil as the material. A total of 1 062 241 high-quality sequences were obtained by high-throughput sequencing,among which a total of 5 882 bacterial taxa (OTUs) were identified. The soil bacterial community structures were similar at the phylum level at different stages after fumigation treatment,but show significant difference at the species level. Metam-sodium fumigation had a great impact on the structure and function of the microbial community, the alpha diversity of bacterial community, the driving factors of community construction, the co-occurring network relationship and the metabolic activity showed an increase or a downward trend at the initial and middle stages after metam-sodium fumigation, but quickly recovered and stabilized in the later period after the application of metam-sodium. The above results indicated that the application of metam-sodium caused a strong “ecological disturbance” to the bacterial community, resulting in phase characteristics of “inhibition-activation-recovery” in the bacterial community. The stability of microbial community was closely related to the resistance and resilience of the community. This study provides a theoretical support for the scientific evaluation of the impact of soil disinfection with metam-sodium on the soil micro-ecosystem and provides clues for the selection of environment-friendly disinfectants.

    Key words

    soil disinfection; metam-sodium; soil microbial community; functional reconstruction

    土壤消毒可以有效控制作物連作條件下的土傳病害和地下害蟲,是保護(hù)農(nóng)田種植行之有效的生產(chǎn)技術(shù)。土壤熏蒸劑是利用其揮發(fā)時產(chǎn)生的氣體來毒殺有害生物的一類農(nóng)藥。然而,由于熏蒸劑作用效果的廣譜性,在施用過程中不可避免地會對土壤中非目標(biāo)微生物的組成及功能活性造成影響,進(jìn)而改變整個土壤生態(tài)系統(tǒng)的結(jié)構(gòu)和功能,尤其是養(yǎng)分轉(zhuǎn)化過程[1]。與其他二硫代氨基甲酸鹽農(nóng)藥一樣,威百畝通過干擾病原菌呼吸酶(如丙酮酸脫氫酶)來發(fā)揮其生物活性。異硫氰酸甲酯是威百畝的降解產(chǎn)物,也是熏蒸過程的關(guān)鍵生物活性物質(zhì)[25]。從施用效果上看,威百畝熏蒸對土傳病害、地下害蟲及雜草的防治效果能夠替代溴甲烷,但其是否是優(yōu)良的環(huán)境友好型熏蒸劑目前尚未見報道。

    土壤熏蒸引起的土壤微生物群落結(jié)構(gòu)的破壞可以視為一種生態(tài)上的“擾動”。土壤熏蒸對土壤生態(tài)系統(tǒng)“擾動”的程度與多種因素有關(guān),主要包括熏蒸劑種類、劑型與施用量[6]、本底微生物群落組成[7]、土壤施肥歷史、土壤理化因子[8]及其質(zhì)地[9]等,并且與施肥策略[10]也有很大的關(guān)系。一方面,威百畝處理造成的生態(tài)“擾動”對土壤微生物群落造成強烈的環(huán)境脅迫,是研究農(nóng)藥施用對土壤微生物群落影響的理想實驗?zāi)P?另一方面,這種特殊的“擾動”可以用于設(shè)計微宇宙試驗,研究“擾動”后土壤微生物生態(tài)系統(tǒng)的重建過程。因此,熏蒸可以作為一種研究手段來深入探討土壤生態(tài)系統(tǒng)的“擾動抗逆恢復(fù)平衡”過程,研究“擾動”對隨后的群落重建和功能恢復(fù)的影響機(jī)制。

    本研究以實驗室威百畝熏蒸處理的土壤為試驗材料,通過高通量測序技術(shù),探究威百畝熏蒸對土壤細(xì)菌群落的影響及 “擾動”后土壤微生物群落重建及功能恢復(fù)情況,以期科學(xué)評價威百畝土壤熏蒸對土壤微生態(tài)系統(tǒng)的影響,為環(huán)境友好型消毒劑的選擇提供理論線索;并為土壤生態(tài)系統(tǒng)的維穩(wěn)機(jī)制提供一定的理論支持。

    1 材料和方法

    1.1 材料

    供試土壤為山西省太原市小店區(qū)黃瓜大棚連作3年以上的土壤。土壤為沙壤土,其理化性質(zhì)為:堿解氮50.08 mg/kg,有機(jī)質(zhì)27.85 g/kg,全氮1.27 g/kg,速效鉀223 mg/kg,有效磷43.96 mg/kg,pH 6.98。

    藥劑:42%威百畝水劑(AS),大連染料化工有限公司。

    1.2 微宇宙的設(shè)置以及取樣設(shè)置

    微宇宙試驗在一個含有300 g土的廣口瓶中進(jìn)行。42%威百畝AS的劑量設(shè)置為70 mg/kg(參考文獻(xiàn)[11],略有改動),共設(shè)置6個重復(fù),以加入等量的蒸餾水為對照,對照設(shè)置3個重復(fù)。試驗時先用滅菌蒸餾水調(diào)整土壤濕度至田間最大持水量的60%,之后將藥劑均勻施加到土壤中,立即使用瓶塞和凡士林封閉瓶口,然后放置在25℃的黑暗條件中恒溫培養(yǎng)7 d。7 d后移除瓶塞,瓶口用透氣塑料膜再次封好,以避免空氣中的雜菌污染土壤。

    移除瓶塞后在微生物的重建過程中定期加入滅菌水來彌補土壤中水分的損失。分別于移除瓶塞后3 d(本研究簡寫為F),15 d(S),40 d(T)從每個廣口瓶中取1 g土壤樣本,盡量保持土壤完整,避免破壞性取樣。部分用于分析土壤物理化學(xué)參數(shù);剩余部分儲存在-80℃冰箱中用于后續(xù)核酸的提取。

    1.3 DNA提取、PCR擴(kuò)增和測序

    使用E.Z.N.A.Soil DNA Kit(Omega,USA)提取土壤DNA。用酶標(biāo)儀Infinite200 PRO(TECAN,Switzerland)測定DNA的純度和濃度。DNA的A260/A280在1.8~2.0時可用于下一步試驗。

    擴(kuò)增體系(25 μL):5 U/μL Plantium Taq 0.5 μL,10×EasyTaq Buffer(Mg2+ Plus)5.0 μL,dNTPs 40 μL,50 μmol/L引物0.5 μL,DNA模板5.0 μL,加入無菌水至25 μL。PCR反應(yīng)條件:95℃ 3 min;95℃ 30 s,62℃30 s,72℃ 30 s,30個循環(huán);72℃ 5 min。

    從每個土壤樣品中提取3個DNA樣品,將其PCR產(chǎn)物充分混勻后送上海美吉生物技術(shù)工程有限公司,通過Illumina Miseq測序平臺,對16S rDNA高變區(qū)V3-V4區(qū)進(jìn)行測序。

    1.4 群落構(gòu)建過程分析

    為了觀察群落構(gòu)建過程的驅(qū)動因素,參考Stegen等[1214]提出一種基于系統(tǒng)發(fā)育結(jié)構(gòu)的分析方法(圖1)。

    1.5 統(tǒng)計分析

    采用單因素方差分析(ANOVA)評估威百畝熏蒸處理后不同階段細(xì)菌群落α多樣性、優(yōu)勢類群、群落構(gòu)建的驅(qū)動因素、共發(fā)生網(wǎng)絡(luò)關(guān)系以及功能預(yù)測的差異。利用非度量多維尺度分析(NMDS)對熏蒸處理后不同階段的細(xì)菌群落進(jìn)行排序分析。共發(fā)生網(wǎng)絡(luò)分析是基于OTU豐度,利用斯皮爾曼等級相關(guān)矩陣(Spearman correlation matrix)評估變量之間的相關(guān)性,篩選出相關(guān)系數(shù)大于0.6且呈顯著相關(guān)的OTUs(否則,剔除)。然后采用false discovery rate(FDR<0.05)校正P值。網(wǎng)絡(luò)分析的參數(shù)運算使用R語言中igragh程序包,并且在Gepih軟件中可視化(https:∥gephi.org/)。物種的共存格局(co-occurrence pattern)是由嵌套程度和離散程度共同決定,因此使用WNODF指數(shù)來評估群落中物種的嵌套程度(degree of nestedness);使用C值來判定群落的聚集/分離(aggregation/segregation)程度。使用PICRUSt2針對樣本中微生物群落的功能組成進(jìn)行功能預(yù)測。

    2 結(jié)果與分析

    2.1 威百畝熏蒸對土壤細(xì)菌群落結(jié)構(gòu)和多樣性的影響

    高通量測序總共獲得1 062 241個高質(zhì)量序列,比對出5 882個細(xì)菌分類單元(OTUs)。測序的覆蓋度指數(shù)為99.19%~99.75%,這表明所比對的序列能夠涵蓋絕大部分土壤樣品中的細(xì)菌序列。

    共有3 439個OTUs(58.46%)是在威百畝熏蒸處理后不同階段(CK、處理后3、15 d和 40 d)所共有的(圖2)。其中,14個(0.23%)、29個(0.49%)、63個(1.07%)和144個(2.44%)OTUs分別只在CK、處理后3、15 d和40 d中檢測出。并且,CK組的OTUs數(shù)相對最小,而處理后40 d(T階段)的OTUs相對最多。

    本研究相對豐度大于0.01%的細(xì)菌門總共比對出13個(圖3)。威百畝熏蒸處理后不同階段,放線菌門Actinobacteria的相對豐度最高(平均相對豐度32.05%),其次是變形菌門Proteobacteria(2224%)和綠彎菌門Chloroflexi(15.80%)。其中,放線菌門在熏蒸處理后3 d的相對豐度最高(33.68%),在CK組的相對豐度最低(26.35%);變形菌門在熏蒸處理后3 d的相對豐度最高(2398%),在15 d的相對豐度最低(21.97%);綠彎菌門在熏蒸處理后40 d的相對豐度最高(1718%),在15 d的相對豐度最低(14.33%)。并且,如圖4所示,放線菌門和Bacteroidota菌群在處理后不同階段的相對豐度有顯著差異(P<0.05),而其他菌群的相對豐度在處理后不同階段無顯著差異(P>0.05)。

    使用非度量多維尺度分析(NMDS)對熏蒸處理后不同階段的處理組與對照組細(xì)菌群落結(jié)構(gòu)(種水平)進(jìn)行排序分析(圖5)。結(jié)果表明,土壤細(xì)菌群落結(jié)構(gòu)在處理后相同階段內(nèi)相似程度高,而在不同階

    段差異較大。組間差異顯著性分析(ADONIS, P<0.01)表明,微生物群落結(jié)構(gòu)在不同階段有顯著

    差異。綜上,在熏蒸處理后不同階段,細(xì)菌群落結(jié)構(gòu)在門水平上較為類似,但是在種水平上差異顯著。

    2.2 熏蒸處理后不同階段的細(xì)菌群落多樣性及其多樣性維持機(jī)制

    如圖6所示,威百畝熏蒸處理后,土壤微生物群落的多樣性受到很大影響。Sobs指數(shù)、ACE指數(shù)和Chao指數(shù)在CK組分別為3 104、4 158和4 200,處理后3 d分別降低為2 678、3 780和3 840;處理后15 d分別升高為3 105、4 144和4 134,且處理后3 d和15 d存在顯著差異(P<0.05);處理后40 d Sobs指數(shù)(3 246)、ACE指數(shù)(4 235)和Chao指數(shù)(4 228)和15 d時基本持平,兩者間無顯著差異(P>0.05)。另一方面,Shannon指數(shù)和Simpson指數(shù),在CK、處理后3、15 d和40 d均沒有顯著差異(P>0.05)。綜上,土壤微生物多樣性(尤其是豐富度)在處理后不同階段呈現(xiàn)出先下降后上升再平穩(wěn)的階段性特征。

    群落多樣性的維持機(jī)制就是群落構(gòu)建的驅(qū)動機(jī)制,即群落構(gòu)建機(jī)制。根據(jù)熏蒸處理后不同階段土壤細(xì)菌群落結(jié)構(gòu)betaNTI值的結(jié)果(圖7),參考Stegen等的群落構(gòu)建劃分方法[12],推斷群落構(gòu)建過程在CK、處理后3 d和40 d更多地受到變量選擇的影響(heterogeneous selection),其betaNTI值(箱線圖中的中位數(shù))大于2(圖7)。并且,在處理后15 d,betaNTI值為1.85,表明細(xì)菌群落在處理后15 d更多地受到隨機(jī)過程的驅(qū)動(|betaNTI|<2)。

    2.3 熏蒸處理后不同階段土壤細(xì)菌群落共發(fā)生網(wǎng)絡(luò)分析

    細(xì)菌群落共發(fā)生網(wǎng)絡(luò)分析(co-occurrence network)主要反映某一環(huán)境條件下微生物群落的共存格局,可以在一定程度上反映群落的種間關(guān)系。通過微生物群落共發(fā)生網(wǎng)絡(luò)分析,相比于熏蒸處理后各階段,CK組的邊數(shù)量(1 648)、節(jié)點數(shù)量(232)、聚集系數(shù)(0.51)、模塊數(shù)(6)和平均度(14.16)都是最小的;而在熏蒸處理后3 d,節(jié)點數(shù)量(282)、聚集系數(shù)(0.56)和模塊數(shù)(16)都是最大的;15 d和40 d的節(jié)點數(shù)量、聚集系數(shù)和模塊數(shù)處于兩者之間。這表明,細(xì)菌群落的共發(fā)生網(wǎng)絡(luò)關(guān)系也具有一定的階段性特征。

    通常將除外圍節(jié)點(peripherals)外的其余3種類型的節(jié)點(connectors、module hubs和network hubs)歸為關(guān)鍵節(jié)點,而網(wǎng)絡(luò)關(guān)系中的關(guān)鍵節(jié)點是群落中起關(guān)鍵作用的核心物種[15]。本研究通過微生物共發(fā)生網(wǎng)絡(luò)的模塊內(nèi)連通度(Zi)和模塊間連通度(Pi)來鑒定核心物種(圖9),結(jié)果表明,熏蒸處理后3 d的核心物種數(shù)量更多,而CK組的核心物種數(shù)量較少。這表明,群落的核心物種在熏蒸后呈現(xiàn)出先上升后下降的階段性特征。

    為了揭示細(xì)菌群落的共存格局(co-occurrence patterns),本研究計算了群落中物種的嵌套指數(shù)WNODF值與C值(表1)。結(jié)果表明,處理后3 d的WNODF值最大(49.23),CK組的WNODF值最小(38.42),這表明群落在處理后3 d的嵌套結(jié)構(gòu)更明顯。另一方面,由于處理后3 d的C值最高(043,P<0.01),表明細(xì)菌群落在處理后3 d具有更顯著的離散格局(segregated pattern)。綜上,細(xì)菌群落的共存格局在處理后也有一定的階段性特征。

    2.4 熏蒸處理后不同階段的代謝活性分析

    通過PICRUSt2對KEGG Pathway level 2進(jìn)行單因素方差分析(表2),篩選出代謝活性前15的微生物群落代謝通路,其中碳水化合物代謝、氨基酸代謝、能量代謝、輔助因子和維生素的代謝、膜運輸、核苷酸代謝、翻譯、脂質(zhì)代謝、復(fù)制和修復(fù)、信號轉(zhuǎn)導(dǎo)、異種生物的生物降解和代謝、降解等在3 d的活性都是最小的。但是,在15 d和40 d又會顯著升高,且和CK組的基本持平。進(jìn)一步推斷可知,細(xì)菌群落的代謝活性在熏蒸處理后不同階段呈現(xiàn)出先下降再上升的階段性特征。

    3 討論

    當(dāng)前我國針對農(nóng)田種植正大力倡導(dǎo)減肥減藥,土壤熏蒸技術(shù)是行之有效的解決手段。然而,土壤熏蒸在極大程度上殺死病原微生物的同時,也不可避免地會顯著改變土壤細(xì)菌[1617]、真菌[18]及線蟲[19]的群落組成,對土壤生態(tài)系統(tǒng)造成“擾動”。

    威百畝作為溴甲烷的替代品之一,因其具有良好的土壤熏蒸效果而在農(nóng)業(yè)領(lǐng)域得到了廣泛應(yīng)用[20]。威百畝熏蒸會對土壤的生物地球化學(xué)循環(huán)、土壤生態(tài)系統(tǒng)產(chǎn)生深遠(yuǎn)影響[11]。本研究結(jié)果表明,威百畝熏蒸處理后不同階段的土壤細(xì)菌群落結(jié)構(gòu)在門水平上較為類似,但是在種水平上差異顯著。Li等[21]研究發(fā)現(xiàn),威百畝熏蒸后土壤中細(xì)菌群落和優(yōu)勢菌群發(fā)生了改變,但是在門水平上樣本之間差異卻不明顯。本研究中,土壤細(xì)菌群落的多樣性(尤其是豐富度)在熏蒸處理后3 d顯著降低,但是在熏蒸處理后15 d和40 d又顯著升高,表現(xiàn)出先下降后上升的階段性特征,這種變化趨勢與卜東欣等[22]的研究結(jié)果類似。這可能是因為威百畝對土壤細(xì)菌群落產(chǎn)生強烈的毒性效應(yīng),可殺死大部分的細(xì)菌類群。但是,少數(shù)耐藥性較強的菌群可存活下來。另外,威百畝打破了原來土壤微生物群落的平衡,使微生物群落形成一個短暫的“生態(tài)真空期”,為這些耐藥性強的細(xì)菌菌群提供了豐富的營養(yǎng)物質(zhì)和生存空間,使其迅速增殖達(dá)到一定的平衡。這就解釋了本研究細(xì)菌群落在熏蒸處理后具有階段性特征的原因。土壤細(xì)菌群落對這種強烈的 “生態(tài)擾動”表現(xiàn)出“抑制激活恢復(fù)”的趨勢[22]。

    本研究得出的群落構(gòu)建機(jī)制也從另一方面佐證了細(xì)菌群落在熏蒸處理后表現(xiàn)出的“抑制激活恢復(fù)”的階段性特征。施用威百畝而造成的“擾動”會導(dǎo)致土壤理化因子和質(zhì)地的劇烈變化,增加了環(huán)境異質(zhì)性,導(dǎo)致細(xì)菌群落的相異性增加,這就是變量選擇對群落的驅(qū)動作用[1214]。本研究結(jié)果也表明,在熏蒸處理后第3天群落的結(jié)構(gòu)、多樣性、共發(fā)生網(wǎng)絡(luò)關(guān)系以及功能屬性等變化都是最大的。而在熏蒸處理后第15天由于耐藥性菌群的大量增殖,隨機(jī)過程(擴(kuò)散、漂變和多樣化)可能占據(jù)主動地位。而在熏蒸處理后40 d,細(xì)菌群落達(dá)到平穩(wěn)狀態(tài),這時又會受到環(huán)境異質(zhì)性以及種內(nèi)排斥等的限制,變量選擇對群落的驅(qū)動作用增加。

    本研究結(jié)果表明:熏蒸處理后3 d時節(jié)點數(shù)量、聚集系數(shù)和模塊數(shù)都是最大的。在微生物共發(fā)生網(wǎng)絡(luò)分析中,模塊化程度越高(modularity),網(wǎng)絡(luò)越密集,表明微生物群落越復(fù)雜[23]。這可能是因為細(xì)菌群落在熏蒸處理后3 d受到變量選擇的驅(qū)動作用,群落的相異性增加,因而此階段的群落共發(fā)生網(wǎng)絡(luò)關(guān)系也就越復(fù)雜。另外,菌群在熏蒸處理后第3天表現(xiàn)出更加顯著的嵌套結(jié)構(gòu)(WNDF值最大)和離散格局(C值最大,且P<0.01),這可能也是變量選擇驅(qū)動的結(jié)果。在隨后的15 d和40 d,群落的節(jié)點數(shù)量、聚集系數(shù)和模塊數(shù)略微降低,可能是變量選擇的驅(qū)動作用減小,群落的相似性增加的結(jié)果。但是,共發(fā)生網(wǎng)絡(luò)關(guān)系的變化趨勢并沒有微生物群落結(jié)構(gòu)、多樣性和功能活性變化那么明顯,可能是因為微生物群落中物種數(shù)極多,代際周期短,并且存在一定程度的功能冗余。

    微生物群落的結(jié)構(gòu)組成可以在時間或者空間上展現(xiàn)出復(fù)雜的變化,例如,在不同的海水層[24]或者亞高山森林的海拔梯度上[25],這些變化最終會影響其所執(zhí)行的群落功能[26]。本研究結(jié)果表明,威百畝熏蒸處理后微生物群落代謝活性呈現(xiàn)出先下降后上升的特征,也表現(xiàn)出“抑制激活恢復(fù)”的階段性特征。并且群落功能活性的變化趨勢也與群落結(jié)構(gòu)的變化大致相同。

    本研究中微生物群落在熏蒸處理后表現(xiàn)出的階段性特征與微生物生態(tài)系統(tǒng)的穩(wěn)定性密切相關(guān)。而微生物生態(tài)系統(tǒng)的穩(wěn)定性又與微生物群落的抵抗力(resistance)和恢復(fù)力(resilience)密切相關(guān)。抵抗力指微生物群落自身具備的抵抗自然或者環(huán)境因子脅迫并維持原來狀態(tài)不變的能力。土壤經(jīng)過威百畝熏蒸處理后,遭受強烈的“生態(tài)擾動”,微生物群落的抵抗力被打破,微生物菌群受到強烈的脅迫。脅迫帶來的選擇壓力可以促進(jìn)初始群落迅速進(jìn)化,而這種迅速進(jìn)化是微生物群落的結(jié)構(gòu)和功能對特定脅迫尤其是長期脅迫產(chǎn)生抵抗力和恢復(fù)力的重要機(jī)制[27],帶來的直接效應(yīng)就是群落中的土壤細(xì)菌群落豐富度、多樣性迅速減小。另一方面,恢復(fù)力是指系統(tǒng)受到自然或人為環(huán)境因子的脅迫后恢復(fù)至原有狀態(tài)的能力[28]。微生物群落較強的適應(yīng)性能夠保證其在變化的環(huán)境中迅速恢復(fù)穩(wěn)定。有研究指出真菌群落和叢枝菌根真菌在不同的環(huán)境條件和農(nóng)業(yè)管理措施下,對水分限制都表現(xiàn)出顯著的抵抗力和恢復(fù)力[29]。

    4 結(jié)論

    本研究以實驗室威百畝熏蒸處理的土壤為材料,通過高通量測序技術(shù),研究威百畝熏蒸處理對土壤細(xì)菌群落結(jié)構(gòu)和功能的影響,也探究了“擾動”后土壤微生物群落重建及功能恢復(fù)的機(jī)制。結(jié)果表明,在威百畝處理后,微生物群落多樣性、結(jié)構(gòu)、群落構(gòu)建機(jī)制、種間關(guān)系和功能都表現(xiàn)出“抑制激活恢復(fù)”的階段性特征。

    參考文獻(xiàn)

    [1] 李青杰, 方文生, 顏冬冬,等. 熏蒸劑對土壤微生物的影響研究進(jìn)展[J]. 農(nóng)藥學(xué)學(xué)報, 2019, 21(S1): 780786.

    [2] 胡洪濤, 朱志剛, 焦忠久,等. 棉隆土壤消毒對高山甘藍(lán)根腫病土壤細(xì)菌群落結(jié)構(gòu)的影響[J]. 中國農(nóng)學(xué)通報, 2020, 36(16): 120127.

    [3] 燕平梅, 李娜, 趙文婧,等. 威百畝熏蒸對土壤硝化作用及amoA型硝化細(xì)菌群落結(jié)構(gòu)的影響[J]. 植物保護(hù)學(xué)報, 2019, 46(6): 13571363.

    [4] 王峰, 呂藝, 劉洋,等. 氯化苦土壤消毒對三七種植后細(xì)菌群落多樣性及烤煙輪作的影響[J]. 昆明理工大學(xué)學(xué)報(自然科學(xué)版), 2019, 44(4): 9096.

    [5] FANG Wensheng, WANG Qiuxia, HAN Dawei, et al. The effects and mode of action of biochar on the degradation of methyl isothiocyanate in soil [J]. Science of the Total Environment, 2016, 565(4): 339345.

    [6] 曹坳程, 張文吉, 劉建華. 溴甲烷土壤消毒替代技術(shù)研究進(jìn)展[J]. 植物保護(hù), 2007,33(1): 1520.

    [7] DICKENS H E, ANDERSON J M. Manipulation of soil microbial community structure in bog and forest soils using chloroform fumigation [J]. Soil Biology & Biochemistry, 1999, 31(14): 20492058.

    [8] DUNGAN R S, IBEKWE A M, YATES S R. Effect of propargyl bromide and 1,3-dichloropropene on microbial communities in an organically amended soil [J]. FEMS Microbiology Ecology, 2003, 43(1): 7587.

    [9] COLLINS H P, ALVA A, BOYDSTON R A, et al. Soil microbial, fungal, and nematode responses to soil fumigation and cover crops under potato production [J]. Biology and Fertility of Soils, 2006, 42(3): 247257.

    [10]GRIFFITHS B S, HALLETT P D, KUAN H L, et al. Biological and physical resilience of soil amended with heavy metal-contaminated sewage sludge [J]. European Journal of Soil Science, 2010, 56(2): 197206.

    [11]燕平梅, 魏愛麗, 喬宏萍,等. 威百畝熏蒸對土壤反硝化作用及nirS型反硝化細(xì)菌群落結(jié)構(gòu)的影響[J]. 植物保護(hù)學(xué)報, 2019, 46(2): 393400.

    [12]STEGEN J C, LIN Xueju, FREDRICKSON J K, et al. Estimating and mapping ecological processes influencing microbial community assembly[J/OL]. Frontiers in Microbiology, 2015, 6: 370. DOI:10.3389/fmicb.2015.00370.

    [13]STEGEN J C, LIN Xueju, KONOPKA A E, et al. Stochastic and deterministic assembly processes in subsurface microbial communities [J]. ISME Journal, 2012, 6(9): 16531664.

    [14]STEGEN J C, LIN Xueju, FREDRICKSON J K, et al. Quantifying community assembly processes and identifying features that impose them [J].ISME Journal, 2013, 7(11): 20692079.

    [15]DENG Ye, JIANG Yihuai, YANG Yunfeng, et al. Molecular ecological network analyses[J/OL]. BMC Bioinformatics, 2012, 13(1): 113. DOI:10.1186/1471-2105-13-113.

    [16]程新勝, 楊建卿. 熏蒸處理對土壤微生物及硝化作用的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2007, 15(6): 5153.

    [17]HOSHINO Y T, MATSUMOTO N. DNA-versus RNA-based denaturing gradient gel electrophoresis profiles of a bacterial community during replenishment after soil fumigation [J]. Soil Biology & Biochemistry, 2007, 39(2): 434444.

    [18]HOSHINO Y T, MATSUMOTO N. Changes in fungal community structure in bulk soil and spinach rhizosphere soil after chemical fumigation as revealed by 18S rDNA PCR-DGGE [J]. Soil Science & Plant Nutrition, 2010, 53(1): 4055.

    [19]WANG K H, MCSORLEY R, KOKALIS-BURELLE N. Effects of cover cropping, solarization, and soil fumigation on nematode communities [J]. Plant & Soil, 2006, 286(2): 229243.

    [20]TRIKY-DOTAN S, AUSTERWEIL M, STEINER B, et al. Generation and dissipation of methyl isothiocyanate in soils following metam sodium fumigation: Impact on verticillium control and potato yield [J]. Plant Disease, 2007, 91(5): 497503.

    [21]LI Jun, HUANG Bin, WANG Qiuxia, et al. Effects of fumigation with metam-sodium on soil microbial biomass, respiration, nitrogen transformation, bacterial community diversity and genes encoding key enzymes involved in nitrogen cycling [J]. Science of the Total Environment, 2016, 598(1): 10271036.

    [22]卜東欣, 張超, 張鑫, 等. 熏蒸劑威百畝對土壤微生物數(shù)量和酶活性的影響[J]. 中國農(nóng)學(xué)通報, 2014, 30(15): 227233.

    [23]OLESEN J M. The modularity of pollination networks[J]. Proceedings of the National Academy of Sciences, 2007, 104(50): 1989119896.

    [24]BARBERN A, FERNNDEZ-GUERRA A, BOHANNAN B J M, et al. Exploration of community traits as ecological markers in microbial metagenomes [J]. Molecular Ecology, 2012, 21(8):19091917.

    [25]SILES J A, MARGESIN R. Abundance and diversity of bacterial, archaeal, and fungal communities along an altitudinal gradient in alpine forest soils: what are the driving factors?[J]. Microbial Ecology, 2016, 72(1): 207220.

    [26]LOUCA S, JACQUES S M S, PIRES A P F, et al. High taxonomic variability despite stable functional structure across microbial communities [J]. Nature Ecology & Evolution, 2016, 1(1): 1517.

    [27]COHAN F M. What are bacterial species?[J]. Annual Review of Microbiology, 2002, 56(1): 457465.

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

    [29]FURZE J R, MARTIN A R, NASIELSKI J, et al. Resistance and resilience of root fungal communities to water limitation in a temperate agroecosystem [J]. Ecology and Evolution, 2017, 7(10): 34433454.

    (責(zé)任編輯:楊明麗)

    精品久久久久久电影网| 丝袜喷水一区| 亚洲成色77777| 男男h啪啪无遮挡| 久久人人爽人人片av| 青春草视频在线免费观看| 国产成人精品福利久久| 3wmmmm亚洲av在线观看| 校园人妻丝袜中文字幕| 成人二区视频| 亚洲中文av在线| 亚洲熟女精品中文字幕| 91精品伊人久久大香线蕉| 另类精品久久| 国产欧美亚洲国产| 精品国产一区二区久久| 久久久久久久精品精品| 亚洲精品国产av蜜桃| 欧美+日韩+精品| 亚洲av欧美aⅴ国产| 各种免费的搞黄视频| videossex国产| 成年av动漫网址| 成年人免费黄色播放视频| 国产精品成人在线| 久久精品久久精品一区二区三区| 三上悠亚av全集在线观看| 亚洲性久久影院| 最新的欧美精品一区二区| √禁漫天堂资源中文www| 亚洲精品日本国产第一区| 久久久a久久爽久久v久久| 国产男女超爽视频在线观看| 亚洲国产av影院在线观看| 卡戴珊不雅视频在线播放| 嘟嘟电影网在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧美日韩卡通动漫| 国产深夜福利视频在线观看| 高清午夜精品一区二区三区| 一边亲一边摸免费视频| 亚洲五月色婷婷综合| 久久婷婷青草| 一级,二级,三级黄色视频| 亚洲欧美一区二区三区国产| 欧美激情 高清一区二区三区| 亚洲av中文av极速乱| 黑人猛操日本美女一级片| a 毛片基地| 成年女人在线观看亚洲视频| 18在线观看网站| 91久久精品国产一区二区三区| 亚洲欧美一区二区三区国产| 99re6热这里在线精品视频| av播播在线观看一区| av天堂久久9| 韩国av在线不卡| 国产不卡av网站在线观看| 一级,二级,三级黄色视频| a级毛片免费高清观看在线播放| 水蜜桃什么品种好| 久久久久久久久久成人| 亚洲天堂av无毛| 免费观看性生交大片5| 丝袜在线中文字幕| 亚洲人成77777在线视频| 国产伦精品一区二区三区视频9| 久久韩国三级中文字幕| 国产亚洲av片在线观看秒播厂| 美女脱内裤让男人舔精品视频| 99久久中文字幕三级久久日本| 一区二区三区乱码不卡18| 国产亚洲av片在线观看秒播厂| 高清午夜精品一区二区三区| 免费av中文字幕在线| 日韩 亚洲 欧美在线| 考比视频在线观看| 视频在线观看一区二区三区| 啦啦啦中文免费视频观看日本| 国产成人一区二区在线| 亚洲欧美成人综合另类久久久| 国产成人a∨麻豆精品| 丝瓜视频免费看黄片| 免费高清在线观看日韩| 免费久久久久久久精品成人欧美视频 | 少妇的逼好多水| 老司机亚洲免费影院| 日韩视频在线欧美| 国产精品偷伦视频观看了| 91久久精品国产一区二区三区| 建设人人有责人人尽责人人享有的| 久久亚洲国产成人精品v| 成人手机av| 人人妻人人添人人爽欧美一区卜| 十八禁高潮呻吟视频| 国产精品久久久久久久电影| 国产亚洲午夜精品一区二区久久| 另类精品久久| 国产有黄有色有爽视频| 男女国产视频网站| 18禁观看日本| 人体艺术视频欧美日本| 男女国产视频网站| 久热久热在线精品观看| .国产精品久久| 亚洲av.av天堂| 久久国内精品自在自线图片| 免费看不卡的av| 黄片无遮挡物在线观看| 能在线免费看毛片的网站| 永久网站在线| 黄片播放在线免费| 国产深夜福利视频在线观看| 欧美精品人与动牲交sv欧美| 国产av国产精品国产| 亚洲精品一区蜜桃| 91午夜精品亚洲一区二区三区| 亚洲情色 制服丝袜| 欧美亚洲日本最大视频资源| 肉色欧美久久久久久久蜜桃| 有码 亚洲区| 下体分泌物呈黄色| 国产欧美亚洲国产| 成人国产av品久久久| 精品国产露脸久久av麻豆| 一本色道久久久久久精品综合| 国产片特级美女逼逼视频| 国产av一区二区精品久久| 亚洲精品,欧美精品| 久久97久久精品| 免费黄频网站在线观看国产| 汤姆久久久久久久影院中文字幕| 国产成人免费无遮挡视频| 日日啪夜夜爽| 亚洲av中文av极速乱| 亚洲欧洲精品一区二区精品久久久 | a级毛色黄片| 赤兔流量卡办理| 精品国产一区二区三区久久久樱花| 国产精品人妻久久久影院| 熟女人妻精品中文字幕| 在线精品无人区一区二区三| 十分钟在线观看高清视频www| 精品一区二区三区视频在线| 十八禁网站网址无遮挡| 亚洲国产精品专区欧美| 少妇人妻精品综合一区二区| 免费av不卡在线播放| 国产男人的电影天堂91| 国产色婷婷99| 超色免费av| 免费大片18禁| 在线观看一区二区三区激情| 国产一区二区三区av在线| 毛片一级片免费看久久久久| 男人添女人高潮全过程视频| 观看av在线不卡| 国语对白做爰xxxⅹ性视频网站| 亚洲丝袜综合中文字幕| 69精品国产乱码久久久| 下体分泌物呈黄色| 久久女婷五月综合色啪小说| 国产亚洲最大av| 久久久久久久大尺度免费视频| 九九在线视频观看精品| 中文字幕最新亚洲高清| 老女人水多毛片| 久久久亚洲精品成人影院| 亚洲精品久久午夜乱码| 男人添女人高潮全过程视频| 搡女人真爽免费视频火全软件| 最后的刺客免费高清国语| 一区二区日韩欧美中文字幕 | 成人国语在线视频| 免费看av在线观看网站| 人妻夜夜爽99麻豆av| 亚洲国产精品一区三区| av在线老鸭窝| 五月伊人婷婷丁香| 欧美xxxx性猛交bbbb| 久久久久国产网址| 激情五月婷婷亚洲| 国产精品.久久久| 亚洲av福利一区| 另类精品久久| av网站免费在线观看视频| 色视频在线一区二区三区| 极品少妇高潮喷水抽搐| 制服诱惑二区| 五月天丁香电影| 在线精品无人区一区二区三| 免费久久久久久久精品成人欧美视频 | av女优亚洲男人天堂| 伦精品一区二区三区| 色吧在线观看| 欧美成人精品欧美一级黄| 你懂的网址亚洲精品在线观看| 久久午夜综合久久蜜桃| 久久国产精品男人的天堂亚洲 | 日本av手机在线免费观看| 乱码一卡2卡4卡精品| 这个男人来自地球电影免费观看 | 国产精品 国内视频| 亚洲欧美日韩卡通动漫| 日韩成人伦理影院| 又粗又硬又长又爽又黄的视频| 国产无遮挡羞羞视频在线观看| 熟女电影av网| 啦啦啦中文免费视频观看日本| 夜夜爽夜夜爽视频| 亚洲美女视频黄频| 国产一级毛片在线| 亚洲精品自拍成人| 国产精品熟女久久久久浪| 国产探花极品一区二区| 免费av中文字幕在线| 高清毛片免费看| 一级,二级,三级黄色视频| 国产精品久久久久久av不卡| 久久毛片免费看一区二区三区| 国产精品秋霞免费鲁丝片| 免费人成在线观看视频色| 观看av在线不卡| 欧美日韩av久久| 性高湖久久久久久久久免费观看| 在线观看美女被高潮喷水网站| 91成人精品电影| 日本av免费视频播放| 日本欧美国产在线视频| 亚洲av二区三区四区| av一本久久久久| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 在现免费观看毛片| 99国产精品免费福利视频| 国产永久视频网站| 国产精品人妻久久久影院| 午夜激情av网站| 日韩熟女老妇一区二区性免费视频| 国国产精品蜜臀av免费| 十八禁高潮呻吟视频| 久久99蜜桃精品久久| 国产精品熟女久久久久浪| 久久久久国产网址| 18禁在线播放成人免费| 一区二区日韩欧美中文字幕 | 午夜激情久久久久久久| 97在线人人人人妻| 亚洲色图 男人天堂 中文字幕 | 亚洲成人一二三区av| 人妻夜夜爽99麻豆av| 国产成人a∨麻豆精品| 91精品一卡2卡3卡4卡| 搡老乐熟女国产| 在线观看www视频免费| 精品少妇黑人巨大在线播放| 欧美日韩在线观看h| 国产亚洲精品久久久com| 国产成人freesex在线| 一级a做视频免费观看| 自线自在国产av| 这个男人来自地球电影免费观看 | av播播在线观看一区| 老司机亚洲免费影院| 看十八女毛片水多多多| 亚洲一级一片aⅴ在线观看| 飞空精品影院首页| 中文欧美无线码| 女性被躁到高潮视频| 五月开心婷婷网| 久久午夜综合久久蜜桃| 九草在线视频观看| 人人澡人人妻人| 亚洲丝袜综合中文字幕| xxx大片免费视频| 乱人伦中国视频| 国产午夜精品久久久久久一区二区三区| 人妻人人澡人人爽人人| 免费高清在线观看日韩| 91精品伊人久久大香线蕉| 久久久欧美国产精品| 精品久久久久久电影网| 麻豆乱淫一区二区| 久久久亚洲精品成人影院| 午夜福利影视在线免费观看| 成人亚洲欧美一区二区av| 最近中文字幕高清免费大全6| 国产精品国产av在线观看| 最近最新中文字幕免费大全7| 18禁动态无遮挡网站| 这个男人来自地球电影免费观看 | 久久免费观看电影| a级毛片黄视频| 免费高清在线观看日韩| 亚洲高清免费不卡视频| 看十八女毛片水多多多| 久久韩国三级中文字幕| 在线精品无人区一区二区三| 青青草视频在线视频观看| 精品人妻在线不人妻| 精品久久蜜臀av无| 在线观看一区二区三区激情| 国产av国产精品国产| 欧美老熟妇乱子伦牲交| 国产亚洲最大av| 久久精品熟女亚洲av麻豆精品| 日韩制服骚丝袜av| 久久精品国产鲁丝片午夜精品| 美女主播在线视频| 人妻系列 视频| 精品人妻偷拍中文字幕| 亚洲国产精品成人久久小说| 国产成人精品无人区| 一区二区三区乱码不卡18| 少妇高潮的动态图| 午夜福利视频精品| 成人亚洲精品一区在线观看| 99久久综合免费| 免费大片黄手机在线观看| 丰满迷人的少妇在线观看| 人妻 亚洲 视频| 国产一区亚洲一区在线观看| 欧美一级a爱片免费观看看| 亚洲av.av天堂| 老熟女久久久| 日韩不卡一区二区三区视频在线| 欧美亚洲日本最大视频资源| 日本av免费视频播放| 国产乱来视频区| 美女大奶头黄色视频| 日韩免费高清中文字幕av| 日韩熟女老妇一区二区性免费视频| 天美传媒精品一区二区| 一个人看视频在线观看www免费| 伦理电影免费视频| 秋霞伦理黄片| √禁漫天堂资源中文www| 黄色欧美视频在线观看| 九色成人免费人妻av| 一级爰片在线观看| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av天美| 久久精品国产自在天天线| 国产成人aa在线观看| 久久久久久久久久人人人人人人| 日韩制服骚丝袜av| 两个人的视频大全免费| 国产精品三级大全| 九九爱精品视频在线观看| 香蕉精品网在线| 精品少妇黑人巨大在线播放| 国产成人精品婷婷| 国产视频内射| 国产av精品麻豆| 久久精品久久精品一区二区三区| 岛国毛片在线播放| 好男人视频免费观看在线| 成人亚洲精品一区在线观看| 一级毛片我不卡| 黑丝袜美女国产一区| 精品国产国语对白av| 国产成人精品在线电影| 亚洲综合精品二区| 人妻夜夜爽99麻豆av| 中文乱码字字幕精品一区二区三区| 夜夜看夜夜爽夜夜摸| 春色校园在线视频观看| 狠狠婷婷综合久久久久久88av| 99久久综合免费| 亚洲中文av在线| 国产成人免费无遮挡视频| 国产一区二区三区综合在线观看 | 午夜免费鲁丝| 尾随美女入室| 一区二区三区四区激情视频| 伊人久久精品亚洲午夜| 国精品久久久久久国模美| 黄色一级大片看看| 丝袜美足系列| 高清午夜精品一区二区三区| 成年人午夜在线观看视频| 国产免费福利视频在线观看| 老司机影院成人| 午夜福利在线观看免费完整高清在| 精品久久久精品久久久| 三级国产精品欧美在线观看| 51国产日韩欧美| 母亲3免费完整高清在线观看 | 在线精品无人区一区二区三| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 国产免费一区二区三区四区乱码| 高清午夜精品一区二区三区| 久久久久精品性色| 亚洲精品,欧美精品| 婷婷色av中文字幕| 人人妻人人澡人人看| av线在线观看网站| 成人漫画全彩无遮挡| kizo精华| 日本av手机在线免费观看| 国产欧美亚洲国产| 成人无遮挡网站| 高清午夜精品一区二区三区| 一本一本综合久久| 久久免费观看电影| 婷婷色av中文字幕| 日韩制服骚丝袜av| 日韩视频在线欧美| 久久久久人妻精品一区果冻| 久久久久国产网址| 国产在视频线精品| 欧美精品人与动牲交sv欧美| 久久久精品区二区三区| 久久精品国产自在天天线| 丝瓜视频免费看黄片| 我的女老师完整版在线观看| 男女无遮挡免费网站观看| 黑人猛操日本美女一级片| 日韩成人伦理影院| 91精品国产国语对白视频| 国产免费现黄频在线看| 日本欧美国产在线视频| 亚洲精品日韩在线中文字幕| 欧美日韩亚洲高清精品| 亚洲伊人久久精品综合| 精品视频人人做人人爽| 9色porny在线观看| 大香蕉久久网| 2021少妇久久久久久久久久久| 欧美精品人与动牲交sv欧美| 免费看av在线观看网站| 精品一区二区免费观看| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲精品一区二区精品久久久 | 国产黄色免费在线视频| 欧美日韩亚洲高清精品| 18禁裸乳无遮挡动漫免费视频| 夜夜爽夜夜爽视频| 99九九线精品视频在线观看视频| 性高湖久久久久久久久免费观看| 久久国内精品自在自线图片| 亚洲av国产av综合av卡| 特大巨黑吊av在线直播| 亚洲av.av天堂| 日韩精品免费视频一区二区三区 | 精品卡一卡二卡四卡免费| 精品熟女少妇av免费看| 亚洲精品第二区| 欧美人与性动交α欧美精品济南到 | 2018国产大陆天天弄谢| 丁香六月天网| 激情五月婷婷亚洲| 岛国毛片在线播放| 一本一本综合久久| 在线观看免费日韩欧美大片 | 大片电影免费在线观看免费| 国产伦理片在线播放av一区| 大片免费播放器 马上看| 大话2 男鬼变身卡| 欧美丝袜亚洲另类| 色婷婷av一区二区三区视频| h视频一区二区三区| 美女xxoo啪啪120秒动态图| 久久鲁丝午夜福利片| 一区二区三区免费毛片| 永久网站在线| 日本黄色片子视频| 亚洲综合色惰| 国产午夜精品一二区理论片| 亚洲色图综合在线观看| 午夜福利在线观看免费完整高清在| 在线免费观看不下载黄p国产| 美女主播在线视频| 天天影视国产精品| 久久人人爽av亚洲精品天堂| 最近中文字幕2019免费版| 欧美日韩av久久| 制服诱惑二区| 午夜福利视频在线观看免费| 国产亚洲最大av| 国产在线免费精品| 人妻人人澡人人爽人人| 汤姆久久久久久久影院中文字幕| 丝瓜视频免费看黄片| 少妇的逼好多水| 亚洲在久久综合| av国产久精品久网站免费入址| 国产女主播在线喷水免费视频网站| 久久99热6这里只有精品| 亚洲国产av新网站| 一区二区日韩欧美中文字幕 | 大香蕉97超碰在线| av播播在线观看一区| 国产乱人偷精品视频| 国产精品欧美亚洲77777| 成年人免费黄色播放视频| 高清欧美精品videossex| 日产精品乱码卡一卡2卡三| 免费黄网站久久成人精品| 久久久国产一区二区| 欧美+日韩+精品| 国产在视频线精品| 黑人欧美特级aaaaaa片| 久久精品人人爽人人爽视色| 熟女人妻精品中文字幕| 最新中文字幕久久久久| 男人添女人高潮全过程视频| 最近中文字幕高清免费大全6| av.在线天堂| 久久99热6这里只有精品| 少妇人妻久久综合中文| 中文字幕久久专区| 最近中文字幕高清免费大全6| 久久人妻熟女aⅴ| 性色av一级| 亚洲av男天堂| 国产不卡av网站在线观看| 亚洲高清免费不卡视频| 亚洲精品自拍成人| 十八禁高潮呻吟视频| 91久久精品国产一区二区三区| 性高湖久久久久久久久免费观看| 黄色视频在线播放观看不卡| 26uuu在线亚洲综合色| 大又大粗又爽又黄少妇毛片口| 色视频在线一区二区三区| 久久人妻熟女aⅴ| 另类精品久久| 777米奇影视久久| 精品少妇黑人巨大在线播放| 国产精品嫩草影院av在线观看| 成人黄色视频免费在线看| 亚洲第一av免费看| 啦啦啦中文免费视频观看日本| 三级国产精品欧美在线观看| 特大巨黑吊av在线直播| 国产精品免费大片| 人妻制服诱惑在线中文字幕| 亚洲国产精品成人久久小说| 又大又黄又爽视频免费| 亚洲精品久久成人aⅴ小说 | 飞空精品影院首页| 又粗又硬又长又爽又黄的视频| 亚洲精品美女久久av网站| 91精品伊人久久大香线蕉| 人体艺术视频欧美日本| 人妻少妇偷人精品九色| 精品国产一区二区久久| 国产不卡av网站在线观看| 高清在线视频一区二区三区| 99九九线精品视频在线观看视频| 肉色欧美久久久久久久蜜桃| 婷婷成人精品国产| 啦啦啦在线观看免费高清www| 久久久久国产网址| 亚洲欧美一区二区三区国产| 亚洲国产精品专区欧美| 国产毛片在线视频| 亚洲av免费高清在线观看| 五月伊人婷婷丁香| 日韩欧美一区视频在线观看| 热re99久久国产66热| 男女边摸边吃奶| 极品少妇高潮喷水抽搐| 丝袜在线中文字幕| 午夜久久久在线观看| av有码第一页| 少妇熟女欧美另类| 日韩大片免费观看网站| 日本黄色日本黄色录像| 国产 精品1| 亚洲怡红院男人天堂| 免费看不卡的av| 我要看黄色一级片免费的| 欧美日韩视频精品一区| 成人综合一区亚洲| 男男h啪啪无遮挡| 蜜臀久久99精品久久宅男| 一个人免费看片子| 亚洲色图综合在线观看| 亚洲精品乱码久久久久久按摩| 精品亚洲成国产av| av国产精品久久久久影院| 久久久久久久久久久久大奶| 国产日韩欧美亚洲二区| 最近中文字幕高清免费大全6| 久久久久久久久久久久大奶| 新久久久久国产一级毛片| 人人澡人人妻人| 天天操日日干夜夜撸| 国产乱人偷精品视频| 观看美女的网站| 母亲3免费完整高清在线观看 | 亚洲高清免费不卡视频| 国产黄频视频在线观看| 美女中出高潮动态图| 精品少妇内射三级| 91aial.com中文字幕在线观看| 一级爰片在线观看| 毛片一级片免费看久久久久| 国产精品成人在线| h视频一区二区三区| 性色avwww在线观看| 婷婷色综合大香蕉| kizo精华| 日韩亚洲欧美综合| 搡女人真爽免费视频火全软件| 色5月婷婷丁香| 国产精品久久久久久精品电影小说| 亚洲伊人久久精品综合| 日本黄色日本黄色录像| 免费黄色在线免费观看| 在线观看三级黄色| 久久国产精品大桥未久av|