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

    沿黃稻區(qū)土壤菌群結(jié)構(gòu)與生態(tài)功能研究

    2024-05-15 21:48:30田瑩瑩郭琦陳星儒孫向輝趙京

    田瑩瑩 郭琦 陳星儒 孫向輝 趙京

    摘? 要:為了解沿黃稻區(qū)不同基質(zhì)土壤的微生物組成,采集了河灘砂質(zhì)土壤(sand soil,SS)和毗鄰稻田土壤(paddy soil,PS).利用16S rRNA基因高通量測(cè)序技術(shù)揭示沿黃稻區(qū)土壤菌群的組成及其潛在的生態(tài)功能.Alpha多樣性分析顯示兩種土壤的菌群豐富度和多樣性指數(shù)無顯著性差異.進(jìn)一步的物種組成分析,發(fā)現(xiàn)unclassified_f? Micrococcaceae、鞘氨醇單胞菌屬(Sphingomonas)、norank_f? noorank_o? Vicinamibacterales、norank_f? Vicinamibacterales在SS和PS中廣布,與土壤基質(zhì)類型無關(guān)豐度,但芽孢桿菌屬(Bacilllus)、硫桿菌屬(Thiobacillus)在PS中的相對(duì)豐度顯著高于SS.基于FAPROTAX數(shù)據(jù)庫(kù)的功能熱圖分析顯示PS(尤其是P2_1、P3_1)中芳香烴降解功能序列的相對(duì)豐度高于SS,結(jié)合菌群相對(duì)豐度結(jié)果推測(cè)類諾卡氏菌屬(Nocardioides)、unclassified_f_Nocardioidaceae是該差異產(chǎn)生的關(guān)鍵.SS中光養(yǎng)微生物豐度顯著高于PS,尤以卟啉桿菌(Porphyrobacter)、unclassified_f_Rhodobacteraceae、紅細(xì)菌(Rhodobacter)為主的不產(chǎn)氧光合細(xì)菌的豐度較高,其中紅細(xì)菌科(Rhodobacteraceae)是參與SS生物膜形成的重要微生物.研究表明土壤類型對(duì)沿黃稻區(qū)菌群組成有一定的影響,并驅(qū)動(dòng)了菌群生態(tài)功能分異的形成.

    關(guān)鍵詞:壤菌群結(jié)構(gòu);生態(tài)功能;稻田土壤;砂質(zhì)土壤

    中圖分類號(hào):S158????? 文獻(xiàn)標(biāo)志碼:A文章編號(hào):1000-2367(2024)03-0139-11

    河南沿黃稻區(qū)每年經(jīng)歷濕地田、旱地的兩熟制管理方式,常年的稻麥輪作方式,人為干擾強(qiáng)烈,土壤有機(jī)質(zhì)含量較高[1].黃河灘地為砂質(zhì)土壤,有機(jī)質(zhì)含量相對(duì)較低,人為干擾少,每年7至10月的黃河主汛期,多數(shù)灘地被季節(jié)性淹沒[2],其他時(shí)間河床裸露干燥,呈旱地特征,并著生以莎草科為主的濕生草本植物.黃河灘地與毗鄰稻田具有相似的干濕交替周期,但二者又有截然不同的土地管理方式,主要是黃河灘地取樣點(diǎn)位于近河岸處,受上游來水影響大,且無耕作歷史.因此,上述兩種土樣可為微生物菌群與生態(tài)功能研究提供良好的材料.

    土壤微生物在生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量轉(zhuǎn)化過程中發(fā)揮著重要作用,其中微生物群落結(jié)構(gòu)與土壤性質(zhì)互為影響.土地利用方式可以改變土壤質(zhì)量和養(yǎng)分循環(huán),又進(jìn)一步影響著土壤微生物組成[3].不穩(wěn)定有機(jī)碳是維持土壤微生物群落的重要功能性環(huán)境因子,并影響著土壤有機(jī)碳變化和貯存[4].河岸土地利用類型的不同也會(huì)導(dǎo)致土壤微生物群落組成差異,其中土壤硝酸鹽(NO-3)和水分含量是參與氮循環(huán)微生物群落的主要驅(qū)動(dòng)因素[5].農(nóng)業(yè)、非農(nóng)業(yè)用地與原生草地相比,其共存的微生物在群落結(jié)構(gòu)和功能上有所不同[6].植被類型不同,如以小麥、玉米替代向日葵單作也會(huì)顯著影響根際微生物群落組成[7].JANGID等[8]研究發(fā)現(xiàn),濕地轉(zhuǎn)變?yōu)檗r(nóng)田后引起細(xì)菌豐度和種類的顯著變化,并確定土地利用變化是微生物群落組成的主要決定因素.微生

    收稿日期:2023-06-02;修回日期:2023-09-08.

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(21906041);河南省科技攻關(guān)項(xiàng)目(232102320100,212102310079);河南省高等學(xué)校重點(diǎn)科研項(xiàng)目(21A610003);河南工學(xué)院博士基金項(xiàng)目(9001/KQ1811).

    作者簡(jiǎn)介(通信作者):田瑩瑩(1988-),女,河南許昌人,河南工學(xué)院講師,博士,研究方向?yàn)槲廴究刂婆c生態(tài)修復(fù)研究,E-mail:tyy@hait.edu.cn.

    引用本文:田瑩瑩,郭琦,陳星儒,等.沿黃稻區(qū)土壤菌群結(jié)構(gòu)與生態(tài)功能研究[J].河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2024,52(3):139-149.(Tian Yingying,Guo Qi,Chen Xingru,et al.Soil microbial community structure and eco-function in the rice zone along the Yellow River[J].Journal of Henan Normal University(Natural Science Edition),2024,52(3):139-149.DOI:10.16366/j.cnki.1000-2367.2023.06.02.0007.)

    物參與特定的土壤生化反應(yīng)[9],反過來影響土壤生態(tài)系統(tǒng)的多功能性[10-12],包括土壤養(yǎng)分循環(huán)、植物生長(zhǎng)、溫室氣體產(chǎn)生等方面.研究表明,氨氧化細(xì)菌在農(nóng)業(yè)生態(tài)系統(tǒng)的硝化過程中起主導(dǎo)作用[5].水稻-蔬菜長(zhǎng)期輪作導(dǎo)致土壤微生物組結(jié)構(gòu)均勻化,微生物的硝化和反硝化能力增加,潛在的N2O排放增加[9].研究顯示,微生物多樣性受土地利用類型改變的影響較小,而微生物組成和結(jié)構(gòu)對(duì)土地利用類型的響應(yīng)相對(duì)更敏感[6,13].在黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展背景下,從機(jī)理上科學(xué)認(rèn)識(shí)沿黃稻區(qū)土壤微生物群落結(jié)構(gòu)與生態(tài)功能有助于更好地理解土壤內(nèi)部物質(zhì)循環(huán)過程與驅(qū)動(dòng)機(jī)制,對(duì)于制定科學(xué)合理的農(nóng)田可持續(xù)發(fā)展模式起著至關(guān)重要的作用.

    本研究以沿黃流域(原陽(yáng)段)黃河灘地和稻田土為研究對(duì)象,利用高通量測(cè)序技術(shù)揭示沿黃稻區(qū)菌群結(jié)構(gòu)、多樣性和差異物種,并基于BugBase、FAPROTAX數(shù)據(jù)庫(kù)預(yù)測(cè)菌群表型和功能.研究結(jié)果為認(rèn)識(shí)和理解沿黃稻區(qū)土壤微生物、功能及元素循環(huán)機(jī)制提供重要參考.

    1? 材料和方法

    1.1? 研究區(qū)域及土壤樣品的采集

    于2021年11月初采用梅花形5點(diǎn)采樣法采集位于河南省原陽(yáng)縣沿黃稻區(qū)表層1 cm的嫩灘砂質(zhì)土(sand soil,SS; 34°55′1.59″ E,113°40′25.01″ N)和稻田土(paddy soil,PS; 35°01′28.442″ E, 113°39′45.714″ N).每個(gè)土壤類型采集3個(gè)點(diǎn),其中SS樣品記為S1_1、S2_1和S3_1,PS樣品記為P1_1、P2_1和P3_1.采樣點(diǎn)S1_1距離河岸約200 m(有植被覆蓋),S2_1距離河岸約50 m(有0.5 cm淤泥層覆蓋),S3_1為近河岸約1 m(受黃河水漲落影響最大).PS采樣點(diǎn)距河岸直線距離約12 km,各點(diǎn)間直線距離100 m.在水稻地上部分收割后采樣,采集時(shí)盡量避開植株根系.所有土壤樣品采集后立即放置在冷藏箱中保存,并及時(shí)運(yùn)回實(shí)驗(yàn)室分析,其中用于土壤理化性質(zhì)分析的土樣置于4 ℃暫存并于48 h內(nèi)分析完畢,用于高通量測(cè)序的樣品在處理前置于-80 ℃冰箱保存.

    1.2? 土壤理化性質(zhì)的測(cè)定

    取適量風(fēng)干土樣,采用水土質(zhì)量比2.5∶1,搖床振蕩30 min,靜置后用便攜式pH計(jì)測(cè)定土壤pH值.土壤中的二價(jià)鐵(FeII)采用鄰菲羅啉比色法測(cè)定,總鐵(FeT)利用鹽酸羥胺還原后采用鄰菲羅啉比色法測(cè)定.土壤總有機(jī)碳(TOC)采用濕式重鉻酸鉀分光光度法測(cè)定.土壤銨態(tài)氮(NH+4-N)、硝態(tài)氮(NO-3-N)、亞硝態(tài)氮(NO-2-N)的測(cè)定利用2 mol/L KCl溶液浸提,其中NH+4-N采用靛酚藍(lán)比色法,NO-3-N采用紫外分光光度法,NO-2-N采用磺胺/鹽酸萘乙二胺分光光度法測(cè)定.土壤全氮(TKN)利用凱氏定氮法測(cè)定[14].

    1.3? DNA的提取,PCR擴(kuò)增與高通量測(cè)序

    按E.Z.N.A.soil DNA kit(Omega Bio-tek,Norcross)試劑盒說明書步驟提取基因組DNA,使用NanoDrop 2000超微量蛋白質(zhì)核酸分析儀(Thermo Fisher Scientific,USA)測(cè)定提取出的DNA濃度和純度.將提取好的DNA樣品利用16S rRNA基因V3-V4區(qū)擴(kuò)增引物進(jìn)行PCR擴(kuò)增,引物序列為338F(5'-ACTCCTACGGGAGGCAGCAG-3')和806R(5'-GGACTACHVGGGTWTCTAAT-3').使用質(zhì)量分?jǐn)?shù)2%瓊脂糖凝膠回收PCR產(chǎn)物,利用AxyPrep DNA Gel Extraction Kit(Axygen Biosciences,Union City)進(jìn)行回收產(chǎn)物純化.

    純化后的產(chǎn)物送至上海美吉生物醫(yī)藥科技有限公司,使用Illumina Miseq PE300平臺(tái)進(jìn)行高通量測(cè)序和分析.使用NEXTFLEX Rapid DNA-Seq Kit進(jìn)行建庫(kù).使用Fastq軟件對(duì)原始測(cè)序序列進(jìn)行質(zhì)控,并用Flash軟件進(jìn)行拼接.使用Qiime2(v2020.2)軟件,利用DADA2算法對(duì)序列進(jìn)行ASV聚類并剔除嵌合體.基于Silva 16S rRNA 數(shù)據(jù)庫(kù)(v138),利用RDP classifier貝葉斯算法對(duì)ASV代表序列進(jìn)行物種分類學(xué)注釋(置信度閾值0.7),獲得物種分類學(xué)注釋結(jié)果.各樣品的16S rRNA基因測(cè)序結(jié)果上傳至國(guó)家微生物科學(xué)數(shù)據(jù)中心(https://nmdc.cn/),登錄號(hào)NMDC40024174~NMDC40024179.

    1.4? 數(shù)據(jù)分析

    基于細(xì)菌群落及土壤理化因子數(shù)據(jù)利用R語(yǔ)言Psych數(shù)據(jù)包繪制相關(guān)性熱圖.基于美吉云平臺(tái),利用ggplot 2包繪制物種豐度圖,利用Kruskal-Wallis秩和檢驗(yàn)、多重比較方差分析檢驗(yàn)菌群的alpha多樣性,利用主坐標(biāo)分析(PCoA)、Bray-Curtis距離進(jìn)行beta多樣性分析,與BugBase數(shù)據(jù)庫(kù)比對(duì)預(yù)測(cè)微生物表型;并與FAPROTAX數(shù)據(jù)庫(kù)比對(duì)預(yù)測(cè)優(yōu)勢(shì)菌屬的功能.

    2? 結(jié)果和分析

    2.1? 沿黃灘地和稻田土主要理化性質(zhì)

    對(duì)比SS和PS的主要理化性質(zhì),發(fā)現(xiàn)除pH外,PS中的FeII、FeT、NH+4-N、NO-3-N、NO-2-N、TKN和TOC的含量均顯著高于SS.PS中pH的平均值為7.90,低于SS的8.53.研究表明,土壤有機(jī)質(zhì)含量有隨pH升高而降低的趨勢(shì),二者間呈極顯著的負(fù)相關(guān)關(guān)系(R =-0.530~-0.332,P<0.001).本研究中,TOC與pH顯著負(fù)相關(guān)(R=-0.881,P=0.02).

    2.2? 沿黃灘地和稻田土中菌群結(jié)構(gòu)

    將16S rRNA基因測(cè)序結(jié)果(去除葉綠體和線粒體)按樣本最小序列數(shù)抽平,構(gòu)建ASV集(ASV_all).將相對(duì)豐度小于0.01的Genus合并,共有60屬.alpha多樣性指數(shù)反應(yīng)群落內(nèi)物種數(shù)量及其相對(duì)豐度,有豐富度指數(shù)(如Chao指數(shù))和物種多樣性指數(shù)(如Shannon、Simpson指數(shù))等.Chao指數(shù)值越高,說明物種數(shù)目較多;Shannon指數(shù)值越高,說明群落多樣性越高;Simpson指數(shù)值越高,則說明群落多樣性越低.分析表明(圖1),SS物種數(shù)目較多而多樣性低,但二者的3個(gè)指數(shù)值均無顯著性差異(P>0.05).對(duì)微生物組beta多樣性以主坐標(biāo)分析(PCoA)進(jìn)行計(jì)算并可視化.SS各點(diǎn)分散,PS各點(diǎn)距離近,表明SS的群落結(jié)構(gòu)差異較大,而PS的菌群組成較為相似,這與SS樣點(diǎn)間生境差異大有關(guān).主成分軸PC1和PC2分別解釋了49.24%和23.06%的結(jié)構(gòu)差異(圖2).

    去除至少在4個(gè)樣本中相對(duì)豐度都小于0.001 7的ASV,得到相對(duì)豐度前35的屬,以此繪制相對(duì)豐度熱圖(圖3(a)),其中前35屬在總序列中的平均相對(duì)豐度為0.439 4±0.030 4.由于S1_1、S2_1、S3_1樣點(diǎn)理化性質(zhì)(如NH+4-N、TOC、pH等)差異較大,使得SS中組成物種的重現(xiàn)性不如PS高.兩種類型土壤中的優(yōu)勢(shì)細(xì)菌在門水平上都分屬于放線菌門(Actinobacteriota)和變形菌門(Proteobacteria),且相對(duì)豐度無差異.屬水平上,菌群的相對(duì)豐度有所不同.Unclassified_f? Micrococcaceae、鞘氨醇單胞菌屬(Sphingomonas)、norank_f? noorank_o? Vicinamibacterales和norank_f? Vicinamibacterales為土壤中的廣布屬,相對(duì)豐度與土壤類型無關(guān).SS中unclassified_f_Planococcaceae、黃桿菌屬(Flavobacterium)、副球菌屬(Paracoccus)、unclassified_f_Comamonadaceae、紅細(xì)菌屬(Rhodobacter)、Nodosilinea_PCC_7104相對(duì)豐度高于PS;PS中類諾卡氏菌屬(Nocardioides)、unclassified_f_Nocardioidaceae、硫桿菌屬(Thiobacillus)、unclassified_f_Intrasporangiaceae、芽孢桿菌屬(Bacillus)、norank_f_Anaerolineaceae、norank_f_Caldilineaceae、norank_f_norank_o_SJA-15、norank_f_norank_o_norank_c_MB-A2-108、Gaiella、norank_f_Sutterelleaceae相對(duì)豐度高于SS(圖3(b)).norank_f_Anaerolineaceae、norank_f_Caldilineaceae、norank_f_norank_o_SJA-15、norank_f_norank_o_norank_c_KD4_96屬綠彎菌門(Chloroflexi),在PP中的豐度顯著高于SS.Chlorolema也為Chloroflexi中的一屬,但僅在S1_1中有較高的相對(duì)豐度.此外,Porphyrobacter、Unclassified_f_Rhodobacteraceae、Rhodobacter作為不產(chǎn)氧光合細(xì)菌,能通過H2、硫化物、FeII等還原物質(zhì)提供電子供體還原CO2進(jìn)行光合自養(yǎng)生長(zhǎng),因而在低有機(jī)質(zhì)的SS中呈現(xiàn)較高的相對(duì)豐度.

    進(jìn)一步地用線性判別(LEfSe)分析確定了組間屬水平上豐度差異顯著的21屬(P<0.05).如附錄圖S1,norank_f_Caldilineaceae、norank_f_norank_o_SJA-15、norank_f_Anaerolineaceae、norank_f_norank_o_norank_c_KD4_96、norank_f_Steroidobacteraceae、norank_f_Sutterellaceae、Thiobacillus等的豐度在PS中顯著高于SS,而unclassified_f_Comamonadaceae、norank_f_Rhizobiales_Incertae_Sedis、Porphyrobacter、unclassified_f_Rhodobacteraceae、Paracoccus、微小桿菌屬(Exiguobacterium)、unclassified_f_Planococcaceae、Flavobacterium等的豐度在PS中顯著低于SS.

    2.3? 細(xì)菌優(yōu)勢(shì)組成物種功能預(yù)測(cè)分析

    BugBase表型預(yù)測(cè)結(jié)果顯示,35個(gè)土壤優(yōu)勢(shì)菌屬涵蓋了全部7個(gè)類型的微生物群落,除厭氧、生物膜形成外的其他功能菌的比例在SS和PS之間無顯著性差異(圖4(a)).選取與生物膜形成和氧需求相關(guān)的類型分析各自組成細(xì)菌的貢獻(xiàn)度,如圖4(b-e)).科水平上,土壤優(yōu)勢(shì)菌中的嚴(yán)格厭氧菌包括norank_o? Vicinamibacterales、Vicinamibacteraceae,其在PS中相對(duì)豐度顯著高于灘地土.鞘脂單胞菌科(Sphingomonadaceae)、類諾卡氏菌科(Nocardioidaceae)、norank_o_norank_c_MB-A2-108、Rhizobiales_Incertae_Sedis為兩種土壤中貢獻(xiàn)度前四的好氧菌.相對(duì)優(yōu)勢(shì)兼性厭氧菌包括微桿菌科(Microbacteriaceae)、嗜氫菌科(Hydrogenophilaceae)、紅細(xì)菌科(Rhodobacteraceae),其中Hydrogenophilaceae在PS中相對(duì)豐度高,Rhodobacteraceae在SS中相對(duì)豐度高,并是參與SS生物膜形成的重要類群.

    為研究土壤優(yōu)勢(shì)菌屬在生物地球化學(xué)循環(huán)中的作用,利用FAPROTAX預(yù)測(cè)確定了12個(gè)功能類群(圖5).它們主要參與土壤碳、硫元素的轉(zhuǎn)化和提供能量,其中化能異養(yǎng)型、需氧化能異養(yǎng)型功能類群占絕對(duì)優(yōu)勢(shì).氧化硫和硫化物的類群在兩種土壤中的比例差異顯著(P<0.01),其他類群無差異.此外,P2_1,P3_1中具PAHs降解功能的細(xì)菌比例相對(duì)升高.S3_1中光養(yǎng)生物占優(yōu)勢(shì),包括光合自養(yǎng)和光合異養(yǎng)兩類,分別為能產(chǎn)生氧氣的藍(lán)藻(主要是浮絲藻屬)和不產(chǎn)生氧氣的光合細(xì)菌(主要是Porphyrobacter、unclassified_f_Rhodobacteraceae、Rhodobacter)貢獻(xiàn),尤以光合異養(yǎng)功能細(xì)菌為主.

    3? 討? 論

    3.1? 沿黃流域土壤細(xì)菌群落結(jié)構(gòu)分析

    本研究中SS與PS的土壤類型不同,孕育了某些特定的微生物類群.SS中unclassified_f_Planococcaceae的優(yōu)勢(shì)序列(ASV2)與動(dòng)球菌屬(Planococcus sp.)相似性高.Planococcus為好氧菌,GAN等[15]從鹽堿土壤中分離得到的Planococcus最適pH為9.0.SS為砂質(zhì)土,透氣性良好,pH值接近9.0,可能是unclassified_f_Planococcaceae相對(duì)豐度高的關(guān)鍵.Paracoccus廣泛存在于多種生境,對(duì)土壤中多環(huán)芳烴(PAHs)、鄰苯二甲酸二甲酯等污染有較好的降解作用[16-17].脫氮副球菌(Paracoccus denitrificans)是Paracoccus的模式種,也是重要的土壤好氧反硝化菌.SS中Paracoccus豐度高于PS,說明SS可能形成獨(dú)特的環(huán)境,較OMs豐富、N源充足的PS更適宜Paracoccus的生長(zhǎng)繁殖.Flavobacterium也是兩種類型土壤中的差異物種.研究發(fā)現(xiàn),F(xiàn)lavobacterium是土壤中魚腥藻屬(Anabaena sp.)PCC7120的重要伴生菌,具有適應(yīng)寡營(yíng)養(yǎng)的特性,可能是SS中相對(duì)豐度高的原因之一[18].

    PS中類諾卡氏菌Nocardioides、Bacillus、Thiobacillus豐度顯著高于SS.土壤中Nocardioides的某些種參與N、S元素循環(huán)[19-20].Bacillus是水稻田常見的組成微生物[21-22].枯草芽孢桿菌(Bacillus subtilis)DBM與礦物成分(高嶺石、針鐵礦)形成的復(fù)合體提高了礦物本身對(duì)銅(Cu)和鉛(Pb)的吸附能力,并保護(hù)了DBM免受重金屬毒害[23].Thiobacillus是自然生境中重要的化能自養(yǎng)菌,主要參與氧化硫化物(H2S,S2O2-3)或單質(zhì)硫?yàn)榱蛩?本研究中Thiobacillus與pH極顯著負(fù)相關(guān)(P<0.01),與FeII、FeT正相關(guān)(P<0.05),表明偏酸性環(huán)境和高含量FeII會(huì)更適宜氧化亞鐵硫桿菌(Thiobacillus ferrooxidans)生存,說明T. ferrooxidans可能是Thiobacillus中的優(yōu)勢(shì)種之一.此外,Thiobacillus與NH+4-N顯著正相關(guān),其優(yōu)勢(shì)序列ASV42與環(huán)境樣本中一反硝化Thiobacillus(登錄號(hào)MG801484.1)相似性(>99%)高,說明ASV42可能具有脫N作用.因此,推測(cè)T. ferrooxidans、T. denitrificans是本研究中Thiobacillus的重要組成物種.此外,Bacillus、Thiobacillus是重要的發(fā)酵型鐵環(huán)原微生物[24],有效改善石灰性土(pH 7.5~8.5)缺鐵狀況[25].土壤中Thiobacillus豐度可直接影響C、N、Fe、S元素轉(zhuǎn)化.Thiobacillus是PS中的優(yōu)勢(shì)菌,因此在稻田干濕周期中Thiobacillus與土壤元素轉(zhuǎn)化的關(guān)系的研究對(duì)稻田的精細(xì)化管理具有重要的指導(dǎo)意義.

    光養(yǎng)微生物有光合產(chǎn)氧和光合不產(chǎn)氧兩類.與PS相比,SS中光合自養(yǎng)微生物總的相對(duì)豐度高,其中又以光合不產(chǎn)氧細(xì)菌的豐度占優(yōu)勢(shì).SS中Symphothece PCC 702、Planthorix為光合產(chǎn)氧細(xì)菌,Chloronema、Rhodobacter和unclassified_f_Rhodobacteraceae為厭氧不產(chǎn)氧光合細(xì)菌.研究發(fā)現(xiàn),Porphyrobacter屬中也包含有好氧不產(chǎn)氧的光合類群[26].SS中Porphyrobacter、Chloronema、Rhodobacter和unclassified_f_Rhodobacteraceae相對(duì)豐度高,主要是因?yàn)楣廨椛渥饔孟?,它們能以H2、硫化物、FeII等還原物質(zhì)提供的電子供體還原CO2光合自養(yǎng)生長(zhǎng),因而在有機(jī)質(zhì)低的灘地環(huán)境中較異養(yǎng)代謝微生物呈現(xiàn)出一定的生存優(yōu)勢(shì).Rhodobacter是稻田土中常見的優(yōu)勢(shì)光合細(xì)菌,大多能隨外界氧氣、有機(jī)質(zhì)濃度的改變進(jìn)行好氧、微需氧或光合厭氧代謝[27].C是植物必需的六種大量元素之首.土壤中的C固存與土壤肥力密切相關(guān).光氧微生物是土壤有機(jī)碳的重要生產(chǎn)者.為提高沿黃稻田有機(jī)碳含量,結(jié)合測(cè)序結(jié)果,可在沿黃稻田干濕周期分別接種好氧Porphyrobacter和厭氧Chloronema、Rhodobacter實(shí)現(xiàn)生物固碳以減少外源有機(jī)肥使用量.

    3.2? 土壤優(yōu)勢(shì)物種的功能及生物地球化學(xué)作用

    土壤菌群組成與其環(huán)境適應(yīng)性有關(guān),并影響著生境內(nèi)元素的生物地球化學(xué)循環(huán).基于BugBase表型預(yù)測(cè)的結(jié)果顯示,優(yōu)勢(shì)菌屬包括了全部的7類表型.嚴(yán)格厭氧功能菌是PS和SS中的重要差異類群,具體為,norank_o? Vicinamibacterales、Vicinamibacteraceae在PS中相對(duì)豐度顯著高于SS.Microbacteriaceae、Hydrogenophilaceae、Rhodobacteraceae是優(yōu)勢(shì)兼性厭氧菌,此外,Rhizobiales_Incertae_Sedis、Rhodobacteraceae與生物膜形成相關(guān),它們?cè)赟S中的相對(duì)豐度均高于PS.基于FAPROTAX數(shù)據(jù)庫(kù)預(yù)測(cè)確定了12個(gè)功能類群,它們主要參與土壤碳、硫等的轉(zhuǎn)化和提供能量.化能異養(yǎng)和需氧化能異養(yǎng)型微生物在總微生物類群中占絕對(duì)優(yōu)勢(shì).氧化硫和硫化物的微生物在兩類土壤中的比例有極顯著差異.SS中光養(yǎng)微生物的相對(duì)豐度普遍高于PS,尤以S3_1(Porphyrobacter、unclassified_f_Rhodobacteraceae、Rhodobacter構(gòu)成)為最高.光合自養(yǎng)微生物的存在提高了光能利用率,供給了部分有機(jī)碳源,對(duì)維持生境內(nèi)穩(wěn)定的菌群結(jié)構(gòu)起重要作用.PS中具硫、硫化物氧化和PAHs降解功能的微生物的比例也相對(duì)較高.T. ferrooxidans、T. denitrificans為Thiobacillus的重要組成物種,通過代謝硫和硫化物影響體系內(nèi)Fe、N元素的轉(zhuǎn)化.微生物降解作用是環(huán)境中PAHs去除的主要途徑[28].與焦化廠污染的荒地相比,農(nóng)田土壤理化特性增加了PAHs降解菌和總細(xì)菌群落的多樣性,與低分子量PAHs相關(guān)的降解菌主要隸屬于變形桿菌,與高分子PAHs相關(guān)的降解菌主要隸屬于放線菌[28].FAPROTAX功能熱圖分析結(jié)果顯示,PS(尤其是P2_1、P3_1)中PAHs降解功能微生物的相對(duì)豐度高于SS.Sphingomonas生境廣,適應(yīng)性強(qiáng),對(duì)復(fù)雜結(jié)構(gòu)的難降解芳香化合物有特異作用,在生物聚合物的生產(chǎn)方面有優(yōu)越性.本研究中Sphingomonas豐度在兩種不同利用方式的土壤中無顯著性差異,說明PS中可能有其他的具芳香烴功能的微生物存在.研究發(fā)現(xiàn),農(nóng)業(yè)土壤中Nocardioides也是具降解PAHs功能的微生物的重要類群之一[29].本研究中Nocardioides、unclassified_f_Nocardioidaceae在P2_1、P3_1中的相對(duì)豐度顯著高于其他樣點(diǎn),與PAHs降解功能的熱圖分析結(jié)果一致,表明Nocardioides、unclassified_f_Nocardioidaceae可能是促使PAHs降解的關(guān)鍵類群.

    unclassified_f? Micrococcaceae為兩種類型土壤中的絕對(duì)優(yōu)勢(shì)屬.本研究中unclassified_f? Micrococcaceae的優(yōu)勢(shì)代表序列ASV1和ASV3與節(jié)桿菌屬(Arthrobacter sp.)相似性達(dá)100%.Arthrobacter菌株嚴(yán)格好氧,具有高效降解有機(jī)污染物和吸附重金屬的能力[30].Arthrobacter sp.ADH-2在10 h內(nèi)對(duì)100 mg/L阿特拉津的降解率為99.9%[31].研究顯示,unclassified_f? Micrococcaceae是健康與發(fā)病煙田微生物的主要差異物種之一,在健康煙田土壤中的相對(duì)豐度是發(fā)病土壤的2倍[32].此外,unclassified_f? Micrococcaceae也是間作甘蔗植株中特有的內(nèi)生菌之一[33].

    生態(tài)系統(tǒng)中存在多元素的耦合作用,因此,參與碳、硫等元素代謝的微生物往往也介導(dǎo)了生境內(nèi)的N轉(zhuǎn)化過程.Paracoccus denitrificans厭氧條件下以H2作為電子供體,利用CO2、HCO-3等無機(jī)碳源,還原硝酸鹽到亞硝酸鹽到氧化氮(NO、N2O)和N2[34].該過程是生物地球化學(xué)氮循環(huán)中重要的反硝化反應(yīng),方程式為:①5H2+2NO-3→N2+4H2O+2OH-;H2+0.35NO-3+0.35H++0.052CO2→0.17N2+1.1H2O+0.010C5H7O2N[34].好氧條件下,P. denitrificans對(duì)溶解性硝酸鹽、亞硝酸鹽的去除以同化作用為主,少部分經(jīng)由反硝化作用去除[35].學(xué)者們也發(fā)現(xiàn)了同時(shí)具有異養(yǎng)硝化-好氧反硝化功能的菌株,為P. denitrificans用于廢水短程脫氮提供參考[36].某些Rhodobacter菌株以NO-3-N為電子受體,有機(jī)碳源為電子供體光合異養(yǎng)生長(zhǎng),是沉積物中一種關(guān)鍵的反硝化菌[37].研究報(bào)道,Rhodobacter是汾河下游水體中nir S型反硝化細(xì)菌群落的3個(gè)主導(dǎo)菌屬之一,在下游某些采樣點(diǎn)的相對(duì)豐度可高于98%,其相對(duì)豐度與NO-3-N負(fù)相關(guān)[38].本研究中Rhodobacter與NO-3-N也顯著負(fù)相關(guān),說明所調(diào)查區(qū)域的Rhodobacter中可能存在反硝化脫氮菌株.T. denitrificans利用元素硫?yàn)殡娮庸w從水生環(huán)境中去除NO-3,是重要的硫自養(yǎng)反硝化菌.溶解性有機(jī)質(zhì)作為一種微生物光敏劑驅(qū)動(dòng)T. denitrificans發(fā)生光電反硝化反應(yīng)[39].

    4? 結(jié)? 論

    本研究利用高通量測(cè)序技術(shù)揭示了沿黃稻區(qū)SS和PS土樣中的菌群結(jié)構(gòu)、多樣性及差異物種,發(fā)現(xiàn)兩種土樣中的菌群豐富度和多樣性指數(shù)差異不顯著,但部分菌群相對(duì)豐度差異顯著.unclassified_f? Micrococcaceae、Sphingomonas為土壤中的廣布種;Bacilllus、Thiobacillus在PS中的相對(duì)豐度顯著高于SS;PS中光合細(xì)菌的相對(duì)豐度則顯著低于SS.特別地,與硫代謝相關(guān)的Thiobacillus,推測(cè)主要由T. denitrificans和T. ferrooxidans組成,而Rhodobacteraceae則是土壤優(yōu)勢(shì)菌屬中重要的促生物膜形成微生物.基于BugBase、FAPROTAX數(shù)據(jù)庫(kù)的菌群表型和功能預(yù)測(cè)結(jié)果顯示,SS和PS均包括了7類表型和12類地化循環(huán)相關(guān)微生物,其中厭氧、氧化硫或硫化物功能的類群在土壤中的相對(duì)豐度差異顯著.

    但本研究樣本量較少,對(duì)SS和PS土壤微生物的差異性認(rèn)識(shí)可能不全面.未來,需要更全面的研究以揭示土壤類型、土地利用方式等影響下的微生物群落結(jié)構(gòu)與功能.

    附錄見電子版(DOI:10.16366/j.cnki.1000-2367.2023.06.02.0007).

    參? 考? 文? 獻(xiàn)

    [1] ??孫建權(quán),胡寧,王書玉,等.水稻旱直播高產(chǎn)栽培技術(shù)初探[J].種業(yè)導(dǎo)刊,2009(10):22-25.

    SUN J Q,HU N,WANG S Y,et al.Preliminary study on high-yield cultivation techniques of dry direct seeding of rice[J].Journal of Seed Industry Guide,2009(10):22-25.

    [2]李榮容,王鵬,韋詩(shī)濤,等.黃河中下游汛期洪水的分期識(shí)別[J].人民黃河,2017,39(11):24-29.

    LI R R,WANG P,WEI S T,et al.Flood season staged in middle and lower Yellow River[J].Yellow River,2017,39(11):24-29.

    [3]LI Y T,XIE X H,ZHU Z J,et al.Land use driven change in soil organic carbon affects soil microbial community assembly in the riparian of Three Gorges Reservoir Region Applied Soil Ecology,2022,176:104467.

    [4]ZHANG H K,F(xiàn)ANG Y Y,ZHANG B G,et al.Land-use-driven change in soil labile carbon affects microbial community composition and function[J].Geoderma,2022,426:116056.

    [5]MAFA-ATTOYE T G,BASKERVILLE M A,OFOSU E,et al.Riparian land-use systems impact soil microbial communities and nitrous oxide emissions in an agro-ecosystem[J].The Science of the Total Environment,2020,724:138148.

    [6]BISSETT A,RICHARDSON A E,BAKER G,et al.Long-term land use effects on soil microbial community structure and function[J].Applied Soil Ecology,2011,51:66-78.

    [7]WEN X Y,DUBINSKY E,WU Y,et al.Wheat,maize and sunflower cropping systems selectively influence bacteria community structure and diversity in their and succeeding crop's rhizosphere[J].Journal of Integrative Agriculture,2016,15(8):1892-1902.

    [8]JANGID K,WILLIAMS M A,F(xiàn)RANZLUEBBERS A J,et al.Land-use history has a stronger impact on soil microbial community composition than aboveground vegetation and soil properties[J].Soil Biology and Biochemistry,2011,43(10):2184-2193.

    [9]SIM J X F,DOOLETTE C L,VASILEIADIS S,et al.Pesticide effects on nitrogen cycle related microbial functions and community composition[J].The Science of the Total Environment,2022,807:150734.

    [10]SHEN J E,TAO Q,DONG Q,et al.Long-term conversion from rice-wheat to rice-vegetable rotations drives variation in soil microbial communities and shifts in nitrogen-cycling through soil profiles[J].Geoderma,2021,404:115299.

    [11]CHEN W Q,WANG J Y,CHEN X A,et al.Soil microbial network complexity predicts ecosystem function along elevation gradients on the Tibetan Plateau[J].Soil Biology and Biochemistry,2022,172:108766.

    [12]TANG Q A,XIA Y Q,TI C P,et al.Partial organic fertilizer substitution promotes soil multifunctionality by increasing microbial community diversity and complexity[J].Pedosphere,2023,33(3):407-420.

    [13]張俊,黃煦旸,彭宇琦.“雙碳”目標(biāo)下黃河流域種植業(yè)碳足跡的時(shí)空異質(zhì)效應(yīng)[J].河南農(nóng)業(yè)大學(xué)學(xué)報(bào),2023,57(5):884-896.

    ZHANG J,HUANG X Y,PENG Y Q.Spatio-temporal heterogeneity of planting carbon footprint in Yellow River Basin under the carbon peaking and carbon neutrality targets[J].Journal of Henan Agricultural University,2023,57(5):884-896.

    [14]榮國(guó)華,周景云,吳鴻宇,等.流動(dòng)分析儀與凱氏定氮儀測(cè)定土壤全氮含量之比較研究[J].水土保持研究,2023,30(1):204-208.

    RONG G H,ZHOU J Y,WU H Y,et al.Comparative study on determination of total nitrogen content in soil by flow analyzer and Kjeldahl nitrogen analyzer[J].Research of Soil and Water Conservation,2023,30(1):204-208.

    [15]GAN L Z,ZHANG H M,TIAN J W,et al.Planococcus salinus sp. nov.,a moderately halophilic bacterium isolated from a saline-alkali soil[J].International Journal of Systematic and Evolutionary Microbiology,2018,68(2):589-595.

    [16]毛健,駱永明,滕應(yīng),等.一株副球菌對(duì)污染土壤中多環(huán)芳烴的降解研究[J].土壤,2009,41(3):448-453.

    MAO J,LUO Y M,TENG Y,et al.Biodegradation of PAHs by Paracoccus aminovorans HPD-2 in contaminated soil[J].Soils,2009,41(3):448-453.

    [17]朱曉慧,劉虹,郭茹鑫,等.黑土微生物群落對(duì)副球菌屬(Paracoccus sp.)QD15-1修復(fù)鄰苯二甲酸二甲酯污染土壤的響應(yīng)[J].環(huán)境污染與防治,2020,42(2):182-186.

    ZHU X H,LIU H,GUO R X,et al.Response of microbial community in black soil to remediation of dimethyl phthalate contaminated soil by Paracoccus sp.QD15-1[J].Environmental Pollution & Control,2020,42(2):182-186.

    [18]荊瑞勇,王麗艷,孫強(qiáng),等.東北旱田土壤中Anabaena伴生細(xì)菌的分離與鑒定[J].微生物學(xué)通報(bào),2020,47(1):130-139.

    JING R Y,WANG L Y,SUN Q,et al.Isolation and identification of bacteria associated with Anabaena from upland soils in northeast China[J].Microbiology China,2020,47(1):130-139.

    [19]CUI Y S,WOO S G,LEE J,et al.Nocardioides daeguensis sp. nov.,a nitrate-reducing bacterium isolated from activated sludge of an industrial wastewater treatment plant[J].International Journal of Systematic and Evolutionary Microbiology,2013,63(Pt_10):3727-3732.

    [20]LACHMAYR K,ODUARAN E.Investigation of Nocardioides sp. metabolic,structural,biochemical properties,and its potential role in the sulfur cycle[J].The FASEB Journal,2015,29(S1):1944-1956.

    [21]曹夢(mèng)雨,桑世飛,王亞男,等.水稻全基因組磷脂酶家族蛋白篩選及其生物信息學(xué)分析[J].河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2022,50(3):143-149.

    CAO M Y,SANG S F,WANG Y N,et al.Genome wide screening and bioinformatics analysis of rice phospholipase family proteins[J].Journal of Henan Normal University(Natural Science Edition),2022,50(3):143-149.

    [22]DING M J,SHANG N J,XIAO Z X,et al.Bacillus aciditolerans sp. nov.,isolated from paddy soil[J].International Journal of Systematic and Evolutionary Microbiology,2019,69(4):1155-1161.

    [23]BAI J,CHAO Y,CHEN Y M,et al.The effect of interaction between Bacillus subtilis DBM and soil minerals on Cu(II)and Pb(II)adsorption[J].Journal of Environmental Sciences,2019,78(4):328-337.

    [24]田瑩瑩,王強(qiáng),趙京,等.黃河灘地和稻田土中鐵還原菌、不產(chǎn)氧光合細(xì)菌分布機(jī)制[J].微生物學(xué)報(bào),2023,63(2):805-820.

    TIAN Y Y,WANG Q,ZHAO J,et al.Distribution mechanism of iron-reducing bacteria and oxygen-producing photosynthetic bacteria in Yellow River beach and paddy soil[J].Acta Microbiologica Sinica,2023,63(2):805-820.

    [25]GHORBANZADEH N,LAKZIAN A,HAGHNIA G H,et al.Isolation and identification of ferric reducing bacteria and evaluation of their roles in iron availability in two calcareous soils[J].Eurasian Soil Science,2014,47(12):1266-1273.

    [26]LI X L,KOBLEK M,F(xiàn)ENG F,et al.Whole-genome sequence of a freshwater aerobic anoxygenic phototroph,Porphyrobacter sp. strain AAP82,isolated from the Huguangyan maar lake in southern China[J].Genome Announcement,2013,1(2):e0007213.

    [27]ZAVALA N,BAEZA L,GONZALEZ S,et al.The effects of different carbon sources on the growth of Rhodobacter sphaeroides[J].Advances in Microbiology,2019,9(8):737-749.

    [28]DU J Q,LIU J X,JIA T,et al.The relationships between soil physicochemical properties,bacterial communities and polycyclic aromatic hydrocarbon concentrations in soils proximal to coking plants[J].Environmental Pollution,2022,298:118823.

    [29]LU C,HONG Y,LIU J A,et al.A PAH-degrading bacterial community enriched with contaminated agricultural soil and its utility for microbial bioremediation[J].Environmental Pollution,2019,251:773-782.

    [30]李娟,CONSTANTINE UWAREMEW,冷艷,等.節(jié)桿菌屬細(xì)菌處理有機(jī)物和重金屬污染物的研究進(jìn)展[J].環(huán)境科學(xué)與技術(shù),2017,40(10):89-97.

    LI J,UWAREMWE C,LENG Y,et al.Progress on the study of biodegradation of organic pollutants and adsorption of heavy metals with Arthrobacter strains[J].Environmental Science & Technology,2017,40(10):89-97.

    [31]韓鵬,洪青,何麗娟,等.阿特拉津降解菌ADH-2的分離、鑒定及其特性研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2009,28(2):406-410.

    HAN P,HONG Q,HE L J,et al.Isolation,identification and characterization of an atrazine-degrading bacteria ADH-2[J].Journal of Agro-Environment Science,2009,28(2):406-410.

    [32]陳乾錦,林書震,李紅麗,等.邵武煙田土壤微生物群落結(jié)構(gòu)變化與煙草青枯病發(fā)生關(guān)系初報(bào)[J].中國(guó)煙草學(xué)報(bào),2019,25(4):64-71.

    CHEN Q J,LIN S Z,LI H L,et al.A preliminary report on relationship between variation of microbial community structure in soil and occurrence of tobacco bacterial wilt in tobacco field in Shaowu[J].Acta Tabacaria Sinica,2019,25(4):64-71.

    [33]肖健,陳思宇,孫妍,等.甘蔗間作不同豆科作物對(duì)甘蔗植株內(nèi)生細(xì)菌多樣性的影響[J].熱帶作物學(xué)報(bào),2021,42(11):3188-3198.

    XIAO J,CHEN S Y,SUN Y,et al.Effect of intercropping with different legume crops on endophytic bacterial diversity of sugarcanes[J].Chinese Journal of Tropical Crops,2021,42(11):3188-3198.

    [34]RICHARDSON D J,SAWERS G,VAN SPANNING R J M.Periplasmic electron-transport systems in bacteria[M].Amsterdam:Elsevier,2013.

    [35]劉道玉,范立民,王瓊,等.脫氮副球菌的好氧反硝化特性及對(duì)養(yǎng)殖水體中氮素的控制[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2012,31(11):2249-2255.

    LIU D Y,F(xiàn)AN L M,WANG Q,et al.The aerobic denitrification characteristics of Paracoccus denitrificans and its nitrogen control for aquaculture water[J].Journal of Agro-Environment Science,2012,31(11):2249-2255.

    [36]李海紅,佟欣宇,宦臣臣,等.高效異養(yǎng)硝化-好氧反硝化菌株TS-1篩選及降解特性[J].應(yīng)用與環(huán)境生物學(xué)報(bào),2020,26(4):791-799.

    LI H H,TONG X Y,HUAN C C,et al.Characterizing the heterotrophically nitrifying and aerobically denitrifying TS-1 strain[J].Chinese Journal of Applied and Environmental Biology,2020,26(4):791-799.

    [37]ZHOU S L,SUN Y E,LI Z X,et al.Characteristics and driving factors of the aerobic denitrifying microbial community in Baiyangdian Lake,Xiong'an new area[J].Microorganisms,2020,8(5):714.

    [38]汪銀龍,馮民權(quán),董向前.汾河下游水體nirS型反硝化細(xì)菌群落組成與無機(jī)氮關(guān)系[J].環(huán)境科學(xué),2019,40(8):3596-3603.

    WANG Y L,F(xiàn)ENG M Q,DONG X Q.Community composition of nir S-type denitrifying bacteria in the waters of the lower reaches of the Fenhe River and its relationship with inorganic nitrogen[J].Environmental Science,2019,40(8):3596-3603.

    [39]HUANG S F,CHEN M,DIAO Y M,et al.Dissolved organic matter acting as a microbial photosensitizer drives photoelectrotrophic denitrification[J].Environmental Science & Technology,2022,56(7):4632-4641.

    Soil microbial community structure and eco-function in the rice zone along the Yellow River

    Tian Yingyinga,b, Guo Qia, Chen Xingrua, Sun Xianghuia,b, Zhao Jinga,b

    (a. Department of Material Science and Engineering; b. Xinxiang City Engineering Research Center for Wastewater Treatment Energy Saving and Emission Reduction, Henan Institute of Technology, Xinxiang 453003, China)

    Abstract: In the present study, samples of riverside sand soil(SS) and adjacent paddy soil(PS) were collected to know the microbial composition of different substrate soils in the rice zone along the Yellow River. The high-throughput sequencing of 16S rRNA gene was used to uncover the soil microbial compositions and to explain microbes' potential eco-function. Alpha diversity analysis showed no significant difference in microbial richness and diversity index between the two soils. Further species analysis showed the ubiquitous distribution of unclassified_f? Micrococcaceae, Sphingomonas, norank_f? noorank_o? Vicinamibacterales and norank_f? Vicinamibacterales, whose abundances were independent of soil types. However, the relative abundances of Bacilllus and Thiobacillus were significantly higher in PS than that in SS. The function heatmap analysis based on FAPROTAX database showed that aromatic compound degrading bacteria in PS(especially P2_1 and P3_1) was higher than that in SS. Combined with the results of bacterial relative abundance, Nocardioides and unclassified_f_Nocardioidaceae might be the key genera contributing to aromatic compound degradation. The phototrophic microbes(PMs) in SS were more abundant than that in PS, which were dominant anoxygenic PMs like Porphyrobacter, unclassified_f_Rhodobacteraceae and Rhodobacter. Moreover, Rhodobacteraceae was likely to participate in biofilm formation in SS. Overall, soil microbial compositions along the Yellow River were to some extent impacted by the soil type, driving the differences of their eco-functions.

    Keywords: soil microbial community structure; eco-function; paddy soil; sand soil

    [責(zé)任編校? 劉洋? 楊浦]

    附? 錄

    午夜视频国产福利| 在线a可以看的网站| h日本视频在线播放| 日韩人妻高清精品专区| 日本色播在线视频| 国产亚洲午夜精品一区二区久久 | 国产精品久久久久久精品古装| 狂野欧美白嫩少妇大欣赏| 中文乱码字字幕精品一区二区三区| 日本一二三区视频观看| 国产亚洲av嫩草精品影院| 在线 av 中文字幕| 午夜激情福利司机影院| 美女高潮的动态| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲高清精品| 欧美极品一区二区三区四区| 久久精品国产a三级三级三级| 精品人妻一区二区三区麻豆| 99久久九九国产精品国产免费| 2021天堂中文幕一二区在线观| 亚洲国产欧美人成| 亚洲欧美日韩另类电影网站 | 免费av观看视频| 边亲边吃奶的免费视频| 亚洲精品国产成人久久av| 精品一区二区三区视频在线| 80岁老熟妇乱子伦牲交| 2022亚洲国产成人精品| 少妇裸体淫交视频免费看高清| 久久久久国产精品人妻一区二区| 大又大粗又爽又黄少妇毛片口| 肉色欧美久久久久久久蜜桃 | 日韩欧美一区视频在线观看 | h日本视频在线播放| 在线观看免费高清a一片| av免费在线看不卡| 777米奇影视久久| 六月丁香七月| 久久久久精品性色| 嘟嘟电影网在线观看| 日韩电影二区| 亚洲成色77777| 欧美日韩在线观看h| 免费电影在线观看免费观看| 青青草视频在线视频观看| 日韩免费高清中文字幕av| 一级av片app| 免费人成在线观看视频色| 亚洲色图av天堂| 2022亚洲国产成人精品| 一区二区av电影网| 久久久精品欧美日韩精品| 一本一本综合久久| 日韩 亚洲 欧美在线| 啦啦啦啦在线视频资源| 最近最新中文字幕大全电影3| .国产精品久久| 亚洲美女视频黄频| 国语对白做爰xxxⅹ性视频网站| 直男gayav资源| 黄色日韩在线| 国产男人的电影天堂91| 51国产日韩欧美| 少妇的逼水好多| 国产黄片美女视频| 国产久久久一区二区三区| 久久国内精品自在自线图片| 欧美xxxx性猛交bbbb| 国产v大片淫在线免费观看| 新久久久久国产一级毛片| 国产免费又黄又爽又色| 97超碰精品成人国产| 日韩av不卡免费在线播放| 欧美少妇被猛烈插入视频| 久久精品国产亚洲网站| www.av在线官网国产| eeuss影院久久| 日日摸夜夜添夜夜爱| 狠狠精品人妻久久久久久综合| 伦理电影大哥的女人| 成人国产麻豆网| 欧美成人一区二区免费高清观看| 中文字幕久久专区| 久久久久久久国产电影| 黄色一级大片看看| 亚洲精品成人久久久久久| 啦啦啦中文免费视频观看日本| 婷婷色av中文字幕| 国产高清不卡午夜福利| 色哟哟·www| 少妇 在线观看| 中国美白少妇内射xxxbb| 国产真实伦视频高清在线观看| 久久99热这里只频精品6学生| 蜜桃亚洲精品一区二区三区| 精品国产露脸久久av麻豆| 五月玫瑰六月丁香| 乱系列少妇在线播放| 午夜视频国产福利| 国产午夜精品久久久久久一区二区三区| 18禁在线播放成人免费| 新久久久久国产一级毛片| 欧美最新免费一区二区三区| 各种免费的搞黄视频| 下体分泌物呈黄色| 91狼人影院| 最新中文字幕久久久久| 免费观看的影片在线观看| 男人添女人高潮全过程视频| 成人黄色视频免费在线看| 国产免费视频播放在线视频| 精品久久久噜噜| 国产乱人视频| 亚洲,一卡二卡三卡| 99久久人妻综合| 精品熟女少妇av免费看| 一级片'在线观看视频| 日韩 亚洲 欧美在线| 国产成人一区二区在线| eeuss影院久久| 一区二区三区免费毛片| 80岁老熟妇乱子伦牲交| 国产精品av视频在线免费观看| videossex国产| 亚洲av成人精品一区久久| 天天躁夜夜躁狠狠久久av| 美女xxoo啪啪120秒动态图| 国产一级毛片在线| 亚洲经典国产精华液单| 九色成人免费人妻av| 天堂俺去俺来也www色官网| 天天一区二区日本电影三级| 国产精品人妻久久久影院| 少妇 在线观看| 亚洲欧美日韩东京热| 美女视频免费永久观看网站| 亚洲成色77777| 亚洲国产欧美在线一区| 男人狂女人下面高潮的视频| 精品国产一区二区三区久久久樱花 | 特级一级黄色大片| 日本wwww免费看| 国产亚洲av嫩草精品影院| 午夜视频国产福利| 亚洲综合色惰| 麻豆精品久久久久久蜜桃| 国产毛片a区久久久久| 18禁裸乳无遮挡动漫免费视频 | 成人黄色视频免费在线看| 国模一区二区三区四区视频| 国产伦理片在线播放av一区| 国产一区亚洲一区在线观看| 久热久热在线精品观看| 成人毛片60女人毛片免费| 听说在线观看完整版免费高清| 美女视频免费永久观看网站| 久久久久久九九精品二区国产| 久久综合国产亚洲精品| 成人美女网站在线观看视频| 亚洲欧美日韩无卡精品| 插阴视频在线观看视频| 伦精品一区二区三区| 精品一区二区免费观看| 水蜜桃什么品种好| 看十八女毛片水多多多| 国产亚洲最大av| 男人狂女人下面高潮的视频| 欧美精品人与动牲交sv欧美| 国产精品久久久久久av不卡| 午夜视频国产福利| 亚洲一级一片aⅴ在线观看| 最近2019中文字幕mv第一页| 国产男女超爽视频在线观看| 少妇裸体淫交视频免费看高清| 乱码一卡2卡4卡精品| 男女啪啪激烈高潮av片| 国模一区二区三区四区视频| 日韩成人伦理影院| av又黄又爽大尺度在线免费看| 成年免费大片在线观看| 国产免费一级a男人的天堂| 综合色丁香网| 亚洲丝袜综合中文字幕| 美女国产视频在线观看| 久热久热在线精品观看| 午夜日本视频在线| 春色校园在线视频观看| 黄色配什么色好看| 一级毛片 在线播放| 亚洲国产最新在线播放| 日韩av免费高清视频| 啦啦啦在线观看免费高清www| 啦啦啦中文免费视频观看日本| 国产中年淑女户外野战色| 男女那种视频在线观看| 白带黄色成豆腐渣| 国产亚洲5aaaaa淫片| 少妇的逼好多水| 亚洲精品国产av成人精品| 欧美高清性xxxxhd video| 特级一级黄色大片| 国产精品福利在线免费观看| 日韩强制内射视频| 熟女人妻精品中文字幕| 精品视频人人做人人爽| 日本免费在线观看一区| 久久久久久久久大av| 丰满少妇做爰视频| 久久久a久久爽久久v久久| 国产有黄有色有爽视频| 精品午夜福利在线看| 老司机影院成人| 精品少妇黑人巨大在线播放| 国产黄a三级三级三级人| 亚洲国产成人一精品久久久| 亚洲精品日韩在线中文字幕| 久久国内精品自在自线图片| 人人妻人人澡人人爽人人夜夜| 亚洲精品国产成人久久av| 国产乱人偷精品视频| 汤姆久久久久久久影院中文字幕| 中文乱码字字幕精品一区二区三区| 久久久久国产网址| 听说在线观看完整版免费高清| 老女人水多毛片| 男女边摸边吃奶| 亚洲精品影视一区二区三区av| 国产成人午夜福利电影在线观看| 中文精品一卡2卡3卡4更新| 国产中年淑女户外野战色| 国产白丝娇喘喷水9色精品| 国产亚洲午夜精品一区二区久久 | 五月伊人婷婷丁香| 国产中年淑女户外野战色| 美女脱内裤让男人舔精品视频| 亚洲av电影在线观看一区二区三区 | 国产亚洲5aaaaa淫片| 亚洲va在线va天堂va国产| 欧美+日韩+精品| 国产国拍精品亚洲av在线观看| 成年av动漫网址| 午夜福利高清视频| 日本欧美国产在线视频| av福利片在线观看| 69人妻影院| 成人欧美大片| www.色视频.com| 免费看av在线观看网站| 日韩一本色道免费dvd| 亚洲,一卡二卡三卡| 人妻 亚洲 视频| 日韩制服骚丝袜av| 香蕉精品网在线| 国产乱人视频| 久久久久网色| 亚洲av二区三区四区| 搡老乐熟女国产| 亚洲人成网站高清观看| 久久久久久久久大av| 又爽又黄无遮挡网站| 国内精品美女久久久久久| 男男h啪啪无遮挡| 尾随美女入室| 国国产精品蜜臀av免费| 国产日韩欧美在线精品| 亚洲伊人久久精品综合| 欧美高清性xxxxhd video| 免费看a级黄色片| 性插视频无遮挡在线免费观看| 麻豆精品久久久久久蜜桃| 成人无遮挡网站| 中文字幕制服av| 亚洲无线观看免费| 少妇人妻一区二区三区视频| 国产av国产精品国产| 99热这里只有精品一区| 99久久九九国产精品国产免费| 99久久中文字幕三级久久日本| 日本一二三区视频观看| 黄色欧美视频在线观看| 国产免费一级a男人的天堂| 男女国产视频网站| 久热久热在线精品观看| 久久鲁丝午夜福利片| 久久久久久伊人网av| 色婷婷久久久亚洲欧美| 国产精品人妻久久久久久| 2022亚洲国产成人精品| 国产成人aa在线观看| av.在线天堂| 丝袜美腿在线中文| 亚洲最大成人中文| 18禁在线播放成人免费| 久久亚洲国产成人精品v| 国产男女内射视频| 免费人成在线观看视频色| 亚洲国产精品专区欧美| 久久精品国产亚洲av天美| 亚洲成人中文字幕在线播放| 一级av片app| 亚洲精品久久久久久婷婷小说| 国产精品熟女久久久久浪| 成年版毛片免费区| 亚洲国产精品成人综合色| 欧美+日韩+精品| 91久久精品国产一区二区三区| 亚洲人成网站高清观看| 国产成年人精品一区二区| 亚洲国产精品国产精品| 少妇丰满av| 一级av片app| 汤姆久久久久久久影院中文字幕| 亚洲精品成人av观看孕妇| 国产精品一区二区性色av| 久久久久九九精品影院| 熟女电影av网| 亚洲欧美清纯卡通| 美女脱内裤让男人舔精品视频| 午夜福利视频精品| 白带黄色成豆腐渣| 亚洲精品影视一区二区三区av| 看免费成人av毛片| 丝袜美腿在线中文| 久久久久久久久久久免费av| 卡戴珊不雅视频在线播放| 人人妻人人爽人人添夜夜欢视频 | 人体艺术视频欧美日本| 成人无遮挡网站| h日本视频在线播放| 亚洲欧美日韩卡通动漫| 少妇高潮的动态图| 亚洲国产精品成人综合色| 天天一区二区日本电影三级| av在线播放精品| 91狼人影院| 我要看日韩黄色一级片| 亚洲欧美成人精品一区二区| 在线观看国产h片| 午夜免费男女啪啪视频观看| 日韩强制内射视频| av在线亚洲专区| 综合色av麻豆| 亚洲国产精品999| av一本久久久久| 日韩电影二区| 丰满乱子伦码专区| 国产成人午夜福利电影在线观看| 成人漫画全彩无遮挡| kizo精华| 深爱激情五月婷婷| 国产成人freesex在线| 午夜亚洲福利在线播放| 天天一区二区日本电影三级| 91久久精品国产一区二区成人| 亚洲国产日韩一区二区| 国产精品久久久久久av不卡| 久久久久久久久久人人人人人人| 能在线免费看毛片的网站| 夫妻性生交免费视频一级片| 插阴视频在线观看视频| 身体一侧抽搐| 亚洲成人久久爱视频| 精品人妻偷拍中文字幕| 国产中年淑女户外野战色| 国产黄色免费在线视频| 亚洲成人一二三区av| 自拍偷自拍亚洲精品老妇| 日本黄色片子视频| 国产精品成人在线| 亚洲av不卡在线观看| 26uuu在线亚洲综合色| 女的被弄到高潮叫床怎么办| 久久精品久久久久久噜噜老黄| 国产91av在线免费观看| 久久久久久伊人网av| 老师上课跳d突然被开到最大视频| 七月丁香在线播放| 精品一区二区三卡| 欧美少妇被猛烈插入视频| 看黄色毛片网站| 久久ye,这里只有精品| 欧美老熟妇乱子伦牲交| 免费观看性生交大片5| 欧美激情久久久久久爽电影| 毛片一级片免费看久久久久| 国产精品久久久久久久电影| 国产片特级美女逼逼视频| 男的添女的下面高潮视频| 欧美xxxx性猛交bbbb| 婷婷色麻豆天堂久久| 久久99热6这里只有精品| 2021天堂中文幕一二区在线观| 免费看av在线观看网站| 免费播放大片免费观看视频在线观看| 国产日韩欧美亚洲二区| 国产亚洲91精品色在线| 国产精品国产av在线观看| 亚洲最大成人av| 国产毛片在线视频| 亚洲色图综合在线观看| 波野结衣二区三区在线| 热99国产精品久久久久久7| 边亲边吃奶的免费视频| 青春草国产在线视频| 欧美日韩视频精品一区| 久久精品国产亚洲av天美| av免费在线看不卡| 日韩一区二区视频免费看| 激情五月婷婷亚洲| 不卡视频在线观看欧美| 小蜜桃在线观看免费完整版高清| av在线天堂中文字幕| 久久久久久久久久成人| 中国国产av一级| 国产精品麻豆人妻色哟哟久久| 午夜福利网站1000一区二区三区| 久久久久久九九精品二区国产| 国产成人免费无遮挡视频| 午夜视频国产福利| xxx大片免费视频| 日日啪夜夜爽| 哪个播放器可以免费观看大片| 97精品久久久久久久久久精品| 九九在线视频观看精品| 一二三四中文在线观看免费高清| 91狼人影院| 丰满乱子伦码专区| 精品一区二区三卡| 99久久中文字幕三级久久日本| 成人鲁丝片一二三区免费| av在线蜜桃| 性插视频无遮挡在线免费观看| 国产伦理片在线播放av一区| 国产av不卡久久| 伊人久久国产一区二区| 少妇人妻一区二区三区视频| 51国产日韩欧美| 永久网站在线| 久久午夜福利片| 色5月婷婷丁香| 一级毛片电影观看| 亚洲三级黄色毛片| 色吧在线观看| 日本三级黄在线观看| 欧美日韩精品成人综合77777| 秋霞在线观看毛片| 国产永久视频网站| 2021少妇久久久久久久久久久| 九九在线视频观看精品| 性插视频无遮挡在线免费观看| 国产毛片a区久久久久| 少妇裸体淫交视频免费看高清| 在线播放无遮挡| 久久久国产一区二区| 国产精品麻豆人妻色哟哟久久| 亚洲国产精品专区欧美| 久久这里有精品视频免费| 欧美日韩亚洲高清精品| 亚洲精品日韩在线中文字幕| 国产淫片久久久久久久久| 在线天堂最新版资源| 久久午夜福利片| av福利片在线观看| 成人漫画全彩无遮挡| 在线播放无遮挡| 久久久精品94久久精品| 国产又色又爽无遮挡免| 亚洲国产欧美人成| 国产伦精品一区二区三区四那| 麻豆国产97在线/欧美| 久久久久久久亚洲中文字幕| 成人国产av品久久久| 免费看不卡的av| 少妇裸体淫交视频免费看高清| 国产男女超爽视频在线观看| 亚洲成人久久爱视频| 中文欧美无线码| 久久精品国产鲁丝片午夜精品| 国产av不卡久久| 伦精品一区二区三区| av国产久精品久网站免费入址| 91久久精品电影网| 欧美 日韩 精品 国产| 18禁裸乳无遮挡动漫免费视频 | 91久久精品国产一区二区成人| 成年女人看的毛片在线观看| 成年免费大片在线观看| 美女国产视频在线观看| 国产成人午夜福利电影在线观看| 欧美成人a在线观看| 国产高潮美女av| 国产精品伦人一区二区| av免费在线看不卡| 欧美老熟妇乱子伦牲交| 欧美日韩国产mv在线观看视频 | 久久ye,这里只有精品| 女人久久www免费人成看片| 国内精品宾馆在线| 深夜a级毛片| videos熟女内射| 成人特级av手机在线观看| 国产成人精品福利久久| 亚洲欧美中文字幕日韩二区| 国产一区二区在线观看日韩| 男人添女人高潮全过程视频| 久久精品久久精品一区二区三区| 亚洲欧美一区二区三区黑人 | 最近的中文字幕免费完整| 欧美日韩精品成人综合77777| 亚洲人成网站高清观看| 观看美女的网站| av.在线天堂| 麻豆精品久久久久久蜜桃| 国产精品久久久久久精品电影| 赤兔流量卡办理| 亚洲aⅴ乱码一区二区在线播放| 寂寞人妻少妇视频99o| 女人被狂操c到高潮| 一本色道久久久久久精品综合| 国产精品不卡视频一区二区| 一区二区三区乱码不卡18| 午夜福利在线在线| 亚洲人成网站在线播| 精品久久国产蜜桃| 久久久久久久久久久免费av| 亚洲人与动物交配视频| 菩萨蛮人人尽说江南好唐韦庄| 少妇的逼水好多| 亚洲精品影视一区二区三区av| 男女边摸边吃奶| 别揉我奶头 嗯啊视频| 成人特级av手机在线观看| 人人妻人人澡人人爽人人夜夜| 最近的中文字幕免费完整| 国产伦精品一区二区三区视频9| 午夜视频国产福利| 一区二区三区免费毛片| 欧美亚洲 丝袜 人妻 在线| 涩涩av久久男人的天堂| 欧美一区二区亚洲| 亚洲欧美一区二区三区黑人 | 韩国高清视频一区二区三区| 97超视频在线观看视频| 日韩国内少妇激情av| 久久国内精品自在自线图片| 免费观看a级毛片全部| 在线观看一区二区三区激情| 少妇裸体淫交视频免费看高清| 欧美xxⅹ黑人| 亚洲av不卡在线观看| 亚洲欧洲日产国产| 成年版毛片免费区| 国产中年淑女户外野战色| 99热网站在线观看| 伦精品一区二区三区| 亚洲av免费高清在线观看| 国产一区二区三区av在线| 精品久久久噜噜| 亚洲欧美日韩卡通动漫| 天天躁夜夜躁狠狠久久av| 亚洲av在线观看美女高潮| 在线观看国产h片| 精品熟女少妇av免费看| 18禁动态无遮挡网站| 亚洲国产欧美人成| 成人亚洲欧美一区二区av| 久久精品久久久久久久性| 亚洲国产欧美在线一区| 五月开心婷婷网| 亚洲欧美一区二区三区国产| 亚洲伊人久久精品综合| 国产伦在线观看视频一区| 老司机影院成人| 日韩不卡一区二区三区视频在线| 18+在线观看网站| 男插女下体视频免费在线播放| 久久国产乱子免费精品| 日韩一区二区视频免费看| videossex国产| 黄片无遮挡物在线观看| 国产探花在线观看一区二区| 大片免费播放器 马上看| 自拍偷自拍亚洲精品老妇| 免费少妇av软件| 免费观看无遮挡的男女| 国产黄色免费在线视频| 99热这里只有是精品50| 亚洲,一卡二卡三卡| 午夜日本视频在线| 日本一本二区三区精品| 亚洲aⅴ乱码一区二区在线播放| 国产黄色免费在线视频| 国产成人精品一,二区| 91久久精品国产一区二区成人| 国产黄色免费在线视频| 99热这里只有是精品50| 在线观看免费高清a一片| 亚洲真实伦在线观看| 欧美97在线视频| 六月丁香七月| 国产午夜精品久久久久久一区二区三区| 青春草视频在线免费观看| 好男人视频免费观看在线| 国产真实伦视频高清在线观看| 午夜免费观看性视频| 最后的刺客免费高清国语| 97人妻精品一区二区三区麻豆| 国产在线男女| 大香蕉97超碰在线| 少妇丰满av| 国产伦理片在线播放av一区| 亚洲国产色片| 黄色日韩在线| 一级av片app|