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

    典型網(wǎng)紋紅土結(jié)構(gòu)土壤細菌群落及其生態(tài)網(wǎng)絡(luò)特征①

    2022-12-09 01:09:48趙曉蕊吳華勇張甘霖
    土壤 2022年5期
    關(guān)鍵詞:網(wǎng)紋紅土群落

    趙曉蕊,吳華勇,張甘霖,3*

    典型網(wǎng)紋紅土結(jié)構(gòu)土壤細菌群落及其生態(tài)網(wǎng)絡(luò)特征①

    趙曉蕊1,2,吳華勇1,2,張甘霖1,2,3*

    (1 土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點實驗室(中國科學(xué)院南京土壤研究所),南京 210008;2 中國科學(xué)院大學(xué),北京 100049;3 中國科學(xué)院流域地理學(xué)重點實驗室(中國科學(xué)院南京地理與湖泊研究所),南京 210008)

    我國亞熱帶地區(qū)均質(zhì)紅土剖面下部的網(wǎng)紋紅土通常被認(rèn)為是第四紀(jì)古氣候和古環(huán)境變化的記錄載體。迄今為止,對土壤細菌群落在土壤結(jié)構(gòu)體尺度上的分異特征及其影響因素了解甚少。作為紅壤關(guān)鍵帶地下結(jié)構(gòu)的重要組成部分,網(wǎng)紋紅土獨特的紅–白相間的土體結(jié)構(gòu)可能引起土壤細菌群落在其上的差異分布。闡明網(wǎng)紋紅土結(jié)構(gòu)土壤細菌群落特征對探索紅壤關(guān)鍵帶深部生態(tài)系統(tǒng)的生物多樣性和多功能性具有重要的生態(tài)學(xué)意義。本研究以典型網(wǎng)紋紅土結(jié)構(gòu)體為研究對象,將其拆分為紅色基質(zhì)和白色網(wǎng)紋兩部分,利用16S rRNA高通量測序分析紅/白部分細菌群落與細菌生態(tài)網(wǎng)絡(luò)的特征。結(jié)果顯示,網(wǎng)紋結(jié)構(gòu)體紅白部分細菌群落組成差異顯著,紅色基質(zhì)部分放線菌門和厚壁菌門等細菌顯著富集;白色網(wǎng)紋部分酸桿菌門、硝化螺旋菌門等細菌的相對豐度更高。網(wǎng)絡(luò)分析結(jié)果表明,白色網(wǎng)紋部分的細菌網(wǎng)絡(luò)節(jié)點和連線較多,聚類系數(shù)較低,說明白色網(wǎng)紋部分細菌之間相互作用關(guān)系更復(fù)雜。典范對應(yīng)分析表明,網(wǎng)紋結(jié)構(gòu)體環(huán)境異質(zhì)性造成了細菌群落組成在網(wǎng)紋紅白部分的差異分布。其中,質(zhì)地是決定細菌群落及其生態(tài)網(wǎng)絡(luò)差異最重要的環(huán)境因子。此外,在網(wǎng)紋的紅色基質(zhì)部分,礦質(zhì)養(yǎng)分的含量也與細菌群落的組成密切相關(guān)。

    紅壤關(guān)鍵帶;網(wǎng)紋紅土結(jié)構(gòu);土壤細菌群落;細菌生態(tài)網(wǎng)絡(luò);影響因素

    土壤微生物是驅(qū)動陸地表層生態(tài)系統(tǒng)即關(guān)鍵帶(Earth’s critical zone)系統(tǒng)生物地球化學(xué)循環(huán)的引擎,是陸地生態(tài)系統(tǒng)的重要組成部分[1]。土壤微生物調(diào)控土壤中的物質(zhì)周轉(zhuǎn)與養(yǎng)分傳遞,影響地表植物多樣性與作物生產(chǎn),驅(qū)動著陸地生態(tài)系統(tǒng)諸多生物地球化學(xué)過程[2-3]。然而,關(guān)于土壤微生物群落特征的研究大多集中在近地表環(huán)境,如表層土壤或耕作層土壤,對于深層土壤微生物的認(rèn)識仍較少[4]。研究表明,盡管深層土壤中的氧氣和養(yǎng)分含量都極為有限,但在該區(qū)域仍生活著大量而多樣的微生物群落[5-6]。此外,深層土壤中的微生物大多具有鐵還原或硫酸鹽還原等獨特的代謝功能,在礦物風(fēng)化與土壤形成等過程中發(fā)揮重要作用[7]。闡明深層土壤中微生物群落特征及其影響因素,對深入探索陸地表層生態(tài)系統(tǒng)的生物多樣性以及生態(tài)系統(tǒng)服務(wù)的多功能性具有重要的生態(tài)學(xué)意義。

    環(huán)境異質(zhì)性顯著影響土壤微生物的組成與多樣性等群落特征。成土環(huán)境和成土過程通常會使土壤內(nèi)部產(chǎn)生差異化的微環(huán)境,并導(dǎo)致土體土壤理化屬性如質(zhì)地、孔隙以及礦質(zhì)養(yǎng)分的異質(zhì)性,這些差異將進一步影響土壤微生物群落的組成與分布。研究表明,即便是土壤環(huán)境中相隔僅幾毫米的微距離,土壤的理化屬性以及土壤微生物的相對豐度和組成等群落特征也大不相同[8-9]。此外,物種之間的相互作用也會影響土壤微生物群落特征。土壤中的微生物不能孤立生活,而總是與其他微生物發(fā)生著共生、寄生、捕食或是競爭等多種相互作用,從而構(gòu)成復(fù)雜的微生物生態(tài)網(wǎng)絡(luò),這些物種之間的相互作用關(guān)系會對群落的物種組成和結(jié)構(gòu)產(chǎn)生不同影響,并最終影響土壤微生物所行使的生態(tài)學(xué)功能。事實上,大多數(shù)生態(tài)學(xué)功能如養(yǎng)分循環(huán)或污染物降解等過程都無法由單一物種完成,而是多個物種之間相互作用關(guān)系的結(jié)果。因此,探究不同土壤環(huán)境中物種之間的相互作用關(guān)系,為更好地理解復(fù)雜生態(tài)學(xué)過程和生態(tài)系統(tǒng)的服務(wù)功能提供了科學(xué)依據(jù)[10]?;诟咄繙y序技術(shù),利用網(wǎng)絡(luò)分析的方法可以間接揭示微生物群落物種之間潛在的相互作用關(guān)系。目前,網(wǎng)絡(luò)分析方法已被成功應(yīng)用于挖掘土壤、湖泊、地下水和活性污泥等生態(tài)系統(tǒng)中微生物的相互作用模式研究[11-12]。

    廣泛分布于我國南方亞熱帶地區(qū)的網(wǎng)紋紅土是紅壤關(guān)鍵帶地下部的典型結(jié)構(gòu)之一,一般出現(xiàn)在均質(zhì)紅土剖面的下方,紅砂巖(基巖)的上方,是第四紀(jì)紅色風(fēng)化殼的重要組成部分。作為第四紀(jì)古氣候和古環(huán)境條件下的特殊產(chǎn)物,網(wǎng)紋紅土中包含著重要的環(huán)境變化信息,在全球環(huán)境變化的研究中具有極高的研究價值[13]。在歷史氣候與復(fù)雜的成土過程的綜合作用下,較為強烈的氧化還原和風(fēng)化淋溶過程不斷發(fā)生,土體中鐵的非均質(zhì)遷移和淋失以及脫硅富鋁化的持續(xù)進行,形成了在野外易于辨認(rèn)的紅白相間的網(wǎng)紋土體結(jié)構(gòu)[14]。網(wǎng)紋紅土的形成和發(fā)育過程記錄了環(huán)境的變化信息,對網(wǎng)紋紅土結(jié)構(gòu)的沉積學(xué)、環(huán)境學(xué)和生態(tài)學(xué)等特征進行解讀可以為第四紀(jì)環(huán)境的研究提供重要的科學(xué)依據(jù)。近年來,研究人員從不同角度對網(wǎng)紋紅土做了大量的科學(xué)研究,包括分析其物質(zhì)來源以解析網(wǎng)紋的成因、形成年代和古氣候指示意義等方面,但并未有對網(wǎng)紋紅土結(jié)構(gòu)的土壤生態(tài)學(xué)特別是土壤微生物學(xué)方面的相關(guān)研究。而由于風(fēng)化和成土過程導(dǎo)致的網(wǎng)紋層獨特的紅白土體很可能會形成獨特的微環(huán)境,并進一步形成差異化的生態(tài)位,從而影響其上微生物群落的組成以及微生物生態(tài)網(wǎng)絡(luò)。

    為揭示獨特的紅白相間的網(wǎng)紋土體結(jié)構(gòu)內(nèi)微生物群落組成、多樣性及生態(tài)網(wǎng)絡(luò)特征,本文以網(wǎng)紋紅土結(jié)構(gòu)體為研究對象,按照土體顏色將網(wǎng)紋結(jié)構(gòu)體拆分成紅色基質(zhì)和白色網(wǎng)紋兩部分,利用微生物16S rRNA高通量測序技術(shù)研究網(wǎng)紋紅土結(jié)構(gòu)體紅/白部分細菌群落的組成與多樣性等特征?;跍y序所獲取的數(shù)據(jù)集,利用CoNet(Co-occurrence network)網(wǎng)絡(luò)構(gòu)建方法挖掘紅/白部分細菌間的相互作用關(guān)系,并結(jié)合網(wǎng)紋土體的理化屬性探討網(wǎng)紋結(jié)構(gòu)體細菌群落特征的影響因素。

    1 材料與方法

    1.1 研究區(qū)概況與樣品采集

    研究區(qū)位于江西省鷹潭市的紅壤關(guān)鍵帶觀測平臺(red soil critical zone observatory, RSCZO)。紅壤關(guān)鍵帶觀測平臺設(shè)于鷹潭市余江縣劉家站墾殖場的孫家小流域(116°53′ E, 28°14′ N),流域總面積51.25 hm2,流域內(nèi)地勢平坦,坡度在3 ~ 5°,海拔在41 ~ 55 m,主要地貌類型為低山丘陵。研究區(qū)氣候類型為中亞熱帶季風(fēng)濕潤氣候,年平均溫度17.5 ℃,年平均降水量約1 795 mm?;诒狙芯康难芯磕康?,在流域內(nèi)選擇了一個典型樣點,該樣點的土地利用方式為水稻田。紅壤關(guān)鍵帶位置以及樣點位置如圖1所示。

    利用液壓旋轉(zhuǎn)鉆機取樣。該樣點網(wǎng)紋土體的深度范圍為1.4 ~ 5.0 m。將網(wǎng)紋樣品等間距(20 cm)分割,放入自封袋置于便攜式冰箱帶回實驗室。采樣時間為2016年4月,本研究共采集網(wǎng)紋樣品18段。在實驗室無菌的環(huán)境中,利用滅菌的塑料鑷子小心將成段網(wǎng)紋樣品的紅色基質(zhì)部分和白色網(wǎng)紋部分剝離開,分別收集紅色和白色部分,獲取樣品共36個。網(wǎng)紋紅土的出露位置、形態(tài)以及拆分后的網(wǎng)紋樣品如圖2所示。

    1.2 理化性質(zhì)測定

    網(wǎng)紋結(jié)構(gòu)體的紅白部分有機碳含量(SOC)用重鉻酸鉀–硫酸消化法測定;全氮(TN)含量用硒粉、硫酸銅、硫酸鉀消化蒸餾–凱氏定氮法測定;顆粒組成用激光粒度儀法測定;礦質(zhì)全量的測定選擇碳酸鋰–硼酸熔融–電感耦合等離子體質(zhì)譜法測定;游離鐵含量的測定用檸檬酸鈉–連二亞硫酸鈉提取,電感耦合等離子質(zhì)譜法測定。

    1.3 DNA提取與高通量測序

    選擇FastDNA Spin Kit for soil 提取試劑盒(MP Biomedicals, Santa Ana, CA, USA),按照操作指南提取網(wǎng)紋紅/白部分微生物總DNA。DNA提取后利用NanoDrop 2000 (Thermo Scientific, Wilmington, DE )在波長為260/280 nm處測定提取的DNA的濃度與質(zhì)量。

    圖1 研究區(qū)(紅壤關(guān)鍵帶)位置與樣點位置

    圖2 網(wǎng)紋紅土的位置、網(wǎng)紋的微形態(tài)和拆分后的網(wǎng)紋結(jié)構(gòu)

    對符合質(zhì)量要求的DNA目標(biāo)片段進行PCR擴增,擴增反應(yīng)的目標(biāo)片段為細菌的16S rRNA 的可變區(qū)V4 ~ V5區(qū)。PCR擴增反應(yīng)所用的引物是帶有barcode的通用引物515F (5’-GTGCCAGCMGCCG CGG-3’)和907R (5’-CCGTCAATTCMTTTRAGTTT- 3’)。PCR反應(yīng)體積共20 μl,其中包括10 ng DNA模板,2.5 μl 2.5 mmol/L的dNTP,兩端引物各0.8 μl,以及0.4 μl的Fast Pfu DNA聚合酶等。PCR反應(yīng)總共進行25個循環(huán),總的反應(yīng)流程為:95 ℃預(yù)變性反應(yīng)3 min,95 ℃變性30 s,55 ℃退火30 s,72 ℃延伸45 s,以及在 72 ℃延伸保溫10 min。

    測序平臺選擇Illumina MiSeq PE250 Platform,通過對原始序列過濾、序列拼接和去除嵌合體等,得到高質(zhì)量序列,并在97% 相似度水平上對操作分類單元(OTU)進行聚類與劃分,獲得基于序列的OTU。

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

    方差分析用SPSS 19.0進行,典范對應(yīng)分析(CCA)、主坐標(biāo)分析(PCoA)以及多元方差分析(PermANOVA)分析均使用R軟件(R version 4.0.0)的Vegan包進行。

    微生物生態(tài)網(wǎng)絡(luò)構(gòu)建方法參考Karoline Faust開發(fā)的基于相關(guān)/相似性的社會復(fù)雜網(wǎng)絡(luò)分析的CoNet網(wǎng)絡(luò)分析方法,在Cytoscape version 3.4.0平臺上運行[15-16]。在相似/相關(guān)性方法中,本研究選擇了Pearson correlation、Spearman correlation、Bray-Curtis dissimilarity和Kullback-Leibler dissimilarity 4種方法。獲取初始網(wǎng)絡(luò)后,利用Benjamini-Hochberg方法對獲得的相關(guān)系數(shù)進行標(biāo)準(zhǔn)化處理,并對相關(guān)關(guān)系的值進行校正,保留校正后<0.05的連線構(gòu)建最終的網(wǎng)絡(luò)。為統(tǒng)計網(wǎng)絡(luò)的拓撲學(xué)屬性差異,本文過濾并提取了每一網(wǎng)紋結(jié)構(gòu)體樣品對應(yīng)的子網(wǎng)絡(luò),利用Network Analyzer分析子網(wǎng)絡(luò)的拓撲特征。

    2 結(jié)果與討論

    2.1 網(wǎng)紋紅土結(jié)構(gòu)土壤理化屬性差異

    網(wǎng)紋土體結(jié)構(gòu)紅白部分的理化屬性如表1和表2所示。礦質(zhì)全量中鐵和錳在網(wǎng)紋紅色部分的含量顯著高于白色部分,硅含量在網(wǎng)紋的白色部分較高。網(wǎng)紋紅土結(jié)構(gòu)體紅白部分理化性質(zhì)的差異與網(wǎng)紋的形成過程密切相關(guān)。在干濕交替的情況下,氧化還原過程是網(wǎng)紋形成的主要成土過程之一。網(wǎng)紋層中的紅色部分傾向于氧化環(huán)境,含有較多高價態(tài)的鐵氧化物如赤鐵礦等,所以該部分土體呈現(xiàn)紅色;在網(wǎng)紋結(jié)構(gòu)體的白色部分,由于水分的作用,在還原的條件下,三價鐵礦物被還原成二價鐵并沿著裂隙不斷遷移,并伴隨著氧氣和養(yǎng)分的消耗,鉀、鈣、鈉、鎂等礦質(zhì)元素的淋失也更強烈。網(wǎng)紋結(jié)構(gòu)體理化屬性的分析也指示紅色基質(zhì)部分可能是一個礦質(zhì)養(yǎng)分相對豐富的微環(huán)境。

    表1 網(wǎng)紋紅土結(jié)構(gòu)紅、白部分理化屬性(g/kg)

    注:表中同列不同小寫字母代表網(wǎng)紋紅土結(jié)構(gòu)紅、白部分間的差異達顯著水平(<0.05);下同。

    表2 網(wǎng)紋紅土結(jié)構(gòu)紅、白部分土壤機械組成(V/V,%)

    2.2 網(wǎng)紋紅土結(jié)構(gòu)土壤細菌群落與組成

    細菌相對豐度分析結(jié)果表明,在網(wǎng)紋的紅色部分,放線菌門(Actinobacteria)、γ-變形菌綱(γ-Proteobacteria)、β-變形菌綱(β-Proteobacteria)、厚壁菌門(Firmicutes)等的相對豐度要顯著高于白色部分。而在網(wǎng)紋的白色部分,酸桿菌門(Acidobacteria)、綠彎菌門(Chloroflexi)、GAL15和硝化螺旋菌門(Nitrospirae)的相對豐度則較高(圖3)。PermANOVA結(jié)果表明,網(wǎng)紋的紅白部分細菌群落組成差異顯著(= 0.001)。

    酸桿菌門是一種典型的寡營養(yǎng)細菌,它們能夠利用幾丁質(zhì)等難分解的碳源作底物,在養(yǎng)分有限的情況下依然具有生存力[17]。此外,硝化螺旋菌門的相對豐度在白色部分也較高,該門類中含有大量參與硫酸鹽還原的寡營養(yǎng)厭氧細菌,因此更適合生活在礦質(zhì)養(yǎng)分含量較低的白色網(wǎng)紋部分[18]。而在紅色基質(zhì)中豐度較高的γ-變形菌綱(γ-Proteobacteria)是典型的富營養(yǎng)細菌。這可能是因為,網(wǎng)紋結(jié)構(gòu)體的白色部分由于長期的還原作用,元素的淋失和消耗都更強烈,為細菌提供了一個養(yǎng)分有限的還原性環(huán)境,因此酸桿菌門等寡營養(yǎng)細菌可以富集在白色網(wǎng)紋中。而對于紅色基質(zhì)部分,相對弱的淋失過程提供給細菌一個礦質(zhì)養(yǎng)分較豐富的生態(tài)環(huán)境,使放線菌門和γ-變形菌綱等細菌可以富集在網(wǎng)紋的紅色部分。

    (圖中*,**,***分別表示在P<0.05,P<0.01,P<0.001水平紅白部分間差異顯著)

    2.3 網(wǎng)紋紅土結(jié)構(gòu)土壤細菌生態(tài)網(wǎng)絡(luò)

    本研究分別構(gòu)建了網(wǎng)紋結(jié)構(gòu)體紅色基質(zhì)和白色網(wǎng)紋部分細菌間相互作用關(guān)系的生態(tài)網(wǎng)絡(luò),網(wǎng)絡(luò)構(gòu)建的結(jié)果如圖4所示。為統(tǒng)計分析紅、白部分細菌生態(tài)網(wǎng)絡(luò)的拓撲參數(shù)是否具有差異,本研究進一步提取了網(wǎng)紋結(jié)構(gòu)體中每個樣品對應(yīng)的子網(wǎng)絡(luò)。

    圖4 網(wǎng)紋紅土結(jié)構(gòu)紅、白部分細菌生態(tài)互作網(wǎng)絡(luò)

    利用Network Analyzer分析提取每個獨立樣品子網(wǎng)絡(luò)的拓撲學(xué)特征,結(jié)果表明,網(wǎng)紋的白色部分細菌生態(tài)互作網(wǎng)絡(luò)節(jié)點數(shù)目和連線數(shù)目都顯著高于紅色部分,細菌之間的負連線數(shù)目也顯著高于紅色部分。而網(wǎng)紋紅色部分的網(wǎng)絡(luò)聚類系數(shù)顯著高于白色部分(圖5)。

    生態(tài)網(wǎng)絡(luò)的拓撲學(xué)性質(zhì)可以在一定程度上反映網(wǎng)絡(luò)所代表的生態(tài)系統(tǒng)中物質(zhì)流、能量流以及信息流的流暢程度等。網(wǎng)絡(luò)的聚類系數(shù)代表網(wǎng)絡(luò)中一個節(jié)點和其他節(jié)點之間的連通程度,聚類系數(shù)越高表示群落對環(huán)境變化的敏感性越強,當(dāng)外界環(huán)境發(fā)生變化時,環(huán)境因子的擾動通常會在較短的時間內(nèi)影響整個網(wǎng)絡(luò)[19]。研究表明,網(wǎng)紋結(jié)構(gòu)體的紅色基質(zhì)部分細菌生態(tài)網(wǎng)絡(luò)的聚類系數(shù)較高,說明該部分細菌之間相互作用關(guān)系更緊密。而在白色網(wǎng)紋部分,較低的聚類系數(shù)則代表細菌生態(tài)網(wǎng)絡(luò)結(jié)構(gòu)比較松散,網(wǎng)絡(luò)內(nèi)細菌之間的相關(guān)關(guān)系對比于紅色部分相對較弱。

    2.4 網(wǎng)紋紅土結(jié)構(gòu)土壤細菌群落的影響因素

    典范對應(yīng)分析結(jié)果(圖6)表明,網(wǎng)紋紅土結(jié)構(gòu)紅、白部分細菌群落的差異是由于紅、白部分的微環(huán)境異質(zhì)性決定的。受第四紀(jì)古氣候的影響,地下水的上下波動和季節(jié)性的干濕交替造成礦質(zhì)元素包括鉀、鈣和鈉等元素的強烈淋溶和鐵、鋁等元素不斷積累。由于網(wǎng)紋結(jié)構(gòu)體一般出現(xiàn)在均質(zhì)紅土剖面的下方,所以能接收到的地表輸入的養(yǎng)分極為有限,因此,礦物所含有的礦質(zhì)養(yǎng)分是微生物生存的生命資源,是顯著影響微生物群落組成的重要因素[20-22]。在本研究中,礦質(zhì)全量中的鐵和錳含量與網(wǎng)紋結(jié)構(gòu)體紅色基質(zhì)部分的細菌群落顯著相關(guān)。鐵是地殼中含量第四高的元素,亞鐵和三價鐵之間的氧化還原過程是紅壤關(guān)鍵帶生態(tài)系統(tǒng)生物地球化學(xué)循環(huán)中的重要過程之一,同時也深刻地影響著紅壤地區(qū)土壤的發(fā)生與演變過程。而土壤和沉積物中亞鐵和三價鐵之間的氧化還原過程常常由土壤細菌主導(dǎo),白色網(wǎng)紋中發(fā)現(xiàn)的、ra和等細菌都是典型的三價鐵還原細菌[23-24]。此外,網(wǎng)紋結(jié)構(gòu)體的機械組成也顯著影響紅、白部分細菌群落的分布。機械組成代表的質(zhì)地狀況是土壤的重要物理屬性之一,同時也代表了土體的水文狀況[25]。而在深部的網(wǎng)紋層土壤中,土體的水文狀況不只與網(wǎng)紋的形成和發(fā)育有關(guān),也決定了養(yǎng)分狀況,因而對細菌群落有顯著影響[26]。

    圖5 網(wǎng)紋紅土結(jié)構(gòu)紅、白部分細菌生態(tài)網(wǎng)絡(luò)的拓撲學(xué)性質(zhì)

    圖6 網(wǎng)紋紅土結(jié)構(gòu)土壤理化屬性與細菌群落結(jié)構(gòu)的典范對應(yīng)分析

    3 結(jié)論

    網(wǎng)紋紅土結(jié)構(gòu)是第四紀(jì)紅色風(fēng)化殼的重要組成部分,對研究第四紀(jì)古氣候和古環(huán)境變遷至關(guān)重要。網(wǎng)紋結(jié)構(gòu)體內(nèi)部紅、白部分細菌群落組成具有顯著差異,這些差異主要是由網(wǎng)紋土體中紅、白部分的環(huán)境條件差異決定。本研究探索性地將網(wǎng)紋結(jié)構(gòu)進行拆分,分別對不同部分的細菌群落特征進行分析,研究結(jié)果對認(rèn)識紅壤關(guān)鍵帶深部生態(tài)系統(tǒng)的生物多樣性和多功能性具有重要的生態(tài)學(xué)意義。

    [1] National Research Council. Basic research opportunities in earth sciences[M]. Washington, DC: National Academies Press, 2001.

    [2] Williams J W, Jackson S T, Kutzbach J E. Projected distributions of novel and disappearing climates by 2100 AD[J]. Proceedings of the National Academy of Sciences of the United States of America, 2007, 104(14): 5738–5742.

    [3] Wilkinson M T, Richards P J, Humphreys G S. Breaking ground: Pedological, geological, and ecological implications of soil bioturbation[J]. Earth-Science Reviews, 2009, 97(1/2/3/4): 257–272.

    [4] Fierer N, Schimel J P, Holden P A. Variations in microbial community composition through two soil depth profiles[J]. Soil Biology and Biochemistry, 2003, 35(1): 167–176.

    [5] Amend J P, Teske A. Expanding frontiers in deep subsurface microbiology[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2005, 219(1/2): 131–155.

    [6] Whitman W B, Coleman D C, Wiebe W J. Prokaryotes: The unseen majority[J]. PNAS, 1998, 95(12): 6578–6583.

    [7] Gleeson D, Mathes F, Farrell M, et al. Environmental drivers of soil microbial community structure and function at the Avon River Critical Zone Observatory[J]. Science of the Total Environment, 2016, 571: 1407–1418.

    [8] 杜雄峰, 厲舒禎, 馮凱, 等. 農(nóng)牧交錯帶草地土壤剖面微生物總量、多樣性和互作網(wǎng)絡(luò)的垂直分布特征[J]. 微生物學(xué)通報, 2020, 47(9): 2789–2806.

    [9] Brockman F J, Murray C J. Subsurface microbiological heterogeneity: Current knowledge, descriptive approaches and applications[J]. FEMS Microbiology Reviews, 1997, 20(3/4): 231–247.

    [10] Xiao X, Liang Y T, Zhou S, et al. Fungal community reveals less dispersal limitation and potentially more connected network than that of bacteria in bamboo forest soils[J]. Molecular Ecology, 2018, 27(2): 550–563.

    [11] Barberán A, Bates S T, Casamayor E O, et al. Using network analysis to explore co-occurrence patterns in soil microbial communities[J]. The ISME Journal, 2012, 6(2): 343–351.

    [12] Ju F, Zhang T. Bacterial assembly and temporal dynamics in activated sludge of a full-scale municipal wastewater treatment plant[J]. The ISME Journal, 2015, 9(3): 683–695.

    [13] Liu C C, Deng C L, Liu Q S. Mineral magnetic studies of the vermiculated red soils in southeast China and their paleoclimatic significance[J]. Paleogeography, Palaeoclimatology, Paleoecology, 2012, 329/330: 173–183.

    [14] Gu Y S, Huang X Y, Zhang W G, et al. Red paleosols development in response to the enhanced east Asia summer monsoon since the mid-Pleistocene in South China: Evidence derived from magnetic properties and molecular fossil records[J]. Journal of Earth Science, 2013, 24(3): 382–396.

    [15] Faust K, Raes J. CoNet app: Inference of biological association networks using Cytoscape[J]. F1000Research, 2016, 5: 1519.

    [16] Smoot M E, Ono K, Ruscheinski J, et al. Cytoscape 2.8: New features for data integration and network visualization[J]. Bioinformatics, 2011, 27(3): 431–432.

    [17] Tsitko I, Lusa M, Lehto J, et al. The variation of microbial communities in a depth profile of an acidic, nutrient-poor boreal bog in southwestern Finland[J]. Open Journal of Ecology, 2014, 4(13): 832–859.

    [18] Attard E, Poly F, Commeaux C, et al. Shifts between- and-like nitrite oxidizers underlie the response of soil potential nitrite oxidation to changes in tillage practices[J]. Environmental Microbiology, 2010, 12(2): 315–326.

    [19] Chow C E T, Kim D Y, Sachdeva R, et al. Top-down controls on bacterial community structure: Microbial network analysis of bacteria, T4-like viruses and protists[J]. The ISME Journal, 2014, 8(4): 816–829.

    [20] Uroz S, Kelly L C, Turpault M P, et al. The mineralosphere concept: Mineralogical control of the distribution and function of mineral-associated bacterial communities[J]. Trends in Microbiology, 2015, 23(12): 751–762.

    [21] Roberts J A. Inhibition and enhancement of microbial surface colonization: The role of silicate composition[J]. Chemical Geology, 2004, 212(3/4): 313–327.

    [22] Carson J K, Gonzalez-Quin?ones V, Murphy D V, et al. Low pore connectivity increases bacterial diversity in soil[J]. Applied and Environmental Microbiology, 2010, 76(12): 3936–3942.

    [23] Lovley D R. Dissimilatory Fe(III) and Mn(IV) reduction[J]. Microbiological Reviews, 1991, 55(2): 259–287.

    [24] Weber K A, Achenbach L A, Coates J D. Microorganisms pumping iron: Anaerobic microbial iron oxidation and reduction[J]. Nature Reviews Microbiology, 2006, 4(10): 752–764.

    [25] Wallis M G, Scotter D R, Horne D J. An evaluation of the intrinsic sorptivity water repellency index on a range of New Zealand soils[J]. Soil Research, 1991, 29(3): 353.

    [26] Girvan M S, Bullimore J, Pretty J N, et al. Soil type is the primary determinant of the composition of the total and active bacterial communities in arable soils[J]. Applied and Environmental Microbiology, 2003, 69(3): 1800–1809.

    Characteristics of Bacterial Communities and Their Co-occurrence Networks in a Typical Plinthic Horizon

    ZHAO Xiaorui1,2, WU Huayong1,2, ZHANG Ganlin1,2,3*

    (1 State Key Laboratory of Soil and Sustainable Agriculture,Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China; 2 University of Chinese Academy of Sciences, Beijing 100049, China; 3 Key Laboratory of Watershed Geographic Sciences, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing 210008, China)

    Plinthite horizon (net-like soil horizon), as a part of the Quaternary red palesols, is believed to be an important record of the Quaternary paleoclimate and paleo-environment changes of subtropical areas in China. As an important part of the subsurface structure of red soil critical zone, plinthite horizon can harbor various bacterial communities, which may play key roles in mineral weathering and soil formation processes. Previous studies usually took the entire soil layer or genetic horizon as a single sample, while the intra-horizon differentiation of bacterial communities and co-occurrence network remain largely unknown. It is hypothesed that the plinthite horizon may affect the distribution of bacterial communities in the red matrix and white veins. Clarifying the characteristics of soil bacterial communities in plinthite horizon and its impact factors would help for exploring the biodiversity and understanding the multi-functionality of red soil ecosystems. By using the high throughput sequencing processes, the plinthite structure were here divided into two parts and the bacterial communities and their ecological networks under red matrix and white veins were analyzed. Results demonstrate that there is difference of bacterial composition between red matrix and white veins. Relative abundance of the phylum of Actinobacteria and Bacteroidetes are significantly higher in the red matrix of the plinthite horizon, while the relative abundance of Acidobacteria and Nitrospirae are significantly higher in the white veins. Bacterial ecological network and its topological features are also different between the red matrix and white veins in the plinthite horizon. Co-occurrence network (CoNet) analysis demonstrate that the ecological network of white veins contains more edges and has a lower network’s clustering coefficient value. CCA analysis demonstrate that bacterial community composition and biodiversity in the plinthite horizon are caused by the environmental heterogeneity between the red matrix and white veins. Major element contents and particle size distribution are mainly impact factors that shape bacterial communities in the red matrix of the plinthite horizon.

    Red soil critical zone; Plinthite horizon; Bacterial community; Co-occurrence network; Impact factors

    S154.36

    A

    10.13758/j.cnki.tr.2022.05.015

    趙曉蕊, 吳華勇, 張甘霖. 典型網(wǎng)紋紅土結(jié)構(gòu)土壤細菌群落及其生態(tài)網(wǎng)絡(luò)特征. 土壤, 2022, 54(5): 986–992.

    國家自然科學(xué)基金項目(41571130051,41771251)和中國科學(xué)院科技服務(wù)網(wǎng)絡(luò)計劃項目(KFJ-STS-ZDTP-039)資助。

    (glzhang@issas.ac.cn)

    趙曉蕊(1986—),女,遼寧錦州人,博士研究生,主要從事紅壤關(guān)鍵帶土壤微生物生態(tài)研究。E-mail: xrzhao@issas.ac.cn

    猜你喜歡
    網(wǎng)紋紅土群落
    網(wǎng)紋甜瓜增加網(wǎng)紋密度的方法
    大學(xué)生牙齦炎齦上菌斑的微生物群落
    合成微生物群落在發(fā)酵食品中的應(yīng)用研究
    種植流翔高鈣網(wǎng)紋瓜 山東農(nóng)民唱響增收曲
    朽木頌
    老友(2020年2期)2020-03-08 14:25:40
    失眠記
    詩潮(2017年4期)2017-12-05 10:16:18
    紅土象湖感化溪
    紅土地(2016年6期)2017-01-15 13:45:58
    如何提高柔印中網(wǎng)紋輥的清洗效率
    今日印刷(2016年2期)2016-03-17 08:25:45
    紅土鎳礦中的鈣、鎂、鐵、鉆、鎳等含量測定
    淺談網(wǎng)紋輥在色彩穩(wěn)定中的重要作用
    免费看日本二区| 亚洲成人久久性| 天天躁日日操中文字幕| 久久99精品国语久久久| av在线亚洲专区| 51国产日韩欧美| 亚洲四区av| 一级二级三级毛片免费看| 99热网站在线观看| 人体艺术视频欧美日本| 国产伦在线观看视频一区| 国产精品久久电影中文字幕| 99riav亚洲国产免费| 综合色av麻豆| 亚洲激情五月婷婷啪啪| 高清在线视频一区二区三区 | 天堂影院成人在线观看| 亚洲精品国产av成人精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 少妇人妻精品综合一区二区 | av卡一久久| 国产成年人精品一区二区| 熟女电影av网| 真实男女啪啪啪动态图| 麻豆乱淫一区二区| 舔av片在线| 日韩欧美三级三区| 一个人免费在线观看电影| 久久久精品欧美日韩精品| 毛片女人毛片| 少妇裸体淫交视频免费看高清| 欧美xxxx黑人xx丫x性爽| 免费不卡的大黄色大毛片视频在线观看 | 一区二区三区免费毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 青青草视频在线视频观看| 久久亚洲精品不卡| 波野结衣二区三区在线| 午夜久久久久精精品| 男插女下体视频免费在线播放| 国产成人精品婷婷| 欧美人与善性xxx| 熟女人妻精品中文字幕| 黄色视频,在线免费观看| 久久精品久久久久久噜噜老黄 | 在线国产一区二区在线| 日韩 亚洲 欧美在线| 天天躁日日操中文字幕| 国产蜜桃级精品一区二区三区| 三级国产精品欧美在线观看| 亚洲性久久影院| 又爽又黄无遮挡网站| 国产久久久一区二区三区| 成人二区视频| 亚洲不卡免费看| 日本一本二区三区精品| 如何舔出高潮| 日韩高清综合在线| 国产伦精品一区二区三区视频9| 黄色日韩在线| 国产乱人偷精品视频| 国产在线男女| www.av在线官网国产| 亚洲成a人片在线一区二区| 中文欧美无线码| 国产大屁股一区二区在线视频| 日韩成人av中文字幕在线观看| 色噜噜av男人的天堂激情| 欧美高清成人免费视频www| 桃色一区二区三区在线观看| 搡女人真爽免费视频火全软件| 成年版毛片免费区| 一卡2卡三卡四卡精品乱码亚洲| 中国美女看黄片| 99在线人妻在线中文字幕| 精品不卡国产一区二区三区| 联通29元200g的流量卡| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区三区四区久久| 久久热精品热| 日日摸夜夜添夜夜添av毛片| 成人性生交大片免费视频hd| 五月玫瑰六月丁香| 亚洲乱码一区二区免费版| 国产精品.久久久| 国产精品麻豆人妻色哟哟久久 | 精品人妻视频免费看| 国产蜜桃级精品一区二区三区| 69av精品久久久久久| 最近2019中文字幕mv第一页| 成人美女网站在线观看视频| 国产精品一区二区三区四区免费观看| 色哟哟哟哟哟哟| 精品久久久久久久久亚洲| 97在线视频观看| 在现免费观看毛片| 日本黄色片子视频| 赤兔流量卡办理| 久久精品91蜜桃| 欧美日韩精品成人综合77777| 婷婷色综合大香蕉| 国国产精品蜜臀av免费| 最近中文字幕高清免费大全6| 可以在线观看的亚洲视频| 久久久久久久久久久丰满| 亚洲,欧美,日韩| 亚洲av中文av极速乱| 欧美精品国产亚洲| 一个人免费在线观看电影| 在现免费观看毛片| 日本撒尿小便嘘嘘汇集6| 亚洲图色成人| 天堂网av新在线| 精品久久久噜噜| 亚洲av电影不卡..在线观看| 久久精品国产99精品国产亚洲性色| 成年免费大片在线观看| 日本黄色片子视频| 欧美成人精品欧美一级黄| 日本黄色片子视频| 老司机福利观看| 狂野欧美白嫩少妇大欣赏| 国产成人aa在线观看| 熟女人妻精品中文字幕| 一本久久中文字幕| 亚洲国产精品成人久久小说 | av国产免费在线观看| 最新中文字幕久久久久| 国产黄a三级三级三级人| 成人特级av手机在线观看| 免费大片18禁| 一个人免费在线观看电影| 一本久久精品| 日韩亚洲欧美综合| 国产精品女同一区二区软件| 好男人视频免费观看在线| 精品久久久久久久末码| 给我免费播放毛片高清在线观看| 国产成人aa在线观看| 免费在线观看成人毛片| 亚洲av免费在线观看| 欧美三级亚洲精品| 亚洲乱码一区二区免费版| 久久久久久大精品| 国产黄色小视频在线观看| 久久国内精品自在自线图片| 色视频www国产| 国产高清三级在线| 老熟妇乱子伦视频在线观看| 国产亚洲5aaaaa淫片| 综合色av麻豆| 国产老妇伦熟女老妇高清| 最近最新中文字幕大全电影3| 免费av不卡在线播放| 成人三级黄色视频| 偷拍熟女少妇极品色| 久久精品国产亚洲网站| 免费人成在线观看视频色| 美女脱内裤让男人舔精品视频 | 搡老妇女老女人老熟妇| 色视频www国产| 尤物成人国产欧美一区二区三区| 亚洲最大成人中文| 亚洲av中文av极速乱| АⅤ资源中文在线天堂| 综合色av麻豆| 亚洲欧美日韩高清专用| 亚洲精品乱码久久久v下载方式| av在线亚洲专区| 欧美+亚洲+日韩+国产| 久久亚洲精品不卡| 观看美女的网站| 亚洲七黄色美女视频| 国内精品宾馆在线| 中文字幕人妻熟人妻熟丝袜美| 国产免费男女视频| 此物有八面人人有两片| 丰满人妻一区二区三区视频av| 女人十人毛片免费观看3o分钟| 久久人人爽人人片av| 国产探花在线观看一区二区| 国产熟女欧美一区二区| 在线观看一区二区三区| 在线观看美女被高潮喷水网站| 一级黄片播放器| 欧美变态另类bdsm刘玥| 亚洲一区二区三区色噜噜| 亚洲av男天堂| 亚洲精品影视一区二区三区av| 午夜福利成人在线免费观看| 在线国产一区二区在线| 男的添女的下面高潮视频| www.色视频.com| 91在线精品国自产拍蜜月| 少妇人妻精品综合一区二区 | 高清午夜精品一区二区三区 | 免费黄网站久久成人精品| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久久久丰满| 九九爱精品视频在线观看| 久久久国产成人精品二区| 久久精品久久久久久久性| 三级毛片av免费| 能在线免费观看的黄片| 国产精华一区二区三区| 午夜爱爱视频在线播放| 国产中年淑女户外野战色| 一进一出抽搐动态| 国产成人aa在线观看| 精品无人区乱码1区二区| 麻豆精品久久久久久蜜桃| 国产亚洲欧美98| 简卡轻食公司| 成人性生交大片免费视频hd| 午夜福利视频1000在线观看| 免费人成视频x8x8入口观看| 久久精品综合一区二区三区| 亚洲七黄色美女视频| 国产老妇伦熟女老妇高清| 免费看日本二区| 国产色爽女视频免费观看| 欧美一级a爱片免费观看看| 人人妻人人看人人澡| 欧美最新免费一区二区三区| 成年女人看的毛片在线观看| 综合色av麻豆| 免费搜索国产男女视频| 99久国产av精品国产电影| 国产一级毛片在线| 又黄又爽又刺激的免费视频.| 国产伦在线观看视频一区| 国产一级毛片七仙女欲春2| 国产三级在线视频| 在线免费观看的www视频| 免费看美女性在线毛片视频| 亚洲精品影视一区二区三区av| 日韩 亚洲 欧美在线| 亚洲美女搞黄在线观看| 午夜免费男女啪啪视频观看| 婷婷精品国产亚洲av| 国产精品久久久久久精品电影| 丰满人妻一区二区三区视频av| 给我免费播放毛片高清在线观看| 成人鲁丝片一二三区免费| 九草在线视频观看| 精品人妻偷拍中文字幕| 偷拍熟女少妇极品色| 国产精品蜜桃在线观看 | 欧美高清成人免费视频www| 在线a可以看的网站| 日韩av不卡免费在线播放| 五月伊人婷婷丁香| 淫秽高清视频在线观看| 国产成人影院久久av| 国产精品久久久久久久久免| 校园春色视频在线观看| 日韩视频在线欧美| 欧美最黄视频在线播放免费| 又爽又黄a免费视频| 精品日产1卡2卡| 日本撒尿小便嘘嘘汇集6| 在线免费十八禁| 久久久久网色| 最后的刺客免费高清国语| 校园人妻丝袜中文字幕| 亚洲av成人av| 日韩精品有码人妻一区| 亚洲激情五月婷婷啪啪| 国产亚洲精品久久久久久毛片| 国产探花极品一区二区| 赤兔流量卡办理| 联通29元200g的流量卡| 九九爱精品视频在线观看| 免费看日本二区| 韩国av在线不卡| 久99久视频精品免费| 午夜精品在线福利| 亚洲精品影视一区二区三区av| 日本免费a在线| 国产黄a三级三级三级人| 国产蜜桃级精品一区二区三区| 亚洲无线在线观看| 亚洲va在线va天堂va国产| 成人二区视频| 美女被艹到高潮喷水动态| 你懂的网址亚洲精品在线观看 | 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 国产黄片美女视频| 又粗又硬又长又爽又黄的视频 | 在线播放无遮挡| 男人舔女人下体高潮全视频| 久久午夜亚洲精品久久| 午夜福利成人在线免费观看| 日韩在线高清观看一区二区三区| 国产免费一级a男人的天堂| 色综合亚洲欧美另类图片| 日日干狠狠操夜夜爽| 韩国av在线不卡| 日韩大尺度精品在线看网址| 日本成人三级电影网站| 日本与韩国留学比较| 亚州av有码| 亚洲五月天丁香| 久久久国产成人精品二区| 天堂av国产一区二区熟女人妻| 黄色配什么色好看| 久久99精品国语久久久| 欧美精品国产亚洲| 国产成人aa在线观看| 亚洲三级黄色毛片| 国产精品1区2区在线观看.| 亚洲欧美精品专区久久| 久久精品综合一区二区三区| 麻豆成人午夜福利视频| 久久久a久久爽久久v久久| h日本视频在线播放| avwww免费| 亚洲七黄色美女视频| 一级av片app| 免费av观看视频| 三级国产精品欧美在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美激情国产日韩精品一区| 18禁黄网站禁片免费观看直播| 少妇裸体淫交视频免费看高清| 最近的中文字幕免费完整| 久久久久久久亚洲中文字幕| 色哟哟哟哟哟哟| 久久婷婷人人爽人人干人人爱| 六月丁香七月| 久久精品国产鲁丝片午夜精品| 久久精品91蜜桃| 床上黄色一级片| 亚洲国产精品合色在线| 亚洲无线观看免费| 男女啪啪激烈高潮av片| 久久综合国产亚洲精品| 成年免费大片在线观看| 三级男女做爰猛烈吃奶摸视频| 看非洲黑人一级黄片| 精品免费久久久久久久清纯| 中出人妻视频一区二区| 亚洲,欧美,日韩| 亚洲美女视频黄频| 六月丁香七月| av黄色大香蕉| 久久久久久久久久久免费av| 亚洲成人中文字幕在线播放| 日本爱情动作片www.在线观看| 蜜桃亚洲精品一区二区三区| kizo精华| 国产不卡一卡二| 久久久久久久久大av| 最近2019中文字幕mv第一页| 日韩大尺度精品在线看网址| 色5月婷婷丁香| 久久精品夜色国产| 观看免费一级毛片| 久久鲁丝午夜福利片| 亚洲av中文av极速乱| 国国产精品蜜臀av免费| 欧美色视频一区免费| av专区在线播放| 欧美+亚洲+日韩+国产| 中国国产av一级| 国产单亲对白刺激| 亚洲第一区二区三区不卡| 日本免费a在线| av在线老鸭窝| 国产一区二区三区av在线 | 欧美最新免费一区二区三区| 国产伦精品一区二区三区视频9| 日本撒尿小便嘘嘘汇集6| 国产老妇女一区| 久久久精品94久久精品| 亚洲人成网站在线播| 身体一侧抽搐| 日本爱情动作片www.在线观看| 国内久久婷婷六月综合欲色啪| 亚洲av中文av极速乱| 哪里可以看免费的av片| 少妇人妻精品综合一区二区 | 国产免费一级a男人的天堂| 亚洲内射少妇av| 亚洲成a人片在线一区二区| 久久久久久久久久成人| 国产精品久久久久久亚洲av鲁大| 我要看日韩黄色一级片| 久久精品国产99精品国产亚洲性色| 亚洲欧美日韩东京热| 久久精品夜夜夜夜夜久久蜜豆| 国产成人影院久久av| av专区在线播放| 2021天堂中文幕一二区在线观| 在线观看美女被高潮喷水网站| 国产精品久久电影中文字幕| 亚洲中文字幕日韩| 国产一区二区亚洲精品在线观看| 99热这里只有是精品在线观看| 国产免费男女视频| 99久久精品热视频| 男人和女人高潮做爰伦理| 中文字幕av在线有码专区| 在线观看一区二区三区| 最近手机中文字幕大全| 99热精品在线国产| 欧美日本亚洲视频在线播放| 99久久无色码亚洲精品果冻| 久久久久久久亚洲中文字幕| 男人舔奶头视频| а√天堂www在线а√下载| 内地一区二区视频在线| 寂寞人妻少妇视频99o| 哪里可以看免费的av片| 麻豆av噜噜一区二区三区| 国产精品不卡视频一区二区| 国产精品一及| 久久精品夜色国产| 国产毛片a区久久久久| 天堂√8在线中文| 美女高潮的动态| 小说图片视频综合网站| 男女边吃奶边做爰视频| 99热6这里只有精品| 特级一级黄色大片| 日韩三级伦理在线观看| 国产麻豆成人av免费视频| 亚洲国产精品合色在线| 亚洲精品色激情综合| 午夜激情福利司机影院| 欧美高清成人免费视频www| 少妇丰满av| 精品一区二区三区视频在线| 亚洲欧美日韩东京热| 黑人高潮一二区| 欧美极品一区二区三区四区| 国产免费一级a男人的天堂| 午夜亚洲福利在线播放| 一本精品99久久精品77| 欧美日韩在线观看h| 最后的刺客免费高清国语| 国产成人午夜福利电影在线观看| 国产蜜桃级精品一区二区三区| 久久人人爽人人片av| 欧美三级亚洲精品| 美女xxoo啪啪120秒动态图| 国产单亲对白刺激| 91久久精品电影网| 亚洲欧美日韩卡通动漫| 美女被艹到高潮喷水动态| 在线观看一区二区三区| 最近的中文字幕免费完整| 国产精品一区二区在线观看99 | 亚洲最大成人手机在线| 成人永久免费在线观看视频| www.色视频.com| 少妇人妻精品综合一区二区 | 午夜精品国产一区二区电影 | 一边摸一边抽搐一进一小说| 亚洲人成网站高清观看| 精华霜和精华液先用哪个| 一本一本综合久久| 伊人久久精品亚洲午夜| 国产精品人妻久久久影院| 久久久国产成人精品二区| 一进一出抽搐gif免费好疼| 91久久精品国产一区二区成人| 亚洲丝袜综合中文字幕| 草草在线视频免费看| 99热全是精品| 国内精品美女久久久久久| 亚洲精品影视一区二区三区av| 久久鲁丝午夜福利片| 3wmmmm亚洲av在线观看| 人妻少妇偷人精品九色| 亚洲欧美日韩东京热| АⅤ资源中文在线天堂| 亚洲精品乱码久久久v下载方式| 九九在线视频观看精品| 久久综合国产亚洲精品| 国产色婷婷99| 亚洲真实伦在线观看| 深爱激情五月婷婷| 久久久精品大字幕| 国产精品av视频在线免费观看| 成人鲁丝片一二三区免费| av在线蜜桃| 好男人视频免费观看在线| 亚洲欧美日韩东京热| 国产一区二区三区在线臀色熟女| 欧美一级a爱片免费观看看| 免费av毛片视频| 国产精品久久久久久久久免| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久末码| 尾随美女入室| 有码 亚洲区| 国产在视频线在精品| 美女大奶头视频| 51国产日韩欧美| 成人特级黄色片久久久久久久| 亚洲av第一区精品v没综合| 久久久欧美国产精品| 亚洲自偷自拍三级| 亚洲第一电影网av| 国产精品,欧美在线| 亚洲av电影不卡..在线观看| 久久婷婷人人爽人人干人人爱| 日韩欧美国产在线观看| 国产高清视频在线观看网站| 一级毛片aaaaaa免费看小| 九色成人免费人妻av| 亚洲不卡免费看| 日韩人妻高清精品专区| 亚洲欧美清纯卡通| 日韩,欧美,国产一区二区三区 | 成人国产麻豆网| 精品无人区乱码1区二区| 国产精品1区2区在线观看.| 国产精品女同一区二区软件| 男女那种视频在线观看| 亚洲成人精品中文字幕电影| 一进一出抽搐动态| 亚洲欧美精品自产自拍| 男女视频在线观看网站免费| 色5月婷婷丁香| 五月玫瑰六月丁香| 国产精品永久免费网站| 2021天堂中文幕一二区在线观| 一本久久精品| 亚洲经典国产精华液单| 舔av片在线| 国产伦精品一区二区三区四那| www.色视频.com| 精品无人区乱码1区二区| 国产成人精品久久久久久| 欧美成人a在线观看| 搞女人的毛片| 国产精品日韩av在线免费观看| 直男gayav资源| 亚洲国产精品sss在线观看| 国产成人aa在线观看| 天堂中文最新版在线下载 | 亚洲精品乱码久久久久久按摩| 日韩高清综合在线| 少妇高潮的动态图| 欧美一区二区亚洲| 小蜜桃在线观看免费完整版高清| 国产精品人妻久久久影院| 国产欧美日韩精品一区二区| 中国国产av一级| 久久久久久久久久久免费av| 一区二区三区免费毛片| 看免费成人av毛片| 欧美性感艳星| 少妇的逼水好多| 床上黄色一级片| 亚洲在线观看片| 亚洲欧美成人综合另类久久久 | 美女 人体艺术 gogo| 国产精品一区二区性色av| 午夜老司机福利剧场| 99久久九九国产精品国产免费| 欧美日韩乱码在线| 成人午夜精彩视频在线观看| 波多野结衣高清作品| 国产 一区 欧美 日韩| 婷婷六月久久综合丁香| 中文字幕免费在线视频6| 亚洲国产精品sss在线观看| 日韩一区二区三区影片| 国产极品天堂在线| 久99久视频精品免费| 日本-黄色视频高清免费观看| 免费搜索国产男女视频| 国产精品乱码一区二三区的特点| 国产视频首页在线观看| 免费av观看视频| 91久久精品国产一区二区三区| 联通29元200g的流量卡| 又粗又硬又长又爽又黄的视频 | 大香蕉久久网| 日日干狠狠操夜夜爽| 国产探花极品一区二区| 久久久精品94久久精品| 久久这里有精品视频免费| 99久久精品热视频| 久久精品国产亚洲av天美| 亚洲精品久久久久久婷婷小说 | 久久99热这里只有精品18| 国产精品久久电影中文字幕| 男人舔女人下体高潮全视频| 中文字幕制服av| 看非洲黑人一级黄片| 亚洲精品456在线播放app| 18禁黄网站禁片免费观看直播| 国产单亲对白刺激| 日韩欧美精品v在线| 成人二区视频| 国产精品福利在线免费观看| 夜夜夜夜夜久久久久| 欧美日韩在线观看h| 国产精品伦人一区二区| 精华霜和精华液先用哪个| 极品教师在线视频| 日韩中字成人| 国产单亲对白刺激| 日本色播在线视频| 晚上一个人看的免费电影| 精品久久久久久久久亚洲| 如何舔出高潮| 99热只有精品国产| 欧美人与善性xxx| 一级毛片我不卡|