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

    藥用植物寧夏枸杞 (Lycium barbarum L.) 土壤細(xì)菌群落演替特征

    2017-06-13 11:05:47張俊華張翼李明
    關(guān)鍵詞:幼齡樹齡老齡

    張俊華,張翼,李明

    藥用植物寧夏枸杞 (Lycium barbarum L.) 土壤細(xì)菌群落演替特征

    張俊華1,張翼2,李明1

    (1 寧夏大學(xué)環(huán)境工程研究院,寧夏銀川 750021;2 寧夏大學(xué)教育學(xué)院,寧夏銀川 750021)

    【目的】研究土層、季節(jié)和樹齡對寧夏枸杞土壤細(xì)菌群落結(jié)構(gòu)的影響,對揭示枸杞種植區(qū)域土壤質(zhì)量變化規(guī)律具有重要意義?!痉椒ā恳詫幭蔫坭街l(xiāng)—寧夏回族自治區(qū)中寧縣為研究區(qū),采用高通量測序分析了不同土層、季節(jié)和樹齡條件下枸杞土壤細(xì)菌生物量及細(xì)菌群落結(jié)構(gòu)和豐度的變化趨勢?!窘Y(jié)果】隨著枸杞樹齡增加,與幼齡 (種植當(dāng)年,< 1 年) 枸杞土壤相比,春季和夏季 0—20 cm 中齡 (6 年) 和老齡 (12 年) 土壤總有機(jī)碳(total organic carbon,TOC) 及 0—40 cm 土壤總有機(jī)氮 (total organic nitrogen,TON) 含量先增加后減少;春季不同樹齡 0—20 cm 土壤微生物生物量碳 (microbial biomass carbon,MBC)和微生物生物量氮(microbial biomass nitrogen,MBN)一直增加,但春季和秋季 20—40 cm 則呈相反趨勢。夏季中齡植株表層土壤 MBC 最高,但MBN 卻最低。隨著樹齡增加,枸杞表層土壤細(xì)菌多樣性普遍呈先增加后減小的趨勢。枸杞土壤變形菌門Proteobacteria、放線菌門 Actinobacteria、擬桿菌門 Bacteroidete 在細(xì)菌群落中占絕對優(yōu)勢。相同季節(jié)老齡土壤Proteobacteria 相對豐度普遍高于幼齡和中齡,幼齡土壤 Planctomycetes 和綠彎菌門 Chloroflexi 豐度高于中齡和老齡。變形桿菌中黃色單胞菌 Xanthomonadales、紅桿菌 Rhodobacterales 和根瘤菌 Rhizobiales 占主導(dǎo)地位,尤其是在秋季中齡和 3 個(gè)季節(jié)的老齡枸杞土壤。土層、季節(jié)和樹齡對土壤微生物數(shù)量影響不同?!窘Y(jié)論】隨著樹齡增加,寧夏枸杞表層土壤微生物生物量和細(xì)菌多樣性呈先增加后降低的趨勢。老齡植株土壤 Xanthomonadales數(shù)量相對最多。土層對枸杞土壤碳分布有極顯著影響,樹齡主要顯著影響土壤氮源和細(xì)菌群落多樣性,季節(jié)對枸杞土壤碳源、微生物量碳氮均有極顯著影響。

    寧夏枸杞;細(xì)菌群落;多樣性;季節(jié);土層;樹齡

    土壤微生物能夠影響許多生態(tài)過程、參與生物地球化學(xué)循環(huán),如有機(jī)質(zhì)和養(yǎng)分的轉(zhuǎn)化與循環(huán)[1],故其多樣性對生態(tài)功能的維持至關(guān)重要。植物種類在養(yǎng)分吸收、累積、分布和回歸中充當(dāng)了一個(gè)重要角色,而土壤微生物在控制植物養(yǎng)分有效性和土壤質(zhì)量方面具有重要作用[2]。藥用植物是對人類健康非常重要的資源。在微生物方面,許多學(xué)者對不同藥用植物根際土壤組分和微生物群體、根際效應(yīng)都做了較深入研究[3–4]。也有學(xué)者對不同耕作或施肥措施下藥用植物土壤有機(jī)質(zhì)、碳氮比、微生物多樣性和功能性指數(shù)、灌木菌根真菌、微生物生物量等的差異進(jìn)行了研究[5–6]。同時(shí),雖然許多證據(jù)顯示植物品種對根際區(qū)系、內(nèi)生菌根真菌有較大影響[7–8],但也有研究表明土壤因素 (尤其是土壤 pH)、種植年限或植物生長階段對微生物多樣性影響較植物品種更顯著[9]。微生物與植物根系防御和植物所需養(yǎng)分的獲取密切相關(guān),因而,有必要證明植物生長年限和發(fā)育階段對微生物的影響[10]。

    枸杞 (Lycium barbarum L.) 是茄科多年生灌木,其果實(shí)、葉片、根系都是傳統(tǒng)中藥。近年來,由于長期密集型管理,寧夏枸杞在 10~15 年就會出現(xiàn)果實(shí)產(chǎn)量減少、品質(zhì)下降的現(xiàn)象。因此,我們假設(shè)寧夏枸杞在生長前期和中期土壤微生物多樣性豐富,但隨著樹齡的增加,出現(xiàn)特有的土壤微生物種群并成為優(yōu)勢種群,或它們的數(shù)量達(dá)到了危害枸杞植株正常生長的程度。在此,本研究通過分析寧夏枸杞不同樹齡、不同季節(jié)、不同土層土壤微生物組成和多樣性,確定不同條件下枸杞土壤微生物細(xì)菌群落多樣性及引起土壤退化的微生物種群。

    1 材料與方法

    1.1 研究區(qū)概況

    寧夏回族自治區(qū)中衛(wèi)市中寧縣,位于東經(jīng) 105°15′~106°05′,北緯 36°49′~37°47′,是中國枸杞之鄉(xiāng),也是寧夏枸杞的主產(chǎn)區(qū)。該地區(qū)年平均氣溫 9.5℃,降水202.1 mm,日平均氣溫 ≥ 10℃ 的天數(shù)約 170 d,積溫 3200℃~3300℃。本研究在中寧縣寧安鎮(zhèn)南橋村選擇樹齡分別為幼齡 (種植當(dāng)年,< 1 年)、中齡(6 年) 和老齡 (12 年) 枸杞地各 3 塊 (每塊面積約666 m2)。該地區(qū)土壤以灌淤土為主, 質(zhì)地為砂壤土。

    1.2 樣品采集

    分別于春季 (枸杞展葉期,4 月 26 日)、夏季 (盛果期,7 月 14 日) 和秋季末期 (落葉期,10 月 25 日)采集土壤樣品。每塊地按 S 形選擇 9 株枸杞樹,在樹冠投影范圍內(nèi),先去除地表枯枝落葉后,在表層(0—20 cm) 和亞表層 (20—40 cm) 取其土壤,挖取具有完整根系的土體,先輕輕抖落大塊不含根系的土壤,然后用力將根表面附著的土壤全部抖落下來,獲得目標(biāo)土壤。每塊地 9 株枸杞分別分層采樣,然后將同一地塊相同層次的土壤充分混勻,將新鮮土樣分成 2 份,一份低溫冷藏帶回實(shí)驗(yàn)室,在–20℃ 下冷凍保存,用于土壤含水量和土壤微生物的測定,另一份帶回室內(nèi)風(fēng)干用于測定土壤 TOC、TON。

    1.3 樣品分析

    1.3.1 土壤樣品理化性狀測定 TOC 用重鉻酸鉀容量法;TON 用酸解法;微生物碳、氮 (MBC, MBN) 用氯仿熏蒸浸提法。

    1.3.2 土壤 DNA 提取 每個(gè)樣品取 0.5 g 土樣提取DNA,本試驗(yàn)采用 FastDNA?SPIN Kit For Soil 試劑盒提取土壤微生物基因組 DNA,按說明書操作步驟進(jìn)行。用 1% 的瓊脂糖凝膠檢驗(yàn) DNA 純度。DNA濃度用 NanoDrop (ND-2000) 檢測。

    1.3.3 illumina HiSeq 測序及數(shù)據(jù)分析 本研究目標(biāo)基因?yàn)榧?xì)菌的 16S rRNA V4 區(qū)基因片段,所用的引物為 515f (GTGCCAGCMGCCGCGGTAA)/806r (GGACTACHVGGGTWTCTAAT)。反應(yīng)體系:15μL Phusion High-Fidelity PCR Master Mix (New England Biolabs, UK),F(xiàn)orward Primers (0.2 μmol/L),Reverse Primers (0.2 μmol/L),10 ngDNA 模板。PCR 反應(yīng)條件:98℃ 1 min,30 × (98℃ 10 s, 50℃ 30 s, 72℃ 60 s),72℃ 5 min。PCR 產(chǎn)物用 2% 瓊脂糖凝膠檢測,選擇位于 200~300 bp 的目標(biāo)條帶,切膠回收純化進(jìn)行下一步的分析。

    用于 PCR 擴(kuò)增的引物外側(cè)都會加入一段 barcode序列進(jìn)行區(qū)分,序列為一段 6 bp 堿基的寡核苷酸序列,每個(gè)樣品用不同 barcode 序列進(jìn)行區(qū)分。Barcode 序列由上海吉美生物醫(yī)藥科技有限公司設(shè)計(jì)并合成。所得 PCR 產(chǎn)物經(jīng)切膠純化后 (Qiagen 膠回收試劑盒),用 NanoDrop 2000 (Thermo Fisher Scientific Inc. USA) 測定濃度,等摩爾混合后采用完成文庫構(gòu)建,安排上機(jī)測序。高通量測序由諾禾致遠(yuǎn) Illumina MiSeq 2500 平臺運(yùn)行 (Beijing, China, http://www.novogene.com/index.php),最后所測得數(shù)據(jù)在 250~300 bp 之間產(chǎn)生。

    序列分析流程如下:MiSeq 測序序列首先需要得到的雙端序列數(shù)據(jù)根據(jù) PE reads 之間的 overlap 關(guān)系,使用軟件 Seqprep 將成對的 reads 拼接 (merge)成一條序列,同時(shí)對 reads 的質(zhì)量和 merge 的效果進(jìn)行質(zhì)控過濾,并根據(jù)序列末端的 box 序列校正序列方向,然后按照 barcode 標(biāo)簽序列識別并區(qū)分樣品得到有效數(shù)據(jù)。

    1.4 數(shù)據(jù)處理

    利用 SAS 中 Duncan 法檢驗(yàn)不同樹齡、季節(jié)和土層枸杞土壤微生物指標(biāo)間的差異顯著性,用MANOVA 法研究土層、季節(jié)和樹齡對土壤理化指標(biāo)和微生物指標(biāo)的交互影響,采用 Pearson 相關(guān)進(jìn)行土壤理化性狀與土壤微生物指標(biāo)間的相關(guān)性分析。

    2 結(jié)果與分析

    2.1 土壤營養(yǎng)元素和微生物生物量

    隨著枸杞樹齡的增加,春季和夏季 0—20 cm 土壤 TOC 先增加后減少,夏季變幅較大 (表 1);秋季表層和春季、夏季亞表層土壤 TOC 有逐漸降低的趨勢;亞表層土壤 TOC 變幅普遍小于表層。不同樹齡土壤 TON 差異顯著;隨著樹齡的增加,0—40 cm 土壤 TON 均顯著增加后顯著降低,中齡土壤 TON 最高,其次為老齡。夏季土壤 TOC 和 TON 普遍高于春季和秋季,秋季最低;表層高于亞表層。春季表層土壤 MBC 和 MBN 隨樹齡的增加而增加,但春季和秋季亞表層則呈相反趨勢;夏季 MBC 隨樹齡的增加先增大后減小,中齡表層土壤 MBC 最高,但MBN 卻最低,亞表層該樹齡兩者均最低。秋季土壤MBC 相對最低。春季和夏季土壤不同樹齡 MBN 差異較小,但秋季中齡土壤 MBN 顯著高于其他 2 個(gè)樹齡。微生物生物量 C/N 變化范圍很大 (0.07~9.19),其中春季亞表層最高,秋季亞表層最低,幼齡樹土壤微生物生物量 C/N 較高 (夏季表層除外)。

    2.2 不同樹齡枸杞土壤細(xì)菌群落多樣性

    可操作分類單元 (operational taxonomic units,OTU) 的豐度初步說明了樣品的物種豐富程度,Chao1指數(shù)和 Shannon 指數(shù)則都指示了群落多樣性、豐富度和均勻度的變化。由表 2 可以看出,春季和夏季 0—20 cm 土壤中齡土壤表層序列數(shù)、OUTs 值、Chao1和 Shannon 指數(shù)最大,老齡土壤序列數(shù)、OUTs 值、Chao1 和 Shannon 指數(shù)普遍最??;該土層秋季老齡土壤 OTUs 和 Chao1 指數(shù)最大。20—40 cm 土壤序列數(shù)和各指數(shù)沒有一致的變化趨勢。夏季土壤 OTUs、Chao1 和 Shannon 指數(shù)大于春季和秋季,表層土壤OTUs、Chao1 和 Shannon 指數(shù)大于亞表層。

    2.3 不同樹齡枸杞土壤細(xì)菌群落相對豐度分析

    從門的角度來計(jì)算不同微生物相對豐度 (圖 1),表層土壤中 Proteobacteria 相對豐度最大,為 35.3%~48.4% (平均為 41.4%);其次為 Actinobacteria,平均豐度為 13.6%;Bacteroidetes、Planctomycetes 和Chloroflexi 平均豐度也大于 7%,而 Nitrospirae 相對最低,只有 0.8%。亞表層土壤中 Proteobacteria 相對豐度依舊最大,平均為 39.6%;其次是 Actinobacteria,相對豐度為 12.3%;Nitrospirae 平均相對豐度有所增加 (1.4%)。夏季表層土壤 Proteobacteria 豐度在三個(gè)季節(jié)中最高,春季和秋季差異很小,Actinobacteria豐度在秋季相對豐度最高,春季和夏季基本持平,此外,秋季 Chloroflexi 和 Planctomycetes 豐度高于春季和夏季,春季 Firmicutes 豐度高于夏季和秋季,其它門類豐度在三個(gè)季節(jié)表層土壤中無明顯差異。夏季亞表層土壤 Proteobacteria 豐度高于春季和秋季,秋季 Actinobacteria、Planctomycetes、Gemmatimonadetes、 Cloroflexi 均高于春季和夏季,春季 Crenarchaeota 高于夏季和秋季。相同季節(jié)老齡土壤 Proteobacteria 普遍豐度高于幼齡和中齡,幼齡土壤 Planctomycetes 和Chloroflexi 豐度高于中齡和老齡,春季幼齡 Bacteroidetes高于中齡和老齡。表層中齡土壤 Firmitutes 豐度最高,老齡最低;亞表層幼齡最高。表層幼齡土壤 Nitrospirae均高于其他兩個(gè)樹齡,其中夏季差異極顯著,中齡土壤 Nitrospirae 豐度最低,亞表層無明顯變化規(guī)律。2.4不同樹齡枸杞土壤群落結(jié)構(gòu)組成分析

    表1 土壤營養(yǎng)元素和微生物指標(biāo)Table1 Values of nutrient elements and microbial biomass for soils in different conditions

    表2 不同季節(jié)和土層各樹齡高通量序列數(shù)、OTUs 和 Shannon 指數(shù)Table2 Summary of high quality sequence numbers, numbers of operational taxonomic units (OTUs, ≥ 97% sequence similarity) and Shannon index in different conditions

    圖1 不同季節(jié)和樹齡土壤在 0—20 cm (左) 和 20—40 cm (右) 細(xì)菌各門類相對豐度圖Fig. 1 The relative abundance of bacterial phylum in 0–20 cm (Left) and 20–40 cm (Right) in the soil of L. barbarum in different seasons

    將微生物屬數(shù)量大于 0.005 作為標(biāo)準(zhǔn)做熱圖 (圖2),表層土壤篩選出 65 個(gè)屬,亞表層篩選出 62 個(gè)屬。表層土壤 Xanthomonadaceae 數(shù)量最多,尤其是在秋季中齡和 3 個(gè)季節(jié)的老齡枸杞土壤;Prielluflaceae含量次之,Candidatus、Flavobacteriaceae、o-iii1-15、Acidimicrobiales、Piscirickettsiaceae、Rhizobiales 和Bacullus 數(shù)量相對也較多。亞表層土壤微生物屬類組成與表層相似,但 Xanthomonadaceae、Flavobacteriaceae、Piscirickettsiaceae 和 Bacullus 較表層有不同程度的減少,Priellulaceae、Candidatus、o-iii1-15 和Rhizobiales 較表層有不同程度的增加。從聚類分析結(jié)果看,表層土壤秋季中齡和老齡土壤微生物數(shù)量及結(jié)構(gòu)相似,相同樹齡土壤春季和夏季微生物數(shù)量及結(jié)構(gòu)相似,亞表層土壤也有類似的變化趨勢。

    2.5 不同因素對土壤微生物量和細(xì)菌多樣性的影響

    從不同因素對土壤微生物指標(biāo)的影響程度來看(表 3),土層對土壤 MBC 和細(xì)菌群落 Shannon 指數(shù)分別有極顯著和顯著影響,季節(jié)對 MBC、MBN、MBC/OM、MBC/TN、細(xì)菌群落 Chao1 和 Shannon指數(shù)均有極顯著或顯著影響;樹齡對 MBN、MBN/ TN 和 Shannon 指數(shù)都有極顯著影響,季節(jié)和樹齡交互對 MBN、MBN/TN 有極顯著影響。土層和季節(jié)、土層和樹齡 2 因子交互及土層、季節(jié)和樹齡 3 因子交互對土壤微生物指標(biāo)均無顯著影響。所以,土層主要影響土壤碳分布,樹齡主要影響土壤氮分布,而季節(jié)對土壤碳和氮都有顯著影響;樹齡能夠極顯著影響細(xì)菌數(shù)量及其多樣性。

    2.6 細(xì)菌群落結(jié)構(gòu) PCA 分析

    圖3 為 0—20 cm (A) 和 20—40 cm (B) 土層不同季節(jié)和不同樹齡土壤微生物群落結(jié)構(gòu)的主成分分析結(jié)果,該結(jié)果可以提供土壤微生物群落結(jié)構(gòu)信息。樣品點(diǎn)距離的遠(yuǎn)近代表了樣品中微生物群落的相似性,距離越遠(yuǎn),差異越大。0—20 cm土層主成分 1 (PC1)和主成分 2 (PC2) 可以解釋 56.58% 的信息,由圖可知,不同季節(jié)幼齡土壤微生物結(jié)構(gòu)與中齡和老齡差異最大,而春季 3 個(gè)樹齡土壤微生物結(jié)構(gòu)與夏季和秋季差異也較大。與表層不同,秋季 3 個(gè)樹齡亞表層土壤微生物結(jié)構(gòu)與春季和夏季差異最大,不同季節(jié)幼齡土壤微生物與其他 2 個(gè)樹齡間差異也較大。從圖可以看出,亞表層土壤季節(jié)差異大于表層。說明樹齡和季節(jié)都可以改變枸杞土壤微生物群落結(jié)構(gòu)。

    3 討論

    3.1 枸杞土壤營養(yǎng)元素與微生物碳氮

    土壤有機(jī)碳 (TOC) 不僅是土壤養(yǎng)分的重要組成部分,而且也是微生物活性的物質(zhì)和能量源泉;土壤有機(jī)氮 (TON) 進(jìn)行礦化,形成的礦質(zhì)氮是作物吸

    圖2 土壤細(xì)菌在不同樹齡枸杞 0—20 cm (上圖) 和 20—40 cm (下圖) 土壤分布情況Fig. 2 Distribution of the selected bacteria in 0–20 cm (Up) and 20–40 cm (Down) soils under different age of L. barbarum[注(Note): 圖例中熱圖顏色表示不同樹齡土壤細(xì)菌屬類相對百分?jǐn)?shù)The heat map colors represent the relative percentages of the bacterial genus within different age of soils with the legend indicated of each panel.]

    收氮素的主要來源;土壤微生物量碳 (MBC) 是土壤有效養(yǎng)分的流動庫和土壤變化的媒介物,它在土壤保持肥力和植物養(yǎng)分吸收方面具有重要作用;而土壤微生物量氮 (MBN) 則是土壤有機(jī)氮的來源[11–12]。隨著寧夏枸杞樹齡的增加,春季和夏季表層土壤TOC 及 0—40 cm 土壤 TON 都先增加后降低,與Ying 等[13]結(jié)論一致,這是因?yàn)橹旋g枸杞處于旺盛生長期,雖然對養(yǎng)分消耗很大,但由于產(chǎn)量高,農(nóng)民對中齡枸杞施用的化肥和有機(jī)肥量大,以保證其高產(chǎn),故 TOC 和 TON 都處于最高值,同時(shí),由于植物根系在土壤中的分布影響有機(jī)碳含量的垂直分布,大量腐爛根系為土壤有機(jī)碳積累提供來源[14],故老齡枸杞土壤 TOC 和 TON 普遍高于幼齡枸杞土壤。春季表層土壤 MBC 和 MBN 隨樹齡的增加而增加[15],但春季和秋季亞表層則呈相反趨勢,這可能是由于當(dāng)?shù)剞r(nóng)民施肥普遍施在 10—20 cm,春季和秋季表層土壤養(yǎng)分豐富,能為微生物提供豐富的食物來源,而亞表層施肥量少,而該層枸杞根系對養(yǎng)分消耗量也較大,所以該層土壤養(yǎng)分缺乏,植株與微生物產(chǎn)生了養(yǎng)分競爭[16],導(dǎo)致主要以氮素為營養(yǎng)來源的微生物量下降,隨著樹齡的增加,MBC 和 MBN 逐漸減少,枸杞土壤環(huán)境逐步惡化[17]。由于夏季盛果期枸杞鮮果 7 天左右采摘一次,中齡植株產(chǎn)量高,對氮素的消耗非常大,故夏季亞表層中齡枸杞 MBC 和MBN 都最小,說明該季亞表層中齡植株氮素和碳源都需要補(bǔ)充。秋季枸杞果實(shí)采摘很少或不再采摘,所以土壤 MBN 大量累積,但中齡和老齡枸杞土壤MBC 均偏低,說明這些樹齡需要大量增加碳源,農(nóng)民需要多施用有機(jī)肥。此外,根系分布及根系分泌物在此過程中也起著積極作用[18–19]。

    表3 季節(jié)、土層、樹齡及其交互作用影響土壤微生物因子的多因素方差分析Table3 F-values form three-way ANOVA testing effects of season, layer, planting age, and their interactions on physicochemical and soil microbial parameters

    圖3 不同樹齡枸杞土壤 0—20 cm (A) 和 20—40 cm (B) 細(xì)菌群落結(jié)構(gòu)的 PCA 分析Fig. 3 Principal components analysis of bacteria communities in 0–20 cm (A) and 20–40 cm (B) in the soil of L. barbarum intercropping under different standing ages

    一般來說,土壤微生物群落的季節(jié)變化與土壤溫度和濕度的季節(jié)變異密切相關(guān)[20],故夏季 TOC、TON 和 MBC 高于春季,秋季最低,這是由于 4 月底當(dāng)?shù)販囟入m然有所升高,但是土溫依舊很低,10月底土壤溫度已經(jīng)下降至 15℃ 左右,這兩個(gè)時(shí)期植物生長緩慢,不利于土壤微生物的生存,而 7 月份,是植株生長茂盛時(shí)期,豐富的植物生物量為土壤提供大量的凋落物,另外發(fā)達(dá)的灌叢植物地下根系,其分泌物和死根是土壤微生物的重要能源物質(zhì),有利于土壤微生物的生長繁殖與土壤養(yǎng)分的積累。但何容等[21]對武夷山不同海拔植被土壤微生物量碳隨季節(jié)的變化研究表明均為早春最大,夏季最小,研究結(jié)果的差異源于海拔高帶來的溫濕度差異而致。因此,影響土壤微生物的因素較多,其生物量的高低仍取決于起關(guān)鍵作用的影響因素[22]。

    3.2 枸杞土壤微生物多樣性

    本研究中隨著樹齡增加,枸杞表層土壤 Chao1和 Shannon 指數(shù)普遍呈先增加后減小的趨勢,說明中齡枸杞土壤微生物群落結(jié)構(gòu)與幼齡和老齡相比最豐富,植株生長至 12 年左右時(shí)多樣性減小,土壤微環(huán)境惡化。但亞表層隨樹齡的增加微生物群落多樣性變化規(guī)律不同,中齡土壤 Chao1 和 Shannon 指數(shù)普遍最低,說明中齡植株土壤在亞表層微生物多樣性較差,這可能是施肥多在表層,亞表層養(yǎng)分含量相對較少,根系對養(yǎng)分吸收與微生物形成競爭,抑制或阻礙了微生物群落的繁殖[15]。

    在門的分類水平上,不同樹齡間活性微生物組成比較一致,但各門相對豐度不同。從細(xì)菌的組成來看,35% 以上的細(xì)菌基因源于 Proteobacteria,在活性微生物中占絕對優(yōu)勢,其他門如 Actinobacteria、Bacteroidetes、Planctomycetes、Chloroflexi、Firmicutes、Acidobacteria、Gemmatimonadetes、Nitrospirae 均為農(nóng)田土壤中的常見門類[23]。研究表明,Acidobacteria 廣泛分布于理化性質(zhì)迥異的各種土壤環(huán)境,是土壤中僅次于 Proteobacteria 的第二大微生物類群,其數(shù)量和組成受 pH 值影響較大[30],但本研究中 Actinobacteria 為該地區(qū)枸杞土壤微生物的第二大類群,這是因?yàn)槠溥m合生長在中性或微堿性的土壤中[24]。李晨華等[25]指出 Actinobacteria 與Proteobacteria 是灰漠土表層土壤中的優(yōu)勢類群,此結(jié)論與本研究基本一致,但該研究指出隨土壤深度的增加,Actinobacteria 相對豐度減少,Proteobacteria相對豐度增加,逐漸成為深層土壤的絕對優(yōu)勢類群,也有相關(guān)研究表明Proteobacteria 相對豐度會隨土壤深度降低,Actinobacteria 相對豐度則在深土層中較高[26]。而本研究發(fā)現(xiàn)隨著深度增加,Actinobacteria與 Proteobacteria 都略有降低。通常情況下放線菌適應(yīng)好氧環(huán)境,而周期性灌溉與研究區(qū)較高的地下水位,使得土壤含水量隨土壤深度的增加而升高,氧氣含量隨之下降,這可能是 Actinobacteria 隨深度其相對豐度降低的一個(gè)原因。有研究指出,Acidobacteria是可以作為較貧瘠土壤環(huán)境的指標(biāo)[25],中齡土壤Acidobacteria 最低,幼齡和老齡都較高,說明幼齡和老齡土壤質(zhì)量都較差。幼齡枸杞表層土壤硝化螺桿菌相對豐度最高,這表明了幼齡枸杞土壤的硝化作用較強(qiáng),從而可以促進(jìn)氧化亞氮的排放[27]。研究區(qū)中 Firmicutes 的絕對優(yōu)勢類群為芽孢桿菌綱(Bacilli),其在抗鹽堿方面具有優(yōu)勢[28],表層中齡樹土壤 Firmitutes 豐度最高,說明該樹齡表層抗鹽堿能力最強(qiáng),老齡抗鹽堿能力最弱。

    本研究變形桿菌中 Xanthomonadales、Rhodobacterales 和 Rhizobiales 占主導(dǎo)地位,尤其是在秋季中齡和 3 個(gè)季節(jié)的老齡枸杞土壤中。Xanthomonadales 是革蘭氏陰性菌,其中的很多物種是植物的病原體,能夠?qū)е罗r(nóng)作物病害[9],其含量與土壤遭受脅迫有關(guān)[29]。雖然表層土壤拮抗微生物L(fēng)ysobacter 比例高于亞表層,對病原菌有一定的制約作用[30],但含量較少,其作用也有限。亞表層Xanthomonadaceae、Bacullus 等減少,Priellulaceae、Rhizobiales 較表層有不同程度的增加,其中Rhizobiales 具有固氮作用[31],說明亞表層土壤較表層土壤略健康。

    3.3 枸杞土壤微生物影響因素

    從不同因素對土壤微生物指標(biāo)的影響程度來看,土層主要對土壤碳源和微生物生物量碳有極顯著影響[32]。季節(jié)對枸杞土壤碳源、微生物生物量碳氮均有極顯著影響,而樹齡主要顯著影響土壤氮源和細(xì)菌群落多樣性,土層、季節(jié)和樹齡 3 因子交互對土壤微生物均無顯著影響。從 PCA 結(jié)果也可以看出,季節(jié)和樹齡都可以顯著影響枸杞土壤微生物結(jié)構(gòu),表層樹齡影響更顯著,但亞表層季節(jié)影響更大。

    4 結(jié)論

    隨著寧夏枸杞樹齡的增加,春季和夏季表層土壤 TOC 及 0—40 cm 土壤 TON 都先增加后降低,春季表層土壤 MBC 和 MBN 隨樹齡的增加而增加,但春季和秋季亞表層則呈相反趨勢。隨著樹齡增加,枸杞表層土壤微生物多樣性普遍呈先增加后減小的趨勢,說明中齡枸杞土壤微生物群落結(jié)構(gòu)最豐富,植株生長至 12 年左右時(shí)多樣性減小,土壤微環(huán)境惡化。但 6 年樹齡植株根系在亞表層微生物多樣性較差,所以建議中齡枸杞施肥深度增加。枸杞土壤Proteobacteria、Actinobacteria、Bacteroidete 在活性微生物中占絕對優(yōu)勢。本研究變形桿菌中Xanthomonadales,Rhodobacterales 和 Rhizobiales 占主導(dǎo)地位,尤其是在秋季中齡和 3 個(gè)季節(jié)的老齡枸杞。不同因素對枸杞土壤細(xì)菌群落演替影響程度不同:季節(jié)對枸杞土壤碳源、微生物生物量碳氮均有極顯著影響,而樹齡主要顯著影響土壤氮源和細(xì)菌群落多樣性,土層、季節(jié)和樹齡 3 因子交互對土壤細(xì)菌群落均無顯著影響。

    [1]He JZ, Li J, Zheng YM. Thoughts on the microbial diversitystability relationship in soil ecosystems[J]. Biodiversity Science, 2013, 21(4): 411–420.

    [2]蔣婧, 宋明華. 植物與土壤微生物在調(diào)控生態(tài)系統(tǒng)養(yǎng)分循環(huán)中的作用[J]. 植物生態(tài)學(xué)報(bào), 2010, 34(8): 979–988. Jiang J, Song MH. Review of the roles of plants and soil microorganisms in regulating ecosystem nutrient cycling[J]. Chinese Journal of Plant Ecology, 2010, 34(8): 979–988.

    [3]Karthikeyan B, Jaleel CA, Lakshmanan GM A, et al. Studies on rhizosphere microbial diversity of some commercially important medicinal plants[J]. Colloids and Surfaces B: Biointerfaces, 2008, 62: 143–145.

    [4]Ahmed EA, Hassan EA, Tobgy KM KE, et al. Evaluation of rhizobacteria of some medicinal plants for plant growth promotion and biological control[J]. Annals of Agricultural Science, 2014, 59(2): 273–280.

    [5]Palese AM, Vignozzi N, Celano G, et al. Influence of soil management on soil physical characteristics and water storage in a mature rainfed olive orchard[J]. Soil & Tillage Research, 2014, 144: 96–109.

    [6]Taguas EV, Gómez JA. Vulnerability of olive orchards under the current CAP (common agricultural policy) regulations on soil erosion: a study case in Southern Spain[J]. Land Use Policy, 2015, 42: 683–694.

    [7]Zubek S, Stefanowicz AM, B?aszkowski J, et al. Arbuscular mycorrhizal fungi and soil microbial communities under contrasting fertilization of three medicinal plants[J]. Applied Soil Ecology, 2012, 59: 106–115.

    [8]Wei M, Zhuge YP, Lou YH, et al. Effects of fertilization on Xanthoceras sorbifolia bunge growth and soil enzyme activities[J]. Journal of Soil and Water Conservation, 2010, 24(2): 237–240.

    [9]Campos SB, Lisboa BB, Flavio AO, et al. Soil suppressiveness and its relations with the microbial community in aBrazilian subtropical agroecosystem under different management systems[J]. Soil Biology and Biochemistry, 2016, 96: 191–197.

    [10]Chaparro JM, Badri DV, Vivanco JM. Rhizosphere microbiome assemblage is affected by plant development[J]. The ISME Journal, 2014, 8(4): 790.

    [11]劉金山, 戴健, 劉洋, 等. 過量施氮對旱地土壤碳、氮及供氮能力的影響[J]. 植物營養(yǎng)與肥料學(xué)報(bào), 2015, 21(1): 112–120. Liu JS, Dai J, Liu Y, et al. Effects of excessive nitrogen fertilization on soil organic carbon and nitrogen and nitrogen supply capacity in dryland[J]. Journal of Plant Nutrition and Fertilizer, 2015, 21(1): 112–120.

    [12]Li R, Liu YN, Cheng WY, et al. Study on remediation of phenanthrene contaminated soil by pulsed dielectric barrier discharge plasma: The role of active species[J]. Chemical Engineering Journal, 2016, 296: 132–140.

    [13]Ying YX, Ding WL, Zhou YQ, et al. Influence of Panax ginseng continuous cropping on metabolic function of soil microbial communities[J]. Chinese Herbal Medicines, 2012, 4(4): 329–334.

    [14]Jobbagy EG, Jackson RB. The vertical distribution of soil organic carbon and its relation to climate and vegetation[J]. Ecological Applications, 2002, 10(2): 423–436.

    [15]王紀(jì)杰. 桉樹人工林土壤質(zhì)量變化特征[D]. 南京: 南京林業(yè)大學(xué)博士學(xué)位論文, 2011. Wang JJ. Variation characteristies of soil quality in EuealyPtus plantation [D]. Nanjing: PhD Dissertation, Nangjing Forestry University, 2011.

    [16]Singh S, Ghoshal N, Singh KP. Variations in soil microbial biomass and crop roots due to differing resource quality inputs in atropical dryland agroecosystem[J]. Soil Biology and Biochemisty, 2007, 39: 76–86.

    [17]Józefowska A, Wós B, Pietrzykowski M. Tree species and soil substrate effects on soil biota during early soil forming stages at afforested mine sites[J]. Applied Soil Ecology, 2016, 102: 70–79.

    [18]Shahzad T, Chenu C, Genet P, et al. Contribution of exudates, arbuscular mycorrhizal fungi and litter depositions to the rhizosphere priming effect induced by grassland species[J]. Soil Biology and Biochemistry, 2015, 80: 146–155.

    [19]Sun BJ, Jia SX, Zhang SX. Tillage, seasonal and depths effects on soil microbial properties in black soil of Northeast China[J]. Soil & Tillage Research, 2016, 155: 421–428.

    [20]張琪. 高寒草原區(qū)不同植被恢復(fù)方式對土壤微生物的影響研究[D]. 蘭州: 蘭州大學(xué)博士學(xué)位論文, 2015. Zhang Q. Influence of vegetation restoration methods on alpine steppe soil microorganisms [D]. Lanzhou: PhD Dissertation, Lanzhou University, 2015.

    [21]何容, 王國兵, 汪家社, 等. 武夷山不同海拔植被土壤微生物量的季節(jié)動態(tài)及主要影響因子[J]. 生態(tài)學(xué)雜志, 2009, 28(3): 394–399. He R, Wang GB, Wang JS, et al. Seasonal variation and its main affecting factors of soil microbial biomass under different vegetations along an elevation gradient in Wuyi Mountains of China[J]. Chinese Journal of Ecology, 2009, 28(3): 394–399.

    [22]Yao HY, Bowman D, Shi W. Seasonal variations of soil microbial biomass and activity in warm- and cool-season turfgrass systems[J]. Soil Biology and Biochemistry, 2011, 43: 1536–1543.

    [23]陳月星. 覆蓋對渭北旱作蘋果園土壤微生物群落特征的影響[D].楊凌: 西北農(nóng)林科技大學(xué)博士學(xué)位論文, 2015. Chen YX. Effects of mulching managements on soil microbial communities in an unirrigated apple orchard of Weibei Loess Plateau [D]. Yangling: PhD Dissertation, Northwest Agriculture and Forestry University, 2015.

    [24]劉方春, 邢尚軍, 馬海林, 等. 干旱生境中接種根際促生細(xì)菌對核桃根際土壤生物學(xué)特征的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2014, 25(5): 1475–1482. Liu FC, Xing SJ, Ma HL, et al. Effects of inoculating plant growthpromoting rhizobacteria on the biological characteristics of walnut(Juglans regia) rhizosphere soil under drought condition[J]. Chinese Journal of Applied Ecology, 2014, 25(5): 1475–1482.

    [25]李晨華, 張彩霞, 唐立松, 等. 長期施肥土壤微生物群落的剖面變化及其與土壤性質(zhì)的關(guān)系[J]. 微生物學(xué)報(bào), 2014, 54(3): 319–329. Li CH, Zhang CX, Tang LS, et al. Effect of long-term fertilizing regime on soil microbial diversity and soil property[J]. Acta Microbiologica Sinica, 2014, 54(3): 319–329.

    [26]肖燁. 三江平原典型濕地類型土壤微生物學(xué)特性對土壤有機(jī)碳的影響[D]. 北京: 中國科學(xué)院大學(xué)博士學(xué)位論文, 2015. Xiao Y. Effects of soil microbiological characteristics on soil organic carbon in the typical wetlands of Sanjiang Plain, Northeast China [D]. Beijing: PhD Dissertation, University of Chinese Academy of Sciences, 2015.

    [27]路丹, 雷靜, 韋燕燕, 等. 短期免耕和壟作對稻田土壤微生物群落及多樣性指數(shù)的影響[J]. 西南農(nóng)業(yè)學(xué)報(bào), 2015, 28(4): 1670–1674. Lu D, Lei J, Wei YY, et al. Effect of microbial community and diversity index of paddy soil under short no-tillage and ridge tillage[J]. Southwest China Journal of Agricultural Sciences, 2015, 28(4): 1670–1674.

    [28]Naushad HS, Gupta RS. Phylogenomics and molecular signatures for species from the plant pathogen-containing order xanthomonadalles[J]. PloS One, 2013, 8(2): e55216.

    [29]汪華. 茶園土壤微生物群落結(jié)構(gòu)對植茶年限、施肥和高溫條件的響應(yīng)研究[D]. 杭州: 浙江大學(xué)博士學(xué)位論文, 2014. Wang H. The response of soil microbial community structure to cultivating age, fertilization and high temperature in tea orchard [D]. Hangzhou: PhD Dissertation, Zhejiang University, 2014.

    [30]Rehman K, Ying Z, Andleeb S, et al. Short term influence of organic and inorganic fertilizer on soil microbial biomass and DNA in summer and spring[J]. Journal of Northeast Agricultural University (English Edition), 2016, 23(1): 20–27.

    [31]Dombrecht B, Marchal K, Vanderleyden J, et al. Prediction and overview of the RpoN-regulon in closely related species of the rhizobiales[J]. Genome Biology, 2002, 3(12): RESEARCH0076.

    [32]Eilers KG, Debenport S, Anderson S, et al. Digging deeper to find unique microbial communities: The strong effect of depth on the structure of bacterial and archaeal communities in soil[J]. Soil Biology and Biochemistry, 2012, 50: 58–65.

    Evolution characteristics of bacteria in soil of medical plant (Lycium barbarum L.)

    ZHANG Jun-hua1, ZHANG Yi2, LI Ming1
    ( 1 Institute of Environmental Engineering, Ningxia University, Yinchuan 750021, China; 2 College of Education, Ningxia University, Yinchuan 750021, China )

    【Objectives】Studying effects of soil layer, seasons and stand ages on soil bacteria of Lycium barbarum L. is very important to reveal trends of soil quality change in cultivated regions of L. barbarum.【Methods】The study was achieved by high throughput sequencing analysis of bacteria community structure and abundance of soil samples from L. barbarum fields in Zhongning County, Ningxia Province, which is the origin of L. barbarum. Three stand ages of L. barbarum were selected, small-age (< 1 year), middle-age (6 year) and old-age (12 year). The soil samples (three plots each) were collected in spring, summer and autumn in surface (0–20 cm) and subsurface (20–40 cm) layers.【Results】Compared to the small-age of soil, the total organic carbon (TOC) of 0–20 cm in spring and summer and total organic nitrogen (TON) of 0–40 cm of the middle-age and old-age were increased and then decreased and the microbial biomass carbon (MBC) andmicrobial biomass nitrogen (MBN) of 0–20 cm in spring were increased from the small-age to old-age, while those of 20–40 cm in spring and autumn were decreased. The MBC of 0–20 cm of the middle-age in summer was highest, while the MBN was lowest. With the increase of stand age, the soil bacteria community diversity of L. barbarum was increased and then decreased. Proteobacteria, Actinobacteria and Bacteroidete were prodominent plyumn. The relative abundance of Proteobacteria of the old-age was higher than that of the smallage and middle-age, Planctomycetes and Chloroflexi of small-age were higher than the middle-age and old-age. Xanthomonadales, Rhodobacterales and Rhizobiales were the prodominent classes, especial for the middle-age in autumn and old-age in three seasons. There were different effects of soil layer, seasons and stand ages on the microbial biomass.【Conclusions】With the increase of stand age, microbial biomass and diversities of soil of L. barbarum have an increased–decreased tendency. Compared to the small-age and middle-age, the soil of the old-age had maximum Xanthomonadales. Carbon of soil was influenced significantly by soil layers in our study, nitrogen and bacterial community diversities were influenced significantly by stand ages and carbon, MBC and MBN were all influenced significantly by seasons.

    Lycium barbarum L.; bacteria community; diversity; season; soil layer; stand age

    2016–09–02 接受日期:2016–11–28

    國家自然科學(xué)基金(41261080);寧夏自治區(qū)環(huán)保專項(xiàng)資金項(xiàng)目資助。

    張俊華(1977—),女,寧夏中衛(wèi)人,博士,研究員,主要從事精準(zhǔn)農(nóng)業(yè)和土壤質(zhì)量提升方面研究。E-mail:zhangjunhua728@163.com

    猜你喜歡
    幼齡樹齡老齡
    幼齡腰果園套種番薯栽培技術(shù)
    基于多光譜圖像的幼齡白木香含鋅量無損診斷
    《老齡科學(xué)研究》(月刊)歡迎訂閱
    老齡蘋果園“三改三減”技術(shù)措施推廣
    煙臺果樹(2021年2期)2021-07-21 07:18:46
    不同樹齡赤霞珠葡萄果實(shí)性狀及品質(zhì)的差異分析
    《老齡工作者之歌》——老齡工作者的心聲
    福利中國(2017年4期)2017-02-06 01:18:47
    巨尾桉不同樹齡及染蟲狀態(tài)桉葉油的成分分析
    高山幼齡果園套種蔬菜技術(shù)
    中國果菜(2016年9期)2016-03-01 01:28:45
    運(yùn)動對老齡腦的神經(jīng)保護(hù)作用研究進(jìn)展
    體育科研(2016年2期)2016-02-28 17:05:53
    閱讀花瓣雨
    91在线精品国自产拍蜜月| 国产精品女同一区二区软件| 狂野欧美激情性xxxx在线观看| 亚洲婷婷狠狠爱综合网| 成年av动漫网址| 久久av网站| 22中文网久久字幕| 免费av中文字幕在线| 性色avwww在线观看| 内地一区二区视频在线| 一个人看视频在线观看www免费| 日韩人妻高清精品专区| 一级毛片电影观看| 色视频在线一区二区三区| 欧美xxⅹ黑人| 久热久热在线精品观看| 最近2019中文字幕mv第一页| 国产精品久久久久成人av| 精品国产国语对白av| 亚洲情色 制服丝袜| av在线播放精品| 久久影院123| 国产综合精华液| 男女边吃奶边做爰视频| 少妇的逼好多水| 亚洲国产精品专区欧美| 美女大奶头黄色视频| 如日韩欧美国产精品一区二区三区 | 一边亲一边摸免费视频| 久久精品国产自在天天线| 汤姆久久久久久久影院中文字幕| 深夜a级毛片| 一级毛片 在线播放| 亚洲精品久久午夜乱码| 女的被弄到高潮叫床怎么办| 三级经典国产精品| av免费在线看不卡| av国产久精品久网站免费入址| 免费大片黄手机在线观看| 久久精品国产亚洲av涩爱| 三级国产精品欧美在线观看| 夫妻性生交免费视频一级片| 亚洲精品aⅴ在线观看| 国内揄拍国产精品人妻在线| 中国国产av一级| 乱系列少妇在线播放| 性色av一级| 亚洲av在线观看美女高潮| 日韩欧美一区视频在线观看 | 久久久久人妻精品一区果冻| 日韩亚洲欧美综合| 秋霞在线观看毛片| 极品教师在线视频| 国产在线男女| 国产免费福利视频在线观看| 最新中文字幕久久久久| 久久人人爽人人爽人人片va| 精品国产国语对白av| 国产欧美日韩一区二区三区在线 | 久久这里有精品视频免费| 欧美精品高潮呻吟av久久| 国产精品国产av在线观看| 在线精品无人区一区二区三| av福利片在线观看| 欧美精品一区二区免费开放| 男女边吃奶边做爰视频| 日韩欧美精品免费久久| 香蕉精品网在线| 午夜免费男女啪啪视频观看| 大陆偷拍与自拍| 成人毛片60女人毛片免费| 久久99蜜桃精品久久| 九九久久精品国产亚洲av麻豆| 久久韩国三级中文字幕| 国产深夜福利视频在线观看| 色94色欧美一区二区| 亚洲伊人久久精品综合| 精品一区二区三区视频在线| 97超碰精品成人国产| 日日啪夜夜撸| 成人18禁高潮啪啪吃奶动态图 | 乱人伦中国视频| 色94色欧美一区二区| www.色视频.com| 亚洲精品国产av成人精品| 欧美精品高潮呻吟av久久| 在线看a的网站| 九色成人免费人妻av| 少妇人妻精品综合一区二区| 免费大片黄手机在线观看| 成人国产麻豆网| 91精品伊人久久大香线蕉| 九九爱精品视频在线观看| 亚洲成色77777| 一级毛片 在线播放| 国产亚洲精品久久久com| 美女中出高潮动态图| 伊人亚洲综合成人网| 国产高清国产精品国产三级| 国产69精品久久久久777片| 一个人看视频在线观看www免费| 日韩伦理黄色片| 大片免费播放器 马上看| 日日啪夜夜爽| 久久久亚洲精品成人影院| 伦理电影免费视频| 成年av动漫网址| 精品国产一区二区三区久久久樱花| 女人精品久久久久毛片| 亚洲人与动物交配视频| 亚洲中文av在线| 日本vs欧美在线观看视频 | 一区二区三区精品91| 久久综合国产亚洲精品| 日本91视频免费播放| 亚洲欧美精品自产自拍| a级毛片免费高清观看在线播放| 麻豆精品久久久久久蜜桃| 伦精品一区二区三区| 69精品国产乱码久久久| 好男人视频免费观看在线| 精品亚洲成a人片在线观看| 日本黄大片高清| 久久久精品94久久精品| 男人狂女人下面高潮的视频| 你懂的网址亚洲精品在线观看| 精品午夜福利在线看| 欧美xxxx性猛交bbbb| 日韩成人伦理影院| 九九久久精品国产亚洲av麻豆| 亚洲四区av| 日韩电影二区| 国产极品粉嫩免费观看在线 | 国产精品国产三级国产专区5o| 国产在线视频一区二区| 亚洲自偷自拍三级| 热re99久久精品国产66热6| 黄色一级大片看看| 国产亚洲91精品色在线| 欧美 亚洲 国产 日韩一| 精品亚洲成国产av| 有码 亚洲区| 色婷婷av一区二区三区视频| 嫩草影院新地址| 免费黄色在线免费观看| 国产爽快片一区二区三区| 国产成人午夜福利电影在线观看| 在线天堂最新版资源| 国产精品福利在线免费观看| 五月玫瑰六月丁香| 男男h啪啪无遮挡| 亚洲精品第二区| 久久久久久久久久久丰满| 免费看光身美女| 菩萨蛮人人尽说江南好唐韦庄| 老司机亚洲免费影院| 亚洲欧美成人精品一区二区| 99久久中文字幕三级久久日本| 亚洲欧美中文字幕日韩二区| 少妇高潮的动态图| 免费av中文字幕在线| 极品教师在线视频| 日韩人妻高清精品专区| 一区二区三区乱码不卡18| a级一级毛片免费在线观看| 91aial.com中文字幕在线观看| 毛片一级片免费看久久久久| 午夜免费男女啪啪视频观看| 国产成人免费观看mmmm| 久久99热6这里只有精品| 九草在线视频观看| √禁漫天堂资源中文www| 国产精品伦人一区二区| 纵有疾风起免费观看全集完整版| 卡戴珊不雅视频在线播放| 午夜影院在线不卡| 日韩精品有码人妻一区| 在线免费观看不下载黄p国产| 精品人妻熟女毛片av久久网站| 七月丁香在线播放| 免费看av在线观看网站| 伊人亚洲综合成人网| av在线播放精品| 久久国产精品男人的天堂亚洲 | 女人久久www免费人成看片| 国产欧美另类精品又又久久亚洲欧美| 国产精品久久久久久精品古装| 最近中文字幕2019免费版| 国国产精品蜜臀av免费| 涩涩av久久男人的天堂| 日韩欧美 国产精品| 国产精品国产三级国产av玫瑰| 国内少妇人妻偷人精品xxx网站| 天天躁夜夜躁狠狠久久av| 我的女老师完整版在线观看| 18禁在线无遮挡免费观看视频| 亚洲成人一二三区av| 亚洲精品一二三| 婷婷色麻豆天堂久久| 自拍欧美九色日韩亚洲蝌蚪91 | 人妻 亚洲 视频| 亚洲欧美一区二区三区黑人 | 国产精品久久久久久久电影| 久久影院123| 久久久久久久大尺度免费视频| 亚洲欧美精品自产自拍| 天堂8中文在线网| 老女人水多毛片| freevideosex欧美| 亚洲综合精品二区| 乱码一卡2卡4卡精品| 久久久久久久大尺度免费视频| 国产免费视频播放在线视频| 夜夜骑夜夜射夜夜干| 日日啪夜夜爽| 亚洲av二区三区四区| 欧美激情极品国产一区二区三区 | 日本猛色少妇xxxxx猛交久久| 国产老妇伦熟女老妇高清| 日韩,欧美,国产一区二区三区| 欧美精品国产亚洲| 极品少妇高潮喷水抽搐| 亚洲天堂av无毛| av不卡在线播放| 免费看光身美女| 日韩av免费高清视频| 午夜精品国产一区二区电影| 亚洲精品,欧美精品| 美女cb高潮喷水在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产精品蜜桃在线观看| 国产白丝娇喘喷水9色精品| 成人毛片60女人毛片免费| 91成人精品电影| 亚洲欧美成人精品一区二区| 国产精品国产三级专区第一集| 免费久久久久久久精品成人欧美视频 | 亚洲内射少妇av| 日韩欧美精品免费久久| 精品人妻一区二区三区麻豆| 三级经典国产精品| 国产伦精品一区二区三区视频9| 精品少妇内射三级| 国产精品人妻久久久久久| 十分钟在线观看高清视频www | 亚洲精品aⅴ在线观看| 黄色毛片三级朝国网站 | 毛片一级片免费看久久久久| 精品少妇黑人巨大在线播放| 国产视频首页在线观看| 亚洲三级黄色毛片| 我的女老师完整版在线观看| 黑人巨大精品欧美一区二区蜜桃 | 简卡轻食公司| 亚洲av男天堂| 天堂8中文在线网| 国产精品一二三区在线看| 免费看不卡的av| 自拍偷自拍亚洲精品老妇| av一本久久久久| 久久久久国产精品人妻一区二区| 国产一区有黄有色的免费视频| 亚洲国产精品一区二区三区在线| 中文精品一卡2卡3卡4更新| 女性被躁到高潮视频| 秋霞伦理黄片| 欧美日韩精品成人综合77777| 亚洲欧美中文字幕日韩二区| 亚洲欧美一区二区三区黑人 | 免费看不卡的av| 久久久久视频综合| 久久精品国产亚洲av天美| 日产精品乱码卡一卡2卡三| 亚洲四区av| 人人妻人人澡人人爽人人夜夜| 9色porny在线观看| 亚洲精华国产精华液的使用体验| 免费久久久久久久精品成人欧美视频 | 中国三级夫妇交换| 九九爱精品视频在线观看| 国产亚洲5aaaaa淫片| 日韩亚洲欧美综合| 亚洲欧美一区二区三区黑人 | 男女啪啪激烈高潮av片| 国产精品久久久久久精品电影小说| 色94色欧美一区二区| 视频中文字幕在线观看| 91精品一卡2卡3卡4卡| 欧美精品国产亚洲| 久久国产精品男人的天堂亚洲 | 精华霜和精华液先用哪个| 亚洲欧美中文字幕日韩二区| 国产女主播在线喷水免费视频网站| 99九九在线精品视频 | 欧美 亚洲 国产 日韩一| 看十八女毛片水多多多| 久久 成人 亚洲| 国产伦在线观看视频一区| 国产精品三级大全| 精品一区二区三卡| 大片电影免费在线观看免费| 国产成人精品一,二区| 国产一区二区在线观看av| 成人二区视频| 亚洲av国产av综合av卡| 久久精品久久精品一区二区三区| 黄片无遮挡物在线观看| 各种免费的搞黄视频| 最后的刺客免费高清国语| 久久久久久久久久久丰满| 亚洲电影在线观看av| 精品久久久久久久久av| 久久99热6这里只有精品| 日韩熟女老妇一区二区性免费视频| 日日撸夜夜添| av女优亚洲男人天堂| 3wmmmm亚洲av在线观看| 综合色丁香网| 乱码一卡2卡4卡精品| 一级,二级,三级黄色视频| 日日摸夜夜添夜夜爱| 久久久久久久久久久丰满| 国产黄片美女视频| 精品午夜福利在线看| 久久久久精品性色| 午夜免费观看性视频| 两个人的视频大全免费| 嫩草影院新地址| 91成人精品电影| 亚洲一区二区三区欧美精品| 一级毛片久久久久久久久女| 妹子高潮喷水视频| 色网站视频免费| 一边亲一边摸免费视频| 色视频www国产| 久久国产精品大桥未久av | 欧美xxxx性猛交bbbb| 日本黄色日本黄色录像| av在线播放精品| 亚洲在久久综合| 美女xxoo啪啪120秒动态图| 80岁老熟妇乱子伦牲交| 97在线视频观看| 久久国产精品男人的天堂亚洲 | 亚洲欧洲国产日韩| 午夜精品国产一区二区电影| 人人妻人人澡人人看| 如何舔出高潮| 99热国产这里只有精品6| 丝袜在线中文字幕| 亚洲精品日韩av片在线观看| 久热久热在线精品观看| 欧美bdsm另类| 国产深夜福利视频在线观看| 久久久久久久久大av| 蜜桃久久精品国产亚洲av| 99热国产这里只有精品6| 女的被弄到高潮叫床怎么办| 国产一区二区在线观看日韩| 纯流量卡能插随身wifi吗| 精品久久久久久久久亚洲| 国产一级毛片在线| 国产熟女午夜一区二区三区 | 国产成人freesex在线| 欧美性感艳星| 精品少妇黑人巨大在线播放| 久久 成人 亚洲| 精品酒店卫生间| 一区二区三区四区激情视频| 日韩熟女老妇一区二区性免费视频| 久久97久久精品| 看免费成人av毛片| av国产精品久久久久影院| 精品少妇黑人巨大在线播放| 久久国产亚洲av麻豆专区| 女性生殖器流出的白浆| 精品久久久精品久久久| 制服丝袜香蕉在线| 人妻少妇偷人精品九色| 这个男人来自地球电影免费观看 | 尾随美女入室| av国产精品久久久久影院| 欧美日本中文国产一区发布| 午夜福利视频精品| 国产av码专区亚洲av| 亚洲va在线va天堂va国产| 国产精品久久久久久久电影| 亚洲国产色片| 午夜日本视频在线| 亚洲成人一二三区av| 国产熟女午夜一区二区三区 | 成人亚洲精品一区在线观看| 亚洲精品国产av成人精品| 亚洲电影在线观看av| 国产精品嫩草影院av在线观看| 午夜福利网站1000一区二区三区| 欧美xxxx性猛交bbbb| 日产精品乱码卡一卡2卡三| 久久精品久久久久久久性| 欧美精品一区二区免费开放| 精品99又大又爽又粗少妇毛片| 草草在线视频免费看| 国产日韩欧美亚洲二区| 2018国产大陆天天弄谢| 97超碰精品成人国产| 亚洲情色 制服丝袜| 亚洲在久久综合| 91久久精品国产一区二区成人| 久久女婷五月综合色啪小说| 99久国产av精品国产电影| 黑丝袜美女国产一区| 一区二区三区四区激情视频| 亚洲综合色惰| 午夜91福利影院| 精品少妇黑人巨大在线播放| 男女边吃奶边做爰视频| 精品人妻一区二区三区麻豆| 日产精品乱码卡一卡2卡三| 一区二区三区免费毛片| 另类亚洲欧美激情| 国产伦精品一区二区三区视频9| 麻豆精品久久久久久蜜桃| 精品熟女少妇av免费看| 色视频在线一区二区三区| 黄色欧美视频在线观看| 成人国产麻豆网| a级毛色黄片| 欧美日韩在线观看h| 美女视频免费永久观看网站| 亚洲综合色惰| 秋霞伦理黄片| 精品少妇久久久久久888优播| 岛国毛片在线播放| 美女cb高潮喷水在线观看| 日韩欧美 国产精品| 麻豆精品久久久久久蜜桃| 欧美日韩精品成人综合77777| 久久久久久久久久久丰满| 亚洲自偷自拍三级| 欧美老熟妇乱子伦牲交| 美女福利国产在线| 肉色欧美久久久久久久蜜桃| 欧美日韩av久久| 一区二区三区精品91| 大香蕉97超碰在线| 国产免费一级a男人的天堂| 乱系列少妇在线播放| 国产黄频视频在线观看| 久久精品国产a三级三级三级| 国产欧美另类精品又又久久亚洲欧美| 国产成人午夜福利电影在线观看| 日本爱情动作片www.在线观看| .国产精品久久| 男人和女人高潮做爰伦理| 日韩av在线免费看完整版不卡| 丰满饥渴人妻一区二区三| 26uuu在线亚洲综合色| 观看免费一级毛片| 国国产精品蜜臀av免费| 国产一区二区三区综合在线观看 | 99久久精品热视频| 亚洲第一av免费看| 18禁在线播放成人免费| 国产黄片美女视频| 亚洲国产精品成人久久小说| 国产 一区精品| 十八禁高潮呻吟视频 | 国产在线视频一区二区| 亚洲精品国产av成人精品| 日韩人妻高清精品专区| 亚洲国产成人一精品久久久| 十分钟在线观看高清视频www | 五月玫瑰六月丁香| 亚洲无线观看免费| 国产免费一区二区三区四区乱码| 22中文网久久字幕| 欧美日韩亚洲高清精品| 精品午夜福利在线看| a级片在线免费高清观看视频| 久久97久久精品| av.在线天堂| 欧美日韩av久久| 美女大奶头黄色视频| 亚洲国产精品成人久久小说| 欧美老熟妇乱子伦牲交| 欧美三级亚洲精品| 卡戴珊不雅视频在线播放| 日韩电影二区| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产亚洲av天美| 久久精品国产自在天天线| 超碰97精品在线观看| 日本与韩国留学比较| 亚洲精品第二区| 亚洲激情五月婷婷啪啪| 国产精品99久久99久久久不卡 | 这个男人来自地球电影免费观看 | 免费av不卡在线播放| 老司机亚洲免费影院| av播播在线观看一区| 十八禁高潮呻吟视频 | www.av在线官网国产| 亚洲精品国产av成人精品| 亚洲精品日韩在线中文字幕| 99热这里只有精品一区| 中国三级夫妇交换| 色5月婷婷丁香| av国产久精品久网站免费入址| 在线观看av片永久免费下载| 免费黄网站久久成人精品| 精品久久久久久久久av| 五月开心婷婷网| av女优亚洲男人天堂| 哪个播放器可以免费观看大片| 国产一区有黄有色的免费视频| 国产日韩一区二区三区精品不卡 | 国产永久视频网站| 三上悠亚av全集在线观看 | 亚洲在久久综合| 99久久精品国产国产毛片| 2022亚洲国产成人精品| 久久精品国产a三级三级三级| 欧美国产精品一级二级三级 | 人妻少妇偷人精品九色| 最新中文字幕久久久久| av专区在线播放| 91久久精品国产一区二区成人| 欧美精品一区二区大全| 亚洲人成网站在线播| 精品久久国产蜜桃| 九草在线视频观看| 水蜜桃什么品种好| 日韩制服骚丝袜av| 亚洲人成网站在线观看播放| 亚洲高清免费不卡视频| 亚洲激情五月婷婷啪啪| 日日啪夜夜爽| 18禁在线播放成人免费| 黄片无遮挡物在线观看| 内地一区二区视频在线| 黑人猛操日本美女一级片| 亚洲av二区三区四区| 少妇人妻精品综合一区二区| 九九久久精品国产亚洲av麻豆| 久久久欧美国产精品| 欧美日韩国产mv在线观看视频| 日产精品乱码卡一卡2卡三| videos熟女内射| 久久人人爽av亚洲精品天堂| 国产精品久久久久久久电影| 男人狂女人下面高潮的视频| 午夜激情久久久久久久| 五月开心婷婷网| 51国产日韩欧美| 日韩制服骚丝袜av| 久久亚洲国产成人精品v| 午夜91福利影院| 自拍偷自拍亚洲精品老妇| 一个人免费看片子| 综合色丁香网| 交换朋友夫妻互换小说| 亚洲图色成人| 欧美性感艳星| 三级国产精品欧美在线观看| 少妇熟女欧美另类| 日本午夜av视频| 亚洲成人手机| 十八禁高潮呻吟视频 | a 毛片基地| 男女边摸边吃奶| 国产精品一二三区在线看| 日日啪夜夜撸| 妹子高潮喷水视频| 日韩中字成人| 九九爱精品视频在线观看| 七月丁香在线播放| 一级毛片aaaaaa免费看小| 丰满乱子伦码专区| 国产午夜精品一二区理论片| 观看美女的网站| 一级毛片久久久久久久久女| 麻豆精品久久久久久蜜桃| 久久精品国产亚洲av天美| 欧美精品一区二区大全| 国产精品久久久久久av不卡| 少妇精品久久久久久久| 一级毛片黄色毛片免费观看视频| 国产黄片视频在线免费观看| 日韩一区二区视频免费看| 亚洲国产精品专区欧美| 蜜桃久久精品国产亚洲av| 亚洲熟女精品中文字幕| 妹子高潮喷水视频| 久久热精品热| 免费黄色在线免费观看| 91午夜精品亚洲一区二区三区| 欧美xxxx性猛交bbbb| 午夜福利在线观看免费完整高清在| 大片免费播放器 马上看| 免费大片黄手机在线观看| av国产精品久久久久影院| 在线播放无遮挡| 亚洲情色 制服丝袜| 在线播放无遮挡| av一本久久久久| 人妻 亚洲 视频| 有码 亚洲区| 自线自在国产av| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久精品电影小说| 亚洲成色77777| 一二三四中文在线观看免费高清| 国产伦在线观看视频一区| 精品久久国产蜜桃| 美女中出高潮动态图|