• <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
    婷婷色综合www| 水蜜桃什么品种好| av线在线观看网站| 亚洲精品美女久久av网站| 中文字幕另类日韩欧美亚洲嫩草| 国产麻豆69| 黑人欧美特级aaaaaa片| 黄网站色视频无遮挡免费观看| 一级爰片在线观看| 大话2 男鬼变身卡| 啦啦啦在线观看免费高清www| 99久久99久久久精品蜜桃| 妹子高潮喷水视频| 男女床上黄色一级片免费看| 日本av免费视频播放| 亚洲综合色网址| 天天躁日日躁夜夜躁夜夜| 国产成人欧美| 啦啦啦 在线观看视频| 欧美日韩成人在线一区二区| 亚洲精品国产色婷婷电影| 超碰成人久久| 成年动漫av网址| 免费观看人在逋| 国产精品秋霞免费鲁丝片| 色视频在线一区二区三区| 最近的中文字幕免费完整| 亚洲av成人不卡在线观看播放网 | 侵犯人妻中文字幕一二三四区| 午夜影院在线不卡| 美女大奶头黄色视频| 国产精品人妻久久久影院| 咕卡用的链子| 免费少妇av软件| 国产熟女欧美一区二区| 最新在线观看一区二区三区 | 久久久久精品性色| 国产成人午夜福利电影在线观看| 黄片无遮挡物在线观看| 黄色毛片三级朝国网站| 高清在线视频一区二区三区| 嫩草影视91久久| 国产伦理片在线播放av一区| 69精品国产乱码久久久| 国产免费一区二区三区四区乱码| 天天躁狠狠躁夜夜躁狠狠躁| 男的添女的下面高潮视频| 在线观看三级黄色| 超碰成人久久| 国产午夜精品一二区理论片| 亚洲一区二区三区欧美精品| 观看美女的网站| 建设人人有责人人尽责人人享有的| 久久毛片免费看一区二区三区| 久久精品国产亚洲av高清一级| 亚洲成国产人片在线观看| 国产在视频线精品| 久久精品熟女亚洲av麻豆精品| 国产免费福利视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 久热爱精品视频在线9| 精品一品国产午夜福利视频| 视频在线观看一区二区三区| www.熟女人妻精品国产| 国产一级毛片在线| 夫妻性生交免费视频一级片| 爱豆传媒免费全集在线观看| 亚洲欧美一区二区三区国产| 国产不卡av网站在线观看| 亚洲一区二区三区欧美精品| 一区二区三区四区激情视频| 国产野战对白在线观看| 精品福利永久在线观看| 黑人欧美特级aaaaaa片| 人人妻,人人澡人人爽秒播 | 国产 一区精品| 午夜免费鲁丝| 熟妇人妻不卡中文字幕| 国产免费福利视频在线观看| √禁漫天堂资源中文www| 97人妻天天添夜夜摸| 80岁老熟妇乱子伦牲交| 香蕉国产在线看| 午夜免费鲁丝| 国产野战对白在线观看| 人体艺术视频欧美日本| 51午夜福利影视在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品国产精品| 国产成人系列免费观看| 极品少妇高潮喷水抽搐| 1024香蕉在线观看| 亚洲第一区二区三区不卡| 天天添夜夜摸| av.在线天堂| 美女主播在线视频| 国产亚洲最大av| 亚洲四区av| 在线亚洲精品国产二区图片欧美| 青春草视频在线免费观看| 天天添夜夜摸| 国产精品成人在线| 精品午夜福利在线看| 中文字幕av电影在线播放| 国产片特级美女逼逼视频| 狂野欧美激情性bbbbbb| 国产无遮挡羞羞视频在线观看| 亚洲av在线观看美女高潮| tube8黄色片| 亚洲精品一二三| 国产高清不卡午夜福利| 亚洲精品自拍成人| 性少妇av在线| 久久久久人妻精品一区果冻| 欧美日韩国产mv在线观看视频| 欧美成人午夜精品| 午夜福利视频在线观看免费| a级毛片在线看网站| 制服人妻中文乱码| 精品亚洲成a人片在线观看| 久久人妻熟女aⅴ| 国产激情久久老熟女| 亚洲欧美中文字幕日韩二区| 十八禁网站网址无遮挡| 久久久久视频综合| 亚洲精品成人av观看孕妇| 男女床上黄色一级片免费看| 国产精品亚洲av一区麻豆 | 欧美日韩一区二区视频在线观看视频在线| 久久久久精品国产欧美久久久 | 久久久久久久久免费视频了| 天天躁夜夜躁狠狠躁躁| 电影成人av| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久久久久免| 男女之事视频高清在线观看 | 一二三四中文在线观看免费高清| 成年动漫av网址| 超色免费av| 国产1区2区3区精品| 国产亚洲av片在线观看秒播厂| 国产女主播在线喷水免费视频网站| 黑人猛操日本美女一级片| 麻豆精品久久久久久蜜桃| 欧美激情极品国产一区二区三区| 亚洲综合色网址| 建设人人有责人人尽责人人享有的| 啦啦啦 在线观看视频| 亚洲免费av在线视频| 国产黄色视频一区二区在线观看| 伦理电影大哥的女人| 国产一区亚洲一区在线观看| 蜜桃国产av成人99| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 美女中出高潮动态图| 超色免费av| 午夜福利视频精品| 久久久精品免费免费高清| 午夜av观看不卡| 精品少妇内射三级| 欧美老熟妇乱子伦牲交| 国产精品亚洲av一区麻豆 | 亚洲美女视频黄频| 少妇 在线观看| 美女高潮到喷水免费观看| 爱豆传媒免费全集在线观看| 欧美 日韩 精品 国产| 美女脱内裤让男人舔精品视频| 午夜激情久久久久久久| 超碰97精品在线观看| 日日撸夜夜添| 男男h啪啪无遮挡| 欧美精品高潮呻吟av久久| 一级,二级,三级黄色视频| 日韩成人av中文字幕在线观看| 高清在线视频一区二区三区| 制服诱惑二区| 九色亚洲精品在线播放| 欧美日本中文国产一区发布| 亚洲国产中文字幕在线视频| 丝瓜视频免费看黄片| 欧美在线一区亚洲| 亚洲精品av麻豆狂野| 国产精品麻豆人妻色哟哟久久| 亚洲av中文av极速乱| 国产高清国产精品国产三级| 在线观看人妻少妇| 国产熟女午夜一区二区三区| a 毛片基地| 亚洲在久久综合| 国产精品亚洲av一区麻豆 | 日韩一卡2卡3卡4卡2021年| 国产黄色视频一区二区在线观看| 丝袜美腿诱惑在线| 校园人妻丝袜中文字幕| 少妇被粗大的猛进出69影院| 成人国产av品久久久| 男的添女的下面高潮视频| 日韩成人av中文字幕在线观看| 无遮挡黄片免费观看| 丝袜人妻中文字幕| 老司机影院成人| 亚洲国产日韩一区二区| 男女高潮啪啪啪动态图| 国产日韩欧美在线精品| 日韩不卡一区二区三区视频在线| 久久精品久久久久久久性| 一级毛片黄色毛片免费观看视频| 国精品久久久久久国模美| 精品视频人人做人人爽| 日本vs欧美在线观看视频| 一级a爱视频在线免费观看| 美女高潮到喷水免费观看| www.精华液| 久久久久久久大尺度免费视频| 久久99一区二区三区| 国产日韩一区二区三区精品不卡| 日韩 欧美 亚洲 中文字幕| 别揉我奶头~嗯~啊~动态视频 | 啦啦啦中文免费视频观看日本| 国产色婷婷99| 亚洲免费av在线视频| 在线观看www视频免费| 成年美女黄网站色视频大全免费| 中文字幕高清在线视频| 亚洲男人天堂网一区| 亚洲av福利一区| 亚洲精品国产av蜜桃| 日韩av不卡免费在线播放| 午夜日韩欧美国产| 捣出白浆h1v1| 中国国产av一级| 国产精品麻豆人妻色哟哟久久| 日本午夜av视频| 亚洲精品日本国产第一区| 九草在线视频观看| 9色porny在线观看| 久热这里只有精品99| 久久久精品国产亚洲av高清涩受| 亚洲精品中文字幕在线视频| 日韩一本色道免费dvd| 王馨瑶露胸无遮挡在线观看| 国产精品久久久人人做人人爽| 亚洲国产av影院在线观看| 亚洲精品第二区| 免费观看人在逋| 久久精品久久精品一区二区三区| 尾随美女入室| 欧美精品av麻豆av| 日韩精品免费视频一区二区三区| 一级毛片 在线播放| 亚洲精品国产区一区二| 十八禁人妻一区二区| 亚洲欧美清纯卡通| 波野结衣二区三区在线| 久久 成人 亚洲| 精品视频人人做人人爽| 久久久久国产一级毛片高清牌| 午夜精品国产一区二区电影| 日本色播在线视频| 韩国高清视频一区二区三区| 久久久亚洲精品成人影院| 欧美乱码精品一区二区三区| 亚洲中文av在线| 久久ye,这里只有精品| 精品亚洲成国产av| 精品一区二区三区四区五区乱码 | 久久精品国产a三级三级三级| 这个男人来自地球电影免费观看 | 91aial.com中文字幕在线观看| 美女国产高潮福利片在线看| 亚洲第一区二区三区不卡| 在线观看免费高清a一片| 午夜免费男女啪啪视频观看| 日韩av免费高清视频| 亚洲国产精品国产精品| 久久久精品国产亚洲av高清涩受| 18禁观看日本| 国产精品一国产av| 亚洲三区欧美一区| 亚洲精品日韩在线中文字幕| 亚洲国产最新在线播放| 久久久久久久久久久免费av| 高清不卡的av网站| 在线天堂中文资源库| 欧美精品av麻豆av| 国产精品国产av在线观看| 在线观看三级黄色| 搡老岳熟女国产| 成人国产麻豆网| 韩国高清视频一区二区三区| 久久99精品国语久久久| 亚洲精品成人av观看孕妇| 高清黄色对白视频在线免费看| 亚洲成人手机| 色94色欧美一区二区| 日本爱情动作片www.在线观看| 免费看不卡的av| 2018国产大陆天天弄谢| 巨乳人妻的诱惑在线观看| 国产精品秋霞免费鲁丝片| 老司机靠b影院| 亚洲一区二区三区欧美精品| 久久久久国产精品人妻一区二区| av网站免费在线观看视频| 亚洲精品国产av成人精品| 大香蕉久久网| 啦啦啦在线免费观看视频4| 亚洲精品一二三| 午夜久久久在线观看| 午夜影院在线不卡| 在线观看三级黄色| 欧美97在线视频| 一本色道久久久久久精品综合| 亚洲美女搞黄在线观看| 久久久久精品久久久久真实原创| 国产精品蜜桃在线观看| 中文字幕人妻熟女乱码| 亚洲熟女毛片儿| 日韩av不卡免费在线播放| 成年人午夜在线观看视频| 欧美日韩综合久久久久久| 久久 成人 亚洲| 捣出白浆h1v1| 日日摸夜夜添夜夜爱| 深夜精品福利| 天堂8中文在线网| 久久天堂一区二区三区四区| 欧美精品av麻豆av| 亚洲第一青青草原| 久久热在线av| 亚洲在久久综合| 黄片播放在线免费| 中文字幕人妻丝袜一区二区 | 免费高清在线观看日韩| 亚洲国产最新在线播放| 香蕉丝袜av| 天天躁日日躁夜夜躁夜夜| 人人妻人人澡人人爽人人夜夜| 99热网站在线观看| 丝袜喷水一区| 日日撸夜夜添| 国产成人系列免费观看| 成人三级做爰电影| 一区福利在线观看| 伦理电影免费视频| 18禁国产床啪视频网站| 久久热在线av| 日本黄色日本黄色录像| 亚洲色图综合在线观看| 性高湖久久久久久久久免费观看| 精品卡一卡二卡四卡免费| 午夜福利免费观看在线| 好男人视频免费观看在线| 女性生殖器流出的白浆| 久久精品国产综合久久久| 亚洲av电影在线进入| 精品少妇久久久久久888优播| 高清视频免费观看一区二区| 亚洲 欧美一区二区三区| 欧美人与善性xxx| 成人18禁高潮啪啪吃奶动态图| 久久毛片免费看一区二区三区| 丰满乱子伦码专区| 免费在线观看完整版高清| 欧美久久黑人一区二区| 亚洲精品美女久久av网站| av国产精品久久久久影院| 国产在线视频一区二区| 中文字幕高清在线视频| 香蕉丝袜av| 一本大道久久a久久精品| 国产av码专区亚洲av| 咕卡用的链子| 男女免费视频国产| 国产黄色视频一区二区在线观看| 欧美亚洲 丝袜 人妻 在线| 一边摸一边做爽爽视频免费| 欧美日韩视频精品一区| 久久人人97超碰香蕉20202| 日韩不卡一区二区三区视频在线| 涩涩av久久男人的天堂| 色综合欧美亚洲国产小说| 欧美在线黄色| 国产亚洲最大av| 1024香蕉在线观看| 精品少妇久久久久久888优播| 在线天堂中文资源库| 黄色毛片三级朝国网站| 国产精品av久久久久免费| 无限看片的www在线观看| 制服丝袜香蕉在线| 日韩 亚洲 欧美在线| 最新的欧美精品一区二区| 麻豆精品久久久久久蜜桃| 亚洲美女黄色视频免费看| 国产精品久久久久久久久免| 亚洲国产av影院在线观看| 国产日韩欧美亚洲二区| 成年女人毛片免费观看观看9 | 在线天堂中文资源库| 卡戴珊不雅视频在线播放| 夫妻性生交免费视频一级片| 老司机亚洲免费影院| 老鸭窝网址在线观看| 一级黄片播放器| 青草久久国产| 日韩制服丝袜自拍偷拍| 亚洲一级一片aⅴ在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品 欧美亚洲| 亚洲精品一二三| 午夜免费鲁丝| 黑丝袜美女国产一区| 亚洲四区av| 久久免费观看电影| 超碰97精品在线观看| 国产极品粉嫩免费观看在线| 极品少妇高潮喷水抽搐| 美女扒开内裤让男人捅视频| 免费在线观看完整版高清| 久久99热这里只频精品6学生| 看十八女毛片水多多多| av不卡在线播放| 精品福利永久在线观看| 在线观看www视频免费| av福利片在线| 丁香六月天网| 国产精品免费大片| 午夜福利,免费看| 日本爱情动作片www.在线观看| 女人久久www免费人成看片| 国产免费一区二区三区四区乱码| 国产成人精品无人区| 国产亚洲av高清不卡| 成人手机av| 中文字幕制服av| 久久这里只有精品19| 无遮挡黄片免费观看| www.精华液| 国产成人一区二区在线| 男男h啪啪无遮挡| 黄频高清免费视频| 少妇 在线观看| 国产 一区精品| 久久精品国产a三级三级三级| 国产一区亚洲一区在线观看| 国产一卡二卡三卡精品 | 少妇被粗大的猛进出69影院| 免费观看性生交大片5| 国产男女超爽视频在线观看| 侵犯人妻中文字幕一二三四区| 久久久久久久久久久久大奶| 亚洲av综合色区一区| 久久久精品国产亚洲av高清涩受| 欧美97在线视频| 成人18禁高潮啪啪吃奶动态图| 九草在线视频观看| 久久精品亚洲熟妇少妇任你| 只有这里有精品99| 亚洲精品一区蜜桃| 多毛熟女@视频| 亚洲精品美女久久久久99蜜臀 | 久久久国产欧美日韩av| 激情五月婷婷亚洲| 美女高潮到喷水免费观看| 婷婷色综合www| 亚洲精品自拍成人| 成年av动漫网址| 777米奇影视久久| 久久天堂一区二区三区四区| 91精品国产国语对白视频| 免费观看人在逋| 亚洲精品第二区| 99热国产这里只有精品6| 国产精品久久久av美女十八| 在线观看国产h片| 亚洲av日韩精品久久久久久密 | 精品国产超薄肉色丝袜足j| 一区二区av电影网| 精品人妻在线不人妻| 久久青草综合色| 免费女性裸体啪啪无遮挡网站| 亚洲精品第二区| 国产极品粉嫩免费观看在线| 久久狼人影院| av国产精品久久久久影院| av.在线天堂| 亚洲熟女毛片儿| 亚洲国产av影院在线观看| 在线观看免费视频网站a站| 黄色怎么调成土黄色| 丝袜美腿诱惑在线| 亚洲欧美激情在线| 精品一区二区三区四区五区乱码 | 悠悠久久av| 18在线观看网站| 成人国语在线视频| 久久99热这里只频精品6学生| 女性生殖器流出的白浆| 亚洲婷婷狠狠爱综合网| 亚洲成人av在线免费| 男男h啪啪无遮挡| 精品国产一区二区三区久久久樱花| 少妇猛男粗大的猛烈进出视频| 欧美最新免费一区二区三区| 国产免费一区二区三区四区乱码| 成人亚洲欧美一区二区av| 咕卡用的链子| 免费观看人在逋| 色婷婷久久久亚洲欧美| 国产成人精品久久久久久| 老司机亚洲免费影院| 欧美精品一区二区大全| 久久久久精品性色| 久久99一区二区三区| 高清欧美精品videossex| 波野结衣二区三区在线| 9191精品国产免费久久| 亚洲第一区二区三区不卡| 欧美成人午夜精品| 国产成人免费观看mmmm| 亚洲,欧美精品.| 一区二区三区乱码不卡18| 婷婷色麻豆天堂久久| 精品亚洲乱码少妇综合久久| 欧美精品av麻豆av| 最近中文字幕2019免费版| 99热全是精品| 成人国产av品久久久| 日本wwww免费看| 在线观看免费视频网站a站| 一区二区三区四区激情视频| 久久99一区二区三区| 美女视频免费永久观看网站| 亚洲精品乱久久久久久| 宅男免费午夜| 王馨瑶露胸无遮挡在线观看| 麻豆av在线久日| 亚洲av在线观看美女高潮| 亚洲精品国产色婷婷电影| 天天添夜夜摸| 国产亚洲一区二区精品| 日韩中文字幕视频在线看片| 亚洲国产日韩一区二区| 亚洲欧美中文字幕日韩二区| 制服诱惑二区| 美女中出高潮动态图| 一二三四中文在线观看免费高清| 久久久亚洲精品成人影院| 亚洲精品一区蜜桃| 成人漫画全彩无遮挡| 国产伦人伦偷精品视频| 如何舔出高潮| 亚洲精品国产一区二区精华液| 亚洲中文av在线| 少妇人妻精品综合一区二区| 一级毛片 在线播放| 国产精品二区激情视频| 男人爽女人下面视频在线观看| 久久精品久久久久久噜噜老黄| h视频一区二区三区| 大片免费播放器 马上看| 汤姆久久久久久久影院中文字幕| 国产亚洲av片在线观看秒播厂| 综合色丁香网| 最近最新中文字幕大全免费视频 | 成人亚洲精品一区在线观看| 国产黄色视频一区二区在线观看| 99热网站在线观看| 国产精品二区激情视频| 亚洲欧美中文字幕日韩二区| 制服诱惑二区| 亚洲精品久久成人aⅴ小说| 美女福利国产在线| 只有这里有精品99| 欧美97在线视频| 毛片一级片免费看久久久久| 十八禁人妻一区二区| 亚洲国产毛片av蜜桃av| 看免费av毛片| 亚洲激情五月婷婷啪啪| 无遮挡黄片免费观看| 啦啦啦中文免费视频观看日本| 久久免费观看电影| 午夜福利视频在线观看免费| 免费观看性生交大片5| 日韩人妻精品一区2区三区| 涩涩av久久男人的天堂| 久久青草综合色| 日韩av在线免费看完整版不卡| 日本91视频免费播放| 一区二区三区四区激情视频| 亚洲精品视频女| 波多野结衣av一区二区av| 亚洲,一卡二卡三卡| 婷婷色麻豆天堂久久| 国产日韩欧美在线精品| 国产成人免费无遮挡视频| 国产精品国产三级专区第一集| 一个人免费看片子| 久久 成人 亚洲| 亚洲精品视频女| 国产日韩欧美在线精品| 精品少妇内射三级| 欧美精品一区二区免费开放| 日韩一区二区视频免费看| av国产久精品久网站免费入址| 国产 一区精品| 亚洲人成网站在线观看播放| 免费在线观看视频国产中文字幕亚洲 | 国产午夜精品一二区理论片| 亚洲自偷自拍图片 自拍| bbb黄色大片|