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

    不同種植模式下鳳丹根際微生物群落結(jié)構(gòu)和多樣性分析

    2022-05-16 23:40:31劉智張虹羅會婷王仲偉湯詩杰
    江蘇農(nóng)業(yè)學(xué)報 2022年2期
    關(guān)鍵詞:鳳丹高通量測序

    劉智 張虹 羅會婷 王仲偉 湯詩杰

    摘要:為探究不同種植模式下鳳丹根際微生物群落結(jié)構(gòu)和多樣性,采用高通量測序技術(shù)對采自櫸樹(Zelkova serrata)下種植(ZS)、娜塔櫟(Quercus nattallii)下種植(QS)、露天種植(OS)的鳳丹根際以及空白土壤(CK)樣品土壤細(xì)菌和真菌進(jìn)行測序。結(jié)果顯示,與空白土壤相比,種植鳳丹提高了根際微生物群落的豐富度和多樣性,不同種植模式下鳳丹根際微生物群落結(jié)構(gòu)有所差異。3種種植模式下細(xì)菌的優(yōu)勢菌門均為變形菌門(Proteobacteria)、酸桿菌門(Acidobacteria)、放線菌門(Actinobacteria)以及疣微菌門(Verrucomicrobia),但不同種植模式下每個門的相對豐度有所差異。所有樣品中總共檢測出12個真菌門,3種種植模式下真菌的優(yōu)勢菌門均為子囊菌門(Ascomycota)和擔(dān)子菌門(Basidiomycota),不同種植模式下各真菌門的相對豐度存在差異;屬水平上真菌群落結(jié)構(gòu)差異較大,ZS組、QS組和OS組的優(yōu)勢菌屬分別為Exophiala、Thelephoraceae_unclassified、Ascomycota_unclassified。測定3種種植模式下鳳丹的生長指標(biāo),發(fā)現(xiàn)林下種植的鳳丹株高、葉長、葉寬均高于露天種植。林下種植模式的細(xì)菌需氧表型、革蘭氏陰性表型、生物膜形成表型以及致病潛力表型高于單一種植模式。本研究從根際微生物角度分析了不同種植模式對鳳丹的影響,為鳳丹林下種植提供參考,為后續(xù)促生菌篩選、菌肥研制提供理論支撐。

    關(guān)鍵詞:鳳丹;根際微生物;高通量測序;微生物多樣性;微生物群落結(jié)構(gòu)

    中圖分類號:S154.3文獻(xiàn)標(biāo)識碼:A文章編號:1000-4440(2022)02-0502-10

    Analysis on community structure and diversity in rhizosphere microorganisms of Paeonia ostii under different planting patterns

    LIU Zhi ZHANG Hong LUO Hui-ting WANG Zhong-wei TANG Shi-jie

    Abstract:In order to explore the rhizosphere microbial community structure and diversity of Paeonia ostii under different planting patterns, high-throughput sequencing technology was used to sequence the rhizosphere soil samples of Paeonia ostii collected from three different planting patterns, namely Zelkova serrata (ZS), Quercus nattallii (QS) and open-air planting (OS) and blank soil(CK). The results showed that compared with the blank soil, planting Paeonia ostii increased the richness and diversity of rhizosphere microbial community, and the rhizosphere microbial community structure of Paeonia ostii was different under different planting patterns. The dominant phyla of bacteria under three planting patterns were Proteobacteria, Acidobacteria, Actinobacteria and Verrucomicrobia, but the relative abundance of each phylum under different planting patterns was different. A total of 12 fungal phyla were detected in all samples. The dominant phyla of fungi under three planting modes were Ascomycota and Basidiomycota, and the relative abundance of each fungal phylum was different under different planting patterns. The fungal community structure was quite different at genus level, and the dominant genera in ZS group, QS group and OS group were Exophiala, Thelephoraceae_unclassified and Ascomycota_unclassified, respectively. The growth indices of Paeonia ostii under three planting patterns were measured, plant height, leaf length and leaf width of Paeonia ostii planted under forest were higher than those planted in open field. The aerobic phenotype, Gram-negative phenotype, forms biofilms phenotype and potentially pathogenic phenotype in cultivation pattrerns under forest were higher than those in single planting mode. In this study, the effects of different planting patterns on Paeonia ostii were analyzed from the perspective of rhizosphere microorganisms, so as to provide reference for understory planting of Paeonia ostii, and provide theoretical support for subsequent screening of probiotics and development of bacterial fertilizer.

    Key words:Paeonia ostii;rhizosphere microorganisms;high-throughput sequencing;microbial diversity;microbial community structure

    鳳丹(Paeonia ostii)屬芍藥科芍藥屬植物,結(jié)實能力強(qiáng),鳳丹籽產(chǎn)量和出油率高。鳳丹籽油富含多種不飽和脂肪酸,其干燥根皮被稱為丹皮,可入藥,花兼具觀賞、食用和藥用價值,在全國范圍內(nèi)得到了廣泛的種植[1-4]。國家發(fā)展和改革委員會等部門于2020年印發(fā)的《關(guān)于科學(xué)利用林地資源促進(jìn)木本糧油和林下經(jīng)濟(jì)高質(zhì)量發(fā)展的意見》指出,到2025年,要新增或改造木本糧油經(jīng)濟(jì)林3.33×106 hm2,木本食用油年產(chǎn)量達(dá)2.5×106 t,林下經(jīng)濟(jì)年產(chǎn)值達(dá)1×1012 元。作為多年生灌木,鳳丹抗逆性強(qiáng)、適應(yīng)性好,可林下種植,與核桃、香椿、油茶等作物已經(jīng)形成了較為成熟的套種模式。林下種植充分利用了土地資源,同時發(fā)揮了林地的林蔭優(yōu)勢,改善了鳳丹的生長環(huán)境,對提高產(chǎn)量、促進(jìn)增收有著積極意義[5]。

    根際是植物根系與土壤微生物相互作用的場所,具體指植物的根表及受根系直接影響的土壤區(qū)域,是地球上最復(fù)雜的生態(tài)系統(tǒng)之一[6-7]。植物的根系分泌物為根際微生物的生長提供營養(yǎng)物質(zhì),大約有10%的凈光合固定碳通過分泌物的形式釋放到根際,從而形成不同的根際微生物群落[8]。不同的根際微生物可以通過不同的方式,比如調(diào)控植物激素水平、影響植物信號的產(chǎn)生與傳導(dǎo)、改善營養(yǎng)獲取方式以及提高植物抗逆性等來影響植物的生長發(fā)育[9-11]。一些根際微生物可以通過形成生物膜來幫助植物抵抗高溫和高鹽脅迫[12],枯草芽孢桿菌(Bacillus subtilis)可以分泌激素類和抗菌素類物質(zhì)來抵抗病原菌[13]。根際微生物群落結(jié)構(gòu)受到植物類型的影響,楊瀟湘等[14]比較油菜和大豆根際土壤微生物群落結(jié)構(gòu)發(fā)現(xiàn),兩者在門水平上的優(yōu)勢菌群相同,但是后者富含具有生防作用的微生物,比如黃桿菌屬(Flavobacterium)、鞘脂單胞菌屬(Sphingomonas)、芽孢桿菌屬(Bacillus)等;大豆根際具有更多的病原細(xì)菌和真菌,比如腸桿菌(Enterobacter)和炭疽菌(Colletotrichum spp.)等。除此之外,土壤環(huán)境、栽培模式、水肥條件、田間管理模式等也會影響植物根際微生物結(jié)構(gòu)。例如劉振香等[15]研究發(fā)現(xiàn),夏玉米田中施用有機(jī)肥可以提高土壤微生物多樣性和豐富度;姜莉莉等[16]研究發(fā)現(xiàn),蘋果園內(nèi)種植黑麥草可提高土壤中鏈格孢菌(Alternaria)和曲霉菌(Aspergillus)的相對豐度,種植白三葉可提高土壤赤霉菌(Gibberella)的相對豐度,而自然生草則會提高青枯菌(Ralstonia)等致病菌的相對豐度。

    相比于傳統(tǒng)土壤微生物研究方法,高通量測序技術(shù)所需樣本量少、能用環(huán)境樣本中直接獲取的總DNA進(jìn)行文庫的構(gòu)建,操作簡單、準(zhǔn)確度高,成為研究土壤微生物的熱門手段[17-19]。前人已經(jīng)研究了不同地區(qū)[20]、不同種植年限[21-22]鳳丹根際微生物多樣性,但對于不同種植模式下鳳丹根際微生物多樣性的研究還未見報道。本研究基于高通量測序技術(shù),研究不同種植模式下鳳丹根際微生物群落結(jié)構(gòu)和多樣性,以期為鳳丹根際微生態(tài)環(huán)境研究提供理論基礎(chǔ),為優(yōu)化油用牡丹栽培技術(shù)提供思路。

    1材料與方法

    1.1試驗材料

    試驗地位于江蘇省中國科學(xué)院植物研究所試驗田(南京市六合區(qū)竹鎮(zhèn)鎮(zhèn),32°33′46″N,118°39′28″E),土壤類型均為黃棕壤,屬亞熱帶季風(fēng)氣候。試驗田內(nèi)設(shè)有櫸樹下種植(ZS)、娜塔櫟下種植(QS)以及露天種植(OS)3種模式,櫸樹和娜塔櫟的株距×行距為4 m×6 m。3種模式下種植的鳳丹均為3年生苗,株距×行距為40 cm×60 cm,種植密度約為1 hm2 42 000株,水肥管理均參照國家林業(yè)局《油用牡丹栽培技術(shù)規(guī)程》[23]進(jìn)行。于2021年8月采樣,每個試驗區(qū)均采用五點采樣法,挖出整株鳳丹植株,采用抖根法采集根際。按照相同方法,在空地上采集地表10~15 cm深度的土壤作為對照組(CK),每組樣品均設(shè)6個生物學(xué)重復(fù)。收集的樣品裝入無菌自封袋,冰盒保存,帶回實驗室備用。土壤理化性質(zhì)見表1。

    1.2鳳丹生長指標(biāo)測定

    生長指標(biāo)于2021年6月底測定,每個小區(qū)隨機(jī)選取10株鳳丹植株,使用卷尺(精度0.1 cm)測量株高,每1株上選取3片健康成熟頂生小葉,使用直尺(精度0.1 cm)測量葉長(葉片基部至頂端長度)和葉寬(葉片最寬處寬度),重復(fù)3次,求平均值。

    1.3根際微生物基因組DNA的提取、擴(kuò)增以及高通量測序

    使用天根DP712試劑盒(天根生化科技有限公司產(chǎn)品)提取根際微生物基因組DNA,采用1%瓊脂糖凝膠電泳檢測其完整性,采用分光光度計檢測其濃度和純度。使用引物F:5′-CCTACGGGNGGCWGCAG-3′,R:5′-GACTACHVGGGTATCTAATCC-3′擴(kuò)增細(xì)菌的16 S rDNA的V3+V4區(qū);使用引物F:5′-GTGARTCATCGAATCTTTG-3′,R:5′-TCCTCCGCTTATTGATATGC-3′擴(kuò)增真菌的ITS2區(qū)域。PCR反應(yīng)體系為:Phusion Hot Start Flex 2×Master Mix 12.5 μl,1 μmol/L的正反向引物各2.5 μl,DNA模板50 ng,加ddH2O至25.0 μl。PCR反應(yīng)條件為:98 ℃ 30 s;98 ℃ 10 s,54 ℃ 30 s,72 ℃ 45 s,35個循環(huán);72 ℃ 10 min。

    PCR產(chǎn)物經(jīng)純化和定量之后構(gòu)建DNA文庫,評估合格的測序文庫(由聯(lián)川生物科技有限公司提供)在Illumina NovaSeq平臺上進(jìn)行測序。

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

    使用軟件Cutadapt(v1.9)去除barcode標(biāo)簽和引物序列,使用Fqtrim(v0.94)和Vsearch(v2.3.4)分別過濾低質(zhì)量序列和嵌合體,利用DADA2進(jìn)行解調(diào),得到特征(feature)序列和特征豐度表格。

    對特征豐度進(jìn)行歸一化后,使用QIIME2計算樣品的α多樣性指數(shù)(Chao1指數(shù)、Observed species指數(shù)、Goods coverage指數(shù)、Shannon指數(shù)、Simpson指數(shù))和β多樣性指數(shù)。使用SILVA數(shù)據(jù)庫(https://www.arb-silva.de/documentation/release138/)注釋特征序列,獲取分類學(xué)信息,使用Excel 2016和SPSS 26處理數(shù)據(jù)。將樣本信息上傳至BugBase(https://bugbase.cs.umn.edu/),預(yù)測細(xì)菌表型。使用R(VennDiagram)繪制特征分布Venn圖,使用R(vegen 2.5.4)繪制PCA圖,使用Origin Pro 2021繪制箱線圖和群落結(jié)構(gòu)圖。

    2結(jié)果與分析

    2.1測序結(jié)果質(zhì)量評估

    原始數(shù)據(jù)處理后,統(tǒng)計高質(zhì)量數(shù)據(jù)。16 S rDNA測序序列個數(shù)為1 717 754,有99.88%的序列長度在400~500 bp,有效數(shù)據(jù)與原始數(shù)據(jù)之比(Valid)均值為85.31%,其中最大值為94.01%,最小值為73.87%;所有樣本的Q30均在80.00%以上,平均值為91.76%。

    ITS2測序序列個數(shù)為1 969 418,有70.10%的序列長度在200~300 bp,29.18%的序列長度在300~400 bp,Vaild最小值為94.26%,Q30均在95%以上。說明,此次測序質(zhì)量良好。

    2.2不同種植模式下鳳丹根際微生物豐富度和多樣性

    與傳統(tǒng)操作分類單元(OTU)相比,DADA2算法可以對擴(kuò)增子數(shù)據(jù)進(jìn)行去噪、修正、去冗余,得到特征數(shù)據(jù),數(shù)據(jù)精確度和物種分辨率大大提升。由表2可見,在16 S rDNA和ITS2測序中,ZS組的特征數(shù)最多,QS組和OS組次之,CK的特征數(shù)最少。

    α多樣性指數(shù)可以用來反映樣本的物種豐富度和多樣性。在16 S rDNA測序中,4組的Observed species指數(shù)之間無顯著差異;4組的Simpson指數(shù)均為1.00。在ITS2測序中,ZS、QS、OS組的Observed species指數(shù)均高于CK; 各組間Simpson指數(shù)和Shannon指數(shù)差異不顯著。

    從圖1可以直觀地看出,每組樣品共有和特有的特征數(shù)目。在16 S rDNA測序中,4組樣品共有的特征數(shù)為839,ZS組、QS組、OS組、CK特有的特征數(shù)分別為:4 631、4 635、5 316和5 239;ITS2測序中,4組樣品共有的特征數(shù)為124,ZS組、QS組、OS組、CK特有的特征數(shù)分別為670、562、582和378。

    從以上結(jié)果可以看出,相比于空白土壤,3種種植模式均提高了鳳丹根際微生物的豐富度和多樣性;3種不同的種植模式中,鳳丹根際細(xì)菌和真菌的豐富度和多樣性也存在差異。

    2.3不同種植模式下鳳丹根際細(xì)菌群落結(jié)構(gòu)

    取鳳丹根際細(xì)菌門水平相對豐度前10位(TOP10)繪制柱狀堆疊圖(圖2A)。在門水平上細(xì)菌群落結(jié)構(gòu)中,變形菌門(Proteobacteria)在4組中豐度最高,ZS組、QS組、OS組均高于CK;其次是酸桿菌門(Proteobacteria),在ZS組和CK中的豐度相當(dāng),均高于QS組和OS組;OS組中放線菌門(Actinobacteria)的豐度大于其他3組;疣微菌門(Verrucomicrobia)在ZS組和QS組中的豐度高于OS組和CK;擬桿菌門(Bacteroidetes)在OS組中的豐度要高于其他3組;綠彎菌門(Chloroflexi)、己科河菌門(Rokubacteria)以及硝化螺旋菌門(Nitrospirae)在CK中的豐度分別為12.9%、4.6%和3.0%,均高于在其他3組中的豐度。

    圖2B是屬水平上鳳丹根際細(xì)菌群落結(jié)構(gòu),屬水平上優(yōu)勢菌多為現(xiàn)階段未分類的細(xì)菌,在ZS組和QS組中,Candidatus_Udaeobacter的相對豐度均為最高,分別為7.80%和6.07%,其次為Subgroup_6_unclassified,相對豐度分別為6.1%和4.8%, 排在第三位的是Rokubacteriales_unclassified,相對豐度分別為3.0%和3.1%。在OS組中Subgroup_6_unclassified的相對豐度最高,之后按照相對豐度高低依次是Rokubacteriales_unclassified、KD4-96_unclassified、Actinobacteria_unclassified。在CK中,按照豐度高低依次是:Rokubacteriales_unclassified、Acidobacteriales_unclassified、Subgroup_6_unclassified、Candidatus_Udaeobacter。

    由此可見,不同的種植模式對鳳丹根際細(xì)菌群落結(jié)構(gòu)是有一定影響的,特別是酸桿菌門、放線菌門和疣微菌門。在屬水平上,Subgroup_6_unclassified、Candidatus_Udaeobacter和Rokubacteriales_unclassified在不同的種植模式下差異較大。

    2.4不同種植模式下鳳丹根際真菌群落結(jié)構(gòu)

    由圖3A可以看出,在ZS組、OS組和CK中子囊菌門(Ascomycota)相對豐度最高,擔(dān)子菌門(Basidiomycota)位居第二,而在QS組中擔(dān)子菌門豐度最高。OS組中Fungi_unclassified、球囊菌門(Glomeromycota)以及接合菌門(Zygomycota)的相對豐度要高于其他3組。毛霉門(Mucoromycota)和被孢霉門(Mortierellomycota)沒有出現(xiàn)在CK中。

    由圖3B可知,ZS組中相對豐度TOP3依次是外瓶霉屬(Exophiala)、叢赤殼屬(Nectria)、Ascomycota_unclassified,其中外瓶霉屬的相對豐度達(dá)到了22.5%;QS組中Thelephoraceae_unclassified的相對豐度最高,為21.3%,其次為叢赤殼屬;OS組中,Ascomycota_unclassified的豐度最高,其次為Fungi_unclassified;而在CK中,相對豐度TOP3為小菇屬(Mycena)(16.6%)、Chaetothyriales_unclassified(12.2%)、鐮刀菌屬(Fusarium)(9.0%)。

    不同種植模式下,鳳丹根際真菌群落結(jié)構(gòu)有著較大的差異。門水平上子囊菌門、擔(dān)子菌門、Fungi_unclassified、接合菌門豐度差異較大。而在屬水平上差異更明顯,每組的優(yōu)勢屬均不同。

    2.5不同種植模式下鳳丹根際微生物β多樣性分析

    對不同種植模式下鳳丹根際細(xì)菌和真菌群落結(jié)構(gòu)進(jìn)行主成分分析,結(jié)果見圖4,根據(jù)樣品代表的點在PCA圖中的距離,可以判斷樣品之間的相似度。同時PCA圖給出了基于Anosim相似性分析計算出的R值,R值越接近于1,說明組間差異越大于組內(nèi)差異。

    從圖4A可以看出,在不同種植模式下鳳丹根際細(xì)菌群落PCA圖中,前2個主成分的累計貢獻(xiàn)率達(dá)到了48.58%。ZS組和QS組樣品對應(yīng)的點在PCA圖中較為集中,且兩組距離較近,說明ZS組和QS組的細(xì)菌群落組成較為相似。OS組中的大部分點與ZS組、QS組分離,CK與ZS組、QS組明顯分離。結(jié)合R=0.729 5可知,不同組別之間差異較大,說明相比于CK,種植鳳丹改變了根際細(xì)菌的群落結(jié)構(gòu),2種林下種植模式的細(xì)菌群落組成較為相似,且與露天種植模式存在差異。

    圖4B是不同種植模式下鳳丹根際真菌群落PCA圖,ZS組、QS組和OS組樣品對應(yīng)的點在PCA圖中距離較近,說明3組的真菌群落結(jié)構(gòu)較為相似,CK與其他3組有一定的分離,結(jié)合R=0.652 9可知,相比于CK,種植鳳丹改變了根際真菌群落結(jié)構(gòu)。

    2.6不同種植模式下鳳丹的生長狀況

    不同種植模式下鳳丹的生長狀況如表3所示,ZS組的株高高于QS組和OS組,其中和OS組差異顯著;ZS組和QS組的葉長均顯著高于OS組;ZS組和QS組的葉寬均顯著高于OS組。

    2.7不同種植模式下鳳丹根際細(xì)菌表型預(yù)測

    細(xì)菌表型預(yù)測結(jié)果(圖5)表明,不同種植模式下鳳丹根際細(xì)菌表型之間存在較大差異。ZS組和QS組的需氧表型顯著高于OS組;可移動元件表型在OS組中最高,顯著高于ZS組和CK,但與QS組沒有顯著差異;兼性厭氧表型表現(xiàn)為ZS組、OS組高于CK, QS組和CK沒有顯著差異;ZS組、QS組的生物膜形成表型要高于CK, OS組與CK無顯著差異;革蘭氏陰性表型在ZS組、QS組和CK中無顯著差異,但ZS組顯著高于OS組;革蘭氏陽性表型則表現(xiàn)為OS組顯著高于ZS組;ZS組和QS組的致病潛力顯著高于CK;ZS組的脅迫耐受低于OS組和CK。厭氧表型在4組中沒有顯著差異。

    3討論

    研究結(jié)果表明,選取合理的種植模式能增強(qiáng)作物根際的交互作用,優(yōu)化根際微生物群落結(jié)構(gòu),對提高土壤肥力、促進(jìn)作物生長發(fā)育有著積極作用[24]。輪作、套作以及連作模式下,丹參根際微生物群落結(jié)構(gòu)存在一定差異,輪作和套作對改善土壤質(zhì)量以及緩解丹參連作障礙有著一定的積極作用[25]。香蕉間作韭菜、假花生和覆蓋稻草處理,顯著提高了香蕉根際放線菌的相對豐度,但是尖孢鐮刀菌(F. oxysporum)的豐度有所下降,降低了香蕉枯萎病的發(fā)生幾率[26]。探究不同種植模式下草莓根際細(xì)菌和真菌群落結(jié)構(gòu),發(fā)現(xiàn)限根種植提高了草莓根際放線菌門、變形菌門、擬桿菌門以及子囊菌門的豐度,降低了綠彎菌門、酸桿菌門和接合菌門的豐度[27]。研究發(fā)現(xiàn),鳳丹與香椿、女貞間作提高了單株地上部分生物量和產(chǎn)量[28],鳳丹-香椿套種模式在定植4年后產(chǎn)量高出單一種植模式51.97%,極大提高了籽粒產(chǎn)量[29]。

    與CK相比,綠彎菌門、己科河菌門以及硝化螺旋菌門在其他處理組中的相對豐度有所降低。綠彎菌多為細(xì)胞絲狀,營養(yǎng)方式多樣,參與CO2的同化過程[30];硝化螺旋菌門可將銨態(tài)氮氧化為硝態(tài)氮,參與土壤氮循環(huán)[31]。3種種植模式下細(xì)菌的優(yōu)勢菌門均為放線菌門,這與張永敢[32]等的研究結(jié)果一致,但是相比于露天種植,林下種植鳳丹根際細(xì)菌酸桿菌門和疣微菌門的相對豐度得到了提高,芽單胞菌門、放線菌門以及擬桿菌門有所下降。相關(guān)研究結(jié)果表明,擬桿菌門在輪作休耕模式中與土壤全鉀含量呈負(fù)相關(guān)[33];而在黃花蒿連作中,隨著種植年限的增加,根際芽單胞菌門和放線菌門的豐度有所降低,說明黃花蒿能選擇性抑制它們的生長繁殖[34]。2種林下種植模式鳳丹根際細(xì)菌群落結(jié)構(gòu)和多樣性也存在一定差異,相比在娜塔櫟下種植鳳丹,在櫸樹下種植鳳丹提高了根際酸桿菌門的豐度,酸桿菌門細(xì)菌能夠在缺氧酸性環(huán)境下生存,抗逆性強(qiáng),對林地中樹葉等植物殘體具有降解功能[35],本研究中ZS組的pH低至3.65,也佐證了這一結(jié)論。

    在所有樣品中,鑒定出的真菌分別屬于12個門,3種種植模式下真菌優(yōu)勢菌門均為子囊菌門和擔(dān)子菌門,且在ZS組中子囊菌門的豐度最高。子囊菌門真菌多為腐生真菌,能分解土壤中的化合物,對土壤養(yǎng)分循環(huán)有著重要的作用,一定程度上可以促進(jìn)作物的生長發(fā)育[36-38];擔(dān)子菌門的一些真菌可以與作物共生形成菌根,增強(qiáng)作物抗性[39]。屬水平上,ZS組叢赤殼屬豐度高于其他3組,被孢霉屬在QS組中豐度最高。王麗娟等[40]從鹿蹄草中分離出1株促生菌,鑒定為叢赤殼屬真菌,該菌株可以顯著提高鹿蹄草的根系活力,提高葉片中葉綠素a的含量、可溶性糖含量以及超氧化物歧化酶活性,從而促進(jìn)鹿蹄草的生長發(fā)育。被孢霉屬真菌主要腐生于土壤,部分種類可以和植物共生,對改善土壤結(jié)構(gòu)有一定作用,能促進(jìn)植物的生長[41]。鐮刀菌屬(Fusarium)和Ceratobasidium作為常見的植物病原菌,對種子的萌發(fā)及植物根莖的生長不利,可能導(dǎo)致作物發(fā)生枯萎病、根腐病等[42-43],而本研究發(fā)現(xiàn)在ZS組中這2種真菌的豐度要低于其他3組,這說明櫸樹下種植鳳丹對抵抗病原菌有一定的作用。

    BugBase是基于測序數(shù)據(jù)進(jìn)行的預(yù)測分析,并不能真實反映細(xì)菌群落的全部表型,但是分析不同種植模式下鳳丹根際的細(xì)菌表型,對于鳳丹的種植仍有一定的理論參考價值。細(xì)菌表型預(yù)測結(jié)果顯示,兼性厭氧表型在CK中最低,這可能是因為隨著鳳丹種植年限增加,連作土壤總孔度減少,從而影響通氣性所導(dǎo)致的[44]。林下種植模式下需氧表型、革蘭氏陰性表型、生物膜形成表型以及致病潛力表型在一定程度上高于單一種植模式,這可能跟變形菌門的豐度升高有關(guān)。變形菌門細(xì)菌在自然界中廣泛存在,革蘭氏陰性細(xì)菌,大多數(shù)都具有脂多糖組成的外膜,其中β-變形菌多為好氧及兼性細(xì)菌,包含很多病原菌,這與楊睿[45]研究貴州喀斯特山區(qū)花椒林根際細(xì)菌表型的結(jié)果相符,隨著變形菌門豐度的升高,具有上述表型的細(xì)菌豐度也隨著升高。研究結(jié)果表明,放線菌門有抗脅迫的能力[46],本研究中QS組和OS組中放線菌門的相對豐度高于其他2組,脅迫耐受表型也表現(xiàn)出相同的變化。細(xì)菌的表型受多種因素影響,在鳳丹的種植中,不同模式下土壤理化性質(zhì)存在差異,環(huán)境異質(zhì)性大,微生物通過表型的適應(yīng)性調(diào)整,以應(yīng)對環(huán)境的改變[47]。

    對果實生長期鳳丹的部分生長指標(biāo)進(jìn)行測定,發(fā)現(xiàn)套種在櫸樹下的鳳丹植株生長狀態(tài)最佳,其株高顯著高于其他2組,這可能和子囊菌門、酸桿菌門中一些有益菌的豐度增加有關(guān)。

    綜上,林下種植鳳丹提高了根際微生物多樣性,不同種植模式下鳳丹根際微生物群落結(jié)構(gòu)存在較大差異,櫸樹下套種鳳丹提高了酸桿菌門、子囊菌門、叢赤殼屬等有益菌的豐度,減少了鐮刀菌屬和Ceratobasidium等致病菌的豐度。測量生長指標(biāo)發(fā)現(xiàn)櫸樹下種植的鳳丹長勢最好,由此推斷,櫸樹下種植鳳丹有助于植株的生長發(fā)育。我們今后還將從土壤代謝組學(xué)、籽粒產(chǎn)量以及鳳丹藥用成分等角度繼續(xù)探索不同種植模式對鳳丹的影響,以期為鳳丹高效栽培技術(shù)以及林農(nóng)復(fù)合經(jīng)營提供新思路。

    參考文獻(xiàn):

    [1]張姍姍,趙凡,魏小豹,等. 鳳丹和紫斑牡丹6個產(chǎn)地種子脂肪酸組分的比較[J]. 中國糧油學(xué)報, 2021, 36(3): 84-90.

    [2]謝一青,吳慶錐,余孟楊,等. 鳳丹南移福建的適應(yīng)性分析[J]. 熱帶作物學(xué)報, 2021,42(2): 592-598.

    [3]王仲偉,王歡利,張文獻(xiàn),等. 江蘇油用牡丹的發(fā)展思考[J]. 江蘇林業(yè)科技, 2016, 43(6): 53-55.

    [4]裴姿琛,胡永宏,劉澤,等. 新型木本油料作物——油用牡丹的開發(fā)和利用[J]. 林業(yè)與生態(tài)科學(xué), 2018, 33(4): 358-363.

    [5]王勇,何舒,熊冰杰,等. 不同栽培模式對人參根際土壤微生物多樣性的影響研究[J]. 中草藥, 2021, 52(17): 5303-5310.

    [6]劉王鎖,李海泉,何毅,等. 根際微生物對植物與土壤交互調(diào)控的研究進(jìn)展[J]. 中國土壤與肥料, 2020(5): 318-327.

    [7]KHAN N, ALI S, SHAHID M , et al. Insights into the interactions among roots, rhizosphere, and rhizobacteria for improving plant growth and tolerance to abiotic stresses: a review[J]. Cells,2021,10(6): 1551

    [8]SONG Y, LI X, YAO S, et al. Correlations between soil metabolomics and bacterial community structures in the pepper rhizosphere under plastic greenhouse cultivation[J]. Science of the Total Environment, 2020, 728: 138439.

    [9]KANG S, KHAN A L, WAQAS M, et al. Plant growth-promoting rhizobacteria reduce adverse effects of salinity and osmotic stress by regulating phytohormones and antioxidants in Cucumis sativus[J]. Journal of Plant Interactions, 2014, 9(1): 673-682.

    [10]YUAN J, ZHAO J, WEN T, et al. Root exudates drive the soil-borne legacy of aboveground pathogen infection[J]. Biomed Central, 2018, 6(1): 156.

    [11]BAL H, NAYAK L,DAS S, et al. Isolation of ACC deaminase producing PGPR from rice rhizosphere and evaluating their plant growth promoting activity under salt stress[J]. Plant Soil,2013, 366: 93-105.

    [12]KIM Y C, LEVEAU J, MCSPADDEN GARDENER B B, et al. The multifactorial basis for plant health promotion by plant-associated bacteria[J]. Applied and Environmental Microbiology, 2011, 77(5): 1548-1555.

    [13]KASIM W A, GAAFAR R M, ABOU-ALI R M, et al. Effect of biofilm forming plant growth promoting rhizobacteria on salinity tolerance in barley[J]. Annals of Agricultural Sciences, 2016,61(2):217-227.

    [14]楊瀟湘,張蕾,黃小琴,等. 基于高通量測序分析大豆和油菜根際微生物群落結(jié)構(gòu)的差異[J]. 應(yīng)用生態(tài)學(xué)報, 2019, 30(7): 2345-2351.

    [15]劉振香,劉鵬,賈緒存,等. 不同水肥處理對夏玉米田土壤微生物特性的影響[J]. 應(yīng)用生態(tài)學(xué)報, 2015, 26(1): 113-121.

    [16]姜莉莉,宮慶濤,武海斌,等. 不同生草處理對蘋果園土壤微生物群落的影響[J]. 應(yīng)用生態(tài)學(xué)報, 2019, 30(10): 3482-3490.

    [17]艾鑠,張麗杰,肖芃穎,等. 高通量測序技術(shù)在環(huán)境微生物領(lǐng)域的應(yīng)用與進(jìn)展[J]. 重慶理工大學(xué)學(xué)報(自然科學(xué)), 2018, 32(9): 111-121.

    [18]王林,李冰,朱健. 高通量測序技術(shù)在人工濕地微生物多樣性研究中的研究進(jìn)展[J]. 中國農(nóng)學(xué)通報, 2016, 32(5): 10-15.

    [19]牛世全,龍洋,李海云,等. 應(yīng)用IlluminaMiSeq高通量測序技術(shù)分析河西走廊地區(qū)鹽堿土壤微生物多樣性[J]. 微生物學(xué)通報, 2017, 44(9): 2067-2078.

    [20]尹原森,馬國勝,曹春燕,等. 不同地區(qū)鳳丹根際土壤微生物功能多樣性分析[J]. 分子植物育種, 2021, 19(20):6918-6926.

    [21]郭麗麗,尹偉倫,郭大龍,等. 油用鳳丹牡丹不同種植時間根際細(xì)菌群落多樣性變化[J]. 林業(yè)科學(xué), 2017, 53(11): 131-141.

    [22]李昱瑩,劉曙光,廉小芳,等. 油用牡丹鳳丹不同種植年限根際真菌群落多樣性變化研究[J]. 基因組學(xué)與應(yīng)用生物學(xué), 2020, 39(4): 1672-1685.

    [23]全國花卉標(biāo)準(zhǔn)化技術(shù)委員會. 油用牡丹栽培技術(shù)規(guī)程: LY/T 2958-2018[S]. 北京: 國家林業(yè)局,2018.

    [24]章家恩,高愛霞,徐華勤,等. 玉米/花生間作對土壤微生物和土壤養(yǎng)分狀況的影響[J]. 應(yīng)用生態(tài)學(xué)報, 2009, 20(7): 1597-1602.

    [25]王悅,楊貝貝,王浩,等. 不同種植模式下丹參根際土壤微生物群落結(jié)構(gòu)變化[J]. 生態(tài)學(xué)報, 2019, 39(13): 4832-4843.

    [26]王麗霞,井濤,殷曉敏,等. 不同種植模式對香蕉根區(qū)土壤養(yǎng)分及可培養(yǎng)微生物數(shù)量的影響[J]. 中國南方果樹, 2020, 49(4): 80-86.

    [27]陳澤斌,陳敏,阮亞男,等. 不同種植模式草莓根際細(xì)菌和真菌群落結(jié)構(gòu)的變化[J]. 昆明學(xué)院學(xué)報, 2019, 41(6): 52-56.

    [28]李娟. 油用牡丹與經(jīng)濟(jì)林木立體栽培模式研究[D]. 鄭州: 鄭州大學(xué), 2019.

    [29]王國霞,羅青,楊玉珍,等. 立體栽培模式對油用牡丹生長和產(chǎn)量的影響[J]. 貴州農(nóng)業(yè)科學(xué), 2018, 46(12): 100-103.

    [30]鮮文東,張瀟橦,李文均. 綠彎菌的研究現(xiàn)狀及展望[J]. 微生物學(xué)報, 2020,60(9):1801-1820.

    [31]劉芯竹. 覆蓋經(jīng)營對毛竹林土壤理化性質(zhì)和細(xì)菌群落多樣性影響[D]. 杭州:浙江農(nóng)林大學(xué), 2021.

    [32]張永敢,趙娟,張玉潔,等. 藥用植物鳳丹(Paeonia suffruticosa)根際土壤細(xì)菌群落16S rRNA基因的ARDRA分析[J]. 生態(tài)學(xué)報, 2016, 36(17): 5564-5574.

    [33]南麗麗,郭全恩,譚杰輝,等. 輪作休耕模式對土壤細(xì)菌群落的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2020, 38(6): 128-134.

    [34]李倩,袁玲,楊水平,等. 連作對黃花蒿生長及土壤細(xì)菌群落結(jié)構(gòu)的影響[J]. 中國中藥雜志, 2016, 41(10): 1803-1810.

    [35]楊安娜,陸云峰,張俊紅,等. 杉木人工林土壤養(yǎng)分及酸桿菌群落結(jié)構(gòu)變化[J]. 林業(yè)科學(xué), 2019, 55(1): 119-127.

    [36]劉婷,肖仲久,李小霞,等. 高通量測序技術(shù)分析茅臺釀酒高粱根際真菌群落特征[J]. 內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2021,43(1): 48-52.

    [37]馬繼瓊,孫一丁,楊奕,等.不同葉稻瘟抗性水稻品種對根際真菌的影響[J]. 江蘇農(nóng)業(yè)科學(xué),2021,49(12):75-80.

    [38]劉震,徐玉鵬,王秀領(lǐng),等.黑龍港苜蓿根際土壤真菌群落結(jié)構(gòu)及多樣性分析[J]. 江蘇農(nóng)業(yè)科學(xué),2021,49(10):197-201.

    [39]薛曉敏,王來平,韓雪平,等. 不同樹盤覆蓋對矮砧蘋果園土壤微生物群落結(jié)構(gòu)和多樣性的影響[J]. 生態(tài)學(xué)報, 2021,41(4): 1528-1536.

    [40]王麗娟,王威威,呂雪,等. 鹿蹄草內(nèi)生真菌的分離鑒定及其對鹿蹄草的生理效應(yīng)[J]. 東北農(nóng)業(yè)大學(xué)學(xué)報, 2013, 44(8): 87-93.

    [41]MICHAL J J, ANATOLI L, ALEXANDRA F, et al. The beneficial root-colonizing fungus Mortierella hyalina promotes the aerial growth of Arabidopsis and activates calcium-dependent responses which restrict Alternaria brassicae-induced disease development in roots.[J]. Molecular Plant-microbe Interactions, MPMI, 2018, 32(3): 351-363.

    [42]張向民. 鐮刀菌屬分類學(xué)研究歷史與現(xiàn)狀[J]. 菌物研究, 2005(2): 63-66.

    [43]唐濤,王帆帆,郭杰,等. 12種生物源殺菌劑對白術(shù)根腐病的防效[J]. 植物保護(hù), 2021, 47(3): 288-293.

    [44]耿貴,楊瑞瑞,於麗華,等. 作物連作障礙研究進(jìn)展[J]. 中國農(nóng)學(xué)通報, 2019, 35(10): 36-42.

    [45]楊睿. 貴州喀斯特山區(qū)花椒林根際土壤微生物群落結(jié)構(gòu)及功能研究[D]. 貴州:貴州師范大學(xué), 2021.

    [46]楊睿,李娟,龍健,等. 貴州喀斯特山區(qū)不同種植年限花椒根際土壤細(xì)菌群落結(jié)構(gòu)特征研究[J]. 生態(tài)環(huán)境學(xué)報, 2021, 30(1): 81-91.

    [47]趙鵬宇,白雪,燕平梅,等. 華北落葉松林土壤細(xì)菌群落結(jié)構(gòu)與表型的環(huán)境異質(zhì)性響應(yīng)[J]. 林業(yè)科學(xué), 2021, 57(7): 101-110.

    (責(zé)任編輯:陳海霞)

    收稿日期:2021-12-16

    基金項目:江蘇省農(nóng)業(yè)科技自主創(chuàng)新基金項目[CX(19)3066] ;江蘇現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)項目[JATS(2021)507]

    作者簡介:劉智(1996-),男,陜西漢中人,碩士研究生,主要從事油用牡丹根際微生物研究。(E-mail)liuzhi2021702@163.com

    通訊作者:王仲偉,(E-mail)w0414115@Sina.cn;湯詩杰,(E-mail)tangshijie69@aliyun.com

    猜你喜歡
    鳳丹高通量測序
    鳳丹形態(tài)及其生理特性的季節(jié)動態(tài)研究
    不同種源地鳳丹后代種子在蘇州萌發(fā)試驗研究
    不同間作模式對田間小氣候特征及鳳丹光合特性和種實性狀的影響
    鳳丹籽油對小鼠H22腫瘤的抑制作用
    鳳丹愈傷組織中丹皮酚含量的測定
    環(huán)狀RNA在疾病發(fā)生中的作用
    川明參輪作對煙地土壤微生物群落結(jié)構(gòu)的影響
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    人參根際真菌群落多樣性及組成的變化
    LncRNAs作為miRNA的靶模擬物調(diào)節(jié)miRNA
    免费在线观看日本一区| 天堂动漫精品| 99久久久亚洲精品蜜臀av| 欧美区成人在线视频| 欧美日韩瑟瑟在线播放| 国产精品乱码一区二三区的特点| 变态另类成人亚洲欧美熟女| 小说图片视频综合网站| 老熟妇乱子伦视频在线观看| 国产高清三级在线| 国产v大片淫在线免费观看| 欧美日韩一级在线毛片| 毛片女人毛片| 午夜福利视频1000在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美一级毛片孕妇| 国产免费男女视频| 无遮挡黄片免费观看| 久久久久久久久久黄片| 国产精品精品国产色婷婷| 深爱激情五月婷婷| 桃红色精品国产亚洲av| 国产在视频线在精品| 成人国产综合亚洲| 亚洲精品日韩av片在线观看 | 热99re8久久精品国产| 国产在线精品亚洲第一网站| 丰满的人妻完整版| 亚洲欧美日韩卡通动漫| 亚洲久久久久久中文字幕| 国产免费男女视频| 日韩 欧美 亚洲 中文字幕| 国产精品自产拍在线观看55亚洲| e午夜精品久久久久久久| 国产激情偷乱视频一区二区| 在线观看av片永久免费下载| 十八禁网站免费在线| 午夜免费激情av| 蜜桃久久精品国产亚洲av| 在线观看舔阴道视频| 国产精品综合久久久久久久免费| 一二三四社区在线视频社区8| 国产成人影院久久av| 亚洲精品乱码久久久v下载方式 | 国产真实乱freesex| 黄色丝袜av网址大全| 欧美激情久久久久久爽电影| 久久这里只有精品中国| 成人无遮挡网站| 免费看日本二区| 亚洲无线观看免费| 国产野战对白在线观看| 搡老妇女老女人老熟妇| 亚洲avbb在线观看| 小说图片视频综合网站| 波多野结衣高清作品| 中出人妻视频一区二区| 欧美绝顶高潮抽搐喷水| 精品国产三级普通话版| 十八禁人妻一区二区| 国内久久婷婷六月综合欲色啪| 亚洲av五月六月丁香网| 丰满的人妻完整版| 国产欧美日韩一区二区精品| 国产亚洲精品综合一区在线观看| 亚洲国产精品999在线| 香蕉丝袜av| xxxwww97欧美| 久久久久国产精品人妻aⅴ院| 亚洲国产精品999在线| 在线观看美女被高潮喷水网站 | 很黄的视频免费| 人人妻人人看人人澡| 嫩草影视91久久| 亚洲av日韩精品久久久久久密| 一进一出抽搐动态| 成人永久免费在线观看视频| 国产精品香港三级国产av潘金莲| 十八禁人妻一区二区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 97超级碰碰碰精品色视频在线观看| 日日夜夜操网爽| 一进一出抽搐gif免费好疼| 天堂动漫精品| 在线看三级毛片| 香蕉av资源在线| 在线免费观看的www视频| 在线免费观看的www视频| 极品教师在线免费播放| 国内毛片毛片毛片毛片毛片| 久久6这里有精品| 亚洲av免费高清在线观看| 国产视频内射| 亚洲天堂国产精品一区在线| 久久久久久久久中文| 黄色日韩在线| 又粗又爽又猛毛片免费看| 男女下面进入的视频免费午夜| 国产成人福利小说| 色播亚洲综合网| 18禁黄网站禁片午夜丰满| 高清在线国产一区| 日韩欧美在线乱码| 国产老妇女一区| 久久久久久九九精品二区国产| 久久久成人免费电影| 中文字幕人成人乱码亚洲影| 欧美精品啪啪一区二区三区| 在线天堂最新版资源| 人妻夜夜爽99麻豆av| 亚洲精品一区av在线观看| 免费av毛片视频| 国产精品自产拍在线观看55亚洲| 51国产日韩欧美| 国产激情欧美一区二区| 精品久久久久久久人妻蜜臀av| 天天一区二区日本电影三级| 国产成人系列免费观看| 国产成人av激情在线播放| 99国产精品一区二区三区| 麻豆久久精品国产亚洲av| 亚洲av二区三区四区| 一级a爱片免费观看的视频| 国产黄a三级三级三级人| 久久久久久久精品吃奶| e午夜精品久久久久久久| 制服人妻中文乱码| 国产精品野战在线观看| av天堂中文字幕网| 亚洲中文字幕日韩| 最好的美女福利视频网| 亚洲av免费在线观看| 午夜亚洲福利在线播放| 毛片女人毛片| 亚洲不卡免费看| 久久久久免费精品人妻一区二区| 搡老妇女老女人老熟妇| 亚洲精品亚洲一区二区| 日韩人妻高清精品专区| 国产蜜桃级精品一区二区三区| 丰满乱子伦码专区| 性色avwww在线观看| 99国产极品粉嫩在线观看| 久久精品国产清高在天天线| 欧美性猛交黑人性爽| 亚洲人与动物交配视频| xxx96com| 久久国产精品影院| 老熟妇仑乱视频hdxx| 五月伊人婷婷丁香| 国产黄a三级三级三级人| 成人无遮挡网站| 亚洲精品456在线播放app | 亚洲av免费在线观看| 日韩av在线大香蕉| 在线看三级毛片| 99热精品在线国产| 国产老妇女一区| 在线观看一区二区三区| 香蕉丝袜av| x7x7x7水蜜桃| 丁香六月欧美| 国产精品综合久久久久久久免费| 亚洲激情在线av| 狠狠狠狠99中文字幕| 18禁黄网站禁片午夜丰满| 两个人看的免费小视频| 免费av不卡在线播放| 亚洲国产精品999在线| 一个人免费在线观看的高清视频| 男人的好看免费观看在线视频| 国产精品久久电影中文字幕| 国产午夜精品论理片| 国产三级中文精品| 99精品欧美一区二区三区四区| 国产精品 国内视频| 精品日产1卡2卡| 成人永久免费在线观看视频| 99国产极品粉嫩在线观看| 欧美性猛交黑人性爽| av女优亚洲男人天堂| 夜夜躁狠狠躁天天躁| 一进一出抽搐gif免费好疼| 麻豆国产97在线/欧美| 最近最新免费中文字幕在线| 日韩欧美精品免费久久 | 老司机在亚洲福利影院| 中文资源天堂在线| 欧美不卡视频在线免费观看| 在线十欧美十亚洲十日本专区| 亚洲人成网站高清观看| 色综合婷婷激情| 国产精品一区二区免费欧美| 国产成人a区在线观看| 国产又黄又爽又无遮挡在线| 级片在线观看| 久久草成人影院| 日本黄色片子视频| 亚洲欧美激情综合另类| 女人被狂操c到高潮| 特大巨黑吊av在线直播| 一级毛片女人18水好多| 三级国产精品欧美在线观看| 国产成年人精品一区二区| 俺也久久电影网| 中文字幕人妻熟人妻熟丝袜美 | 亚洲av熟女| 精品午夜福利视频在线观看一区| 国产精品99久久久久久久久| 亚洲精品在线观看二区| 噜噜噜噜噜久久久久久91| 在线免费观看不下载黄p国产 | 亚洲五月天丁香| 99热这里只有精品一区| 婷婷精品国产亚洲av在线| 亚洲精华国产精华精| 一级毛片高清免费大全| 色吧在线观看| 欧美成人免费av一区二区三区| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 看黄色毛片网站| 国产免费av片在线观看野外av| 国产成人aa在线观看| 免费无遮挡裸体视频| 国产一区二区亚洲精品在线观看| 亚洲不卡免费看| 此物有八面人人有两片| 午夜福利视频1000在线观看| 成人国产综合亚洲| 国内毛片毛片毛片毛片毛片| 天堂av国产一区二区熟女人妻| 国内揄拍国产精品人妻在线| 首页视频小说图片口味搜索| 午夜福利高清视频| 国产精华一区二区三区| 黄色日韩在线| 国产精品精品国产色婷婷| 九色国产91popny在线| 精品人妻1区二区| 久久久久九九精品影院| 国产精品一区二区三区四区免费观看 | 美女大奶头视频| 亚洲av免费高清在线观看| 亚洲最大成人中文| 乱人视频在线观看| 99久久精品一区二区三区| 精品久久久久久,| 真人一进一出gif抽搐免费| 色老头精品视频在线观看| 国产国拍精品亚洲av在线观看 | 国产三级在线视频| 国产麻豆成人av免费视频| 亚洲欧美日韩高清在线视频| 国产高清激情床上av| 亚洲人成伊人成综合网2020| 三级国产精品欧美在线观看| 黄片大片在线免费观看| 免费av不卡在线播放| 18禁在线播放成人免费| 午夜视频国产福利| 99在线人妻在线中文字幕| 日韩欧美在线二视频| 亚洲精品一区av在线观看| 精品欧美国产一区二区三| 国产精品一区二区三区四区免费观看 | 日韩亚洲欧美综合| 久久香蕉精品热| 午夜亚洲福利在线播放| 中文字幕人妻熟人妻熟丝袜美 | 女人被狂操c到高潮| 成人av一区二区三区在线看| 18禁裸乳无遮挡免费网站照片| 熟妇人妻久久中文字幕3abv| 欧美性猛交╳xxx乱大交人| 人妻丰满熟妇av一区二区三区| 一a级毛片在线观看| 97超级碰碰碰精品色视频在线观看| 国产一区二区三区视频了| 国产精品 国内视频| 国产成+人综合+亚洲专区| 女生性感内裤真人,穿戴方法视频| 中文字幕人妻熟人妻熟丝袜美 | 精品久久久久久久久久免费视频| 一进一出好大好爽视频| 国内久久婷婷六月综合欲色啪| 无遮挡黄片免费观看| 99精品久久久久人妻精品| 国产精品1区2区在线观看.| 在线观看美女被高潮喷水网站 | 又黄又爽又免费观看的视频| 国产黄色小视频在线观看| 国产欧美日韩精品一区二区| 亚洲熟妇中文字幕五十中出| 国产成人a区在线观看| 婷婷六月久久综合丁香| 内地一区二区视频在线| 午夜a级毛片| 国产成人欧美在线观看| 午夜影院日韩av| 亚洲人与动物交配视频| 麻豆久久精品国产亚洲av| 亚洲精品久久国产高清桃花| 在线十欧美十亚洲十日本专区| 色综合站精品国产| 国内精品久久久久精免费| 久久亚洲精品不卡| 免费看日本二区| 中文字幕av成人在线电影| 午夜免费激情av| 又黄又爽又免费观看的视频| 性色av乱码一区二区三区2| 成人一区二区视频在线观看| 国产精品一区二区三区四区久久| 久久国产精品人妻蜜桃| 国产精品久久久久久精品电影| www国产在线视频色| 亚洲 国产 在线| 一本久久中文字幕| 又紧又爽又黄一区二区| 国模一区二区三区四区视频| 国产精品一区二区三区四区免费观看 | 免费观看的影片在线观看| 亚洲最大成人中文| 亚洲内射少妇av| 少妇的丰满在线观看| 色老头精品视频在线观看| 日韩人妻高清精品专区| 香蕉av资源在线| 久久久久国内视频| 日本一本二区三区精品| 欧美+亚洲+日韩+国产| 日本免费a在线| 国产综合懂色| 久久久久亚洲av毛片大全| 国内毛片毛片毛片毛片毛片| 国产精品久久电影中文字幕| 一区二区三区激情视频| 国产成人福利小说| 国产亚洲欧美98| 国产精品99久久99久久久不卡| 亚洲精品粉嫩美女一区| 欧美日韩亚洲国产一区二区在线观看| 亚洲一区二区三区色噜噜| 日本a在线网址| 国产精品一区二区免费欧美| 亚洲美女黄片视频| 午夜福利视频1000在线观看| 夜夜躁狠狠躁天天躁| 又粗又爽又猛毛片免费看| 欧美成人性av电影在线观看| 日韩欧美国产在线观看| 老熟妇仑乱视频hdxx| 91在线精品国自产拍蜜月 | 18禁黄网站禁片免费观看直播| 黄片大片在线免费观看| 久久久久免费精品人妻一区二区| 国产免费av片在线观看野外av| 好男人在线观看高清免费视频| netflix在线观看网站| 搡女人真爽免费视频火全软件 | 久久6这里有精品| 日日摸夜夜添夜夜添小说| 宅男免费午夜| 99久国产av精品| 亚洲专区中文字幕在线| 一区二区三区免费毛片| 中文字幕av在线有码专区| 好男人电影高清在线观看| 制服丝袜大香蕉在线| 国产成人啪精品午夜网站| 久久久久亚洲av毛片大全| 久久精品91无色码中文字幕| 中文字幕高清在线视频| av国产免费在线观看| 久久九九热精品免费| 特级一级黄色大片| 老汉色av国产亚洲站长工具| 国产精品久久久久久久久免 | 国产午夜精品论理片| 日韩欧美三级三区| 亚洲第一欧美日韩一区二区三区| 国产一区二区三区视频了| 69人妻影院| 免费人成在线观看视频色| 日本 av在线| aaaaa片日本免费| 中文亚洲av片在线观看爽| 免费看美女性在线毛片视频| 久久久精品大字幕| 日本 av在线| 日本黄色片子视频| 国产高清视频在线播放一区| 色精品久久人妻99蜜桃| 国产成+人综合+亚洲专区| 长腿黑丝高跟| 国产视频内射| 少妇熟女aⅴ在线视频| 亚洲18禁久久av| 中出人妻视频一区二区| ponron亚洲| 亚洲五月天丁香| 亚洲无线观看免费| АⅤ资源中文在线天堂| 内地一区二区视频在线| 精品人妻1区二区| 国产精品永久免费网站| 色在线成人网| 欧美日韩一级在线毛片| 一级黄色大片毛片| 18禁美女被吸乳视频| 美女被艹到高潮喷水动态| 香蕉av资源在线| 欧美成人免费av一区二区三区| 亚洲av成人不卡在线观看播放网| 亚洲精品久久国产高清桃花| 国产精品亚洲美女久久久| 久久久久久国产a免费观看| 久久伊人香网站| 欧美+亚洲+日韩+国产| 欧美成狂野欧美在线观看| 亚洲中文字幕一区二区三区有码在线看| 久9热在线精品视频| 成人国产一区最新在线观看| 亚洲第一欧美日韩一区二区三区| 国产免费av片在线观看野外av| 国产aⅴ精品一区二区三区波| 国产高清视频在线观看网站| 午夜免费男女啪啪视频观看 | 黄色日韩在线| 亚洲欧美日韩高清专用| 国产av在哪里看| 欧美一级毛片孕妇| 91麻豆精品激情在线观看国产| 夜夜躁狠狠躁天天躁| 琪琪午夜伦伦电影理论片6080| 中文资源天堂在线| 99久久无色码亚洲精品果冻| 国产精品久久久久久久电影 | 欧美+亚洲+日韩+国产| 给我免费播放毛片高清在线观看| 人人妻人人看人人澡| 在线免费观看的www视频| 日本免费一区二区三区高清不卡| 久久久精品欧美日韩精品| 亚洲无线观看免费| 久久久成人免费电影| 亚洲第一欧美日韩一区二区三区| 日本免费一区二区三区高清不卡| 美女cb高潮喷水在线观看| av天堂中文字幕网| 成人午夜高清在线视频| 亚洲中文日韩欧美视频| 中文字幕人妻丝袜一区二区| 亚洲国产精品久久男人天堂| 亚洲自拍偷在线| aaaaa片日本免费| 97碰自拍视频| 精品熟女少妇八av免费久了| 人妻夜夜爽99麻豆av| 午夜免费激情av| 欧美日韩瑟瑟在线播放| 欧美性感艳星| 国产亚洲精品综合一区在线观看| 午夜福利免费观看在线| 一区二区三区国产精品乱码| 99精品欧美一区二区三区四区| 一本精品99久久精品77| 一a级毛片在线观看| 人妻久久中文字幕网| 一a级毛片在线观看| 99久久综合精品五月天人人| 国产真实伦视频高清在线观看 | 亚洲熟妇中文字幕五十中出| 中出人妻视频一区二区| 免费看a级黄色片| 三级毛片av免费| 免费无遮挡裸体视频| 69人妻影院| 国产精品乱码一区二三区的特点| 91在线精品国自产拍蜜月 | 免费大片18禁| 99久久精品一区二区三区| 波多野结衣高清作品| 真实男女啪啪啪动态图| 久久久成人免费电影| 亚洲第一电影网av| 草草在线视频免费看| 亚洲国产中文字幕在线视频| 一a级毛片在线观看| 校园春色视频在线观看| 天堂动漫精品| 国产免费男女视频| 日本一二三区视频观看| 人人妻人人澡欧美一区二区| 高清在线国产一区| 九九热线精品视视频播放| 国产精品一区二区三区四区免费观看 | 亚洲av成人精品一区久久| 综合色av麻豆| 精品福利观看| 别揉我奶头~嗯~啊~动态视频| 日韩成人在线观看一区二区三区| 制服人妻中文乱码| 最近视频中文字幕2019在线8| 国产黄片美女视频| 女警被强在线播放| 一a级毛片在线观看| 一级毛片女人18水好多| 91在线精品国自产拍蜜月 | 日本黄色片子视频| 国产私拍福利视频在线观看| 在线播放无遮挡| 99久久无色码亚洲精品果冻| 久久久久久九九精品二区国产| 亚洲欧美精品综合久久99| 亚洲,欧美精品.| 亚洲欧美日韩高清在线视频| АⅤ资源中文在线天堂| 啦啦啦观看免费观看视频高清| 亚洲精品粉嫩美女一区| 欧美不卡视频在线免费观看| 欧美日韩福利视频一区二区| 色噜噜av男人的天堂激情| 国内精品久久久久精免费| 一个人免费在线观看电影| 欧美三级亚洲精品| 国产单亲对白刺激| 天堂动漫精品| 成人高潮视频无遮挡免费网站| 嫩草影院入口| 在线天堂最新版资源| 毛片女人毛片| 国内精品美女久久久久久| 麻豆国产av国片精品| 日韩欧美 国产精品| 黄色片一级片一级黄色片| 亚洲精品亚洲一区二区| 成年免费大片在线观看| 美女大奶头视频| 亚洲人成网站在线播| 国产免费av片在线观看野外av| 最新美女视频免费是黄的| 久久久久久久久久黄片| 精品久久久久久久久久免费视频| 国产国拍精品亚洲av在线观看 | 成人国产一区最新在线观看| aaaaa片日本免费| 丁香欧美五月| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产清高在天天线| 亚洲人成网站高清观看| 有码 亚洲区| 人人妻人人看人人澡| 免费大片18禁| 免费无遮挡裸体视频| 女警被强在线播放| 亚洲精品日韩av片在线观看 | 日本成人三级电影网站| 老司机午夜十八禁免费视频| 国产精品一及| 欧美绝顶高潮抽搐喷水| 最近在线观看免费完整版| 国产高潮美女av| 国产视频内射| 精华霜和精华液先用哪个| 欧美日韩亚洲国产一区二区在线观看| 男女下面进入的视频免费午夜| 亚洲avbb在线观看| 国产精品av视频在线免费观看| 精品人妻一区二区三区麻豆 | 国产一区二区亚洲精品在线观看| 99久久精品一区二区三区| 国产91精品成人一区二区三区| 亚洲天堂国产精品一区在线| 99在线人妻在线中文字幕| 精品久久久久久久人妻蜜臀av| 一个人免费在线观看的高清视频| 亚洲国产日韩欧美精品在线观看 | 国产激情欧美一区二区| 国产私拍福利视频在线观看| 又黄又粗又硬又大视频| 别揉我奶头~嗯~啊~动态视频| 搡女人真爽免费视频火全软件 | 久久久久久久久久黄片| 在线十欧美十亚洲十日本专区| 可以在线观看的亚洲视频| 九九久久精品国产亚洲av麻豆| 欧美乱码精品一区二区三区| 午夜福利视频1000在线观看| eeuss影院久久| 99久久精品一区二区三区| 亚洲精品久久国产高清桃花| 少妇人妻一区二区三区视频| 欧美最黄视频在线播放免费| 国产精品 国内视频| 男女午夜视频在线观看| 中文字幕久久专区| 中亚洲国语对白在线视频| 最好的美女福利视频网| 高潮久久久久久久久久久不卡| 99国产精品一区二区蜜桃av| 神马国产精品三级电影在线观看| 亚洲在线自拍视频| 免费人成视频x8x8入口观看| 人人妻,人人澡人人爽秒播| 日本一本二区三区精品| 亚洲黑人精品在线| 亚洲真实伦在线观看| 中文字幕久久专区| 18美女黄网站色大片免费观看| 日韩成人在线观看一区二区三区| 丁香六月欧美| 嫁个100分男人电影在线观看| 国产一区二区在线观看日韩 |