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

    典型設施菜地中土壤微生物代謝功能多樣性

    2022-03-11 07:14:34苑學霞楊貴華王文博孫晨曦鄔元娟方麗萍
    農業(yè)工程學報 2022年24期
    關鍵詞:安丘壽光菜地

    苑學霞,張 勇,楊貴華,王文博,孫晨曦,鄔元娟,方麗萍

    ·農業(yè)生物環(huán)境與能源工程·

    典型設施菜地中土壤微生物代謝功能多樣性

    苑學霞1,張 勇2,楊貴華3,王文博1,孫晨曦1,鄔元娟1,方麗萍1

    (1. 山東省農業(yè)科學院農業(yè)質量標準與檢測技術研究所/山東省食品質量與安全檢測技術重點實驗室,濟南 250100;2. 山東省國土空間規(guī)劃院,濟南 250014;3. 濟南市農業(yè)技術推廣服務中心,濟南 250099)

    土壤微生物代謝功能多樣性是維持土壤生態(tài)系統(tǒng)健康的關鍵。為評價設施蔬菜種植對土壤微生物代謝功能多樣性的影響,該研究采用Biolog-Eco微平板法,研究2個典型設施蔬菜種植市、不同種植年限設施菜地中土壤微生物代謝功能多樣性;采用相關分析和冗余分析分別研究了與碳源利用相關的細菌群落及影響因子。結果表明:兩地設施菜地土壤中平均顏色變化率(Average Well Color Development,AWCD)、Shannon指數(shù)、Simpson指數(shù)和McIntosh指數(shù)隨著種植年限的增加而降低,但安丘種植14 a和壽光種植10 a菜地土壤不符合此規(guī)律。兩地間土壤微生物對碳源利用的差異大于種植年限導致的差異。在安丘土壤中,16種碳源(分屬糖類、氨基酸類、羧酸類、聚合物類、酚類和胺類)與不同門細菌顯著相關(<0.05);是與碳源相關種類最多的細菌,與9種碳源顯著相關(<0.05)。在壽光土壤中,11種碳源(分屬糖類、氨基酸類、羧酸類和聚合物類)與不同門細菌顯著相關(<0.05);是與碳源相關種類最多的細菌,與4種碳源顯著相關(<0.05)。在安丘土壤中,Cd對土壤微生物碳源利用有顯著負作用(<0.01),是影響土壤微生物碳源利用的最強環(huán)境因子,有機質(Organic Matter,OM)對土壤微生物碳源利用有顯著正作用(<0.01)。在壽光土壤中,Zn、OM、Cd對土壤微生物碳源利用有顯著負作用(<0.05),As、pH值對土壤微生物碳源利用有顯著正作用(<0.05)。綜上,兩地之間土壤微生物代謝功能多樣性、與之相關的細菌群落及影響因子均有差異,因此,不同地區(qū)應因地制宜,采取不同耕作措施以改善土壤微生態(tài)環(huán)境,保障設施菜地土壤健康。

    土壤;微生物;設施菜地;栽培年限;Biolog-Eco

    0 引 言

    設施蔬菜種植因其突破了氣候對蔬菜生長的限制、高產出、高收益等優(yōu)勢,成為蔬菜種植體系中極為重要的一部分。中國設施蔬菜產業(yè)發(fā)展迅猛,其栽培面積已達386萬hm2 [1]。設施蔬菜長期種植經(jīng)常引起土壤質量的變化,其主要原因包括以下兩個方面:一是由于設施菜地內封閉性、可控性、缺少雨水淋洗、復種指數(shù)高等因素,其土壤環(huán)境與露地土壤環(huán)境顯著不同[2-3];二是設施蔬菜栽培過程中化學品的大量施用和土地利用強度增加導致的土壤質量問題,具體包括:大量的化肥施用導致的土壤酸化、次生鹽漬化及營養(yǎng)元素失衡[4-5],畜禽糞便、農膜、農藥等農業(yè)投入品施用導致的污染物(重金屬、抗生素、有機氯農藥、鄰苯二甲酸酯等[6])增加。

    土壤微生物多樣性和功能的時空演變規(guī)律及其驅動機制是近年來土壤微生物學主要研究方向之一[7]。土壤微生物是土壤生態(tài)系統(tǒng)的重要組成部分,其通過調節(jié)碳、氮、磷、硫等元素水平,直接或間接地促進養(yǎng)分循環(huán)和肥力維持[8],在維護陸地生態(tài)系統(tǒng)健康方面發(fā)揮著重要作用。土壤微生物多樣性通過提高土壤微生物功能而促進作物生長[9],維持土壤的可持續(xù)利用。由于對環(huán)境因素和土壤養(yǎng)分狀況高度敏感,土壤微生物多樣性被廣泛用作土壤健康的重要評價指標[10-11]。土壤微生物對不同碳源代謝功能多樣性是土壤微生物多樣性變化的重要表征[12]。Biolog-Eco微平板因其靈敏度高、分辨力強、最大限度保留微生物群落原有的代謝特征等特點,而被廣泛應用于表征農田[12]、園地[13]、森林[14]、草地[15]、濕地[16]生態(tài)系統(tǒng)中土壤微生物代謝功能多樣性變化。不同的耕作方式和種植年限均會影響土壤微生物代謝功能多樣性。賈鵬麗等[17]發(fā)現(xiàn)玉米地中土壤微生物的Shannon指數(shù)、McIntosh指數(shù)、Simpson指數(shù)均高于水稻、土豆、大豆地。顧美英等[18]對Biolog碳代謝功能進行分析,發(fā)現(xiàn)隨著核桃定植年限的增加,土壤中細菌活性和Shannon指數(shù)先升高后降低,在10~15 a達到最高。

    關于設施菜地種植對土壤微生物影響的研究發(fā)現(xiàn),設施黃瓜連作顯著影響土壤細菌群落結構[19]。對安徽地區(qū)設施菜地進行研究,發(fā)現(xiàn)不同種植年限(3、5、7、10 a)顯著影響土壤中細菌多樣性、豐富度和群落組成[20]。設施菜地中土壤微生物群落變化主要源于栽培過程中土壤理化性質、重金屬等的變化。土壤中大多數(shù)細菌與As、Cd、Cr、Ni和Pb呈顯著正相關,而與Cu、Mn和Zn呈負相關[21]。pH值、Ni[22]、Pb[3]和空間距離[23]等均是引起土壤微生物群落變化的主導因素。以上研究多采用針對核酸的高通量測序方法,研究重點為微生物分類多樣性和遺傳多樣性。但微生物碳代謝活性是了解土壤碳素循環(huán)和儲存的核心[24],因此,有必要從土壤微生物碳代謝角度對設施菜地土壤中微生物代謝功能的變化,尤其是對不同種植年限和地點導致的差異及影響因素進行深入研究。

    山東省是中國重要的設施蔬菜生產地,其設施菜地種植面積占全國總面積的50%以上[25]。本文選取山東省2個典型設施蔬菜生產市(壽光、安丘),采用Biolog-Eco微平板法,研究不同種植年限、不同區(qū)域設施菜地中土壤微生物代謝功能多樣性的差異,進一步分析與碳源利用相關的細菌群落及利用碳源的主要影響因子。本研究將設施蔬菜栽培作用下土壤微生物對碳的代謝功能和細菌群落、土壤中非生物物質聯(lián)系起來,以期為深入了解設施菜地長期連續(xù)種植下土壤質量變化、促進設施菜地土壤的可持續(xù)利用提供科學依據(jù)。

    1 材料與方法

    1.1 研究地點及樣品采集

    研究地點位于山東省安丘市和壽光市。通過調查,在安丘市選擇5個蔬菜大棚:大棚集中在300 m范圍內,大棚結構、種植作物等基本一致,種植年限分別為3 a(AQ1)、7 a(AQ2)、10 a(AQ3)、14 a(AQ4)和16 a(AQ5)。土壤類型為潮土。以豬糞、牛糞為基肥(用量約為7.5 kg/m2),種植過程中根據(jù)生長期沖施化學肥料。種植作物為黃瓜、西紅柿,種植時間從9月下旬至來年6月中旬。在壽光市選擇6個蔬菜大棚:大棚集中在600 m范圍內,大棚結構、種植作物等基本一致,種植年限分別為3 a(SG1)、5 a(SG2)、7 a(SG3)、10 a(SG4)、15 a(SG5)和17 a(SG6)。土壤類型為褐土。以稻殼雞糞為基肥(用量約為8.3 kg/m2),種植過程中根據(jù)生長期沖施化學肥料。種植作物為茄子,種植時間從8月中旬至來年6月底。

    土壤取樣深度為0~20 cm。每個大棚設3個重復:由東向西一條直線上均勻分布3個區(qū)域,每個區(qū)域采用梅花采樣法取樣,將土壤現(xiàn)場混合均勻,去除土樣中植物根系、石塊等雜質,然后將其置于4 ℃車載冰箱中運回實驗室。按四分法分為三份:一份置于4 ℃冰箱保存,用于微生物代謝功能多樣性分析;一份風干后用瑪瑙棒研磨,過篩(孔徑0.25 mm),4 ℃避光保存,用于土壤pH值、有機質(Organic Matter,OM)和重金屬測定;最后一份置于-80 ℃超低溫冰箱保存,用于細菌群落分析。

    1.2 微生物代謝功能多樣性測定

    微生物代謝功能多樣性測定采用Biolog-Eco微平板法[12]。將10 g新鮮土壤置于100 mL 0.85%的無菌NaCl溶液中,在200 r/min下避光振蕩1 h。用0.85%的無菌NaCl溶液再稀釋100倍。將稀釋液接種到Biolog Eco微平板中,每個微孔接種150L。微平板置于(25±1) ℃下恒溫培養(yǎng)7 d,每24 h在BIOLOG EmaxTM(Biolog,Hayward,CA,USA )上讀取590 nm和750 nm 波長下的吸光值(上述操作均在無菌條件下進行)。

    土壤微生物對碳源的利用強度采用每孔平均顏色變化率(Average Well Color Development,AWCD)來描述,計算見式(1)。

    式中C為第孔的吸光值;為對照孔的吸光值;為涉及的碳源的個數(shù)。

    用Shannon、Simpson和McIntosh 3個指數(shù)表征土壤微生物對碳源利用的功能多樣性。其中Shannon指數(shù)用于評價微生物的豐富度,計算見式(2);Simpson指數(shù)用于評價土壤微生物群落中最常見物種的優(yōu)勢度,計算見式(3);McIntosh指數(shù)用于評價群落物種的均一性,計算見式(4)。

    式中P是第孔相對吸光值(C-)與相對吸光值總和的比;n是第孔相對吸光值。

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

    采用SPSS 17.0 Duncan多重比較判斷處理間顯著性差異(<0.05)。采用SPSS 17.0進行主成分分析(Principal Componets Analysis,PCA)?;赟PSS 17.0 Spearman相關性分析結果,采用Gephi 0.9.2繪制相關性網(wǎng)絡關系圖。采用Canoco 5.0對微生物對碳源的利用情況進行冗余分析(Redundancy Analysis,RDA)。

    2 結果與分析

    2.1 設施菜地中土壤微生物對碳源的利用

    AWCD與土壤中能利用單一碳源的微生物的數(shù)量和種類有關,可反映土壤微生物的密度和活性[26]。AWCD的值越大,土壤微生物的密度和活性越高,土壤微生物對碳源的總利用強度越高。AWCD隨培養(yǎng)時間變化見圖 1。由圖1可知,除AQ4外,安丘其他菜地土壤的AWCD隨著種植年限的增加而降低。除SG4外,壽光其他菜地土壤的AWCD差異如下:SG1和SG2(兩者之間差異不顯著)>SG3和SG5(兩者之間差異不顯著)>SG6。綜上,設施菜地中土壤微生物對碳源的總利用強度隨著種植年限的增加而降低。但AQ4和SG4并不遵循此規(guī)律:AQ4土壤的AWCD均顯著高于AQ2、AQ3和AQ5,SG4土壤的AWCD為壽光土壤中最低。

    Biolog-Eco微平板有31種碳源,根據(jù)化學基團分為6類:糖類、氨基酸類、羧酸類、聚合物類、酚類和胺類[24]。為分析設施菜地中土壤微生物對6類碳源的利用情況,分別計算6類碳源AWCD隨培養(yǎng)時間變化,結果見圖 2。由圖2可知,在安丘土壤中:AQ1中6類碳源的AWCD相對較高。AQ5中6類碳源的AWCD相對較低。在壽光土壤中:SG1和SG2中糖類、羧酸類、聚合物類和胺類碳源的AWCD相對較高,但SG1中氨基酸類和酚類碳源的AWCD均為較低。其他處理各類碳源AWCD沒有明顯規(guī)律。表明兩地土壤對各分類碳源的利用偏好不同,較為復雜。綜合分析培養(yǎng)168 h后6類碳源AWCD,在6類碳源中,氨基酸類碳源AWCD最高(0.791~1.211,平均值1.014),酚類碳源AWCD最低(0.251~0.900,平均值0.562)。表明土壤微生物對氨基酸類碳源的利用能力最強,對酚類碳源的利用能力最弱。

    注:AQ1~AQ5分別代表安丘種植年限3、7、10、14和16 a設施菜地;SG1~SG6分別代表壽光種植年限3、5、7、10、15和17 a設施菜地。下同。

    圖2 每類碳源利用平均顏色變化率

    2.2 設施菜地中土壤微生物功能多樣性指數(shù)

    不同的多樣性指數(shù)表征土壤微生物功能多樣性的不同方面。由表1知,除AQ4外,安丘菜地土壤中Shannon指數(shù)隨著種植年限的增加而降低,其中AQ1的Shannon指數(shù)顯著高于AQ2和AQ3,AQ2和AQ3的Shannon指數(shù)顯著高于AQ5(<0.05)。在壽光,除SG4外,SG1和SG2的Shannon指數(shù)顯著高于SG3和SG5,SG3和SG5的Shannon指數(shù)顯著高于SG6(<0.05)。Simpson指數(shù)、McIntosh指數(shù)與Shannon指數(shù)的變化規(guī)律相一致,只在統(tǒng)計顯著性上略有不同。表明設施菜地中土壤微生物的豐富度、群落中最常見物種的優(yōu)勢度、群落物種的均一性隨種植年限的增加而降低,但AQ4和SG4并不遵循此規(guī)律。綜合分析兩地之間各指數(shù)差異,發(fā)現(xiàn)安丘土壤中Shannon指數(shù)和Simpson指數(shù)均顯著高于壽光(<0.05),而安丘土壤中McIntosh指數(shù)卻顯著低于壽光(<0.05)。表明安丘土壤中微生物的豐富度、群落中最常見物種的優(yōu)勢度高于壽光,但群落物種的均一性低于壽光。

    表1 設施菜地土壤中微生物碳源利用多樣性指數(shù)

    注:同地區(qū)中同列不同字母表示差異顯著(<0.05);AQ,安丘;SG,壽光。下同。

    Note: Different letters in the same column of one county means significant difference (<0.05). AQ, Anqiu; SG, Shouguang. Same as below.

    2.3 設施菜地中土壤微生物群落代謝多樣性

    主成分分析可直觀顯示不同土壤微生物對31種碳源的利用情況。對培養(yǎng)168 h 后AWCD數(shù)據(jù)標準化變換后進行典型變量分析,提取前2個主成分因子作荷載圖,結果見圖3。由圖3知,第一主成分(PC1)方差貢獻率為28.83%,第二主成分(PC2)方差貢獻率為11.63%。兩地之間土壤微生物對碳源的利用存在明顯的空間分異。安丘土壤主要分布在第1、第2象限,壽光土壤主要分布在第3、第4象限。在同一個地方不同種植年限設施菜地土壤之間有一定的分異,但其差異小于兩地之間的差異。表明不同的土壤類型、種植模式等對土壤微生物群落代謝多樣性的影響大于種植年限的影響。

    2.4 與碳源利用相關的細菌群落

    不同設施菜地土壤中細菌群落結果見文獻[22]。相對豐度很低的微生物群落其生態(tài)功能相對有限,故本文中選擇優(yōu)勢細菌群落(相對豐度大于1%[27])(門水平,13種),分析其與碳源利用的相關性,結果見圖4。結果發(fā)現(xiàn),在安丘土壤中,6類16種碳源與不同門細菌顯著相關(<0.05),表明以上16種碳源可被土壤細菌利用;在壽光土壤中,糖類、氨基酸類、羧酸類和聚合物類等4類11種碳源與不同門細菌顯著相關(<0.05),表明以上11種碳源可被土壤細菌利用。

    圖3 設施菜地土壤中微生物碳源利用的主成分分析

    注:只顯示顯著正相關(P<0.05)。

    由圖4知,在安丘土壤中,、、和等4個門細菌與不同碳源顯著相關(<0.05)。其中與9種碳源顯著相關(<0.05),是與碳源相關種類最多的細菌,表明此地可利用的碳源種類最多。在壽光土壤中,、、等9個門細菌與不同碳源顯著相關(<0.05)。其中與4種碳源顯著相關(<0.05),是與碳源相關種類最多的細菌,表明此地可利用的碳源種類最多。

    2.5 土壤微生物利用碳源的影響因子

    不同設施菜地土壤的理化性質和重金屬含量結果見文獻[22]。對兩地土壤理化性質、重金屬和微生物對31種碳源的利用情況進行RDA分析,以評估土壤微生物利用碳源的主要影響因子,結果見圖5,其中箭頭連線的長度和夾角分別代表作用的大小和正負。由圖5可知,在安丘土壤中,RDA1和RDA2可分別解釋30.92%和16.26%的土壤微生物對碳源的利用情況。Cd對土壤微生物碳源利用的負作用顯著(<0.01),是影響土壤微生物碳源利用的最強環(huán)境因子,對總效應的貢獻率為23.4%;OM對土壤微生物碳源利用的正作用顯著(<0.01)。在壽光土壤中,RDA1和RDA2分別可解釋31.02%和9.97%的土壤微生物對碳源的利用情況。Zn、OM、Cd對土壤微生物碳源利用的負作用顯著(<0.05),As、pH值對土壤微生物碳源利用的正作用顯著(<0.05)。

    注:1.β-甲基-D-葡萄糖苷;2.D-半乳糖酸γ-內酯;3.D-木糖;4.i-赤藻糖醇;5.D-甘露醇;6.N-乙酰-D-葡萄糖胺;7.D-纖維二糖;8.α-D-葡萄糖-1-磷酸;9.α-D-乳糖;10.D, L-α-磷酸甘油;11.L-精氨酸;12.L-天冬酰胺;13.L-苯丙氨酸;14.L-絲氨酸;15.L-蘇氨酸;16.甘氨酰-L-谷氨酸;17.丙酮酸甲酯;18.D-半乳糖醛酸;19.γ-羥丁酸;20.D-葡糖胺酸;21.衣康酸;22.α-丁酮酸;23.D-蘋果酸;24.吐溫40;25.吐溫80;26.α-環(huán)式糊精;27.肝糖;28.2-羥基苯甲酸;29.4-羥基苯甲酸;30.苯乙胺;31.腐胺。

    3 討 論

    Biolog-Eco微板上的31種碳源多為植物根系分泌物的主要成分,是根際土壤微生物有效碳源的代表[2]。在東北不同土地利用類型下土壤微生物對不同碳源的利用強度不同,種植玉米、水稻、土豆、大豆和荒地土壤中微生物利用最強的碳源種類均有不同[17]。云南森林中土壤微生物對氨基酸類和羧酸類碳源的利用強度高,對酚類和胺類的利用強度低,并且受季節(jié)的影響[14]。本研究發(fā)現(xiàn)土壤微生物對氨基酸類碳源的利用能力最強,對酚類碳源的利用能力最弱。由此可知,土壤類型、植被、季節(jié)等均可影響土壤微生物對不同碳源的利用強度。

    連續(xù)多年的蔬菜種植可引發(fā)嚴重的土壤問題,進而導致土壤微生物多樣性減少[2]。高通量測序結果發(fā)現(xiàn)設施蔬菜種植對土壤微生物群落的影響:Tian等[28]發(fā)現(xiàn)隨著設施蔬菜種植年限的延長,細菌中大部分屬的相對豐度在第3年略有增加后逐漸下降;Song等[29]發(fā)現(xiàn)真菌的多樣性隨設施蔬菜種植年限的增加而降低,但細菌的多樣性沒有顯著變化。土壤中微生物群落與代謝功能密切相關[30]。設施蔬菜連作過程中,人為干擾和自毒性物質引起根際土壤微生物群落的改變,進而導致土壤微生物代謝功能的失調[31]。趙輝等[32]研究發(fā)現(xiàn),設施蔬菜種植5、7 a后,AWCD和3個多樣性指數(shù)均隨種植年限的增加而逐漸降低,與本研究結果一致。但本研究中AQ4和SG4并不符合其他年限菜地土壤中代謝功能多樣性變化的規(guī)律,具體原因有待進一步研究。

    本研究發(fā)現(xiàn)兩地之間土壤微生物對碳源利用具有差異,可能有幾個方面原因:土壤類型是土壤微生物豐度和多樣性的主要驅動因素[33],研究發(fā)現(xiàn)土壤類型對土壤微生物群落結構的影響甚至大于重金屬Cd的影響[34],本研究中兩地土壤分別屬于潮土和褐土,這可能是兩地之間土壤微生物代謝功能多樣性差異的原因之一;由于不同種類糞肥中含有的動物腸道菌、抗生素含量及營養(yǎng)成分等不同,施用牛糞、雞糞和豬糞對土壤微生物多樣性的影響也截然不同[35-36],本研究中壽光施用稻殼雞糞,而安丘施用牛糞和豬糞。因此,糞肥種類也可能是導致兩地土壤微生物代謝功能多樣性差異原因;兩地種植的作物不同對土壤微生物代謝功能多樣性不同也有貢獻。不同作物產生的植物凋落物的成分和數(shù)量[37]、根系分泌物的成分和數(shù)量[38]、植物根系的特征[39]等不同。由于植物為土壤微生物提供了大量的碳源(特別是可溶性碳源),其對微生物代謝活性及多樣性起決定性作用[40];另外,連續(xù)單作是減少土壤微生物多樣性的主要原因之一[2],本研究中壽光為茄子單作,而安丘是黃瓜和西紅柿輪作。綜上,本研究中兩地之間土壤微生物對于碳源利用的差異主要源于土壤類型、糞肥種類、種植作物及種植方式的差異。

    糖類、氨基酸類、羧酸類等根系分泌物的主要功能是為土微生物提供營養(yǎng)物質[41]。本研究結果進一步證實了這一點。安丘土壤中細菌可利用16種碳源(6類)、壽光土壤中細菌可利用11種碳源(4類)。在細菌利用碳源種類數(shù)量方面:在安丘土壤中,豐度為第二[22],是利用碳源種類最多的細菌,可以利用9種碳源。在壽光土壤中,雖然豐度較低[22],但是可利用4種碳源,是此地利用碳源種類最多的細菌。另外,兩地中雖然是豐度最高的細菌[22],但利用碳源的種類都不是最多。因此,細菌的豐度的高低與可利用碳源的種類數(shù)量多少并不一致。這個結果的另一個原因也可能是,本研究所選擇的碳源種類只有31種,所覆蓋的碳源種類有限。

    土壤微生物主要利用土壤中的可溶性有機碳進行碳代謝,因此土壤中可溶性有機碳是影響微生物碳源利用的決定性因素[40]。本研究發(fā)現(xiàn)OM是碳源利用的重要影響因子,但是在兩地中作用截然不同,在安丘土壤中OM對土壤微生物碳源利用的影響是正作用,而在壽光土壤中其影響是負作用。這可能是因為OM的成分復雜,不同土壤類型中能溶解的碳的種類和濃度可能存在差異,這一點還有待進一步研究。重金屬對土壤微生物種群具有抑制或毒性作用,其作用取決于重金屬類型和濃度[3,42],而連續(xù)多年設施蔬菜栽培可引起土壤中重金屬的積累[6,43]。在重金屬污染的土壤中,具有更高金屬耐受性的微生物比其他微生物生長得更好,導致微生物物種替換或多樣性改變。研究發(fā)現(xiàn)重金屬污染對土壤中真菌的多樣性抑制大于細菌,土壤中微生物碳儲存和對碳的利用模式隨之改變[42]。本研究中重金屬Cd(安丘、壽光)、Zn(壽光)對土壤微生物對碳源的利用強度的負作用,進一步證實了重金屬(特別是Cd)對土壤微生物代謝功能多樣性的抑制作用。另外,多項研究發(fā)現(xiàn)pH值顯著影響土壤微生物群落[20,22],本研究在壽光土壤中發(fā)現(xiàn)pH值顯著影響土壤微生物對碳源的利用。因此,pH值可能通過影響土壤中微生物群落而影響其對碳源的代謝功能。

    4 結 論

    本研究從不同種植年限和地點2個方面,初步揭示了設施菜地連續(xù)種植對土壤微生物代謝功能多樣性的影響,得到如下結論:

    1)除AQ4(安丘種植14 a)和SG4(壽光種植10 a)外,其他設施菜地土壤中平均顏色變化率、Shannon指數(shù)、Simpson指數(shù)和McIntosh指數(shù)隨著種植年限的增加而降低。

    2)兩地導致的土壤微生物對碳源利用的差異顯著,大于種植年限導致的差異。在安丘土壤中,是與碳源相關種類最多的細菌,與9種碳源顯著相關。在壽光土壤中,是與碳源相關種類最多的細菌,與4種碳源顯著相關。

    3)在安丘土壤中,Cd對土壤微生物碳源利用有顯著負作用,是影響安丘土壤微生物碳源利用的最強環(huán)境因子,有機質對土壤微生物碳源利用有顯著正作用。在壽光土壤中,Zn、OM、Cd對土壤微生物碳源利用有顯著負作用,As、pH值對土壤微生物碳源利用有顯著正作用。設施蔬菜長期種植引起的土壤質量下降及重金屬(特別是Cd)對土壤微生物的負作用應引起重視。

    [1] 成臣,易艷紅,廖雅汶,等. 贛南地區(qū)設施菜地土壤全硒含量分布特征[J]. 農業(yè)工程學報,2022,38(9):333-338.

    Cheng Chen, Yi Yanhong, Liao Yawen, et al. Distribution characteristics of the total selenium content in greenhouse vegetable soil in southern Jiangxi Province of China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(9): 333-338. (in Chinese with English abstract)

    [2] Shen W S, Hu M C, Qian D, et al. Microbial deterioration and restoration in greenhouse-based intensive vegetable production systems[J]. Plant and Soil, 2021, 463: 1-18.

    [3] Sun J, Pan L, Li Z, et al. Comparison of greenhouse and open field cultivations across China: Soil characteristics, contamination and microbial diversity[J]. Environmental Pollution, 2018, 243: 1509-1516.

    [4] Li J, Wan X, Liu X, et al. Changes in soil physical and chemical characteristics in intensively cultivated greenhouse vegetable fields in North China[J]. Soil and Tillage Research, 2019, 195: 104366.

    [5] 蔡祖聰. 我國設施栽培養(yǎng)分管理中待解的科學和技術問題[J]. 土壤學報,2019,56(1):36-43.

    Cai Zucong. Scientific and technological issures of nutrient management under greenhouse cultivation in China[J]. Acta Pedologica Sinica, 2019, 56(1): 36-43. (in Chinese with English abstract)

    [6] Kianpoor Kalkhajeh Y, Huang B, Hu W, et al. Environmental soil quality and vegetable safety under current greenhouse vegetable production management in China[J]. Agriculture, Ecosystems & Environment, 2021, 307: 107230.

    [7] 陳香,李衛(wèi)民,劉勤. 基于文獻計量的近30年國內外土壤微生物研究分析[J]. 土壤學報,2020,57(6):1458-1470.

    Chen Xiang, Li Weimin, Liu Qin. Bibliometric-based analysis of researches on soil microbes at home and abroad in the past 30 years[J]. Acta Pedologica Sinica, 2020, 57(6): 1458-1470. (in Chinese with English abstract)

    [8] Li L, Li T, Meng H, et al. Effects of seven-year fertilization reclamation on bacterial community in a coal mining subsidence area in Shanxi, China[J]. International Journal of Environmental Research and Public Health, 2021, 18(23): 12504.

    [9] Shu X, He J, Zhou Z, et al. Organic amendments enhance soil microbial diversity, microbial functionality and crop yields: A meta-analysis[J]. Science of the Total Environment, 2022, 829: 154627.

    [10] 王宇峰,孟會生,李廷亮,等. 培肥措施對復墾土壤微生物碳氮代謝功能多樣性的影響[J]. 農業(yè)工程學報,2020,36(24):81-90.

    Wang Yufeng, Meng Huisheng, Li Tingliang, et al. Effects of fertilization regime on the functional diversity of microbial carbon and nitrogen metabolism in reclaimed soil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(24): 81-90. (in Chinese with English abstract)

    [11] Wang C N, Wu R L, Li Y Y, et al. Effects of pesticide residues on bacterial community diversity and structure in typical greenhouse soils with increasing cultivation years in Northern China[J]. Science of the Total Environment, 2020, 710: 136321.

    [12] He W, Ye W, Sun M, et al. Soil phosphorus availability and stoichiometry determine microbial activity and functional diversity of fluvo-aquic soils under long-term fertilization regimes[J]. Journal of Soils and Sediments, 2022, 22(4): 1214-1227.

    [13] Zhu Z, Bai Y, Lv M, et al. Soil fertility, microbial biomass, and microbial functional diversity responses to four years fertilization in an apple orchard in North China[J]. Horticultural Plant Journal, 2020, 6(4): 223-230.

    [14] 何中聲,陳佳嘉,朱靜,等. 戴云山南坡不同海拔森林土壤微生物功能多樣性特征及影響因素[J]. 生態(tài)學報,2022,42(9):3504-3515.

    He Zhongsheng, Chen Jiajia, Zhu Jing, et al. Characteristics of microbial functional diversity and its influecing factors of forest soils at different elevations on the southern slope of Daiyun Mountain[J]. Acta Ecologica Sinica, 2022, 42(9): 3504-3515. (in Chinese with English abstract)

    [15] Furtak K, Grz?dziel J, Ga??zka A, et al. Prevalence of unclassified bacteria in the soil bacterial community from floodplain meadows (fluvisols) under simulated flood conditions revealed by a metataxonomic approachss[J]. Catena, 2020, 188: 104448.

    [16] Yin Y, Yan Z. Variations of soil bacterial diversity and metabolic function with tidal flat elevation gradient in an artificial mangrove wetland[J]. Science of The Total Environment, 2020, 718: 137385.

    [17] 賈鵬麗,馮海艷,李淼. 東北黑土區(qū)不同土地利用方式下農田土壤微生物多樣性[J]. 農業(yè)工程學報,2020,36(20):171-178.

    Jia Pengli, Feng Haiyan, Li Miao. Soil microbial diverstiy of black soil under different land use patterns in northeast China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(20): 171-178. (in Chinese with English abstract)

    [18] 顧美英,徐萬里,馬凱,等. 不同定植年限核桃園土壤細菌群落多樣性及碳代謝功能特征[J]. 生態(tài)學雜志,2021,40(7):2045-2056.

    Gu Meiying, Xu Wanli, Ma Kai, et al. Soil bacterial communtiy diversity and carbon source emtabolism function in walnut orchard with different stand ages[J]. Chinese Journal of Ecology, 2021, 40(7): 2045-2056. (in Chinese with English abstract)

    [19] Zhao Y N, Mao X X, Zhang M S, et al. Response of soil microbial communities to continuously mono-cropped cucumber under greenhouse conditions in a calcareous soil of north China[J]. Journal of Soils and Sediments, 2020, 20(5): 2446-2459.

    [20] Zhang J, Wang P, Tian H, et al. Pyrosequencing-based assessment of soil microbial community structure and analysis of soil properties with vegetable planted at different years under greenhouse conditions[J]. Soil and Tillage Research, 2019, 187: 1-10.

    [21] Deng W, Zhang A, Chen S, et al. Heavy metals, antibiotics and nutrients affect the bacterial community and resistance genes in chicken manure composting and fertilized soil[J]. Journal of Environmental Management, 2020, 257: 109980.

    [22] Yuan X, Zhang Y, Sun C, et al. Profile of bacterial community and antibiotic resistance genes in typical vegetable greenhouse soil[J]. International Journal of Environmental Research and Public Health, 2022, 19(13): 7742.

    [23] Zhang K, Adams J M, Shi Y, et al. Environment and geographic distance differ in relative importance for determining fungal community of rhizosphere and bulk soil[J]. Environmental Microbiology, 2017, 19(9): 3649-3659.

    [24] Liu Z D, Song Y Y, Ma X Y, et al. Deep soil microbial carbon metabolic function is important but often neglected: A study on the Songnen Plain reed wetland, Northeast China[J/OL]. Fundamental Research, 2022. [2022-12-20]. https://doi.org/10.1016/j.fmre.2022.06.013.

    [25] Zhao X, Wang J, Zhu L, et al. Field-based evidence for enrichment of antibiotic resistance genes and mobile genetic elements in manure-amended vegetable soils[J]. Science of the Total Environment, 2019, 654: 906-913.

    [26] Choi K H, Dobbs F C. Comparison of two kinds of Biolog microplates (GN and ECO) in their ability to distinguish among aquatic microbial communities[J]. Journal of Microbiological Methods, 1999, 36(3): 203-213.

    [27] Yang Y, Chen Q, Yu W, et al. Estimating soil bacterial abundance and diversity in the Southeast Qinghai-Tibet Plateau[J]. Geoderma, 2022, 416: 115807.

    [28] Tian T, Chen Z, Tian Y Q, et al. Microbial diversity in solar greenhouse soils in Round-Bohai Bay-Region, China: The influence of cultivation year and environmental condition[J]. Environmental Science and Pollution Research, 2017, 24(29): 23236-23249.

    [29] Song Y, Xu M, Li X N, et al. Long-term plastic greenhouse cultivation changes soil microbial community structures: A case study[J]. Journal of Agricultural and Food Chemistry, 2018, 66(34): 8941-8948.

    [30] Song Y, Li X N, Yao S, et al. Correlations between soil metabolomics and bacterial community structures in the pepper rhizosphere under plastic greenhouse cultivation[J]. Science of the Total Environment, 2020, 728: 138439.

    [31] Li Y, Chi J, Ao J, et al. Effects of different continuous cropping years on bacterial community and diversity of cucumber rhizosphere soil in solar-greenhouse[J]. Current Microbiology, 2021, 78(6): 2380-2390.

    [32] 趙輝,王喜英,徐仕強,等. 貴州武陵片區(qū)不同種植年限設施菜地土壤微生物群落的結構和功能多樣性[J]. 河南農業(yè)科學,2021,50(1):81-91.

    Zhao Hui, Wang Xiying, Xu Shiqiang, et al. Soil microbial communtiy structure and functional diversity in vegetable greenhouse for differetn plantin years in Wuling area of Guizhou province[J]. Journal of Henan Agricultural Sciences, 2021, 50(1): 81-91. (in Chinese with English abstract)

    [33] Sadet-Bourgeteau S, Houot S, Karimi B, et al. Microbial communities from different soil types respond differently to organic waste input[J]. Applied Soil Ecology, 2019, 143: 70-79.

    [34] Liu H, Yang Y, Yang Y, et al. Dynamics of fungal and bacterial communities in different types of soil ageing with different dosages of cadmium[J]. Ecotoxicology and Environmental Safety, 2022, 242: 113860.

    [35] Han X M, Hu H W, Chen Q L, et al. Antibiotic resistance genes and associated bacterial communities in agricultural soils amended with different sources of animal manures[J]. Soil Biology & Biochemistry, 2018, 126: 91-102.

    [36] Yuan X, Zhang Y, Fan L, et al. Temporal dynamics of antibiotic resistance genes in vegetable greenhouse soils following different manure applications[J]. Journal of Soil Science and Plant Nutrition, 2022, 22: 5144-5158.

    [37] Xu Z, Yu G, Zhang X, et al. Divergence of dominant factors in soil microbial communities and functions in forest ecosystems along a climatic gradient[J]. Biogeosciences, 2018, 15: 1217-1228.

    [38] Williams A, Langridge H, Straathof A L, et al. Root functional traits explain root exudation rate and composition across a range of grassland species[J]. Journal of Ecology, 2022, 110(1): 21-33.

    [39] Spitzer C M, Lindahl B, Wardle D A, et al. Root trait–microbial relationships across tundra plant species[J]. New Phytologist, 2021, 229(3): 1508-1520.

    [40] Hernández-Cáceres D, Stokes A, Angeles-Alvarez G, et al. Vegetation creates microenvironments that influence soil microbial activity and functional diversity along an elevation gradient[J]. Soil Biology and Biochemistry, 2022, 165: 108485.

    [41] 吳林坤,林向民,林文雄. 根系分泌物介導下植物-土壤-微生物互作關系研究進展與展望[J]. 植物生態(tài)學報,2014,38(3):298-310.

    Wu Linkun, Lin Xiangmin, Lin Wenxiong. Advances and perspective in research on plant-soil-microbe interactions mediated by root exudates[J]. Chinese Journal of Plant Ecology, 2014, 38(3): 298-310. (in Chinese with English abstract)

    [42] Xu Y, Seshadri B, Bolan N, et al. Microbial functional diversity and carbon use feedback in soils as affected by heavy metals[J]. Environment International, 2019, 125: 478-488.

    [43] 苑學霞,方麗萍,張?zhí)?,? 不同年限設施菜地土壤中重金屬和抗生素污染特征[J]. 生態(tài)環(huán)境學報,2020,29(8):1669-1674.

    Yuan Xuexia, Fang Liping, Zhang Taiping, et al. Characteristics of heavy metals and antibiotics pollution in vegetabel greenhouses after different culivation years[J]. Ecology and Environmental Sciences, 2020, 29(8): 1669-1674. (in Chinese with English abstract)

    Soil microbial metabolism functional diversity in typical facility vegetable fields

    Yuan Xuexia1, Zhang Yong2, Yang Guihua3, Wang Wenbo1, Sun Chenxi1, Wu Yuanjuan1, Fang Liping1

    (1.-,,250100,;2.,250014,;3.,250099,)

    Metabolic functional diversity of soil microorganisms can greatly contribute to the health of the soil ecosystems in facility vegetable field. However, the soil quality is seriously deteriorated in recent years, due to the frequency input of organic and chemical fertilizers, pesticides, and fungicides. The Biolog-Eco method has been widely used to demonstrate the carbon metabolic activity of culturable microorganisms. The ecological status of microbial communities can also be characterized by the metabolic functional diversity in diverse environments. In this study, a systematic investigation was carried out to clarify the impacts of cultivation periods and geographical location on the functional diversity of soil microbial metabolism using a Biolog-Eco microplate. Correlation analysis and Redundancy Analysis (RDA) were used to reveal the relative bacterial phyla and effect factors. The soil samples were collected from 11 facility vegetable fields in the two typical vegetable cultivation counties (Anqiu and Shouguang) of China. The key influencing factors were then identified for the microorganisms and carbon source utilization. The results showed that the Average Well Color Development (AWCD), Simpson index, McIntosh index, and Shannon index decreased with the increase of cultivation years, except for the AQ4 (cultivation history of 14 a in Anqiu), and SG4 (cultivation history of 10 a in Shouguang). The Principal Component Analysis (PCA) result demonstrated that the difference in the soil microbial utilization of carbon sources between two counties was much more significant than those among the soils with different cultivation histories in the same county. It infers that the geographical location played a more important role than the cultivation years. Correlation analysis showed that 16 different carbon sources (six types of carbon source involved) were closely related to the different bacterial phyla (<0.05), andof soil bacterial phylum was positively related to the most kinds of carbon sources (9 kinds) in the Anqiu soil. In Shouguang soil, 11 different carbon sources (four types: carbohydrates, amino acids, carboxylic acids, and polymers) were closely related to the different bacterial phyla (<0.05), andpositively related to most kinds of carbon sources (4 kinds). RDA showed that the Cd was the top environmental factor, negatively affecting the soil microbial carbon sources utilization in Anqiu soil (<0.01), while the positive effect of organic matter was significant (<0.01). Zn, organic matter, and Cd were negatively dominated the utilization of microbial carbon sources (<0.05) in Shouguang, whereas, the effects of As and pH were significantly positive (<0.05). This finding can provide direct evidence that the significant decrease in soil microbial functional diversity was caused by the long-term cultivation of facility vegetable, different relative bacteria and effect factors caused by location in the north of China. It is urgent to take different measures according to different situation to improve the soil micro-environment for the healthy environments in the vegetable fields.

    soils; microorganism; facility vegetable field; cultivation years; Biolog-Eco

    10.11975/j.issn.1002-6819.2022.24.016

    X172

    A

    1002-6819(2022)-24-0145-08

    苑學霞,張勇,楊貴華,等. 典型設施菜地中土壤微生物代謝功能多樣性[J]. 農業(yè)工程學報,2022,38(24):145-152.doi:10.11975/j.issn.1002-6819.2022.24.016 http://www.tcsae.org

    Yuan Xuexia, Zhang Yong, Yang Guihua, et al. Soil microbial metabolism functional diversity in typical facility vegetablefields[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(24): 145-152. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.24.016 http://www.tcsae.org

    2022-10-10

    2022-11-16

    國家自然科學基金項目(41907030)

    苑學霞,研究員,研究方向為農業(yè)環(huán)境安全與風險評估。Email:xuexiayuan@163.com

    猜你喜歡
    安丘壽光菜地
    小白兔進菜地
    菜地游樂園
    菜地的周長
    菜地里的狂歡
    山東安丘倉方早生桃生產管理技術探討
    山東壽光農產品批發(fā)市場最新蔬菜行情
    壽光大蔥
    安丘蜜桃標準化栽培技術規(guī)程
    中國果菜(2015年2期)2015-03-11 20:01:07
    安丘董家莊漢畫像石墓主人之謎
    大眾考古(2014年2期)2014-06-26 08:29:32
    壽光溫室野蛞蝓的發(fā)生與防治
    長江蔬菜(2014年1期)2014-03-11 15:09:57
    狂野欧美激情性bbbbbb| 国产精品亚洲av一区麻豆 | 亚洲一级一片aⅴ在线观看| 亚洲熟女毛片儿| 亚洲av日韩在线播放| 亚洲欧美一区二区三区久久| 最近中文字幕2019免费版| 午夜免费鲁丝| 精品人妻熟女毛片av久久网站| 多毛熟女@视频| 久久久精品免费免费高清| 久久久久久久精品精品| 国产一卡二卡三卡精品 | 巨乳人妻的诱惑在线观看| 国产男女内射视频| 国产成人精品福利久久| 亚洲情色 制服丝袜| 精品酒店卫生间| 在线观看免费高清a一片| 亚洲av日韩在线播放| 亚洲精品第二区| 如日韩欧美国产精品一区二区三区| 国产日韩欧美视频二区| 国产极品天堂在线| 麻豆精品久久久久久蜜桃| 老鸭窝网址在线观看| 国产精品久久久久久精品电影小说| 欧美在线一区亚洲| 超碰成人久久| 成人毛片60女人毛片免费| 国产精品一区二区精品视频观看| 纵有疾风起免费观看全集完整版| 97人妻天天添夜夜摸| 女人久久www免费人成看片| 高清不卡的av网站| 免费看不卡的av| 国产免费视频播放在线视频| 黑丝袜美女国产一区| 国产精品一区二区在线不卡| 国产淫语在线视频| 亚洲国产精品一区三区| 亚洲中文av在线| av网站在线播放免费| 国产成人一区二区在线| 国产黄色视频一区二区在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲成色77777| av在线老鸭窝| 久久久精品免费免费高清| 国产淫语在线视频| 男人添女人高潮全过程视频| 秋霞在线观看毛片| 免费在线观看视频国产中文字幕亚洲 | 国产精品av久久久久免费| 国产一卡二卡三卡精品 | 啦啦啦 在线观看视频| 在线亚洲精品国产二区图片欧美| 咕卡用的链子| 久久影院123| 亚洲av电影在线进入| 青春草视频在线免费观看| 久久久久久人妻| 一二三四在线观看免费中文在| 欧美日韩av久久| 成人18禁高潮啪啪吃奶动态图| 看非洲黑人一级黄片| 交换朋友夫妻互换小说| 日韩,欧美,国产一区二区三区| 男女无遮挡免费网站观看| 丝瓜视频免费看黄片| 欧美成人午夜精品| 亚洲四区av| 亚洲国产欧美网| 国产成人免费无遮挡视频| 青青草视频在线视频观看| 丰满乱子伦码专区| 丰满乱子伦码专区| 国产亚洲午夜精品一区二区久久| 亚洲国产欧美在线一区| 久久久精品国产亚洲av高清涩受| 汤姆久久久久久久影院中文字幕| 国产免费现黄频在线看| 看十八女毛片水多多多| 色网站视频免费| 国产有黄有色有爽视频| 国产福利在线免费观看视频| 成人亚洲欧美一区二区av| 日韩欧美一区视频在线观看| 亚洲av综合色区一区| 国产日韩欧美视频二区| 免费看av在线观看网站| 婷婷成人精品国产| 丝袜脚勾引网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品国产av成人精品| 日日摸夜夜添夜夜爱| 91精品国产国语对白视频| 一区二区三区四区激情视频| 国产精品一区二区在线不卡| www.精华液| 亚洲精品国产色婷婷电影| 欧美激情高清一区二区三区 | 男女免费视频国产| 亚洲人成77777在线视频| 男女无遮挡免费网站观看| 97精品久久久久久久久久精品| 日韩成人av中文字幕在线观看| 国产激情久久老熟女| 欧美日韩成人在线一区二区| videosex国产| 精品酒店卫生间| 久久久国产欧美日韩av| 欧美日韩精品网址| av在线播放精品| 国产精品 欧美亚洲| 免费人妻精品一区二区三区视频| 久久久久精品人妻al黑| 精品国产一区二区三区四区第35| 黄色 视频免费看| 午夜日韩欧美国产| √禁漫天堂资源中文www| 日本wwww免费看| 亚洲一区二区三区欧美精品| 久久久久久免费高清国产稀缺| 制服人妻中文乱码| 少妇人妻 视频| 亚洲激情五月婷婷啪啪| 啦啦啦 在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 在线免费观看不下载黄p国产| 久久精品国产综合久久久| 中文精品一卡2卡3卡4更新| av在线观看视频网站免费| 在线天堂最新版资源| 精品酒店卫生间| 精品免费久久久久久久清纯 | 中文欧美无线码| 黄色一级大片看看| 久久久久久人妻| 久热爱精品视频在线9| 又黄又粗又硬又大视频| 制服人妻中文乱码| 操美女的视频在线观看| 少妇人妻久久综合中文| 亚洲国产欧美一区二区综合| 国产在线视频一区二区| 巨乳人妻的诱惑在线观看| 亚洲美女搞黄在线观看| 国产精品嫩草影院av在线观看| 黑人欧美特级aaaaaa片| 精品一区在线观看国产| 日本wwww免费看| 欧美亚洲 丝袜 人妻 在线| 成人毛片60女人毛片免费| 久久女婷五月综合色啪小说| 午夜福利乱码中文字幕| 亚洲国产精品国产精品| 黄色 视频免费看| 女性被躁到高潮视频| 天美传媒精品一区二区| 青春草亚洲视频在线观看| 亚洲成国产人片在线观看| 国产一区二区三区av在线| 久久人妻熟女aⅴ| 成人亚洲精品一区在线观看| 久久婷婷青草| 日本91视频免费播放| 亚洲伊人色综图| 韩国av在线不卡| 亚洲综合精品二区| 激情视频va一区二区三区| 纯流量卡能插随身wifi吗| 免费少妇av软件| 宅男免费午夜| 国产精品一区二区在线观看99| 中文字幕人妻丝袜制服| 免费高清在线观看视频在线观看| 免费在线观看黄色视频的| 久久久国产一区二区| 亚洲精品久久午夜乱码| 不卡av一区二区三区| 婷婷成人精品国产| 久久婷婷青草| 欧美激情极品国产一区二区三区| 90打野战视频偷拍视频| 伊人亚洲综合成人网| av在线观看视频网站免费| 日本av免费视频播放| 亚洲国产av新网站| 日韩av不卡免费在线播放| 又大又爽又粗| av免费观看日本| 成人国产av品久久久| 精品一品国产午夜福利视频| 欧美97在线视频| h视频一区二区三区| 国产又色又爽无遮挡免| 国产成人一区二区在线| 美女高潮到喷水免费观看| 国产精品av久久久久免费| 高清黄色对白视频在线免费看| 午夜日本视频在线| 午夜av观看不卡| 国产一区二区在线观看av| 成人国语在线视频| 搡老乐熟女国产| 婷婷色综合www| 999精品在线视频| tube8黄色片| 久久天堂一区二区三区四区| 女人久久www免费人成看片| 国产精品久久久久久人妻精品电影 | 亚洲国产毛片av蜜桃av| 岛国毛片在线播放| 免费久久久久久久精品成人欧美视频| 成人国语在线视频| 精品一区在线观看国产| 大片免费播放器 马上看| 麻豆精品久久久久久蜜桃| 十八禁高潮呻吟视频| 亚洲精品一二三| 成人亚洲精品一区在线观看| 久久99一区二区三区| 国产国语露脸激情在线看| 日本猛色少妇xxxxx猛交久久| 丝袜美足系列| 国产日韩欧美视频二区| 国产又爽黄色视频| 欧美激情极品国产一区二区三区| 十八禁人妻一区二区| 久久久久视频综合| xxxhd国产人妻xxx| www日本在线高清视频| 国精品久久久久久国模美| 国产精品熟女久久久久浪| 女人爽到高潮嗷嗷叫在线视频| 亚洲综合精品二区| 婷婷色综合www| 亚洲成人手机| 欧美久久黑人一区二区| 菩萨蛮人人尽说江南好唐韦庄| 欧美最新免费一区二区三区| 肉色欧美久久久久久久蜜桃| 美女脱内裤让男人舔精品视频| 国产精品av久久久久免费| 国产精品久久久久久人妻精品电影 | 亚洲美女搞黄在线观看| 国产爽快片一区二区三区| 国产淫语在线视频| 一区福利在线观看| 日韩精品免费视频一区二区三区| 欧美人与善性xxx| 99re6热这里在线精品视频| 岛国毛片在线播放| av有码第一页| 欧美最新免费一区二区三区| 天天躁夜夜躁狠狠久久av| 午夜福利网站1000一区二区三区| 亚洲国产欧美一区二区综合| 91aial.com中文字幕在线观看| 十八禁人妻一区二区| 欧美人与性动交α欧美精品济南到| 日本色播在线视频| 日韩伦理黄色片| 狂野欧美激情性xxxx| 一级片'在线观看视频| 国产高清国产精品国产三级| 久久久久久久精品精品| 国精品久久久久久国模美| 国产精品一区二区精品视频观看| 久久久久国产精品人妻一区二区| 一边摸一边做爽爽视频免费| 国产爽快片一区二区三区| 亚洲国产精品国产精品| 国产成人午夜福利电影在线观看| 久久亚洲国产成人精品v| 捣出白浆h1v1| 欧美xxⅹ黑人| 久久久亚洲精品成人影院| 最黄视频免费看| 免费不卡黄色视频| 少妇的丰满在线观看| 青春草亚洲视频在线观看| 日本av手机在线免费观看| 国产精品久久久久成人av| 国产精品熟女久久久久浪| 97精品久久久久久久久久精品| 国产欧美日韩一区二区三区在线| 国产在线一区二区三区精| a级片在线免费高清观看视频| 人妻一区二区av| 国产高清不卡午夜福利| 国产福利在线免费观看视频| 视频在线观看一区二区三区| 色综合欧美亚洲国产小说| 欧美另类一区| 日日爽夜夜爽网站| 最近手机中文字幕大全| 亚洲国产精品成人久久小说| 在线亚洲精品国产二区图片欧美| 人人妻,人人澡人人爽秒播 | 国语对白做爰xxxⅹ性视频网站| 免费日韩欧美在线观看| 欧美亚洲日本最大视频资源| 精品亚洲乱码少妇综合久久| 男人添女人高潮全过程视频| 岛国毛片在线播放| 免费女性裸体啪啪无遮挡网站| 久久精品久久久久久噜噜老黄| 各种免费的搞黄视频| 黄片播放在线免费| 好男人视频免费观看在线| 五月天丁香电影| 人人妻,人人澡人人爽秒播 | 日韩一区二区三区影片| 青春草视频在线免费观看| 男女午夜视频在线观看| 精品少妇黑人巨大在线播放| 91aial.com中文字幕在线观看| 天堂8中文在线网| 丰满饥渴人妻一区二区三| 国产不卡av网站在线观看| 最近中文字幕2019免费版| 婷婷色综合大香蕉| av在线播放精品| 欧美日韩亚洲高清精品| 99re6热这里在线精品视频| 最黄视频免费看| 日本爱情动作片www.在线观看| 麻豆乱淫一区二区| 亚洲人成网站在线观看播放| 亚洲精品视频女| 国产精品 欧美亚洲| 人妻人人澡人人爽人人| 国产深夜福利视频在线观看| 久久精品人人爽人人爽视色| 别揉我奶头~嗯~啊~动态视频 | 精品一区二区三卡| 国产免费福利视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美一区二区三区久久| 老司机靠b影院| 啦啦啦 在线观看视频| av免费观看日本| 久久精品久久精品一区二区三区| 欧美精品亚洲一区二区| 国产成人欧美| 亚洲视频免费观看视频| 欧美精品人与动牲交sv欧美| 国产男人的电影天堂91| 欧美人与性动交α欧美软件| 韩国精品一区二区三区| 岛国毛片在线播放| av福利片在线| 久久97久久精品| 国产成人欧美| 欧美日本中文国产一区发布| 波多野结衣一区麻豆| 久久久久久久精品精品| 欧美亚洲日本最大视频资源| 黄色一级大片看看| 最黄视频免费看| 美女福利国产在线| 亚洲美女搞黄在线观看| 欧美最新免费一区二区三区| 一级毛片 在线播放| 亚洲欧美清纯卡通| 国产精品香港三级国产av潘金莲 | 最新的欧美精品一区二区| 欧美久久黑人一区二区| 国产 精品1| 欧美老熟妇乱子伦牲交| av在线观看视频网站免费| 我的亚洲天堂| 麻豆乱淫一区二区| 黄色视频不卡| 免费不卡黄色视频| 久久久久网色| 欧美老熟妇乱子伦牲交| 制服诱惑二区| 天美传媒精品一区二区| 亚洲精品国产av蜜桃| 人人妻人人澡人人爽人人夜夜| 少妇人妻久久综合中文| 国产老妇伦熟女老妇高清| 国产成人午夜福利电影在线观看| 免费人妻精品一区二区三区视频| 啦啦啦中文免费视频观看日本| 日日啪夜夜爽| av网站免费在线观看视频| 久久人人爽av亚洲精品天堂| 亚洲精品国产av蜜桃| 国产精品av久久久久免费| 高清av免费在线| 另类精品久久| 久久精品人人爽人人爽视色| 中文字幕色久视频| 麻豆乱淫一区二区| 国产一区二区激情短视频 | 999久久久国产精品视频| 蜜桃在线观看..| 一级片免费观看大全| 国产亚洲午夜精品一区二区久久| www日本在线高清视频| 另类亚洲欧美激情| 国产女主播在线喷水免费视频网站| 成人影院久久| 最近2019中文字幕mv第一页| 亚洲,欧美精品.| 亚洲av国产av综合av卡| 丰满乱子伦码专区| 午夜激情av网站| 日韩精品免费视频一区二区三区| 国产福利在线免费观看视频| 精品午夜福利在线看| 国产精品熟女久久久久浪| 亚洲成人免费av在线播放| 精品国产露脸久久av麻豆| 一区二区av电影网| 丰满少妇做爰视频| 97在线人人人人妻| 亚洲第一区二区三区不卡| 极品人妻少妇av视频| 午夜日韩欧美国产| 亚洲av成人不卡在线观看播放网 | 亚洲国产精品一区三区| 亚洲精品第二区| 国产在线免费精品| 观看av在线不卡| 国语对白做爰xxxⅹ性视频网站| 高清在线视频一区二区三区| 国产精品av久久久久免费| 七月丁香在线播放| 日韩免费高清中文字幕av| 国产精品麻豆人妻色哟哟久久| 国产av码专区亚洲av| 各种免费的搞黄视频| 五月天丁香电影| 80岁老熟妇乱子伦牲交| 亚洲精品自拍成人| 我的亚洲天堂| 中文字幕制服av| 亚洲欧美一区二区三区黑人| 一本大道久久a久久精品| 老司机靠b影院| 国产一区二区激情短视频 | 韩国高清视频一区二区三区| 一本一本久久a久久精品综合妖精| 最新在线观看一区二区三区 | 免费在线观看黄色视频的| 国产成人a∨麻豆精品| 精品亚洲成a人片在线观看| 免费观看av网站的网址| 又粗又硬又长又爽又黄的视频| 丝袜脚勾引网站| videosex国产| 亚洲免费av在线视频| 亚洲一区中文字幕在线| 秋霞伦理黄片| 在线观看三级黄色| 黑人巨大精品欧美一区二区蜜桃| 亚洲av国产av综合av卡| 男人爽女人下面视频在线观看| av福利片在线| 日日爽夜夜爽网站| 免费黄频网站在线观看国产| 欧美xxⅹ黑人| 黄片播放在线免费| 久久人人97超碰香蕉20202| 你懂的网址亚洲精品在线观看| 少妇的丰满在线观看| 午夜激情av网站| 欧美老熟妇乱子伦牲交| 中文字幕人妻熟女乱码| 日韩免费高清中文字幕av| 精品一品国产午夜福利视频| 成人国产麻豆网| 丝袜美足系列| 日韩精品免费视频一区二区三区| 亚洲av中文av极速乱| av国产精品久久久久影院| 久久久久久人妻| 久久精品国产亚洲av高清一级| 综合色丁香网| 亚洲精品一区蜜桃| 麻豆av在线久日| 日本爱情动作片www.在线观看| 人人妻,人人澡人人爽秒播 | 999精品在线视频| 亚洲成人av在线免费| 丰满乱子伦码专区| 欧美97在线视频| 久久这里只有精品19| 丝袜脚勾引网站| 精品国产乱码久久久久久男人| 亚洲欧洲精品一区二区精品久久久 | 老司机影院成人| 午夜福利免费观看在线| 久久人人爽av亚洲精品天堂| 精品国产一区二区三区四区第35| 少妇人妻 视频| 热re99久久精品国产66热6| 最近最新中文字幕大全免费视频 | 一二三四中文在线观看免费高清| 精品人妻在线不人妻| 亚洲av男天堂| 如何舔出高潮| 91国产中文字幕| 侵犯人妻中文字幕一二三四区| 国产熟女欧美一区二区| 亚洲精品av麻豆狂野| 久久毛片免费看一区二区三区| 女人精品久久久久毛片| 777久久人妻少妇嫩草av网站| 亚洲成色77777| 伦理电影免费视频| 美女午夜性视频免费| 七月丁香在线播放| 久久久久久人妻| 成人黄色视频免费在线看| 一级黄片播放器| 日韩一本色道免费dvd| 老司机亚洲免费影院| 日韩熟女老妇一区二区性免费视频| 日本猛色少妇xxxxx猛交久久| 国产爽快片一区二区三区| 在线观看免费日韩欧美大片| 美女国产高潮福利片在线看| 亚洲成色77777| 人人妻人人添人人爽欧美一区卜| avwww免费| 成人国语在线视频| 久久久久久人妻| 亚洲色图 男人天堂 中文字幕| 亚洲欧美精品综合一区二区三区| 天堂中文最新版在线下载| 色播在线永久视频| 老汉色av国产亚洲站长工具| 欧美精品av麻豆av| 久久婷婷青草| 国产精品人妻久久久影院| 高清在线视频一区二区三区| 国产国语露脸激情在线看| 国产精品久久久久久精品古装| av在线观看视频网站免费| 精品少妇一区二区三区视频日本电影 | 成年动漫av网址| 久久这里只有精品19| 国产无遮挡羞羞视频在线观看| 亚洲一区二区三区欧美精品| 中文乱码字字幕精品一区二区三区| 久久鲁丝午夜福利片| 精品一区在线观看国产| 巨乳人妻的诱惑在线观看| bbb黄色大片| 三上悠亚av全集在线观看| 日本色播在线视频| 老汉色∧v一级毛片| 香蕉国产在线看| 亚洲欧美精品综合一区二区三区| 纵有疾风起免费观看全集完整版| 亚洲自偷自拍图片 自拍| 美女脱内裤让男人舔精品视频| 热re99久久国产66热| 波野结衣二区三区在线| netflix在线观看网站| 老汉色av国产亚洲站长工具| 我的亚洲天堂| 亚洲,欧美,日韩| 精品第一国产精品| 国产精品国产av在线观看| 久久久久精品国产欧美久久久 | 在线精品无人区一区二区三| 制服人妻中文乱码| 多毛熟女@视频| 亚洲国产看品久久| 中国三级夫妇交换| 免费黄频网站在线观看国产| 性少妇av在线| www.自偷自拍.com| 午夜老司机福利片| 自线自在国产av| 国产国语露脸激情在线看| 九九爱精品视频在线观看| 在线观看免费午夜福利视频| 999久久久国产精品视频| 日韩伦理黄色片| 国产爽快片一区二区三区| 一区二区av电影网| 伊人久久国产一区二区| 日本一区二区免费在线视频| 丝瓜视频免费看黄片| 又大又黄又爽视频免费| 久久精品久久久久久久性| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久久久久免| 欧美日韩成人在线一区二区| 少妇 在线观看| 婷婷色麻豆天堂久久| 久久久久国产精品人妻一区二区| 精品福利永久在线观看| 一级a爱视频在线免费观看| 欧美精品一区二区免费开放| 三上悠亚av全集在线观看| 精品一区二区免费观看| 精品一区二区三卡| 午夜久久久在线观看| 最近2019中文字幕mv第一页| 在线看a的网站| 青春草国产在线视频| 性少妇av在线| 另类亚洲欧美激情| 搡老乐熟女国产| 日本黄色日本黄色录像| xxxhd国产人妻xxx|