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

    甘青青蘭根際土壤理化性質(zhì)及微生物群落結(jié)構(gòu)特征分析

    2024-06-16 01:43:23張二豪劉盼盼何萍簡閱徐雨婷陳誠欣祿亞洲蘭小中索朗桑姆
    關(guān)鍵詞:微生物群落土壤理化性質(zhì)高通量測序

    張二豪 劉盼盼 何萍 簡閱 徐雨婷 陳誠欣 祿亞洲 蘭小中 索朗桑姆

    摘要:為研究甘青青蘭根際土壤理化性質(zhì)與根際土壤細(xì)菌和真菌群落結(jié)構(gòu)組成特征,以西藏自治區(qū)工布江達(dá)縣(GB)、卡諾區(qū)(KR)和洛隆縣(LL)的甘青青蘭根際土壤微生物為研究對象,采用Illumina MiSeq高通量測序技術(shù)對根際土壤細(xì)菌和真菌群落結(jié)構(gòu)組成進(jìn)行分析,同時(shí)測定根際土壤的理化性質(zhì),并與根際土壤核心微生物菌群進(jìn)行相關(guān)性分析。結(jié)果表明,不同地區(qū)甘青青蘭根際土壤理化性質(zhì)差異顯著。甘青青蘭根際土壤共獲得3 900個(gè)細(xì)菌OTUs和1 990個(gè)真菌OTUs,不同地區(qū)間土壤微生物多樣性存在明顯差異。在門水平,GB、KR 和LL 樣品的優(yōu)勢細(xì)菌門均為放線菌門(Actinobacteria),GB 樣品的優(yōu)勢真菌門是擔(dān)子菌門(Basidiomycota),KR和LL樣品的優(yōu)勢真菌門是子囊菌門(Ascomycota);在屬水平,不同地區(qū)甘青青蘭根際土壤優(yōu)勢微生物菌群存在明顯差異。主成分分析顯示,不同地區(qū)甘青青蘭根際土壤細(xì)菌和真菌群落結(jié)構(gòu)組成差異較大。核心微生物菌群分析表明,甘青青蘭根際土壤核心細(xì)菌菌群有257個(gè)屬;核心真菌菌群有102個(gè)屬。相關(guān)性分析表明,核心微生物菌群的改變與根際土壤理化因子間存在不同程度的相關(guān)性,且土壤全鉀、全磷、全氮、速效鉀和速效氮含量是影響甘青青蘭根際土壤微生物群落結(jié)構(gòu)組成的重要因子。綜上所述,不同地區(qū)甘青青蘭根際土壤微生物群落結(jié)構(gòu)組成存在顯著差異,且根際土壤微生物群落結(jié)構(gòu)的變化與土壤理化因子密切相關(guān),為甘青青蘭人工種植及篩選有益微生物菌群提供了理論依據(jù)。

    關(guān)鍵詞:甘青青蘭;土壤理化性質(zhì);高通量測序;微生物群落;核心微生物;相關(guān)性分析

    doi:10.13304/j.nykjdb.2022.0741

    中圖分類號:S154.3 文獻(xiàn)標(biāo)志碼:A 文章編號:1008‐0864(2024)03‐0201‐13

    甘青青蘭(Dracocephalum tanguticum Maxim.)又名“唐古特青蘭”,是唇形科(Lamiaceae)青蘭屬(Dracocephalum)的一種多年生草本植物,主要分布在海拔1 900~4 600 m的西藏東南、青海東部和四川西部等地。甘青青蘭作為一種藏醫(yī)藥臨床常用藏藥材,具有清熱解毒、止咳化痰、和胃舒肝等功效[1]。研究表明,甘青青蘭含有豐富的藥用化學(xué)成分,如黃酮及黃酮苷類、萜類、揮發(fā)油、多糖和無機(jī)元素等[2-4],具有抗缺氧、降血糖、抗菌和消炎等藥理活性[5-8]。

    藥用植物的品質(zhì)與植物品種、根際土壤微生物、土壤理化性質(zhì)及氣候等因素密切相關(guān)[9‐10]。根際土壤是指受植物根系及根系分泌物影響的離根軸數(shù)毫米的微區(qū)域,根際土壤微生物是根際土壤微生態(tài)的重要組成部分,在物質(zhì)循環(huán)、能量轉(zhuǎn)換、植物生長發(fā)育、抗病和抗逆境等方面扮演重要角色[11]。研究表明,紫萁小菇(Mycena osmundicola)和蜜環(huán)菌(Armillariella mellea)能促進(jìn)天麻(Gastrodia elata Bl.)種子的萌發(fā)[12];假單胞菌屬(Pseudomonas)、醋桿菌屬(Acrtobacter)、腸桿菌屬(Enterobacter)和芽孢桿菌屬(Bacillus)等能促進(jìn)藥用植物對土壤養(yǎng)分的吸收和利用,進(jìn)而促進(jìn)植物的生長發(fā)育[11]。Kim等[13]從人參根際分離的解淀粉芽孢桿菌(Bacillus amyloliquefaciens)AK‐0能抑制人參銹腐病菌(Cylindrocarpon destructans)的繁殖。Egamberdieva 等[14] 研究表明根瘤菌(Mesorhizobium sp.)和假單胞菌(Pseudomonasextremorientalis)能增強(qiáng)甘草對鹽的耐受性。研究表明,叢枝菌根真菌能促進(jìn)藥用植物中酚類、萜類和生物堿的合成和積累[15];枯草芽孢桿菌(Bacillus subtilis)能提高黨參的生物量[16];印度梨狀芽孢桿菌能促進(jìn)藥用植物黃花蒿中青蒿素的合成[17]。根際土壤微生物群落結(jié)構(gòu)的變化與土壤理化性質(zhì)密切相關(guān)。研究表明,土壤養(yǎng)分含量、pH和含水率是影響珍珠豬毛菜根際土壤細(xì)菌群落結(jié)構(gòu)組成的重要因子[18];土壤pH及氮、磷含量等影響喜馬拉雅紫茉莉根際土壤細(xì)菌的群落結(jié)構(gòu)及多樣性[19]。因此,研究藥用根際土壤微生物群落結(jié)構(gòu)組成及根際土壤理化性質(zhì),對藥材種植及分離篩選益生微生物類群具有重要意義。

    甘青青蘭作為一種傳統(tǒng)藏藥材,其相關(guān)研究主要集中在其藥用化學(xué)成分及藥理學(xué)活性等方面,而有關(guān)根際土壤微生物群落結(jié)構(gòu)組成及土壤理化性質(zhì)的研究尚未見報(bào)道。為進(jìn)一步探討甘青青蘭根際土壤微生物群落結(jié)構(gòu)組成差異與土壤理化性質(zhì)的關(guān)聯(lián)性,本研究采用高通量測序技術(shù)結(jié)合土壤理化性質(zhì)測定,對3個(gè)地區(qū)的甘青青蘭根際土壤微生物群落結(jié)構(gòu)組成進(jìn)行分析,旨在揭示其核心微生物類群與土壤理化性質(zhì)間的關(guān)系,以期為甘青青蘭人工種植和益生微生物的篩選提供理論依據(jù)。

    1 材料與方法

    1.1 樣本采集

    2020 年8 月于西藏昌都市卡若區(qū)(Karuo,KR)、昌都市洛隆縣(Luolong,LL)和林芝市工布江達(dá)縣(Gongbujiangda,GB)分別采集2年生甘青青蘭根際土壤樣品(表1)。在每個(gè)采樣點(diǎn),根據(jù)隨機(jī)采樣原則,隨機(jī)選擇9株長勢良好的甘青青蘭植株,采用“抖根法”收集根際土壤樣品,每個(gè)采樣點(diǎn)3個(gè)重復(fù),并過直徑為2 mm的過濾篩以除去殘枝落葉,裝入無菌袋并混合均勻,置于液氮中并帶回實(shí)驗(yàn)室,-80 ℃保存?zhèn)溆谩?/p>

    1.2 根際土壤理化性質(zhì)分析

    參照《土壤農(nóng)化分析》[20]操作方法測定根際土壤pH、總磷(total phosphorus,TP)、速效磷(available phosphorus, AP)、總鉀(total potassium,TK)、速效鉀(available potassium,AK)、總氮(totalnitrogen,TN)、速效氮(available nitrogen,AN)和土壤有機(jī)質(zhì)(soil organic matter,SOM)含量。

    1.3 根際土壤樣品總DNA 提取、PCR 擴(kuò)增及Illumina 測序

    根據(jù)土壤基因組DNA提取試劑盒(TIANGEN)操作說明提取根際土壤細(xì)菌和真菌總DNA,1%瓊脂糖凝膠電泳檢測DNA質(zhì)量,用紫外分光光度計(jì)NanoDrop ND-2000(Thermo Scientific, Wilmington,USA)測定DNA 含量和純度。利用通用引物338F(5-ACTCCTACGGGAGGCAGCA-3)和806R(5-GGACTACHVGGGTWTCTAAT-3)擴(kuò)增細(xì)菌16S rRNA 基因V3~V4 區(qū),利用通用引物ITS1F(CTTGGTCATTTAGAGGAAGTAA)和ITS2R(GCTGCGTTCTTCATCGATGC)擴(kuò)增真菌內(nèi)源轉(zhuǎn)錄間隔區(qū)(internally transcribed spacer,ITS)。PCR體系為20 μL,包含2×Taq PCR MasterMix 10 μL、10 mmol·L-1正反向引物各1 μL、基因組DNA模板1 μL、ddH2O 7 μL。PCR 程序:94 ℃ 5 min;94 ℃30 s,55 ℃ 30 s,72 ℃ 1 min,35 個(gè)循環(huán);72 ℃10 min,4 ℃保存。2%瓊脂糖凝膠電泳檢測PCR擴(kuò)增產(chǎn)物,膠回收擴(kuò)增產(chǎn)物后委托上海派森諾科技股份有限公司進(jìn)行高通量測序。

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

    利用QIIME2(2019.4)對Illumina MiSeq 測序數(shù)據(jù)進(jìn)行去噪、拼接和質(zhì)控;在97%相似度水平下,利用Vsearch(v2.13.4)對高質(zhì)量序列進(jìn)行OTU(operational taxonomic unit)聚類分析;基于Greengenes 數(shù)據(jù)庫(Release 13.8,http://greengenes.secondgenome. com/)和UNITE 數(shù)據(jù)庫(Release8.0,https://unite.ut.ee/)分別對根際土壤細(xì)菌和真菌OTUs 代表序列進(jìn)行物種分類學(xué)注釋;利用QIIME2(2019.4)和R 語言進(jìn)行分類學(xué)組成、α 多樣性和β多樣性分析。其中α多樣性包括Chao1指數(shù)、Simpson 指數(shù)、Shannon 指數(shù)和Shannoneven指數(shù);β 多樣性分析包括主成分分析(principalcomponent analysis,PCA)和樣本層級聚類。基于R語言進(jìn)行相關(guān)性Heatmap分析;利用SPSS 25.0軟件對土壤理化性質(zhì)和alpha多樣性指數(shù)進(jìn)行單因素方差分析。

    2 結(jié)果與分析

    2.1 樣品根際土壤理化性質(zhì)分析

    由圖1可知,KR 樣品pH 呈偏堿性,而LL 和GB樣品呈偏酸性;LL樣品中的SOM含量顯著高于GB 樣品,但與KR 樣品無顯著差異,且KR 與GB間無顯著差異;GB、KR和LL樣品中TN和TP含量差異顯著,且KR 中含量最高,GB中含量最低,KR樣品中AN和AP含量顯著高于GB,但KR樣品與LL樣品中AN含量無顯著差異,GB和LL樣品中AP含量無顯著差異;GB、KR和LL樣品中TK含量差異顯著,且LL樣品中最高,KR次之,GB最低;KR和LL樣品中AK含量最高,且顯著高于GB樣品,但樣品KR和LL間無顯著差異。由此表明,不同地區(qū)甘青青蘭根際土壤理化性質(zhì)差異顯著。

    2.2 樣品序列統(tǒng)計(jì)及alpha 多樣性分析

    高通量測序結(jié)果(表2)表明,3個(gè)地區(qū)甘青青蘭根際土壤細(xì)菌共產(chǎn)生48 763條有效序列,其中,GB、KR和LL分別產(chǎn)生15 316、16 714和16 733條有效序列,序列平均長度376 bp,測序覆蓋度為97.11%~97.78%;根際土壤真菌共產(chǎn)生143 744條有效序列,其中,GB、KR 和LL 分別產(chǎn)生42 048、51 312和50 384條有效序列,序列平均長度260 bp,測序覆蓋度為99.67%~99.85%。由此說明,測序結(jié)果能基本反應(yīng)樣品真實(shí)微生物群落結(jié)構(gòu)組成。

    在97%相似度水平下,3個(gè)地區(qū)甘青青蘭根際土壤共獲得3 900個(gè)細(xì)菌OTUs和1 990個(gè)真菌OTUs,其中,GB、KR、LL樣品分別注釋到的細(xì)菌OTUs 和真菌OTUs 分別為1 149、1 473、1 278 和489、883和618個(gè)。韋恩圖(圖2)分析顯示,3個(gè)地區(qū)共有的細(xì)菌OTUs和真菌OTUs分別為667和91個(gè),分別占總數(shù)的17.10%和4.57%;GB、KR和LL樣品特有的細(xì)菌OTUs分別為190、417和280個(gè),特有的真菌OTUs分別為242、556和347個(gè)。由此說明,不同地區(qū)甘青青蘭根際土壤微生物群落結(jié)構(gòu)組成存在差異。

    α多樣性結(jié)果(圖3)表明,KR樣品根際土壤細(xì)菌和真菌的豐富度指數(shù)最高,顯著高于LL 和GB,GB樣品根際土壤細(xì)菌和真菌的豐富度指數(shù)最低,顯著低于LL樣品;KR樣品根際土壤細(xì)菌多樣性指數(shù)最高,顯著高于GB樣品,但與LL樣品無顯著差異;LL 樣品根際土壤真菌多樣性指數(shù)最高,但與KR樣品無顯著差異,GB樣品多樣性指數(shù)最低,且顯著低于LL和KR樣品。以上結(jié)果表明,不同地區(qū)甘青青蘭根際土壤微生物α多樣性存在明顯差異。

    2.3 樣品微生物種群歸類分析

    由表3可知,根際土壤細(xì)菌共注釋到27門、74綱、177目、269科、492屬,其中GB樣品注釋到26門、62綱、133目、197科、360屬;KR樣品注釋到26門、68綱、159目、237科、415屬;LL樣品注釋到23門、55綱、125目、193科、332屬。根際土壤真菌共注釋到11門、38綱、99目、213科、387屬,其中GB樣品注釋到10門、28綱、66目、132科、190屬;KR樣品注釋到9門、31綱、83目、163科、276屬;LL樣品注釋到8門、29綱、62目、137科、225屬。不同種源地甘青青蘭根際土壤微生物在各分類等級上的數(shù)量存在明顯差異,說明不同種源地甘青青蘭根際土壤微生物群落結(jié)構(gòu)組成存在差異。

    2.4 樣品微生物群落結(jié)構(gòu)組成分析

    在門水平,3個(gè)地區(qū)甘青青蘭根際土壤細(xì)菌和真菌群落結(jié)構(gòu)組成如圖4所示。對于細(xì)菌群落,放線菌門(Actinobacteriota)、變形菌門(Proteobacteria)、厚壁菌門(Firmicutes)和綠彎菌門(Chloroflexi)的相對豐度較高。其中,在GB樣品中,相對豐度>1%的有3個(gè)菌門,分別是放線菌門(45.41%)、變形菌門(19.20%)和厚壁菌門(30.63%);在KR 和LL 樣品中,相對豐度>1%的有4個(gè)菌門,分別是放線菌門(57.84%和67.44%)、變形菌門(28.92%和16.69%)、厚壁菌門(3.84% 和3.57%)和綠彎菌門(3.09% 和6.06%);說明3個(gè)地區(qū)甘青青蘭根際土壤優(yōu)勢細(xì)菌門類型基本相同,但相對豐度差異較大。對于真菌群落,GB 樣品中的優(yōu)勢真菌門為擔(dān)子菌門(Basidiomycota,85.43%)、子囊菌門(Ascomycota,10.01%)和被孢霉門(Mortierellomycota,2.84%);KR和LL 樣品中的優(yōu)勢真菌門是子囊菌門(67.86%和75.25%)、擔(dān)子菌門(19.89% 和15.89%)、被孢霉門(10.83% 和7.17%)和unclassified_k_Fungi(1.21%和1.10%),說明3個(gè)地區(qū)甘青青蘭根際土壤中的優(yōu)勢真菌門組成相似,但相對豐度差異較大。

    在屬水平,3個(gè)地區(qū)甘青青蘭根際土壤細(xì)菌和真菌群落結(jié)構(gòu)組成如圖5所示。在相對豐度>1%條件下,GB樣品有5個(gè)細(xì)菌屬,其中能準(zhǔn)確分類的有4 個(gè),分別是膨脹芽孢桿菌屬(Tumebacillus, 8.46%) 、短芽孢桿菌屬(Brevibacillus, 8.43%) 、土壤紅桿菌屬(Solirubrobacter,7.22%)和芽孢桿菌屬(Bacillus,6.03%);KR和LL樣品各有3個(gè)細(xì)菌屬,其中能準(zhǔn)確分類的有2 個(gè),分別是土壤紅桿菌屬(Solirubrobacter,6.05% 和8.37%)和芽孢桿菌屬(2.56%和2.52%),說明KR和LL樣品中的優(yōu)勢細(xì)菌屬相同,但與GB樣品存在明顯差異。在相對豐度>1%條件下,GB樣品有6個(gè)真菌屬,分別是濕傘屬(Hygrocybe,78.31%) 、產(chǎn)油菌屬(Solicoccozyma,6.40%)、被孢霉屬(Mortierella,2.84%)、毛殼菌屬(Chaetomium,1.80%)、錐毛殼屬(Coniochaeta,1.29%)和青霉菌屬(Penicillium,1.16%);KR樣品中有15個(gè)真菌屬,其中能準(zhǔn)確分類的有11個(gè),分別是產(chǎn)油菌屬(15.22%)、被孢霉屬(10.83%)、赤霉菌屬(Gibberella,7.74%)、寡囊盤菌屬(Thelebolus,6.51%) 、蠟盤菌屬(Rutstroemia,6.03%)、青霉菌屬(3.83%)、枝孢瓶霉屬(Cladophialophora,3.29%)、錐毛殼屬(3.14%)、Knufia(2.20%)、Dactylonectria(1.70%)和外瓶柄霉屬(Exophiala,1.28%);LL樣品中有15個(gè)真菌屬,其中能準(zhǔn)確分類的有9 個(gè),分別是Montagnula(14.75%)、被孢霉屬(7.17%)、背芽突霉屬(Cadophora,6.13%)、枝孢瓶霉屬(5.62%)、瓶霉屬(Phialophora,4.19%)、鐮刀菌屬(Fusarium,2.95%)、Knufia(2.73%)、產(chǎn)油菌屬(2.32%)和錐毛殼屬(1.92%)。以上結(jié)果表明,不同種源地甘青青蘭根際土壤優(yōu)勢真菌屬存在明顯差異。

    2.5 樣品微生物β 多樣性分析

    為了分析不同樣品間根際土壤細(xì)菌和真菌群落結(jié)構(gòu)組成差異,對甘青青蘭根際土壤微生物進(jìn)行基于非加權(quán)組平均法(unweighted pair groupmethod with arithmetic mean,UPGMA)的樣本層級聚類和主成分分析(principal component analysis,PCA),結(jié)果如圖6和圖7所示。由圖6可知,每個(gè)地區(qū)根際土壤細(xì)菌和真菌單獨(dú)聚為1個(gè)類群,相對而言,KR和LL樣品距離較近。由圖7可知,主成分1(PC1)和主成分2(PC2)對根際土壤細(xì)菌的差異貢獻(xiàn)率分別為33.09% 和26.08%,累計(jì)貢獻(xiàn)率為59.17%;對根際土壤真菌的差異貢獻(xiàn)率分別為87.64%和9.2%,累計(jì)貢獻(xiàn)率為96.84%,能較好地區(qū)分不同地區(qū)甘青青蘭根際土壤樣品,說明不同地區(qū)甘青青蘭根際土壤細(xì)菌和真菌群落結(jié)構(gòu)組成差異較大。

    2.6 樣品微生物核心菌群分析

    甘青青蘭根際土壤細(xì)菌和真菌在屬分類水平韋恩圖見圖8。由圖8可知,甘青青蘭根際土壤細(xì)菌共注釋到492個(gè)屬,其中核心細(xì)菌菌群有257個(gè)屬,占總細(xì)菌屬的52.24%,能準(zhǔn)確分類的核心細(xì)菌菌群有133個(gè)屬;甘青青蘭根際土壤真菌共注釋到387個(gè)屬,其中核心真菌菌群有102個(gè)屬,占總真菌屬的26.36%,能準(zhǔn)確分類的核心真菌菌群有72個(gè)屬。

    為了分析根際土壤理化因子與核心微生物菌群間的相互關(guān)系,基于Pearson算法對根際土壤理化因子與相對豐度前20的核心微生物菌群進(jìn)行相關(guān)性分析,結(jié)果如圖9所示。TK和AK與紅游動菌屬(Rhodoplanes)、Gaiella 屬和假諾卡氏菌屬(Pseudonocardia)等8個(gè)屬的相對豐度呈顯著正相關(guān),而與芽孢桿菌屬(Bacillus)、膨脹芽孢桿菌屬、短芽孢桿菌屬、鏈霉菌屬(Streptomyces)和類芽孢桿菌屬(Paenibacillus)的相對豐度呈顯著負(fù)相關(guān);SOM 與紅游動菌屬、Gaiella 和假諾卡氏菌屬等8個(gè)屬的相對豐度呈顯著正相關(guān),而與芽孢桿菌屬、鏈霉菌屬和假諾卡氏菌屬的相對豐度呈顯著負(fù)相關(guān);TP 與假諾卡氏菌屬、微小桿菌屬(Microvirga)、芽球菌屬(Blastococcus)和節(jié)桿菌屬(Arthrobacter)的相對豐度呈顯著正相關(guān),而與固定桿菌屬(Conexibacter)和慢生根瘤菌屬(Bradyrhizobium)的相對豐度呈顯著負(fù)相關(guān);AP和pH與微小桿菌屬和芽球菌屬的相對豐度呈顯著正相關(guān),而與固定桿菌屬呈顯著負(fù)相關(guān);TN和AN與假諾卡氏菌屬、芽球菌屬和節(jié)桿菌屬等呈顯著正相關(guān),而與固定桿菌屬、 慢生根瘤菌屬和類芽孢桿菌屬等呈顯著負(fù)相關(guān)。錘舌菌屬(Leohumicola) 、Fusicolla、乳突赤殼屬(Thelonectria)、錐毛殼屬和外瓶柄霉屬與AP、pH、TN、TP和AN呈顯著正相關(guān);枝孢瓶霉屬、Knufia和異莖點(diǎn)霉屬(Paraphoma)與TK、AK 和SOM 呈顯著正相關(guān);鐮刀菌屬、瓶霉屬和背芽突霉屬與TK 和SOM 呈顯著正相關(guān);濕傘屬與AK、SOM、TK、TN和AN呈顯著負(fù)相關(guān);毛殼菌屬與AK、TK、和AN呈顯著負(fù)相關(guān)。以上結(jié)果表明,根際土壤理化性質(zhì)影響根際土壤核心微生物菌群的相對豐度。

    2.7 根際土壤微生物多樣性指數(shù)與土壤理化因子的相關(guān)性分析

    由圖10可知,AK、TP、TN 和AN 與根際土壤細(xì)菌的Chao1指數(shù)、香農(nóng)指數(shù)和Shannoneven均勻度指數(shù)呈顯著正相關(guān);TK與根際土壤細(xì)菌的香農(nóng)指數(shù)和Shannoneven 均勻度指數(shù)呈顯著正相關(guān);pH與根際土壤細(xì)菌的Chao1指數(shù)呈顯著正相關(guān);辛普森指數(shù)與根際土壤的TK、AK、TN和AN含量呈顯著負(fù)相關(guān)。根際土壤真菌的香農(nóng)指數(shù)和Shannoneven均勻度指數(shù)與TK、AK、TN和AN含量呈顯著正相關(guān);Chao1 指數(shù)與AP、pH、TN、TP 和AN 呈顯著正相關(guān);辛普森指數(shù)與AK、SOM、TK、TN和AN呈顯著負(fù)相關(guān)。以上結(jié)果表明,根際土壤中AK、TK、TP、TN和AN含量是影響甘青青蘭根際土壤微生物群落多樣性的重要因子。

    3 討論

    微生物作為根際土壤物質(zhì)循環(huán)和能量轉(zhuǎn)換的主要驅(qū)動因子,在植物根際土壤微生態(tài)中發(fā)揮重要作用。植物類型及根際土壤理化因子在很大程度上影響根際土壤微生物群落的結(jié)構(gòu)組成。研究表明,土壤pH、TK、TN、AN、TP和AP是影響喜馬拉雅紫茉莉根際土壤細(xì)菌群落組成和多樣性的重要因子[19]。魏艷晨等[18]研究表明,土壤養(yǎng)分含量、pH和含水率影響珍珠豬毛菜根際土壤細(xì)菌多樣性。李茜等[21]研究發(fā)現(xiàn),根際土壤理化性質(zhì)與不同栽培年限人參根際土壤細(xì)菌群落結(jié)構(gòu)的改變關(guān)系密切。李毳等[22‐23]研究發(fā)現(xiàn),土壤環(huán)境因子和植物種類是影響藥用植物根際土壤微生物群落結(jié)構(gòu)組成和多樣性的重要因素。本研究表明,3個(gè)種源地甘青青蘭根際土壤理化性質(zhì)及根際土壤微生物多樣性存在顯著差異,且根際土壤的AK、TK、TP、TN和AN含量是影響甘青青蘭根際土壤微生物群落多樣性的重要因子。Burns等[24]研究發(fā)現(xiàn),植物種類、土壤理化因子和空間因素對土壤微生物群落多樣性有顯著影響。綜上所述,根際土壤理化因子和空間因素是導(dǎo)致不同地區(qū)甘青青蘭根際土壤微生物群落多樣性差異的關(guān)鍵因素。

    3個(gè)地區(qū)甘青青蘭根際土壤的優(yōu)勢細(xì)菌類群是放線菌門、變形菌門、厚壁菌門和綠彎菌門。研究表明,土壤優(yōu)勢細(xì)菌類群主要集中在變形菌門、放線菌門和厚壁菌門等10個(gè)門[25],與本研究結(jié)果一致。放線菌作為一類重要的生防菌,能產(chǎn)生多種抗生素,進(jìn)而增強(qiáng)植物的抗病能力[26];變形菌門和厚壁菌門具有固氮、土壤修復(fù)和物質(zhì)降解等功能,有助于提高藥用植物的質(zhì)量和抗逆境能力[27‐28]。這3個(gè)菌門的細(xì)菌在甘青青蘭根際土壤中的富集可能參與了甘青青蘭的生長發(fā)育、抗逆和藥用成分合成等。李潔等[29]研究表明,白首烏根際土壤中的優(yōu)勢真菌門是子囊菌門、被孢霉門、擔(dān)子菌門、球囊菌門(Glomeromycota)、類原生動物門(Rozellomycota)和梳霉門(Kickxellomycota),與本研究結(jié)果存在差異。本研究結(jié)果表明,子囊菌門和擔(dān)子菌門是甘青青蘭根際土壤中的優(yōu)勢真菌門,說明植物類型是影響根際土壤真菌類型的關(guān)鍵因素。研究表明,子囊菌門和擔(dān)子菌門是參與土壤有機(jī)質(zhì)分解的主要真菌群落,且與植物根際氮循環(huán)、植物菌根形成及植物生長發(fā)育息息相關(guān)[30]。由此說明,子囊菌門和擔(dān)子菌門的富集有助于甘青青蘭的生長發(fā)育。

    本研究結(jié)果表明,在屬水平,3個(gè)地區(qū)的甘青青蘭根際土壤優(yōu)勢微生物菌群差異明顯,這與樣本層級聚類和主成分分析結(jié)果一致,說明土壤環(huán)境因子和氣候類型等因素是影響根際土壤微生物群落組成的關(guān)鍵因素[24]。核心微生物菌群是一類影響宿主健康和生態(tài)系統(tǒng)功能的關(guān)鍵物種,在植物生長發(fā)育和物質(zhì)代謝調(diào)控中發(fā)揮重要作用[31]。本研究將3個(gè)地區(qū)甘青青蘭根際土壤中共有的微生物確定為核心菌屬,從根際土壤微生物種類和分布規(guī)律來看,甘青青蘭根際土壤核心細(xì)菌菌群種類明顯多于核心真菌類群。芽孢桿菌屬、假諾卡氏菌屬、鏈霉菌屬、慢生根瘤菌屬、鞘氨醇單胞菌屬和節(jié)桿菌屬等是甘青青蘭根際土壤細(xì)菌中的核心菌屬;被孢霉屬、赤霉菌屬、背芽突霉屬、亞隔孢殼屬、青霉菌屬、鐮刀菌屬和毛殼菌屬等是根際土壤真菌中的核心菌屬。研究表明,芽孢桿菌有促進(jìn)植物生長和增加藥用成分合成的作用[32‐33];假諾卡氏菌屬中的許多成員能夠產(chǎn)生抗生素、維生素和酶制劑活性物質(zhì)[34];鞘氨醇單胞菌屬能促進(jìn)大豆的生長和抗?jié)B透脅迫[35];節(jié)桿菌屬和鏈霉菌屬能增強(qiáng)植物的抗逆境脅迫[36‐39];被孢霉屬具有產(chǎn)不飽和脂肪酸和抗植物病原菌的潛在功能;毛殼菌屬在自然界中廣泛存在,具有潛在生防功能[40];青霉菌屬具有抑制植物病害的功效,并能促進(jìn)櫻桃和玉米等植物的生長[41‐42]。以上結(jié)果表明,甘青青蘭根際土壤富含豐富的、具有潛在生物學(xué)功能的核心微生物菌群。根際土壤理化因子與根際土壤核心微生物菌群相關(guān)分析表明,根際土壤理化因子與核心微生物菌群存在一定程度的相關(guān)性,為甘青青蘭的人工種植提供了理論依據(jù)。綜上所述,甘青青蘭根際土壤富含豐富的微生物資源,研究不同地區(qū)甘青青蘭根際土壤微生物群落結(jié)構(gòu)組成,分析核心微生物菌群與土壤理化因子間的相互關(guān)系,對充分挖掘益生微生物資源和指導(dǎo)甘青青蘭人工種植具有重要理論意義。

    參考文獻(xiàn)

    [1] 中華人民共和國衛(wèi)生部藥典委員會.衛(wèi)生部藥品標(biāo)準(zhǔn)·藏

    藥分冊[M].北京:人民衛(wèi)生出版社,1995.

    [2] LI Q, LIU Y Q, HAN L F, et al .. Chemical constituents and

    quality control of two Dracocephalum species based on highperformance

    liquid chromatographic fingerprints coupled with

    tandem mass spectrometry and chemometrics [J]. J. Separaton

    Sci., 2016, 39(21):4075-4081.

    [3] 左馬怡,楊超,田倩,等. 青蘭屬植物甘青青蘭化學(xué)成分研

    究[J].云南民族大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,24(2):101-103.

    ZUO M Y, YANG C, TIAN Q, et al .. Chemical constituents of

    Dracocephalum tanguticum Maxim. of genus Dracocephalum [J].

    J. Yunnan Minzu Univ. (Nat. Sci.), 2015, 24(2):101-103.

    [4] 劉建英,劉玉梅.青蘭屬植物的化學(xué)成分及藥理作用研究進(jìn)

    展[J].食品科學(xué),2012,33(13):314-319.

    LIU J Y, LIU Y M. Research advance in chemical compositions

    and pharmacological effects of Dracocephalum [J]. Food Sci.,

    2012, 33(13):314-319.

    [5] 謝建鋒,朱林燕,孔子銘,等.唐古特青蘭總黃酮的提取及其

    體外抗氧化活性的研究[J]. 華西藥學(xué)雜志,2015,30(4):

    422-424.

    XIE J F, ZHU L Y, KONG Z M, et al .. Study on the extraction

    of total flavonoids from Dracocephalumn tanguticum and its

    activity in vitro [J]. West China J. Pharm. Sci., 2015, 30(4):

    422-424.

    [6] 李永芳,李延斌,楊梅,等.唐古特青蘭對低氧性肺動脈高壓

    大鼠ET-1和NO的影響[J].中成藥,2016,38(10):2260-2262.

    [7] YAO H, ZHOU L C, TANG L L, et al .. Protective effects of

    luteolin-7-O-glucoside against starvation-induced injury

    through upregulation of autophagy in H9c2 Cells [J]. Biosci.

    Trends., 2017, 11(5):557-564.

    [8] 胡君茹,姜華,李喜香.對藏藥甘青青蘭4種提取物抑菌作用

    的研究[J].西部中醫(yī)藥,2014,27(6):10-12.

    HU J R, JIANG H, LI X X. Study on antibacterial activity of four

    kinds of extracts from Tibetan medicine Ganqing Qinglan [J].

    Western. J. Tradition Chin. Med., 2014, 27(6):10-12.

    [9] 楊毅,田侃,倪新興,等.中藥材品質(zhì)影響因素實(shí)證研究[J].中

    藥材,2016,39(6):1251-1256.

    [10] 羅守楊,田春杰,張加凡,等.根際微生物對藥材道地性的影

    響及應(yīng)用前景展望[J].時(shí)珍國醫(yī)國藥,2022,33(4):948-950.

    [11] CARDONE L, CASTRONUOVO D, PERNIOLA M, et al ..

    Evaluation of corm origin and climatic conditions on saffron

    (Crocus sativus L.) yield and quality [J]. J. Sci. Food Agric.,

    2019, 99(13):5858-5869.

    [12] 徐錦堂,郭順星,范黎,等.天麻種子與小菇屬真菌共生萌發(fā)

    的研究[J].菌物學(xué)報(bào),2001(01):137-141.

    XU J T, GUO S X, FAN L, et al .. Symbiotic germination

    between Gasterodia elata and fungal species of Mycena [J].

    Mycosystema, 2001(01):137-141.

    [13] KIM Y S, BALARAJU K, JEON Y. Biological characteristics

    of Bacillus amyloliquefaciens AK‐0 and suppression of ginseng

    root rot caused by Cylindrocarpon destructans [J]. J. Appl.

    Microbiol., 2017, 122(1):166-179.

    [14] EGAMBERDIEVA D, LI L, LINDSTROM K, et al .. A

    synergistic interaction between salt-tolerant Pseudomonas and

    Mesorhizobium strains improves growth and symbiotic

    performance of liquorice (Glycyrrhiza uralensis Fish.) under salt

    stress [J]. Appl. Microbiol. Biotechnol., 2016, 100(6):2829-2841.

    [15] KHALEDIYAN N, WEISANY W, SCHENK P M. Arbuscular

    mycorrhizae and rhizobacteria improve growth, nutritional

    status and essential oil production in Ocimum basilicum and

    Satureja hortensis [J/OL]. Ind. Crops Prod., 2020, 160:113163

    [2022-08-10]. https://doi.org/10.1016/j.indcrop.2020.113163.

    [16] RUI H R, CHEN Y L, PARE P W, et al .. Beneficial soil

    microbe promotes seed germination, plant growth and

    photosynthesis in herbal crop Codonopsis pilosula [J]. Crop

    Pasture Sci., 2016, 67(1):91-98.

    [17] ARORA M, SAXENA P, CHOUDHARY D K, et al .. Dual

    symbiosis between Piriformospora indica and Azotobacter

    chroococcum enhances the artemisinin content in Artemisia

    annua L. [J]. World J. Microbiol. Biotechnol., 2016, 32(2):1-10.

    [18] 魏艷晨,陳吉祥,王永剛,等.荒漠植物珍珠豬毛菜根際土壤

    細(xì)菌多樣性與土壤理化性質(zhì)相關(guān)性分析[J].中國農(nóng)業(yè)科技

    導(dǎo)報(bào),2022,24(5):209-217.

    WEI Y C, CHEN J X, WANG Y G, et al .. Analysis of bacterial

    diversity in the rhizosphere soil of Salsola passerina and its

    correlation with the soil physical and chemical properties [J]. J.

    Agric Sci Technol., 2022, 24(5):209-217.

    [19] 張二豪,趙潤東,祿亞洲,等.不同種源地喜馬拉雅紫茉莉根

    際土壤理化性質(zhì)及細(xì)菌群落結(jié)構(gòu)組成分析[J].山東農(nóng)業(yè)科

    學(xué),2022,54(4): 82-89.

    ZHANG E H, ZHAO R D, LU Y Z, et al .. Soil physiochemical

    properties and microbial community structure of mirabilis

    himalaica rhizosphere from different provenances [J].

    Shandong Agric Sci., 2022, 54(4): 82-89.

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

    1999:34-108.

    [21] 李茜,徐瑞蔓,陳迪,等.不同栽培年限人參根際土壤細(xì)菌群

    落與土壤理化性質(zhì)和酶活性的相關(guān)性[J].植物營養(yǎng)與肥料

    學(xué)報(bào),2022,28(2):313-324.

    L X, XU R M, CHEN D, et al .. Correlation of bacterial

    community with soil physicochemical properties and enzyme

    activities in rhizosphere soil under different cultivation years of

    Ginseng (Panax ginseng C. A. Mey.) [J]. J. Plant Nutr. Fert.,

    2022, 28(2):313-324.

    [22] 李毳,劉怡,劉晉仙.藥用植物根際細(xì)菌群落多樣性驅(qū)動因素

    分析[J].生態(tài)環(huán)境學(xué)報(bào),2020,29(10):1988-1993.

    LI C, LIU Y, LIU J X. Analysis of driving factors of rhizosphere

    bacterial community diversityin three genuine medicine plants [J].

    Ecol. Environ. Sci., 2020, 29(10):1988-1993.

    [23] 李毳,劉怡.晉東南3種道地藥材植物根際真菌群落特性[J].

    生態(tài)環(huán)境學(xué)報(bào),2019,28(7):1388-1393.

    LI C, LIU Y. Characteristics of rhizosphere fungal community

    of three genuine medicinal plants in southeast Shanxi province

    [J]. Ecol. Environ. Sci., 2019, 28(7):1388-1393.

    [24] BURNS J H, ANACKER B L, STRAUSS S Y, et al .. Soil

    microbial community variation correlates most strongly with

    plant species identity, followed by soil chemistry, spatial

    location and plant genus [J]. AoB Plants, 2015, 7(3):207-240.

    [25] PHILIPPOT L, ANDERSSON S G E, BATTIN T J, et al .. The

    ecological coherence of high bacterial taxonomic ranks [J]. Nat.

    Rev. Microbiol., 2010, 8(7):523-529.

    [26] MENENDEZ E, CARRO L. Actinobacteria and their role as

    plant probiotics//giri b, prasad r, wu qs, varma a. biofertilizers

    for sustainable agriculture and environment [J]. Biofert. Sust.

    Agric. Environ., 2019, 55:333-351.

    [27] 宋兆齊,王莉,劉秀花,等.云南4處酸性熱泉中的變形菌門細(xì)

    菌多樣性[J].河南農(nóng)業(yè)大學(xué)學(xué)報(bào),2016,50(3):376-382.

    SONG Z Q, WANG L, LIU X Y, et al .. The diversities of

    proteobacteria in four acidic hot springs in Yunnan [J]. J.

    Henan Agric. Univ., 2016, 50(3):376-382.

    [28] 郭鳳仙,劉越,唐麗,等.藥用植物根際微生物研究現(xiàn)狀與展

    望[J].中國農(nóng)業(yè)科技導(dǎo)報(bào),2017,19(5):12-21.

    GUO F X, LIU Y, TANG L, et al .. Research status and

    prospect on rhizosphere microbiome of medicinal plants [J]. J.

    Agric. Sci. Technol., 2017, 19(5):12-21.

    [29] 李潔,林瑩,徐美玉,等.泰山白首烏根際土壤真菌多樣性分

    析[J].中國農(nóng)業(yè)科技導(dǎo)報(bào),2022,24(6):70-81.

    LI J, LIN Y, XU M Y, et al .. Sequencing analysis of fungal

    diversity in rhizosphere soil of Cynanchum bungei Decne [J]. J.

    Agric. Sci. Technol., 2022, 24(6):70-81.

    [30] YELLE D J, RALPH J, LU F, et al .. Evidence for cleavage of

    lignin by a brown rot basidiomycete [J]. Environ. Microbiol.,

    2010, 10(7):1844-1849.

    [31] 許國琪,劉怡萱,曹鵬熙,等. 青藏高原冰川棘豆(Oxytropis

    glacialis)內(nèi)生菌核心微生物組的界定及其互作網(wǎng)絡(luò)分析[J].

    微生物學(xué)通報(bào),2020,47(9):2746-2758.

    XU G Q, LIU Y X, CAO P X, et al .. Core microflora and

    endophytic interaction network of Oxytropis glacialis in Qinghai-

    Tibet plateau [J]. Microbiol. China, 2020, 47(9):2746-2758.

    [32] LIU Y H, WEI Y Y, MOHAMAD O, et al .. Diversity,

    community distribution and growth promotion activities of

    endophytes associated with halophyte Lycium ruthenicum

    Murr [J/OL]. 3 Biotech., 2019, 9(4):144 [2022-08-10]. https://

    doi.org/10.1007/s13205-019-1678-8.

    [33] SONG X, WU H, YIN Z, et al.. Endophytic bacteria isolated from

    panax ginseng improves ginsenoside accumulation in adventitious

    ginseng root culture [J/OL]. Molecules, 2017, 22(6):837 [2022-08-

    10]. https://doi.org/10.3390/molecules22060837.

    [34] 呂志堂,劉志恒.假諾卡氏菌科的分類學(xué)研究進(jìn)展[J].微生

    物學(xué)通報(bào),1999(3):210-215.

    [35] ASAF S, KHAN A L, KHAN M A, et al .. Osmoprotective

    functions conferred to soybean plants via inoculation with

    Sphingomonas sp. LK11 and exogenous trehalose [J/OL].

    Microbiol. Res., 2017:135 [2022-08-10]. https://doi.org/10.1016/

    j.micres.2017.08.009..

    [36] NIU X G, SONG L C, XIAO Y N, et al .. Drought-tolerant plant

    growth-promoting rhizobacteria associated with foxtail millet in

    a semi-arid agroecosystem and their potential in alleviating

    drought stress [J/OL]. Front. Microbiol., 2018, 8:2580 [2022-

    08-10]. https://doi.org/10.3389/fmicb.2017.02580.

    [37] DIAS M P, BASTOS M S, XAVIER V B, et al .. Plant growth

    and resistance promoted by Streptomyces spp. in tomato [J].

    Plant Physiol. Biochem., 2017, 118(9):479-493.

    [38] 金航. 絲狀真菌被孢霉生長群體感應(yīng)的初步研究[D].武漢:

    華中科技大學(xué),2016.

    JIN H. A preliminary research on the quorum-sensing of

    filamentous Mortierella [D]. Wuhan: Huazhong University of

    Science and Technology, 2016.

    [39] 于幫紅,徐慧,張傳博.一株冬蟲夏草來源真菌的鑒定及抑菌

    活性檢測[J].中國實(shí)驗(yàn)方劑學(xué)雜志,2016,22(8):36-40

    YU B H, XU H, ZHANG C B. Identification of fungi from

    Cordyceps sinensis and its antibacterial activity [J]. Chin. J.

    Exp. Tradit. Med. Form., 2016, 22(8):36-40.

    [40] 黃曉麗.毛殼菌(Chaetomium spp.)生防菌株的篩選及其生物

    防治機(jī)制研究[D].雅安:四川農(nóng)業(yè)大學(xué),2009.

    HUANG X L. Screens biocontrol strains from Chaetomium spp.

    and studies on biocontrol mechanisms of Chaetomium spp. [D].

    Yaan: Sichuan Agricultural University, 2009.

    [41] LIU X, CUI X, JI D, et al .. Luteolin-induced activation of the

    phenylpropanoid metabolic pathway contributes to quality

    maintenance and disease resistance of sweet cherry [J/OL].

    Food Chem., 2020, 342(1): 128309 [2022-08-10]. https://doi.

    org/10.1016/j.foodchem.2020.128309.

    [42] WANG X Y, ZHANG X H, YANG M X, et al.. Multi-Site

    evaluation of accumulated temperature and rainfall for maize yield

    and disease in Loess plateau [J/OL]. Agriculture, 2021, 11(4):373

    [2022-08-10]. https://doi.org/10.3390/agriculture11040373.

    (責(zé)任編輯:張冬玲)

    猜你喜歡
    微生物群落土壤理化性質(zhì)高通量測序
    微生物群落變化與妊娠結(jié)局關(guān)系的研究進(jìn)展
    川明參輪作對煙地土壤微生物群落結(jié)構(gòu)的影響
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    人參根際真菌群落多樣性及組成的變化
    LncRNAs作為miRNA的靶模擬物調(diào)節(jié)miRNA
    土壤調(diào)理劑對植煙土壤微生物群落的影響
    馬纓丹入侵對根際和非根際土壤理化性質(zhì)的影響
    林下養(yǎng)殖蚯蚓對鹽堿地土壤和苗木生長的影響
    稻蟹共作模式下稻蟹產(chǎn)出與土壤理化性質(zhì)的研究
    內(nèi)蒙古典型草原土壤有機(jī)碳與土壤理化性質(zhì)的關(guān)系
    国产 一区精品| 色视频www国产| 国产精品久久久久久精品电影小说 | 色吧在线观看| 一级毛片久久久久久久久女| 亚洲高清免费不卡视频| 少妇 在线观看| 亚洲成人一二三区av| 如何舔出高潮| tube8黄色片| 亚洲国产精品999| 国产在线免费精品| 亚洲av欧美aⅴ国产| 久久久亚洲精品成人影院| 91aial.com中文字幕在线观看| 久久99热这里只有精品18| 久久久久久久久久久免费av| 十八禁网站网址无遮挡 | 欧美zozozo另类| 插阴视频在线观看视频| 中文乱码字字幕精品一区二区三区| 99热这里只有是精品在线观看| 中文字幕亚洲精品专区| 欧美成人精品欧美一级黄| 亚洲不卡免费看| 久久久a久久爽久久v久久| 久久久精品免费免费高清| 99国产精品免费福利视频| 老司机影院毛片| 亚洲欧美精品自产自拍| 精华霜和精华液先用哪个| 色吧在线观看| 国产免费视频播放在线视频| 国产视频内射| 高清视频免费观看一区二区| 国产成人精品久久久久久| 免费人妻精品一区二区三区视频| 国产成人91sexporn| 多毛熟女@视频| 久久青草综合色| 亚洲精品日本国产第一区| 成人二区视频| 国产一区有黄有色的免费视频| 97超碰精品成人国产| 三级经典国产精品| av在线观看视频网站免费| 九九久久精品国产亚洲av麻豆| 精品午夜福利在线看| 日韩强制内射视频| 欧美xxxx黑人xx丫x性爽| 18禁裸乳无遮挡免费网站照片| 丝袜脚勾引网站| 成年免费大片在线观看| 少妇人妻精品综合一区二区| 精品人妻一区二区三区麻豆| 中国美白少妇内射xxxbb| 看十八女毛片水多多多| 欧美精品人与动牲交sv欧美| 26uuu在线亚洲综合色| 成人特级av手机在线观看| 高清不卡的av网站| 春色校园在线视频观看| 欧美日韩国产mv在线观看视频 | 80岁老熟妇乱子伦牲交| 久久国产精品大桥未久av | 欧美xxxx性猛交bbbb| 嘟嘟电影网在线观看| 丝瓜视频免费看黄片| 日韩欧美一区视频在线观看 | 夜夜爽夜夜爽视频| 久久婷婷青草| 高清午夜精品一区二区三区| 少妇丰满av| 亚洲一级一片aⅴ在线观看| 深爱激情五月婷婷| 日韩免费高清中文字幕av| 人妻系列 视频| 极品少妇高潮喷水抽搐| 一级毛片aaaaaa免费看小| 国产成人精品福利久久| 国产黄频视频在线观看| 一级毛片电影观看| 麻豆成人午夜福利视频| 亚洲精品一区蜜桃| 国产一区二区三区综合在线观看 | 少妇高潮的动态图| 欧美丝袜亚洲另类| 国产爽快片一区二区三区| 国产色爽女视频免费观看| 国产伦在线观看视频一区| 日本-黄色视频高清免费观看| 国产 一区精品| 免费久久久久久久精品成人欧美视频 | 1000部很黄的大片| 一本—道久久a久久精品蜜桃钙片| 一区二区av电影网| 搡老乐熟女国产| 在线观看免费高清a一片| 少妇人妻精品综合一区二区| 久久精品久久久久久久性| tube8黄色片| 亚洲久久久国产精品| 国产成人一区二区在线| av在线观看视频网站免费| 五月玫瑰六月丁香| 能在线免费看毛片的网站| 久久人人爽人人片av| 99热这里只有是精品50| 制服丝袜香蕉在线| 亚洲精品色激情综合| 亚洲精品久久午夜乱码| 久久久久久久久久久丰满| 亚洲不卡免费看| av在线老鸭窝| 国产有黄有色有爽视频| 日本vs欧美在线观看视频 | 国产精品秋霞免费鲁丝片| 国产精品女同一区二区软件| 国产高潮美女av| 99视频精品全部免费 在线| 日本午夜av视频| 国产伦精品一区二区三区视频9| 嫩草影院新地址| 国产成人aa在线观看| 97超碰精品成人国产| 五月伊人婷婷丁香| 亚洲四区av| 亚洲丝袜综合中文字幕| 精华霜和精华液先用哪个| 嫩草影院入口| h视频一区二区三区| freevideosex欧美| 国产精品久久久久久av不卡| 日韩视频在线欧美| 国产真实伦视频高清在线观看| 亚洲色图综合在线观看| 99热这里只有是精品50| 国产午夜精品久久久久久一区二区三区| a级毛片免费高清观看在线播放| 成人特级av手机在线观看| 亚洲精品日本国产第一区| 国产 一区 欧美 日韩| 久久女婷五月综合色啪小说| 亚洲av成人精品一二三区| 久久人人爽av亚洲精品天堂 | 九九爱精品视频在线观看| 久久亚洲国产成人精品v| 在线天堂最新版资源| 亚洲精品日韩av片在线观看| 成人特级av手机在线观看| 国产在视频线精品| 国产欧美日韩一区二区三区在线 | 高清午夜精品一区二区三区| 久久人妻熟女aⅴ| 亚洲怡红院男人天堂| 青青草视频在线视频观看| 乱系列少妇在线播放| 国产高清三级在线| 精品国产露脸久久av麻豆| 成人国产麻豆网| 成人18禁高潮啪啪吃奶动态图 | 久久久久久久久久久免费av| 一区二区三区精品91| 在线观看三级黄色| 欧美精品亚洲一区二区| 欧美97在线视频| 大陆偷拍与自拍| 亚洲av成人精品一二三区| 人人妻人人添人人爽欧美一区卜 | 韩国高清视频一区二区三区| 熟女av电影| 熟女电影av网| 日韩大片免费观看网站| 国产精品嫩草影院av在线观看| 高清在线视频一区二区三区| 熟女电影av网| 一边亲一边摸免费视频| 久久久久精品性色| 精品久久久噜噜| 久久久久久久久久成人| 精品久久久久久久久亚洲| 草草在线视频免费看| 日韩成人av中文字幕在线观看| 能在线免费看毛片的网站| 亚洲欧美一区二区三区国产| 欧美高清成人免费视频www| 看十八女毛片水多多多| 亚洲综合色惰| 性色avwww在线观看| 美女cb高潮喷水在线观看| 久久精品国产亚洲av天美| 国产在线视频一区二区| 一区二区三区精品91| 狂野欧美激情性xxxx在线观看| 国产成人a∨麻豆精品| 免费观看的影片在线观看| 涩涩av久久男人的天堂| 亚洲电影在线观看av| 高清日韩中文字幕在线| 大香蕉97超碰在线| 欧美97在线视频| 欧美激情极品国产一区二区三区 | av福利片在线观看| 亚洲成人中文字幕在线播放| 国产精品国产三级国产av玫瑰| 丰满少妇做爰视频| 老女人水多毛片| 午夜福利在线观看免费完整高清在| 色婷婷av一区二区三区视频| 能在线免费看毛片的网站| 伊人久久国产一区二区| 免费av不卡在线播放| 久久精品久久久久久噜噜老黄| 久久久久久久精品精品| 久久久久久久大尺度免费视频| 香蕉精品网在线| 91精品一卡2卡3卡4卡| 看十八女毛片水多多多| 最近最新中文字幕免费大全7| 国产精品一及| 搡女人真爽免费视频火全软件| 黄片无遮挡物在线观看| 亚洲欧美日韩另类电影网站 | av不卡在线播放| 午夜老司机福利剧场| 青春草国产在线视频| videos熟女内射| 午夜免费男女啪啪视频观看| 欧美最新免费一区二区三区| 2022亚洲国产成人精品| 亚洲国产精品专区欧美| 国产亚洲最大av| 免费大片18禁| 午夜福利影视在线免费观看| 性色av一级| 精品人妻一区二区三区麻豆| 国产淫语在线视频| 男女边摸边吃奶| 少妇的逼好多水| 久久久久网色| 国产淫语在线视频| 欧美精品人与动牲交sv欧美| 亚洲国产欧美人成| 下体分泌物呈黄色| 国产成人一区二区在线| 亚洲av.av天堂| 狂野欧美白嫩少妇大欣赏| 国产成人精品福利久久| 草草在线视频免费看| 插阴视频在线观看视频| 亚洲人成网站高清观看| av网站免费在线观看视频| 视频中文字幕在线观看| 国产日韩欧美亚洲二区| 久久久亚洲精品成人影院| 色视频在线一区二区三区| 直男gayav资源| 国产精品一区二区性色av| 日韩强制内射视频| 18禁动态无遮挡网站| h日本视频在线播放| 26uuu在线亚洲综合色| 亚洲精品一区蜜桃| 麻豆成人av视频| 亚洲av男天堂| 下体分泌物呈黄色| 大香蕉久久网| 亚洲国产欧美在线一区| 亚洲精品456在线播放app| 久久久久久久久久久免费av| 成人国产av品久久久| 这个男人来自地球电影免费观看 | 久久精品国产亚洲av涩爱| 久久国产乱子免费精品| 97热精品久久久久久| 夜夜爽夜夜爽视频| av在线观看视频网站免费| 舔av片在线| 中文字幕亚洲精品专区| 两个人的视频大全免费| 舔av片在线| av在线app专区| 日韩欧美 国产精品| 国产精品欧美亚洲77777| 亚洲人与动物交配视频| 91狼人影院| 国产精品成人在线| 亚洲人与动物交配视频| 亚洲av欧美aⅴ国产| 夜夜骑夜夜射夜夜干| 网址你懂的国产日韩在线| 一级毛片久久久久久久久女| 国产精品不卡视频一区二区| 97在线视频观看| 一个人免费看片子| 色婷婷av一区二区三区视频| 直男gayav资源| 亚洲精品日韩av片在线观看| 国产精品久久久久久精品古装| 91久久精品电影网| 99久国产av精品国产电影| av福利片在线观看| 国产免费一区二区三区四区乱码| 亚洲色图av天堂| 黄色怎么调成土黄色| kizo精华| 日韩强制内射视频| 免费看不卡的av| 又粗又硬又长又爽又黄的视频| 又爽又黄a免费视频| 国产一区二区在线观看日韩| 亚洲国产成人一精品久久久| 少妇被粗大猛烈的视频| 在线亚洲精品国产二区图片欧美 | 肉色欧美久久久久久久蜜桃| 我要看日韩黄色一级片| 熟妇人妻不卡中文字幕| 久久久久网色| 精品亚洲成国产av| 肉色欧美久久久久久久蜜桃| 色婷婷av一区二区三区视频| 亚洲美女搞黄在线观看| 国产伦精品一区二区三区视频9| 久热久热在线精品观看| 美女主播在线视频| 91精品国产九色| 国产极品天堂在线| 51国产日韩欧美| 人妻夜夜爽99麻豆av| 午夜老司机福利剧场| 一个人看的www免费观看视频| 久久久久久久精品精品| 亚洲丝袜综合中文字幕| 高清av免费在线| 少妇精品久久久久久久| 一本—道久久a久久精品蜜桃钙片| 伊人久久国产一区二区| av国产免费在线观看| 欧美3d第一页| freevideosex欧美| 国产片特级美女逼逼视频| 在线天堂最新版资源| av卡一久久| 人妻制服诱惑在线中文字幕| 在线观看免费高清a一片| 欧美激情极品国产一区二区三区 | 国产成人午夜福利电影在线观看| 国产精品一区二区在线不卡| 尾随美女入室| 嫩草影院新地址| 国产探花极品一区二区| 免费观看无遮挡的男女| 人妻制服诱惑在线中文字幕| 国产成人a∨麻豆精品| www.av在线官网国产| 日本av手机在线免费观看| 久久亚洲国产成人精品v| 我要看日韩黄色一级片| 亚洲国产精品999| 91久久精品国产一区二区三区| 精品一区二区三卡| 国产精品欧美亚洲77777| 国产伦精品一区二区三区四那| 七月丁香在线播放| 日本欧美视频一区| 久久久国产一区二区| 青春草视频在线免费观看| 久久久久视频综合| 国产色婷婷99| 久久久精品94久久精品| 午夜视频国产福利| 日韩强制内射视频| 亚洲无线观看免费| 黄色怎么调成土黄色| 亚洲最大成人中文| 亚洲高清免费不卡视频| 91狼人影院| 国产成人免费观看mmmm| 精品亚洲成国产av| 精品午夜福利在线看| 夜夜爽夜夜爽视频| 色哟哟·www| av视频免费观看在线观看| 黑人猛操日本美女一级片| 男女无遮挡免费网站观看| 亚洲精品乱久久久久久| 国产淫语在线视频| 一级毛片电影观看| 国产高清三级在线| 午夜激情久久久久久久| 国产又色又爽无遮挡免| 国国产精品蜜臀av免费| 亚洲性久久影院| 国产精品一区二区在线不卡| 高清视频免费观看一区二区| 亚洲,一卡二卡三卡| 亚洲中文av在线| 一二三四中文在线观看免费高清| 最近中文字幕2019免费版| 美女cb高潮喷水在线观看| 大话2 男鬼变身卡| 亚洲欧美日韩另类电影网站 | av女优亚洲男人天堂| 在线观看av片永久免费下载| 少妇 在线观看| 国产精品偷伦视频观看了| 国产一级毛片在线| 亚洲av成人精品一二三区| 尤物成人国产欧美一区二区三区| 国产成人免费观看mmmm| 联通29元200g的流量卡| 精品一区在线观看国产| 亚洲国产精品999| av在线播放精品| 色综合色国产| 精品久久久久久久末码| 日本-黄色视频高清免费观看| 亚洲欧美清纯卡通| 交换朋友夫妻互换小说| 狂野欧美激情性bbbbbb| 欧美精品一区二区免费开放| 欧美成人午夜免费资源| 一本色道久久久久久精品综合| 欧美97在线视频| 在线观看三级黄色| 97超碰精品成人国产| 偷拍熟女少妇极品色| 伦精品一区二区三区| 亚洲一级一片aⅴ在线观看| 午夜激情福利司机影院| 美女脱内裤让男人舔精品视频| 天堂8中文在线网| 欧美变态另类bdsm刘玥| 久久99热这里只有精品18| 麻豆精品久久久久久蜜桃| av免费观看日本| 久久精品国产鲁丝片午夜精品| 国产一区亚洲一区在线观看| 免费看不卡的av| 一级毛片久久久久久久久女| 波野结衣二区三区在线| 国产片特级美女逼逼视频| 这个男人来自地球电影免费观看 | 视频中文字幕在线观看| 欧美最新免费一区二区三区| 亚洲国产色片| 亚洲精品,欧美精品| 狠狠精品人妻久久久久久综合| 久久久成人免费电影| 又黄又爽又刺激的免费视频.| 精品午夜福利在线看| a级一级毛片免费在线观看| 久久影院123| 日韩,欧美,国产一区二区三区| 男男h啪啪无遮挡| 你懂的网址亚洲精品在线观看| 在线 av 中文字幕| 久久久亚洲精品成人影院| 婷婷色av中文字幕| 久久99热这里只有精品18| 国产免费又黄又爽又色| 久久人人爽av亚洲精品天堂 | 久久久久网色| 亚洲美女黄色视频免费看| 狠狠精品人妻久久久久久综合| 亚洲一级一片aⅴ在线观看| 国产色爽女视频免费观看| 最后的刺客免费高清国语| 免费看日本二区| 久久99蜜桃精品久久| 2021少妇久久久久久久久久久| 成人亚洲精品一区在线观看 | 免费在线观看成人毛片| 日韩欧美 国产精品| 国产老妇伦熟女老妇高清| 日本-黄色视频高清免费观看| 国模一区二区三区四区视频| 女性生殖器流出的白浆| av国产精品久久久久影院| 91狼人影院| 噜噜噜噜噜久久久久久91| 成人午夜精彩视频在线观看| 久久97久久精品| 人妻制服诱惑在线中文字幕| 婷婷色综合www| 在线观看av片永久免费下载| 中文字幕亚洲精品专区| 国产男人的电影天堂91| 欧美激情国产日韩精品一区| av国产久精品久网站免费入址| 寂寞人妻少妇视频99o| 天堂俺去俺来也www色官网| 一区二区av电影网| 久久精品熟女亚洲av麻豆精品| 视频中文字幕在线观看| 免费观看av网站的网址| 性色avwww在线观看| 老司机影院毛片| av国产久精品久网站免费入址| 欧美激情极品国产一区二区三区 | 熟妇人妻不卡中文字幕| 亚洲综合色惰| 欧美高清成人免费视频www| 美女福利国产在线 | 最后的刺客免费高清国语| 亚洲激情五月婷婷啪啪| 亚洲精品乱码久久久v下载方式| 2018国产大陆天天弄谢| 欧美激情国产日韩精品一区| 成人无遮挡网站| 国产精品.久久久| 偷拍熟女少妇极品色| 中国美白少妇内射xxxbb| 一级a做视频免费观看| 免费看日本二区| 久久热精品热| 久久青草综合色| 国产伦在线观看视频一区| 日本爱情动作片www.在线观看| 不卡视频在线观看欧美| 欧美另类一区| 大又大粗又爽又黄少妇毛片口| 又粗又硬又长又爽又黄的视频| 国产精品99久久久久久久久| 美女福利国产在线 | 免费观看无遮挡的男女| 1000部很黄的大片| 亚洲av不卡在线观看| 一级片'在线观看视频| 日本黄色日本黄色录像| 亚洲国产色片| 国内揄拍国产精品人妻在线| av在线app专区| 青春草亚洲视频在线观看| 欧美日韩综合久久久久久| 久久人人爽人人片av| 日本vs欧美在线观看视频 | 国产日韩欧美在线精品| 国产精品女同一区二区软件| 日本-黄色视频高清免费观看| 免费看av在线观看网站| 亚洲国产欧美在线一区| 久久久久国产网址| 精品人妻熟女av久视频| 99久久精品热视频| 熟女av电影| 成人无遮挡网站| 免费人妻精品一区二区三区视频| 国产色爽女视频免费观看| 欧美日韩视频高清一区二区三区二| 天堂中文最新版在线下载| 美女高潮的动态| 视频区图区小说| 秋霞伦理黄片| 男男h啪啪无遮挡| 亚洲一区二区三区欧美精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品国产色婷婷电影| 成人亚洲欧美一区二区av| 国产中年淑女户外野战色| 一级毛片 在线播放| 欧美+日韩+精品| av不卡在线播放| 黄片wwwwww| 小蜜桃在线观看免费完整版高清| 精品久久久精品久久久| 深夜a级毛片| 一级爰片在线观看| 久久6这里有精品| 成人毛片a级毛片在线播放| 丝瓜视频免费看黄片| 久久影院123| 午夜福利在线观看免费完整高清在| 欧美成人精品欧美一级黄| 国产精品三级大全| 成人亚洲欧美一区二区av| 国产精品国产三级国产专区5o| 免费不卡的大黄色大毛片视频在线观看| 日韩国内少妇激情av| 乱系列少妇在线播放| 色婷婷av一区二区三区视频| 日韩欧美一区视频在线观看 | 能在线免费看毛片的网站| 在线看a的网站| 在线免费观看不下载黄p国产| 一区在线观看完整版| 伊人久久国产一区二区| 日本猛色少妇xxxxx猛交久久| 另类亚洲欧美激情| 舔av片在线| 五月开心婷婷网| 99久久精品国产国产毛片| 天天躁夜夜躁狠狠久久av| 成年免费大片在线观看| 最黄视频免费看| 国产亚洲av片在线观看秒播厂| 五月玫瑰六月丁香| 久久久欧美国产精品| 久久99热这里只频精品6学生| 人人妻人人澡人人爽人人夜夜| 99久国产av精品国产电影| 亚洲av国产av综合av卡| 国产欧美日韩一区二区三区在线 | 亚洲久久久国产精品| 亚洲成人手机| 91久久精品电影网| 国产一级毛片在线| 成人国产av品久久久| av在线蜜桃| av线在线观看网站| 天堂8中文在线网| 国产永久视频网站| av在线老鸭窝| 国精品久久久久久国模美| 汤姆久久久久久久影院中文字幕|