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

    濱海鹽漬化土壤中氨氧化微生物豐度和多樣性特征①

    2022-02-09 09:51:16宋延靜張曉黎
    土壤 2022年6期

    宋延靜,張曉黎,付 嬈,李 萌,王 潔,馬 蘭*

    濱海鹽漬化土壤中氨氧化微生物豐度和多樣性特征①

    宋延靜1,2,3,張曉黎2,付 嬈1,3,李 萌1,3,王 潔1,馬 蘭1,3*

    (1 山東省蠶業(yè)研究所,山東煙臺 264002; 2 中國科學(xué)院海岸帶環(huán)境過程與生態(tài)修復(fù)重點(diǎn)實(shí)驗(yàn)室(煙臺海岸帶研究所),山東煙臺 264003;3 山東省農(nóng)業(yè)科學(xué)院黃河三角洲現(xiàn)代農(nóng)業(yè)研究院,山東東營 257091)

    為探究濱海鹽土不同鹽度梯度下氨氧化微生物的豐度和多樣性特征,利用土壤化學(xué)和分子生態(tài)學(xué)手段(定量PCR、T-RFLP)對萊州灣南岸及黃河口4個河口斷面(黃河、白浪河、堤河、膠萊河)14個不同鹽度(6.4‰ ~ 51.1‰)鹽漬化土壤樣品的氨氧化古菌(AOA)和氨氧化細(xì)菌(AOB)的硝化潛勢、豐度及多樣性進(jìn)行了分析。結(jié)果發(fā)現(xiàn):土壤硝化潛勢在高鹽度(34.7‰ ~51.1‰)條件下被顯著抑制,主要受土壤鹽度、pH和NO3–-N水平顯著影響;AOA-基因豐度比AOB-高出兩個數(shù)量級,在中鹽度時豐度最高(9.92×106copies/g土),在高鹽度時受到顯著抑制(5.28×106copies/g土,<0.05);AOB-基因豐度受鹽度的影響,低鹽度時顯著高于中、高鹽度條件;然而AOA和AOB的多樣性和群落結(jié)構(gòu)受鹽度梯度影響不大。相關(guān)分析表明,硝化潛勢與AOA和AOB豐度均無顯著相關(guān)性,而與AOA/AOB比值以及AOA的Shannon指數(shù)顯著負(fù)相關(guān)。由此可見,濱海鹽土中,鹽度的波動對土壤氨氧化活性和功能微生物都會產(chǎn)生劇烈的影響,鹽度和pH造成的土壤氮素有效性的變化可能是影響濱海鹽土硝化活性和氨氧化微生物豐度及群落組成的關(guān)鍵因素。

    濱海鹽土;硝化潛勢;氨氧化細(xì)菌和古菌;鹽度梯度;群落結(jié)構(gòu)

    土壤鹽漬化是造成濱海土壤質(zhì)量下降的重要原因之一。濱海鹽土高鹽度、高pH的特征對土壤肥力和微生物的多樣性及生態(tài)功能均有不利影響[1]。濱海鹽土氮的有效性低、保水保肥性差,嚴(yán)重制約區(qū)域的土壤地力。解析濱海鹽土氮循環(huán)微生物的活性、多樣性及分布規(guī)律,能夠?yàn)槔斫鉃I海土壤生態(tài)系統(tǒng)功能退化提供基礎(chǔ)數(shù)據(jù),對改善土壤性質(zhì),提高肥力,探索并最終形成鹽漬化生態(tài)系統(tǒng)修復(fù)的技術(shù)體系具有重要意義。

    大量研究表明,鹽分會抑制氮的礦化作用。但也有部分研究顯示,輕度鹽漬化會促進(jìn)氮的礦化,而重度鹽漬化會抑制氮的礦化,尤其是當(dāng)土壤中含鹽量達(dá)到0.5% ~ 1.0% 時,土壤的硝化作用會被完全抑制[2-4]。硝化過程受微生物驅(qū)動,其中氨氧化過程是其限速步驟,由化能自養(yǎng)的氨氧化細(xì)菌(AOB)和氨氧化古菌(AOA)共同完成。這兩種菌都含有編碼催化氨氧化第一步反應(yīng)的氨單加氧酶基因()。利用氨單加氧酶α亞基基因()作為分子標(biāo)記可研究AOB和AOA在環(huán)境中的多樣性、群落組成及豐度特征[5]。鹽度也是影響土壤氨氧化微生物功能和多樣性的重要因子,高鹽和低鹽環(huán)境中,AOA和AOB的群落組成也差別迥異[6]。Sahan等[7]證實(shí)韋斯特謝爾德河口的AOB多樣性在低鹽區(qū)高于高鹽區(qū)。Mosier和Francis[8]在對舊金山灣沉積物AOA的分析中發(fā)現(xiàn),來自低鹽河口區(qū)的AOA形成了獨(dú)特的低鹽群。濱海鹽土除了具有明顯的鹽度梯度外,還有營養(yǎng)貧瘠、pH偏高等特征。然而目前針對這種土壤的AOA和AOB活性、組成、豐度及其環(huán)境調(diào)控機(jī)制尚不清楚。

    萊州灣沿岸海水入侵面積超過4 300 km2,鹽漬化面積達(dá)2 100 km2,西南岸鹽漬化程度最高,土壤自沿海到陸地呈現(xiàn)明顯的梯度鹽漬化特征。前期研究發(fā)現(xiàn),該環(huán)境鹽度對藍(lán)細(xì)菌的組成和分布具有強(qiáng)烈影響[9],但鹽度梯度對氮循環(huán)微生物的影響并不清楚。由此,本研究擬利用地球化學(xué)、生物化學(xué)與微生物分子生態(tài)學(xué)交叉手段,研究萊州灣濱海鹽漬土不同鹽度梯度下AOA和AOB的多樣性、群落組成、活性特征及其對環(huán)境因子的響應(yīng),以為深入研究鹽漬化土壤氮素循環(huán)特征及其功能類群的適應(yīng)機(jī)制提供理論依據(jù)。

    1 材料與方法

    1.1 樣品采集及處理

    本研究在萊州灣西南岸(119°09′ ~ 119°36′ E,36°56′ ~ 37°46′ N)的黃河口(YR)、白浪河口(BR)、堤河口(DR)及膠萊河口(JR)布設(shè)4個樣帶斷面、14個站位(圖1)進(jìn)行樣品采集。其中,黃河口斷面5個站位(YR01 ~ YR05),白浪河口斷面2個站位(DR01和DR02),堤河口斷面3個站位(DR01 ~ DR03),膠萊河口斷面4個站位(JR01 ~ JR04)。

    樣品采集時間為2014年7月,五點(diǎn)取樣法取表層0 ~ 5 cm土樣混勻。樣品冰上保存轉(zhuǎn)移到實(shí)驗(yàn)室,一部分4℃保鮮并于48 h內(nèi)進(jìn)行氨氧化活性測定,另一部分–80℃ 冷凍,用于測定土壤環(huán)境因子和提取土壤DNA。

    1.2 土壤理化因子的測定

    土壤pH測定采用電位法;土壤含水量測定采用烘干法;土壤鹽度通過測定水土質(zhì)量比為1∶5土壤溶液的電導(dǎo)率來獲得。土壤經(jīng)冷凍干燥后,2.5 g干土中加入25 ml 2 mol/L KCl溶液浸提,振蕩1 h,過濾,利用連續(xù)流動分析儀測定土壤銨態(tài)氮(NH4+-N)、硝態(tài)氮(NO3–-N)和亞硝態(tài)氮(NO2–-N)含量。土壤全氮(TN)、總有機(jī)碳(TOC)含量利用元素分析儀測定;土壤粒徑(GS)采用激光粒度儀分析。

    1.3 土壤硝化潛勢(nitrification potential,NP)的測定

    采用氯酸鹽抑制法[10],以每小時每克干土產(chǎn)生的NO2–-N的量來表征。

    1.4 土壤總DNA的提取與測定

    稱取0.5 ~ 1.0 g新鮮土樣,準(zhǔn)確記錄數(shù)值,通過土壤含水量換算為干土質(zhì)量,以便后續(xù)基因拷貝數(shù)的計(jì)算。利用FastDNA?SPIN?(MP Biomedical,美國)土壤DNA提取試劑盒按說明書要求提取DNA。提取的DNA用微量紫外分光光度計(jì)(NanoDrop 2000c,美國)測定濃度及純度,分裝后于–80℃保存。

    1)AOB和AOA豐度分析。采用熒光定量PCR技術(shù)測定AOB-(amoA1F/amoA2R)和AOA-(Arch-amoA-for/Arch-amoA-rev)基因的拷貝數(shù)。計(jì)算質(zhì)??截悢?shù)時,分別稀釋AOA和AOB質(zhì)粒至10–1~ 10–8,以拷貝數(shù)的對數(shù)為橫坐標(biāo),CT值為縱坐標(biāo)建立標(biāo)準(zhǔn)曲線。每個樣品3次重復(fù),并設(shè)不加模板的反應(yīng)管為陰性對照。所有PCR反應(yīng)的擴(kuò)增效率均在90% 以上,相關(guān)系數(shù)(2)均大于98%。

    圖1 萊州灣西南岸河口采樣站位圖

    2)AOB和AOA群落結(jié)構(gòu)分析。以土壤總DNA為模板,采用末端限制性片斷長度多態(tài)性方法(T-RFLP)分析鹽度梯度土壤中AOB(amoA1F/ amoA2R)和AOA(Arch-amoAF/Arch-amoAR)的群落結(jié)構(gòu)。其中每對引物的前引物5'端標(biāo)記有熒光基團(tuán)6-羧基熒光素(FAM)。AOA-內(nèi)切酶為RsaI和HhaI,AOB-內(nèi)切酶為RsaI和MspI,酶切產(chǎn)物送生工生物工程(上海)股份有限公司進(jìn)行毛細(xì)管電泳檢測。

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

    利用SPSS統(tǒng)計(jì)軟件中的單因素方差分析(ANOVA)檢驗(yàn)鹽度梯度對硝化潛勢、AOB和AOA的豐度及多樣性的影響,并利用最小顯著差異法(LSD)比較不同鹽度間的差異;利用Spearman相關(guān)性分析檢驗(yàn)硝化潛勢、AOB和AOA的豐度及多樣性與環(huán)境因子間的關(guān)系。通過群落生態(tài)分析軟件PRIMER-6對AOA和AOB的群落結(jié)構(gòu)進(jìn)行非度量多維尺度(NMDS)排序,并利用ANOSIM檢驗(yàn)鹽度梯度對AOB和AOA群落結(jié)構(gòu)的影響。利用CANOCO軟件中的典范對應(yīng)分析(CCA或RDA)探索氨氧化微生物群落結(jié)構(gòu)與環(huán)境變量之間的關(guān)系。

    2 結(jié)果與分析

    2.1 土壤理化因子

    本研究獲得的所有樣品可劃分為低鹽度(鹽度平均值6.8‰,變幅6.4‰ ~ 7.6‰)、中鹽度(鹽度平均值16.44‰,變幅12.5‰ ~ 23.9‰)、高鹽度(鹽度平均值44.47‰,變幅34.7‰ ~ 51.3‰)3個梯度[11], 其中低鹽度包括站位YR05、JR01、DR01,中鹽度包括站位YR01、YR04、BR01、BR02、DR03、JR02、JR03、JR04,高鹽度包括YR02、DR02、YR03。pH范圍為7.35 ~ 7.58,屬于弱堿性土壤。無機(jī)氮營養(yǎng)鹽中以NO3–-N為主,其中,低、中、高鹽度土壤NO3–-N含量平均值分別為6.05、26.37和23.25 mg/kg。NH4+-N和NO2–-N含量不同鹽度梯度間差別不大。全氮含量范圍為0.59 ~ 0.88 g/kg,總有機(jī)碳含量5.77 ~ 7.19 g/kg(表1)。

    2.2 土壤硝化潛勢

    從圖2中可以看出,高鹽度土壤硝化潛勢平均為0.12 μg/(g·h),與低鹽度(0.28 μg/(g·h))和中鹽度(0.33 μg/(g·h))相比顯著降低(<0.01),說明高鹽度顯著抑制了土壤的硝化潛勢。

    表1 土壤樣品理化因子(n=3)

    注:d0.5表示平均粒徑,指粒度累計(jì)分布中50% 所對應(yīng)的直徑。

    (柱圖上方不同小寫字母表示組間差異顯著(P<0.05),下同)

    2.3 基于amoA基因熒光定量PCR的AOB和AOA豐度

    萊州灣濱海鹽土中AOB的豐度為104~ 105copies/g土,AOA的豐度為106copies/g土(圖3)。所有樣品中AOA-基因豐度高出AOB-約2個數(shù)量級,說明萊州灣濱海鹽土中的氨氧化微生物群落以AOA為主。

    AOA-基因豐度在中鹽度土壤中的平均值最高,為9.92×106copies/g土,其次是低鹽度土壤(8.45×106copies/g土),高鹽度土壤中最低(5.28×106copies/g土),其中高鹽度土壤中AOA-的豐度顯著低于低、中鹽度土壤(<0.05,圖3A)。

    圖3 不同鹽度梯度間AOA(A)和AOB(B)amoA基因豐度的比較

    AOB-基因豐度在鹽度梯度間存在極顯著性差異(<0.01,圖3B),其中低鹽度土壤AOB-基因豐度為3.48×105copies/g土,顯著高于中鹽度土壤(4.90×104copies/g土)和高鹽度土壤(4.16×104copies/g土)。

    2.4 基于T-RFLP分析的土壤AOA和AOB多樣性和群落結(jié)構(gòu)

    除樣品JR05的AOA-基因沒有獲得理想的結(jié)果外,其他樣品經(jīng)T-RFLP分析發(fā)現(xiàn),AOA-和AOB-的T-RF數(shù)目在低、中、高鹽度土壤中均無顯著性差異(表2)。AOA的均勻度指數(shù)和Shannon指數(shù)在低、高鹽度土壤中差異不顯著(>0.05),但均顯著高于中鹽度土壤(<0.05)(表2)。

    NMDS分析結(jié)果顯示,所有樣品的AOA類群可分為2個群組:A組(BR01、BR02、DR03、YR04、JR01、YR02和YR03)和B組(JR04、JR02、JR03、YR01和DR02)。低鹽度樣品DR01的AOA群落結(jié)構(gòu)與其他樣品差異較大(圖4A)。ANOSIM檢驗(yàn)(表3)發(fā)現(xiàn),AOA的群落結(jié)構(gòu)不受鹽度梯度影響(0.23)。AOB的NMDS分析結(jié)果(圖4B)顯示,所有樣品被分為3個群組:A組包括DR01、DR02、DR03、JR01和BR01;B組包括YR01、YR02、YR03、YR05、JR03、和BR02;C組包括YR04和JR04。

    表2 基于T-RFLP分析的濱海鹽土AOA和AOB的多樣性指數(shù)

    注:鹽度梯度間的總體差異水平用值指征;同列不同小寫字母表示組間差異顯著(<0.05)。

    圖4 AOA(A)和AOB(B)amoA基因群落結(jié)構(gòu)的非度量多維標(biāo)度(NMDS)分析

    表3 鹽度梯度間AOA和AOB amoA基因群落結(jié)構(gòu)的ANOSIM分析

    注:值用來檢驗(yàn)組間(兩組或多組)差異是否顯著大于組內(nèi)差異,<0.05時表示組間差異大于組內(nèi)差異。

    2.5 土壤AOA和AOB豐度、多樣性和群落結(jié)構(gòu)及環(huán)境因子間的關(guān)系

    從表4可以看出,萊州灣濱海鹽土的硝化潛勢與鹽度(= –0.58)和pH(= –0.55)顯著負(fù)相關(guān)(<0.05),而與NO3–-N水平顯著正相關(guān)(=0.37,<0.05)。AOA的豐度與土壤pH顯著正相關(guān)(=0.39,<0.05),而與土壤全氮(= –0.36)、總有機(jī)碳(= –0.37)、NH4+-N(= –0.45)和NO2–-N含量(= –0.33)顯著負(fù)相關(guān)(<0.05)。AOB的豐度與所有環(huán)境因子的相關(guān)性均不顯著(>0.05)。另外,AOA/AOB比值與土壤pH顯著正相關(guān)(=0.36,<0.05),而與全氮含量(= –0.39,<0.05)顯著負(fù)相關(guān)。AOA和AOB的Shannon指數(shù)與所有環(huán)境因子均無顯著相關(guān)性(>0.05)。典范對應(yīng)分析發(fā)現(xiàn),AOA和AOB的群落結(jié)構(gòu)與所有環(huán)境因子均無顯著性相關(guān)(>0.05)。

    將硝化潛勢與微生物指標(biāo)(豐度、AOA/AOB比值、Shannon指數(shù))進(jìn)行Spearman相關(guān)分析發(fā)現(xiàn),硝化潛勢與AOA和AOB的基因豐度并無顯著相關(guān)(>0.05),但與AOA/AOB比值存在顯著負(fù)相關(guān)關(guān)系(= –0.35,=0.04);另外,硝化潛勢與AOA的Shannon指數(shù)也呈顯著負(fù)相關(guān)關(guān)系(= –0.68,=0.008)。

    表4 AOA和AOB豐度、多樣性指數(shù)與環(huán)境因子間的Spearman相關(guān)性(ρ值)

    注:為Spearman相關(guān)系數(shù);AOA/AOB為AOA-與AOB-豐度的比值;*、**分別表示相關(guān)性達(dá)0.05和0.01顯著水平。

    3 討論

    3.1 萊州灣西南岸濱海鹽土特性

    萊州灣河口區(qū)土壤具有明顯的鹽度梯度特征,鹽度范圍為6.4‰ ~ 51.3‰,跨度較大,總體沿入海方向呈升高趨勢。土壤屬弱堿性,肥力較低,無機(jī)氮(NH4+-N、NO3–-N和NO2–-N總和)含量在4.33 ~ 132.98 mg/kg,遠(yuǎn)遠(yuǎn)低于耕地的氮水平,有機(jī)碳含量也低于一般土壤水平[12]。

    3.2 濱海鹽土硝化潛力的特征及影響因素

    萊州灣濱海鹽土的潛在硝化速率在0.11 ~ 0.51 μg/(g·h),處于大多數(shù)河口區(qū)硝化速率的低值區(qū)[13-14],可能的原因是這些河口區(qū)的研究對象為沉積物,其中氮含量相對較高,而本研究區(qū)域土壤NH4+-N含量非常低,平均僅為0.14 mg/kg,因此限制了硝化活性。將所有樣品按鹽度高(34.7‰ ~ 51.3‰)、中(12.5‰ ~ 23.9‰)、低(6.4‰ ~ 7.6‰)分組發(fā)現(xiàn),土壤硝化潛勢并不是簡單地隨鹽度升高而降低,而是中鹽度最高,但與低鹽度差異不顯著,高鹽度最低(圖2)。該結(jié)果與以往的報道一致,比如Colne河口[14]、Randers Fjord河口[15]和Douro河口[16]。鹽度對硝化潛勢的影響機(jī)制還不明確。有研究認(rèn)為,鹽度升高會降低土壤對NH4+-N的吸附能力[17],使氨氧化反應(yīng)的底物減少,因此限制了氨氧化活性。在低鹽度情況下,土壤NH4+-N可利用性較高,高濃度NH4+-N可能會抑制很多氨氧化微生物類群的活性[13],特別對AOA具有抑制作用[18];中鹽度土壤的NH4+-N水平可能更適合氨氧化活性;而高鹽度土壤中NH4+-N可利用性低,又限制了氨氧化活性[12]。同時有研究認(rèn)為,鹽度變化會影響氨氧化微生物的生理,改變氨氧化微生物的群落結(jié)構(gòu),進(jìn)而影響土壤氨氧化活性[19]。

    Spearman相關(guān)性分析顯示(表5),除鹽度外,硝化潛勢還與pH顯著負(fù)相關(guān)。萊州灣西南岸土壤pH范圍為7.11 ~ 7.94,屬于弱堿性土壤。大量研究表明,土壤pH是影響硝化活性的關(guān)鍵因素,pH 7.5 ~ 8.5是土壤硝化活性的最佳pH范圍,而酸性和堿性條件都會抑制土壤的硝化活性[20-21]。一方面可能是堿性和酸性條件下產(chǎn)生的自由NH3和亞硝酸對氨氧化微生物細(xì)胞的毒性作用限制了硝化活性[22];另一方面可能是酸性條件下,NH3被離子化變成NH4+,增加了氨氧化的能量需求[23]。按照這些假設(shè),本研究中的pH范圍應(yīng)該比較適合土壤的硝化活性,然而結(jié)果卻是硝化活性隨土壤pH的升高而降低,推測其可能原因是,該區(qū)域土壤堿度升高是海水入侵導(dǎo)致的土壤鹽漬化所致,因此,這種負(fù)相關(guān)可能是鹽度和pH共同作用的結(jié)果。

    3.3 濱海鹽土氨氧化微生物的豐度特征和環(huán)境響應(yīng)

    萊州灣濱海鹽土所有樣品中氨氧化微生物以AOA為主,這與國內(nèi)外大多數(shù)河口環(huán)境的研究結(jié)果相一致[13,24]。當(dāng)然,也有例外,在Douro河口[16]和Cochin河口[25],AOB的豐度顯著高于AOA,較高的NH4+-N水平被認(rèn)為可能是造成這種特殊現(xiàn)象的主要原因。AOA-基因豐度在中鹽度土壤中較高,高鹽度土壤中最低(圖3A),而AOB豐度對鹽度的響應(yīng)趨勢與AOA不同,在中、高鹽度土壤中均被顯著抑制(圖3B)。AOA和AOB對鹽度梯度的響應(yīng)與很多河口類似[26],即相對于AOB,AOA對鹽度的耐受能力更高。這可能是由于AOA和AOB對NH4+-N親和力的不同所導(dǎo)致的[27],即AOB偏愛高NH4+-N環(huán)境,低鹽促進(jìn)土壤對NH4+-N的吸附,因而有利于AOB的生長;而AOA對NH4+-N的親和力較高,偏愛低NH4+-N環(huán)境,中鹽度土壤的NH4+-N水平可能比較適合AOA生長。當(dāng)然,也有可能是AOA和AOB細(xì)胞對鹽度的生理反應(yīng)不同從而導(dǎo)致了這種生境選擇[19]。AOA豐度與土壤pH顯著正相關(guān),與已有的認(rèn)知“AOA在酸性環(huán)境中占主導(dǎo),而AOB在中性和堿性土壤中占優(yōu)勢”相悖[20,23]。在萊州灣西南岸鹽土中,土壤pH在氨氧化微生物的適宜范圍內(nèi),此時pH的影響比較弱,而有限的氮底物則可能是氨氧化微生物生長的限制因素,pH越高,越利于NH4+-N向NH3的轉(zhuǎn)化,因而促進(jìn)了AOA的生長。此外,AOA豐度與土壤肥力顯著負(fù)相關(guān)(表5)。一般情況下,AOA在肥力低的環(huán)境中貢獻(xiàn)較大,主要是因?yàn)锳OA對底物氨的親和力遠(yuǎn)高于AOB[27],因此在低肥力土壤中更具競爭優(yōu)勢。白浪河口有效氮含量最低,因此,該區(qū)域AOA豐度明顯高于其他河口,而AOB豐度則相反。

    3.4 濱海鹽土氨氧化微生物群落的分布特征

    鹽度被認(rèn)為是控制河口區(qū)AOB和AOA群落結(jié)構(gòu)和多樣性最主要的因素。Bernhard等[13]發(fā)現(xiàn)在低、中鹽度樣品中,AOB主要隸屬于亞硝化螺菌屬相關(guān)類群以及亞硝化單胞菌屬,而在高鹽度環(huán)境中,AOB群落結(jié)構(gòu)演替為以亞硝化螺菌屬為主。另外,在韋斯特謝爾德河口,低鹽度站點(diǎn)AOB的多樣性較高鹽度站點(diǎn)高[28]。在對舊金山灣沉積物AOA的分析中發(fā)現(xiàn),來自低鹽度河口區(qū)的AOA形成了獨(dú)特的低鹽度類群[6]。然而本研究中,鹽度對AOB和AOA的多樣性和群落結(jié)構(gòu)均無顯著影響。

    3. 5 硝化潛勢與氨氧化微生物的關(guān)系

    硝化潛勢與AOA或AOB的豐度均無顯著相關(guān)性,而與AOA/AOB比值顯著負(fù)相關(guān)(表5)。這說明在萊州灣西南岸濱海鹽土中,雖然AOA是優(yōu)勢氨氧化菌,但其豐度變化并不能很好地解釋氨氧化活性的變化。該結(jié)果與之前本課題組在該區(qū)域堿土中(鹽度0.52‰,pH 8.0)添加生物質(zhì)炭后,硝化潛勢與AOA豐度顯著相關(guān)的結(jié)果不一致[29],暗示了環(huán)境因素(比如鹽度及其造成的NH4+-N有效性差異,以及生物質(zhì)炭添加干預(yù))是調(diào)控該區(qū)域?yàn)I海鹽土硝化潛勢的主要因素。而且,即使該環(huán)境中AOA在豐度上占優(yōu)勢,這些AOA也不一定表現(xiàn)氨氧化活性。另外,AOA/AOB比值一般受土壤NH4+-N水平調(diào)控,隨NH4+-N水平的升高而降低。在本研究貧瘠的土壤中,NH4+-N恰恰是硝化活性的限制因素,因此硝化潛勢與AOA/AOB比值間的負(fù)相關(guān)關(guān)系,有可能是其對氮素需求的間接表現(xiàn)。

    4 結(jié)論

    萊州灣西南岸濱海鹽土是非常脆弱的生態(tài)系統(tǒng),鹽度高、pH高、營養(yǎng)匱乏。在這種環(huán)境中,鹽度的波動對土壤氨氧化活性和功能微生物都會產(chǎn)生劇烈的影響。總體來說,中鹽度(12.5‰ ~ 23.9‰)比較適宜土壤硝化活性,而高鹽度會顯著抑制硝化活性。AOA和AOB對鹽度波動響應(yīng)并不一致:AOA偏愛中鹽度環(huán)境,而AOB則更喜低鹽度環(huán)境。鹽度和pH造成的土壤氮素有效性的變化可能是影響濱海鹽土硝化活性和微生物豐度及群落組成的核心因素。

    [1] Dendooven L, Alcántara-Hernández R J, Valenzuela- Encinas C, et al. Dynamics of carbon and nitrogen in an extreme alkaline saline soil: A review[J]. Soil Biology and Biochemistry, 2010, 42(6): 865–877.

    [2] 張先富, 李卉, 洪梅, 等. 蘇打鹽堿土對氮轉(zhuǎn)化的影響[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版), 2012, 42(4): 1145–1150.

    [3] Reddy N, Crohn D M. Effects of soil salinity and carbon availability from organic amendments on nitrous oxide emissions[J]. Geoderma, 2014, 235/236: 363–371.

    [4] 李紅強(qiáng), 姚榮江, 楊勁松, 等. 鹽漬化對農(nóng)田氮素轉(zhuǎn)化過程的影響機(jī)制和增效調(diào)控途徑[J]. 應(yīng)用生態(tài)學(xué)報, 2020, 31(11): 3915–3924.

    [5] Gubry-Rangin C, Nicol G W, Prosser J I. Archaea rather than bacteria control nitrification in two agricultural acidic soils[J]. FEMS Microbiology Ecology, 2010, 74(3): 566–574.

    [6] Bernhard A E, Donn T, Giblin A E, et al. Loss of diversity of ammonia-oxidizing bacteria correlates with increasing salinity in an estuary system[J]. Environmental Microbiology, 2005, 7(9): 1289–1297.

    [7] Sahan E, Muyzer G. Diversity and spatio-temporal distribution of ammonia oxidizing archaea and bacteria in sediments of the Westerschelde estuary. FEMS Microbiology Ecology, 2008, 64: 175–186.

    [8] Mosier A C, Francis C A. Relative abundance and diversity of ammonia-oxidizing Archaea and bacteria in the San Francisco Bay Estuary[J]. Environmental Microbiology, 2008, 10(11): 3002–3016.

    [9] 李寒, 張曉黎, 郭曉紅, 等. 濱海鹽漬化土壤中藍(lán)細(xì)菌多樣性及分布[J]. 微生物學(xué)通報, 2015, 42(5): 957–967.

    [10] Kurola J, Salkinoja-Salonen M, Aarnio T, et al. Activity, diversity and population size of ammonia-oxidizing bacteria in oil-contaminated landfarming soil[J]. FEMS Microbiology Letters, 2005, 250(1): 33–38.

    [11] В.А.科夫達(dá) , В.С.穆娜朵娃 , В.В.葉戈洛夫 , В.П.施特爾戈洛夫 , 祝壽泉. 按鹽漬化程度和性質(zhì)并結(jié)合植物耐鹽性而擬定的土壤分類[J]. 土壤, 1962(3): 59–62.

    [12] 李二煥, 沈俊, 鞠靖, 等. 蘇北濱海鹽土區(qū)河流入??谕寥利}分及養(yǎng)分分布特征[J]. 中國水土保持科學(xué), 2016, 14(1): 79–88.

    [13] Bernhard A E, Tucker J, Giblin A E, et al. Functionally distinct communities of ammonia-oxidizing bacteria along an estuarine salinity gradient[J]. Environmental Microbiology, 2007, 9(6): 1439–1447.

    [14] Li J L, Nedwell D B, Beddow J, et al.Gene abundances and nitrification potential rates suggest that benthic ammonia-oxidizing bacteria and not Archaea dominate N cycling in the Colne Estuary, United Kingdom[J]. Applied and Environmental Microbiology, 2015, 81(1): 159–165.

    [15] Meyer R L, Kj?r T, Revsbech N P. Use of NOx- microsensors to estimate the activity of sediment nitrification and NOx- consumption along an estuarine salinity, nitrate, and light gradient[J]. Aquatic Microbial Ecology, 2001, 26: 181–193.

    [16] Magalh?es C M, Joye S B, Moreira R M, et al. Effect of salinity and inorganic nitrogen concentrations on nitrification and denitrification rates in intertidal sediments and rocky biofilms of the Douro River Estuary, Portugal[J]. Water Research, 2005, 39(9): 1783–1794.

    [17] Rysgaard S, Thastum P, Dalsgaard T, et al. Effects of salinity on NH4+adsorption capacity, nitrification, and denitrification in Danish estuarine sediments[J]. Estuaries, 1999, 22(1): 21.

    [18] Hatzenpichler R, Lebedeva E V, Spieck E, et al. A moderately thermophilic ammonia-oxidizing crenarchaeote from a hot spring[J]. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(6): 2134–2139.

    [19] Bollmann A, Laanbroek H J. Influence of oxygen partial pressure and salinity on the community composition of ammonia-oxidizing bacteria in the Schelde Estuary[J]. Aquatic Microbial Ecology, 2002, 28: 239–247.

    [20] Jiang X J, Hou X Y, Zhou X, et al. pH regulates key players of nitrification in paddy soils[J]. Soil Biology and Biochemistry, 2015, 81: 9–16.

    [21] Sahrawat K L. Factors affecting nitrification in soils[J]. Communications in Soil Science and Plant Analysis, 2008, 39(9/10): 1436–1446.

    [22] Halling-Sorensen B, Jorgensen S E. Studies in environmental science[M]. Elsevier, 1993: 55–118.

    [23] Nicol G W, Leininger S, Schleper C, et al. The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing Archaea and bacteria[J]. Environmental Microbiology, 2008, 10(11): 2966–2978.

    [24] Yao H Y, Campbell C D, Chapman S J, et al. Multi-factorial drivers of ammonia oxidizer communities: Evidence from a national soil survey[J]. Environmental Microbiology, 2013, 15(9): 2545–2556.

    [25] Puthiya Veettil V, Abdulaziz A, Chekidhenkuzhiyil J, et al. Bacterial domination over Archaea in ammonia oxidation in a monsoon-driven tropical estuary[J]. Microbial Ecology, 2015, 69(3): 544–553.

    [26] Bernhard A E, Landry Z C, Blevins A, et al. Abundance of ammonia-oxidizing Archaea and bacteria along an estuarine salinity gradient in relation to potential nitrification rates[J]. Applied and Environmental Microbiology, 2010, 76(4): 1285–1289.

    [27] 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.

    [28] Sahan E, Muyzer G. Diversity and spatio-temporal distribution of ammonia-oxidizing Archaea and Bacteria in sediments of the Westerschelde Estuary[J]. FEMS Microbiology Ecology, 2008, 64(2): 175–186.

    [29] Song Y J, Zhang X L, Ma B, et al. Biochar addition affected the dynamics of ammonia oxidizers and nitrification in microcosms of a coastal alkaline soil[J]. Biology and Fertility of Soils, 2014, 50(2): 321–332.

    Patterns of Activity and Community of Ammonia Oxidizers Along Salinity Gradient in Salinized Soils

    SONG Yanjing1, 2, 3, ZHANG Xiaoli2, FU Rao1,3, LI Meng1,3, WANG Jie1, MA Lan1,3*

    (1 Sericultural Research Institute of Shandong Province, Yantai, Shandong 264002, China; 2 Key Laboratory of Coastal Environmental Processes and Ecological Remediation, Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Yantai, Shandong 264003, China; 3 Institute of Modern Agriculture on Yellow River Delta, Shandong Academy of Agricultural Sciences, Dongying, Shandong 257091, China)

    To explore the effects of the salinity gradient on the abundance and diversity of ammonia-oxidizing microorganisms in coastal salinized soils, salinized soil samples with different salinities (6.4‰–51.3‰) were collected from four estuary transects (Yellow River, Bailang River, Di River, and Jiaolai River) in the southwest coast of Laizhou Bay. The methods of soil chemistry and molecular ecology (qPCR and T-RFLP) were applied to analyze the activities, abundance and community structures of ammonia-oxidizing archaea (AOA) and bacteria (AOB). The results showed that soil nitrification potential (NP) was significantly inhibited at high salinity (34.7‰–51.3‰) and it was significantly correlated with soil salinity, pH, and NO3–-N level. Abundance of AOA-gene was two orders of magnitude higher than that of AOB-gene, which was highest at medium salinity (9.92×106copies/g soil), and was inhibited at high salinity (5.28×106copies/g soil,<0.05). Abundance of AOB-was correlated with soil salinity, which was highest at low salinity and significantly reduced at medium and high salinities. However, the diversities and community structures of both AOA and AOB were little affected by the salinity gradient based on the T-RFLP data. Furthermore, correlation analysis showed that NP was not significantly correlated with abundance of AOA or AOB, but negatively significantly correlated with the ratio of AOA/AOB and the Shannon index of AOA. In conclusion, the ammonia oxidation activity and its functional microorganisms in coastal saline soils were affected intensively by the fluctuation of salinity. The changes of soil nitrogen availability caused by salinity and pH may be the key factors affecting nitrification activity, microbial abundance and community composition in coastal saline soils.

    Coastal salinized soil; Nitrification potential; Ammonia-oxidizing bacteria and archaea; Salinity gradient; Community structure

    X172

    A

    10.13758/j.cnki.tr.2022.06.009

    宋延靜, 張曉黎, 付嬈, 等. 濱海鹽漬化土壤中氨氧化微生物豐度和多樣性特征. 土壤, 2022, 54(6): 1157–1164.

    國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2019YFD1002703)、山東省自然科學(xué)基金項(xiàng)目(ZR2021QD145)和山東省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)科技創(chuàng)新工程項(xiàng)目(CXGC2022A37)資助。

    通訊作者(cysmalan@shandong.cn)

    宋延靜(1985—),山東煙臺人,博士,助理研究員,主要從事土壤微生物生態(tài)、氮循環(huán)及鹽堿灘涂生物修復(fù)技術(shù)研究。E-mail: yjsong1214@163.com

    久久国产精品男人的天堂亚洲 | 人体艺术视频欧美日本| 亚洲av二区三区四区| 我的老师免费观看完整版| 欧美日韩精品成人综合77777| a级毛片免费高清观看在线播放| 成人无遮挡网站| 特大巨黑吊av在线直播| 内射极品少妇av片p| 寂寞人妻少妇视频99o| 亚洲国产av新网站| 中文字幕av电影在线播放| 三级国产精品片| 黄色欧美视频在线观看| 国产女主播在线喷水免费视频网站| 亚洲国产精品专区欧美| 午夜福利视频精品| 老熟女久久久| 亚洲经典国产精华液单| 日本vs欧美在线观看视频 | 黄色日韩在线| 搡女人真爽免费视频火全软件| 我要看日韩黄色一级片| 色视频www国产| av视频免费观看在线观看| 大话2 男鬼变身卡| 精品久久久久久久久av| 日韩 亚洲 欧美在线| 国产免费又黄又爽又色| 国产成人精品无人区| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久av不卡| 亚洲精品乱码久久久v下载方式| 伊人久久精品亚洲午夜| 欧美bdsm另类| 久久午夜福利片| 美女视频免费永久观看网站| 黄色一级大片看看| 欧美日韩在线观看h| 久久久精品94久久精品| 亚洲欧美一区二区三区国产| 色5月婷婷丁香| 男人舔奶头视频| 亚洲精品乱久久久久久| 亚洲国产欧美在线一区| 又粗又硬又长又爽又黄的视频| 国产男人的电影天堂91| 精品一区二区三区视频在线| 99热这里只有是精品50| 中文天堂在线官网| 91久久精品国产一区二区三区| 水蜜桃什么品种好| 一本久久精品| 99视频精品全部免费 在线| 黄色一级大片看看| 久久鲁丝午夜福利片| 成人毛片a级毛片在线播放| 国产男人的电影天堂91| 免费在线观看成人毛片| 亚洲欧美日韩东京热| 中文欧美无线码| 欧美性感艳星| 三上悠亚av全集在线观看 | 国产精品一区二区在线不卡| 97精品久久久久久久久久精品| h视频一区二区三区| 纵有疾风起免费观看全集完整版| 噜噜噜噜噜久久久久久91| 少妇人妻精品综合一区二区| 中文字幕久久专区| 大陆偷拍与自拍| 久久国产精品大桥未久av | 男女国产视频网站| 国产成人免费观看mmmm| 国产亚洲最大av| 亚洲高清免费不卡视频| 最近2019中文字幕mv第一页| 亚洲美女搞黄在线观看| 伦精品一区二区三区| 久久精品国产亚洲av涩爱| 国产精品国产av在线观看| 欧美+日韩+精品| www.av在线官网国产| 亚洲成色77777| 国产黄色免费在线视频| 99久久精品热视频| 国国产精品蜜臀av免费| 国产成人免费无遮挡视频| 嘟嘟电影网在线观看| 中文资源天堂在线| 3wmmmm亚洲av在线观看| 国产高清有码在线观看视频| 日本91视频免费播放| 国产精品久久久久久久电影| 久久久精品94久久精品| 欧美成人精品欧美一级黄| 国产欧美亚洲国产| 精品久久久久久电影网| 国产精品久久久久久精品电影小说| 日本欧美视频一区| 一区二区三区精品91| 噜噜噜噜噜久久久久久91| 丰满饥渴人妻一区二区三| 91精品一卡2卡3卡4卡| 亚洲av国产av综合av卡| 久久韩国三级中文字幕| 欧美97在线视频| 边亲边吃奶的免费视频| 国产精品欧美亚洲77777| 亚洲国产精品成人久久小说| 欧美高清成人免费视频www| 亚洲精品日韩在线中文字幕| 色视频在线一区二区三区| 欧美成人精品欧美一级黄| 久久久国产精品麻豆| 色婷婷av一区二区三区视频| 最近的中文字幕免费完整| 18+在线观看网站| 国产精品偷伦视频观看了| 久久久久人妻精品一区果冻| 亚洲成人一二三区av| 女人久久www免费人成看片| 精品一区二区免费观看| 国产美女午夜福利| 国产精品一区二区三区四区免费观看| xxx大片免费视频| 亚洲av免费高清在线观看| h视频一区二区三区| av在线观看视频网站免费| 91在线精品国自产拍蜜月| av线在线观看网站| 男女啪啪激烈高潮av片| 免费观看性生交大片5| 自线自在国产av| 成年美女黄网站色视频大全免费 | 日韩免费高清中文字幕av| 亚洲第一av免费看| 如日韩欧美国产精品一区二区三区 | 久久免费观看电影| 免费观看无遮挡的男女| 免费人成在线观看视频色| 自拍偷自拍亚洲精品老妇| 最近中文字幕高清免费大全6| 久久精品国产鲁丝片午夜精品| 大话2 男鬼变身卡| 亚洲欧美日韩卡通动漫| 青春草视频在线免费观看| 在线观看免费日韩欧美大片 | 成年av动漫网址| 成人亚洲欧美一区二区av| xxx大片免费视频| 国产精品国产三级国产专区5o| 欧美日本中文国产一区发布| 中文精品一卡2卡3卡4更新| 狂野欧美激情性bbbbbb| 亚洲中文av在线| 亚洲精品乱久久久久久| 极品教师在线视频| 22中文网久久字幕| 99久久综合免费| 黄色欧美视频在线观看| 日本色播在线视频| 日韩免费高清中文字幕av| 另类亚洲欧美激情| 视频区图区小说| 精品一区二区三区视频在线| 伦精品一区二区三区| 国产伦精品一区二区三区视频9| 成年人午夜在线观看视频| 欧美老熟妇乱子伦牲交| 国产精品国产三级国产av玫瑰| 熟女人妻精品中文字幕| 国产精品一区二区在线观看99| 精品午夜福利在线看| 亚洲精品第二区| 久久鲁丝午夜福利片| 国产精品国产三级国产av玫瑰| 最近的中文字幕免费完整| av有码第一页| 日韩强制内射视频| 极品少妇高潮喷水抽搐| 欧美日韩在线观看h| 性色avwww在线观看| 国产精品久久久久久久久免| 国产欧美日韩精品一区二区| 国产av一区二区精品久久| 在线观看人妻少妇| 黄色日韩在线| 日韩av免费高清视频| 在线观看www视频免费| .国产精品久久| 一级毛片电影观看| 涩涩av久久男人的天堂| 国产精品无大码| 18禁在线无遮挡免费观看视频| 男人添女人高潮全过程视频| 26uuu在线亚洲综合色| 亚洲色图综合在线观看| 久久这里有精品视频免费| 最黄视频免费看| 久久国内精品自在自线图片| 一个人看视频在线观看www免费| 久久久久久久久大av| 国产精品99久久99久久久不卡 | 欧美激情极品国产一区二区三区 | 国产精品久久久久成人av| 亚洲欧美日韩另类电影网站| videossex国产| 天堂俺去俺来也www色官网| 91成人精品电影| 国产亚洲一区二区精品| 97在线人人人人妻| 亚洲av在线观看美女高潮| 丝袜喷水一区| 久久久国产欧美日韩av| 国产成人精品一,二区| 99国产精品免费福利视频| 久久人妻熟女aⅴ| 视频区图区小说| 91久久精品国产一区二区成人| 国产无遮挡羞羞视频在线观看| 18+在线观看网站| 汤姆久久久久久久影院中文字幕| 国产成人精品婷婷| 国产精品久久久久久精品电影小说| 久热这里只有精品99| 成人毛片a级毛片在线播放| 国产精品免费大片| 国产黄色视频一区二区在线观看| 久久精品熟女亚洲av麻豆精品| 国内少妇人妻偷人精品xxx网站| 夫妻午夜视频| 亚洲图色成人| 69精品国产乱码久久久| 国产中年淑女户外野战色| 久久久久久久久久人人人人人人| 亚洲,一卡二卡三卡| 成人午夜精彩视频在线观看| 一本色道久久久久久精品综合| 另类精品久久| 在线看a的网站| 国产黄片视频在线免费观看| 在线播放无遮挡| 七月丁香在线播放| 中文字幕亚洲精品专区| 日本与韩国留学比较| 人妻人人澡人人爽人人| 亚洲中文av在线| 免费av中文字幕在线| 十八禁高潮呻吟视频 | 特大巨黑吊av在线直播| 国产精品福利在线免费观看| 乱系列少妇在线播放| 黑丝袜美女国产一区| 成人18禁高潮啪啪吃奶动态图 | 少妇精品久久久久久久| 国产永久视频网站| 免费观看的影片在线观看| 黄色欧美视频在线观看| 亚洲精品,欧美精品| 夜夜骑夜夜射夜夜干| 内射极品少妇av片p| 精品久久国产蜜桃| 中文字幕久久专区| 亚洲欧美成人综合另类久久久| 精品国产国语对白av| 欧美成人午夜免费资源| 免费观看av网站的网址| 2018国产大陆天天弄谢| 久久99热这里只频精品6学生| 欧美亚洲 丝袜 人妻 在线| 有码 亚洲区| 少妇 在线观看| 欧美高清成人免费视频www| 一区二区三区精品91| 国产在视频线精品| 国产一区二区三区综合在线观看 | 免费人妻精品一区二区三区视频| 嫩草影院新地址| 免费黄频网站在线观看国产| 久久97久久精品| 自线自在国产av| av网站免费在线观看视频| 免费av中文字幕在线| 丰满饥渴人妻一区二区三| 国产高清不卡午夜福利| 国产精品人妻久久久影院| 美女cb高潮喷水在线观看| 国产视频内射| 国产在线一区二区三区精| 人人妻人人澡人人爽人人夜夜| 日本色播在线视频| 丰满迷人的少妇在线观看| 国产精品国产三级国产专区5o| 亚洲人成网站在线观看播放| 成人毛片a级毛片在线播放| 欧美精品国产亚洲| 久久国内精品自在自线图片| 老熟女久久久| 日本免费在线观看一区| 欧美另类一区| 另类精品久久| 欧美最新免费一区二区三区| 亚洲av男天堂| av卡一久久| 欧美一级a爱片免费观看看| 两个人的视频大全免费| 色哟哟·www| 天堂俺去俺来也www色官网| 精品酒店卫生间| 久久久久久久久大av| 人妻少妇偷人精品九色| 亚洲在久久综合| 各种免费的搞黄视频| 国产精品蜜桃在线观看| 最黄视频免费看| 亚洲av综合色区一区| 亚洲精品久久久久久婷婷小说| 国产片特级美女逼逼视频| 夜夜看夜夜爽夜夜摸| 91精品国产国语对白视频| 一个人看视频在线观看www免费| 我要看日韩黄色一级片| 人人妻人人看人人澡| 有码 亚洲区| 夜夜爽夜夜爽视频| 桃花免费在线播放| 涩涩av久久男人的天堂| 18禁裸乳无遮挡动漫免费视频| 国产乱人偷精品视频| 看非洲黑人一级黄片| 成人二区视频| 国产一区二区三区av在线| 亚洲精品成人av观看孕妇| 永久免费av网站大全| 一本一本综合久久| 天天操日日干夜夜撸| 老司机影院成人| 九九在线视频观看精品| 大片免费播放器 马上看| 精品酒店卫生间| 91在线精品国自产拍蜜月| 久久久久久久久久成人| 欧美日韩视频高清一区二区三区二| 国产午夜精品一二区理论片| 国产精品三级大全| 18+在线观看网站| 男女国产视频网站| 中文资源天堂在线| 日本av手机在线免费观看| 亚洲性久久影院| 男女国产视频网站| 国产伦理片在线播放av一区| 久久久久久久精品精品| 久久久久久久久久成人| 日日爽夜夜爽网站| av福利片在线观看| 春色校园在线视频观看| 国产成人精品无人区| 精品一区二区三区视频在线| 日韩一本色道免费dvd| 日韩 亚洲 欧美在线| 国产精品一区二区性色av| 一本一本综合久久| 天天躁夜夜躁狠狠久久av| 国产片特级美女逼逼视频| 寂寞人妻少妇视频99o| 老司机亚洲免费影院| 大话2 男鬼变身卡| 国产黄色视频一区二区在线观看| 22中文网久久字幕| 国产日韩欧美视频二区| 街头女战士在线观看网站| 麻豆乱淫一区二区| 街头女战士在线观看网站| 热re99久久精品国产66热6| 91aial.com中文字幕在线观看| 免费看av在线观看网站| 日本av免费视频播放| 精品久久久精品久久久| 黄色配什么色好看| 国产免费视频播放在线视频| 人人妻人人看人人澡| 99热网站在线观看| 久久久a久久爽久久v久久| 亚洲精品色激情综合| 日韩视频在线欧美| 91久久精品国产一区二区三区| 精品一品国产午夜福利视频| 色哟哟·www| 国产午夜精品一二区理论片| 一级毛片久久久久久久久女| 国产在线免费精品| 插阴视频在线观看视频| 亚洲欧美精品自产自拍| 最近的中文字幕免费完整| 一级,二级,三级黄色视频| 少妇熟女欧美另类| 国产一区二区三区av在线| 欧美日韩亚洲高清精品| 美女国产视频在线观看| 少妇精品久久久久久久| 夫妻午夜视频| 伦精品一区二区三区| av国产精品久久久久影院| 国产成人freesex在线| 日本爱情动作片www.在线观看| 亚洲国产日韩一区二区| 日韩欧美精品免费久久| 午夜福利网站1000一区二区三区| 欧美日韩av久久| 日日啪夜夜撸| 久久久久久久大尺度免费视频| 少妇人妻 视频| 丁香六月天网| 啦啦啦在线观看免费高清www| 国产免费又黄又爽又色| 欧美另类一区| 亚洲欧美日韩东京热| 99久久综合免费| 亚洲电影在线观看av| 国产淫片久久久久久久久| 久久久久久久久久久久大奶| a级毛色黄片| 五月伊人婷婷丁香| 亚洲国产精品国产精品| 日韩av免费高清视频| 亚洲精华国产精华液的使用体验| 久久97久久精品| 黑丝袜美女国产一区| 日本欧美国产在线视频| 最新的欧美精品一区二区| 嘟嘟电影网在线观看| 国产精品秋霞免费鲁丝片| 国产高清三级在线| 一本—道久久a久久精品蜜桃钙片| 亚洲av成人精品一区久久| 亚州av有码| 青春草视频在线免费观看| 国产视频首页在线观看| 国产伦理片在线播放av一区| 91在线精品国自产拍蜜月| 老司机影院毛片| 一本一本综合久久| 亚洲一区二区三区欧美精品| 天堂俺去俺来也www色官网| 精品国产乱码久久久久久小说| 日日啪夜夜爽| 精品午夜福利在线看| 日韩一区二区视频免费看| 日韩强制内射视频| 欧美+日韩+精品| 美女内射精品一级片tv| 国产成人精品久久久久久| 51国产日韩欧美| 久久99蜜桃精品久久| 成年av动漫网址| 一区二区三区精品91| 国产一区亚洲一区在线观看| 国产乱人偷精品视频| 国产亚洲欧美精品永久| 亚洲精华国产精华液的使用体验| av国产精品久久久久影院| 精品国产一区二区三区久久久樱花| 在线 av 中文字幕| 一区二区三区精品91| 蜜桃久久精品国产亚洲av| 99久久中文字幕三级久久日本| 黄色欧美视频在线观看| 青春草亚洲视频在线观看| 女人久久www免费人成看片| 最后的刺客免费高清国语| 成人毛片a级毛片在线播放| 在线观看www视频免费| 九九久久精品国产亚洲av麻豆| 日本欧美视频一区| 久久午夜综合久久蜜桃| 亚洲中文av在线| 国产精品熟女久久久久浪| 国产精品国产三级专区第一集| 亚洲国产欧美日韩在线播放 | 久热久热在线精品观看| 亚洲国产av新网站| 久久狼人影院| 偷拍熟女少妇极品色| 国产成人精品福利久久| 国内精品宾馆在线| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久精品电影小说| 大话2 男鬼变身卡| 日本爱情动作片www.在线观看| 国产精品久久久久久精品电影小说| 国产有黄有色有爽视频| 亚洲国产最新在线播放| 国产在线免费精品| av一本久久久久| 国产精品国产av在线观看| 亚洲精品国产色婷婷电影| 国产精品国产av在线观看| 噜噜噜噜噜久久久久久91| 久久久久久久久久成人| 国产男人的电影天堂91| 老女人水多毛片| 亚洲一级一片aⅴ在线观看| 国产精品蜜桃在线观看| 国内揄拍国产精品人妻在线| 亚洲av不卡在线观看| 亚洲天堂av无毛| 妹子高潮喷水视频| 2022亚洲国产成人精品| 亚洲人与动物交配视频| 日韩人妻高清精品专区| 男人和女人高潮做爰伦理| 午夜福利在线观看免费完整高清在| 好男人视频免费观看在线| 人人妻人人爽人人添夜夜欢视频 | 青青草视频在线视频观看| 老女人水多毛片| 日本与韩国留学比较| 亚洲欧美中文字幕日韩二区| 伊人久久精品亚洲午夜| 在线观看国产h片| 国产午夜精品一二区理论片| 国产国拍精品亚洲av在线观看| 国产乱人偷精品视频| av在线观看视频网站免费| 国产成人aa在线观看| 久久久久久久久久人人人人人人| 国产免费视频播放在线视频| 女性生殖器流出的白浆| 99久久综合免费| 国产av精品麻豆| 欧美精品亚洲一区二区| 婷婷色综合大香蕉| 男女无遮挡免费网站观看| 亚洲色图综合在线观看| 永久网站在线| 亚洲国产欧美日韩在线播放 | 内射极品少妇av片p| av不卡在线播放| 久久国产乱子免费精品| 又黄又爽又刺激的免费视频.| 亚洲国产欧美在线一区| 高清毛片免费看| 成人毛片60女人毛片免费| 日日撸夜夜添| 丝瓜视频免费看黄片| www.色视频.com| 欧美精品高潮呻吟av久久| 欧美精品国产亚洲| 午夜视频国产福利| 久久影院123| 国产在线视频一区二区| 美女xxoo啪啪120秒动态图| 中国美白少妇内射xxxbb| 国产综合精华液| 亚洲国产日韩一区二区| 亚洲国产最新在线播放| 亚洲精品,欧美精品| 久久久久久久亚洲中文字幕| 在线观看美女被高潮喷水网站| 欧美日韩综合久久久久久| 丰满饥渴人妻一区二区三| 亚洲自偷自拍三级| 内射极品少妇av片p| 国产老妇伦熟女老妇高清| 一区在线观看完整版| 欧美一级a爱片免费观看看| 久久6这里有精品| 欧美精品一区二区大全| 亚洲精品456在线播放app| 九九久久精品国产亚洲av麻豆| 国产高清三级在线| 狂野欧美激情性xxxx在线观看| 亚洲,欧美,日韩| 国产成人a∨麻豆精品| 精品人妻熟女毛片av久久网站| 91久久精品国产一区二区成人| 全区人妻精品视频| 中文乱码字字幕精品一区二区三区| 亚洲av国产av综合av卡| 色视频www国产| 亚洲精品视频女| 最近的中文字幕免费完整| 汤姆久久久久久久影院中文字幕| 在线观看www视频免费| 69精品国产乱码久久久| 中文字幕亚洲精品专区| 欧美日韩视频精品一区| 久热久热在线精品观看| 日本黄大片高清| 午夜激情福利司机影院| 亚洲激情五月婷婷啪啪| 亚洲精品中文字幕在线视频 | 国产亚洲欧美精品永久| 久久av网站| 精品少妇内射三级| 制服丝袜香蕉在线| 狂野欧美白嫩少妇大欣赏| 日本午夜av视频| 亚洲图色成人| 日韩av免费高清视频| 高清视频免费观看一区二区| 91精品国产国语对白视频| 高清黄色对白视频在线免费看 | 国产精品国产三级国产av玫瑰| 国产欧美亚洲国产| 麻豆乱淫一区二区| 午夜激情久久久久久久| 国产高清不卡午夜福利| 国产一级毛片在线| 色哟哟·www| 亚洲电影在线观看av| 另类精品久久| 国产精品三级大全|