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

    模擬氮沉降對(duì)貝加爾針茅草原土壤氮轉(zhuǎn)化微生物的影響

    2019-10-23 06:01:06劉紅梅張海芳張艷軍楊殿林
    關(guān)鍵詞:古菌銨態(tài)氮拷貝數(shù)

    劉紅梅,張海芳,秦 潔,王 慧,張艷軍,楊殿林

    (農(nóng)業(yè)農(nóng)村部環(huán)境保護(hù)科研監(jiān)測(cè)所,農(nóng)業(yè)農(nóng)村部產(chǎn)地環(huán)境污染防控重點(diǎn)實(shí)驗(yàn)室/天津市農(nóng)業(yè)環(huán)境與農(nóng)產(chǎn)品安全重點(diǎn)實(shí)驗(yàn)室,天津300191)

    大氣氮沉降增加呈現(xiàn)全球化趨勢(shì),我國(guó)是繼歐洲、北美之后的第3大氮沉降區(qū),大氣氮沉降遠(yuǎn)高于全球平均水平[1]。大氣氮沉降借助其對(duì)土壤碳固定[2]、植物氮素利用等的直接或間接作用[3],極大地干預(yù)了生態(tài)系統(tǒng)碳蓄積和氮素重新分配過程[4-5]。高氮沉降導(dǎo)致土壤氮轉(zhuǎn)化過程發(fā)生變化[6]、土壤可利用性氮含量增加[7]、氮循環(huán)微生物數(shù)量及組成發(fā)生變化[8],進(jìn)而對(duì)生態(tài)系統(tǒng)產(chǎn)生影響。

    土壤氮轉(zhuǎn)化主要過程包括生物固氮作用、氨化作用、硝化作用和反硝化作用,這些過程都由土壤微生物所驅(qū)動(dòng)。氮沉降進(jìn)入土壤后的一系列轉(zhuǎn)化離不開與氮素轉(zhuǎn)化的相關(guān)微生物參與。近年來(lái)的研究表明,施氮肥顯著影響固氮菌、硝化菌和反硝化菌群落結(jié)構(gòu)和豐度[8]。適量氮肥添加提高固氮菌豐度,促進(jìn)土壤微生物的固氮功能[9];高濃度無(wú)機(jī)氮添加抑制固氮菌生長(zhǎng)[10],從而抑制固氮微生物生長(zhǎng)。Ning等[11]研究發(fā)現(xiàn)氮相關(guān)功能基因?qū)Φ砑铀俾时憩F(xiàn)出不同的敏感性,氨氧化細(xì)菌AOB-amoA基因比氨氧化古菌AOA-amoA基因?qū)Φ砑痈用舾小M寥腊毖趸?xì)菌豐度隨氮添加量增加而增加,而氨氧化古菌豐度則無(wú)明顯變化[12]。氮肥種類和氮肥施用量會(huì)影響土壤氨氧化微生物組成和豐度[13]。受氮素限制的黃土高原鹽堿草地,施氮后顯著地提高了微生物生物量氮轉(zhuǎn)化速率[14]。在亞熱帶森林、溫帶森林中研究發(fā)現(xiàn),N2O排放、CH4吸收與AOA、AOB群落豐度分別呈現(xiàn)正、負(fù)相關(guān)關(guān)系,氨氧化菌群落動(dòng)態(tài)能夠解釋土壤CH4吸收和N2O排放之間的消長(zhǎng)作用[15]。

    目前,在氮沉降增加的背景下,針對(duì)干旱和半干旱草原土壤氮轉(zhuǎn)化功能基因豐度對(duì)氮沉降增加的響應(yīng)特征還不夠明確。為此,本研究通過施氮模擬氮沉降增加,研究不同氮添加水平對(duì)涉及土壤氮循環(huán)中的固氮nifH基因、氨氧化細(xì)菌AOB-amoA基因、氨氧化古菌AOA-amoA基因和反硝化基因nirK豐度的影響,通過分析氮轉(zhuǎn)化功能基因豐度與土壤理化因子之間的相關(guān)性關(guān)系,研究氮添加對(duì)微生物介導(dǎo)的固氮、氨氧化和反硝化過程的影響及其反饋,為深入認(rèn)識(shí)草地生態(tài)系統(tǒng)氮循環(huán)及其對(duì)氮沉降增加的響應(yīng)機(jī)制提供參考。

    1 材料與方法

    1.1 研究區(qū)域概況與樣地設(shè)置

    研究區(qū)域概況、試驗(yàn)樣地設(shè)置及試驗(yàn)前植被類型和基本理化性質(zhì)詳見參考文獻(xiàn)[16]。試驗(yàn)樣地位于內(nèi)蒙古貝加爾針茅草原(48°27′~48°35′N,119°35′~119°41′E)。于2010年開始進(jìn)行圍欄試驗(yàn),氮添加水平設(shè)置參考國(guó)際上同類研究試驗(yàn)[17-18]。氮添加處理設(shè)8個(gè)水平,依次為:0、15、30、50、100、150、200、300 kg N·hm-2·a-1,分別記為:對(duì)照N0,低氮(N15、N30和N50),高氮(N100、N150、N200和N300)。

    氮素為NH4NO3,每年6月中旬和7月中旬進(jìn)行氮添加。共設(shè)置8個(gè)處理水平,4次重復(fù)。小區(qū)面積8 m×8 m,小區(qū)間設(shè)置2 m隔離帶,重復(fù)間設(shè)置5 m隔離帶。

    1.2 樣品采集與處理

    土壤樣品采集于2015年8月中旬,用土鉆在各小區(qū)采集0~10 cm土壤樣品,每小區(qū)采集10個(gè)點(diǎn),混合均勻。去除土壤入侵物和根系,將土壤樣品分為兩份,一份保存于冰盒中帶回實(shí)驗(yàn)室,存于-20℃冰柜中;一份土壤樣品于室內(nèi)自然風(fēng)干。

    1.3 測(cè)定方法

    土壤理化性質(zhì)測(cè)定:土壤pH采用玻璃電極法(土水比1∶2.5),土壤總有機(jī)碳測(cè)定采用重鉻酸鉀外加熱法,土壤全氮用凱氏定氮法,土壤銨態(tài)氮和硝態(tài)氮含量采用氯化鉀溶液提取-流動(dòng)分析儀(QC8000)測(cè)定[19]。土壤微生物生物量氮采用氯仿熏蒸-K2SO4提取-TOC儀(Multi N/C 3000)測(cè)定[20]。土壤可溶性氮以測(cè)定土壤微生物生物量氮時(shí)未熏蒸土壤的全氮含量表示,土壤可溶性有機(jī)氮以土壤可溶性氮減去銨態(tài)氮和硝態(tài)氮的差值表示[21]。

    土壤DNA提取:稱取0.25 g在冰柜中保存的土壤樣品,用PowerSoil?DNA Isolationkit試劑盒,按照說(shuō)明書提取土壤總DNA。用Nanodrop2000超微量分光光度計(jì)和1%瓊脂糖電泳檢測(cè)DNA質(zhì)量。將通過質(zhì)量檢測(cè)的DNA樣品于-20℃冰箱保存待用。

    標(biāo)準(zhǔn)品配制:用目的基因引物進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體系為50 μL,包括上下游引物各1 μL,模板1 μL,2×Taq MasterMix 25 μL,滅菌超純水22 μL。擴(kuò)增條件為94℃5 min預(yù)變性后;94℃30 s,55℃30 s,72℃ 30 s,共30個(gè)循環(huán);最后1輪循環(huán)完成后再72℃延伸10 min。反應(yīng)完畢后,用1%瓊脂糖凝膠電泳檢查擴(kuò)增結(jié)果,目的片段進(jìn)行瓊脂糖凝膠回收。將回收PCR產(chǎn)物連接至pMD18-T載體,轉(zhuǎn)化至大腸桿菌DH5α感受態(tài)細(xì)胞中,經(jīng)Amp+、IPTG和X-gal的LB平板篩選陽(yáng)性克隆,送去生工生物有限公司進(jìn)行序列測(cè)定分析。

    標(biāo)準(zhǔn)曲線制作:通過做預(yù)實(shí)驗(yàn),將4種提取測(cè)序正確的質(zhì)粒標(biāo)準(zhǔn)品從101~105倍進(jìn)行10倍稀釋,每個(gè)梯度取2 μL做模板建立標(biāo)準(zhǔn)曲線。

    定量PCR檢測(cè):土壤固氮基因nifH、氨氧化細(xì)菌AOB、氨氧化古菌AOA和反硝化細(xì)菌定量擴(kuò)增引物見表1。擴(kuò)增反應(yīng)體系總體積為20 μL,包括2×GoT-aq?qPCR Master Mix 10 μL,10 μmol·L-1的上下游引物各0.5 μL,DNA模板2 μL(1~10 ng),滅菌超純水7 μL。將加好樣的96-PCR板置于熒光定量PCR儀進(jìn)行擴(kuò)增反應(yīng),每個(gè)樣品3個(gè)重復(fù)。擴(kuò)增反應(yīng)條件為95℃預(yù)變性30 s;95℃ 變性5 s,60℃退火40 s,72℃延伸30 s,40個(gè)循環(huán)。

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

    應(yīng)用SPSS 17.0統(tǒng)計(jì)軟件進(jìn)行單因素方差分析(One-way ANOVA)、Pearson相關(guān)性分析。

    2 結(jié)果與分析

    2.1 氮沉降對(duì)土壤化學(xué)性質(zhì)的影響

    連續(xù)6年不同氮添加處理下,土壤化學(xué)性質(zhì)變化見表2。高氮添加(N100、N150、N200和N300)土壤有機(jī)碳、硝態(tài)氮和微生物生物量氮含量高于或顯著高于低氮添加(N15、N30和N50)和對(duì)照N0。7個(gè)氮添加處理的土壤pH均低于或顯著低于對(duì)照N0,且高氮添加pH顯著低于低氮添加。土壤全氮含量在不同氮添加處理間無(wú)顯著差異。

    2.2 氮沉降對(duì)土壤氮轉(zhuǎn)化功能基因豐度的影響

    固氮微生物功能基因nifH豐度在各氮添加處理中的變化范圍在5.26×105~1.43×106copies·g-1soil之間(圖1A)。隨著氮添加水平的升高,nifH豐度表現(xiàn)為先升高后降低的趨勢(shì),其中N50處理是對(duì)照N0的2.39倍。N15、N30、N50、N100、N150和 N200處理nifH豐度高于或顯著高于對(duì)照N0。N300處理nifH豐度低于對(duì)照N0,但無(wú)顯著差異(P>0.05)。

    表1 定量PCR擴(kuò)增引物Table 1 Amplification primer of quantitative PCR

    表2 不同氮添加處理下土壤化學(xué)性質(zhì)Table 2 Soil chemical properties in different nitrogen addition treatments

    氨氧化細(xì)菌AOB-amoA基因豐度在各氮添加處理中的變化范圍為1.33×107~2.48×107copies·g-1soil之間(圖1B)。N15、N30和N50處理AOB-amoA基因豐度與對(duì)照相比無(wú)顯著差異(P>0.05),高氮添加(N100、N150、N200和N300)AOB-amoA基因豐度均顯著高于對(duì)照N0(P<0.05)。表明氮沉降增加使AOB-amoA基因豐度提高。氨氧化古菌AOA-amoA基因豐度在各氮添加處理中的變化范圍為2.56×106~3.97×106copies·g-1soil之間(圖1C)。N30處理AOA-amoA基因豐度最高,但與對(duì)照相比無(wú)顯著差異,其他6個(gè)氮添加處理均低于對(duì)照N0。土壤中AOB的基因拷貝數(shù)遠(yuǎn)高于AOA基因拷貝數(shù),且高氮添加顯著提高了AOB/AOA比值(圖1E)。

    圖1 不同氮處理的功能基因豐度Figure 1 Functional genes in different nitrogen treatments

    土壤反硝化功能基因nirK豐度在不同氮添加處理中的變化范圍為 2.88×107~7.53×107copies·g-1soil之間(圖1D)。N15處理nirK基因豐度低于對(duì)照N0,但無(wú)顯著差異(P>0.05);N30、N50和N100處理nirK基因豐度均顯著高于對(duì)照N0(P<0.05);高氮添加處理(N150、N200和N300)nirK基因豐度顯著低于對(duì)照N0(P<0.05),顯著提高了氨氧化/反硝化豐度比值(圖1F)。

    表3 功能基因拷貝數(shù)與化學(xué)性質(zhì)之間的相關(guān)性分析Table 3 Correlation analysis between functional gene copies and soil chemical properties

    2.3 土壤氮素轉(zhuǎn)化功能基因豐度與土壤因子的相關(guān)性

    土壤氮轉(zhuǎn)化功能基因拷貝數(shù)與土壤化學(xué)性質(zhì)相關(guān)性分析見表3。土壤nifH基因拷貝數(shù)與總有機(jī)碳、硝態(tài)氮和微生物生物量氮含量呈顯著負(fù)相關(guān)(P<0.05),與土壤銨態(tài)氮含量呈極顯著正相關(guān)(P<0.01)。土壤AOB-amoA基因拷貝數(shù)與有機(jī)碳、硝態(tài)氮、可溶性有機(jī)氮和微生物生物量氮含量呈極顯著正相關(guān)(P<0.01),與土壤pH呈極顯著負(fù)相關(guān)(P<0.01)。土壤AOA-amoA基因拷貝數(shù)與土壤pH呈極顯著正相關(guān)(P<0.01),與土壤有機(jī)碳、硝態(tài)氮和微生物生物量氮含量呈顯著負(fù)相關(guān)(P<0.05)。土壤nirK基因拷貝數(shù)與土壤pH和銨態(tài)氮含量呈極顯著正相關(guān)關(guān)系(P<0.01),與土壤有機(jī)碳、硝態(tài)氮、可溶性有機(jī)氮和微生物生物量氮含量呈顯著負(fù)相關(guān)關(guān)系(P<0.05)。AOB/AOA與土壤有機(jī)碳、硝態(tài)氮、可溶性有機(jī)氮和微生物生物量氮呈極顯著正相關(guān)(P<0.01),與pH呈極顯著負(fù)相關(guān)(P<0.01)。氨氧化/反硝化與土壤有機(jī)碳、硝態(tài)氮、可溶性有機(jī)氮和微生物生物量氮呈極顯著正相關(guān)(P<0.01),與pH和銨態(tài)氮含量呈顯著負(fù)相關(guān)(P<0.05)。

    3 討論

    3.1 土壤固氮功能基因?qū)Φ两淀憫?yīng)的綜合分析

    土壤固氮微生物是一類重要的氮轉(zhuǎn)化功能微生物,其能夠?qū)⒖諝庵械牡€原成氨供給植物吸收利用。固氮微生物數(shù)量與輸入土壤的外源氮素密切相關(guān),施氮量過高會(huì)降低固氮菌的豐度并抑制固氮酶的活性,并對(duì)固氮菌群落組成產(chǎn)生顯著影響[26]。本研究中,隨著氮添加水平的升高,nifH豐度表現(xiàn)為先升高后降低趨勢(shì),N300處理降低了nifH基因相對(duì)豐度,其他6個(gè)氮添加處理均高于或顯著高于對(duì)照N0。即低于200 kg N·hm-2·a-1時(shí)更有利于固氮菌的生長(zhǎng),而高于300 kg N·hm-2·a-1時(shí)反而抑制其生長(zhǎng)。這與劉彩霞等[27]得出的氮沉降(銨態(tài)氮)高于60 kg N·hm-2·a-1抑制杉木林固氮功能微生物生長(zhǎng)結(jié)論一致。Zhang等[7]研究也表明,高濃度無(wú)機(jī)氮添加抑制固氮菌生長(zhǎng)。前人研究認(rèn)為,氮肥輸入降低固氮菌的豐度是由于土壤酸化或高氮含量造成的[28]。

    土壤固氮菌豐度與土壤有機(jī)碳含量、pH、有效磷含量等其他化學(xué)因子密切相關(guān)。本研究相關(guān)分析表明,nifH基因豐度與土壤有機(jī)碳含量呈顯著負(fù)相關(guān)性。由于生物固氮過程需要消耗大量能量,而該過程在很大程度上依賴于土壤中的有機(jī)碳含量。高氮添加雖然提高了土壤有機(jī)碳含量,但高氮添加導(dǎo)致高的速效氮含量反而會(huì)抑制土壤的生物固氮過程。本研究表明,nifH基因豐度與銨態(tài)氮含量呈極顯著正相關(guān),與土壤硝態(tài)氮含量呈顯著負(fù)相關(guān)。表明氮添加量增加引起土壤可利用性氮含量變化是引起固氮功能基因豐度發(fā)生變化的重要驅(qū)動(dòng)因素。Wang等[29]研究表明,長(zhǎng)期施肥引起小麥生長(zhǎng)季土壤銨態(tài)氮含量改變,顯著影響固氮菌豐度,不利于土壤微生物的固氮功能。本研究中,低氮添加促進(jìn)土壤的生物固氮過程,但過高氮添加反而會(huì)抑制土壤的生物固氮過程。表明未來(lái)連續(xù)高氮沉降可能不利于貝加爾針茅草原土壤微生物的固氮功能。

    3.2土壤氨氧化功能基因?qū)Φ两淀憫?yīng)的綜合分析

    硝化作用是氮素轉(zhuǎn)化的關(guān)鍵過程,由含氨單加氧酶基因(amoA)、氨氧化細(xì)菌(AOB)和氨氧化古菌(AOA)共同驅(qū)動(dòng)氨氧化過程是硝化作用的限速步驟,其數(shù)量和結(jié)構(gòu)組成變化可反映土壤氮素變化[30]。目前氮肥或氮添加對(duì)土壤硝化作用的影響研究大多是針對(duì)這一步驟開展的。一般情況下,氨是微生物進(jìn)行硝化作用的直接底物,酸性土壤中銨根離子分解處理的氨較少,從而抑制AOB生長(zhǎng)[31]。本研究,低氮添加(N15、N30和N50)對(duì)氨氧化細(xì)菌AOB相對(duì)豐度無(wú)顯著影響,而高氮添加顯著提高了氨氧化細(xì)菌豐度。說(shuō)明連續(xù)6 a高氮添加,促進(jìn)了AOB的生長(zhǎng)繁殖。這與Li等[32]對(duì)內(nèi)蒙古草原10 a氮添加試驗(yàn)AOB基因豐度變化結(jié)果一致。Ai等[33]在對(duì)小麥玉米輪作潮土研究表明,氮肥添加可顯著提高AOB的數(shù)量。其可能原因是土壤中存在對(duì)NH+4有較高親和力的未知的AOB[34]或者是某些AOB的生長(zhǎng)不以NH+4作為唯一底物[35]。與AOB豐度變化情況不同,除N30處理與對(duì)照相比無(wú)顯著差異外,其余6種氮添加處理均降低了土壤氨氧化古菌基因豐度,說(shuō)明長(zhǎng)期氮沉降增加會(huì)導(dǎo)致土壤氨氧化古菌豐度降低。長(zhǎng)期施氮肥導(dǎo)致南方紅壤中氨氧化古菌豐度降低[36],這表明盡管草原暗栗鈣土和南方紅壤在pH、有機(jī)質(zhì)和碳氮水平上存在較大差異,但長(zhǎng)期氮添加和長(zhǎng)期施用氮肥對(duì)這兩種土壤中氨氧化古菌的影響存在一致性。在本研究中,各氮添加處理土壤AOB-amoA基因拷貝數(shù)比AOA-amoA基因拷貝數(shù)高一個(gè)數(shù)量級(jí),可見本試驗(yàn)環(huán)境下更有利于AOB生長(zhǎng)。這與Yao等[37]在茶園酸性土壤AOA主導(dǎo)氨氧化作用的研究結(jié)果不一致。前人研究表明,AOB和AOA生長(zhǎng)對(duì)土壤含氮量響應(yīng)不同,AOB對(duì)高氮環(huán)境較為適應(yīng),而AOA傾向于低氮環(huán)境[38]。本試驗(yàn)中高氮添加處理土壤硝態(tài)氮含量顯著高于低氮添加和對(duì)照處理,從而促進(jìn)了AOB生長(zhǎng)。但對(duì)于這一原因推測(cè)的合理性還需要進(jìn)一步深入研究。

    大量研究表明,多種土壤環(huán)境因子,包括pH、溫度、土壤類型和利用方式等都會(huì)影響氨氧化微生物的群落組成和數(shù)量,從而使氮添加對(duì)土壤氨氧化微生物影響存在很大不確定性。本研究氨氧化細(xì)菌和氨氧化古菌與土壤理化因子的相關(guān)分析表明,影響氨氧化功能基因相對(duì)豐度的土壤理化因子為pH、有機(jī)碳、硝態(tài)氮和微生物生物量氮含量。這與Hayden等[39]對(duì)澳大利亞農(nóng)業(yè)與草地氮轉(zhuǎn)化功能基因驅(qū)動(dòng)因子研究結(jié)果一致。相關(guān)分析表明,AOB基因拷貝數(shù)與土壤有機(jī)碳含量呈極顯著正相關(guān),說(shuō)明連續(xù)6 a氮添加提高了貝加爾針茅草原土壤有機(jī)碳含量是促進(jìn)氨氧化細(xì)菌生長(zhǎng)繁殖的重要因素。一些研究表明,氨氧化細(xì)菌適合在氮素含量豐富的土壤中生存,且其數(shù)量與土壤pH呈顯著負(fù)相關(guān)[14,40]。本研究中,氨氧化細(xì)菌基因拷貝數(shù)與土壤pH呈極顯著負(fù)相關(guān)。這與其研究結(jié)果一致。原因可能是,隨著氮添加水平的增加,土壤pH降低,使得土壤中耐酸性AOB數(shù)量增多[41]。氨氧化古菌與土壤化學(xué)因子的相關(guān)性分析表明,氨氧化古菌基因拷貝數(shù)與土壤pH呈極顯著正相關(guān),說(shuō)明貝加爾針茅草原土壤中AOA基因豐度隨著pH降低而顯著降低。高氮添加增加了土壤中硝態(tài)氮含量,相應(yīng)的AOB基因拷貝數(shù)增加,而AOA基因拷貝數(shù)則減少,AOB/AOA比值提高,且隨著氮添加水平的增加,提高的幅度明顯增加。表明連續(xù)6 a高氮添加,改變了氨氧化微生物相對(duì)豐度,氨氧化微生物對(duì)氨氧化的相對(duì)貢獻(xiàn)也可能發(fā)生改變。但氨氧化微生物對(duì)氨氧化作用的相對(duì)貢獻(xiàn)率是否發(fā)生改變,還需采用同位素示蹤技術(shù)進(jìn)一步研究。相關(guān)分析表明,AOA基因拷貝數(shù)與土壤硝態(tài)氮含量呈極顯著負(fù)相關(guān)關(guān)系,這與Di等[38]的AOA適合于低氮環(huán)境的研究結(jié)論一致。本研究中長(zhǎng)期連續(xù)氮添加導(dǎo)致的土壤化學(xué)性質(zhì)變化,尤其是土壤pH降低導(dǎo)致土壤中AOA豐度明顯降低,可能降低土壤中AOA的氨氧化功能,進(jìn)一步影響貝加爾針茅草原土壤的硝化作用。

    3.3 土壤反硝化功能基因?qū)Φ两淀憫?yīng)的綜合分析

    反硝化作用是將硝酸根還原為N2O或N2的過程,是由許多厭氧或兼性厭氧微生物參與的一系列酶催化過程。這些酶主要由功能基因硝酸還原酶narG/napA、亞硝酸還原酶nirK/nirS、一氧化氮還原酶norB和氧化亞氮還原酶nosZ編碼。不同反硝化功能基因型對(duì)環(huán)境因子的影響存在較大差異。其中,亞硝酸還原酶nirK/nirS催化亞硝酸鹽還原成NO是微生物驅(qū)動(dòng)反硝化作用關(guān)鍵步驟。已有研究表明,nirK型反硝化細(xì)菌對(duì)施肥種類和施肥量都十分敏感[42]。連續(xù)6 a氮添加顯著影響了nirK相對(duì)豐度,N30、N50和N100處理顯著提高了nirK相對(duì)豐度,N150、N200和N300處理顯著降低了nirK相對(duì)豐度。表明,高氮添加(N150、N200和N300)對(duì)nirK豐度有明顯的抑制作用。

    環(huán)境因子如土壤pH、溫度、有機(jī)碳、硝態(tài)氮含量是影響反硝化作用的重要因素。本研究相關(guān)分析表明,nirk相對(duì)豐度與土壤pH、有機(jī)碳、硝態(tài)氮、銨態(tài)氮、可溶性有機(jī)氮和微生物生物量氮含量均呈顯著相關(guān)性(表3)。Xie等[43]研究認(rèn)為土壤有機(jī)碳和硝酸鹽含量是驅(qū)動(dòng)青藏高原草甸反硝化功能基因nirK豐度的重要土壤因子。Jahangir等[44]研究表明,反硝化活性的高低與提供給細(xì)菌生長(zhǎng)的有機(jī)質(zhì)的質(zhì)量和數(shù)量密切相關(guān)。土壤反硝化作用最適pH范圍為7.0~8.0,在此范圍內(nèi)反硝化作用最好[45]。在本研究中,高氮添加(N150、N200和N300)處理pH顯著低于對(duì)照N0,降低了土壤pH,抑制了反硝化功能基因nirK豐度,并由此可能對(duì)反硝化作用產(chǎn)生影響。一般情況下,土壤AOB基因豐度增加會(huì)提高反硝化作物的底物硝酸鹽含量,進(jìn)而促進(jìn)反硝化微生物驅(qū)動(dòng)的反硝化作用,造成氮素流失。但本研究中,高氮(N150、N200和N300)顯著降低了反硝化細(xì)菌nirK基因豐度,而硝態(tài)氮含量卻顯著提高。說(shuō)明高氮添加驅(qū)動(dòng)的土壤氨氧化作用高于驅(qū)動(dòng)的土壤的反硝化作用。土壤nirK基因豐度與土壤硝態(tài)氮含量呈極顯著負(fù)相關(guān)驗(yàn)證了這一觀點(diǎn)。

    上述研究表明,貝加爾針茅草原土壤反硝化活性受土壤化學(xué)性質(zhì)特別是pH、有機(jī)碳、硝態(tài)氮、銨態(tài)氮、可溶性有機(jī)氮和微生物生物量氮含量的顯著影響。此外,由于不同反硝化過程所涉及的酶不同,同時(shí)攜帶不同反硝化基因的微生物種群對(duì)環(huán)境因子響應(yīng)的差異,因此,揭示不同氮添加水平對(duì)反硝化功能微生物基因豐度的影響,需要對(duì)不同類型的反硝化微生物開展深入研究。

    4 結(jié)論

    (1)隨著氮添加水平的增加,固氮細(xì)菌nifH基因豐度呈現(xiàn)先升高后降低趨勢(shì)。低于200 kg N·hm-2·a-1時(shí)促進(jìn)固氮菌生長(zhǎng)。影響固氮菌nifH主要環(huán)境因子是土壤有機(jī)碳、銨態(tài)氮、硝態(tài)氮和微生物生物量氮含量。

    (2)高氮添加(N100、N150、N200和N300)提高了AOB基因豐度,降低了AOA基因豐度,并顯著提高了AOB/AOA,表明氮添加量增加可能使AOB在貝加爾針茅草原土壤硝化過程中占主導(dǎo)作用。影響AOB、AOA基因豐度的主要環(huán)境因子是土壤pH、有機(jī)碳、硝態(tài)氮和微生物生物量氮含量。

    (3)高氮添加(N150、N200和N300)降低了反硝化細(xì)菌nirK基因豐度。氮添加量的增加促進(jìn)了AOB主導(dǎo)的氨氧化過程,而反硝化微生物nirK豐度的降低提高了氨氧化產(chǎn)物的積累,繼而提高了土壤中的硝態(tài)氮含量。影響反硝化細(xì)菌nirK的主要環(huán)境因子是土壤pH、有機(jī)碳、硝態(tài)氮、銨態(tài)氮、可溶性有機(jī)氮和微生物生物量氮含量。

    猜你喜歡
    古菌銨態(tài)氮拷貝數(shù)
    不同pH和氧氣條件下土壤古菌與海洋古菌的競(jìng)爭(zhēng)適應(yīng)機(jī)制*
    變油為氣,“榨干”廢棄油田
    大自然探索(2022年5期)2022-07-11 03:10:33
    海洋古菌
    線粒體DNA拷貝數(shù)變異機(jī)制及疾病預(yù)測(cè)價(jià)值分析
    不同質(zhì)地土壤銨態(tài)氮吸附/解吸特征
    不同鹽堿化土壤對(duì)NH+4吸附特性研究
    胎兒染色體組拷貝數(shù)變異與產(chǎn)前超聲異常的相關(guān)性分析
    有機(jī)質(zhì)對(duì)城市污染河道沉積物銨態(tài)氮吸附-解吸的影響*
    DNA序列拷貝數(shù)變化決定黃瓜性別
    東太平洋海隆深海熱液區(qū)沉積物古菌多樣性分析
    欧美 日韩 精品 国产| 亚洲欧美色中文字幕在线| 综合色丁香网| 又大又黄又爽视频免费| 22中文网久久字幕| 免费看光身美女| 亚洲第一av免费看| 欧美 日韩 精品 国产| 99热全是精品| 老司机影院成人| 午夜老司机福利剧场| 欧美国产精品一级二级三级| 久久久久精品久久久久真实原创| 一级,二级,三级黄色视频| 99视频精品全部免费 在线| 大香蕉久久网| 赤兔流量卡办理| 日本欧美视频一区| 纯流量卡能插随身wifi吗| 国产精品一区www在线观看| 免费观看无遮挡的男女| 久久人人爽人人爽人人片va| www.色视频.com| 亚洲av中文av极速乱| 亚洲欧美一区二区三区黑人 | videos熟女内射| 久久久久久人妻| 99热全是精品| 国产精品熟女久久久久浪| 中文字幕av电影在线播放| 美女xxoo啪啪120秒动态图| 色吧在线观看| av福利片在线| 男女边摸边吃奶| 欧美日韩综合久久久久久| 蜜桃在线观看..| 日韩制服丝袜自拍偷拍| 亚洲国产精品成人久久小说| 在线观看免费高清a一片| 久久ye,这里只有精品| av免费观看日本| 亚洲av国产av综合av卡| 最近最新中文字幕大全免费视频 | 国产午夜精品一二区理论片| 丝袜喷水一区| 丝袜在线中文字幕| 老司机亚洲免费影院| 久久人人爽av亚洲精品天堂| 国产麻豆69| 中文欧美无线码| 22中文网久久字幕| 永久网站在线| 女性生殖器流出的白浆| 尾随美女入室| h视频一区二区三区| 新久久久久国产一级毛片| 国产深夜福利视频在线观看| 精品人妻一区二区三区麻豆| 丰满饥渴人妻一区二区三| 精品国产国语对白av| 两个人看的免费小视频| 成人国产av品久久久| 精品人妻熟女毛片av久久网站| 久久狼人影院| 国产在线一区二区三区精| 韩国精品一区二区三区 | 十八禁网站网址无遮挡| 精品人妻熟女毛片av久久网站| 成人免费观看视频高清| 尾随美女入室| 深夜精品福利| 久久精品熟女亚洲av麻豆精品| 国产亚洲精品第一综合不卡 | 又大又黄又爽视频免费| 最新的欧美精品一区二区| 成年人免费黄色播放视频| 中国美白少妇内射xxxbb| 另类亚洲欧美激情| 午夜福利乱码中文字幕| 亚洲精品久久成人aⅴ小说| 伊人久久国产一区二区| 999精品在线视频| 午夜老司机福利剧场| 欧美日韩av久久| 成人综合一区亚洲| 久久免费观看电影| 国产高清国产精品国产三级| 成人二区视频| 日韩精品免费视频一区二区三区 | 亚洲综合色惰| 9色porny在线观看| 亚洲精品中文字幕在线视频| www.av在线官网国产| 91aial.com中文字幕在线观看| 中国国产av一级| 国产黄色视频一区二区在线观看| 少妇的逼水好多| 丝袜人妻中文字幕| 观看美女的网站| 欧美亚洲日本最大视频资源| 岛国毛片在线播放| av免费在线看不卡| av女优亚洲男人天堂| 男女无遮挡免费网站观看| 国产免费福利视频在线观看| 精品一区在线观看国产| 中文字幕制服av| 欧美激情极品国产一区二区三区 | 熟妇人妻不卡中文字幕| 国产精品一区二区在线观看99| 亚洲国产av新网站| 精品卡一卡二卡四卡免费| 三上悠亚av全集在线观看| 九草在线视频观看| 日本vs欧美在线观看视频| 精品卡一卡二卡四卡免费| videos熟女内射| 在线天堂最新版资源| 亚洲一级一片aⅴ在线观看| av不卡在线播放| 日韩视频在线欧美| 欧美日韩成人在线一区二区| 亚洲欧美日韩卡通动漫| 在线观看人妻少妇| 一级片免费观看大全| 国产亚洲最大av| 亚洲五月色婷婷综合| av视频免费观看在线观看| 超色免费av| 91在线精品国自产拍蜜月| 久久久久久人人人人人| 99视频精品全部免费 在线| 少妇的丰满在线观看| 免费日韩欧美在线观看| 亚洲性久久影院| 观看av在线不卡| 精品少妇内射三级| 久久精品熟女亚洲av麻豆精品| 人成视频在线观看免费观看| www.熟女人妻精品国产 | 少妇被粗大猛烈的视频| 午夜91福利影院| 久久久久久久亚洲中文字幕| 国产精品成人在线| 久久韩国三级中文字幕| 欧美97在线视频| 男的添女的下面高潮视频| 人人澡人人妻人| 久久久久久人人人人人| 午夜福利在线观看免费完整高清在| 又粗又硬又长又爽又黄的视频| 久久精品夜色国产| 国产精品.久久久| 亚洲av.av天堂| av不卡在线播放| 午夜免费男女啪啪视频观看| 亚洲国产精品专区欧美| 欧美成人午夜精品| 国产av精品麻豆| 久久久国产一区二区| 麻豆乱淫一区二区| 老司机影院毛片| 亚洲国产欧美在线一区| 99热全是精品| 中文字幕另类日韩欧美亚洲嫩草| 久久影院123| 久热久热在线精品观看| 成人毛片a级毛片在线播放| 丰满迷人的少妇在线观看| 久久久久久久国产电影| 亚洲精品aⅴ在线观看| 国产日韩一区二区三区精品不卡| 新久久久久国产一级毛片| 人成视频在线观看免费观看| 91精品伊人久久大香线蕉| 久久久久久久亚洲中文字幕| 久久99蜜桃精品久久| 啦啦啦视频在线资源免费观看| 飞空精品影院首页| 最近中文字幕2019免费版| 免费av中文字幕在线| 热re99久久精品国产66热6| 久久精品人人爽人人爽视色| 精品久久蜜臀av无| 欧美日韩亚洲高清精品| 永久免费av网站大全| 人人澡人人妻人| 国产乱人偷精品视频| 欧美精品一区二区大全| 啦啦啦视频在线资源免费观看| 亚洲国产毛片av蜜桃av| 日本av免费视频播放| 一级毛片我不卡| av视频免费观看在线观看| 日本午夜av视频| 99九九在线精品视频| 男女无遮挡免费网站观看| 国产精品一国产av| 春色校园在线视频观看| 91精品三级在线观看| 少妇人妻久久综合中文| 国产精品国产三级国产专区5o| 久久久久久久久久人人人人人人| 欧美精品av麻豆av| 久久久久久久亚洲中文字幕| 99热国产这里只有精品6| 亚洲四区av| 色5月婷婷丁香| 精品一区二区免费观看| 美女视频免费永久观看网站| 边亲边吃奶的免费视频| 久久综合国产亚洲精品| 夜夜爽夜夜爽视频| 欧美精品一区二区免费开放| 日本欧美国产在线视频| 色94色欧美一区二区| √禁漫天堂资源中文www| 蜜桃国产av成人99| 精品熟女少妇av免费看| 91在线精品国自产拍蜜月| 免费少妇av软件| 国产精品偷伦视频观看了| 大片免费播放器 马上看| 免费少妇av软件| 亚洲国产精品一区二区三区在线| 大陆偷拍与自拍| 九九在线视频观看精品| 国产日韩一区二区三区精品不卡| 又黄又爽又刺激的免费视频.| 22中文网久久字幕| 国产男人的电影天堂91| 最近2019中文字幕mv第一页| 黄色一级大片看看| 国产视频首页在线观看| 黄色 视频免费看| 国产精品一区二区在线观看99| 黑人猛操日本美女一级片| 国产精品国产三级国产专区5o| 高清视频免费观看一区二区| 成人国语在线视频| 亚洲精品中文字幕在线视频| 中文字幕人妻丝袜制服| 久久热在线av| 久久精品国产鲁丝片午夜精品| 美女脱内裤让男人舔精品视频| 久久人人爽人人爽人人片va| 三上悠亚av全集在线观看| 精品国产露脸久久av麻豆| 久久人人爽人人爽人人片va| 99热这里只有是精品在线观看| 国产成人精品在线电影| av免费在线看不卡| 狂野欧美激情性bbbbbb| 五月开心婷婷网| 日韩三级伦理在线观看| 麻豆乱淫一区二区| 香蕉丝袜av| 丰满乱子伦码专区| 亚洲三级黄色毛片| av福利片在线| av女优亚洲男人天堂| 国产白丝娇喘喷水9色精品| 久久久国产一区二区| 男人舔女人的私密视频| 亚洲精品aⅴ在线观看| 国产综合精华液| 亚洲第一av免费看| 国产一区二区三区av在线| 丝袜美足系列| 成年动漫av网址| 日本与韩国留学比较| 嫩草影院入口| 久久鲁丝午夜福利片| 一级毛片 在线播放| 亚洲欧美一区二区三区国产| 人人妻人人爽人人添夜夜欢视频| 女人精品久久久久毛片| 亚洲精品456在线播放app| 国产熟女欧美一区二区| 青春草亚洲视频在线观看| 国产高清三级在线| 日本av免费视频播放| 高清欧美精品videossex| 亚洲,欧美,日韩| 精品一区二区三区四区五区乱码 | 久久精品国产亚洲av天美| 色网站视频免费| 日韩制服骚丝袜av| 曰老女人黄片| 久久精品人人爽人人爽视色| 免费av中文字幕在线| 国产亚洲精品第一综合不卡 | 麻豆乱淫一区二区| 久久久国产精品麻豆| 中文精品一卡2卡3卡4更新| 中文字幕制服av| 大陆偷拍与自拍| 国产日韩一区二区三区精品不卡| 午夜视频国产福利| 女人精品久久久久毛片| 欧美亚洲日本最大视频资源| 日产精品乱码卡一卡2卡三| 国产精品久久久久久av不卡| 人人妻人人爽人人添夜夜欢视频| 2018国产大陆天天弄谢| 亚洲国产毛片av蜜桃av| 久久精品久久久久久噜噜老黄| 老司机亚洲免费影院| 最近中文字幕2019免费版| 久久ye,这里只有精品| 乱人伦中国视频| 亚洲美女黄色视频免费看| 欧美日韩国产mv在线观看视频| 久久精品熟女亚洲av麻豆精品| 亚洲伊人久久精品综合| 久久精品久久久久久久性| 狠狠婷婷综合久久久久久88av| 在线观看免费高清a一片| 青春草视频在线免费观看| 69精品国产乱码久久久| 国产一级毛片在线| 男女下面插进去视频免费观看 | 国产乱人偷精品视频| 日韩人妻精品一区2区三区| 大香蕉97超碰在线| 欧美最新免费一区二区三区| 久久国产亚洲av麻豆专区| 成年动漫av网址| 日本午夜av视频| 99热6这里只有精品| 久久久久久久久久久免费av| 一级片'在线观看视频| av在线老鸭窝| 26uuu在线亚洲综合色| 人妻系列 视频| 熟女电影av网| 一二三四中文在线观看免费高清| 亚洲精华国产精华液的使用体验| 十分钟在线观看高清视频www| 老司机影院成人| 欧美激情极品国产一区二区三区 | 国产精品一二三区在线看| 久久久久久久久久久久大奶| 卡戴珊不雅视频在线播放| 日韩精品免费视频一区二区三区 | 精品国产国语对白av| av有码第一页| av天堂久久9| 免费少妇av软件| 一区在线观看完整版| 亚洲精品第二区| 欧美最新免费一区二区三区| 成人18禁高潮啪啪吃奶动态图| 人妻 亚洲 视频| 免费在线观看黄色视频的| 水蜜桃什么品种好| 欧美精品人与动牲交sv欧美| 人妻人人澡人人爽人人| 欧美日韩一区二区视频在线观看视频在线| 9色porny在线观看| 免费人成在线观看视频色| 搡女人真爽免费视频火全软件| 97人妻天天添夜夜摸| 色网站视频免费| 国产免费一区二区三区四区乱码| 国产成人精品一,二区| 亚洲精品久久久久久婷婷小说| 中文字幕精品免费在线观看视频 | 高清av免费在线| 男人舔女人的私密视频| av免费观看日本| 亚洲欧美一区二区三区国产| 中文字幕人妻丝袜制服| 老司机影院成人| 亚洲人成77777在线视频| 黄色毛片三级朝国网站| 亚洲,欧美,日韩| 亚洲情色 制服丝袜| 日本欧美视频一区| 亚洲av成人精品一二三区| 黄色 视频免费看| 午夜福利在线观看免费完整高清在| 日本爱情动作片www.在线观看| 国产午夜精品一二区理论片| 黄片播放在线免费| 欧美激情 高清一区二区三区| 777米奇影视久久| 久久国产精品男人的天堂亚洲 | 一级黄片播放器| 寂寞人妻少妇视频99o| videossex国产| 久久久久久人人人人人| av在线观看视频网站免费| 99九九在线精品视频| 日韩成人伦理影院| 久久久国产一区二区| 久久精品熟女亚洲av麻豆精品| 亚洲国产av新网站| 国产色爽女视频免费观看| 国产极品粉嫩免费观看在线| 香蕉精品网在线| 日本色播在线视频| 天堂俺去俺来也www色官网| 亚洲第一区二区三区不卡| 最近手机中文字幕大全| 少妇高潮的动态图| 咕卡用的链子| 女人久久www免费人成看片| 美女xxoo啪啪120秒动态图| 高清视频免费观看一区二区| 激情五月婷婷亚洲| 丝袜喷水一区| 免费久久久久久久精品成人欧美视频 | xxxhd国产人妻xxx| 另类亚洲欧美激情| 亚洲精品国产色婷婷电影| 国产又色又爽无遮挡免| 免费久久久久久久精品成人欧美视频 | 伊人久久国产一区二区| 亚洲精品色激情综合| 亚洲国产精品国产精品| 看非洲黑人一级黄片| 国产欧美日韩一区二区三区在线| 91aial.com中文字幕在线观看| 亚洲精品456在线播放app| 亚洲激情五月婷婷啪啪| 久久国内精品自在自线图片| 成人亚洲欧美一区二区av| 丰满少妇做爰视频| 亚洲中文av在线| 精品卡一卡二卡四卡免费| 国产一区二区三区综合在线观看 | 国产探花极品一区二区| www.熟女人妻精品国产 | 日韩不卡一区二区三区视频在线| 欧美成人午夜精品| av.在线天堂| a级毛片黄视频| av.在线天堂| 久久久久精品性色| 欧美xxⅹ黑人| 99热6这里只有精品| 成人毛片60女人毛片免费| 最后的刺客免费高清国语| 制服丝袜香蕉在线| 美国免费a级毛片| 天堂中文最新版在线下载| videossex国产| 久久久国产精品麻豆| 国语对白做爰xxxⅹ性视频网站| 99国产精品免费福利视频| 黄色毛片三级朝国网站| 97超碰精品成人国产| 国产视频首页在线观看| 久久毛片免费看一区二区三区| 中国国产av一级| 久久久欧美国产精品| 亚洲av福利一区| 97在线人人人人妻| 各种免费的搞黄视频| 免费人成在线观看视频色| 精品久久蜜臀av无| kizo精华| 两个人看的免费小视频| 观看av在线不卡| 美女国产高潮福利片在线看| 婷婷色麻豆天堂久久| 久久精品夜色国产| 日韩av不卡免费在线播放| 建设人人有责人人尽责人人享有的| 亚洲美女黄色视频免费看| 热re99久久精品国产66热6| 最新中文字幕久久久久| 久久人人97超碰香蕉20202| 国产精品秋霞免费鲁丝片| 国产一区有黄有色的免费视频| 亚洲综合色惰| 久久精品aⅴ一区二区三区四区 | 最近最新中文字幕免费大全7| 国产爽快片一区二区三区| 久久99精品国语久久久| 热re99久久精品国产66热6| 寂寞人妻少妇视频99o| 国产黄色视频一区二区在线观看| 男的添女的下面高潮视频| 久久久精品区二区三区| 日本av免费视频播放| 国产无遮挡羞羞视频在线观看| 亚洲欧美成人综合另类久久久| 久久久久久久大尺度免费视频| 国产国语露脸激情在线看| 久久国产精品男人的天堂亚洲 | 又黄又爽又刺激的免费视频.| 18禁裸乳无遮挡动漫免费视频| 国产又色又爽无遮挡免| 9色porny在线观看| 九九在线视频观看精品| 2021少妇久久久久久久久久久| 蜜臀久久99精品久久宅男| 一边摸一边做爽爽视频免费| 免费人妻精品一区二区三区视频| 精品99又大又爽又粗少妇毛片| 99热国产这里只有精品6| 国产成人欧美| 日本vs欧美在线观看视频| videos熟女内射| 亚洲国产精品专区欧美| 国产又爽黄色视频| 亚洲综合色网址| 久久精品熟女亚洲av麻豆精品| 亚洲国产欧美在线一区| 久久久久精品人妻al黑| 久久99热这里只频精品6学生| 国产av精品麻豆| 一本久久精品| 国产女主播在线喷水免费视频网站| 激情视频va一区二区三区| 黄片无遮挡物在线观看| 亚洲欧美色中文字幕在线| 亚洲欧洲日产国产| 国产又色又爽无遮挡免| 久久人人爽av亚洲精品天堂| 亚洲丝袜综合中文字幕| 黑人猛操日本美女一级片| 深夜精品福利| 黄色毛片三级朝国网站| 人人妻人人澡人人爽人人夜夜| 国产精品一区www在线观看| 亚洲一区二区三区欧美精品| 亚洲av在线观看美女高潮| 少妇人妻久久综合中文| a级片在线免费高清观看视频| 青青草视频在线视频观看| 国产有黄有色有爽视频| 男人爽女人下面视频在线观看| 永久网站在线| 曰老女人黄片| 欧美日韩国产mv在线观看视频| 色94色欧美一区二区| 人妻人人澡人人爽人人| 成人手机av| av线在线观看网站| 久久久久人妻精品一区果冻| 国产av精品麻豆| 国产亚洲精品久久久com| 如日韩欧美国产精品一区二区三区| 中文字幕亚洲精品专区| 国产国拍精品亚洲av在线观看| 人妻 亚洲 视频| 18禁在线无遮挡免费观看视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品日本国产第一区| 亚洲av中文av极速乱| av在线app专区| 大片电影免费在线观看免费| 黄色配什么色好看| 午夜福利影视在线免费观看| 在线观看免费日韩欧美大片| 香蕉精品网在线| 国产男女内射视频| 97超碰精品成人国产| 一区二区三区四区激情视频| 亚洲一码二码三码区别大吗| 在线免费观看不下载黄p国产| 亚洲av男天堂| 80岁老熟妇乱子伦牲交| 久久韩国三级中文字幕| 亚洲国产精品一区二区三区在线| 97人妻天天添夜夜摸| 精品久久久久久电影网| 亚洲国产毛片av蜜桃av| av在线播放精品| 熟女电影av网| 国产69精品久久久久777片| 日韩不卡一区二区三区视频在线| 久热久热在线精品观看| 男女边吃奶边做爰视频| 国产探花极品一区二区| 国产 精品1| 日本av免费视频播放| 在线观看一区二区三区激情| 少妇的逼好多水| 欧美最新免费一区二区三区| 最近中文字幕高清免费大全6| 精品亚洲成a人片在线观看| 十分钟在线观看高清视频www| 午夜福利,免费看| 观看美女的网站| 99re6热这里在线精品视频| 日韩欧美精品免费久久| 久久久久久久久久久久大奶| 少妇 在线观看| av片东京热男人的天堂| 亚洲人成网站在线观看播放| 中文字幕精品免费在线观看视频 | 黄色毛片三级朝国网站| 日韩av在线免费看完整版不卡| 99热这里只有是精品在线观看| 日本色播在线视频| 日韩视频在线欧美| 久久午夜综合久久蜜桃| 免费女性裸体啪啪无遮挡网站| 亚洲av免费高清在线观看| 22中文网久久字幕| 国产一级毛片在线| 日本vs欧美在线观看视频| 国产伦理片在线播放av一区| 国产精品久久久久久av不卡| 免费女性裸体啪啪无遮挡网站| 亚洲欧美一区二区三区国产| 欧美亚洲日本最大视频资源| 久久精品国产亚洲av涩爱| 亚洲欧洲精品一区二区精品久久久 | 又黄又粗又硬又大视频| 一区二区三区精品91|