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

    輪作模式對植煙土壤酶活性及真菌群落的影響

    2016-07-21 01:25:46陳丹梅段玉琪楊宇虹
    生態(tài)學(xué)報(bào) 2016年8期
    關(guān)鍵詞:多樣性輪作真菌

    陳丹梅, 段玉琪,楊宇虹, 晉 艷, 袁 玲,*

    1 西南大學(xué)資源環(huán)境學(xué)院,重慶 400716 2 云南省煙草農(nóng)業(yè)科學(xué)研究院,昆明 650031

    ?

    輪作模式對植煙土壤酶活性及真菌群落的影響

    陳丹梅1, 段玉琪2,楊宇虹2, 晉艷2, 袁玲1,*

    1 西南大學(xué)資源環(huán)境學(xué)院,重慶400716 2 云南省煙草農(nóng)業(yè)科學(xué)研究院,昆明650031

    摘要:種植模式顯著影響土壤理化生物學(xué)性質(zhì),并與作物土傳真菌病害的發(fā)生密切相關(guān)。試驗(yàn)選擇云南省具有代表性的紅壤,設(shè)置烤煙-休閑-玉米(T-B-M)、烤煙-油菜-玉米(T-C-M)、烤煙-油菜-水稻(T-C-R)和烤煙-苕子-水稻(T-V-R)等4種輪作模式,利用化學(xué)、酶學(xué)分析及454高通量測序技術(shù),研究了土壤養(yǎng)分、酶活性及真菌群落結(jié)構(gòu),旨在為烤煙的合理輪作提供科學(xué)依據(jù)。經(jīng)16a的不同輪作種植后,土壤pH變化于5.6—6.4之間,仍然適合種植水稻、玉米、油菜和烤煙等多種作物。在T-V-R處理的土壤中,烤煙產(chǎn)量、煙葉產(chǎn)值和上中等煙比例最高,有機(jī)質(zhì)比原初提高45.11%,堿解氮、有效磷和微生物量碳氮顯著增加,蔗糖酶、脫氫酶、脲酶、磷酸酶和過氧化氫酶活性也顯著高于其它輪作模式,說明T-V-R改善了土壤生態(tài)環(huán)境,促進(jìn)了微生物繁衍,數(shù)量增加,活性增強(qiáng)。土壤真菌的18S rDNA讀數(shù)依次為13097(T-B-M)、11345(T-C-M)、12939(T-C-R)和13763(T-V-R),分別代表530、378、395和581種(類)的真菌,由子囊菌門、擔(dān)子菌門、接合菌門、壺菌門和尚待鑒定的真菌等構(gòu)成,其中尚待鑒定的真菌屬種和子囊菌門占絕大部分。在不同輪作處理的土壤中,前15種優(yōu)勢真菌的豐富度變化于29.46%—62.86%之間,優(yōu)勢菌株的相似性極低,說明土壤真菌的種群結(jié)構(gòu)因輪作模式不同而異。T-V-R處理土壤中的真菌多樣性指數(shù)最高,優(yōu)勢度指數(shù)最低,說明T-V-R輪作的土壤適合多種真菌的繁殖生長,種群數(shù)量增加。多種真菌共同存在,互相制約,可防止病原真菌過度繁殖,降低作物發(fā)生真菌病害幾率。從作物產(chǎn)量和產(chǎn)值、土壤有機(jī)質(zhì)、養(yǎng)分和真菌種群結(jié)構(gòu)看,T-V-R優(yōu)于其它3種輪作模式,值得推廣應(yīng)用。

    關(guān)鍵詞:輪作;土壤;養(yǎng)分;真菌;多樣性

    同一作物或近緣作物長期連續(xù)種植會(huì)造成土壤養(yǎng)分異常積累或過度消耗,病原微生物迅速繁衍,化感物質(zhì)積累,土壤微生物種群結(jié)構(gòu)失衡,作物產(chǎn)量品質(zhì)降低和病害嚴(yán)重發(fā)生等[1]。因此,多種作物輪作輔以休閑一般優(yōu)于種植單一作物,研究不同輪作模式對土壤物理、化學(xué)、生物學(xué)性質(zhì)的影響,可為合理輪作提供科學(xué)依據(jù)。

    在我國熱帶和亞熱帶地區(qū),通常采用一年三熟或兩年五熟的高強(qiáng)度土壤利用模式,需要大量施用化肥,導(dǎo)致土壤有機(jī)質(zhì)減少,結(jié)構(gòu)破壞,病害嚴(yán)重[2]。馬鈴薯長期連作后,根際土壤中的細(xì)菌數(shù)顯著低于其他種植模式,但真菌則相反,作物發(fā)生真菌病害的風(fēng)險(xiǎn)提高[3]。玉米連作年限不同,根際微生物的數(shù)量、種群和多樣性也不一樣。隨著連作時(shí)間延長,它們的數(shù)量和種群減少,多樣性逐漸降低[4]。烤煙長期連作,提高黑脛病和青枯病的發(fā)病率,降低煙葉產(chǎn)量和質(zhì)量[5]。相反,輪作則能有效中斷病原菌寄主,降低作物病害發(fā)生率[5- 7]。值得注意的是,輪作模式不同,對土壤物理、化學(xué)、生物學(xué)性質(zhì)的影響也不一樣。連續(xù)5a實(shí)施水稻-油菜輪作,土壤有機(jī)質(zhì)含量顯著高于水稻-小麥和水稻-黑麥草輪作[8]。在我國東北黑土上,實(shí)施黃瓜連作、小麥-黃瓜和毛苕子-黃瓜輪作,輪作顯著增加根際細(xì)菌種類,但減少結(jié)瓜后期的真菌種群,尤以毛苕子-黃瓜輪作最為顯著[9]。微生物是土壤的重要組成部分,釋放多種酶類,參與土壤有機(jī)質(zhì)降解、腐殖質(zhì)合成和養(yǎng)分循環(huán)等生物化學(xué)過程[10]。土壤微生物數(shù)量和種群隨種植植物的種類不同而發(fā)生變化,作物根系分泌物和凋落物對土壤微生物的種群結(jié)構(gòu)產(chǎn)生重要影響[11],并與土傳病害密切相關(guān)[12]。

    在微生物生態(tài)學(xué)研究中,利用Roche Genome Sequencer FLX測序平臺(tái)進(jìn)行454高通量測序是重要的分子生物學(xué)手段。該方法通過測定微生物保守性較高的16S rDNA和18S rDNA序列,并與基因庫中的已知序列進(jìn)行對比,能準(zhǔn)確、快速、靈敏地確定待測微生物屬種。自454高通量測序技術(shù)應(yīng)用以來,已檢測出數(shù)百至數(shù)千種土壤微生物,是傳統(tǒng)培養(yǎng)方法所獲得微生物數(shù)量的數(shù)十倍甚至數(shù)百倍。此外,該方法還能檢測出不可培養(yǎng)的微生物[13]。云南省是我國烤煙主產(chǎn)區(qū),烤煙輪作比例達(dá)80%以上,煙地夏季主要輪作玉米和水稻,分別占種植面積的64.7%和27.6%,冬季一般種植油菜和綠肥或休閑[14]。選擇科學(xué)的輪作模式有益于減輕連作障礙,是保持土地生產(chǎn)力和保障整個(gè)農(nóng)業(yè)生產(chǎn)健康和持續(xù)發(fā)展的重要措施[15]。因此,本研究以云南省煙草農(nóng)業(yè)科學(xué)研究院的長期輪作試驗(yàn)為對象,采用化學(xué)、酶學(xué)方法及454高通量測序技術(shù),揭示了云南煙區(qū)主要輪作模式對作物產(chǎn)量和產(chǎn)值、土壤養(yǎng)分、酶活性及真菌種群結(jié)構(gòu)的影響,為合理選擇烤煙輪作模式奠定基礎(chǔ)。

    1材料與方法

    1.1試驗(yàn)地概況

    云南省煙草農(nóng)業(yè)科學(xué)研究院研究和試驗(yàn)基地位于北緯24°14′,東經(jīng) 102°30′,海拔1680 m,年均溫度15.9℃,年降雨量918mm,雨季(4—9月)降雨量占全年的79.5%,年日照時(shí)數(shù)2072h。供試土壤為云南省典型、具有代表性的砂質(zhì)紅壤,基礎(chǔ)土壤的pH6.4,有機(jī)質(zhì)10.70 g/kg,全氮0.54 g/kg,全磷0.11 g/kg,全鉀6.43 g/kg,有效氮82.0 mg/kg,有效磷 9.01 mg/kg,有效鉀160.0 mg/kg。

    1.2試驗(yàn)處理

    試驗(yàn)始于1998年,設(shè)置云南普遍采用的4種種植模式,包括烤煙-休閑-玉米(T-B-M)、烤煙-油菜-玉米(T-C-M)、烤煙-油菜-水稻(T-C-R)和烤煙-苕子-水稻(T-V-R)等。小區(qū)面積28m2,重復(fù)3次,隨機(jī)區(qū)組排列??緹?、玉米和水稻一般在5月初移栽,苕子和油菜秋末播種。根據(jù)當(dāng)?shù)卮筇锷a(chǎn)情況,烤煙施用純氮105kg/hm2,P2O574 kg/hm2,K2O 263 kg/hm2;油菜施用純氮90kg/hm2,P2O530 kg/hm2;玉米和水稻施用純氮112.5kg/hm2,P2O540 kg/hm2;苕子施用純氮75kg/hm2。此外,水稻、玉米和油菜秸稈全部切碎還田;綠肥收割多次,未徹底收割的殘樁全部翻壓入土。在烤煙種植季節(jié),基肥由N∶P2O5∶K2O=10∶10∶25烤煙專用肥提供,施氮量占施肥總量的70%;剩余30%的氮素由烤煙專用追肥(N∶P2O5∶K2O = 10∶0∶25)提供,在移栽后7—10d和30—40d均分兩次穴施。在玉米、油菜及苕子種植季節(jié),化學(xué)氮磷分別由尿素和過磷酸鈣提供。其中,玉米和油菜70%的氮和全部磷肥作基肥,剩余30%的氮分別在油菜云苔期和玉米大喇叭口期作追肥穴施,苕子全施基肥。

    1.3土樣采集與分析

    在不同輪作模式中,栽種或休閑是其中的重要環(huán)節(jié),它們均可能改變土壤理化生物學(xué)性質(zhì),進(jìn)而影響后季作物-烤煙生長發(fā)育和產(chǎn)量品質(zhì)。因此,試驗(yàn)于2014年春季油菜和苕子旺長期,采集0—20cm的耕層土壤(含休閑處理),揀去雜物。部分土壤立即液氮冷凍備測微生物量碳氮和真菌18S rDNA序列,另取部分土壤風(fēng)干,常規(guī)分析土壤pH、有機(jī)質(zhì)和堿解氮、有效磷和有效鉀[16]。微生物量碳氮采用氯仿熏蒸-0.5mol/L K2SO4提取,用K2Cr2O7氧化法測碳和腚酚藍(lán)比色法測氮[17]。分別用3,5-二硝基水楊酸比色法、磷酸苯二鈉比色法、腚酚藍(lán)比色法、高錳酸鉀滴定容量法和TTC比色法測定土壤蔗糖酶、磷酸酶、脲酶、過氧化氫酶和脫氫酶活性[18]。

    提取土壤樣品中的真菌18S rDNA,接著采用真菌18S區(qū)段中V5—V7區(qū)的通用引物0817F(5′- TTAGCATGGAATAATRRAATAGGA-3′)和1196R(5′- TCTGGACCTGGTGAGTTTCC-3′)進(jìn)行擴(kuò)增。再參照454高通量測序方法,純化、定量和均一化真菌18S rDNA,送上海美吉生物科技有限公司利用Roche Genome Sequencer FLX測序平臺(tái)進(jìn)行測序[19]。然后,對有效序列進(jìn)行去雜、修剪、去除嵌合體序列等過濾處理,得到優(yōu)化序列,通過聚類分析形成分類單元(operational taxonomic units, OTUs),采用BLAST程序?qū)Ρ菺enBank(http://ncbi.nlm.nih.gov)中的已知序列,根據(jù)97%的相似度確定18S rRNA基因序列對應(yīng)的真菌屬(種)名稱。

    于2014年烤煙收獲季節(jié),分葉位采收煙葉,三段式烘烤,按煙草行業(yè)標(biāo)準(zhǔn)分級(jí)定價(jià),記錄上、中、下等煙葉產(chǎn)量,均價(jià)和產(chǎn)值。

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

    利用土壤真菌屬(種)數(shù)(OTUs)和18S rDNA序列數(shù)(Reads)計(jì)算土壤真菌的種群特征值,包括多樣性指數(shù)和優(yōu)勢度指數(shù)。

    Shannon-Wiener多樣性指數(shù)H的計(jì)算公式為:

    式中,Sobs為得到的OTUs的總數(shù)量;ni為真菌i的OTUs數(shù)量;N為得到的真菌的總數(shù)量。

    Simpson優(yōu)勢度指數(shù)D的計(jì)算公式為:

    式中,Sobs為得到的OTUs的總數(shù)量;ni為真菌i的OTUs數(shù)量;N為得到的真菌的總數(shù)量[20- 21]。

    真菌豐富度為某種真菌的18S rDNA讀數(shù)占真菌18S rDNA總讀數(shù)的百分?jǐn)?shù)[22- 23]。

    試驗(yàn)數(shù)據(jù)用Excel進(jìn)行基本計(jì)算,SPSS16.0軟件進(jìn)行統(tǒng)計(jì)分析,差異顯著性水平為P≤ 0.05。

    2試驗(yàn)結(jié)果

    2.1煙葉產(chǎn)量、產(chǎn)值和上中等煙比例

    在T-V-R處理的土壤中,煙葉的產(chǎn)量、產(chǎn)值和上中等煙比例都為最高,煙葉產(chǎn)量和產(chǎn)值相比T-C-M處理分別增加了22.88%和25.30%;與T-C-R處理相比,上中等煙比例提高了3.62%(表1)。

    表1 不同輪作模式對烤煙產(chǎn)量,產(chǎn)值和上中等煙比例的影響

    T-B-M:烤煙-休閑-玉米輪作 flue-cured tobacco-blank-maize crop rotation;T-C-M:烤煙-油菜-玉米輪作 flue-cured tobacco-canola-maize crop rotation;T-C-R:烤煙-油菜-水稻輪作 flue-cured tobacco-canola-rice crop rotation;T-V-R:烤煙-苕子-水稻輪作 flue-cured tobacco-vetch-rice crop rotation; 表中數(shù)據(jù)為平均數(shù)±標(biāo)準(zhǔn)差;在同一列中,有不同小寫字母者表示差異顯著,(P≤ 0.05)

    2.2土壤pH、有機(jī)質(zhì)與有效養(yǎng)分

    在不同輪作的土壤中,pH、有機(jī)質(zhì)和有效養(yǎng)分含量都不同程度的增加或無顯著變化(表2)。與原始土壤相比,T-B-M和T-C-M處理使土壤pH值顯著增加,最高為6.4,T-V-R和T-C-R處理對土壤的pH值無顯著影響;T-C-M、T-C-R和T-V-R分別使有機(jī)質(zhì)增加25.49%、16.23%和45.11%,T-B-M則無顯著變化;T-V-R顯著增加土壤堿解氮含量,為原始土壤的1.32倍;4種輪作處理均提高有效磷含量,為原始土壤的1.12—3.67倍;T-B-M的有效鉀含量最高,其余3種輪作的有效鉀含量與原始土壤相似。

    表2 不同輪作模式對土壤pH、有機(jī)質(zhì)和有效養(yǎng)分的影響

    表中數(shù)據(jù)為平均數(shù)±標(biāo)準(zhǔn)差;在同一列中,有不同小寫字母者表示差異顯著,(P≤ 0.05)

    2.3土壤酶活

    不同輪作模式對土壤酶活的影響各不相同(表3)。在T-V-R處理的土壤中,蔗糖酶、磷酸酶、脫氫酶、過氧化氫酶以及脲酶活性最高;在T-B-M、T-C-M和T-C-R處理的土壤中,蔗糖酶和脲酶活性分別變化于36.07—39.05和0.63—0.65之間,處理之間無顯著差異;在T-C-M處理的土壤中,磷酸酶和脫氫酶活性低于其余3種輪作處理;T-B-M和T-C-R處理的土壤中,過氧化氫酶活性無顯著差異,但顯著低于T-C-M和T-V-R。

    表3 不同種植模式對土壤酶活性的影響

    表中數(shù)據(jù)為平均數(shù)±標(biāo)準(zhǔn)差;在同一列中,有不同小寫字母者表示差異顯著(P≤ 0.05)

    2.4微生物生物量

    圖1可見,在T-V-R處理的土壤中,微生物碳含量最高,為170.0mg/kg;T-C-M和T-C-R處理微生物碳含量次之但無顯著差異,為117.2—119.4mg/kg;T-B-M最低,僅為61.9mg/kg。微生物氮含量的變化趨勢則有所不同,T-V-R和T-C-M含量最高,達(dá)到8.9—9.4mg/kg,其次為T-B-M,T-C-R含量最低。實(shí)施不同輪作處理后,土壤微生物量碳氮比變化于7.5—19.1之間。

    圖1 種植模式對土壤微生物碳氮含量的影響Fig.1 Influence of rotation treatments on microbial C and N in soilT-B-M:烤煙-休閑-玉米輪作 flue-cured tobacco-blank-maize crop rotation;T-C-M:烤煙-油菜-玉米輪作 flue-cured tobacco-canola-maize crop rotation;T-C-R:烤煙-油菜-水稻輪作 flue-cured tobacco-canola-rice crop rotation;T-V-R:烤煙-苕子-水稻輪作 flue-cured tobacco-vetch-rice crop rotation

    圖2 不同輪作處理的土壤中,真菌稀釋性曲線 Fig.2 Fungi rarefaction curves in soil under different rotation treatments

    2.5土壤真菌2.5.1真菌稀釋曲線

    隨機(jī)抽取測序樣品中的18S rDNA讀數(shù)(reads),以真菌屬(種)數(shù)(OTUs)為縱坐標(biāo),18S rDNA讀數(shù)為橫坐標(biāo),獲得稀釋曲線(圖2)[24]。結(jié)果表明,抽樣讀數(shù)大約在1500以下時(shí),細(xì)菌種(屬)類數(shù)迅速增加;讀數(shù)在1500—4000之間,細(xì)菌種(屬)數(shù)緩慢增加;讀數(shù)超過4000之后,其種(屬)類數(shù)的增長逐漸趨于平緩。但是,細(xì)菌稀釋曲線因輪作處理不同而異,T-V-R最高,T-B-M次之,T-C-R和T-C-M最低,且?guī)缀踔丿B。

    2.5.2真菌門類

    在T-B-M、T-C-M、T-C-R和T-V-R處理的土壤中,454-高通量測序分別獲得了13097,11345,12939和13763個(gè)18S rDNA 序列,分別代表530、378、395和581種真菌(OTUs),歸屬于子囊菌門(Ascomycota)、壺菌門(Chytridiomycota)、擔(dān)子菌門(Basidiomycota)、接合菌門(Zygomycota)和尚待鑒定的類型 (表4)。其中,T-B-M、T-C-M和T-C-R 3個(gè)處理中待定真菌OTUs數(shù)量最多,占總OTUs數(shù)的46.03%—53.40%,其次是子囊菌門,比例為34.53%—37.22%;T-V-R處理中子囊菌OTUs數(shù)量最多,占總數(shù)的45.78%,其次是未知類型,占總數(shù)的39.76%。4個(gè)處理中,接合菌門、擔(dān)子菌門和壺菌門合計(jì)約占總量的12.08%—15.61%。

    表4 不同輪作處理土壤中真菌OTUs數(shù)量

    *表中數(shù)據(jù)為平均數(shù);在同一行中,有不同小寫字母者表示差異顯著(P≤ 0.05)

    2.5.3優(yōu)勢真菌

    在各輪作處理的土壤中,前15種優(yōu)勢真菌的豐富度合計(jì)占總量的29.46%—62.86%,尚待鑒定的真菌豐富度高達(dá)6.21%—45.75%。其中,在T-C-M和T-C-R處理中,待定菌株最多,高達(dá)9株,累計(jì)豐富度分別為45.75%和41.09%;其次是T-B-M處理,待定真菌有7株,豐富度為21.92%;在T-V-R處理中,待定菌株數(shù)量最少,僅有5株,豐富度為6.21%(表5)。

    在前15種優(yōu)勢菌株中,4種處理之間真菌屬(種)的相似性極低,無共同存在的真菌屬(種),且每種處理特有的真菌種類均超過9種(表5)。此外,待定真菌-4(unclassified fungi-4)共同存在于T-B-M、T-C-R和T-V-R 3個(gè)處理中;枝孢菌(Cladosporium)同時(shí)存在于T-C-R和T-V-R處理中;T-B-M和T-V-R處理共有的真菌是尖孢鐮刀菌(Fusarium_oxysporum)、接合菌-1(Zygomycota-1)、鐮刀菌(Fusarium)和鏈格孢菌(Alternaria);T-C-M和T-C-R處理共有的真菌為待定真菌-7(unclassified fungi-7)、待定真菌-10(unclassified fungi-10)、待定真菌-12(unclassified fungi-12)、待定真菌-14(unclassified fungi-14)。

    2.5.4真菌群落特征值

    在T-V-R的土壤中,真菌多樣性指數(shù)最高,高達(dá)5.34;T-B-M次之,為4.97;T-C-M和T-C-R最低,變化于4.10—4.28之間。但是,優(yōu)勢度指數(shù)T-C-R最高,達(dá)到0.04;T-V-R最低,僅0.01(表6)。

    3討論

    有機(jī)質(zhì)是土壤的重要成分,顯著影響土壤物理、化學(xué)、生物學(xué)性質(zhì)。秸稈還田既可補(bǔ)充土壤有機(jī)質(zhì)和作物需要的養(yǎng)分,又能改良土壤理化生物學(xué)性質(zhì)[25]。在16年的長期輪作過程中,每年都有大量的秸稈還田,土壤有機(jī)質(zhì)和有效養(yǎng)分含量提高或無顯著變化,說明在高強(qiáng)度輪作種植條件下,輔以秸稈還田或冬季休閑均可保持或提高土壤生產(chǎn)力。此外,各輪作處理的有效磷含量遠(yuǎn)遠(yuǎn)超過原始土壤,故適量減施磷肥很有必要。土壤pH值變化于5.6—6.4,都處于適合多種作物生長的pH值范圍內(nèi)[26]。綜合作物產(chǎn)量、產(chǎn)值、土壤pH、有機(jī)質(zhì)和有效養(yǎng)分的狀況,以T-V-R種植模式最佳。

    土壤有機(jī)質(zhì)是微生物的碳源和氮源[27]。在輪作處理的土壤中,每年都實(shí)施秸稈還田,各種作物秸稈腐解所釋放的有機(jī)質(zhì)組成具有多樣性,可滿足多種微生物,尤其是異養(yǎng)型微生物-真菌對碳源和養(yǎng)分的需要。在實(shí)施T-V-R輪作的處理中,水旱輪作,嫌/好氣交替,創(chuàng)造了適合多種微生物繁衍的不同土壤環(huán)境;加之烤煙、水稻和苕子的近緣性小,有機(jī)成分差異較大,可滿足不同微生物的碳源和營養(yǎng)需要,有益于它們的繁殖生長。此外,油菜為直根系,而苕子為須根系,后者根系發(fā)達(dá),能分泌大量的單(多)糖、氨基酸和有機(jī)酸等碳水化合物,向根際微生物提供碳源,促進(jìn)其生長繁殖[28]。因此,在T-V-R土壤中,微生物量碳氮最高,18S rNDA讀數(shù)最大,真菌屬(種)類最多,土壤酶活性最強(qiáng)。脫氫酶、脲酶和磷酸酶活性與土壤微生物量呈顯著正相關(guān),推測與土壤微生物的分泌作用有關(guān),類似前人研究結(jié)果[29- 30]。

    表5 在不同輪作土壤中,15種優(yōu)勢真菌的豐富度

    表6各輪作處理土壤真菌群落特征值

    Table 6Characteristic values of fungal communities in soil under different rotation treatments

    種植模式Rotationmodals多樣性指數(shù)Diversityindex優(yōu)勢度指數(shù)DominantindexT-B-M4.97b0.02cT-C-M4.28c0.03bT-C-R4.10c0.04aT-V-R5.34a0.01d

    表中數(shù)據(jù)為平均數(shù);在同一列中,有不同小寫字母者表示差異顯著,(P≤ 0.05)

    在不同輪作處理的土壤中,微生物碳氮比相差2.5倍以上,意味著土壤微生物的組成結(jié)構(gòu)各不相同。本文采用454高通量測序方法,揭示了真菌從門類到屬(種)的組成、優(yōu)勢菌株、豐富度和多樣性特征等。結(jié)果表明,稀釋性曲線呈迅速增加-緩慢增加-逐漸平緩的變化趨勢,說明18S rDNA的樣本量和測序數(shù)據(jù)合理,準(zhǔn)確地反映了樣本中的真菌數(shù)及組成情況;稀釋曲線的峰值高低不一,意味著輪作處理影響了土壤真菌的多樣性。一般認(rèn)為,多樣性指數(shù)表示生物群落中的物種多寡,數(shù)值愈大表示群落中的物種越豐富;優(yōu)勢度指數(shù)越大,生物群落內(nèi)的奇異度越高,優(yōu)勢種群突出[31]。生物種群豐富和多樣性指數(shù)高是生態(tài)環(huán)境健康穩(wěn)定的重要表現(xiàn)[32]。在健康穩(wěn)定的生態(tài)環(huán)境中,生物多樣性指數(shù)較高,優(yōu)勢度指數(shù)較低。T-V-R顯著提高真菌群落的多樣性指數(shù),但降低優(yōu)勢度指數(shù)。因此,T-V-R處理的土壤生態(tài)環(huán)境較好,使之適合多種真菌的繁殖生長,種群數(shù)量增加。多種真菌共同存在,互相制約,可防止某些病原真菌過度繁殖,抑制作物病害的發(fā)生。值得注意的是,454高通量測序表明,在不同種植模式的土壤中,真菌的豐富度因輪作處理和它們的種類不同而異,存在378—581種(屬)的真菌,遠(yuǎn)遠(yuǎn)超過了目前的常規(guī)培養(yǎng)、磷脂脂肪酸和PCR-DGGE能達(dá)到的水平。此外,用形態(tài)和生理生化方法鑒定微生物需要豐富的經(jīng)驗(yàn),程序復(fù)雜,耗時(shí)長,容易出現(xiàn)誤判。說明454高通量測序檢測靈敏度高,鑒定準(zhǔn)確性可靠,更能深入研究土壤真菌的組成和群落結(jié)構(gòu)。本項(xiàng)研究表明,18S rDNA讀數(shù)和真菌OTUs以T-V-R最高,T-C-M最低,說明T-V-R土壤中的真菌數(shù)和屬(種)類最多。

    除了待定真菌外,子囊菌占各處理真菌的絕大多數(shù),都為40%左右甚至超過40%以上。經(jīng)16a不同輪作后,4種輪作處理的土壤中,前15種優(yōu)勢菌株的相似性極低,說明長期輪作極大地改變了土壤中真菌的組成及類型。需要指出的是,土壤真菌參與土壤有機(jī)質(zhì)礦化,分解進(jìn)入土壤的作物殘?bào)w和有機(jī)肥中的有機(jī)質(zhì),供給植物營養(yǎng)[33]。但是,大多數(shù)真菌同時(shí)也為植物病原菌,引起多種真菌病害。例如,子囊菌引起根腐、莖腐、果(穗)腐和枝枯等;鐮刀菌普遍為植物真菌性病害;葡萄肉坐菌幾乎能危害所有果樹,如香蕉蕉腐病、桃樹流膠病等;座囊菌是我國主要的果樹病菌,引起香蕉葉斑病,梨葉斑病,葡萄黑痘病等;鏈格孢菌95%以上兼性寄生于植物,引起多種植物病害,如小麥葉枯病,玉米大斑病,茄子早疫病,白菜黑斑病等;炭角菌廣泛分布于全球各地,在熱帶和亞熱帶地區(qū)尤為常見,能引起桑、茶、櫟、蘋果根腐病和黑腐病等[34]。從防治作物土傳真菌病害的角度看,合理輪作是完全可能和必要的。

    4結(jié)論

    總之,4種輪作方式均有益于保持或提高土壤肥力和土地生產(chǎn)力。從作物產(chǎn)量、產(chǎn)值、土壤有機(jī)質(zhì)、有效養(yǎng)分、酶活性、微生物量、真菌種群結(jié)構(gòu)和多樣性等方面看,烤煙-苕子-水稻輪作最佳,是值得推廣的一種種植模式。

    參考文獻(xiàn)(References):

    [1]王飛, 李世貴, 徐鳳花, 顧金剛. 連作障礙發(fā)生機(jī)制研究進(jìn)展. 中國土壤與肥料, 2013, (5): 6- 13.

    [2]張北贏, 陳天林, 王兵. 長期施用化肥對土壤質(zhì)量的影響. 中國農(nóng)學(xué)通報(bào), 2010, 26(11): 182- 187.

    [3]劉高遠(yuǎn), 郭天文, 譚雪蓮, 馬得禎, 張國宏, 陳光榮. 不同栽培方式下馬鈴薯土壤微生物區(qū)系的動(dòng)態(tài)變化. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào), 2014, 33(4): 19- 24.

    [4]邢會(huì)琴, 肖占文, 閆吉智, 馬建倉, 孟嫣. 玉米連作對土壤微生物和土壤主要養(yǎng)分的影響. 草業(yè)科學(xué), 2011, 28(10): 1777- 1780.

    [5]晉艷, 楊宇虹, 段玉琪, 孔光輝. 烤煙輪作、連作對煙葉產(chǎn)量質(zhì)量的影響. 西南農(nóng)業(yè)學(xué)報(bào), 2004, 17(增刊): 267- 271.

    [6]Zhang W, Long X Q, Huo X D,Chen Y F, Lou K. 16S rRNA-based PCR-DGGE analysis of actinomycete communities in fields with continuous cotton cropping in Xinjiang, China. Soil Microbiology, 2013, 66(2): 385- 393.

    [7]鄒雨坤, 張靜妮, 楊殿林, 陳秀蓉, 張?zhí)烊? 趙建寧, 趙帥. 不同利用方式下羊草草原土壤生態(tài)系統(tǒng)微生物群落結(jié)構(gòu)的PLFA分析. 草業(yè)學(xué)報(bào), 2011, 20(4): 27- 33.

    [8]王先華. 稻田不同輪作方式對培肥地力的作用. 耕作與栽培, 2002, (6): 9- 10.

    [9]于高波, 吳鳳芝, 周新剛. 小麥、毛苕子與黃瓜輪作對土壤微生態(tài)環(huán)境及產(chǎn)量的影響. 土壤學(xué)報(bào), 2011, 48(1): 175- 184.

    [10]Kennedy A C, Smith K L. Soil microbial diversity and the sustainability of agricultural soils. Plant and Soil, 1995, 170(1): 75- 86.

    [11]張鼎華, 葉章發(fā), 李寶福. 杉木、馬尾松輪作對林地土壤肥力和林木生長的影響. 林業(yè)科學(xué), 2011, 37(5): 10- 15.

    [12]袁龍剛, 張軍林, 張朝陽, 別國勇. 連作對辣椒根際土壤微生物區(qū)系影響的初步研究. 陜西農(nóng)業(yè)科學(xué), 2006, (2): 49- 50.

    [13]段曌, 肖煒, 王永霞, 賴泳紅, 崔曉龍. 454測序技術(shù)在微生物生態(tài)學(xué)研究中的應(yīng)用. 微生物學(xué)雜志, 2011, 31(5): 76- 81.

    [14]李天福, 王樹會(huì), 王彪, 邵巖. 云南烤煙輪作現(xiàn)狀與對策分析. 煙草農(nóng)業(yè)科學(xué), 2006, 2(2): 216- 222.

    [15]郁飛燕, 李友軍. 輪作對土壤肥力的影響. 科技傳播, 2011, (16): 60- 60.

    [16]楊劍虹, 王成林, 代亨林. 土壤農(nóng)化分析與環(huán)境監(jiān)測. 北京: 中國大地出版社, 2008: 26- 75.

    [17]Frosteg?rd ?, Tunlid A, B??th E. Phospholipid fatty acid composition, biomass, and activity of microbial communities from two soil types experimentally exposed to different heavy metals. Applied and Environmental Microbiology, 1993, 59(11): 3605- 3617.

    [18]關(guān)松蔭. 土壤酶及其研究法. 北京: 農(nóng)業(yè)出版社, 1986: 274- 339.

    [19]Herlemann D P R, Lundin D, Labrenz M, Jürgens K, Zheng A L, Aspeborg H, Andersson A F. Metagenomic de novo assembly of an aquatic representative of the verrucomicrobial class spartobacteria. MBIO, 2013, 4(3): e00569- 12.

    [20]Krebs C J. Ecological Methodology. New York:Benjamin Cummings, 1989: 654- 654.

    [21]Schloss P D, Gevers D, Westcott S L. Reducing the effects of PCR amplification and sequencing artifacts on 16S rDNA-based studies. PLoS ONE, 2011, 6(12): e27310.

    [22]Roesch L F W, Fulthorpe R R, Riva A, Casella G, Hadwin A K M, Kent A D, Daroub S H, Camargo F A O, Farmerie W G, Triplett E W. Pyrosequencing enumerates and contrasts soil microbial diversity. The ISME Journal, 2007, 1(4): 283- 290.

    [23]Lim Y M, Kim B K, Kim C, Jung H S, Kim B S, Lee J H, Chun J. Assessment of soil fungal communities using pyrosequencing. The Journal of Microbiology, 2010, 48(3): 284- 289.

    [24]Amato K R, Yeoman C J, Kent A, Righini N, Carbonero F, Estrada A, Gaskins H R, Stumpf R M, Yildirim S, Torralba M, Gillis M, Wilson B A, Nelson K E, White B A, Leigh S R. Habitat degradation impacts black howler monkey (Alouattapigra) gastrointestinal microbiomes. The ISME Journal, 2013, 7: 1344- 1353.

    [25]徐國偉, 李帥, 趙永芳, 陳明燦, 李友軍. 秸稈還田與施氮對水稻根系分泌物及氮素利用的影響研究. 草業(yè)學(xué)報(bào), 2014, 23(2): 140- 146.

    [26]李念勝, 王樹聲. 土壤pH值與烤煙質(zhì)量. 中國煙草, 1986, (2): 12- 14.

    [27]Katarina H. Soil microbial community structure in relation to vegetation management on former agricultural land. Soil Biology & Biochemistry, 2002, 34(9): 1299- 1307.

    [28]焦彬. 中國綠肥. 北京: 農(nóng)業(yè)出版社, 1986: 97- 110.

    [29]袁亮. 設(shè)施栽培土壤微生物量和酶活性的變化規(guī)律及其與土壤肥力的關(guān)系[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2007.

    [30]單貴蓮, 初曉輝, 羅富成, 馬玉寶, 李臨杭, 陳功. 圍封年限對典型草原土壤微生物及酶活性的影響. 草原與草坪, 2012, 32(1): 1- 6.

    [31]Hunter P R, Gaston M A. Numerical index of the discriminatory ability of typing systems: An application of Simpson′s index of diversity. Journal of Clinical Microbiology, 1988, 26(11): 2465- 2466.

    [32]高東, 何霞紅. 生物多樣性與生態(tài)系統(tǒng)穩(wěn)定性研究進(jìn)展. 生態(tài)學(xué)雜志, 2010, 29(12): 2507- 2513.

    [33]徐嚴(yán). 若爾蓋高原濕地土壤真菌的初步研究[D]. 成都: 四川農(nóng)業(yè)大學(xué), 2008.

    [34]康振生. 我國植物真菌病害的研究現(xiàn)狀及發(fā)展策略. 植物保護(hù), 2010, 36(3): 9- 12.

    Influence of crop rotation on enzyme activities and fungal communities in flue-cured tobacco soil

    CHEN Danmei1, DUAN Yuqi2, YANG Yuhong2, JIN Yan2, YUAN Ling1,*

    1CollegeofResourcesandEnvironment,SouthwestUniversity,Chongqing400716,China2YunnanTobaccoAgricultureResearchInstitute,Kunming650031,China

    Abstract:Rotation practices have significant influences on the chemical, physical, and biological properties of soils, as well as on the occurrence of crop diseases caused by soil fungi. A long-term rotation experiment was set up in a red soil area of the Yunnan province with four rotation models: flue-cured tobacco-blank-maize (T-B-M), flue-cured tobacco-canola-maize (T-C-M), flue-cured tobacco-canola-rice (T-C-R), and flue-cured tobacco-vetch-rice (T-V-R). Chemical and enzymatic analysis and 454 pyrosequencing were used to study changes in nutrients, enzymes activities, and fungal colony structures in order to provide a scientific basis for the rational use of rotation of flue-cured tobacco. After 16 years of crop rotation, the soil pH changed from 5.6 to 6.4, which is still suitable for the cultivation of flue-cured tobacco, canola, rice, and maize. The yield, output value, and quality of flue-cured tobacco were increased significantly, and organic matter was increased by 45.11% under the condition of T-V-R treatment. Similar changes were also found in available nitrogen and phosphorus, soil microbial biomass carbon and nitrogen, and activities of sucrase, dehydrogenase, urease, phosphatase, and catalase. This suggests improvement of the soil ecological environment, promotion of microbial reproduction, and increases in microbial numbers and activities. The reads of fungal 18S rDNA sequences were 13097 in the soil with T-B-M, 11345 with T-C-M, 12939 with T-C-R, and 13763 with T-V-R, representing 530378395, and 581 fungal genera or species, respectively. All fungi belonged to Ascomycota, Basidiomycota, Zygomycota, Chytridiomycota, or unidentified groups, with the majority belonging to Ascomycota. The abundance of the top 15 predominant fungi accounted for 29.46%—62.86% of the total abundance, and less similarity was found in the tested soils, indicating variation in fungal colony structure in the soils with variable rotation models. T-V-R showed the highest Shannon diversity index but the lowest Simpson index of fungal communities, implying that the soils with T-V-R were favorable for the reproduction and growth of multiple fungi and increment of their groups. The co-presence of a variety of fungi in a soil could result in reciprocal inhibition and prevent their over-reproduction, thereby decreasing the possibility of the crops being infected by pathogenic fungi. With respect to crop yield, organic matter in soils, nutrients, and fungal community structure, T-V-R was superior to the other three rotation models, which could be encouraged in the tobacco-growing areas in Yunnan Province.

    Key Words:rotation; soil; nutrient; fungus; diversity

    基金項(xiàng)目:云南省煙草公司項(xiàng)目(2013YN11); 國家煙草專賣局項(xiàng)目(110201302016); 科技部“973”項(xiàng)目(2013CB127405)

    收稿日期:2014- 11- 05; 網(wǎng)絡(luò)出版日期:2015- 08- 24

    *通訊作者

    Corresponding author.E-mail: lingyuanh@aliyun.com

    DOI:10.5846/stxb201411052183

    陳丹梅, 段玉琪,楊宇虹, 晉艷, 袁玲.輪作模式對植煙土壤酶活性及真菌群落的影響.生態(tài)學(xué)報(bào),2016,36(8):2373- 2381.

    Chen D M, Duan Y Q, Yang Y H, Jin Y, Yuan L.Influence of crop rotation on enzyme activities and fungal communities in flue-cured tobacco soil.Acta Ecologica Sinica,2016,36(8):2373- 2381.

    猜你喜歡
    多樣性輪作真菌
    農(nóng)民瓜菜“接茬輪作”增效益
    高等大型真菌與人類
    科學(xué)(2020年2期)2020-08-24 07:56:56
    真菌造房子
    輪作休耕,讓土地歇一會(huì)兒
    耕地輪作效益調(diào)查報(bào)告
    新時(shí)期群文輔導(dǎo)工作的特征
    海洋微生物次生代謝的生物合成機(jī)制
    科技資訊(2016年19期)2016-11-15 10:39:12
    舞蹈表演的表現(xiàn)形式多樣性研究
    人間(2016年27期)2016-11-11 16:27:23
    水磨地區(qū)蕨類植物多樣性調(diào)查分析
    艾滋病合并侵襲性真菌感染的診治
    人人妻人人澡人人爽人人夜夜| 91aial.com中文字幕在线观看| 亚洲欧美成人综合另类久久久| √禁漫天堂资源中文www| 9191精品国产免费久久| 飞空精品影院首页| 高清视频免费观看一区二区| freevideosex欧美| 精品久久国产蜜桃| 性色avwww在线观看| 蜜桃国产av成人99| 91精品三级在线观看| 午夜视频国产福利| 欧美精品人与动牲交sv欧美| 欧美人与性动交α欧美精品济南到 | 校园人妻丝袜中文字幕| 国产精品国产三级国产专区5o| 日韩视频在线欧美| 美女脱内裤让男人舔精品视频| 中文字幕人妻熟女乱码| 亚洲色图综合在线观看| 黄色 视频免费看| 亚洲精华国产精华液的使用体验| 午夜日本视频在线| 九色成人免费人妻av| 中文字幕精品免费在线观看视频 | 日韩中文字幕视频在线看片| av不卡在线播放| 蜜桃国产av成人99| 欧美精品高潮呻吟av久久| 巨乳人妻的诱惑在线观看| 丁香六月天网| 午夜福利乱码中文字幕| 深夜精品福利| 人人妻人人澡人人爽人人夜夜| 国产高清不卡午夜福利| 寂寞人妻少妇视频99o| 天堂中文最新版在线下载| 免费观看av网站的网址| 一级片免费观看大全| 高清黄色对白视频在线免费看| 午夜福利网站1000一区二区三区| freevideosex欧美| 又大又黄又爽视频免费| 少妇的丰满在线观看| 国产成人精品婷婷| 久久影院123| 美女xxoo啪啪120秒动态图| 国产精品女同一区二区软件| 少妇猛男粗大的猛烈进出视频| av卡一久久| videossex国产| 欧美国产精品va在线观看不卡| av黄色大香蕉| 免费女性裸体啪啪无遮挡网站| 亚洲av成人精品一二三区| 亚洲色图综合在线观看| 日韩制服丝袜自拍偷拍| 亚洲成色77777| 久久久久精品久久久久真实原创| 在线观看免费高清a一片| 亚洲精品456在线播放app| 亚洲,欧美,日韩| 大香蕉久久成人网| 亚洲av电影在线观看一区二区三区| 最近2019中文字幕mv第一页| 久久亚洲国产成人精品v| 丝袜人妻中文字幕| 国产精品女同一区二区软件| 久久久久久人人人人人| 久久精品人人爽人人爽视色| www.色视频.com| 大片免费播放器 马上看| 国产深夜福利视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 爱豆传媒免费全集在线观看| 精品一区二区三区视频在线| 久久婷婷青草| 国产精品国产三级专区第一集| 国产高清国产精品国产三级| 蜜桃国产av成人99| 亚洲成色77777| 精品久久久久久电影网| 亚洲欧美清纯卡通| 久久久国产欧美日韩av| 成人毛片a级毛片在线播放| 亚洲欧洲日产国产| 国产精品一二三区在线看| 久久这里有精品视频免费| 美女大奶头黄色视频| 亚洲内射少妇av| 黄色视频在线播放观看不卡| 国产不卡av网站在线观看| 日产精品乱码卡一卡2卡三| 久久影院123| 亚洲成人一二三区av| 久久久久久久久久成人| 草草在线视频免费看| 美女内射精品一级片tv| 免费少妇av软件| 国产成人精品无人区| 国产激情久久老熟女| 最近2019中文字幕mv第一页| 久久精品国产亚洲av涩爱| 一区在线观看完整版| 亚洲精品av麻豆狂野| 自线自在国产av| 欧美日韩综合久久久久久| 日日撸夜夜添| 水蜜桃什么品种好| 亚洲久久久国产精品| 99热国产这里只有精品6| 丝袜美足系列| 日本欧美视频一区| 一区二区三区精品91| 青青草视频在线视频观看| 香蕉精品网在线| 五月开心婷婷网| 久久这里只有精品19| 99九九在线精品视频| 国产亚洲精品第一综合不卡 | 一边摸一边做爽爽视频免费| 亚洲国产日韩一区二区| 少妇人妻久久综合中文| av女优亚洲男人天堂| 80岁老熟妇乱子伦牲交| 美女主播在线视频| av电影中文网址| 精品国产一区二区三区四区第35| 2022亚洲国产成人精品| 国产成人免费无遮挡视频| av线在线观看网站| 交换朋友夫妻互换小说| 人人澡人人妻人| 久久久久久人人人人人| 国产成人精品久久久久久| 国产男女内射视频| 国产精品国产三级国产专区5o| 亚洲精品美女久久久久99蜜臀 | 国产男女内射视频| 丝袜在线中文字幕| 你懂的网址亚洲精品在线观看| 一级片'在线观看视频| 水蜜桃什么品种好| 国产一区二区激情短视频 | 黄色毛片三级朝国网站| av在线app专区| 寂寞人妻少妇视频99o| 一级毛片电影观看| 18禁动态无遮挡网站| 波野结衣二区三区在线| 国产精品三级大全| 国产欧美日韩一区二区三区在线| 国产有黄有色有爽视频| 久久ye,这里只有精品| 久热久热在线精品观看| 久久人妻熟女aⅴ| 有码 亚洲区| 韩国高清视频一区二区三区| 欧美日韩综合久久久久久| 国产国拍精品亚洲av在线观看| 少妇人妻精品综合一区二区| 日韩熟女老妇一区二区性免费视频| 色吧在线观看| 三级国产精品片| 在线天堂中文资源库| 我的女老师完整版在线观看| 黄色视频在线播放观看不卡| 女人久久www免费人成看片| 国产精品国产三级国产av玫瑰| av女优亚洲男人天堂| 人人妻人人爽人人添夜夜欢视频| 精品国产一区二区久久| a 毛片基地| 国产成人一区二区在线| 日韩大片免费观看网站| 热re99久久精品国产66热6| 成人免费观看视频高清| 国产精品三级大全| 国产深夜福利视频在线观看| 成人毛片a级毛片在线播放| 欧美日韩视频高清一区二区三区二| 九九爱精品视频在线观看| 免费观看在线日韩| 青青草视频在线视频观看| 国产在线一区二区三区精| 亚洲国产精品一区二区三区在线| 亚洲精品久久午夜乱码| 涩涩av久久男人的天堂| 欧美精品一区二区免费开放| 一本大道久久a久久精品| 日韩在线高清观看一区二区三区| 久久久久久伊人网av| 午夜激情av网站| 久久精品国产亚洲av涩爱| 最近最新中文字幕免费大全7| 51国产日韩欧美| 国产欧美亚洲国产| 国产精品三级大全| 国产乱人偷精品视频| 天天操日日干夜夜撸| 日本爱情动作片www.在线观看| 亚洲av在线观看美女高潮| 一级片免费观看大全| 一区二区三区四区激情视频| 午夜福利乱码中文字幕| 亚洲国产精品专区欧美| 99热这里只有是精品在线观看| 波野结衣二区三区在线| 久久人人爽av亚洲精品天堂| 精品一区二区三卡| 免费播放大片免费观看视频在线观看| 亚洲精品成人av观看孕妇| 黄色一级大片看看| av一本久久久久| 中文字幕免费在线视频6| 男女免费视频国产| 久久国产亚洲av麻豆专区| 成人综合一区亚洲| 草草在线视频免费看| 亚洲欧美一区二区三区国产| 夫妻性生交免费视频一级片| 国产精品女同一区二区软件| 国产淫语在线视频| 亚洲av.av天堂| 99热6这里只有精品| 国产精品一区二区在线观看99| √禁漫天堂资源中文www| 天天躁夜夜躁狠狠躁躁| 一区二区三区精品91| 亚洲综合精品二区| 欧美日本中文国产一区发布| 亚洲欧洲日产国产| 免费黄频网站在线观看国产| 成人手机av| 人妻人人澡人人爽人人| 亚洲,欧美,日韩| 国产精品人妻久久久影院| 热re99久久精品国产66热6| 欧美xxⅹ黑人| av国产精品久久久久影院| 国产免费又黄又爽又色| 少妇人妻 视频| 一本久久精品| 久久狼人影院| 一区二区av电影网| 超碰97精品在线观看| 免费女性裸体啪啪无遮挡网站| 色视频在线一区二区三区| 欧美精品高潮呻吟av久久| 免费观看a级毛片全部| 五月天丁香电影| 国产白丝娇喘喷水9色精品| 亚洲欧洲日产国产| 国产免费又黄又爽又色| 日韩伦理黄色片| 侵犯人妻中文字幕一二三四区| 亚洲av电影在线进入| 青春草亚洲视频在线观看| 99久久人妻综合| 欧美日韩成人在线一区二区| 99久久精品国产国产毛片| 成年人免费黄色播放视频| 亚洲内射少妇av| 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| 国产精品99久久99久久久不卡 | 亚洲国产精品999| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久精品古装| 成人毛片60女人毛片免费| 毛片一级片免费看久久久久| 国产成人午夜福利电影在线观看| 日韩av免费高清视频| 欧美成人午夜免费资源| 欧美bdsm另类| 免费看av在线观看网站| 日韩一本色道免费dvd| 婷婷色综合www| 欧美+日韩+精品| 色哟哟·www| 精品亚洲成国产av| 五月伊人婷婷丁香| 国产精品久久久久久精品古装| 伊人亚洲综合成人网| videossex国产| 国产在线视频一区二区| 亚洲国产精品一区三区| 蜜桃在线观看..| av又黄又爽大尺度在线免费看| 国产爽快片一区二区三区| 蜜臀久久99精品久久宅男| 精品一区二区三区视频在线| 国产免费福利视频在线观看| 欧美人与性动交α欧美软件 | 老司机亚洲免费影院| 最黄视频免费看| 少妇人妻精品综合一区二区| 美女大奶头黄色视频| 亚洲精品国产av成人精品| 国产成人精品福利久久| 久久97久久精品| 亚洲欧美成人精品一区二区| 97在线视频观看| 一个人免费看片子| 色婷婷av一区二区三区视频| 亚洲av在线观看美女高潮| 日本午夜av视频| 少妇熟女欧美另类| 91午夜精品亚洲一区二区三区| 欧美 亚洲 国产 日韩一| 亚洲经典国产精华液单| 热re99久久国产66热| 高清视频免费观看一区二区| 青青草视频在线视频观看| 赤兔流量卡办理| 国产一区亚洲一区在线观看| av在线老鸭窝| 少妇被粗大猛烈的视频| 亚洲国产色片| 免费在线观看完整版高清| 日韩伦理黄色片| av电影中文网址| 久热这里只有精品99| 免费观看a级毛片全部| 一级毛片我不卡| 亚洲国产看品久久| 亚洲天堂av无毛| 午夜91福利影院| 国产日韩欧美在线精品| 一区在线观看完整版| 久久精品国产综合久久久 | 免费不卡的大黄色大毛片视频在线观看| 99久久精品国产国产毛片| 日韩欧美一区视频在线观看| 久久久久精品久久久久真实原创| 欧美日韩av久久| 久久久久久久久久成人| 五月玫瑰六月丁香| 香蕉精品网在线| 国产精品久久久久久久电影| 9色porny在线观看| 精品亚洲成国产av| 国产精品无大码| 毛片一级片免费看久久久久| 少妇的逼水好多| 日本色播在线视频| 一本久久精品| 99国产精品免费福利视频| 极品人妻少妇av视频| 欧美日韩av久久| 国产在线视频一区二区| 一级片'在线观看视频| 妹子高潮喷水视频| 高清黄色对白视频在线免费看| 精品第一国产精品| 韩国精品一区二区三区 | 一级黄片播放器| 一级片免费观看大全| 男人爽女人下面视频在线观看| 日本-黄色视频高清免费观看| 精品国产国语对白av| 国产在线免费精品| 热re99久久国产66热| av有码第一页| 久久国产亚洲av麻豆专区| 免费大片黄手机在线观看| 国产成人av激情在线播放| 欧美老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 成年女人在线观看亚洲视频| 欧美精品国产亚洲| 久久精品aⅴ一区二区三区四区 | 色视频在线一区二区三区| 国产伦理片在线播放av一区| 欧美日韩视频高清一区二区三区二| 中国美白少妇内射xxxbb| 伦理电影免费视频| 18禁国产床啪视频网站| 精品第一国产精品| 捣出白浆h1v1| 久久韩国三级中文字幕| 色网站视频免费| 99九九在线精品视频| 成人亚洲欧美一区二区av| 男男h啪啪无遮挡| 国产一区二区激情短视频 | 妹子高潮喷水视频| 99香蕉大伊视频| 国产成人欧美| 日韩精品有码人妻一区| 亚洲国产精品国产精品| 国产成人免费观看mmmm| 国产精品久久久av美女十八| 高清av免费在线| h视频一区二区三区| 有码 亚洲区| 一区二区日韩欧美中文字幕 | av在线app专区| 日日撸夜夜添| 日韩视频在线欧美| 91久久精品国产一区二区三区| 午夜视频国产福利| 国产精品女同一区二区软件| 久久久久国产网址| 国产免费一区二区三区四区乱码| 狠狠婷婷综合久久久久久88av| 高清欧美精品videossex| 狠狠婷婷综合久久久久久88av| 久久国内精品自在自线图片| 男女下面插进去视频免费观看 | 两性夫妻黄色片 | 男女啪啪激烈高潮av片| 国产乱来视频区| 日产精品乱码卡一卡2卡三| 激情视频va一区二区三区| 日本与韩国留学比较| 另类精品久久| 搡女人真爽免费视频火全软件| 日韩制服丝袜自拍偷拍| 免费不卡的大黄色大毛片视频在线观看| 亚洲第一av免费看| 欧美国产精品一级二级三级| 午夜福利,免费看| 久久青草综合色| 大香蕉久久网| 美国免费a级毛片| 成年美女黄网站色视频大全免费| 国产精品99久久99久久久不卡 | 在线天堂最新版资源| av不卡在线播放| 欧美精品人与动牲交sv欧美| 国产精品一二三区在线看| 一区在线观看完整版| 男人舔女人的私密视频| 日韩大片免费观看网站| 日韩一区二区视频免费看| 国产精品一区二区在线不卡| 国产在线视频一区二区| 成人无遮挡网站| 全区人妻精品视频| 亚洲欧美中文字幕日韩二区| 亚洲精品,欧美精品| 黄色 视频免费看| 亚洲av成人精品一二三区| 精品亚洲成国产av| 多毛熟女@视频| 国产成人午夜福利电影在线观看| 最近中文字幕高清免费大全6| 黄色配什么色好看| 亚洲精品国产av成人精品| 成人毛片60女人毛片免费| 免费黄色在线免费观看| 日韩中字成人| 久久久a久久爽久久v久久| 国产男人的电影天堂91| 熟女av电影| 久久青草综合色| 久久女婷五月综合色啪小说| 国产精品嫩草影院av在线观看| 一级,二级,三级黄色视频| a级毛色黄片| 五月伊人婷婷丁香| 一个人免费看片子| 伦理电影免费视频| 国产成人精品福利久久| 久久国内精品自在自线图片| 这个男人来自地球电影免费观看 | 欧美国产精品va在线观看不卡| 午夜福利视频精品| 成年人免费黄色播放视频| 一区二区三区四区激情视频| 肉色欧美久久久久久久蜜桃| 亚洲欧美一区二区三区国产| 日本午夜av视频| videos熟女内射| 2018国产大陆天天弄谢| 高清视频免费观看一区二区| 国产免费一区二区三区四区乱码| 肉色欧美久久久久久久蜜桃| 秋霞在线观看毛片| 国产一区二区在线观看日韩| 国产高清三级在线| 国产片内射在线| 99re6热这里在线精品视频| 91精品伊人久久大香线蕉| 午夜福利在线观看免费完整高清在| 女性被躁到高潮视频| 91午夜精品亚洲一区二区三区| 999精品在线视频| 国产国语露脸激情在线看| 午夜影院在线不卡| 久久午夜综合久久蜜桃| 久久女婷五月综合色啪小说| 下体分泌物呈黄色| 五月伊人婷婷丁香| 91成人精品电影| 国产激情久久老熟女| 午夜免费男女啪啪视频观看| 久久青草综合色| 一级毛片黄色毛片免费观看视频| 日韩 亚洲 欧美在线| 婷婷色麻豆天堂久久| 精品99又大又爽又粗少妇毛片| 各种免费的搞黄视频| 国产精品人妻久久久久久| 亚洲精品av麻豆狂野| 精品熟女少妇av免费看| 一级黄片播放器| 亚洲欧洲国产日韩| 久久这里有精品视频免费| xxxhd国产人妻xxx| 天堂中文最新版在线下载| 三级国产精品片| 日本wwww免费看| 久久这里有精品视频免费| 看免费av毛片| 一级毛片黄色毛片免费观看视频| 韩国av在线不卡| 又大又黄又爽视频免费| 成年女人在线观看亚洲视频| 亚洲 欧美一区二区三区| 天堂中文最新版在线下载| 黄网站色视频无遮挡免费观看| 精品一区二区三区四区五区乱码 | 成人黄色视频免费在线看| 在线观看人妻少妇| 两个人看的免费小视频| 久久久久久久久久人人人人人人| 少妇 在线观看| 最新中文字幕久久久久| 日韩电影二区| 日本爱情动作片www.在线观看| 免费女性裸体啪啪无遮挡网站| 韩国高清视频一区二区三区| 最近最新中文字幕免费大全7| 蜜桃在线观看..| 王馨瑶露胸无遮挡在线观看| 日韩av在线免费看完整版不卡| 中文字幕免费在线视频6| 久久久国产一区二区| 国产一区二区三区综合在线观看 | 18在线观看网站| 亚洲精品一区蜜桃| 中文字幕最新亚洲高清| 黑丝袜美女国产一区| 一边摸一边做爽爽视频免费| 99热网站在线观看| 宅男免费午夜| 成人无遮挡网站| 亚洲精品久久久久久婷婷小说| 韩国精品一区二区三区 | 看免费成人av毛片| 亚洲精品日本国产第一区| 69精品国产乱码久久久| 一级爰片在线观看| 亚洲 欧美一区二区三区| 男女下面插进去视频免费观看 | 亚洲精品色激情综合| 欧美激情国产日韩精品一区| 大香蕉97超碰在线| 香蕉国产在线看| 男女午夜视频在线观看 | 免费高清在线观看日韩| 亚洲精品第二区| 亚洲国产av影院在线观看| 国产又色又爽无遮挡免| 久久精品国产综合久久久 | 久久精品国产a三级三级三级| 丰满乱子伦码专区| 久久女婷五月综合色啪小说| 韩国av在线不卡| 日韩大片免费观看网站| 国产探花极品一区二区| 亚洲国产色片| 一区二区av电影网| 亚洲精品一二三| 亚洲国产最新在线播放| 大香蕉久久成人网| 国产又爽黄色视频| 久久国产亚洲av麻豆专区| 亚洲激情五月婷婷啪啪| 精品亚洲成a人片在线观看| 天堂8中文在线网| 欧美xxⅹ黑人| 中文字幕免费在线视频6| www.熟女人妻精品国产 | 熟女人妻精品中文字幕| 国产国语露脸激情在线看| 国产成人午夜福利电影在线观看| 一区二区三区乱码不卡18| 日本与韩国留学比较| 国产xxxxx性猛交| 少妇猛男粗大的猛烈进出视频| 婷婷色麻豆天堂久久| tube8黄色片| 国产av国产精品国产| 91午夜精品亚洲一区二区三区| 日日啪夜夜爽| 人人妻人人添人人爽欧美一区卜| 永久免费av网站大全| 亚洲美女搞黄在线观看| 中国国产av一级| 国产成人精品久久久久久| 中文字幕免费在线视频6| 亚洲欧洲精品一区二区精品久久久 | 卡戴珊不雅视频在线播放| 日韩一本色道免费dvd| videos熟女内射| 久久精品久久精品一区二区三区| 永久免费av网站大全| 国产在线一区二区三区精| 欧美日韩视频高清一区二区三区二| 午夜福利网站1000一区二区三区| 日韩制服丝袜自拍偷拍|