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

    毛烏素沙地苔蘚結(jié)皮對沙化土壤性質(zhì)和細菌群落的影響*

    2021-11-15 05:21:52張雨虹張韶陽張樹煇柳宇鑫李靜泉
    土壤學報 2021年6期
    關鍵詞:結(jié)皮苔蘚沙地

    張雨虹,張韶陽,張樹煇,柳宇鑫,趙 吉,李靜泉?

    毛烏素沙地苔蘚結(jié)皮對沙化土壤性質(zhì)和細菌群落的影響*

    張雨虹,張韶陽,張樹煇,柳宇鑫,趙 吉,李靜泉?

    (內(nèi)蒙古大學生態(tài)與環(huán)境學院,內(nèi)蒙古環(huán)境污染控制與廢物資源化重點實驗室,呼和浩特 010021)

    生物土壤結(jié)皮能有效防風固沙,對沙漠地區(qū)生態(tài)恢復意義重大,但目前缺乏對生物結(jié)皮層與同厚度裸沙表層之間以及它們各自下方土壤之間的比較與分析。采用熒光定量PCR和高通量測序技術,對毛烏素沙地苔蘚結(jié)皮層與裸沙表層,結(jié)皮下方土壤與裸沙下層土壤細菌16S rRNA基因豐度和群落結(jié)構(gòu)進行分析。結(jié)合土壤理化因子數(shù)據(jù),分析苔蘚結(jié)皮對沙地細菌群落多樣性的影響。對苔蘚結(jié)皮層和結(jié)皮下方土壤,以及裸沙地的相應表層和下層土壤,進行樣品采集,通過綜合分析更細致研究生物結(jié)皮對土壤細菌群落及其生存環(huán)境的影響。結(jié)果表明:相比于裸沙表層,苔蘚結(jié)皮顯著提高了土壤速效養(yǎng)分、全氮和有機質(zhì)的水平,結(jié)皮下方土壤中速效養(yǎng)分、全磷和有機質(zhì)含量均高于裸沙下層。結(jié)皮層及其下方土壤中粗粉砂和粉粒的含量顯著高于裸沙表層和下層,表明苔蘚結(jié)皮顯著改善了沙地土壤的理化性質(zhì)。有效磷、速效鉀、有效氮、黏粒、粉粒、粗砂是影響沙地土壤細菌群落組成的重要環(huán)境因子。苔蘚結(jié)皮層中16S rRNA基因拷貝數(shù)最高,結(jié)皮層及其下方土壤中的細菌豐度顯著高于裸沙表層和下層。多樣性分析顯示,結(jié)皮下方土壤中細菌多樣性最高。在苔蘚結(jié)皮層中,發(fā)現(xiàn)分類學地位尚未明確的難培養(yǎng)屬unclassifiled_f__Micromonosporaceae和norank_c__Cyanobacteria,以及(假諾卡氏菌屬)的相對豐度均顯著高于裸沙表層;在結(jié)皮下方土壤中,發(fā)現(xiàn)分類學地位尚未明確的難培養(yǎng)屬norank_c__Acidobacteria和(紅色桿菌屬)的相對豐度均顯著高于裸沙下層,這些差異顯著的物種對穩(wěn)定沙地土壤結(jié)構(gòu)具有重要作用。因此,苔蘚結(jié)皮的形成對沙地原有表層及下層土壤的細菌群落產(chǎn)生顯著影響,利于生物固沙。研究結(jié)果為風沙治理和荒漠生態(tài)恢復提供了重要的微生物學理論依據(jù)。

    苔蘚結(jié)皮;土壤理化性質(zhì);細菌群落;多樣性;對比分析

    第5次全國監(jiān)測數(shù)據(jù)顯示,我國沙化土地面積達到172.12萬km2,占國土面積的17.93%[1]。沙化土地主要集中在我國北方地區(qū),對當?shù)亟?jīng)濟和社會的發(fā)展,人們的生活以及生態(tài)平衡影響巨大。長久以來,我國一直致力于沙化土地生態(tài)恢復治理,尤其在沙漠與人類生活區(qū)域交界處,進行了大量的沙漠治理,阻止了沙漠對人類生存環(huán)境的進一步侵襲。在已使用的治沙方法中,生物土壤結(jié)皮(biological soil crusts)是其中效果明顯的一種。它是土壤碳氮循環(huán)的重要組成部分,不僅能提高土壤的穩(wěn)定性、抗風蝕和保持水分的能力,還能增加土壤肥力,提高植物吸收營養(yǎng)元素的能力,在沙漠生態(tài)系統(tǒng)中占有重要地位[2-6]。生物土壤結(jié)皮組成復雜,具有耐干旱和高溫等特點,其中尤屬苔蘚結(jié)皮的抗機械干擾、固沙和蓄水能力最強。

    微生物是土壤活性的重要組成部分,國內(nèi)外學者對生物土壤結(jié)皮的微生物群落特性進行了研究。Maier等[7]對西班牙塔伯納斯盆地的生物土壤結(jié)皮進行研究發(fā)現(xiàn),藍藻門(Cyanobacteria)在結(jié)皮層中豐度高而在結(jié)皮下方土壤中豐度低;相反地,結(jié)皮下方土壤中酸桿菌門(Acidobacteria)、疣微菌門(Verrucomicrobia)、芽單胞菌門(Gemmatimonadetes)、浮霉菌門(Planctomycetes)、裝甲菌門(Armatimonadetes)的含量明顯更高??屏_拉多高原地衣結(jié)皮中真菌多樣性明顯高于之前藍藻結(jié)皮中測定的結(jié)果;和采樣地對應的細菌相比真菌的生物量更少,多樣性更低;真菌多樣性和人為干擾呈負相關,和地衣結(jié)皮覆蓋呈正相關[8]。騰格里沙漠苔蘚結(jié)皮中的微生物生物量和微生物磷脂脂肪酸濃度高于藍藻-地衣結(jié)皮,生物土壤結(jié)皮的發(fā)育年齡與微生物生物量和微生物磷脂脂肪酸濃度呈正相關,并且顯著影響0~20 cm土壤層的微生物生物量碳和氮;生物土壤結(jié)皮在沙丘地帶的定殖和發(fā)育有益于植被恢復區(qū)的土壤微生物特性和土壤質(zhì)量[9]。賈麗娟等[10]研究發(fā)現(xiàn),與周圍臨近裸土相比,中央戈壁石下生物土壤結(jié)皮細菌群落結(jié)構(gòu)和多樣性明顯不同,生物結(jié)皮和臨近裸土中的細菌最優(yōu)勢菌門分別為Cyanobacteria和Actinobacteria,生物結(jié)皮中的優(yōu)勢屬為、和Cyanobacteria綱的1個未知屬,臨近裸土中的優(yōu)勢屬為、和3個屬于Actinobacteria門和Acidobacteria門的未知屬;與臨近裸土相比,生物結(jié)皮中有效磷、速效氮和葉綠素的含量均明顯升高,而pH則略微降低。目前有關生物結(jié)皮對土壤微生物群落影響的研究,多集中在生物結(jié)皮本身,或?qū)⑸锝Y(jié)皮與周圍裸地進行橫向比較,或與結(jié)皮下層土壤間進行縱向比較,或?qū)⒉煌愋徒Y(jié)皮(藻結(jié)皮、地衣結(jié)皮和苔蘚結(jié)皮)之間進行比較。故樣品采集多集中于生物結(jié)皮層、結(jié)皮下方土壤以及裸地土壤,而作為對照的裸地土壤并未區(qū)分表層和下層,缺乏對生物結(jié)皮層與同厚度裸地表層之間以及它們各自下方土壤之間的微生物群落組成與多樣性方面的細致比較與分析。本研究采集毛烏素沙地苔蘚結(jié)皮層和結(jié)皮下方土壤,以及裸沙地的相應表層和下層土壤,通過綜合分析更加細致地研究了生物結(jié)皮對土壤細菌群落及其生存環(huán)境的影響。研究結(jié)果為固沙治理和荒漠生態(tài)恢復提供理論依據(jù)。

    1 材料與方法

    1.1 土壤樣品采集

    樣地選在毛烏素沙地東南緣神木縣圪丑溝(38°53′44″~38°53′49″N,109°53′43″~109°53′46″E)。選取3塊油蒿周圍地面長有苔蘚結(jié)皮的樣地,分別采集結(jié)皮層(0~1.5 cm)和結(jié)皮層下方沙土(1.5~11.5 cm)的樣品;并且選取3塊同一區(qū)域的裸沙地作為對照,分別采集表層(0~1.5 cm)和下層(1.5~11.5 cm)的土樣。每塊樣地的大小為0.5 m×0.5 m,用鏟子沿對角線取3次混勻后作為一個樣品。苔蘚

    結(jié)皮3塊樣地呈直線分布,相互之間距離為10~15 m(裸沙樣地情況一樣)。不同土壤類型的命名分別為結(jié)皮層(Biological soil crusts,BC),結(jié)皮下方土壤(Layer beneath biological soil crusts,LBBC),裸沙表層土壤(Surface soil,SS),裸沙下層土壤(Layer beneath surface soil,LBSS)。采集的土壤樣品于采樣途中低溫保存,后續(xù)用于微生物多樣性測定的樣品于–80℃下保存;用于理化性質(zhì)測定的樣品自然風干后過60目篩。

    1.2 土壤理化性質(zhì)測定

    采用pH計()法測定土壤pH,碳酸氫鈉浸提—鉬銻抗分光光度法()測定土壤有效磷(AP),重鉻酸鉀容量法測定土壤有機質(zhì)(OM),堿解擴散法測定土壤有效氮(AN),硝酸鈉浸提—四苯硼鈉比濁法測定土壤速效鉀(AK),測定土壤全氮()、全磷。

    土壤粒度測定:稱量3 g土樣于100 mL燒杯中,加入10% H2O220 mL,于通風櫥內(nèi)控溫電熱板上加熱,溫度為80~120℃(最佳溫度90℃),H2O2耗盡時繼續(xù)加入H2O2,直至不再有氣泡。再加入10% HCl 10 mL,耗盡時繼續(xù)加入HCl至不再有氣泡。加滿蒸餾水靜置一夜,抽濾出水,除過量酸,如此反復數(shù)次,直至呈中性。抽濾出去水后,加入10 mL(NaPO3)6溶液,超聲震蕩10 min,用激光粒度儀(Mastersizer2000激光粒度分析儀)測定土壤粒度。具體指標包括粗砂(CS)、細砂(FS)、粗粉砂(CST)、粉粒(ST)和黏粒(CY)。

    1.3 土壤樣品16S rRNA基因擴增子高通量測序

    土壤樣本采集后送上海美吉生物醫(yī)藥科技有限公司進行16S rRNA基因擴增子高通量測序。具體過程為:采用通過引物338F(5′-ACTCCTACGGGA GGCAGCAG-3′)和 806R(5′-GGACTACHVGGG TWTCTAAT-3′)[11]對V3-V4可變區(qū)進行PCR擴增。 PCR產(chǎn)物用2%瓊脂糖凝膠回收純化,用Tris-HCl 洗脫后,用2%瓊脂糖電泳檢測。利用QuantiFluorTM-ST(Promega,USA)進行定量檢測。根據(jù) Illumina MiSeq 平臺(Illumina,SanDiego,USA)標準操作規(guī)程將純化后的擴增片段構(gòu)建 PE 2×300 的文庫。運用Illumina 公司的Miseq PE300 平臺進行測序。測序完成后,數(shù)據(jù)上傳美吉云平臺。登錄平臺(https://cloud.majorbio.com/)后,先進行基礎分析以獲得OUT豐度表,具體過程為:使用Trimmomatic軟件對原始測序序列進行質(zhì)量控制,并用FLASH軟件進行拼接。使用UPARSE 軟件(version 7.1 http://drive5.com/uparse/),以97%的相似度為依據(jù),對序列進行OTU聚類;并用UCHIME 軟件提出嵌合體。利用RDP classifier(http://rdp.cme.msu.edu/)對每條序列進行物種分類注釋,比對Silva數(shù)據(jù)庫(SSU123),將比對閾值設為70%。隨后在美吉云平臺進行交互分析,包括在OTU水平進行的排序分析和VPA分析,在門水平和屬水平進行的細菌群落組成分析,在屬水平進行的組間物種差異分析、相關性熱圖分析和共現(xiàn)性網(wǎng)絡分析。采用Spearman相關系數(shù)進行環(huán)境因子與物種之間的相關性Heatmap分析、Network網(wǎng)絡分析以及兩組組間物種差異顯著性Student’s T檢驗。

    1.4 熒光定量PCR擴增

    熒光定量PCR采用同上高通量測序的引物。20 μL反應體系為:10 μL 2 X ChamQ SYBR Color qPCR Master Mix,上下游引物(5 μmol·L–1)各0.8 μL,2 μL模板,0.4 μL 50 X ROX Reference Dye 1,6 μL ddH2O。擴增程序為:95℃預變性3 min;95℃變性5 s,58℃退火30 s,72℃延伸1 min。使用ABI7300型熒光定量PCR儀(Applied Biosystems,USA)進行擴增。每個樣品設置3個重復,以土壤干重計算最終基因豐度。

    1.5 試驗數(shù)據(jù)統(tǒng)計分析

    通過Excel2003和SPSS19.0對試驗數(shù)據(jù)進行統(tǒng)計分析。

    2 結(jié) 果

    2.1 結(jié)皮對土壤理化性質(zhì)的影響

    BC樣本組的各項土壤化學指標均總體上高于其他三個土壤樣本組(表1)。BC中有效磷、速效鉀、有效氮含量均顯著高于SS中的含量,可見結(jié)皮層能夠有效富集這些速效養(yǎng)分;相比于LBSS,LBBC中有效磷、速效鉀、有效氮的含量有升高趨勢。BC中全氮和有機質(zhì)的含量顯著高于SS中的含量,LBBC中全磷和有機質(zhì)含量高于LBSS。在4個樣本組中,BC的pH最低,其土壤略呈酸性。

    土壤結(jié)皮顯著改變了沙地土壤顆粒組成(表2),其中BC和LBBC相比于SS和LBSS有黏粒出現(xiàn),其組成比率分別為2.75%和0.18%;BC和LBBC中的粉粒含量分別為7.57%和5.04%,粗粉砂含量分別為42.28%和48.77%,2個指標均顯著高于SS和LBSS;而細砂的含量則顯著低于SS和LBSS,可見BC和LBBC的土壤樣本中顆粒粒徑有明顯變小的趨勢。生物結(jié)皮的形成顯著增加了土壤中細顆粒物含量。

    表1 不同樣本組的土壤化學性質(zhì)

    注:1)平均值±標準誤差,同列不同字母表示試驗處理間差異顯著(Duncan method,< 0.05)。2)BC,結(jié)皮層;LBBC,結(jié)皮下方土壤;SS,裸沙表層土壤;LBSS,裸沙下層土壤。下同。Note:1)Means ± standard error,different letters in the same row mean the significant difference between treatments (Duncan method,< 0.05). 2)BC,Biological soil crusts;LBBC,Layer beneath biological soil crusts;SS,Surface soil;LBSS,Layer beneath surface soil. AP,available phosphorus;AK,available potassium;AN,available nitrogen;TP,total phosphorus;TN,total nitrogen;OM,organic matter. The same below.

    表2 不同樣本組的土壤顆粒組成

    2.2 結(jié)皮對土壤細菌16S rRNA基因豐度的影響

    熒光定量PCR試驗結(jié)果顯示,BC、LBBC、SS、LBSS四個土壤樣本組的16S rRNA基因豐度的均值分別為1.79×1010copies·g–1干土、2.43×109copies·g–1干土、1.99×107copies·g–1干土、4.55×106copies·g–1干土。將基因豐度取對數(shù)作圖(圖1),可見BC和LBBC樣本組的細菌豐度均顯著高于SS和LBSS樣本組,并且BC樣本組也顯著高于LBBC樣本組。

    2.3 結(jié)皮對土壤細菌多樣性的影響

    LBBC的香農(nóng)指數(shù)和辛普森指數(shù)值在四組樣本中最大,其微生物多樣性在四組樣本中最高。LBSS的香農(nóng)指數(shù)和辛普森指數(shù)值居中。BC和SS的多樣性指數(shù)值非常接近,均低于LBBC和LBSS土壤樣本組(表3)。

    表3 不同樣本組的多樣性指數(shù)

    2.4 結(jié)皮對土壤細菌群落組成的影響

    在門(phylum)水平上繪制細菌群落組成柱形圖(圖2),將相對豐度小于0.01歸為others(圖中顯示相對豐度前13的細菌門)。各樣本組相對豐度排在前5的細菌門分別為:BC樣本組中,Actinobacteria(38.52%)、Proteobacteria(26.01%)、Acidobacteria(8.59%)、Cyanobateria(8.12%)、Chloroflexi(7.66%);SS樣本組中,Actinobacteria(32.22%)、Proteobacteria(27.18%)、Bacteroidetes(15.62%)、Acidobacteria(8.08%)、Chloroflexi(6.10%);LBBC樣本組中,Actinobacteria(31.01%)、Proteobacteria(23.72%)、Acidobacteria(16.20%)、Chloroflexi(13.43%)、Gemmatimonadetes(4.18%);LBSS樣本組中,Actinobacteria(33.40%)、Proteobacteria(30.84%)、Acidobacteria(10.00%)、Chloroflexi(8.94%)、Bacteroidetes(5.92%)。

    在屬(genus)水平上繪制細菌群落組成柱形圖(圖3),將相對豐度小于0.04歸為others(圖中顯示相對豐度前15的細菌屬)。各樣本組相對豐度排在前5的細菌屬分別為:BC樣本組中,unclassifiled_f__Micromonosporaceae(11.58%)、norank_c__Cyanobacteria(7.13%)、RB41(4.35%)、unclassified_o__Sphigomonoadales(3.18%)、(3.16%);SS樣本組中,(10.22%)、norank_o__Frankiales(4.72%)、RB41(4.25%)、(3.91%)、(3.22%);LBBC樣本組中,norank_c__Acidobacteria(7.59%)、RB41(5.30%)、(3.38%)、norank_c__Actinobacteria(2.84%)、(2.77%);LBSS樣本組中,(3.84%)、norank_o__Acidimicrobiales(3.54%)、RB41(3.51%)、(3.00%)、unclassified_f__ Micrococcaceae(2.70%)。

    2.5 結(jié)皮對土壤細菌群落組成組間差異的影響

    在PCA(Principal component analysis)分析中 PC1軸和PC2軸的貢獻率分別為 41.23%和27.73%,兩者貢獻率之和為68.96%。由于4個樣本組在PCA分析中的排序結(jié)果與RDA分析(圖4)相比沒有差異,故PCA分析排序結(jié)果參見圖4。BC在PC1軸上與其他三個土壤樣本組能夠區(qū)分明晰,表示BC與其他三個土壤樣本組的細菌群落組成存在明顯差異。BC與SS之間間隔最遠,兩者之間群落組成差異最大;LBBC與LBSS的樣本點分布分別呈兩個方向,說明在BC的影響下,LBBC的細菌群落組成與LBSS相比出現(xiàn)差異。

    2.6 環(huán)境因子與土壤細菌群落組成之間的關系

    圖4a為土壤化學性質(zhì)與細菌群落組成之間的RDA(Redundancy analysis)分析圖,RDA1軸和RDA2軸的解釋度分別為38.88%和16.63%,兩者之和的解釋度為55.51%。各環(huán)境因子對樣本細菌群落組成的影響程度為有效磷、速效鉀、有效氮的影響最大;有機質(zhì)、全磷和全氮影響次之;pH影響最小。圖4b為土壤顆粒組成與細菌群落組成之間的RDA分析圖,RDA1軸的解釋度為39.73%,RDA2軸的解釋度為25.58%,兩者之和的解釋度為65.31%。土壤顆粒組成對樣本細菌群落組成的影響程度為粉砂、黏粒、粗砂的影響最大,細砂次之,粗粉砂最小。

    在圖4a和圖4b中樣本點均分布在兩個區(qū)域,Ⅰ區(qū)為RDA1正半軸部分,SS和LBSS樣本點在此分布;Ⅱ區(qū)為RDA1負半軸部分,BC樣本點在此分布;LBBC樣本點在Ⅰ區(qū)和Ⅱ區(qū)都有分布。在圖4a中,Ⅰ區(qū)分布的SS和LBSS樣本組的細菌群落組成除與pH呈正相關外,與其他因子均呈負相關;Ⅱ區(qū)分布的BC樣本組的細菌群落組成除與pH值呈負相關外,與其他因子均呈正相關。在圖4b中,Ⅰ區(qū)分布的SS和LBSS樣本組的細菌群落組成與粗粉砂、粗砂、粉粒、黏粒呈負相關,與細砂呈正相關;Ⅱ區(qū)分布的BC樣本組的細菌群落組成與粗粉砂、粗砂、粉粒、黏粒呈正相關,與細砂呈負相關。

    對與BC樣本組細菌群落組成呈正相關的土壤化學因子(有效磷、速效鉀、有效氮、全磷、全氮、有機質(zhì))和土壤顆粒組成(粗砂、粗粉砂、粉粒、黏粒)之間進行VPA分析,結(jié)果顯示土壤顆粒組成的單獨解釋度為0.075,高于土壤化學因子的單獨解釋度0.057。對與BC樣本組細菌群落組成呈負相關的土壤化學因子(pH)和土壤顆粒組成(細砂)之間進行VPA分析,結(jié)果顯示細砂的單獨解釋度為0.238,高于pH的單獨解釋度0.078。綜合來看,在微生物組成影響方面,土壤顆粒組成的影響權(quán)重高于土壤化學因子。

    2.7 組間物種差異

    選取屬水平上相對豐度前15的物種進行組間差異性分析。圖5a為BC和SS樣本組之間物種差異分析??煽闯鯞C與SS比較中在unclassifiled_ f__Micromonosporaceae(屬于小單孢菌科的1個未知屬),norank_c__Cyanobacteria(屬于藍藻綱的1個未知屬)以及(假諾卡氏菌屬)方面存在極顯著差異,且均為BC中的相對豐度高于SS中的相對豐度。由圖5b可知,LBBC與LBSS兩個樣本組之間在(紅色桿菌屬)存在極顯著差異,在norank_c__Acidobacteria(屬于酸桿菌綱的1個未知屬)存在顯著差異,且均為LBBC中的相對豐度高于LBSS中的相對豐度。

    2.8 微生物種類與環(huán)境因子的相關性熱圖

    選取屬水平上相對豐度前18的物種與土壤化學性質(zhì)繪制相關性Heatmap圖。如圖6a所示,與SS相比在BC中相對豐度有顯著優(yōu)勢的unclassifiled_f__Micromonosporaceae、norank_c__ Cyanobacteria和均與AP 、OM、AK、AN、TP、TN呈正相關,都與pH呈負相關。其中,unclassifiled_f__Micromonosporaceae與TN呈極顯著正相關,與AP呈顯著正相關;norank_ c__Cyanobacteria與AP、AK、TP、TN均呈極顯著正相關;與AP呈極顯著正相關,與AK、AN呈顯著正相關。與LBSS相比在LBBC中相對豐度有顯著優(yōu)勢的norank_c__Acidobacteria,除與OM呈正相關外,與其他測定的土壤養(yǎng)分,包括AP 、AK、AN、TP、TN,均呈負相關;除與TN呈負相關外,與AP、OM、AK、AN、TP均呈正相關。

    選取屬水平上相對豐度前18的物種與土壤顆粒組成繪制相關性Heatmap圖。如圖6b所示,與SS和LBSS分別相比在BC和LBBC中相對豐度有顯著優(yōu)勢的5個屬(、unclassifiled_ f__Micromonosporaceae、norank_c__Cyanobacteria、、norank_c__Acidobacteria)均與FS(裸沙地的主要顆粒組成)呈負相關,其中與FS呈顯著負相關;并且該5個屬與ST(結(jié)皮土壤中含量顯著增多)均呈正相關。此外,、norank_c__Cyanobacteria、、norank_c__Acidobacteria這4個屬與CST(同樣是結(jié)皮土壤中含量顯著增多)均呈正相關,而且與CST呈極顯著正相關。此外,、unclassifiled_f__Micromonosporaceae、norank_c__Cyanobacteria、這4個屬與CY呈正相關,并且與CY呈顯著正相關。

    2.9 共現(xiàn)性網(wǎng)絡

    苔蘚結(jié)皮的形成影響B(tài)C和LBBC中細菌的群落組成,與裸沙地表層和下層的對照分別相比其相對豐度具有顯著優(yōu)勢的物種中,norank_c__ Acidobacteria在共現(xiàn)性網(wǎng)絡(圖7)中與其他細菌之間的聯(lián)系最密切,連接度(the degree of connections)最高,且多為正相關作用關系;unclassifiled_f__ Micromonosporaceae、norank_c__Cyanobacteria和與其他細菌之間的聯(lián)系程度次之,且unclassifiled_f__Micromonosporaceae和多為負相關作用關系,norank_c__ Cyanobacteria的正、負相關作用關系參半;與其他細菌之間的聯(lián)系程度最小。與SS相比在BC中相對豐度有顯著優(yōu)勢的unclassifiled_ f__Micromonosporaceae、norank_c__Cyanobacteria和三者之間互為正相關作用關系,與LBSS相比在LBBC中相對豐度有顯著優(yōu)勢的norank_c__Acidobacteria和之間亦為正相關作用關系。此外,放線菌門(藍色結(jié)點所示)在苔蘚結(jié)皮區(qū)域土壤中占有非常重要的地位。

    3 討 論

    3.1 苔蘚結(jié)皮對沙化土壤物理和化學性質(zhì)的影響

    苔蘚結(jié)皮層所測定的各項土壤養(yǎng)分指標均高于裸沙表層、結(jié)皮下方沙土和裸沙下層。從水平方向上看,生物土壤結(jié)皮的速效養(yǎng)分含量顯著高于裸沙表層,這是因為苔蘚結(jié)皮中的微生物及植物能夠有效地富集速效養(yǎng)分。生物土壤結(jié)皮的pH值略低于7,在4個樣本組中最低,可能是結(jié)皮層中的微生物能夠分泌酸性物質(zhì)所致[12-13]。生物結(jié)皮的形成使得土壤化學性質(zhì)發(fā)生改變,結(jié)皮中的微生物群落結(jié)構(gòu)與環(huán)境化學存在密切聯(lián)系[14]。當生物結(jié)皮將養(yǎng)分富集后,植被灌叢可以有效利用這些養(yǎng)分而得到生長,而且植物可以減少干擾并提供遮陰條件,促進生物結(jié)皮的生長[15],并且不同植被類型對生物結(jié)皮的發(fā)育影響不同[16]。在這樣的相互促進下,植物和結(jié)皮共同生長,起到固沙的作用。結(jié)皮下方土壤的有效磷、有效氮、速效鉀含量均高于裸沙下層土壤,可見在結(jié)皮層的影響下,結(jié)皮下方土壤中的速效養(yǎng)分含量有所增加。在土壤粒度組成方面,苔蘚結(jié)皮層與結(jié)皮下方土壤均有黏粒出現(xiàn),且粉粒和粗粉砂的含量分別高于裸沙表層和裸沙下層,說明生物土壤結(jié)皮的土壤粒徑小于裸沙土壤。大量粗粉砂的出現(xiàn)是生物結(jié)皮形成的前提條件和物質(zhì)基礎,而結(jié)皮層的發(fā)育又促進了粗粉砂風化,故而逐漸有了黏粒的出現(xiàn)[13]。生物結(jié)皮層具有明顯富集細顆粒物的能力,并隨著生物結(jié)皮層形成年份的延長逐漸增加,這種能力對改善沙土土壤理化性質(zhì),對沙土發(fā)育演變過程具有積極作用[17-18]。

    3.2 苔蘚結(jié)皮對沙化土壤細菌群落的影響

    細菌群落組成柱形圖分析結(jié)果顯示,在門水平上,結(jié)皮層樣本中Actinobacteria和Cyanobateria的相對豐度分別為38.52%和8.12%,明顯高于裸沙表層樣本;在屬水平上,結(jié)皮層樣本中unclassifiled_ f__Micromonosporaceae(屬于小單孢菌科的1個未知屬,歸于Actinobacteria門)和norank_c__Cyanobacteria(屬于藍藻綱的1個未知屬,歸于Cyanobacteria門)的相對豐度分別為11.58%和7.13%,也明顯高于裸沙表層樣本??梢姛o論從門水平還是屬水平上分析苔蘚結(jié)皮優(yōu)勢菌群的結(jié)果一致,其對生物固沙具有重要作用。結(jié)皮層的樣本點最為集中,說明其組內(nèi)樣本之間細菌群落組成差異最?。辉诮Y(jié)皮層影響下,結(jié)皮下方土壤組內(nèi)樣本之間細菌群落組成差異小于裸沙下層土壤。說明苔蘚結(jié)皮的形成可使土壤細菌群落組成組內(nèi)之間差異變小。苔蘚結(jié)皮的形成使結(jié)皮層及其下方土壤中的細菌豐度均顯著高于裸沙表層和下層,并且結(jié)皮層中的細菌豐度最高。在四組樣本中,結(jié)皮下方土壤細菌多樣性最高,而非結(jié)皮層土壤細菌多樣性最高。唐凱[19]對渾善達克沙地生物土壤結(jié)皮研究發(fā)現(xiàn)結(jié)皮層的細菌多樣性低于結(jié)皮下方土壤。從共現(xiàn)性網(wǎng)絡分析可知,在苔蘚結(jié)皮層相對豐度最高的unclassifiled_f__Micromonosporaceae(屬于小單孢菌科的1個未知屬)與其他細菌之間多為負相關作用關系,可能與小單孢菌科的微生物能夠分泌抗生素有關[20];而結(jié)皮下方土壤中相對豐度最高的norank_c__Acidobacteria(屬于酸桿菌綱的1個未知屬)與其他細菌之間多為正相關作用關系。土壤細菌群落中優(yōu)勢菌群與其他細菌之間的正、負相關作用關系可能會影響土壤細菌群落的多樣性。苔蘚結(jié)皮層的形成對沙地原有表層細菌群落產(chǎn)生重要影響,unclassifiled_f__Micromonosporaceae和norank_ c__Cyanobacteria的相對豐度明顯高于裸沙表層;并且在土壤理化性質(zhì)方面使其各項土壤養(yǎng)分指標均明顯高于裸沙表層,但對pH影響較小。從Heatmap分析可以看出,pH對4個樣本組中的細菌種類沒有顯著影響,而有效磷、速效鉀、全磷、全氮等因子對結(jié)皮樣本組中的優(yōu)勢細菌unclassifiled_f__ Micromonosporaceae和norank_c__Cyanobacteria具有顯著或極顯著影響。說明結(jié)皮的形成對細菌群落的影響大于pH。

    在之前的研究中,缺少對裸沙地進行表層和下層的分層取樣,故而在進行結(jié)皮與裸沙的對比分析時針對性不足。而從細菌群落組成分析中可知,裸沙表層和裸沙下層土壤之間存在明顯差異,將對照沙地進行分層取樣,對比分析才更準確。在毛烏素沙地苔蘚結(jié)皮的影響下,不同深度土壤細菌群落中的一些物種發(fā)生顯著變化。且unclassifiled_f__ Micromonosporaceae和norank_c__Cyanobacteria在苔蘚結(jié)皮中的相對豐度位于前兩位,norank_ c__Acidobacteria在結(jié)皮下方土壤中的相對豐度排在首位,進一步說明苔蘚結(jié)皮的形成對結(jié)皮層土壤及其下方土壤細菌群落中優(yōu)勢物種的形成具有重要影響。叢枝菌根真菌在生物土壤結(jié)皮中廣泛存在,對于荒漠生態(tài)系統(tǒng)具有重要作用,其群落多樣性與生物結(jié)皮之間關聯(lián)密切[21]。在后續(xù)研究中可考慮對裸沙和生物結(jié)皮進行分層取樣,更加細致地對比分析結(jié)皮形成對荒漠土壤叢枝菌根真菌群落的影響。

    3.3 苔蘚結(jié)皮及其下方土壤中豐度顯著增加的細菌的固沙能力

    與裸沙表層和下層土壤分別相比,在毛烏素沙地苔蘚結(jié)皮和其下方土壤中相對豐度顯著增加的物種,在固沙過程中應具有重要作用。其中,unclassifiled_f__Micromonosporaceae、和同屬于放線菌門(Actinobacteria),Actinobacteria是常見的土壤寄居者,有極強的黏附能力,可以成為貯菌源,其分泌出來的黏液可以將沙土黏結(jié)[22],其絲狀菌體也利于結(jié)皮層結(jié)構(gòu)的穩(wěn)定。并且Actinobacteria是土壤團聚體的重要組成部分,負責穩(wěn)定土壤團聚體的結(jié)構(gòu)[23-24]。細菌norank_c__Cyanobacteria屬于藍藻門(Cyanobacteria),Cyanobacteria在苔蘚結(jié)皮中占有重要地位[25-26],生物結(jié)皮演替早期階段的藻結(jié)皮具有顯著高于地衣結(jié)皮和苔蘚結(jié)皮的固氮能力,有助于荒漠生態(tài)系統(tǒng)進行原始養(yǎng)分積累;在維管植物大規(guī)模出現(xiàn)后,可以繼續(xù)固定大氣中的氮素,不斷積累土壤有機質(zhì),為周圍的維管植物和土壤微生物提供養(yǎng)分來源,促進了荒漠生態(tài)系統(tǒng)的正向演替[27-28]。此外,掃描電子顯微鏡觀察發(fā)現(xiàn)絲狀藍藻與胞外多聚分泌物形成一個錯綜復雜的網(wǎng)絡,能夠在土壤表面結(jié)合、附著礦物顆粒,增強了生物結(jié)皮的結(jié)構(gòu)穩(wěn)定性[29]。細菌norank_c__Acidobacteria屬于酸桿菌門(Acidobacteria),同Actinobacteria一樣,Acidobacteria也是土壤團聚體的重要組成部分,在穩(wěn)定土壤團聚體的結(jié)構(gòu)方面發(fā)揮著重要作用[23-24]。

    4 結(jié) 論

    沙地苔蘚結(jié)皮能有效富集速效養(yǎng)分和全效養(yǎng)分,并細化沙粒。有效磷、速效鉀、有效氮、黏粒、粉粒、粗砂是影響沙地苔蘚結(jié)皮區(qū)域細菌群落組成的重要環(huán)境因子。苔蘚結(jié)皮層中16S rRNA基因拷貝數(shù)最高,結(jié)皮下方土壤中微生物多樣性最高。在研究生物結(jié)皮的過程中,對照沙地也要分層取樣,微生物群落對比分析才更準確。苔蘚結(jié)皮的形成對結(jié)皮層土壤及其下方土壤細菌群落中優(yōu)勢物種的形成具有重要影響,這些物種對穩(wěn)定沙地土壤結(jié)構(gòu)具有重要作用。因此,苔蘚結(jié)皮的形成對沙地原有表層及下層土壤的細菌群落產(chǎn)生顯著影響,利于生物固沙。

    [1] Zhou Y,Yang X C,Jin Y X,et al. Classification of the desertification control models in North China[J]. Journal of Desert Research,2020,40(3):106—114. [周穎,楊秀春,金云翔,等. 中國北方沙漠化治理模式分類[J]. 中國沙漠,2020,40(3):106—114.]

    [2] Gao L Q,Bowker M A,Xu M X,et al. Biological soil crusts decrease erodibility by modifying inherent soil properties on the Loess Plateau,China[J]. Soil Biology & Biochemistry,2017,105:49—58.

    [3] Li X R,Tan H J,Hui R,et al. Researches in biological soil crust of China:A review[J]. Chinese Science Bulletin,2018,63(23):2320—2334. [李新榮,譚會娟,回嶸,等. 中國荒漠與沙地生物土壤結(jié)皮研究[J]. 科學通報,2018,63(23):2320—2334.]

    [4] Young M H,F(xiàn)enstermaker L F,Belnap J. Monitoring water content dynamics of biological soil crusts[J]. Journal of Arid Environments,2017,142:41—49.

    [5] Zhao L N,Liu Y B,Wang Z R,et al. Bacteria and fungi differentially contribute to carbon and nitrogen cycles during biological soil crust succession in arid ecosystems[J]. Plant and Soil,2020,447(1/2):379—392.

    [6] Chiquoine L P,Abella S R,Bowker M A. Rapidly restoring biological soil crusts and ecosystem functions in a severely disturbed desert ecosystem[J]. Ecological Applications,2016,26(4):1260—1272.

    [7] Maier S,Schmidt T S B,Zheng L J,et al. Analyses of dryland biological soil crusts highlight lichens as an important regulator of microbial communities[J]. Biodiversity and Conservation,2014,23(7):1735—1755.

    [8] Bates S T,Nash T H III,Sweat K G III,et al. Fungal communities of lichen-dominated biological soil crusts:Diversity,relative microbial biomass,and their relationship to disturbance and crust cover[J]. Journal of Arid Environments,2010,74(10):1192—1199.

    [9] Liu Y M,Li X R,Xing Z S,et al. Responses of soil microbial biomass and community composition to biological soil crusts in the revegetated areas of the Tengger Desert[J]. Applied Soil Ecology,2013,65:52—59.

    [10] Jia L J,Tang K,Lan H Q,et al. Diversity and structure of hypolithic bacteria community of Zhongyang Gobi[J]. Microbiology China,2018,45(12):2603—2613. [賈麗娟,唐凱,蘭慧青,等. 中央戈壁石下生物土壤結(jié)皮中細菌群落結(jié)構(gòu)和多樣性[J]. 微生物學通報,2018,45(12):2603—2613.]

    [11] Zhang B. Characteristics of coupled ultrasonic and alkaline pretreatment on enhancing anaerobic digestion of high-solid sludge[D]. Tianjin:Tianjin University,2017. [張博. 超聲耦合堿解預處理促進高固污泥厭氧消化的工藝特性研究[D]. 天津:天津大學,2017. ]

    [12] Kakeh J,Gorji M,Sohrabi M,et al. Effects of biological soil crusts on some physicochemical characteristics of rangeland soils of Alagol,Turkmen Sahra,NE Iran[J]. Soil & Tillage Research,2018,181:152—159.

    [13] Yang J Z,Bu C F,Zhang X C. Primary study on development characteristics of biological soil crust in Mu Us Sandy Land in Northern Shaanxi[J]. Journal of Soil and Water Conservation,2009,23(6):162—165,189. [楊建振,卜崇峰,張興昌. 陜北毛烏素沙地生物結(jié)皮發(fā)育特征的初步研究[J]. 水土保持學報,2009,23(6):162—165,189.]

    [14] Swenson T L,Karaoz U,Swenson J M,et al. Linking soil biology and chemistry in biological soil crust using isolate exometabolomics[J]. Nature Communications,2018,9:19.

    [15] Ghiloufi W,Chaieb M. Vascular plant diversity associated with biological soil crusts:Insights from Mediterranean arid ecosystem[J]. African Journal of Ecology,2017,55(2):252—255.

    [16] Wu Y S,Yin R P,Tian X M,et al. Development characteristics of biological crusts under artificial vegetation in southern Mu Us sandy land[J]. Journal of Desert Research,2018,38(2):339—344. [吳永勝,尹瑞平,田秀民,等. 毛烏素沙地南緣人工植被區(qū)生物結(jié)皮發(fā)育特征[J]. 中國沙漠,2018,38(2):339—344.]

    [17] Yan D R,Huang H G,Zhang S N,et al. Nutrients and particle composition characteristics in moss biological crusts[J]. Journal of Arid Land Resources and Environment,2018,32(10):111—116. [閆德仁,黃海廣,張勝男,等. 沙漠苔蘚生物結(jié)皮層養(yǎng)分及顆粒組成特征[J]. 干旱區(qū)資源與環(huán)境,2018,32(10):111—116.]

    [18] He F L,Guo C X,Wu H,et al. Effect of biological soil crust succession on soil texture,nutrient contents,and microbial populations of dune surfaces at the edge of the Minqin Oasis[J]. Acta Ecologica Sinica,2017,37(18):6064—6073. [何芳蘭,郭春秀,吳昊,等. 民勤綠洲邊緣沙丘生物土壤結(jié)皮發(fā)育對淺層土壤質(zhì)地、養(yǎng)分含量及微生物數(shù)量的影響[J]. 生態(tài)學報,2017,37(18):6064—6073.]

    [19] Tang K. Community structure and diversity of bacteria in biological soil crusts of Hunshandake deserts[D]. Huhhot:Inner Mongolia Agricultural University,2016. [唐凱. 渾善達克沙地生物土壤結(jié)皮細菌多樣性和群落結(jié)構(gòu)研究[D]. 呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學,2016. ]

    [20] Adnani N,Vazquez-Rivera E,Adibhatla S,et al. Investigation of interspecies interactions between marine Micromonosporaceae for identification of novel antibiotics[J]. Planta Medica,2015,81(11):913.

    [21] Qi J H,Liu Y B,Li X R,et al. AMF diversity analysis of lichen and moss biocrusts in Shapotou region[J]. Acta Pedologica Sinica,2020,57(4):986—994. [漆婧華,劉玉冰,李新榮,等. 沙坡頭地區(qū)地衣和蘚類結(jié)皮叢枝菌根真菌多樣性研究[J]. 土壤學報,2020,57(4):986—994.]

    [22] Cong L S. Investigation of bacterial community structure in biological soil crusts of Inner Mongolia deserts[D]. Huhhot:Inner Mongolia Agricultural University,2010. [叢立雙. 內(nèi)蒙古荒漠生物結(jié)皮中細菌群落分析[D]. 呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學,2010. ]

    [23] Mummey D,Holben W,Six J,et al. Spatial stratification of soil bacterial populations in aggregates of diverse soils[J]. Microbial Ecology,2006,51(3):404—411.

    [24] Zheng H,Wang X,Luo X X,et al. Biochar-induced negative carbon mineralization priming effects in a coastal wetland soil:Roles of soil aggregation and microbial modulation[J]. Science of the Total Environment,2018,610/611:951—960.

    [25] Maier S,Tamm A,Wu D M,et al. Photoautotrophic organisms control microbial abundance,diversity,and physiology in different types of biological soil crusts[J]. The ISME Journal,2018,12(4):1032—1046.

    [26] Chilton A M,Neilan B A,Eldridge D J. Biocrust morphology is linked to marked differences in microbial community composition[J]. Plant and Soil,2018,429(1/2):65—75.

    [27] Chamizo S,Cantón Y,Miralles I,et al. Biological soil crust development affects physicochemical characteristics of soil surface in semiarid ecosystems[J]. Soil Biology & Biochemistry,2012,49:96—105.

    [28] Zhang B C,Li R H,Xiao P,et al. Cyanobacterial composition and spatial distribution based on pyrosequencing data in the Gurbantunggut Desert,Northwestern China[J]. Journal of Basic Microbiology,2016,56(3):308—320.

    [29] Malam Issa O,Le Bissonnais Y,Défarge C,et al. Role of a cyanobacterial cover on structural stability of sandy soils in the Sahelian part of western Niger[J]. Geoderma,2001,101(3/4):15—30.

    Effect of Moss Crust on Sandy Soil Properties and Bacterial Community in Mu Us Sandy Land

    ZHANG Yuhong, ZHANG Shaoyang, ZHANG Shuhui, LIU Yuxin, ZHAO Ji, LI Jingquan?

    (School of Ecology and Environment of Inner Mongolia University, Inner Mongolia Key Laboratory of Environmental Pollution Prevention and Waste Resource Recycle, Huhhot 010021, China)

    【Objective】Biological soil crusts can be effective in wind-preventing and sand-fixing. This is of great significance to ecological restoration in desert areas. Nevertheless, there is no literature showing a comparative analysis between the biological soil crusts and the bare sand surface layer with the same thickness or for soils beneath them. 【Method】The abundance of 16S rRNA gene and the structure of bacterial communities of the moss crust layer, bare sand surface layer, and soil layers beneath them in Mu Us sandy land was analyzed by fluorescent quantitative PCR and high-throughput sequencing. Based on the data of soil physical and chemical factors, the effect of moss crust on the diversity of bacterial community in sandy land was analyzed. Also, the moss crust layer and soil layer beneath it, and the corresponding surface and underlying layers of bare sand were collected and the effects of biological crusts on soil bacterial communities and their living environments was analyzed more accurately through a comprehensive analysis. 【Result】Compared with bare sand surface layer, moss crust significantly increased the contents of soil available nutrients, total nitrogen, and organic matter. Also, the contents of available nutrients, total phosphorus, and organic matter in the soil beneath moss crust were higher than those in the soil beneath bare sand surface layer. The contents of coarse silt and silt in moss crust and its underlying soil were significantly higher than those in bare sand surface layer and its underlying soil. The results showed that moss crust significantly improved the physical and chemical properties of sandy soils. Additionally, the available phosphorus, potassium, nitrogen, clay, silt, and coarse sand were important factors affecting bacterial community composition of sandy soil. The moss crust layer had the highest 16S rRNA gene copies and the bacterial abundance of moss crust and its underlying soil was significantly higher than that of bare sand surface layer and its underlying soil. The results of diversity analysis showed that the soil beneath moss crust had the highest bacterial diversity. Futhermore, the relative abundances of unclassifiled_f__Micromonosporaceae, norank_c__Cyanobacteria, andinmoss crust were significantly higher than those in the bare sand surface layer. Also, the relative abundances of norank_c__ Acidobacteria andin the soil beneath moss crust were significantly higher than those in the soil beneath bare sand surface layer. These genera with significant differences play an important role in stabilizing the soil structure of sandy land. 【Conclusion】The formation of moss crust has a significant effect on the bacterial community of the original surface and underlying soils in sandy land. This is beneficial for biological sand-fixation. This study provides an important theoretical basis for sandstorm control and desert ecological restoration in Microbiology.

    Moss crust; Physical and chemical properties of soil; Bacterial community; Diversity; Comparative analysis

    S154.36

    A

    10.11766/trxb202012060221

    張雨虹,張韶陽,張樹煇,柳宇鑫,趙吉,李靜泉.毛烏素沙地苔蘚結(jié)皮對沙化土壤性質(zhì)和細菌群落的影響[J]. 土壤學報,2021,58(6):1585–1597.

    ZHANG Yuhong,ZHANG Shaoyang,ZHANG Shuhui,LIU Yuxin,ZHAO Ji,LI Jingquan. Effect of Moss Crust on Sandy Soil Properties and Bacterial Community in Mu Us Sandy Land[J]. Acta Pedologica Sinica,2021,58(6):1585–1597.

    *國家重點研發(fā)計劃項目(2017YFC0506703)、內(nèi)蒙古自治區(qū)重大基礎研究開放課題(30500-515330304)、內(nèi)蒙古自治區(qū)自然科學基金項目(2018MS03074)資助 Supported by the National Key Research and Development Program of China (No. 2017YFC0506703),the Opening Project for Major Basic Research of Inner Mongolia Autonomous Region (No. 30500-515330304),and the Natural Science Foundation of Inner Mongolia Autonomous Region (No. 2018MS03074)

    Corresponding author,E-mail:ljq198327@163.com

    張雨虹(1995—),女,內(nèi)蒙古呼和浩特人,碩士研究生,主要從事土壤微生物研究。E-mail:zhangyuhong950512@163.com

    2020–12–06;

    2021–03–03;

    2021–07–01

    (責任編輯:盧 萍)

    猜你喜歡
    結(jié)皮苔蘚沙地
    The mysterious “glacier mice”
    NaCl與Na2SO4風沙土鹽結(jié)皮力學特性比較分析?
    “小矮人”苔蘚的大用途
    軍事文摘(2021年18期)2021-12-02 01:28:36
    能鉆過柔軟沙地的蛇形機器人
    軍事文摘(2021年18期)2021-12-02 01:27:58
    黃土高原不同發(fā)育階段生物結(jié)皮的導水和持水特征
    呼倫貝爾沙地實現(xiàn)良性逆轉(zhuǎn)
    沙地迷宮
    黃土高原典型生物結(jié)皮對坡面產(chǎn)流產(chǎn)沙過程的影響
    風滾草
    苔蘚能“對付”空氣污染嗎?
    亚洲图色成人| 一级毛片我不卡| 成人黄色视频免费在线看| h视频一区二区三区| 亚洲综合精品二区| 精品一区在线观看国产| 夫妻午夜视频| 黄片无遮挡物在线观看| 成人免费观看视频高清| 岛国毛片在线播放| 亚洲国产毛片av蜜桃av| 亚洲三区欧美一区| 国产老妇伦熟女老妇高清| 日本欧美国产在线视频| 天堂8中文在线网| 精品第一国产精品| 色视频在线一区二区三区| 一级片'在线观看视频| 亚洲国产欧美在线一区| 一级,二级,三级黄色视频| 欧美精品一区二区大全| 一级毛片黄色毛片免费观看视频| 大陆偷拍与自拍| 在线天堂最新版资源| 亚洲男人天堂网一区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人一区二区在线| 亚洲成av片中文字幕在线观看 | 免费人妻精品一区二区三区视频| 国产精品一国产av| 久久久久久久久久久久大奶| 日韩av在线免费看完整版不卡| 日产精品乱码卡一卡2卡三| 日韩av在线免费看完整版不卡| 永久免费av网站大全| 亚洲欧美精品自产自拍| 成年美女黄网站色视频大全免费| 中文字幕另类日韩欧美亚洲嫩草| 国产av一区二区精品久久| 国产日韩欧美亚洲二区| 午夜激情av网站| 曰老女人黄片| 欧美日韩精品网址| 国产乱来视频区| 麻豆精品久久久久久蜜桃| 超碰97精品在线观看| 熟妇人妻不卡中文字幕| av线在线观看网站| videos熟女内射| 老女人水多毛片| videos熟女内射| 欧美97在线视频| 久久毛片免费看一区二区三区| a级片在线免费高清观看视频| 九草在线视频观看| 免费观看av网站的网址| 纵有疾风起免费观看全集完整版| 只有这里有精品99| av免费观看日本| 久久精品国产自在天天线| 美女国产视频在线观看| 欧美国产精品一级二级三级| 人人妻人人添人人爽欧美一区卜| 91精品伊人久久大香线蕉| 国产精品久久久久久av不卡| av电影中文网址| 在线免费观看不下载黄p国产| 妹子高潮喷水视频| 亚洲国产成人一精品久久久| 欧美成人午夜免费资源| 一级爰片在线观看| 午夜免费男女啪啪视频观看| 日韩在线高清观看一区二区三区| 亚洲在久久综合| 久久精品人人爽人人爽视色| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区欧美精品| 婷婷色综合大香蕉| 一本色道久久久久久精品综合| 中文字幕人妻丝袜一区二区 | 一区二区日韩欧美中文字幕| 久热这里只有精品99| 18+在线观看网站| 国产精品.久久久| 欧美日韩一区二区视频在线观看视频在线| 午夜福利在线观看免费完整高清在| 精品人妻偷拍中文字幕| 如日韩欧美国产精品一区二区三区| 久久久久精品久久久久真实原创| 中文字幕色久视频| 成人影院久久| 精品卡一卡二卡四卡免费| www日本在线高清视频| 婷婷色av中文字幕| videos熟女内射| 精品人妻在线不人妻| 久久精品人人爽人人爽视色| 亚洲国产精品一区二区三区在线| av在线观看视频网站免费| 亚洲欧洲国产日韩| 亚洲内射少妇av| 日本av手机在线免费观看| 毛片一级片免费看久久久久| 99久国产av精品国产电影| 日韩中文字幕欧美一区二区 | 多毛熟女@视频| 久久99热这里只频精品6学生| 日本欧美国产在线视频| 亚洲精品中文字幕在线视频| 免费黄频网站在线观看国产| 免费看不卡的av| 女人久久www免费人成看片| 寂寞人妻少妇视频99o| 91国产中文字幕| 中国三级夫妇交换| 精品第一国产精品| 夜夜骑夜夜射夜夜干| 国产日韩欧美视频二区| 999精品在线视频| 在线观看人妻少妇| 欧美变态另类bdsm刘玥| av有码第一页| 久久久亚洲精品成人影院| 亚洲av国产av综合av卡| 日本wwww免费看| 男女国产视频网站| 好男人视频免费观看在线| 久久人人爽人人片av| 久久精品国产a三级三级三级| 男人添女人高潮全过程视频| 90打野战视频偷拍视频| 国产黄色免费在线视频| 少妇人妻 视频| 国产日韩一区二区三区精品不卡| av卡一久久| 久久久久精品性色| 中文字幕制服av| 亚洲国产看品久久| 免费av中文字幕在线| 少妇的逼水好多| 亚洲在久久综合| 狠狠婷婷综合久久久久久88av| 精品酒店卫生间| 久久免费观看电影| 日韩一卡2卡3卡4卡2021年| 久久影院123| 欧美 亚洲 国产 日韩一| 女人被躁到高潮嗷嗷叫费观| 日韩一区二区视频免费看| 免费大片黄手机在线观看| 边亲边吃奶的免费视频| 亚洲欧美成人综合另类久久久| 国产精品一区二区在线不卡| 夜夜骑夜夜射夜夜干| 国产一级毛片在线| 午夜免费观看性视频| 丝袜美腿诱惑在线| 国产一区二区三区av在线| 亚洲综合色网址| 日日摸夜夜添夜夜爱| 在线免费观看不下载黄p国产| 熟妇人妻不卡中文字幕| 精品一品国产午夜福利视频| 久久久久久久久久久免费av| 国产精品国产三级专区第一集| 一本色道久久久久久精品综合| 少妇猛男粗大的猛烈进出视频| 国产乱来视频区| 亚洲第一av免费看| 波野结衣二区三区在线| 日本av免费视频播放| 亚洲欧洲日产国产| 国产精品久久久久成人av| 热re99久久精品国产66热6| 亚洲国产毛片av蜜桃av| 看十八女毛片水多多多| 亚洲av国产av综合av卡| 女人精品久久久久毛片| 午夜日韩欧美国产| 熟妇人妻不卡中文字幕| 国产精品免费大片| 天天躁日日躁夜夜躁夜夜| 成年av动漫网址| 十八禁网站网址无遮挡| 久久精品亚洲av国产电影网| 男女高潮啪啪啪动态图| 水蜜桃什么品种好| 美女午夜性视频免费| 亚洲欧美精品自产自拍| 可以免费在线观看a视频的电影网站 | 精品人妻一区二区三区麻豆| 成人亚洲精品一区在线观看| av福利片在线| 久久97久久精品| 国产精品久久久久久久久免| 91精品三级在线观看| 91aial.com中文字幕在线观看| 丰满乱子伦码专区| 又黄又粗又硬又大视频| 涩涩av久久男人的天堂| 国产高清不卡午夜福利| 国产麻豆69| 热re99久久国产66热| 精品99又大又爽又粗少妇毛片| 国产一区二区三区综合在线观看| 久久99蜜桃精品久久| 啦啦啦中文免费视频观看日本| 丝袜在线中文字幕| 中文字幕人妻丝袜一区二区 | 亚洲精品国产一区二区精华液| 美女xxoo啪啪120秒动态图| 夜夜骑夜夜射夜夜干| av线在线观看网站| 伦精品一区二区三区| 男的添女的下面高潮视频| 各种免费的搞黄视频| 美女福利国产在线| 青春草国产在线视频| 亚洲精品久久成人aⅴ小说| 亚洲第一区二区三区不卡| 99香蕉大伊视频| 亚洲 欧美一区二区三区| 老汉色∧v一级毛片| 国产又色又爽无遮挡免| 国产成人精品无人区| 免费不卡的大黄色大毛片视频在线观看| 啦啦啦视频在线资源免费观看| 晚上一个人看的免费电影| 欧美日韩精品网址| 侵犯人妻中文字幕一二三四区| 国产乱人偷精品视频| 美女高潮到喷水免费观看| 亚洲av国产av综合av卡| 男女高潮啪啪啪动态图| 一本大道久久a久久精品| 日日啪夜夜爽| 久久鲁丝午夜福利片| 1024视频免费在线观看| 一级片'在线观看视频| 国产精品熟女久久久久浪| 精品亚洲成国产av| videosex国产| 亚洲美女黄色视频免费看| 久久ye,这里只有精品| 九色亚洲精品在线播放| 久久国产精品男人的天堂亚洲| 国产有黄有色有爽视频| 欧美日韩一级在线毛片| 精品一区在线观看国产| 亚洲精品日本国产第一区| 亚洲国产毛片av蜜桃av| 欧美日韩综合久久久久久| 日本wwww免费看| 丰满少妇做爰视频| 久久久a久久爽久久v久久| 亚洲成色77777| 精品人妻偷拍中文字幕| 日韩三级伦理在线观看| 色婷婷久久久亚洲欧美| 极品少妇高潮喷水抽搐| 免费女性裸体啪啪无遮挡网站| 日本vs欧美在线观看视频| 久久精品国产综合久久久| 999久久久国产精品视频| 看免费成人av毛片| 黄片无遮挡物在线观看| 日韩成人av中文字幕在线观看| 晚上一个人看的免费电影| 欧美+日韩+精品| 成人黄色视频免费在线看| 成人漫画全彩无遮挡| 水蜜桃什么品种好| 久热这里只有精品99| 两个人免费观看高清视频| 999精品在线视频| av片东京热男人的天堂| 热99久久久久精品小说推荐| 搡老乐熟女国产| 成人手机av| 韩国高清视频一区二区三区| 中文字幕制服av| 女人久久www免费人成看片| 婷婷色综合大香蕉| 国产免费视频播放在线视频| 国产在线一区二区三区精| 中文字幕人妻丝袜制服| 爱豆传媒免费全集在线观看| 国产老妇伦熟女老妇高清| 午夜免费观看性视频| 免费人妻精品一区二区三区视频| 日本爱情动作片www.在线观看| 狠狠精品人妻久久久久久综合| 成人免费观看视频高清| 国产精品秋霞免费鲁丝片| 亚洲成人手机| 精品一区二区三区四区五区乱码 | 亚洲成人一二三区av| 欧美黄色片欧美黄色片| 亚洲人成电影观看| 精品第一国产精品| 亚洲成人一二三区av| 国产一区二区三区综合在线观看| 久久这里只有精品19| 不卡av一区二区三区| 蜜桃在线观看..| 婷婷色综合大香蕉| 一级毛片我不卡| 80岁老熟妇乱子伦牲交| 国产 一区精品| 精品国产国语对白av| 丝袜喷水一区| 1024视频免费在线观看| 久久99热这里只频精品6学生| 高清在线视频一区二区三区| 黄频高清免费视频| av一本久久久久| 香蕉精品网在线| 国产白丝娇喘喷水9色精品| 久久久久视频综合| 欧美精品亚洲一区二区| 亚洲精品一二三| 久久久久精品久久久久真实原创| 色播在线永久视频| 天天影视国产精品| 成人毛片a级毛片在线播放| 亚洲情色 制服丝袜| 久久久久久人人人人人| 国产av码专区亚洲av| 大码成人一级视频| 国产精品一二三区在线看| 久久久久久久久久久免费av| 大片电影免费在线观看免费| 中文乱码字字幕精品一区二区三区| 国产av一区二区精品久久| 亚洲熟女精品中文字幕| 一区二区三区四区激情视频| 女人高潮潮喷娇喘18禁视频| 国产精品 国内视频| 丰满迷人的少妇在线观看| 久久精品国产a三级三级三级| 制服丝袜香蕉在线| 久久久久人妻精品一区果冻| 香蕉精品网在线| 男女免费视频国产| 啦啦啦视频在线资源免费观看| 精品视频人人做人人爽| 一区在线观看完整版| 99久久精品国产国产毛片| 1024香蕉在线观看| 亚洲av成人精品一二三区| 最近2019中文字幕mv第一页| 久久鲁丝午夜福利片| 亚洲国产欧美在线一区| 激情视频va一区二区三区| 精品国产乱码久久久久久男人| 亚洲国产精品一区三区| 亚洲少妇的诱惑av| 男女边摸边吃奶| 国产1区2区3区精品| 国产色婷婷99| 国产精品无大码| 一区二区三区乱码不卡18| 欧美激情 高清一区二区三区| av福利片在线| 老司机影院成人| 欧美日本中文国产一区发布| 搡女人真爽免费视频火全软件| 国产 精品1| 国产男女超爽视频在线观看| √禁漫天堂资源中文www| 中文字幕制服av| 婷婷色麻豆天堂久久| 99久久综合免费| av视频免费观看在线观看| 视频区图区小说| 欧美日韩av久久| 日韩,欧美,国产一区二区三区| 日韩三级伦理在线观看| 热re99久久精品国产66热6| 欧美激情高清一区二区三区 | 日本爱情动作片www.在线观看| 久久久久久久亚洲中文字幕| 国产精品无大码| 老司机亚洲免费影院| 亚洲,欧美精品.| 大码成人一级视频| 婷婷色av中文字幕| 一级毛片黄色毛片免费观看视频| 少妇的丰满在线观看| 亚洲av综合色区一区| 老鸭窝网址在线观看| 国产免费一区二区三区四区乱码| 亚洲,欧美精品.| 啦啦啦在线观看免费高清www| 热re99久久精品国产66热6| 亚洲第一区二区三区不卡| 91精品伊人久久大香线蕉| 国产午夜精品一二区理论片| 成人国语在线视频| 性色av一级| 亚洲欧美精品自产自拍| 久久青草综合色| 国产精品香港三级国产av潘金莲 | 日韩欧美一区视频在线观看| 26uuu在线亚洲综合色| 国产片特级美女逼逼视频| 欧美激情极品国产一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 亚洲综合精品二区| 国产又色又爽无遮挡免| 成人国产av品久久久| 国产成人免费观看mmmm| 99热全是精品| 精品少妇黑人巨大在线播放| 亚洲色图 男人天堂 中文字幕| 精品久久久精品久久久| 精品一区二区三卡| 天天躁夜夜躁狠狠久久av| 亚洲国产欧美网| 深夜精品福利| 国产精品 欧美亚洲| 日本91视频免费播放| 人人妻人人澡人人看| 精品少妇一区二区三区视频日本电影 | 亚洲五月色婷婷综合| 在线天堂最新版资源| 啦啦啦中文免费视频观看日本| 热99久久久久精品小说推荐| 成人国语在线视频| 在线亚洲精品国产二区图片欧美| 国产精品 国内视频| 亚洲欧美一区二区三区久久| 欧美日本中文国产一区发布| 国产精品 欧美亚洲| 美女中出高潮动态图| 日韩制服骚丝袜av| 国产精品久久久久久久久免| 视频在线观看一区二区三区| 不卡av一区二区三区| 国产高清不卡午夜福利| 最近中文字幕高清免费大全6| 午夜福利一区二区在线看| 啦啦啦在线观看免费高清www| 国产乱人偷精品视频| 两性夫妻黄色片| 久久精品夜色国产| 男女免费视频国产| 精品卡一卡二卡四卡免费| 亚洲成av片中文字幕在线观看 | 狂野欧美激情性bbbbbb| 大陆偷拍与自拍| 中文字幕av电影在线播放| 亚洲,欧美精品.| www.精华液| 国产一区有黄有色的免费视频| 99热国产这里只有精品6| 国产精品国产三级国产专区5o| 黄频高清免费视频| 大片免费播放器 马上看| 十八禁网站网址无遮挡| 一级毛片电影观看| 国产免费视频播放在线视频| 国产精品不卡视频一区二区| 777米奇影视久久| 久久精品国产a三级三级三级| 男人操女人黄网站| 免费黄色在线免费观看| 亚洲欧美精品自产自拍| 色播在线永久视频| 日本黄色日本黄色录像| av一本久久久久| 国产又色又爽无遮挡免| 免费久久久久久久精品成人欧美视频| 丝袜喷水一区| 亚洲av男天堂| 黑人巨大精品欧美一区二区蜜桃| 免费在线观看黄色视频的| 在线观看www视频免费| 午夜日本视频在线| 97精品久久久久久久久久精品| 成人毛片60女人毛片免费| 一区二区三区乱码不卡18| 国产成人精品在线电影| 免费观看性生交大片5| 国产毛片在线视频| 中文字幕最新亚洲高清| 国产午夜精品一二区理论片| av网站在线播放免费| 一级爰片在线观看| 十分钟在线观看高清视频www| 韩国精品一区二区三区| 久久精品国产自在天天线| 妹子高潮喷水视频| 亚洲精品视频女| 下体分泌物呈黄色| 国产成人精品久久久久久| 精品卡一卡二卡四卡免费| 国产精品.久久久| 高清不卡的av网站| 久久久久精品人妻al黑| 久久综合国产亚洲精品| 国产又色又爽无遮挡免| 免费av中文字幕在线| 在线观看免费视频网站a站| 亚洲av福利一区| 在线观看一区二区三区激情| 日本vs欧美在线观看视频| 日韩制服丝袜自拍偷拍| 最黄视频免费看| 中文字幕制服av| 99热国产这里只有精品6| 18禁动态无遮挡网站| 69精品国产乱码久久久| 超碰97精品在线观看| 制服诱惑二区| 另类精品久久| 久久精品亚洲av国产电影网| 秋霞在线观看毛片| 电影成人av| 最近2019中文字幕mv第一页| 日韩人妻精品一区2区三区| 亚洲欧美一区二区三区国产| 国产精品 国内视频| 最近2019中文字幕mv第一页| 2022亚洲国产成人精品| 侵犯人妻中文字幕一二三四区| 亚洲第一区二区三区不卡| 又黄又粗又硬又大视频| 中文字幕人妻丝袜制服| 亚洲欧美清纯卡通| 国产精品av久久久久免费| 午夜日韩欧美国产| 国产精品av久久久久免费| 一级毛片我不卡| 狠狠婷婷综合久久久久久88av| av电影中文网址| 狠狠婷婷综合久久久久久88av| 夜夜骑夜夜射夜夜干| 成年美女黄网站色视频大全免费| 免费日韩欧美在线观看| 人成视频在线观看免费观看| 水蜜桃什么品种好| 欧美亚洲 丝袜 人妻 在线| 黄网站色视频无遮挡免费观看| 久久99精品国语久久久| 久久久久久免费高清国产稀缺| 亚洲色图 男人天堂 中文字幕| 亚洲精品中文字幕在线视频| 国产熟女午夜一区二区三区| 亚洲国产成人一精品久久久| 69精品国产乱码久久久| 天天躁夜夜躁狠狠躁躁| 国产麻豆69| 丁香六月天网| 在线观看三级黄色| 精品一品国产午夜福利视频| 丰满乱子伦码专区| 两个人免费观看高清视频| 久久久久人妻精品一区果冻| 免费在线观看完整版高清| 色视频在线一区二区三区| 国产成人精品福利久久| 在线观看免费高清a一片| 国产男女超爽视频在线观看| 高清欧美精品videossex| 精品国产乱码久久久久久男人| 日本欧美国产在线视频| 国产成人精品无人区| 最近2019中文字幕mv第一页| 国产极品天堂在线| 国产精品久久久久久精品古装| 99香蕉大伊视频| 高清在线视频一区二区三区| 少妇的丰满在线观看| 久久精品aⅴ一区二区三区四区 | 国产精品久久久久成人av| 在线看a的网站| 亚洲精品成人av观看孕妇| 大香蕉久久成人网| 中文字幕制服av| 一本—道久久a久久精品蜜桃钙片| 国产精品麻豆人妻色哟哟久久| 久久久国产一区二区| 成人午夜精彩视频在线观看| 91久久精品国产一区二区三区| 人人妻人人爽人人添夜夜欢视频| 校园人妻丝袜中文字幕| 午夜福利在线免费观看网站| 我要看黄色一级片免费的| 美女大奶头黄色视频| 妹子高潮喷水视频| 日韩精品有码人妻一区| 韩国精品一区二区三区| 国产精品亚洲av一区麻豆 | 免费少妇av软件| 男女高潮啪啪啪动态图| 日韩在线高清观看一区二区三区| 在线观看美女被高潮喷水网站| 伦理电影大哥的女人| 精品卡一卡二卡四卡免费| 亚洲欧美一区二区三区黑人 | 久久免费观看电影| 精品国产超薄肉色丝袜足j| 26uuu在线亚洲综合色| 亚洲精品一区蜜桃| 高清av免费在线| 国产毛片在线视频| 亚洲成人手机| 肉色欧美久久久久久久蜜桃| 97人妻天天添夜夜摸| 大香蕉久久成人网| 老熟女久久久| 欧美日韩亚洲国产一区二区在线观看 | 女性被躁到高潮视频| 亚洲av国产av综合av卡|