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

    圍墾對(duì)崇明東灘濕地全程氨氧化微生物的影響①

    2021-01-12 02:29:22宋怡軒張?chǎng)卫?/span>賈仲君張耀鴻
    土壤 2020年6期
    關(guān)鍵詞:拷貝數(shù)灘涂硝化

    宋怡軒,金 銳,張?chǎng)卫?,?潔,賈仲君,張耀鴻*

    圍墾對(duì)崇明東灘濕地全程氨氧化微生物的影響①

    宋怡軒1,2,金 銳1,張?chǎng)卫?,張 潔1,賈仲君2,張耀鴻1*

    (1南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心/江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室,南京 210044;2土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點(diǎn)實(shí)驗(yàn)室(中國科學(xué)院南京土壤研究所),南京 210008)

    全程氨氧化細(xì)菌(complete ammonia oxidizers,Comammox)的發(fā)現(xiàn)根本改變了學(xué)術(shù)界對(duì)硝化過程的認(rèn)識(shí),但其地理分異規(guī)律及對(duì)氮轉(zhuǎn)化過程的貢獻(xiàn)仍不清楚。本研究選擇長江口崇明東灘不同圍墾年限(0、27、51、86 a)稻田表層耕作土壤,采用好氧培養(yǎng)試驗(yàn)測(cè)定土壤硝化潛力;通過標(biāo)靶功能基因?qū)崟r(shí)熒光定量硝化微生物的數(shù)量變異特征,包括全程氨氧化細(xì)菌(Comammox)、氨氧化細(xì)菌(ammonia-oxidizing bacteria,AOB)和氨氧化古菌(ammonia-oxidizing archaea,AOA)。結(jié)果表明,與圍墾0 a的自然灘涂濕地相比,圍墾27、51、86 a的水稻土凈硝化速率從N 2.24 mg/(kg·d)分別增加至N 19.3、11.6和11.4 mg/(kg·d),增幅高達(dá)5.1 倍~ 8.7倍。AOA的豐度與圍墾年限顯著正相關(guān)。自然灘涂濕地中AOA和AOB的數(shù)量分別為0.34×107copies/g和1.14×107copies/g,圍墾86 a后增幅最高可達(dá)27.9倍。自然灘涂濕地中Comammox Clade A和Clade B的基因拷貝數(shù)高于圍墾稻田土壤,且Comammox Clade A隨著圍墾年限增加其豐度顯著增加。統(tǒng)計(jì)分析發(fā)現(xiàn),AOB與土壤硝化速率顯著正相關(guān),可能在圍墾水稻土氨氧化過程中發(fā)揮了重要作用;而Comammox Clade A和Clade B與土壤總有機(jī)碳(TOC)、銨含量(NH4+)呈顯著負(fù)相關(guān)關(guān)系,可能更適應(yīng)于營養(yǎng)貧瘠的灘涂自然濕地土壤。

    圍墾;全程氨氧化細(xì)菌;氨氧化

    氨氧化作用作為土壤氮素轉(zhuǎn)化的關(guān)鍵過程之一, 是硝化過程的限速步驟。一直以來,人們普遍認(rèn)為硝化作用由兩類不同的化能自養(yǎng)型細(xì)菌:氨氧化細(xì)菌(ammonia-oxidizing bacteria, AOB)和亞硝酸鹽氧化細(xì)菌(nitrite-oxidizing bacteria, NOB)完成的兩步氮循環(huán)過程[1-2]。2004年,化能自養(yǎng)型氨氧化古菌(ammonia- oxidizing archaea, AOA)被發(fā)現(xiàn),隨后的許多研究表明,這類微生物在土壤和水生系統(tǒng)的硝化過程中發(fā)揮著重要作用,豐富了驅(qū)動(dòng)好氧氨氧化過程的微生物類型。從熱力學(xué)和動(dòng)力學(xué)角度來看,相比于經(jīng)典分步硝化,獨(dú)立完成氨氧化和亞硝酸鹽氧化的全程氨氧化過程產(chǎn)生的能量更多,在產(chǎn)率較高而生長速率較低的條件下具有更大優(yōu)勢(shì)[3-4]。2015年,Daims等[3]和van Kessel等[5]發(fā)現(xiàn)了硝化螺菌屬(spp.)中有微生物能全程完成硝化作用,即一步將氨氮轉(zhuǎn)化為硝態(tài)氮,這類微生物統(tǒng)稱為全程氨氧化細(xì)菌(complete ammonia oxidizers, Comammox)。該細(xì)菌的發(fā)現(xiàn)使人們對(duì)氮循環(huán)過程有了一個(gè)全新的認(rèn)識(shí)。通過分析宏基因組學(xué)數(shù)據(jù)發(fā)現(xiàn)Comammox細(xì)菌與硝化螺旋菌屬spp.中其他微生物生境相似,廣泛分布于土壤、淡水濕地、地?zé)釡厝⑽鬯幚韽S和飲用水處理系統(tǒng)中[3, 5-6]。然而,Comammox在濱海濕地生態(tài)系統(tǒng)中的分布特征和生態(tài)學(xué)重要性卻鮮有報(bào)道。

    崇明東灘位于長江口崇明島的最東端,是長江口規(guī)模最大、發(fā)育最完善的河口型潮汐灘涂濕地。入海徑流攜帶的泥沙在河口附近沉積,塑造了廣闊的潮灘濕地。隨著經(jīng)濟(jì)的發(fā)展,人類對(duì)土地的需求越來越迫切,灘涂圍墾改農(nóng)田成為解決土地問題的有效途徑。一方面,圍墾區(qū)的土壤來自于灘涂淤積的泥沙,土壤發(fā)育時(shí)間相對(duì)較短,其土壤本底較為均一;另一方面,圍墾后的濕地不再受潮汐影響,與近海間的物質(zhì)能量交換基本消失,特別是水稻種植、翻耕、施肥等農(nóng)業(yè)管理措施,使圍墾區(qū)生境發(fā)生顯著變化,形成了獨(dú)特的濕地類型:圍墾區(qū)稻田濕地,特別是大量的農(nóng)業(yè)活性氮進(jìn)入稻田生態(tài)系統(tǒng)后,使其成為氮循環(huán)的熱區(qū)。因此,開展濱海河口濕地土壤氨氧化過程及其相關(guān)功能微生物的研究,有助于深入認(rèn)識(shí)圍墾造田和農(nóng)業(yè)氮輸入對(duì)氮關(guān)鍵轉(zhuǎn)化過程的作用機(jī)理,具有重要的理論意義與科學(xué)價(jià)值。

    基于此,本文以典型濱海河口濕地——崇明東灘為研究區(qū),采用空間代替時(shí)間的研究方法,研究圍墾造田條件下濱海濕地土壤的好氧氨氧化過程,探討其空間變異特征及主要影響因素;通過微生物分子生物學(xué)研究方法,分析濕地土壤中全程氨氧化微生物(Clade A 和 Clade B)和經(jīng)典氨氧化微生物(AOA 和 AOB)的豐度變化,評(píng)估這些微生物功能基因的生態(tài)位分化及其對(duì)硝化作用的相對(duì)貢獻(xiàn),以期為長江口濕地環(huán)境保護(hù)與生態(tài)治理提供重要的參考依據(jù)。

    1 材料和方法

    1.1 土壤樣品采集

    采樣地點(diǎn)為上海市崇明島(121°09′ ~ 121°54′E,31°27′ ~ 31°51′N),屬于典型的亞熱帶季風(fēng)氣候,終年溫?zé)?,降水充沛,年均?5.3 ℃,年降水量1 003.7 mm[7]。本研究選取東灘濕地保護(hù)區(qū)灘涂土壤為圍墾 0 a土壤作為對(duì)照,其他3個(gè)水稻土的圍墾年限分別為27、51、86 a。圍墾植稻不同年限的實(shí)驗(yàn)樣點(diǎn)根據(jù)Cui等[8]的參考文獻(xiàn)選取。為了使采樣點(diǎn)能更準(zhǔn)確反映不同的土壤發(fā)育年限,更具代表性,在每個(gè)樣區(qū)內(nèi)以S形設(shè)置6個(gè)采樣點(diǎn),各采樣點(diǎn)間距約為10 ~ 15 m,用土鉆取表層0 ~ 20 cm鮮土,并將該6個(gè)采樣點(diǎn)的土壤均勻混合為1個(gè)混合樣本,于保溫箱加冰袋保存運(yùn)回實(shí)驗(yàn)室冷凍保存?zhèn)溆谩?/p>

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

    土壤理化性質(zhì)的分析測(cè)定主要參考《土壤農(nóng)業(yè)化學(xué)分析方法》[9]。土壤全氮含量采用半微量凱氏定量法測(cè)定。土壤總有機(jī)碳采用濃硫酸–重鉻酸鉀消煮,硫酸亞鐵滴定法測(cè)定。土壤NH4+-N和NO– 3-N用2 mol/L KCl溶液浸提后,采用AA3流動(dòng)分析儀測(cè)定;亞硝態(tài)氮采用重氮反應(yīng)–比色法測(cè)定。土壤pH采用水土比為2.5︰1提取水溶液后,用數(shù)字酸度計(jì)(PHSJ- 6L,上海儀電科學(xué)儀器股份有限公司)測(cè)定;土壤電導(dǎo)率用電導(dǎo)法測(cè)定,采用水土為5︰1提取水溶液后用電導(dǎo)率儀(DDSJ-319L,上海儀電科學(xué)儀器股份有限公司)測(cè)定。

    硝化速率的測(cè)定采用好氣培養(yǎng)法[10],取沉積物土樣10 g(干基),加15 ml去離子水,將土樣分散均勻,形成懸濁液。稱取總重量并記錄。25 ℃避光預(yù)培養(yǎng)7 d以便激活土壤微生物。預(yù)培養(yǎng)結(jié)束后,將每個(gè)土樣分為兩組:第一組的3個(gè)錐形瓶,直接測(cè)定其中的NO– 3-N含量(0)。第二組的6個(gè)錐形瓶中,用移液槍加入(NH4)2SO4溶液使N最終濃度達(dá)到200 mg/L,再根據(jù)重量法補(bǔ)加去離子水,使之去離子水總量為20 ml。然后在恒溫培養(yǎng)箱中25 ℃避光培養(yǎng)7 d后,取出錐形瓶,測(cè)定其中的NO– 3-N 含量(= 7)。

    硝化速率的計(jì)算公式為:

    式中:為硝化速率(N, mg/(kg·d));(NO– 3)和(NO– 3)0分別為7 d和0 d 時(shí)NO– 3-N含量。

    1.3 土壤DNA提取

    利用Fast DNA Spin kit for soil提取試劑盒 (MPBiomedicals, USA) 提取土壤樣品中的總DNA,取0.5 g土壤樣品,按試劑盒提供的試劑和步驟提取DNA,F(xiàn)astPrep?FP120 細(xì)胞破碎儀以速度 6 m/s,時(shí)間 40 s進(jìn)行細(xì)胞破碎。取部分DNA提取液用分光光度計(jì)(NanoDrop -1000 UV-Vis)測(cè)定DNA濃度和純度( OD260 /OD280和OD260 /OD230)。土壤DNA保存于–80 ℃冰箱待用。

    1.4 實(shí)時(shí)熒光定量PCR

    在C1000TMReal-Time System擴(kuò)增儀上進(jìn)行實(shí)時(shí)熒光定量PCR擴(kuò)增。測(cè)定氨氧化古菌(AOA)、氨氧化細(xì)菌(AOB)和全程氨氧化細(xì)菌(Comammox Clade A與總Comammox)的基因拷貝數(shù)。AOB、AOA、總Comammox和Comammox Clade A的基因定量 PCR 所用擴(kuò)增引物和擴(kuò)增條件見表1[6,11-14]。反應(yīng)體系均為20 μl,包括DNA樣品1 μl、Taq DNA聚合酶10 μl、前后引物各0.5 μl、無菌水8 μl。

    表 1 AOA, AOB, Comammox Clade A和總Comammox 的amoA基因引物序列及其擴(kuò)增條件

    提取AOA、AOB、Comammox Clade A和總Comammox的基因的重組質(zhì)粒,并通過測(cè)序驗(yàn)證,再用分光光度計(jì)( NanoDrop ND-1000 UV-Vis)測(cè)定質(zhì)粒濃度,并用無菌水將質(zhì)粒稀釋6 ~ 8個(gè)梯度,用于制作定量PCR的標(biāo)準(zhǔn)曲線,根據(jù)標(biāo)準(zhǔn)曲線計(jì)算出目的基因的拷貝數(shù)。Comammox Clade B基因的拷貝數(shù)為總Comammox 與Comammox Clade A的差值。

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

    用IBM Statistics SPSS 16.0軟件進(jìn)行統(tǒng)計(jì)分析,通過單因素方差分析(one-way ANOVA)和Pearson分析土壤理化性質(zhì)、氨氧化微生物豐度的差異性以及相關(guān)性檢驗(yàn),顯著性水平= 0.05。

    2 結(jié)果與分析

    2.1 不同圍墾年限濕地土壤的理化性質(zhì)變化特征

    崇明島濕地不同圍墾年限土壤理化性質(zhì)如表2所示。長期植稻顯著提升了水稻土無機(jī)氮含量,降低了硫酸鹽含量。圍墾稻田土壤NO– 3含量為11.79 ~ 18.67 mg/kg,是未圍墾灘涂濕地土壤的1.8倍 ~ 2.9倍;圍墾稻田土壤NH4+含量為14.48 ~ 23.37 mg/kg,比未圍墾灘涂濕地土壤高出1.1倍 ~ 2.7倍;且圍墾51 a稻田土壤的NH4+和NO– 3含量均顯著高于其他3個(gè)采樣點(diǎn)。土壤NO– 2含量在不同采樣點(diǎn)之間沒有顯著差異。土壤SO2– 4 濃度與EC值的變化趨勢(shì)一致,也隨著圍墾年限增加而明顯減小,且未圍墾灘涂濕地土壤顯著高于圍墾稻田。圍墾稻田的土壤總有機(jī)碳(TOC)介于15.29 ~ 18.34 g/kg,顯著高于未圍墾的自然灘涂濕地,且隨圍墾年限增加有顯著增加的趨勢(shì)。土壤全氮(TN)含量介于0.84 ~ 1.12 g/kg,4個(gè)采樣點(diǎn)之間沒有顯著差異。未圍墾灘涂濕地土壤的pH為7.86,顯著高于圍墾稻田。土壤電導(dǎo)率隨著圍墾年限的增加顯著減小,且未圍墾灘涂濕地土壤的EC值是圍墾稻田的3.2倍 ~ 6.9倍。

    表2 不同圍墾年限土壤的理化性質(zhì)

    注:同一列數(shù)據(jù)小寫字母不同表示不同圍墾年限處理間差異達(dá)<0.05顯著水平,下表同。

    2.2 不同圍墾年限濕地土壤的硝化速率

    如圖1所示,圍墾0 a的自然灘涂濕地硝化速率最低,僅為N 2.22 mg/(kg·d),圍墾27、51、86 a后水稻土硝化速率分別為N 19.3、11.6和11.4 mg/(kg·d),增幅高達(dá)8.7倍、5.2倍和5.1倍。其中,圍墾27 a稻田土壤的硝化速率顯著高于圍墾51 a和86 a土壤,分別高出66% 和67%。

    圖1 不同圍墾年限土壤的硝化速率

    2.3 氨氧化細(xì)菌和古菌的豐度特征

    由圖2可知,圍墾0、27、51、86 a的土壤AOA的基因拷貝數(shù)(以每克干土計(jì))分別為0.34×107、1.10×107、4.07×107、9.95×107copies/g,隨著圍墾年限增加而顯著增大(< 0.05)。圍墾27、51、86 a的稻田土壤AOB的基因拷貝數(shù)(以每克干土計(jì))為8.43×107、6.96×107、6.69×107copies/g,分別是圍墾0 a灘涂濕地土壤的7.4倍、6.1倍和5.9倍。相關(guān)性分析發(fā)現(xiàn),4個(gè)采樣點(diǎn)的AOB基因拷貝數(shù)與土壤硝化速率呈顯著正相關(guān)(= 0.96*)。

    2.4 全程氨氧化細(xì)菌的豐度特征

    近年來國外研究者發(fā)現(xiàn)了具備單獨(dú)將NH4+氧化為NO– 3的可純培養(yǎng)細(xì)菌,均屬于亞硝酸鹽氧化細(xì)菌中的硝化螺菌屬,被稱為全程氨氧化細(xì)菌。由圖3可知,圍墾0、27、51、86 a土壤的Comammox Clade A 拷貝數(shù)分別為1.04×107、0.14×107、0.64×107、1.00×107copies/g。其中,圍墾稻田土壤的Clade A 拷貝數(shù)隨著圍墾年限增加而顯著增大,但均低于圍墾0 a的灘涂濕地土壤。Comammox Clade B在4個(gè)采樣點(diǎn)中的豐度特征與Clade A變化特征有所不同,其拷貝數(shù)分別為5.73×106、1.78×106、0.48× 106、2.04×106copies/g;其中,圍墾稻田土壤的Clade B 拷貝數(shù)均顯著低于圍墾0 a的灘涂濕地土壤。

    圖2 不同圍墾年限濕地土壤AOA和AOB amoA功能基因拷貝數(shù)

    2.5 不同氨氧化微生物之間的豐度比值

    由表3可見,崇明島圍墾年限0、27、51、86 a的土壤中,不同氨氧化微生物之間的豐度比值A(chǔ)OA/AOB介于0.13 ~ 1.49,表明這兩類微生物的功能基因拷貝數(shù)大致在一個(gè)數(shù)量級(jí)水平;Clade A/Clade B介于0.77 ~ 13.35,表明Clade A的功能基因與Clade B的變化范圍較大;AOA、AOB與Clade A 的功能基因拷貝數(shù)的比值介于0.33 ~ 61.43。AOA、AOB與Clade B 的功能基因拷貝數(shù)的比值介于0.60 ~ 146.18。

    圖3 不同圍墾年限Comammox amoA功能基因拷貝數(shù)

    表 3 不同氨氧化微生物的豐度比值

    2.6 氨氧化微生物豐度與土壤理化性質(zhì)的相關(guān)性

    對(duì)圍墾0、27、51、86 a的濕地土壤理化性質(zhì)與氨氧化微生物豐度進(jìn)行了相關(guān)性分析。由表4可見,AOB基因拷貝數(shù)與NH4+、NO– 3呈正相關(guān),且與NO– 3含量達(dá)顯著水平;而Clade A和Clade B均與NH4+、NO– 3呈負(fù)相關(guān),且與NH4+含量都達(dá)到顯著水平。推測(cè)全程氨氧化微生物(Clade A、Clade B)更適應(yīng)于營養(yǎng)貧瘠條件,而AOB可能在活性氮源充足條件下主導(dǎo)硝化作用。TOC、NO– 2與AOA、Clade A、Clade B之間呈負(fù)相關(guān)關(guān)系,而與AOB呈正相關(guān),且TOC與AOB、Clade A、Clade B的相關(guān)性達(dá)到顯著水平。4類氨氧化微生物的基因拷貝數(shù)均與pH呈正相關(guān),且與AOA的相關(guān)性顯著。SO2– 4、EC與AOB豐度呈顯著負(fù)相關(guān),與Clade B表現(xiàn)為顯著正相關(guān)。

    表 4 氨氧化微生物豐度與土壤基本理化性質(zhì)的相關(guān)性

    注:* 表示相關(guān)性達(dá)< 0.05顯著水平。

    3 討論

    本試驗(yàn)以長江口崇明島東灘濕地為對(duì)象來研究典型氨氧化微生物(AOB、AOA)與全程氨氧化細(xì)菌(Comammox)的分布及豐度特征,對(duì)了解長江口濕地生態(tài)系統(tǒng)中硝化作用微生物分布及數(shù)量特征具有重要意義。本研究中,不同圍墾年限(0、27、51、86 a)的濕地土壤中均檢測(cè)到全程氨氧化細(xì)菌,證實(shí)了全程氨氧化細(xì)菌存在于長江口崇明島圍墾稻田和自然灘涂濕地中。國外研究者基于宏基因組數(shù)據(jù)庫比對(duì)功能基因,推測(cè)全程氨氧化細(xì)菌廣泛存在于稻田土壤[3,5]、濕地[3]、淡水生態(tài)系統(tǒng)、工程環(huán)境[6]等。Wang等[2]發(fā)現(xiàn)不同地區(qū)自來水處理廠的樣本中Comammox豐度遠(yuǎn)高于AOA和AOB的的豐度,達(dá)到2倍 ~ 5倍。不同圍墾年限的AOA、AOB與全程氨氧化細(xì)菌(Clade A或Clade B)的功能基因拷貝數(shù)的比值遠(yuǎn)遠(yuǎn)大于1(表4),表明在此生境中經(jīng)典氨氧化微生物比全程氨氧化細(xì)菌在數(shù)量上更具優(yōu)勢(shì)。本研究的采樣區(qū)域?yàn)楹j懡诲e(cuò)帶,其全程氨氧化細(xì)菌兩個(gè)分支Clade A與Clade B 的功能基因的檢測(cè)結(jié)果,對(duì)于揭示長江入海河口區(qū)域中是否存在全程氨氧化細(xì)菌可提供重要的參考數(shù)據(jù)。在本研究中,Clade A基因拷貝數(shù)與AOA、AOB的數(shù)量級(jí)相同,均為107copies/g,且Clade A與Clade B比值范圍為0.8 ~ 13.4,說明在長江口圍墾稻田濕地與自然灘涂濕地中全程氨氧化細(xì)菌主要以Comammox Clade A居多。萬琪慧等[15]在紫色水稻土中也檢測(cè)到Comammox Clade A和Clade B,且Clade A/Clade B比值范圍為16 ~ 34,表明兩個(gè)分支以Clade A為主,本試驗(yàn)結(jié)果與此一致。新近,國外研究者在意大利農(nóng)田土壤中檢測(cè)到了兩分支的Comammox豐度值均高于104genecopies/ng(以DNA計(jì)),其中Comammox Clade A的基因拷貝數(shù)比Comam-mox Clade B高出2 ~ 3個(gè)數(shù)量級(jí)[6],暗示了這些生境中Comammox Clade A在氮硝化過程中可能發(fā)揮著重要的作用。

    通常認(rèn)為,氮素含量豐富,且中性至堿性的土壤有利于AOB生長繁殖,施加氮肥的稻田土壤往往以AOB起主導(dǎo)作用進(jìn)行硝化過程。鐘文輝等人[16]研究發(fā)現(xiàn)長期施用氮肥提高氨氧化細(xì)菌數(shù)量可能是施氮土壤硝化率和硝化勢(shì)較高的原因之一。本試驗(yàn)中,相關(guān)性分析發(fā)現(xiàn)AOB與土壤TOC、NO– 3呈顯著正相關(guān);隨著圍墾年限增加,硝化速率和AOB的基因拷貝數(shù)的變化趨勢(shì)一致,兩者呈顯著正相關(guān),表明AOB對(duì)氨氧化過程起主導(dǎo)作用。黃蓉等[17]在提高酸性紅壤的pH后,發(fā)現(xiàn)AOB提高2個(gè)數(shù)量級(jí),而硝化速率也顯著增加,本試驗(yàn)結(jié)果與此一致。Comammox的基因拷貝數(shù)均明顯低于AOB,且與硝化速率呈顯著負(fù)相關(guān),暗示了Comammox對(duì)硝化速率的相對(duì)貢獻(xiàn)較小。兩類全程氨氧化細(xì)菌即Comammox Clade A與Comammox Clade B有相似規(guī)律,未開墾灘涂濕地的基因拷貝數(shù)普遍高于圍墾稻田土壤,暗示這類微生物在灘涂自然濕地的氮轉(zhuǎn)化過程中可能起著重要的作用。

    大量的研究結(jié)果發(fā)現(xiàn),AOA比AOB對(duì)氨具有更高的親和力,因此當(dāng)氨濃度很低時(shí),AOA對(duì)氨的競(jìng)爭(zhēng)遠(yuǎn)勝于AOB。Martens等[18]發(fā)現(xiàn)在海洋中的AOA可以通過調(diào)整細(xì)胞機(jī)制在低于10 nmol/L的底物濃度環(huán)境中進(jìn)行生長。最近的研究發(fā)現(xiàn)Comammox對(duì)氨的親和力大于AOA。本研究中,Comammox分支A和B(Comammox Clade A與Comammox Clade B)的A基因拷貝數(shù)與TOC、NH4+濃度均呈顯著負(fù)相關(guān)(表4),推測(cè)在沒有圍墾的自然灘涂濕地中,TOC、活性氮含量少的貧營養(yǎng)環(huán)境可能更有利于Comam-mox的生長繁殖。

    綜上所述,長江口濕地中檢測(cè)到ComammoxA基因,且與AOA和AOB共同生存,可能存在對(duì)底物和碳源的競(jìng)爭(zhēng)。如要精確計(jì)算三者對(duì)表觀硝化速率的相對(duì)貢獻(xiàn),需要進(jìn)一步考察其RNA反轉(zhuǎn)錄、蛋白質(zhì)等代謝水平[19],進(jìn)一步從宏基因組、轉(zhuǎn)錄組、蛋白質(zhì)組、代謝組等方面探究其生理生化特性,從而可以準(zhǔn)確評(píng)估其實(shí)際對(duì)硝化作用的貢獻(xiàn)率。

    4 結(jié)論

    1)本試驗(yàn)在4種圍墾年限的稻田土壤中檢測(cè)到全程氨氧化細(xì)菌(Comammox)A基因,表明Comammox廣泛分布在該圍墾地區(qū),且以分支Clade A為主。

    2)未開墾灘涂濕地中Comammox Clade A 和Clade BA基因普遍高于圍墾稻田土壤。

    3)半程氨氧化細(xì)菌(AOB)與硝化速率隨圍墾年限的變化一致,可能對(duì)硝化速率貢獻(xiàn)最大。

    [1] Nunes-Alves C. Do it yourself nitrification[J]. Nature Reviews Microbiology, 2016, 14(2): 61.

    [2] Wang Y L, Ma L P, Mao Y P, et al. Comammox in drinking water systems[J]. Water Research, 2017, 116: 332–341.

    [3] Daims H, Lebedeva E V, Pjevac P, et al. Complete nitrification by Nitrospira bacteria[J]. Nature, 2015, 528(7583): 504–509.

    [4] Costa E, Pérez J, Kreft J U. Why is metabolic labour divided in nitrification?[J]. Trends in Microbiology, 2006, 14(5): 213–219.

    [5] van Kessel M A H J, Speth D R, Albertsen M, et al. Complete nitrification by a single microorganism[J]. Nature, 2015, 528(7583): 555–559.

    [6] Pjevac P, Schauberger C, Poghosyan L, et al. AmoA- targeted polymerase chain reaction primers for the specific detection and quantification of comammox Nitrospira in the environment[J]. bioRxiv, 2017, DOI:10.1101/096891.

    [7] Compilation Committee of Chongming Annals. Chongming county annals[M]. Shanghai: Shanghai People’s Publishing House, 1989: 42.

    [8] Cui J, Liu C, Li Z L, et al. Long-term changes in topsoil chemical properties under centuries of cultivation after reclamation of coastal wetlands in the Yangtze Estuary, China[J]. Soil and Tillage Research, 2012, 123: 50–60.

    [9] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國農(nóng)業(yè)科技出版社, 2000.

    [10] Chu H Y, Fujii T, Morimoto S, et al. Community structure of ammonia-oxidizing bacteria under long-term application of mineral fertilizer and organic manure in a sandy loam soil[J]. Applied and Environmental Microbiology, 2007, 73(2): 485–491.

    [11] Beman J M, Francis C A. Diversity of ammonia-oxidizing Archaea and bacteria in the sediments of a hypernutrified subtropical estuary: Bahía del Tóbari, Mexico[J]. Applied and Environmental Microbiology, 2006, 72(12): 7767– 7777.

    [12] Rotthauwe J H, Witzel K P, Liesack W. The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia- oxidizing populations[J]. Applied and Environmental Microbiology, 1997, 63(12): 4704–4712.

    [13] Yu C D, Hou L J, Zheng Y L, et al. Evidence for complete nitrification in enrichment culture of tidal sediments and diversity analysis of clade a comammox Nitrospira in natural environments[J]. Applied Microbiology and Biotechnology, 2018, 102(21): 9363–9377.

    [14] Wang M Y, Huang G H, Zhao Z R, et al. Newly designed primer pair revealed dominant and diverse comammox amoA gene in full-scale wastewater treatment plants[J]. Bioresource Technology, 2018, 270: 580–587.

    [15] 萬琪慧, 王書玲, 趙偉燁, 等. 重慶紫色水稻土中“全程”和“半程”氨氧化微生物的垂直分異[J]. 微生物學(xué)報(bào), 2019, 59(2): 291-302.

    [16] 鐘文輝, 蔡祖聰, 尹力初, 等. 種植水稻和長期施用無機(jī)肥對(duì)紅壤氨氧化細(xì)菌多樣性和硝化作用的影響[J]. 土壤學(xué)報(bào), 2008, 45(1): 105–111.

    [17] 黃蓉, 張金波, 鐘文輝, 等. 土地利用方式對(duì)萬木林土壤氨氧化微生物豐度的影響[J]. 土壤, 2012, 44(4): 581–587.

    [18] Martens-Habbena W, Berube P M, Urakawa H, et al. Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria[J]. Nature, 2009, 461(7266): 976–979.

    [19] Liu S, Shen L D, Lou L P, et al. Spatial distribution and factors shaping the niche segregation of ammonia- oxidizing microorganisms in the Qiantang River, China[J]. Applied and Environmental Microbiology, 2013, 79(13): 4065–4071.

    Effect of Reclamation on Abundance of Comammox Microorganism in Chongming Wetland

    SONG Yixuan1,2, JIN Rui1, ZHANG Xinlei1, ZHANG Jie1, JIA Zhongjun2, ZHANG Yaohong1*

    (1 Nanjing University of Information Science and Technology, Collaborative Innovation Center for Forecast and Evaluation of Meteorological Disasters/Jiangsu Provincial Key Laboratory of Agricultural Meteorology, Nanjing 210044, China; 2 State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China)

    The discovery of comammox broadens our understanding of the nitrification process, but little is known about its natural geographical distribution and its contribution to nitrogen conversion process. In this study, the surface soils of the paddy fields with different reclamation years (0, 27, 51 and 86 a) were selected in the Chongming wetland of the Yangtze River estuary, soil nitrification rates and its influencing factors were determined by aerobic incubation, and the target functional genewas analyzed by real-time fluorescent quantitative PCR to analyze the quantitative variation characteristics of Comammox, ammoxidox bacteria (AOB) and ammonia archaea (AOA). The results showed that the net nitrification rates of paddy soils with 40, 75 and 120 a of reclamation were N 19.26, 11.63 and 11.43 mg/(kg·d), respectively, which were significantly 5.1-8.7 times higher than that of natural tidal wetland (<0.05). Thegene copy number of soil ammonia archaea AOA was between 0.34×107and 9.95×107copies/g, which increased significantly with the increase of reclamation year. Thegene copy number of soil ammonia bacteria AOB ranged from 1.14×107to 8.43×107copies/g, which was significantly positively correlated with soil nitrification rate (= 0.96*,<0.05), indicating that AOB might play a dominated role in nitrification process. The Comammox Clade A and Clade B copy number of the soils in the paddy fields (27, 51 and 86 a) were lower than that of tidal flat wetland soil (no reclamation), and Clade A increased significantly with the increase of the reclamation year. Comammox Clade A/Clade B was significantly negatively correlated with soil total organic carbon (TOC) and ammonium content (NH4+), indicating that the ammonia-oxidizing microbes (Clade A, Clade B) are more suitable for nutritional infertility. In summary, Comammox bacteria is widely present in the coastal wetland soils of the Yangtze River estuary, which possibly may play an important role in nitrogen transformation process of tidal flat natural wetlands.

    Reclamation; Comammox; Ammonia oxidation

    S154

    A

    10.13758/j.cnki.tr.2020.06.014

    宋怡軒, 金銳, 張?chǎng)卫? 等. 圍墾對(duì)崇明東灘濕地全程氨氧化微生物的影響. 土壤, 2020, 52(6): 1196–1202.

    國家自然科學(xué)基金項(xiàng)目(41671247, 41103039),江蘇省自然科學(xué)基金項(xiàng)目( BK20171455) 和江蘇省大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃項(xiàng)目(201810300094X, 201710300018)資助。

    (yhzhang@nuist.edu.cn)

    宋怡軒(1995—),女,四川南充人,碩士研究生,主要從事農(nóng)業(yè)應(yīng)對(duì)氣候變化研究。E-mail: 2287784164@qq.com

    猜你喜歡
    拷貝數(shù)灘涂硝化
    中國發(fā)起外來灘涂草治理戰(zhàn)
    英語世界(2023年10期)2023-11-17 09:18:22
    霞浦灘涂
    線粒體DNA拷貝數(shù)變異機(jī)制及疾病預(yù)測(cè)價(jià)值分析
    胎兒染色體組拷貝數(shù)變異與產(chǎn)前超聲異常的相關(guān)性分析
    MBBR中進(jìn)水有機(jī)負(fù)荷對(duì)短程硝化反硝化的影響
    海邊灘涂軟淤地基擋潮閘基礎(chǔ)設(shè)計(jì)
    厭氧氨氧化與反硝化耦合脫氮除碳研究Ⅰ:
    DNA序列拷貝數(shù)變化決定黃瓜性別
    沿海灘涂圍墾施工的技術(shù)要點(diǎn)
    精品(2015年9期)2015-01-23 01:35:32
    海水反硝化和厭氧氨氧化速率同步測(cè)定的15N示蹤法及其應(yīng)用
    少妇的丰满在线观看| 亚洲国产欧美网| 日日摸夜夜添夜夜添小说| 成人国产综合亚洲| 白带黄色成豆腐渣| 欧美日韩国产亚洲二区| 国产精品亚洲美女久久久| 黑人欧美特级aaaaaa片| svipshipincom国产片| 内地一区二区视频在线| 亚洲久久久久久中文字幕| 欧美成人免费av一区二区三区| 免费在线观看成人毛片| 午夜激情欧美在线| 午夜精品一区二区三区免费看| 免费看美女性在线毛片视频| 国产国拍精品亚洲av在线观看 | 午夜亚洲福利在线播放| 国产精品一区二区三区四区久久| 99精品欧美一区二区三区四区| 一a级毛片在线观看| 99久久精品热视频| 午夜免费观看网址| 久久香蕉国产精品| 真人一进一出gif抽搐免费| 级片在线观看| 欧美日韩综合久久久久久 | 欧美日韩瑟瑟在线播放| 国产精品香港三级国产av潘金莲| 69人妻影院| 亚洲人成伊人成综合网2020| 国产伦在线观看视频一区| 亚洲精品亚洲一区二区| 18禁在线播放成人免费| 99久国产av精品| 好看av亚洲va欧美ⅴa在| 亚洲精品在线观看二区| 天天添夜夜摸| www国产在线视频色| 欧美日韩国产亚洲二区| av欧美777| 亚洲av熟女| 国产v大片淫在线免费观看| 亚洲一区二区三区不卡视频| 国产成人系列免费观看| 亚洲av熟女| 成人亚洲精品av一区二区| 一个人看视频在线观看www免费 | 男女那种视频在线观看| 人妻夜夜爽99麻豆av| 国产欧美日韩精品亚洲av| 美女cb高潮喷水在线观看| 99久久成人亚洲精品观看| 国产黄片美女视频| 久久中文看片网| 国产亚洲av嫩草精品影院| 中文在线观看免费www的网站| 欧美日韩综合久久久久久 | 精品国产亚洲在线| 有码 亚洲区| 麻豆国产av国片精品| 怎么达到女性高潮| av中文乱码字幕在线| 国产精品爽爽va在线观看网站| 国产美女午夜福利| 一级毛片高清免费大全| 国产午夜精品久久久久久一区二区三区 | 午夜老司机福利剧场| svipshipincom国产片| 岛国在线观看网站| 国产高潮美女av| 欧美+亚洲+日韩+国产| 午夜影院日韩av| 午夜福利欧美成人| 国产三级在线视频| 亚洲精品成人久久久久久| 亚洲精品美女久久久久99蜜臀| 嫩草影院精品99| 久久亚洲真实| 国内毛片毛片毛片毛片毛片| 18美女黄网站色大片免费观看| 亚洲人成伊人成综合网2020| 99久久九九国产精品国产免费| 色噜噜av男人的天堂激情| 免费人成视频x8x8入口观看| 亚洲国产精品久久男人天堂| 不卡一级毛片| 亚洲无线在线观看| e午夜精品久久久久久久| 老司机深夜福利视频在线观看| 亚洲成人精品中文字幕电影| 色播亚洲综合网| 国产成人啪精品午夜网站| 国产精品1区2区在线观看.| 国产精品久久久人人做人人爽| 18禁裸乳无遮挡免费网站照片| 久久国产精品人妻蜜桃| 国产三级中文精品| 久久6这里有精品| 蜜桃久久精品国产亚洲av| 一个人看的www免费观看视频| 久久99热这里只有精品18| 亚洲五月天丁香| 中文字幕人妻熟人妻熟丝袜美 | 麻豆成人午夜福利视频| 亚洲国产日韩欧美精品在线观看 | 国产色爽女视频免费观看| 精品一区二区三区视频在线观看免费| 黄色视频,在线免费观看| 免费看日本二区| 色老头精品视频在线观看| 久久精品国产亚洲av香蕉五月| 亚洲五月婷婷丁香| 亚洲成av人片在线播放无| 午夜福利18| 18禁美女被吸乳视频| 精品99又大又爽又粗少妇毛片 | 99国产精品一区二区三区| 啦啦啦免费观看视频1| 男女视频在线观看网站免费| 精品国内亚洲2022精品成人| 国产一区二区在线观看日韩 | 国产精品乱码一区二三区的特点| 国产私拍福利视频在线观看| 亚洲人成伊人成综合网2020| 亚洲国产精品合色在线| 99热只有精品国产| 亚洲人与动物交配视频| 一本综合久久免费| 色在线成人网| 国产成人影院久久av| 国产成人av激情在线播放| 91字幕亚洲| 国产国拍精品亚洲av在线观看 | 毛片女人毛片| 99久久精品国产亚洲精品| 精品日产1卡2卡| 亚洲精品乱码久久久v下载方式 | 高清日韩中文字幕在线| 免费大片18禁| 国产美女午夜福利| e午夜精品久久久久久久| 最近最新中文字幕大全电影3| 亚洲最大成人中文| 中文字幕人成人乱码亚洲影| 亚洲无线在线观看| 免费观看的影片在线观看| av片东京热男人的天堂| 99热这里只有是精品50| 级片在线观看| 欧美高清成人免费视频www| 精品一区二区三区视频在线观看免费| 1000部很黄的大片| 精品国产三级普通话版| 19禁男女啪啪无遮挡网站| 男人的好看免费观看在线视频| 国产色爽女视频免费观看| 精品人妻偷拍中文字幕| 搡老妇女老女人老熟妇| 亚洲av成人精品一区久久| 又粗又爽又猛毛片免费看| 叶爱在线成人免费视频播放| 在线播放国产精品三级| 日本五十路高清| 亚洲av第一区精品v没综合| 欧美另类亚洲清纯唯美| 成人国产综合亚洲| 国产精品av视频在线免费观看| 国产三级黄色录像| 精品一区二区三区视频在线 | 99国产精品一区二区蜜桃av| 真人做人爱边吃奶动态| 国产伦精品一区二区三区视频9 | 精品国产亚洲在线| 俄罗斯特黄特色一大片| 久久久久久久亚洲中文字幕 | 在线观看美女被高潮喷水网站 | 亚洲成av人片在线播放无| 亚洲片人在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲最大成人中文| 亚洲无线在线观看| 国产在线精品亚洲第一网站| 久久欧美精品欧美久久欧美| 午夜免费激情av| 成人精品一区二区免费| 精品国产美女av久久久久小说| 麻豆成人午夜福利视频| 久久国产精品影院| 免费观看人在逋| 国产综合懂色| 成人三级黄色视频| 亚洲欧美日韩高清专用| 网址你懂的国产日韩在线| 老汉色av国产亚洲站长工具| 法律面前人人平等表现在哪些方面| 免费在线观看影片大全网站| 亚洲欧美一区二区三区黑人| 国产精品电影一区二区三区| 国产蜜桃级精品一区二区三区| 国产成人av教育| 老熟妇乱子伦视频在线观看| 亚洲av成人av| 桃色一区二区三区在线观看| 999久久久精品免费观看国产| 欧美性猛交╳xxx乱大交人| 欧美国产日韩亚洲一区| 国产午夜精品久久久久久一区二区三区 | 69人妻影院| 在线观看日韩欧美| 亚洲在线自拍视频| 女同久久另类99精品国产91| 国产av在哪里看| 嫩草影院精品99| 久久久久性生活片| 久久精品亚洲精品国产色婷小说| 国产精品影院久久| 美女黄网站色视频| 精品免费久久久久久久清纯| 美女 人体艺术 gogo| 真实男女啪啪啪动态图| 网址你懂的国产日韩在线| 亚洲片人在线观看| 母亲3免费完整高清在线观看| 校园春色视频在线观看| 亚洲一区二区三区色噜噜| 观看美女的网站| 精品久久久久久久久久免费视频| 久久精品夜夜夜夜夜久久蜜豆| 成人性生交大片免费视频hd| 欧美不卡视频在线免费观看| 黄色视频,在线免费观看| 国产 一区 欧美 日韩| 欧美又色又爽又黄视频| 亚洲欧美日韩高清专用| 久久香蕉精品热| 一个人免费在线观看电影| 99热这里只有是精品50| 精品国产美女av久久久久小说| aaaaa片日本免费| 九九在线视频观看精品| 亚洲国产精品成人综合色| 99精品在免费线老司机午夜| 午夜免费观看网址| 法律面前人人平等表现在哪些方面| 亚洲av成人不卡在线观看播放网| 男人和女人高潮做爰伦理| 欧美绝顶高潮抽搐喷水| 淫秽高清视频在线观看| 看免费av毛片| 国产精品,欧美在线| 首页视频小说图片口味搜索| 成人国产一区最新在线观看| 欧美色视频一区免费| 网址你懂的国产日韩在线| 天天添夜夜摸| 小蜜桃在线观看免费完整版高清| 99在线人妻在线中文字幕| 日本三级黄在线观看| 91字幕亚洲| 午夜福利视频1000在线观看| 黄色日韩在线| 两个人视频免费观看高清| 国产伦精品一区二区三区四那| 草草在线视频免费看| 免费在线观看亚洲国产| 久久久久久九九精品二区国产| 日韩av在线大香蕉| 亚洲无线观看免费| 床上黄色一级片| 亚洲精品乱码久久久v下载方式 | xxx96com| 免费av毛片视频| e午夜精品久久久久久久| 精品国产亚洲在线| 性色avwww在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | av天堂在线播放| 国产精华一区二区三区| 欧美最新免费一区二区三区 | 成人三级黄色视频| 深爱激情五月婷婷| 成年免费大片在线观看| 国产精品99久久久久久久久| 国产高清视频在线播放一区| 色吧在线观看| 久久性视频一级片| 看免费av毛片| 欧美一区二区精品小视频在线| 十八禁网站免费在线| 国产成人av教育| 成人高潮视频无遮挡免费网站| 日韩欧美一区二区三区在线观看| 国产精品嫩草影院av在线观看 | 嫩草影院精品99| 又爽又黄无遮挡网站| 欧美绝顶高潮抽搐喷水| 久久这里只有精品中国| 欧美成狂野欧美在线观看| 黄色丝袜av网址大全| 又紧又爽又黄一区二区| 美女黄网站色视频| 久久久国产成人免费| 看黄色毛片网站| 午夜福利18| 日韩欧美精品免费久久 | 窝窝影院91人妻| 免费av毛片视频| 午夜福利成人在线免费观看| 国产一区二区在线观看日韩 | 好男人电影高清在线观看| 亚洲五月婷婷丁香| 久久精品91蜜桃| 一边摸一边抽搐一进一小说| 欧美日韩综合久久久久久 | 国产欧美日韩一区二区三| 中出人妻视频一区二区| 在线视频色国产色| 亚洲第一欧美日韩一区二区三区| 美女免费视频网站| 蜜桃亚洲精品一区二区三区| 久久久久亚洲av毛片大全| 国产淫片久久久久久久久 | 国产午夜精品久久久久久一区二区三区 | 久久久久久久久中文| 五月玫瑰六月丁香| 两个人看的免费小视频| 亚洲av美国av| 精品久久久久久久久久久久久| 午夜福利欧美成人| 国产精品久久电影中文字幕| 国产成人系列免费观看| 老司机福利观看| 97人妻精品一区二区三区麻豆| 在线看三级毛片| 欧美一区二区精品小视频在线| 男人和女人高潮做爰伦理| 亚洲中文字幕日韩| 欧美成狂野欧美在线观看| 熟女电影av网| 一本久久中文字幕| 不卡一级毛片| 成人性生交大片免费视频hd| 伊人久久大香线蕉亚洲五| 在线观看免费午夜福利视频| 久久精品人妻少妇| 亚洲人成伊人成综合网2020| 欧美日本亚洲视频在线播放| 成年女人毛片免费观看观看9| 精品人妻1区二区| 国产中年淑女户外野战色| av在线蜜桃| 青草久久国产| 国产高潮美女av| 亚洲国产色片| 亚洲精品亚洲一区二区| 97超级碰碰碰精品色视频在线观看| 波多野结衣高清无吗| 少妇人妻精品综合一区二区 | av欧美777| 99热这里只有精品一区| 国产不卡一卡二| 国产视频内射| 日韩欧美精品v在线| e午夜精品久久久久久久| 日本五十路高清| 日韩有码中文字幕| 亚洲欧美日韩东京热| 一级a爱片免费观看的视频| 国产精品久久视频播放| 少妇的逼水好多| 国产成年人精品一区二区| 三级毛片av免费| 老司机深夜福利视频在线观看| 88av欧美| 国产日本99.免费观看| 黄色日韩在线| 国内毛片毛片毛片毛片毛片| 尤物成人国产欧美一区二区三区| 日本五十路高清| 久久性视频一级片| 99国产精品一区二区三区| 麻豆国产97在线/欧美| 日本在线视频免费播放| 免费看a级黄色片| 天天躁日日操中文字幕| 国产欧美日韩精品一区二区| 五月伊人婷婷丁香| 亚洲五月婷婷丁香| 欧美+日韩+精品| 哪里可以看免费的av片| 欧美中文日本在线观看视频| 国产免费一级a男人的天堂| 久久亚洲真实| 久久国产精品影院| 男人和女人高潮做爰伦理| 我的老师免费观看完整版| 最近在线观看免费完整版| 欧美一区二区精品小视频在线| bbb黄色大片| 嫩草影视91久久| aaaaa片日本免费| 女人高潮潮喷娇喘18禁视频| www国产在线视频色| 色噜噜av男人的天堂激情| 日本三级黄在线观看| 亚洲激情在线av| 两个人看的免费小视频| 国产精品日韩av在线免费观看| 一区二区三区国产精品乱码| 中文字幕熟女人妻在线| 欧美绝顶高潮抽搐喷水| 中亚洲国语对白在线视频| 久久香蕉精品热| 国产亚洲av嫩草精品影院| 日韩欧美 国产精品| 欧美日韩福利视频一区二区| 国产三级中文精品| 99久国产av精品| 国产高清激情床上av| www.熟女人妻精品国产| 美女高潮喷水抽搐中文字幕| 国产探花在线观看一区二区| 黄色视频,在线免费观看| 伊人久久精品亚洲午夜| 欧美在线一区亚洲| 欧美三级亚洲精品| 无人区码免费观看不卡| 色在线成人网| www.熟女人妻精品国产| 亚洲熟妇熟女久久| 亚洲午夜理论影院| 亚洲欧美一区二区三区黑人| 国产精品1区2区在线观看.| 亚洲五月婷婷丁香| 97碰自拍视频| 蜜桃亚洲精品一区二区三区| 一进一出好大好爽视频| 18美女黄网站色大片免费观看| 欧美日韩综合久久久久久 | 中文字幕人妻丝袜一区二区| 成人精品一区二区免费| 国产伦人伦偷精品视频| 国产精品1区2区在线观看.| 久久九九热精品免费| 亚洲成人久久性| 亚洲成人久久爱视频| 国产野战对白在线观看| 蜜桃久久精品国产亚洲av| 精品乱码久久久久久99久播| 国产伦精品一区二区三区四那| 中文字幕熟女人妻在线| 婷婷亚洲欧美| 成人无遮挡网站| 日韩中文字幕欧美一区二区| 给我免费播放毛片高清在线观看| 在线观看av片永久免费下载| 免费人成视频x8x8入口观看| 精品一区二区三区av网在线观看| 国产精品av视频在线免费观看| 在线看三级毛片| 91久久精品电影网| 嫁个100分男人电影在线观看| 亚洲国产精品合色在线| 最近视频中文字幕2019在线8| 99精品在免费线老司机午夜| 色综合亚洲欧美另类图片| 少妇的丰满在线观看| 婷婷精品国产亚洲av在线| 精品欧美国产一区二区三| 又紧又爽又黄一区二区| 欧美日韩乱码在线| 免费在线观看日本一区| 草草在线视频免费看| 亚洲男人的天堂狠狠| 久久精品人妻少妇| 久久久国产成人免费| 日本成人三级电影网站| 香蕉丝袜av| 午夜日韩欧美国产| 久久久久久国产a免费观看| 麻豆一二三区av精品| 我要搜黄色片| 婷婷精品国产亚洲av| 午夜精品在线福利| 女人高潮潮喷娇喘18禁视频| 在线观看日韩欧美| 欧美三级亚洲精品| 两人在一起打扑克的视频| 午夜a级毛片| 最近最新中文字幕大全电影3| 久久久国产精品麻豆| 国产毛片a区久久久久| 美女cb高潮喷水在线观看| 波多野结衣高清无吗| 内地一区二区视频在线| АⅤ资源中文在线天堂| 熟妇人妻久久中文字幕3abv| 日韩高清综合在线| 婷婷精品国产亚洲av在线| 成人亚洲精品av一区二区| 99热只有精品国产| 亚洲一区二区三区不卡视频| 一区二区三区免费毛片| 国内少妇人妻偷人精品xxx网站| 亚洲精品456在线播放app | 精品人妻1区二区| 亚洲avbb在线观看| 天堂√8在线中文| 国产成+人综合+亚洲专区| 最近最新免费中文字幕在线| 国产免费一级a男人的天堂| 少妇的丰满在线观看| 97超视频在线观看视频| 欧美日本视频| 俺也久久电影网| av天堂在线播放| 首页视频小说图片口味搜索| 宅男免费午夜| 日日摸夜夜添夜夜添小说| 国内毛片毛片毛片毛片毛片| 国产视频一区二区在线看| 久久午夜亚洲精品久久| 成人特级黄色片久久久久久久| 国产淫片久久久久久久久 | 他把我摸到了高潮在线观看| 精品国内亚洲2022精品成人| 欧美性感艳星| 亚洲 欧美 日韩 在线 免费| 亚洲一区高清亚洲精品| 国产精品免费一区二区三区在线| 色综合婷婷激情| 欧美黑人巨大hd| 国产乱人视频| 欧美中文综合在线视频| 久久久久久大精品| 免费av不卡在线播放| 熟女人妻精品中文字幕| 久久这里只有精品中国| 在线国产一区二区在线| 一区二区三区免费毛片| 女人高潮潮喷娇喘18禁视频| 国产色爽女视频免费观看| 欧美一区二区亚洲| 夜夜看夜夜爽夜夜摸| 性欧美人与动物交配| 亚洲18禁久久av| 在线十欧美十亚洲十日本专区| 亚洲av免费高清在线观看| 午夜福利在线在线| www国产在线视频色| 草草在线视频免费看| 51国产日韩欧美| 国内精品久久久久精免费| 久久中文看片网| 精品人妻一区二区三区麻豆 | 色老头精品视频在线观看| 99国产综合亚洲精品| 亚洲一区二区三区不卡视频| 听说在线观看完整版免费高清| 脱女人内裤的视频| 亚洲av电影在线进入| 午夜日韩欧美国产| 亚洲性夜色夜夜综合| 亚洲最大成人中文| 天天躁日日操中文字幕| 三级男女做爰猛烈吃奶摸视频| 亚洲精华国产精华精| 在线天堂最新版资源| 亚洲成av人片免费观看| 麻豆久久精品国产亚洲av| 欧美日韩福利视频一区二区| 国内少妇人妻偷人精品xxx网站| 国产成人aa在线观看| 欧美又色又爽又黄视频| 人妻久久中文字幕网| 亚洲精品成人久久久久久| 国产单亲对白刺激| 一区福利在线观看| 久久伊人香网站| 国产欧美日韩一区二区三| 国产高清三级在线| 亚洲av成人av| 亚洲中文日韩欧美视频| 香蕉丝袜av| 欧美日韩精品网址| 在线观看午夜福利视频| 1000部很黄的大片| 国产成人aa在线观看| 亚洲真实伦在线观看| 校园春色视频在线观看| 欧美黄色片欧美黄色片| 欧美最黄视频在线播放免费| 国产一区二区激情短视频| 亚洲国产欧美网| 91在线精品国自产拍蜜月 | 99国产精品一区二区蜜桃av| 波野结衣二区三区在线 | 国产高清三级在线| 国产一区二区激情短视频| 亚洲一区二区三区不卡视频| 一个人免费在线观看电影| 国产毛片a区久久久久| 桃红色精品国产亚洲av| 91字幕亚洲| 国产精品久久久久久久久免 | 欧美+亚洲+日韩+国产| 蜜桃亚洲精品一区二区三区| 国产91精品成人一区二区三区| 国产激情欧美一区二区| 久久精品91无色码中文字幕| 亚洲成人精品中文字幕电影| 99精品在免费线老司机午夜| 久久久久国内视频| 午夜影院日韩av|