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

    不同茬口和施氮水平對(duì)南瓜根際細(xì)菌碳分解潛力的影響①

    2021-03-23 07:49:22邱虎森楊慎驕周新國(guó)
    土壤 2021年1期
    關(guān)鍵詞:茬口根際南瓜

    邱虎森,楊慎驕*,周新國(guó)

    不同茬口和施氮水平對(duì)南瓜根際細(xì)菌碳分解潛力的影響①

    邱虎森1,2,楊慎驕1,2*,周新國(guó)1,2

    (1 中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)田灌溉研究所,河南新鄉(xiāng) 453002;2 河南商丘農(nóng)田生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站,河南商丘 476000)

    為探明根際細(xì)菌群落多樣性及有機(jī)碳分解功能基因?qū)Ω鞔胧┑捻憫?yīng)特征,以收獲小麥和西瓜后種植的南瓜根際土壤為研究對(duì)象,以兩種施氮水平(純氮施入量分別為1.0、2.0 g/kg)為處理,結(jié)合高通量測(cè)序和PICRUSt功能預(yù)測(cè)解析不同茬口和施氮水平下南瓜根際細(xì)菌群落結(jié)構(gòu)和參與有機(jī)碳分解代謝細(xì)菌的關(guān)鍵功能基因豐度的差異。結(jié)果表明:各處理南瓜根際土壤細(xì)菌Shannon多樣性指數(shù)和細(xì)菌群落結(jié)構(gòu)無差異;所有處理中變形菌門、擬桿菌門、放線菌門、藍(lán)藻細(xì)菌門相對(duì)豐度之和達(dá)到90%。冗余分析發(fā)現(xiàn),小麥茬土壤中硝態(tài)氮的含量、西瓜茬土壤中銨態(tài)氮和硝態(tài)氮的含量對(duì)南瓜根際細(xì)菌群落結(jié)構(gòu)的作用高于其他理化指標(biāo)。高氮水平下,與西瓜茬相比,小麥茬南瓜根際土壤變形菌門的相對(duì)豐度降低,但擬桿菌門和放線菌門的相對(duì)豐度升高(<0.05)。無論施氮水平高低,與西瓜茬相比,小麥茬南瓜根際土壤分解半纖維素和淀粉類物質(zhì)諸如α-淀粉酶和β-半乳糖苷酶等產(chǎn)酶基因豐度提高。因此,施氮水平并未影響不同茬口土壤南瓜根際土壤微環(huán)境和細(xì)菌種群結(jié)構(gòu)和多樣性的穩(wěn)定性;但高氮處理下,茬口類型影響參與有機(jī)碳分解的關(guān)鍵細(xì)菌門類相對(duì)豐度。與西瓜茬土壤相比,小麥茬土壤可能更能促進(jìn)南瓜根際土壤有機(jī)碳的分解。

    根際土壤;細(xì)菌群落多樣性;有機(jī)碳分解代謝;施氮水平

    根際是植物與土壤物質(zhì)和能量交換的活躍區(qū)域[1],植物通過根系分泌物為根際微生物的生長(zhǎng)提供基質(zhì)和信號(hào)[2-4],影響著微生物生長(zhǎng)的生態(tài)位、細(xì)菌的形態(tài)和生理等特征[5]。而受根系分泌物形態(tài)的影響,根際微生物則能通過改變根際養(yǎng)分的有效性,影響植物對(duì)養(yǎng)分的吸收,協(xié)助抵御病原菌侵染和改變對(duì)土壤有機(jī)碳的分解速率[6-9]。土壤微生物是有機(jī)碳分解轉(zhuǎn)化的驅(qū)動(dòng)者,其中,細(xì)菌作為根際最豐富的微生物,它們對(duì)土壤有機(jī)碳分解轉(zhuǎn)化的影響不可忽視[10-11]。

    氮作為植物生長(zhǎng)所需的大量元素,氮素水平成為土壤生產(chǎn)力的主要限制性因素。綜合前人研究,施氮量對(duì)土壤細(xì)菌群落的組成和多樣性影響各異。Li等[12]研究發(fā)現(xiàn),種植玉米的砂質(zhì)黏壤土細(xì)菌群落的多樣性和豐富度隨著氮肥用量的增加而減少;而也有研究發(fā)現(xiàn)施氮量對(duì)常年耕種的農(nóng)田作物根際土壤細(xì)菌組成有顯著影響,但對(duì)細(xì)菌多樣性沒有影響[13-14]。究其原因可能歸為兩個(gè)方面:①土壤類型或者同種土壤的不同種植制度都會(huì)引起土壤理化性質(zhì)的差異和土壤細(xì)菌群落結(jié)構(gòu)的改變,進(jìn)而造成細(xì)菌種群對(duì)氮肥的響應(yīng)各異[15];②土壤細(xì)菌群落結(jié)構(gòu)的差異致其驅(qū)動(dòng)土壤氮循環(huán)過程的不同,進(jìn)而影響土壤中氮的形態(tài),而氮形態(tài)可以通過影響植物根系沉積物形態(tài)來調(diào)節(jié)根際土壤細(xì)菌的活性和群落組成[16]。基于反饋調(diào)節(jié)機(jī)制,由根系分泌物引起的根際土壤細(xì)菌群落結(jié)構(gòu)和多樣性的差異也會(huì)直接影響土壤有機(jī)碳的分解代謝[17]。大量研究表明,細(xì)菌對(duì)有機(jī)碳形態(tài)具有選擇利用性[18-19],相對(duì)于含–CH3的物質(zhì),含–COOH的有機(jī)物質(zhì)更容易被細(xì)菌礦化;革蘭氏陰性細(xì)菌菌群優(yōu)先選擇易分解態(tài)有機(jī)質(zhì),而革蘭氏陽性細(xì)菌菌群對(duì)難分解態(tài)有機(jī)質(zhì)的利用更有優(yōu)勢(shì)[8,20]。由于參與有機(jī)碳分解代謝的細(xì)菌種群復(fù)雜多樣,一種細(xì)菌可能含有參與有機(jī)碳分解的多種關(guān)鍵酶基因,單從細(xì)菌種群結(jié)構(gòu)的水平解釋土壤有機(jī)碳的分解過程遠(yuǎn)遠(yuǎn)不夠[21]。而基于PICRUSt功能預(yù)測(cè)有助于深入挖掘參與有機(jī)質(zhì)分解過程中的關(guān)鍵功能基因豐度的差異。綜上,土壤性質(zhì)和施氮水平調(diào)控細(xì)菌群落結(jié)構(gòu)和參與有機(jī)質(zhì)分解代謝的微生物過程,因此,深入了解不同茬口土壤和施氮水平對(duì)下季作物根際土壤細(xì)菌參與有機(jī)碳分解代謝過程關(guān)鍵功能基因豐度的影響,對(duì)改善種植制度和提高土壤有機(jī)碳的固定能力具有指導(dǎo)性意義。

    基于上述目的,本研究以收獲小麥和西瓜后種植的南瓜根際土壤為研究對(duì)象,以施氮水平為處理,通過南瓜盆栽試驗(yàn)研究根際土壤細(xì)菌群落結(jié)構(gòu)對(duì)不同茬口土壤和施氮水平的響應(yīng),重點(diǎn)解析參與有機(jī)碳分解代謝細(xì)菌的關(guān)鍵功能基因表達(dá)特征,為深刻理解根際土壤細(xì)菌群落參與有機(jī)碳分解代謝的生理生態(tài)過程提供一定的理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)設(shè)計(jì)

    南瓜盆栽試驗(yàn)在河南商丘農(nóng)田生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站(34°35.22′ N,1115°34.52′ E)日光溫室內(nèi)開展。供試土壤采自于站內(nèi)農(nóng)田,為西瓜茬土壤和小麥茬土壤,土壤質(zhì)地為砂壤土。采用多點(diǎn)采樣法分別采集西瓜茬和小麥茬表層0 ~ 20 cm土樣,去除動(dòng)植物殘?bào)w及碎石,風(fēng)干后過5 mm篩,充分混合備用。供試土壤基本理化指標(biāo)如表1所示。

    表1 供試土壤基本理化性質(zhì)

    供試白籽南瓜(日本鐵甲)經(jīng)消毒、漂洗以及催芽后選擇健康且長(zhǎng)勢(shì)一致的幼苗移栽到溫室中裝有10.00 kg風(fēng)干土的盆缽(直徑25 cm′高20 cm)。在南瓜苗移栽定植后(每盆2棵南瓜苗),按照相當(dāng)于純氮施入量分別為1.00和2.00 g/kg施入復(fù)合肥料(N∶P2O5∶K2O = 15∶15∶15)。按照正交試驗(yàn)設(shè)計(jì)設(shè)置4個(gè)處理,分別為:西瓜茬土+低氮、西瓜茬土+高氮、小麥茬土+低氮和小麥茬土+高氮,每個(gè)處理10個(gè)重復(fù)。在南瓜生長(zhǎng)期間保持溫室白天溫度范圍為18 ~ 30 ℃,夜間溫度低于15 ℃。

    1.2 樣品采集與測(cè)定

    1.2.1 樣品采集 在種植30 d后,先將南瓜植株從土壤中輕輕拔出,抖掉根系上的松散土,然后將與根系緊密結(jié)合的土壤刷下作為根際土。根際土混合均勻后,經(jīng)過液氮速凍處理,保存在–80 ℃冰箱,用于細(xì)菌群落多樣性分析。收集盆缽中剩余的土壤,混合均勻后用于土壤理化指標(biāo)的測(cè)定。

    1.2.2 樣品測(cè)定 采用Fast DNA?SPIN試劑盒(Qbiogene Inc., USA)提取土壤總DNA,DNA濃度和純度利用NanoDrop2000進(jìn)行檢測(cè),利用1% 瓊脂糖凝膠電泳檢測(cè)DNA提取質(zhì)量。在設(shè)計(jì)測(cè)序引物時(shí),在序列中添加bar-code序列以區(qū)分各個(gè)樣品測(cè)序數(shù)據(jù)。用338F(5′-ACTCCTACGGGAGGCAGCAG-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)引物對(duì)V3 ~ V4可變區(qū)進(jìn)行PCR擴(kuò)增,擴(kuò)增程序?yàn)椋?5 ℃預(yù)變性3 min,27個(gè)循環(huán)(95 ℃變性30 s,55 ℃退火30 s, 72 ℃延伸30 s),最后72 ℃延伸10 min(PCR儀:ABI GeneAmp? 9700型)。擴(kuò)增體系為20 μl,4 μl 5×FastPfu緩沖液,2 μl 2.5 mmol/L dNTPs,0.8 μl引物(5 μmmol/L),0.4 μl FastPfu聚合酶;10 ng DNA模板。使用2% 瓊脂糖凝膠回收PCR產(chǎn)物,利用AxyPrep DNA Gel Extraction Kit(Axygen Biosciences, Union City, CA, USA)進(jìn)行純化,Tris-HCl洗脫,2% 瓊脂糖電泳檢測(cè)。利用QuantiFluor?-ST(Promega, USA)進(jìn)行檢測(cè)定量。根據(jù)定量結(jié)果和測(cè)序量要求,取PCR產(chǎn)物構(gòu)建測(cè)序文庫。構(gòu)建好的文庫在Illumina Miseq PE300測(cè)序平臺(tái)測(cè)序。

    基于《土壤農(nóng)化分析》方法[22],稱取新鮮土樣,按照水土比10∶1(/),采用0.01 mol/L KCl浸提,AA3連續(xù)流動(dòng)分析儀測(cè)定銨態(tài)氮和硝態(tài)氮含量;有效磷采用NaHCO3浸提–鉬銻抗比色法測(cè)定;速效鉀采用NH4OAc浸提–火焰光度法測(cè)定;土壤水分含量采用105 ℃ 烘干稱重法測(cè)定。

    1.3 數(shù)據(jù)分析

    數(shù)據(jù)處理在Excel 2007中完成。Shannon多樣性指數(shù)利用QIME 1.80計(jì)算。細(xì)菌群落功能基因預(yù)測(cè)基于PICRUSt分析實(shí)現(xiàn)[23],通過PICRUSt(PICRUSt軟件存儲(chǔ)了greengene ID對(duì)應(yīng)的KEGG Ortholog(KO)信息)對(duì)OTU豐度表進(jìn)行標(biāo)準(zhǔn)化,并去除16S 標(biāo)記基因在物種基因組中的拷貝(copy)數(shù)目的影響;然后通過每個(gè)OTU 對(duì)應(yīng)的greengene ID,獲得OTU對(duì)應(yīng)的(KO)信息;并計(jì)算各功能基因?qū)?yīng)的KO豐度,該分析可在http://huttenhower.sph.harvard.edu/galaxy/中在線實(shí)現(xiàn)?;赗 3.4.3的vegan包,采用非參數(shù)多元統(tǒng)計(jì)比對(duì)組間細(xì)菌群落結(jié)構(gòu)差異;基于SPSS 16.0配對(duì)T檢驗(yàn)研究組間土壤理化性質(zhì)的差異;基于STAMP的Welch T檢驗(yàn)對(duì)比組間微生物種屬差異和組間參與氮循環(huán)的功能基因相對(duì)豐度差異;采用雙因素方差分析揭示不同茬口土壤和施氮水平對(duì)南瓜根際細(xì)菌Shannon多樣性以及微生物群落結(jié)構(gòu)的影響?;贑anoco 4.5冗余分析揭示土壤理化性質(zhì)對(duì)根際細(xì)菌群落結(jié)構(gòu)的作用。所有作圖在Origin 8.5中完成。

    2 結(jié)果與分析

    2.1 土壤理化性質(zhì)的差異

    配對(duì)T檢驗(yàn)的結(jié)果發(fā)現(xiàn)(表2),不同茬口土壤和施氮水平并未影響土壤銨態(tài)氮和硝態(tài)氮的含量,且土壤中以硝態(tài)氮含量為主。但不同茬口土壤顯著影響低氮處理?xiàng)l件下有效磷的含量,有效磷在小麥茬土壤中的含量約為在西瓜茬土壤中的1.7倍(<0.05)。在西瓜茬土壤中,高氮處理土壤有效磷的含量約為低氮處理土壤的2倍(<0.05)。在小麥茬土壤中,高氮處理土壤速效鉀的含量為426.01 mg/kg,約為低氮處理土壤的1.5倍(<0.05)。

    表2 土壤理化性質(zhì)差異

    注:表中不同小寫字母代表同種茬口土壤條件下不同施氮水平之間差異在<0.05水平顯著;不同大寫字母代表同種施氮水平下不同茬口土壤之間差異在<0.05水平顯著。

    2.2 細(xì)菌種群多樣性及群落結(jié)構(gòu)的差異

    基于圖1的結(jié)果發(fā)現(xiàn),不同茬口土壤和施氮水平并未改變根際土壤細(xì)菌Shannon多樣性指數(shù),且對(duì)細(xì)菌群落結(jié)構(gòu)的影響不顯著(表3,>0.05)。在門水平上對(duì)土壤中細(xì)菌群落結(jié)構(gòu)分析發(fā)現(xiàn),所有處理細(xì)菌群落結(jié)構(gòu)中相對(duì)豐度較高的群體為變形菌門(Proteo-bacteria)、擬桿菌門(Bacterodidetes)、放線菌門(Actinobacteria)和藍(lán)藻細(xì)菌門(Cyanobacteria),四者之和約達(dá)到90%(圖2A)。與西瓜茬土壤高氮處理相比,小麥茬土壤高氮處理根際變形菌門的相對(duì)豐度降低,但擬桿菌門和放線菌門的相對(duì)豐度均升高(<0.05)。在屬水平上,相對(duì)豐度前20位的菌屬主要屬于變形菌門、擬桿菌門和綠彎菌門(Chloroflexi),其中變形菌門假單胞菌屬()和擬桿菌門黃桿菌屬()為各處理的優(yōu)勢(shì)菌屬(圖2B)。在西瓜茬土壤中,與低氮處理相比,高氮處理降低了氫噬胞菌屬(的相對(duì)豐度(<0.05)。與西瓜茬土壤低氮處理相比,小麥茬土壤低氮處理明顯提高了假單胞菌屬的相對(duì)豐度(<0.05)。與西瓜茬土壤高氮處理相比,小麥茬土壤高氮處理明顯提高了假黃單胞菌屬()、根瘤菌屬()、黃桿菌屬()、新鞘氨醇桿菌屬()和節(jié)桿菌屬()的相對(duì)豐度(圖2B,<0.05)。冗余分析的結(jié)果發(fā)現(xiàn),土壤銨態(tài)氮、硝態(tài)氮、有效磷、速效鉀以及含水量對(duì)小麥茬土壤和西瓜茬土壤根際細(xì)菌群落的解釋率分別為98.8% 和98.2%,其中小麥茬土壤中硝態(tài)氮的作用大于其他指標(biāo),而在西瓜茬土壤中銨態(tài)氮和硝態(tài)氮的作用相對(duì)更大(圖3)。

    2.3 參與碳分解代謝的關(guān)鍵功能基因豐度差異

    基于圖4和Welch’s T檢驗(yàn)的結(jié)果表明,在同種茬口土壤條件下,施氮水平對(duì)參與有機(jī)碳分解代謝的關(guān)鍵功能基因豐度的影響不大;但在相同施肥水平下,參與有機(jī)碳分解代謝的關(guān)鍵功能基因豐度受不同茬口土壤的影響顯著。在低氮處理下,小麥茬土壤顯著提高了根際含有α-淀粉酶基因(K01176)和β-半乳糖苷酶基因的豐度(<0.05)。在高氮處理下,小麥茬土壤顯著提高了根際含有β-半乳糖苷酶基因的豐度(<0.05)。

    3 討論

    3.1 土壤理化性質(zhì)對(duì)根際土壤細(xì)菌群落結(jié)構(gòu)及多樣性的影響

    本研究結(jié)果發(fā)現(xiàn),不同茬口和施氮水平并未改變根際土壤細(xì)菌Shannon多樣性指數(shù)和細(xì)菌群落結(jié)構(gòu),說明南瓜經(jīng)過為期1個(gè)月的生長(zhǎng),土壤根際細(xì)菌群落組成已經(jīng)相當(dāng)穩(wěn)定。隨著培養(yǎng)時(shí)間的延長(zhǎng),不同處理間影響根際細(xì)菌群落結(jié)構(gòu)及多樣性的關(guān)鍵因子差異逐漸消失,根際微環(huán)境的穩(wěn)定,導(dǎo)致不同處理根際細(xì)菌所處的根際土壤微環(huán)境較為一致,并未引起不同處理間細(xì)菌群落多樣性的改變[24-25]。冗余分析的結(jié)果發(fā)現(xiàn),土壤銨態(tài)氮、硝態(tài)氮、有效磷、速效鉀以及含水量對(duì)小麥茬土壤和西瓜茬土壤根際細(xì)菌群落的解釋率分別為98.8% 和98.2%,其中小麥茬土壤中硝態(tài)氮的作用大于其他指標(biāo),而在西瓜茬土壤中銨態(tài)氮和硝態(tài)氮的作用相對(duì)更大,說明可利用氮源含量的高低在調(diào)控土壤根際細(xì)菌種群結(jié)構(gòu)中起著至關(guān)重要的作用[12]。而各處理間南瓜根際土壤銨態(tài)氮與硝態(tài)氮含量差異不顯著,這也是導(dǎo)致不同處理間細(xì)菌群落結(jié)構(gòu)和多樣性較穩(wěn)定的原因之一。基于表2的結(jié)果,小麥茬與西瓜茬南瓜根際土壤有效磷的含量存在差異,但這種差異不足以改變細(xì)菌群落結(jié)構(gòu)組成,也側(cè)面說明土壤有效磷不是調(diào)控細(xì)菌群落結(jié)構(gòu)的限制因子。

    3.2 根際土壤參與有機(jī)碳分解代謝細(xì)菌功能菌群的變化

    在門水平上對(duì)土壤中細(xì)菌群落結(jié)構(gòu)分析發(fā)現(xiàn),所有處理細(xì)菌群落結(jié)構(gòu)中相對(duì)豐度較高的群體為變形菌門(Proteobacteria)、擬桿菌門(Bacterodidetes)和放線菌門(Actinobacteria),三者之和約達(dá)到90%,說明這3大門類的細(xì)菌種群具有較強(qiáng)的環(huán)境適應(yīng)能力,無論是土壤環(huán)境的改變還是人為活動(dòng)的影響,都不會(huì)影響細(xì)菌三大門類的主導(dǎo)地位[13,26],但環(huán)境的差異會(huì)影響到其相對(duì)豐度的變化。大量研究認(rèn)為,土壤可利用碳源含量及形態(tài)的改變都會(huì)影響參與碳分解的微生物種類和生物量[8,18]。變形菌門在易分解有機(jī)碳組分的分解中占主導(dǎo)地位,擬桿菌門種群可以表達(dá)糖基水解酶參與纖維素和半纖維素分解的程度[27-28],而放線菌門屬于主要參與難分解有機(jī)碳降解的門類[29]。在本研究中,高氮處理下,小麥茬南瓜根際土壤變形菌門的相對(duì)豐度低于西瓜茬土壤,但擬桿菌門和放線菌門的相對(duì)豐度高于西瓜茬土壤。這從側(cè)面說明在高氮處理下,小麥茬南瓜根際土壤中易利用碳源的含量可能低于西瓜茬南瓜根際土壤[30],但參與難分解有機(jī)碳組分微生物種群的增加可能會(huì)加速小麥茬南瓜根際碳后期的分解速率。

    基于圖4和Welch’s T檢驗(yàn)的結(jié)果,施氮水平對(duì)同種茬口南瓜根際土壤參與有機(jī)碳分解的關(guān)鍵功能基因豐度無顯著影響,說明氮素含量雖然在細(xì)菌群落結(jié)構(gòu)中的作用較大,但并不是限制有機(jī)碳分解的因子。但無論是在低氮處理還是在高氮處理,小麥茬土壤β-半乳糖苷酶基因(K01190)的豐度均顯著高于西瓜茬土壤;另外,在低氮處理下,與西瓜茬南瓜根際土壤相比,小麥茬南瓜根際土壤顯著提高了含有α-淀粉酶基因(K01176)豐度。研究表明β-半乳糖苷酶基因是編碼β-半乳糖苷酶的基因,β-半乳糖苷酶是分解半纖維素的關(guān)鍵酶類,而α-淀粉酶基因編碼形成α-淀粉酶,可以水解淀粉內(nèi)部的α-1,4-糖苷鍵,水解產(chǎn)物為糊精、低聚糖和單糖[30]??傊@兩種基因都是促進(jìn)復(fù)雜有機(jī)物質(zhì)分解轉(zhuǎn)化成簡(jiǎn)單有機(jī)物質(zhì)的關(guān)鍵,因而小麥茬土壤根際細(xì)菌可能會(huì)加速有機(jī)質(zhì)分解過程[26]。

    4 結(jié)論

    1)不同茬口土壤和施氮水平并未改變南瓜根際土壤細(xì)菌Shannon多樣性指數(shù)和群落結(jié)構(gòu),可利用氮源含量在調(diào)節(jié)南瓜根際土壤細(xì)菌群落結(jié)構(gòu)中的作用至關(guān)重要。

    2)高氮處理下,小麥茬土壤根際細(xì)菌群落中參與難分解有機(jī)物質(zhì)的種群豐度提升;且無論施氮高低,與西瓜茬土壤相比,小麥茬土壤的南瓜根際參與有機(jī)質(zhì)分解的關(guān)鍵基因豐度顯著提升,因此,與小麥茬土壤相比,西瓜茬土壤更有利于有機(jī)碳的固定。

    [1] Yang Z P, Yang W P, Li S C, et al. Variation of bacterial community diversity in rhizosphere soil of sole-cropped versus intercropped wheat field after harvest[J]. PLoS One, 2016, 11(3): e0150618.

    [2] Chaparro J M, Badri D V, Vivanco J M. Rhizosphere microbiome assemblage is affected by plant develop-ment [J]. The ISME Journal, 2014, 8(4): 790–803.

    [3] Huang X F, Chaparro J M, Reardon K F, et al. Rhizosphere interactions: root exudates, microbes, and microbial communities[J]. Botany, 2014, 92(4): 267–275.

    [4] 王曉婷, 陳瑞蕊, 井忠旺, 等. 水稻和小麥根際效應(yīng)及細(xì)菌群落特征的比較研究[J]. 土壤學(xué)報(bào), 2019, 56(2): 443–453.

    [5] Ibekwe A M, Poss J A, Grattan S R, et al. Bacterial diversity in cucumber () rhizosphere in response to salinity, soil pH, and boron[J]. Soil Biology and Biochemistry, 2010, 42(4):567–575.

    [6] Morrissey J P, Dow J M, Mark G L, et al. Are microbes at the root of a solution to world food production?[J]. EMBO Reports, 2004, 5(10): 922–926.

    [7] Mendes R, Kruijt M, de Bruijn I, et al. Deciphering the rhizosphere microbiome for disease-suppressive bacteria[J]. Science, 2011, 332(6033): 1097–1100.

    [8] Bird J A, Herman D J, Firestone M K. Rhizosphere priming of soil organic matter by bacterial groups in a grassland soil[J]. Soil Biology and Biochemistry, 2011, 43(4): 718–725.

    [9] 楊珍, 戴傳超, 王興祥, 等. 作物土傳真菌病害發(fā)生的根際微生物機(jī)制研究進(jìn)展[J]. 土壤學(xué)報(bào), 2019, 56(1): 12–22.

    [10] Six J, Frey S D, Thiet R K, et al. Bacterial and fungal contributions to carbon sequestration in agroecosystems[J]. Soil Science Society of America Journal, 2006, 70(2): 555–569.

    [11] Prashar P, Kapoor N, Sachdeva S. Rhizosphere: its structure, bacterial diversity and significance[J]. Reviews in Environmental Science and Bio, 2014, 13(1): 63–77.

    [12] Li J G, Shen M C, Hou J F, et al. Effect of different levels of nitrogen on rhizosphere bacterial community structure in intensive monoculture of greenhouse lettuce[J]. Scientific Reports, 2016, 6: 25305.

    [13] Fierer N, Bradford M A, Jackson R B. Toward an ecological classification of soil bacteria[J]. Ecology, 2007, 88(6): 1354–1364.

    [14] Wang R Q, Xiao Y P, Lv F, et al. Bacterial community structure and functional potential of rhizosphere soils as influenced by nitrogen addition and bacterial wilt disease under continuous sesame cropping[J]. Applied Soil Ecology, 2018, 125:117–127.

    [15] Souza R C, Hungria M, Cant?o M E, et al. Metagenomic analysis reveals microbial functional redundancies and specificities in a soil under different tillage and crop-management regimes[J]. Applied Soil Ecology, 2015, 86:106–112.

    [16] Giagnoni L, Pastorelli R, Mocali S, et al. Availability of different nitrogen forms changes the microbial communities and enzyme activities in the rhizosphere of maize lines with different nitrogen use efficiency[J]. Applied Soil Ecology, 2016, 98:30–38.

    [17] Wang X J, Tang C X, Severi J, et al. Rhizosphere priming effect on soil organic carbon decomposition under plant species differing in soil acidification and root exudation[J]. New Phytologist, 2016, 211(3): 864–873.

    [18] Qiu H S, Zheng X D, Ge T, et al. Weaker priming and mineralisation of low molecular weight organic substances in paddy than in upland soil[J]. European Journal of Soil Biology, 2017, 83:9–17.

    [19] 邱虎森, 蘇以榮, 劉杰云, 等. 易利用態(tài)有機(jī)物質(zhì)對(duì)水稻土甲烷排放的激發(fā)作用[J]. 2018, 50(3): 537–542.

    [20] Derrien D, Plain C, Courty P E, et al. Does the addition of labile substrate destabilise old soil organic matter?[J]. Soil Biology and Biochemistry, 2014, 76:149–160.

    [21] 趙帆, 趙密珍, 王鈺, 等. 龐夫花基于高通量測(cè)序研究草莓根際微生物群落結(jié)構(gòu)和多樣性[J]. 2019, 51(1): 51–60

    [22] 鮑士旦. 土壤農(nóng)化分析[M]. 3版. 北京: 中國(guó)農(nóng)業(yè)出版社, 2000.

    [23] Langille M G I, Zaneveld J, Caporaso J G, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences[J]. Nature Biotechnology, 2013, 31(9): 814–821.

    [24] Zhalnina K, Louie K B, Hao Z, et al. Dynamic root exudate chemistry and microbial substrate preferences drive patterns in rhizosphere microbial community assembly[J]. Nature Microbiology, 2018, 3(4): 470–480.

    [25] Brimecombe M J, De Leij F A A M, Lynch J M. Rhizodeposition and microbial populations//Pinton R, Veranini Z, Nannipieri P. The rhizosphere biochemistry and organic substances at the soil–plant interface. New York, USA: Taylor & Francis Group, 2007: 73–98.

    [26] Zheng W, Zhao Z Y, Gong Q L, et al. Effects of cover crop in an apple orchard on microbial community composition, networks, and potential genes involved with degradation of crop residues in soil[J]. Biology and Fertility of Soils, 2018, 54(6): 743–759.

    [27] Wegner C E, Liesack W. Microbial community dynamics during the early stages of plant polymer breakdown in paddy soil[J]. Environmental Microbiology, 2016, 18(9): 2825–2842.

    [28] Wang G Z, Wang Y R, Yang P L, et al. Molecular detection and diversity of xylanase genes in alpine tundra soil[J]. Applied Microbiology and Biotechnology, 2010, 87(4): 1383–1393.

    [29] Sonia M T, Hafedh B, Abdennaceur H, et al. Studies on the ecology of actinomycetes in an agricultural soil amended with organic residues: II. Assessment of enzymatic activities of Actinomycetales isolates[J]. World Journal of Microbiology and Biotechnology, 2011, 27(10): 2251– 2259.

    [30] Trivedi P, Anderson I C, Singh B K. Microbial modulators of soil carbon storage: integrating genomic and metabolic knowledge for global prediction[J]. Trends in Microbiology, 2013, 21(12): 641–651.

    Effects of Different Crop Rotation and Nitrogen Fertilization on Carbon Decomposition Potential of Pumpkin Rhizosphere Bacteria

    QIU Husen1,2, YANG Shenjiao1,2*, ZHOU Xinguo1,2

    (1 Farmland Irrigation Research Institute, Chinese Academy of Agricultural Sciences, Xinxiang, Henan 453002, China; 2 National Agro-ecological System Observation and Research Station of Shangqiu, Shangqiu, Henan 476000, China)

    In order to understand better the responses of rhizosphere bacterial diversity and organic carbon decomposition metabolism genes to tillage measures, rhizosphere soils of pumpkin rotated with wheat or watermelon were taken as test objects, and two nitrogen levels (N, 1.0 and 2.0 g/kg) were designed, the rhizosphere bacterial community structure and abundance of organic carbon decomposition metabolism genes were analyzed based on high-throughput sequencing technologies and PICRUSt analysis. The results showed that no significant difference was found in bacterial Shannon diversity index and community structure between different treatments. Proteobacteria, Bacterodidetes, Actinobacteria and Cyanobacteria were four dominant phyla in rhizosphere bacterial community structures under all treatments, with the sum of relative abundance reached 90%. Redundancy analysis showed that nitrate nitrogen in soils rotated with wheat, ammonium nitrogen and nitrate nitrogen in soils rotated with watermelon, played more important role in regulating bacterial community structure than those of other physiochemical indexes. At high nitrogen level, the relative abundance of Proteobacteriain soils rotated with wheat was lower than that in soils rotated with watermelon, but the relative abundance of Bacterodidetes and Actinobacteria was higher than that in soils rotated with watermelon (<0.05). The abundances of enzymes gene expression such as α-amylase and β-galactosidase, important for the decomposition of hemicellulose and starch, were higher in soils rotated with wheat compared to that in soils rotated with watermelon. In conclusion, crop rotation and nitrogen application influence little on the microenvironment of pumpkin rhizosphere soil and the stability of bacterial community structure and diversity. At high nitrogen level, soil types rotated with different crops affect the relative abundance of bacterial phylum participating in organic carbon decomposition. The decomposition rate of organic carbon in soils rotated with wheat is higher than that in soils rotated with watermelon.

    Rhizosphere soil; Bacterial diversity; Organic carbon decomposition metabolism; Nitrogen fertilization

    S154.1

    A

    10.13758/j.cnki.tr.2021.01.018

    邱虎森, 楊慎驕, 周新國(guó). 不同茬口和施氮水平對(duì)南瓜根際細(xì)菌碳分解潛力的影響. 土壤, 2021, 53(1): 133–139.

    國(guó)家自然科學(xué)基金項(xiàng)目(51609246)和河南省重點(diǎn)研發(fā)與推廣專項(xiàng)(192102110087)資助。

    (yangshenjiao@hotmail.com)

    邱虎森(1987—),男,河南項(xiàng)城人,博士,助理研究員,主要從事土壤生物與元素循環(huán)研究。E-mail: qiuhusen2008@163.com

    猜你喜歡
    茬口根際南瓜
    根際微生物對(duì)植物與土壤交互調(diào)控的研究進(jìn)展
    蕎麥輪作倒茬對(duì)比試驗(yàn)
    茬口的熊
    玉米‖花生茬口對(duì)冬小麥旗葉光化學(xué)活性的影響
    南瓜燈
    南瓜
    大灰狼(2018年10期)2018-10-31 14:47:16
    黃花蒿葉水提物對(duì)三七根際尖孢鐮刀菌生長(zhǎng)的抑制作用
    促植物生長(zhǎng)根際細(xì)菌HG28-5對(duì)黃瓜苗期生長(zhǎng)及根際土壤微生態(tài)的影響
    南瓜車開動(dòng)啦
    鹽脅迫下白蠟根際微域的營(yíng)養(yǎng)元素狀況
    伦理电影大哥的女人| 国产精品国产三级国产av玫瑰| 少妇熟女欧美另类| 草草在线视频免费看| 欧美日韩视频精品一区| 男女边摸边吃奶| 欧美日韩综合久久久久久| 看免费成人av毛片| 91精品一卡2卡3卡4卡| 久久狼人影院| 日本欧美视频一区| 亚洲国产毛片av蜜桃av| 精品人妻一区二区三区麻豆| 又爽又黄a免费视频| 亚洲三级黄色毛片| av视频免费观看在线观看| 国产成人aa在线观看| 日韩欧美精品免费久久| 午夜91福利影院| 老司机影院成人| 亚洲综合精品二区| 成人漫画全彩无遮挡| 偷拍熟女少妇极品色| 黑人猛操日本美女一级片| 国产亚洲欧美精品永久| 各种免费的搞黄视频| 视频区图区小说| 亚洲av男天堂| 亚洲一级一片aⅴ在线观看| 久久久久久久久大av| 日韩不卡一区二区三区视频在线| 久久99蜜桃精品久久| 精品久久久久久久久亚洲| 日日撸夜夜添| 少妇猛男粗大的猛烈进出视频| 亚洲精品中文字幕在线视频 | 十八禁网站网址无遮挡 | 高清毛片免费看| 内射极品少妇av片p| 亚洲精品自拍成人| 久久综合国产亚洲精品| 极品教师在线视频| 亚洲欧美日韩东京热| 国产午夜精品久久久久久一区二区三区| 免费av中文字幕在线| 一本一本综合久久| 欧美97在线视频| 久久久久久久久久人人人人人人| 菩萨蛮人人尽说江南好唐韦庄| 久久 成人 亚洲| 久热久热在线精品观看| 成人亚洲精品一区在线观看| 我的女老师完整版在线观看| 成年人午夜在线观看视频| 国产高清三级在线| 亚洲av欧美aⅴ国产| 大片电影免费在线观看免费| 老司机亚洲免费影院| 丝袜喷水一区| 交换朋友夫妻互换小说| 少妇被粗大猛烈的视频| 99久久中文字幕三级久久日本| 777米奇影视久久| 日韩三级伦理在线观看| 亚洲综合色惰| 日本猛色少妇xxxxx猛交久久| 色婷婷av一区二区三区视频| 精品久久久久久电影网| av免费在线看不卡| 日韩伦理黄色片| 日韩大片免费观看网站| 国产男女内射视频| 日韩强制内射视频| 看非洲黑人一级黄片| 黑丝袜美女国产一区| 成人亚洲欧美一区二区av| 人体艺术视频欧美日本| 亚洲精品国产成人久久av| 嫩草影院入口| 色视频www国产| √禁漫天堂资源中文www| 高清欧美精品videossex| 在线观看人妻少妇| 天天操日日干夜夜撸| 高清av免费在线| 国产欧美另类精品又又久久亚洲欧美| 水蜜桃什么品种好| 亚洲久久久国产精品| 如日韩欧美国产精品一区二区三区 | 久久这里有精品视频免费| 涩涩av久久男人的天堂| 你懂的网址亚洲精品在线观看| 精品一区二区免费观看| 春色校园在线视频观看| 亚洲三级黄色毛片| 最新的欧美精品一区二区| 久久影院123| 青春草国产在线视频| 熟女av电影| av播播在线观看一区| 免费大片18禁| 乱码一卡2卡4卡精品| 夜夜看夜夜爽夜夜摸| av又黄又爽大尺度在线免费看| 女性被躁到高潮视频| 日韩电影二区| 免费人妻精品一区二区三区视频| 国产精品久久久久久av不卡| av在线app专区| 欧美日韩一区二区视频在线观看视频在线| 又粗又硬又长又爽又黄的视频| 亚洲国产精品成人久久小说| 不卡视频在线观看欧美| 最黄视频免费看| 少妇猛男粗大的猛烈进出视频| 新久久久久国产一级毛片| 一级毛片黄色毛片免费观看视频| av又黄又爽大尺度在线免费看| 人妻少妇偷人精品九色| 亚洲欧美日韩东京热| 亚洲国产av新网站| 国产片特级美女逼逼视频| 精品视频人人做人人爽| 黄片无遮挡物在线观看| 91久久精品国产一区二区成人| 久久精品国产亚洲av天美| 老司机影院毛片| 一本色道久久久久久精品综合| 色婷婷久久久亚洲欧美| 亚洲在久久综合| 久久精品国产自在天天线| 亚洲成人手机| 欧美精品高潮呻吟av久久| 精品人妻偷拍中文字幕| 国产av国产精品国产| 丝瓜视频免费看黄片| 天堂中文最新版在线下载| 国产男女内射视频| 国产精品久久久久久精品古装| 国产在线免费精品| 少妇人妻久久综合中文| 99热全是精品| 精品国产一区二区三区久久久樱花| 日日撸夜夜添| www.色视频.com| 欧美一级a爱片免费观看看| 精品亚洲成a人片在线观看| 赤兔流量卡办理| .国产精品久久| 国产黄色免费在线视频| 国产亚洲5aaaaa淫片| 国产伦精品一区二区三区视频9| xxx大片免费视频| 久久午夜综合久久蜜桃| 日本爱情动作片www.在线观看| 国产精品熟女久久久久浪| 国产成人aa在线观看| 色视频www国产| 青春草视频在线免费观看| 亚洲欧美成人精品一区二区| 亚洲第一区二区三区不卡| 精品国产一区二区久久| 啦啦啦视频在线资源免费观看| av在线app专区| a级毛片在线看网站| 在线观看美女被高潮喷水网站| 久久久久久久久久久丰满| 欧美精品一区二区大全| 亚洲精品久久午夜乱码| 新久久久久国产一级毛片| 精品人妻偷拍中文字幕| 大片免费播放器 马上看| 免费少妇av软件| 最新的欧美精品一区二区| 免费观看无遮挡的男女| 国产精品国产av在线观看| 哪个播放器可以免费观看大片| videossex国产| 午夜免费观看性视频| 国产精品无大码| 99九九在线精品视频 | 欧美精品国产亚洲| 欧美日韩亚洲高清精品| 美女大奶头黄色视频| 视频中文字幕在线观看| 精品视频人人做人人爽| 精品一区二区免费观看| 久久人人爽人人片av| 80岁老熟妇乱子伦牲交| 日日撸夜夜添| 欧美成人精品欧美一级黄| 美女脱内裤让男人舔精品视频| 男女边吃奶边做爰视频| 自线自在国产av| 久久久久久久久久久久大奶| 久久午夜福利片| 秋霞伦理黄片| 美女国产视频在线观看| 熟女人妻精品中文字幕| 国产白丝娇喘喷水9色精品| 欧美精品国产亚洲| 国产在线一区二区三区精| 精品人妻偷拍中文字幕| 夜夜骑夜夜射夜夜干| 赤兔流量卡办理| 99久久精品热视频| 国产黄色视频一区二区在线观看| 老司机影院毛片| 狂野欧美激情性xxxx在线观看| 精品亚洲乱码少妇综合久久| 成人亚洲精品一区在线观看| 日韩熟女老妇一区二区性免费视频| 久久久a久久爽久久v久久| 国产伦理片在线播放av一区| 黄色一级大片看看| 欧美精品人与动牲交sv欧美| 看非洲黑人一级黄片| 久久影院123| 国产成人a∨麻豆精品| 欧美xxxx性猛交bbbb| 欧美精品高潮呻吟av久久| 99国产精品免费福利视频| 纵有疾风起免费观看全集完整版| 少妇丰满av| 日韩强制内射视频| 麻豆精品久久久久久蜜桃| 亚洲精品一二三| 九九久久精品国产亚洲av麻豆| 国产片特级美女逼逼视频| 岛国毛片在线播放| 欧美性感艳星| 在线观看免费视频网站a站| 亚洲欧洲国产日韩| 日本黄色片子视频| 欧美日韩一区二区视频在线观看视频在线| 黄色怎么调成土黄色| 在线观看一区二区三区激情| 国内揄拍国产精品人妻在线| 免费黄色在线免费观看| 国产成人精品久久久久久| 欧美日韩视频高清一区二区三区二| 免费黄频网站在线观看国产| 欧美日韩精品成人综合77777| 国产精品嫩草影院av在线观看| 国产探花极品一区二区| 国产精品欧美亚洲77777| 九九爱精品视频在线观看| 男女啪啪激烈高潮av片| 99久国产av精品国产电影| 久久精品国产亚洲av涩爱| 久久人人爽人人爽人人片va| av在线app专区| 少妇人妻一区二区三区视频| 久久精品国产亚洲网站| 在线观看av片永久免费下载| 99久久精品国产国产毛片| 看免费成人av毛片| 人人妻人人添人人爽欧美一区卜| 欧美变态另类bdsm刘玥| 777米奇影视久久| 狠狠精品人妻久久久久久综合| 亚洲激情五月婷婷啪啪| 亚洲欧美中文字幕日韩二区| 国产av码专区亚洲av| av国产精品久久久久影院| av视频免费观看在线观看| 99久久中文字幕三级久久日本| 美女视频免费永久观看网站| 国产精品国产三级专区第一集| 免费观看的影片在线观看| 大香蕉久久网| 国产亚洲最大av| 亚洲图色成人| 亚洲欧美精品自产自拍| 交换朋友夫妻互换小说| 成人亚洲欧美一区二区av| 极品人妻少妇av视频| 男男h啪啪无遮挡| 精品亚洲成国产av| 日本-黄色视频高清免费观看| 精品久久久久久久久av| 最近2019中文字幕mv第一页| 国产在线男女| 欧美亚洲 丝袜 人妻 在线| 亚洲内射少妇av| 国产精品99久久99久久久不卡 | 97在线视频观看| 国产成人aa在线观看| 赤兔流量卡办理| 内地一区二区视频在线| 成人黄色视频免费在线看| 麻豆成人午夜福利视频| 99精国产麻豆久久婷婷| 午夜福利视频精品| 少妇高潮的动态图| 嫩草影院新地址| 午夜免费观看性视频| 国产精品国产三级国产av玫瑰| 亚洲精品国产av蜜桃| 中文欧美无线码| 国产 精品1| 国产av码专区亚洲av| 少妇猛男粗大的猛烈进出视频| 亚洲性久久影院| 国产精品三级大全| 男人添女人高潮全过程视频| 国产亚洲91精品色在线| 欧美精品高潮呻吟av久久| 99久久精品热视频| 久久久久久久久久久丰满| 欧美日韩精品成人综合77777| 青春草国产在线视频| 美女大奶头黄色视频| 大香蕉97超碰在线| 亚洲欧美一区二区三区黑人 | av有码第一页| 国产综合精华液| a级一级毛片免费在线观看| 亚洲自偷自拍三级| 成人美女网站在线观看视频| 精品午夜福利在线看| 涩涩av久久男人的天堂| 草草在线视频免费看| 女的被弄到高潮叫床怎么办| 最黄视频免费看| 亚洲成人手机| 成人免费观看视频高清| 亚洲欧美日韩卡通动漫| 精品酒店卫生间| 日本av手机在线免费观看| a级片在线免费高清观看视频| 精品午夜福利在线看| 观看免费一级毛片| 国产一区亚洲一区在线观看| 黄色欧美视频在线观看| 国产男人的电影天堂91| 男人添女人高潮全过程视频| 爱豆传媒免费全集在线观看| 精品久久久噜噜| 日日撸夜夜添| 一个人免费看片子| 日日撸夜夜添| 欧美成人精品欧美一级黄| 国产一区二区三区av在线| 秋霞在线观看毛片| 日本与韩国留学比较| 国产精品国产三级国产专区5o| 亚洲精品乱码久久久久久按摩| 亚洲精品日本国产第一区| 国产av码专区亚洲av| 亚洲国产最新在线播放| 人体艺术视频欧美日本| 插阴视频在线观看视频| 亚洲欧美精品自产自拍| 久久久久精品性色| 亚洲精品自拍成人| 韩国高清视频一区二区三区| 亚洲欧美精品专区久久| 超碰97精品在线观看| 国产午夜精品久久久久久一区二区三区| 免费在线观看成人毛片| 日日啪夜夜撸| 成人二区视频| 观看av在线不卡| 欧美成人精品欧美一级黄| 免费观看a级毛片全部| 久热久热在线精品观看| 看非洲黑人一级黄片| 中文字幕制服av| 哪个播放器可以免费观看大片| 一级av片app| 国产精品福利在线免费观看| 久热久热在线精品观看| 婷婷色综合大香蕉| 高清毛片免费看| 国产伦精品一区二区三区视频9| 高清毛片免费看| 亚洲怡红院男人天堂| 蜜桃久久精品国产亚洲av| 少妇被粗大猛烈的视频| 精品亚洲乱码少妇综合久久| 日韩三级伦理在线观看| av专区在线播放| 欧美另类一区| 亚洲精华国产精华液的使用体验| 高清av免费在线| 国产乱人偷精品视频| 欧美日韩视频精品一区| 2021少妇久久久久久久久久久| 在线播放无遮挡| 亚洲精品乱久久久久久| 在线亚洲精品国产二区图片欧美 | 91在线精品国自产拍蜜月| 内射极品少妇av片p| 日本爱情动作片www.在线观看| 人妻制服诱惑在线中文字幕| 一二三四中文在线观看免费高清| 国产在线免费精品| 99久久综合免费| 国产精品一区二区在线观看99| 成人无遮挡网站| 国产熟女欧美一区二区| 国产精品嫩草影院av在线观看| 亚洲精品aⅴ在线观看| 国产 精品1| 久久久午夜欧美精品| 麻豆乱淫一区二区| 国产一区二区三区av在线| 18禁动态无遮挡网站| 久久99蜜桃精品久久| 男人和女人高潮做爰伦理| 久久久精品免费免费高清| 国产极品天堂在线| 中文资源天堂在线| 交换朋友夫妻互换小说| 一级片'在线观看视频| 精品久久久精品久久久| 免费人成在线观看视频色| 汤姆久久久久久久影院中文字幕| 久久精品久久精品一区二区三区| 91成人精品电影| 亚洲精品久久久久久婷婷小说| 国精品久久久久久国模美| 亚洲精品日韩在线中文字幕| 免费观看无遮挡的男女| 国产伦精品一区二区三区四那| 色吧在线观看| 免费黄网站久久成人精品| 亚洲精品日韩av片在线观看| 一本大道久久a久久精品| 欧美高清成人免费视频www| 日韩制服骚丝袜av| 天堂俺去俺来也www色官网| 欧美日韩综合久久久久久| 欧美日韩视频高清一区二区三区二| 午夜免费观看性视频| 大码成人一级视频| 久久久久久人妻| 久久99蜜桃精品久久| 国产精品99久久99久久久不卡 | 另类精品久久| 日韩av在线免费看完整版不卡| 日韩免费高清中文字幕av| 黄色一级大片看看| 午夜免费鲁丝| 另类亚洲欧美激情| 欧美精品国产亚洲| 国产高清三级在线| 99热网站在线观看| 日韩视频在线欧美| 久久综合国产亚洲精品| 老司机影院毛片| 女性被躁到高潮视频| 天天躁夜夜躁狠狠久久av| 国产日韩欧美视频二区| 国产极品粉嫩免费观看在线 | 中国美白少妇内射xxxbb| 不卡视频在线观看欧美| av视频免费观看在线观看| 777米奇影视久久| 在线观看一区二区三区激情| 18禁裸乳无遮挡动漫免费视频| 免费黄网站久久成人精品| 看非洲黑人一级黄片| 丝袜在线中文字幕| 免费观看a级毛片全部| 另类精品久久| av播播在线观看一区| 麻豆乱淫一区二区| 日本-黄色视频高清免费观看| 亚洲欧洲精品一区二区精品久久久 | 午夜福利,免费看| 亚洲av电影在线观看一区二区三区| 成人午夜精彩视频在线观看| 国产精品免费大片| 又黄又爽又刺激的免费视频.| 一区二区av电影网| 三上悠亚av全集在线观看 | 国产视频内射| 欧美xxⅹ黑人| 精品卡一卡二卡四卡免费| 人人妻人人爽人人添夜夜欢视频 | 免费看日本二区| 边亲边吃奶的免费视频| 国产精品熟女久久久久浪| 熟女电影av网| 精品卡一卡二卡四卡免费| 三级经典国产精品| 免费大片18禁| 99re6热这里在线精品视频| 精品少妇久久久久久888优播| 在线精品无人区一区二区三| 我要看日韩黄色一级片| 在线播放无遮挡| 久久人人爽av亚洲精品天堂| 人妻 亚洲 视频| 丰满乱子伦码专区| 日韩一区二区视频免费看| 黑人高潮一二区| 精品亚洲成a人片在线观看| 国产伦精品一区二区三区视频9| 亚洲不卡免费看| 少妇人妻一区二区三区视频| 成人二区视频| 免费高清在线观看视频在线观看| 亚洲欧美一区二区三区国产| 国产熟女欧美一区二区| 精品一区二区三区视频在线| a级毛片在线看网站| 曰老女人黄片| 午夜福利网站1000一区二区三区| av国产精品久久久久影院| 日本91视频免费播放| 99久久精品热视频| 国产日韩欧美在线精品| 国内揄拍国产精品人妻在线| 老熟女久久久| 国内揄拍国产精品人妻在线| 免费黄网站久久成人精品| 久久久久久久久大av| 国产深夜福利视频在线观看| 午夜日本视频在线| 黄色视频在线播放观看不卡| 男女免费视频国产| 国产美女午夜福利| 边亲边吃奶的免费视频| 国产中年淑女户外野战色| 国产极品天堂在线| 韩国高清视频一区二区三区| 国产欧美日韩一区二区三区在线 | 午夜视频国产福利| a级毛片免费高清观看在线播放| 少妇的逼水好多| 国产亚洲91精品色在线| 久久久欧美国产精品| 人人妻人人添人人爽欧美一区卜| 精品人妻熟女毛片av久久网站| 最近中文字幕高清免费大全6| 久久久久久久久久久免费av| 亚洲,一卡二卡三卡| 久久狼人影院| 日韩精品有码人妻一区| 久久久久久久国产电影| 婷婷色综合大香蕉| 亚洲图色成人| 黑人高潮一二区| 在线观看免费视频网站a站| 激情五月婷婷亚洲| 伦精品一区二区三区| 91精品国产国语对白视频| 国产有黄有色有爽视频| 亚洲国产色片| 久久热精品热| 午夜91福利影院| 国产熟女欧美一区二区| 国内精品宾馆在线| 我的老师免费观看完整版| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| 国产中年淑女户外野战色| 亚洲美女黄色视频免费看| 日本黄色片子视频| 免费高清在线观看视频在线观看| 狠狠精品人妻久久久久久综合| 国产精品人妻久久久影院| 观看免费一级毛片| 99热网站在线观看| 少妇精品久久久久久久| 精华霜和精华液先用哪个| 亚洲av在线观看美女高潮| 日韩大片免费观看网站| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 欧美成人午夜免费资源| 亚洲av福利一区| 国产伦理片在线播放av一区| 久久精品久久精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 91在线精品国自产拍蜜月| 大陆偷拍与自拍| 国产淫片久久久久久久久| 久久99蜜桃精品久久| 亚洲人成网站在线播| 亚洲精品色激情综合| 国产午夜精品一二区理论片| 久久青草综合色| av福利片在线观看| 夜夜骑夜夜射夜夜干| 我的老师免费观看完整版| 99热6这里只有精品| 国产黄色免费在线视频| 中文字幕制服av| 国产亚洲91精品色在线| 人人妻人人看人人澡| 日日摸夜夜添夜夜爱| 日本-黄色视频高清免费观看| 99热这里只有是精品在线观看| 少妇人妻久久综合中文| 久热这里只有精品99| 国产黄片视频在线免费观看| 国产有黄有色有爽视频| 麻豆精品久久久久久蜜桃| 亚洲国产精品专区欧美| 国产成人午夜福利电影在线观看| 91精品伊人久久大香线蕉| 国产永久视频网站| 如何舔出高潮| 一边亲一边摸免费视频| 黄色配什么色好看| 国产精品熟女久久久久浪| 久久国产乱子免费精品| 在线观看人妻少妇| 国产精品女同一区二区软件| 你懂的网址亚洲精品在线观看| 色视频www国产|