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

    甘肅野生羊肚菌根際土壤真菌群落與環(huán)境因子相互關(guān)系

    2022-04-11 11:29:16趙玉卉楊阿麗魏甲乾冉永紅陳志榮路等學(xué)
    微生物學(xué)雜志 2022年1期
    關(guān)鍵詞:子囊菌根羊肚

    趙玉卉, 郭 瑞, 楊阿麗, 秦 鵬, 金 輝,魏甲乾, 冉永紅, 陳志榮, 路等學(xué)*

    (1.甘肅省科學(xué)院 生物研究所,甘肅 蘭州 730000;2.中國科學(xué)院 蘭州化學(xué)物理研究所 中國科學(xué)院西北特色植物資源化學(xué)重點實驗室,甘肅 蘭州 730000;3.定西市農(nóng)業(yè)科學(xué)研究院,甘肅 定西 743000;4.溫州醫(yī)科大學(xué) 生物醫(yī)學(xué)工程學(xué)院,浙江 溫州 325000)

    羊肚菌(Morchellaspp.)屬子囊菌門(Ascomycota)盤菌綱(Pezizomycetes)盤菌目(Pezizales)羊肚菌科(Moechellaceae)羊肚菌屬(Morchella),因其外形似羊肚而得名[1]。羊肚菌富含蛋白質(zhì)、大量人體必需氨基酸,多種維生素和高含量的鐵、鋅以及多種礦質(zhì)元素[2-3],具有抗射線、抗腫瘤、抗氧化、抗疲勞作用以及降血壓、降血脂、調(diào)節(jié)免疫力的功效[4]。羊肚菌經(jīng)濟(jì)價值和科研價值高,人工栽培一直是研究的熱點。2012年后,我國羊肚菌大田種植技術(shù)取得重要進(jìn)展,在全國范圍內(nèi)得到推廣,種植面積陡增。經(jīng)過10多年的發(fā)展,出菇不穩(wěn)定、不出菇、產(chǎn)量不高等問題依然普遍存在,且有新的問題不斷出現(xiàn),如病蟲害[5-8]和連作障礙[9-10]等問題。究其原因是羊肚菌基礎(chǔ)理論的研究不夠,羊肚菌的生活史、菌核的機(jī)理、分生孢子的作用等有待進(jìn)一步研究;羊肚菌是一種土生菌,不覆土羊肚菌無法出菇,可見土壤中含有特定刺激其結(jié)實的有益因子[11]。土壤中除了營養(yǎng)物質(zhì),豐富的土壤微生物對其生長有一定的影響。已有研究者對羊肚菌內(nèi)生細(xì)菌、根際土壤細(xì)菌多樣性等進(jìn)行了研究[12-15],結(jié)果表明羊肚菌根際土壤的細(xì)菌多樣性指數(shù)和豐度比非根際土壤的高。關(guān)于羊肚菌根際土壤真菌群落方面的研究較少,何培新等[16]對羊肚菌內(nèi)生真菌進(jìn)行了研究,結(jié)果表明粗柄羊肚菌的內(nèi)生真菌中,毛殼屬是3個發(fā)生地的共有內(nèi)生真菌,其他內(nèi)生真菌均為發(fā)生地特異種類。有研究顯示真菌群落中,高豐度的孔球孢屬(Gilmaniella)似乎與室內(nèi)紅褐羊肚菌(Morchellarufobrunnea)的結(jié)實有關(guān),而優(yōu)勢菌群頭束霉屬(Cephalotrichum)與不出菇有關(guān)[17]。土壤理化性質(zhì)與土壤酶是影響土壤微生物群落數(shù)量、結(jié)構(gòu)的主要因素之一[18],但羊肚菌根際土壤理化性質(zhì)及酶活性與真菌群落相互關(guān)系的研究未見相關(guān)報道。本研究采用高通量測序技術(shù),對甘肅不同地區(qū)羊肚菌根際土壤真菌群落組成、多樣性等結(jié)構(gòu)特征進(jìn)行系統(tǒng)研究,并結(jié)合土壤理化性質(zhì)和酶活性,分析羊肚菌根際微生態(tài)各要素之間的相互關(guān)系,以期揭示根際微生態(tài)系統(tǒng)對羊肚菌生態(tài)適應(yīng)的內(nèi)在機(jī)制,并為羊肚菌的人工種植及生態(tài)功能研究提供參考。

    1 材料與方法

    1.1 材料

    1.1.1 樣品來源 根據(jù)羊肚菌在甘肅省不同區(qū)域的分布情況,設(shè)置9個采樣點,在相同經(jīng)緯度分別設(shè)置DS、DC采樣點和LC、LK采樣點。原因是DS、LC采樣點的羊肚菌屬黑色系,而DC、LK采樣點的羊肚菌屬黃色系,其他采樣點羊肚菌均為黑色系,具體信息見表1。每個采樣點內(nèi)設(shè)置3個10 m×10 m樣方,作為平行重復(fù),每個樣方內(nèi)選取長勢良好的羊肚菌,用干凈的鐵鏟采集子實體及其半徑5 cm、深度0~5 cm的土壤,從距離羊肚菌根際最近的位置用無菌毛刷收集土壤,混合土樣,低溫帶回實驗室,放入-80 ℃冰箱冷凍備用。

    1.1.2 試劑與儀器 土壤DNA提取試劑盒E.Z.N.ATM Mag-Bind Soil DNA Kit(OMEGA公司);土壤酶活試劑盒(蘇州科銘生物技術(shù)有限公司)。pH計(PHS-3C,上海越平);紫外分光光度計(TU-1950,北京普析通用儀器有限責(zé)任公司);高頻紅外碳硫測定儀(COREY-200,德陽市科瑞儀器設(shè)備廠); 自動定氮儀(K12,上海晟聲自動化分析儀器有限公司);電感耦合等離子體直讀光譜儀(ICAP-7400,賽默飛世爾科技有限公司); 原子熒光光譜儀(AFS-2000,北京海光儀器有限公司);離子色譜儀(Aquion,賽默飛世爾科技有限公司); 火焰原子吸收儀(Z-2000,日本日立公司)。

    表1 野生羊肚菌采樣地信息表

    1.2 方法

    1.2.1 土壤理化性質(zhì)及酶活性測定 土壤pH值采用酸度計測定[19],以去除CO2的純水為浸提劑,V水∶V土=2.5∶1;土壤溫度使用土壤溫度計測定;采用重鉻酸鹽法[20]測定土壤有機(jī)質(zhì)含量;采用重鉻酸鉀硫酸消化法[21]測定土壤全氮含量;全碳使用高頻紅外碳硫測定儀測定;銨態(tài)氮采用自動定氮儀測定;硝態(tài)氮采用離子色譜儀測定;有效磷采用紫外分光光度計測定;速效鉀采用火焰原子吸收儀測定;土壤微量元素(Fe、Ca、Mg)含量采用電感耦合等離子體直讀光譜儀測定;Se測定采用原子熒光光譜儀測定。土壤含水量采用烘干法測定,用0.1 g 精度的天平稱取土樣重量,記作土樣濕質(zhì)量重M,105 ℃烘干土樣6~8 h至恒重,記作土樣干質(zhì)量重Ms。土壤含水量(%)=((烘干前鋁盒及土樣質(zhì)量-烘干后鋁盒及土樣質(zhì)量)/(烘干后鋁盒及土樣質(zhì)量-烘干空鋁盒質(zhì)量))×100%。土壤多酚氧化酶活性采用沒食子素比色法測定;土壤過氧化物酶活性采用鄰苯三酚比色法測定;土壤脲酶活性采用靛酚藍(lán)比色法測定;脫氫酶活性采用三苯基四氮唑氯化物(TTC) 比色法測定;土壤堿性磷酸酶活性采用磷酸苯二鈉比色法測定;芳基硫酸酯酶活性采用對-硝基苯酚比色法測定;谷氨酰胺酶活性采用氨比色法測定;β-葡萄糖苷酶活性采用對-硝基苯酚比色法測定。

    1.2.2 土壤樣品總DNA提取與Mi-Seq測序 對采集到的羊肚菌根際土壤樣品,利用土壤DNA提取試劑盒E.Z.N.ATM Mag-Bind Soil DNA Kit對土壤總基因組DNA進(jìn)行提取,采用瓊脂糖凝膠電泳檢測DNA的純度和濃度。擴(kuò)增區(qū)域為ITS1-2,通用引物ITS1F:CCCTACACGACGCTCTTCCG-ATCTN (barcode) CTTGGTCATTTAGAGGAAGTA-

    A,ITS2-R引物:GTGACTGGAGTTCCTTGGCACC-CGAGAATTCCAGCTGCGTTCTTCATCGATGC。擴(kuò)增體系:15 μL×Taq2 master Mix,1 μL Bar- PCR primer F(10 μmol/L), 1 μL primerR(10 μmol/L),Genomic DNA 10~20 ng,加H2O補(bǔ)足至30 μL。擴(kuò)增條件為94 ℃ 3 min,94 ℃ 30 s、45 ℃ 20 s、65 ℃ 30 s 5個循環(huán),94 ℃ 20 s、55 ℃ 20 s、72 ℃ 30 s 20個循環(huán),72 ℃ 5 min。第二輪擴(kuò)增引入Illuminn橋式PCR兼容引物。擴(kuò)增體系:15 μL 2×Taqmaster Mix,1 μL Bar-PCR primerF(10 μmol/L), 1 μL primerR(10 μmol/L),PCR products(上一輪)20 ng,加H2O補(bǔ)足至30 μL。擴(kuò)增條件為95 ℃ 3 min,94 ℃ 20 s、55 ℃ 20 s,72 ℃ 30 s 5個循環(huán),72 ℃ 5 min。PCR擴(kuò)增產(chǎn)物經(jīng)瓊脂糖凝膠電泳檢測后回收純化,應(yīng)用Illumina MiSeq 2×300 bp平臺進(jìn)行宏基因組測序。測序服務(wù)和生物信息分析工作由上海生工生物工程有限公司完成。

    1.2.3 生物信息學(xué)分析 Illumina Mi-seqTM得到的原始圖像數(shù)據(jù)文件經(jīng)CASAVA堿基識別 (Base Calling) 分析轉(zhuǎn)化為原始測序序列(Sequenced Reads),稱之為 Raw Data或Raw Reads。Miseq測序序列中含有barcode序列,以及測序時加入的引物和接頭序列,首先需要去除引物接頭序列,再根據(jù)PE reads之間的overlap關(guān)系,將成對的reads 拼接 (merge) 成一條序列,然后按照barcode標(biāo)簽序列識別并區(qū)分樣品得到各樣本數(shù)據(jù),最后對各樣本數(shù)據(jù)的質(zhì)量進(jìn)行質(zhì)控過濾,得到各樣本有效數(shù)據(jù)(clean reads)。用Usearch 5.2.236軟件對所有樣品的全部clean reads序列進(jìn)行OTUs聚類(97%的一致性)。用Qiime 中的 BLAST方法與 RDP數(shù)據(jù)庫進(jìn)行物種注釋分析,OTUs 的豐度初步表明樣品的物種豐度。獲得的數(shù)據(jù)提交GenBank,登錄號:DS(SRR11905350、SRR11905361、SRR11905362), DC(SRR11905340、SRR11905341、SRR11905342), TB (SRR11905337、SRR11905338、SRR11905339), LC(SRR11905336、SRR11905359、SRR11905360), LK(SRR11905356、SRR11905357、 SRR11905358), ZQ(SRR11905353、SRR11905354、SRR11905355), DB(SRR11905349、SRR11905351、SRR11905352), LQ(SRR11905346、SRR11905347、SRR11905348), LX(SRR11905343、SRR11905344、SRR11905345)。

    1.2.4 數(shù)據(jù)分析 采用QIIME計算的指數(shù),利用軟件R進(jìn)行稀釋度曲線繪制;利用Mothur 軟件進(jìn)行單個樣品的α多樣性分析(Alpha Diversity),

    基于OTU的結(jié)果,對Shannon、Chao1、ACE指數(shù)進(jìn)行計算,用以進(jìn)行生物多樣性分析。用RDP classifier分析方法,基于Bergey′s taxonomy,采用Na?ve Bayesian assignment算法對每條序列在不同層級水平上計算其分配到此rank中的概率值。用Qiime 中的 BLAST方法與 Unite數(shù)據(jù)庫進(jìn)行物種注釋分析,各組樣品在不同水平上的分類比較柱形圖、單個樣品的群落分布柱形圖使用R軟件獲得。

    2 結(jié)果與分析

    2.1 野生羊肚菌根際土壤理化性質(zhì)及酶活性

    由表2可知,ZQ、DB、LQ、LX采樣點羊肚菌根際土壤呈酸性,DS、DC、TB、LC、LK采樣點羊肚菌根際土壤呈堿性,pH范圍6.46~8.19。大量元素隨海拔升高有逐漸增加的趨勢,DB采樣點含量最高,DC采樣點含量最低,各采樣點微量元素變化各異。9個采樣點中,土壤多酚氧化酶(PPO)、過氧化物酶(POD)、脲酶(UE)、芳基硫酸酯酶(ASF)隨海拔升高而降低;而β-葡萄糖苷酶(β-GC)、脫氫酶(DHA)、堿性磷酸酶(AKP)隨著海拔升高而升高;β-葡萄糖苷酶的酶活性在DB采樣點顯著高于其他采樣點;谷氨酰胺酶(GLS)在各采樣點變化不大。

    表2 野生羊肚菌采樣點根際土壤理化性質(zhì)及酶活性

    續(xù)表2

    2.2 野生羊肚菌根際土壤真菌群落豐度

    9個羊肚菌根際土壤樣品,以大于97%相似性聚類得到21 846個真菌OTU,統(tǒng)計結(jié)果見表3。各樣品測序覆蓋度在98.92%以上,樣品稀釋曲線顯示,隨著測序樣品數(shù)目的增加,9個樣品的OTUs稀釋曲線趨于平緩,表明本實驗測序的數(shù)據(jù)量漸進(jìn)合理,能較全面地反映測序樣品的真菌群落組成,數(shù)據(jù)量的增多對于發(fā)現(xiàn)新的OTU數(shù)目的貢獻(xiàn)變小(圖1)。

    圖1 野生羊肚菌根際土壤真菌OTU稀疏度曲線圖Fig.1 Curve of fungi OTU thinning in rhizosphere soil of wild morels

    表3 野生羊肚菌根際土壤測序結(jié)果及Alpha多樣性分析

    2.3 不同采樣點真菌多樣性及相關(guān)性分析

    Shannon指數(shù)反應(yīng)群落的多樣性受樣品群落中物種豐富度和物種均勻度的影響,Chao1指數(shù)和ACE指數(shù)反映樣品中群落的豐富度。9個采樣點羊肚菌根際土壤真菌Alpha分析結(jié)果見表3,Shannon指數(shù)分析結(jié)果為DB>LC>DS>TB>LX >ZQ>DC>LQ>LK; Chao1指數(shù)分析結(jié)果為LK>ZQ>LC>DB>TB>DS>LX>DC>LQ;ACE指數(shù)分析結(jié)果為ZQ>LC>DB>TB>DS>LX>LQ>DC>LK;DB、LC采樣點羊肚菌根際土壤真菌多樣性最高,Shannon指數(shù)分別為5.01和4.61; LK、ZQ 、LC采樣點根際土壤真菌豐富度最高。

    在97%相似度下得到每個樣品的OTU個數(shù),利用Python繪制韋恩圖,展示出多樣品間相同的和特有的OTU數(shù)以及重疊情況,如圖2所示,9個采樣點共有的OTU數(shù)32個,特有的OTU數(shù)較多。

    2.4 野生羊肚菌根際土壤真菌群落組成分析

    不同采樣點羊肚菌根際土壤在門水平的物種組成結(jié)果見圖3,9個采樣點中共檢測到18個門的真菌。采樣點DS、TB、LC、LK、ZQ、DB、LQ、LX優(yōu)勢菌門為子囊菌門,DC采樣點優(yōu)勢菌門為擔(dān)子菌門,其次是子囊菌門;除DC采樣點,其他采樣點子囊菌門的占比均在45%以上,LK、LQ子囊菌門的占比達(dá)到96%以上;未分類真菌在LK采樣點占比最低為0.38%,在LC采樣點占比最高為23.38%;在所有采樣點中擔(dān)子菌和子囊菌的占比最高的是LQ采樣點,為98.13%,最低的是ZQ采樣點,為60.94%。

    圖2 不同采樣點真菌多樣性的相關(guān)性分析Fig.2 The similarity analysis of fungal diversity among samples

    圖3 野生羊肚菌根際土壤真菌在門水平上的分類比較Fig.3 Taxonomic comparison of fungal in rhizosphere soil of wild morels at phylum level

    不同采樣點羊肚菌根際土壤在綱水平的物種組成結(jié)果見圖4,9個采樣點中共檢測到55個綱的真菌。DS、DC、TB、LC采樣點共有的優(yōu)勢菌綱是座囊菌綱(Dothideomycetes)和糞殼菌綱(Sordariomycetes);DC、ZQ、DB、LQ、LX采樣點的共有優(yōu)勢菌綱是傘菌綱(Agarico-mycetes);LK、LQ采樣點共有的優(yōu)勢菌綱是占比很高的盤菌綱(Pezizomycetes);DB、LQ、LX、LC采樣點共有的優(yōu)勢菌綱是錘舌菌綱(Leotiomycetes);ZQ采樣點獨有的優(yōu)勢菌綱是散囊菌綱(Euroti-omycetes)。

    圖4 野生羊肚菌根際土壤真菌在綱水平上的分類比較Fig.4 Taxonomic comparison of fungal in rhizosphere soil of wild morels at class level

    不同采樣點羊肚菌根際土壤在科水平的物種組成結(jié)果見圖5,9個采樣點中共檢測到350個科的真菌。在科水平,不同采樣點的優(yōu)勢菌科差異較大,DS、TB、LC采樣點,存在共同的優(yōu)勢菌科亞隔孢殼科小雙腔菌科(Didymellaceae);所有采樣點中,均存在一定量的unclassified,說明羊肚菌根際土壤存在一些未知的真菌;DS采樣點的優(yōu)勢菌科為小雙腔菌科(Didymellaceae)(23.03%)、暗殼腔菌科(Phaeosphaeriaceae)(10.63%)、unclassified_Fungi(9.75%)、unclassified(5.96%);DC采樣點的優(yōu)勢菌科為未分類的木耳目(unclassified_Auriculariales)(40.54%)、unclassified (7.65%)、格孢腔菌科(Sporormiaceae)(7.26%)、綠核蟲草(Cordycipitaceae)(5.21%);TB采樣點優(yōu)勢菌科為unclassified(11.96%)、小雙腔菌科(Didymellaceae)(10.15%)、球腔菌科(Mycosphaerellaceae)(9.49%),枝孢霉科(Cladosporiaceae)(7.66%);LC采樣點的優(yōu)勢菌科為 unclassified_Fungi(17.11%)、unclassified_Leotiomycetes(13.58%)、小雙腔菌科(Didymellaceae)(8.56%)、unclassified (7.27%);LK采樣點的優(yōu)勢菌科為盤菌科(Pezizaceae) (94.58%)、unclassified_Rozellomycota(2.31%);ZQ采樣點的優(yōu)勢菌科為曲霉科(Aspergillaceae)(19.18%)、unclassified(14.07%)、unclassified_GS11(13.13%);蠟傘科(Hygrophoraceae)(12.36%);DB采樣點的優(yōu)勢菌科為蠟殼耳科(Sebacinaceae)(16.36%)、unclassified(14.96%)、晶杯菌科(Hyaloscyphaceae)(8.34%)、unclassified_Fungi(7.10%);LQ采樣點的優(yōu)勢菌科為羊肚菌科(Morchellaceae)(89.96%)、莢胞腔菌科(Sporormiaceae)(0.98%);LX采樣點的優(yōu)勢菌科為柔膜菌科(Helotiaceae)(18.97%)、unclassified(7.86%)、unclassified_Ascomycota(7.81%)、unclassified_Leotiomycetes(6.82%)。

    不同采樣點羊肚菌根際土壤在屬水平的物種組成結(jié)果見圖6,9個采樣點中共檢測到775個屬的真菌。DS采樣點Paraboeremia(15.56%)、unclassified_phaeosphaeriaceae未分類的暗殼腔菌科(10.45%)、unclassified_Fungi(9.75%)、Nothophoma(6.25%);DC未分類的木耳屬類unclassified_Auriculariales(40.54%)、unclassified(7.65%)、光黑殼屬(Preussia)(7.08%)、白僵菌屬(Beauveria)(5.13%);TB采樣點未分類的木耳屬類unclassified_Auriculariales(20.31%)、Paraboeremia(7.78%)、unclassified_Fungi(6.86%)、unclassified(6.81%);LC采樣點未分類的真菌unclassified_Fungi(17.11%)、未分類的錘舌菌綱(unclassified_Leotiomycetes)(13.58%)、unclassified(6.27%)、被孢霉屬(Mortierella)(5.68%);LK采樣點盤菌屬(Peziza)(94.58%)、未分類的羅茲菌門unclassified_Rozellomycota(2.31%);ZQ采樣點青霉屬(Penicillium)(19.17%)、unclassified(14.07%)、unclassified_GS11(13.13%)、濕傘屬(Hygrocybe)(12.31%);DB采樣點蠟殼菌屬(Sebacina)(16.35%)、unclassified(14.96%)、未分類的真菌門類unclassified_Fungi(7.10%)、Gyoerffyella(6.24%);LQunclassified_Morchellaceae(89.34%)、LXunclassified_Helotiaceae(17.51%)、unclassified(7.86%)、unclassified_Ascomycota(7.81%)。

    圖5 野生羊肚菌根際土壤真菌在科水平上的分類比較Fig.5 Taxonomic comparison of fungal in rhizosphere soil of wild morels at family level

    圖6 野生羊肚菌根際土壤真菌在屬水平上的分類比較Fig.6 Taxonomic comparison of fungal in rhizosphere soil of wild morels at genus level

    2.5 野生羊肚菌土壤真菌群落結(jié)構(gòu)與土壤環(huán)境因子相互關(guān)系

    對9個采樣點野生羊肚菌根際土壤真菌群落進(jìn)行主成分分析,結(jié)果見圖7。真菌群落受主坐標(biāo)成分PCA1與PCA2的影響分別達(dá)到67%和16%。DS、LC、DB、DC、LX、TB采樣點真菌群落在PCA2軸上,ZQ采樣點距離PCA2軸很近,表明這6個采樣點的真菌群落結(jié)構(gòu)相似,受主成分PCA2的影響;LQ 采樣點在PCA1軸上,距離PCA2軸較遠(yuǎn),受主成分PCA1的影響較大;LK采樣點距離PCA2軸較遠(yuǎn),受主成分PCA1的影響較大;LQ和LK兩個采樣點的真菌群落與其他6個采樣點的真菌群落差異較大,LQ和LK采樣點真菌群落的差異最大。

    通過RDA分析,探明土壤環(huán)境因子與羊肚菌根際土壤真菌群落的相關(guān)關(guān)系,結(jié)果見圖8。第一、二排序軸累計解釋率分別為 55.2%和29.7%, 累計解釋率達(dá)到84.9%,說明第一、二排序軸能較好地反映真菌群落與土壤理化性質(zhì)、酶活性之間的相關(guān)關(guān)系。主要的真菌類群中,子囊菌門與速效鉀、谷氨酰胺酶、脲酶、硝態(tài)氮、堿性磷酸酶箭頭方向一致,說明子囊菌與速效鉀、谷氨酰胺酶、脲酶、硝態(tài)氮、堿性磷酸酶和脫氫酶正相關(guān);Rozellomycota和Olpidiomycota與硒、鐵、有效磷、氨態(tài)氮、脫氫酶呈正相關(guān),與鈣、多酚氧化酶、鎂、pH呈負(fù)相關(guān);壺菌門、unclassified_Fungi、Mortierellomycota與脫氫酶、堿性磷酸酶、硝態(tài)氮呈負(fù)相關(guān)。子囊菌門在LK與LQ采樣點優(yōu)勢明顯;LC、DS采樣點相聚較近,真菌菌群結(jié)構(gòu)相近;DC、DB采樣點相距較近,真菌菌群結(jié)構(gòu)相近;ZQ采樣點的真菌群落組成與其它采樣點差異較大。

    圖7 野生羊肚菌土壤真菌群結(jié)構(gòu)PCA分析Fig.7 PCA analysis of soil fungal community structure of wild morels

    圖8 野生羊肚菌根際土壤門水平真菌群落與環(huán)境因子的相關(guān)性RDA分析Fig.8 RDA analysis of the correlation between the fungal community and environmental factors in the rhizosphere soil phylum of wild morels紅色箭頭代表土壤酶活和理化性質(zhì)在水平面上的相對位置,箭頭與排序軸的夾角不同,則其余弦值不同,即二者的相關(guān)性強(qiáng)度也不同。夾角越小,相關(guān)性越大;箭頭的長度越長,代表該環(huán)境因子發(fā)揮的作用越大;藍(lán)色箭頭代表門水平上的物種分布,箭頭越長代表物種在樣品中的影響越大;黑色圓點代表不同采樣點。各真菌編號代表: F1:子囊菌門;F2:擔(dān)子菌門;F3:unclassified;F4:unclassified_Fungi;F5:Rozellomycota;F6:Mortierellomycota;F7:Cercozoa;F8:Glomeromycota;F9:壺菌門;F10:Olpidiomycota;AN:氨態(tài)氮;AP:有效磷;QAP:速效鉀;N:硝態(tài)氮;SH:含水率;ST:土壤溫度;TC:全碳;OM:有機(jī)質(zhì);PPO:多酚氧化酶;POD:過氧化物酶;UE:脲酶;ASF:芳基硫酸酯酶;β-GC:β-葡萄糖苷酶;DHA:脫氫酶;AKP:堿性磷酸酶;GLS:谷氨酰胺酶The red arrow represent the relative position of soil physical and chemical properties and enzyme activity on the horizontal plane, the angle between the arrow and the sort axis is different, the correlation strength is also different, the smaller the angle, the greater the correlation; The longer the length of the arrow, the greater the effect of the environmental factor; The blue arrow represents the species distribution at the class level, and the longer the arrow, the greater the impact of the species in the sample.F1:Ascomycota;F2:Basidiomycota;F3:unclassified;F4:unclassified_Fungi;F5:Rozellomycota;F6:Mortierellomycota;F7:Cercozoa;F8:Glomeromycota;F9:Chytridiomycota;F10:Olpidiomycota;AN:Ammonia nitrogen;AP:Available phosphorus;QAP:Rapidly available potassium;N:Nitrate;SH:Soil humidity;ST:Soil temperature;TC:Total carbon;OM:Organic material;PPO:Polyphenol oxidase;POD:Peroxidase;UE:Urease;ASF:Arylsulfatase;β-GC:β-glucosidase;DHA:Dehydrogenase;AKP:Alkaline phosphatase;GLS:glutaminase

    3 討 論

    土壤中微生物的種類、群落結(jié)構(gòu)對生態(tài)系統(tǒng)的結(jié)構(gòu)和功能具有重要的影響,土壤微生物與植物的互作已有很多的研究,而關(guān)于大型真菌與微生物互作的研究很少。本研究通過對甘肅不同地區(qū)羊肚菌根際真菌群落組成進(jìn)行系統(tǒng)分析,發(fā)現(xiàn)子囊菌門和擔(dān)子菌門在9個采樣地的真菌類群中均占優(yōu)勢地位,有研究顯示[22-23]子囊菌門是真菌豐度最高的真菌類群,在偏堿與偏酸性的土壤中均能大量生存,對土壤的適應(yīng)性較強(qiáng),這與本研究的結(jié)果一致;在綱、科、屬水平土壤真菌群落組成差異明顯,這與閻威文等[24]的研究結(jié)果一致。閻威文等的研究顯示野生羊肚菌根際土壤真菌優(yōu)勢菌屬為被孢霉屬(Mortierella),我們的研究結(jié)果中,只有在LC采樣點被孢霉屬是優(yōu)勢真菌群落之一。毋校輝[25]對四川8個地區(qū)羊肚菌種植根際土壤的高通量測序結(jié)果顯示有兩個采樣點的優(yōu)勢真菌類群為被孢霉屬,閻威文等[24]研究顯示栽培羊肚菌根際土壤真菌優(yōu)勢菌屬為紅菇屬(Russula)。由此可見,無論是野生羊肚菌還是栽培羊肚菌,其根際土壤真菌群落的差異均較大,另外,根際土壤的采集時間不同,真菌群落組成也有差異。

    研究表明,土壤微生物多樣性以及群落組成結(jié)構(gòu)受土壤理化性質(zhì)的影響顯著[26]。本研究中,羊肚菌在偏酸與偏堿的環(huán)境均可生長,pH值隨海拔的增加而減小。RDA分析結(jié)果顯示,除pH值外,有效磷、Ca、Se含量與真菌群落相關(guān)性較大,其他土壤理化性質(zhì)對根際土壤真菌群落的影響很??;尤其是土壤中的大量元素與真菌群落相關(guān)性小,其原因有待進(jìn)一步分析。有研究表明土壤微生物群落的空間變異因土壤化學(xué)性質(zhì)的不同而異[27],各地區(qū)土壤化學(xué)性質(zhì)存在差異,致使其土壤原始菌群不同,使各地羊肚菌生長受侵?jǐn)_程度不同,說明了羊肚菌栽培中土壤選擇的重要性。

    土壤酶是生態(tài)系統(tǒng)物質(zhì)循環(huán)以及能量轉(zhuǎn)化過程中重要的參與者[28],其主要來源為微生物。脫氫酶發(fā)揮的作用較大,脫氫酶能酶促脫氫反應(yīng),起著氫的中間傳遞體的作用;多酚氧化酶與優(yōu)勢真菌群落結(jié)構(gòu)呈負(fù)相關(guān),能將土壤中芳香族化合物氧化成醌,是表征對土壤中污染物凈化能力的指標(biāo)之一[29];磷酸酶是表征土壤碳、氮、磷等元素循環(huán)的因素之一,堿性磷酸酶與優(yōu)勢真菌群落子囊菌有顯著的相關(guān)性,其他酶活性對真菌群落的影響不大。通過比較,認(rèn)為土壤酶活性對真菌群落的影響較理化性質(zhì)影響大。

    羊肚菌人工種植的過程中,羊肚菌利用的碳源主要來源于外援營養(yǎng)袋,羊肚菌菌絲在地下如何代謝、與微生物的相互作用及如何積累營養(yǎng)物質(zhì)支持出菇還不清楚,進(jìn)一步開展野生羊肚菌菌絲生長期、原基形成期、子囊果成熟期土壤微生物組成與環(huán)境因子的相關(guān)性研究,能夠為人工種植羊肚菌提供理論參考。

    猜你喜歡
    子囊菌根羊肚
    紅心曲來源Thermoascus aurantiacus QH-1 的生長代謝特征
    羊肚菌“飛天記”
    軍事文摘(2023年18期)2023-10-31 08:10:18
    小小羊肚菌 拓寬致富路
    現(xiàn)代溫室羊肚菌栽培技術(shù)
    破壁方式對冠突散囊菌有性繁殖體提取物抑菌活性的影響
    外生菌根真菌菌劑的制備及保存研究
    園林科技(2020年2期)2020-01-18 03:28:26
    粗糙脈孢菌7種子囊型歸類教學(xué)探究
    遺傳(2019年11期)2019-11-28 12:03:42
    子囊霉素生物合成的研究進(jìn)展
    羊肚參的生藥學(xué)研究
    重金屬污染土壤的生物修復(fù)——菌根技術(shù)的應(yīng)用
    美女被艹到高潮喷水动态| 久久精品国产亚洲网站| 99热这里只有精品一区| 欧美人与善性xxx| ponron亚洲| 男人舔女人下体高潮全视频| 国产精品爽爽va在线观看网站| 亚洲精品国产av成人精品| 特大巨黑吊av在线直播| 99热精品在线国产| 欧美不卡视频在线免费观看| 国产真实伦视频高清在线观看| 亚洲内射少妇av| 性色avwww在线观看| 人妻久久中文字幕网| 在线国产一区二区在线| ponron亚洲| 国产成年人精品一区二区| 一个人观看的视频www高清免费观看| 韩国av在线不卡| 久久久精品大字幕| 国产白丝娇喘喷水9色精品| 久久人人精品亚洲av| 一边亲一边摸免费视频| 在线免费十八禁| 一本精品99久久精品77| 国产成人一区二区在线| 精品人妻视频免费看| 国产精品日韩av在线免费观看| 久久午夜福利片| 特大巨黑吊av在线直播| 久久久久久久久久久免费av| 久久婷婷人人爽人人干人人爱| 在线观看美女被高潮喷水网站| 欧美xxxx性猛交bbbb| 亚洲在线观看片| 欧美区成人在线视频| 国产一区二区在线观看日韩| 久久韩国三级中文字幕| 国产久久久一区二区三区| 内射极品少妇av片p| 此物有八面人人有两片| 淫秽高清视频在线观看| 国产亚洲欧美98| 少妇的逼水好多| 国内精品美女久久久久久| 亚洲综合色惰| 亚洲va在线va天堂va国产| 亚洲图色成人| 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| а√天堂www在线а√下载| a级毛色黄片| 麻豆一二三区av精品| 精品一区二区三区人妻视频| 亚洲国产欧美人成| 日韩一区二区三区影片| 国产欧美日韩精品一区二区| 久久久久久九九精品二区国产| 丰满人妻一区二区三区视频av| 午夜福利在线观看吧| 免费黄网站久久成人精品| 国产成人freesex在线| 国产三级在线视频| a级毛片a级免费在线| 简卡轻食公司| 精品久久久久久久人妻蜜臀av| 两个人的视频大全免费| 色综合站精品国产| 精品无人区乱码1区二区| 久久久欧美国产精品| 久久欧美精品欧美久久欧美| a级毛片a级免费在线| 国产中年淑女户外野战色| 国产精品野战在线观看| 亚洲真实伦在线观看| 久久久精品欧美日韩精品| 三级毛片av免费| 亚洲av中文字字幕乱码综合| 日韩av不卡免费在线播放| 长腿黑丝高跟| 免费在线观看成人毛片| 久久鲁丝午夜福利片| 日本色播在线视频| 好男人视频免费观看在线| 亚洲图色成人| 国产女主播在线喷水免费视频网站 | 国产淫片久久久久久久久| 欧美一级a爱片免费观看看| 欧美三级亚洲精品| 亚洲国产欧美人成| 午夜激情欧美在线| 天天一区二区日本电影三级| 日本色播在线视频| 欧美xxxx黑人xx丫x性爽| 国产精品永久免费网站| 国产免费一级a男人的天堂| 亚洲av成人av| 日本撒尿小便嘘嘘汇集6| 尾随美女入室| 久久婷婷人人爽人人干人人爱| 97热精品久久久久久| 级片在线观看| 午夜福利高清视频| 午夜亚洲福利在线播放| 久久精品国产鲁丝片午夜精品| 日本成人三级电影网站| 黑人高潮一二区| 亚洲欧美日韩高清专用| 91精品一卡2卡3卡4卡| 人人妻人人澡欧美一区二区| 两个人视频免费观看高清| 免费黄网站久久成人精品| 在线观看66精品国产| 亚洲av.av天堂| 九九热线精品视视频播放| 99国产精品一区二区蜜桃av| 国产日本99.免费观看| 高清毛片免费看| 青春草视频在线免费观看| 亚洲,欧美,日韩| 小说图片视频综合网站| 五月伊人婷婷丁香| 深夜精品福利| 亚洲aⅴ乱码一区二区在线播放| 久久久久九九精品影院| 一进一出抽搐动态| 亚洲欧美中文字幕日韩二区| 99在线视频只有这里精品首页| 两个人视频免费观看高清| 中国国产av一级| 精品一区二区免费观看| 国产不卡一卡二| 久久国产乱子免费精品| 国产成年人精品一区二区| 在线免费观看的www视频| 3wmmmm亚洲av在线观看| 亚洲成人中文字幕在线播放| 人体艺术视频欧美日本| 久久99精品国语久久久| 又粗又硬又长又爽又黄的视频 | 亚洲四区av| 亚洲美女搞黄在线观看| 中文字幕av在线有码专区| av天堂中文字幕网| 国产乱人视频| 久久久久久大精品| 久久久国产成人精品二区| 欧美激情久久久久久爽电影| 中文精品一卡2卡3卡4更新| 亚洲欧美中文字幕日韩二区| 欧美+日韩+精品| av国产免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 高清在线视频一区二区三区 | 国内少妇人妻偷人精品xxx网站| 色综合站精品国产| 一级av片app| 网址你懂的国产日韩在线| 美女黄网站色视频| 桃色一区二区三区在线观看| 欧美日韩一区二区视频在线观看视频在线 | 日本免费a在线| 成人亚洲精品av一区二区| 插逼视频在线观看| 中国国产av一级| 成人亚洲精品av一区二区| 插逼视频在线观看| 免费大片18禁| 日本一本二区三区精品| 亚洲乱码一区二区免费版| 一区二区三区四区激情视频 | 久久久久久国产a免费观看| 成人鲁丝片一二三区免费| 久久久久久久久久久免费av| 亚洲自偷自拍三级| 亚洲熟妇中文字幕五十中出| 成人鲁丝片一二三区免费| 国语自产精品视频在线第100页| 日韩欧美三级三区| 99热6这里只有精品| 日日摸夜夜添夜夜爱| 久久99热6这里只有精品| 亚洲色图av天堂| 99久久精品热视频| 精品人妻熟女av久视频| 老熟妇乱子伦视频在线观看| 亚洲欧洲日产国产| 亚洲精品乱码久久久v下载方式| 夜夜看夜夜爽夜夜摸| 欧美性猛交黑人性爽| 国产v大片淫在线免费观看| 熟女电影av网| a级毛片免费高清观看在线播放| 欧美成人精品欧美一级黄| 99视频精品全部免费 在线| 国产黄a三级三级三级人| 99热这里只有是精品在线观看| 国产精品国产高清国产av| 国产精品,欧美在线| 国产精品美女特级片免费视频播放器| 一个人看的www免费观看视频| 波多野结衣高清作品| 日韩一区二区三区影片| 国产成人精品久久久久久| 欧美一区二区亚洲| 国产成人精品一,二区 | 日本一本二区三区精品| 日韩中字成人| 精品无人区乱码1区二区| 免费无遮挡裸体视频| 亚洲av成人av| 99久国产av精品| 婷婷精品国产亚洲av| 精品人妻视频免费看| 久久久精品94久久精品| 日本欧美国产在线视频| 国产日韩欧美在线精品| 欧美三级亚洲精品| 日日撸夜夜添| 亚洲四区av| 日韩欧美国产在线观看| 九九在线视频观看精品| 99久久精品热视频| 亚洲精品国产av成人精品| 成年av动漫网址| 成人永久免费在线观看视频| 国产三级在线视频| 一区二区三区四区激情视频 | 欧美日韩在线观看h| 人体艺术视频欧美日本| 午夜a级毛片| 女人十人毛片免费观看3o分钟| 免费在线观看成人毛片| 久久精品国产亚洲av涩爱 | 久久人人爽人人片av| 天天一区二区日本电影三级| 嫩草影院精品99| 尤物成人国产欧美一区二区三区| 国产精品无大码| 日本欧美国产在线视频| 岛国毛片在线播放| 青青草视频在线视频观看| 级片在线观看| 亚洲自拍偷在线| 亚洲欧美日韩高清专用| 亚洲人成网站在线播放欧美日韩| 日日干狠狠操夜夜爽| 国产欧美日韩精品一区二区| 久久久久久久久久黄片| 九九爱精品视频在线观看| 在线观看av片永久免费下载| 欧美成人a在线观看| 久久久久久久久中文| 久久久成人免费电影| 成年免费大片在线观看| 精品一区二区三区人妻视频| 午夜老司机福利剧场| 亚洲性久久影院| 国产男人的电影天堂91| 一进一出抽搐动态| 97超碰精品成人国产| 亚洲自拍偷在线| 午夜福利在线观看吧| 久久久精品大字幕| 2022亚洲国产成人精品| av卡一久久| 老师上课跳d突然被开到最大视频| 成人毛片60女人毛片免费| 久久精品国产清高在天天线| 国产精品嫩草影院av在线观看| 国产在视频线在精品| 精品熟女少妇av免费看| 国产精品,欧美在线| 美女cb高潮喷水在线观看| 国内久久婷婷六月综合欲色啪| 久久久久久久久久久免费av| 色综合亚洲欧美另类图片| 亚洲欧美清纯卡通| 国产在线精品亚洲第一网站| 寂寞人妻少妇视频99o| 国产一区二区激情短视频| 美女 人体艺术 gogo| 日日摸夜夜添夜夜添av毛片| 黄片wwwwww| 久久久久九九精品影院| 在线a可以看的网站| 久久综合国产亚洲精品| 91精品一卡2卡3卡4卡| 又粗又硬又长又爽又黄的视频 | 亚洲真实伦在线观看| 精品欧美国产一区二区三| 老熟妇乱子伦视频在线观看| av在线播放精品| 在线国产一区二区在线| 国产精品三级大全| 99精品在免费线老司机午夜| 亚洲国产欧洲综合997久久,| 国产精品永久免费网站| 婷婷六月久久综合丁香| 欧美极品一区二区三区四区| 欧美日韩精品成人综合77777| 一进一出抽搐动态| 国内精品美女久久久久久| 国产精品久久视频播放| 国产黄片美女视频| 夜夜爽天天搞| 黄色一级大片看看| 一级二级三级毛片免费看| 青春草视频在线免费观看| 久久精品国产清高在天天线| 精品人妻一区二区三区麻豆| 亚洲综合色惰| 高清在线视频一区二区三区 | 精品午夜福利在线看| 在线播放国产精品三级| 99久久精品热视频| 亚洲国产精品国产精品| 草草在线视频免费看| 久久久a久久爽久久v久久| 久久99精品国语久久久| 桃色一区二区三区在线观看| av在线观看视频网站免费| 日韩 亚洲 欧美在线| 日本黄大片高清| 亚洲精品久久国产高清桃花| 九九在线视频观看精品| 毛片一级片免费看久久久久| 国产精品乱码一区二三区的特点| 午夜亚洲福利在线播放| 成人永久免费在线观看视频| 日本免费a在线| 国产亚洲精品av在线| 欧美一区二区国产精品久久精品| 国产精品野战在线观看| 国产精品人妻久久久影院| 中文在线观看免费www的网站| 久久精品国产自在天天线| 亚洲一区高清亚洲精品| 91aial.com中文字幕在线观看| 最近的中文字幕免费完整| 色哟哟哟哟哟哟| 久久国内精品自在自线图片| 亚洲一区二区三区色噜噜| 午夜免费激情av| 亚洲色图av天堂| 亚洲欧美日韩卡通动漫| 三级经典国产精品| 18禁在线无遮挡免费观看视频| www.av在线官网国产| 国产黄色小视频在线观看| 99久久久亚洲精品蜜臀av| 一本久久中文字幕| 午夜福利在线观看免费完整高清在 | 亚洲国产日韩欧美精品在线观看| 嫩草影院精品99| 午夜免费男女啪啪视频观看| a级毛片a级免费在线| 赤兔流量卡办理| 少妇的逼好多水| 联通29元200g的流量卡| 97在线视频观看| 成人欧美大片| 91久久精品国产一区二区三区| 看片在线看免费视频| 国产国拍精品亚洲av在线观看| 十八禁国产超污无遮挡网站| 日韩三级伦理在线观看| 成人一区二区视频在线观看| 欧美一区二区国产精品久久精品| 国产伦精品一区二区三区四那| 亚洲熟妇中文字幕五十中出| 可以在线观看毛片的网站| 中文精品一卡2卡3卡4更新| 日本在线视频免费播放| 日韩欧美精品v在线| 毛片女人毛片| 欧美激情久久久久久爽电影| or卡值多少钱| 国产成人a区在线观看| 美女大奶头视频| 亚洲av免费高清在线观看| 一个人看视频在线观看www免费| 亚州av有码| 成人二区视频| 国产乱人视频| 久久国内精品自在自线图片| 看十八女毛片水多多多| 青青草视频在线视频观看| 99热这里只有是精品50| avwww免费| 久久精品国产亚洲网站| 两个人的视频大全免费| 18禁裸乳无遮挡免费网站照片| 欧美+日韩+精品| 久久精品影院6| 欧洲精品卡2卡3卡4卡5卡区| av天堂在线播放| 毛片女人毛片| 在线观看66精品国产| 青青草视频在线视频观看| 大又大粗又爽又黄少妇毛片口| 国内揄拍国产精品人妻在线| 国产一级毛片在线| 日韩精品有码人妻一区| 国产在视频线在精品| 亚洲av免费在线观看| 免费观看精品视频网站| 国产精品爽爽va在线观看网站| 欧美日韩国产亚洲二区| 免费搜索国产男女视频| 插逼视频在线观看| 乱系列少妇在线播放| 国产精品久久久久久亚洲av鲁大| 男人狂女人下面高潮的视频| 久久99热这里只有精品18| 尤物成人国产欧美一区二区三区| 亚洲熟妇中文字幕五十中出| 国产成人午夜福利电影在线观看| 少妇裸体淫交视频免费看高清| 亚洲国产精品合色在线| 免费人成视频x8x8入口观看| 精品熟女少妇av免费看| 热99re8久久精品国产| 亚洲欧美精品专区久久| 日本三级黄在线观看| 97超视频在线观看视频| 看片在线看免费视频| 偷拍熟女少妇极品色| 黄色一级大片看看| 99久久精品国产国产毛片| 国产女主播在线喷水免费视频网站 | 九九热线精品视视频播放| 国产午夜精品论理片| 成人午夜高清在线视频| 美女脱内裤让男人舔精品视频 | 欧美区成人在线视频| 欧美色视频一区免费| 男女啪啪激烈高潮av片| 如何舔出高潮| 在线免费观看的www视频| 99精品在免费线老司机午夜| 免费观看人在逋| 国语自产精品视频在线第100页| 男人的好看免费观看在线视频| 欧美精品一区二区大全| 桃色一区二区三区在线观看| 婷婷六月久久综合丁香| 99精品在免费线老司机午夜| 我要搜黄色片| 在线观看免费视频日本深夜| 亚洲最大成人手机在线| 亚洲av中文av极速乱| 又爽又黄无遮挡网站| 草草在线视频免费看| 国产成人精品一,二区 | 国产老妇女一区| 在现免费观看毛片| 99久久无色码亚洲精品果冻| 国产白丝娇喘喷水9色精品| 毛片女人毛片| 免费看a级黄色片| 国产精品久久电影中文字幕| 亚洲精品色激情综合| 中文字幕av在线有码专区| 亚洲七黄色美女视频| 国产亚洲精品久久久久久毛片| 国产人妻一区二区三区在| 黄片无遮挡物在线观看| 色尼玛亚洲综合影院| 一边亲一边摸免费视频| 伊人久久精品亚洲午夜| 日韩高清综合在线| 男女啪啪激烈高潮av片| 一卡2卡三卡四卡精品乱码亚洲| 99热网站在线观看| 国产午夜精品一二区理论片| 色播亚洲综合网| 欧美三级亚洲精品| 国产视频内射| 欧美日韩精品成人综合77777| 老女人水多毛片| 国产一区二区亚洲精品在线观看| 久久鲁丝午夜福利片| 禁无遮挡网站| 国产极品天堂在线| 日本在线视频免费播放| 国产精品日韩av在线免费观看| 不卡视频在线观看欧美| 国产亚洲精品久久久com| 日产精品乱码卡一卡2卡三| 亚洲美女搞黄在线观看| 国产精品三级大全| 少妇高潮的动态图| 午夜精品在线福利| 国产真实伦视频高清在线观看| 天堂影院成人在线观看| 嫩草影院精品99| 国产精品1区2区在线观看.| 成人av在线播放网站| 亚洲国产精品sss在线观看| 老女人水多毛片| 高清日韩中文字幕在线| 午夜免费男女啪啪视频观看| 亚洲高清免费不卡视频| 又粗又爽又猛毛片免费看| 成熟少妇高潮喷水视频| 91午夜精品亚洲一区二区三区| 欧美极品一区二区三区四区| 一边亲一边摸免费视频| 九草在线视频观看| 1024手机看黄色片| 婷婷色综合大香蕉| 联通29元200g的流量卡| 欧美性猛交╳xxx乱大交人| 一级黄色大片毛片| 狂野欧美白嫩少妇大欣赏| 婷婷色综合大香蕉| 99热精品在线国产| 日日摸夜夜添夜夜爱| 悠悠久久av| 精品日产1卡2卡| av又黄又爽大尺度在线免费看 | 日本黄色视频三级网站网址| 美女高潮的动态| 国产精品国产三级国产av玫瑰| 美女被艹到高潮喷水动态| 欧美zozozo另类| 精品一区二区三区人妻视频| 99久国产av精品| 国产成人精品婷婷| 日本熟妇午夜| av女优亚洲男人天堂| 亚洲精品日韩av片在线观看| 中文字幕久久专区| 欧美在线一区亚洲| 欧美一区二区精品小视频在线| 校园春色视频在线观看| 九色成人免费人妻av| 少妇熟女aⅴ在线视频| 99在线视频只有这里精品首页| 日韩中字成人| 国产激情偷乱视频一区二区| 91久久精品国产一区二区成人| 国产爱豆传媒在线观看| 国产91av在线免费观看| 亚洲无线在线观看| 亚洲人成网站在线播放欧美日韩| 欧美丝袜亚洲另类| 免费电影在线观看免费观看| 日韩成人伦理影院| 久久九九热精品免费| av免费观看日本| 国产精品电影一区二区三区| av卡一久久| 久久久午夜欧美精品| 乱码一卡2卡4卡精品| 青春草亚洲视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美日韩高清在线视频| 一夜夜www| 亚洲va在线va天堂va国产| 国产成人一区二区在线| 一本一本综合久久| 国产不卡一卡二| 精品久久久久久久人妻蜜臀av| 欧美成人一区二区免费高清观看| av天堂在线播放| 久久久久久久亚洲中文字幕| 女人十人毛片免费观看3o分钟| 国产乱人视频| 国产大屁股一区二区在线视频| 日韩成人av中文字幕在线观看| 边亲边吃奶的免费视频| 99久久人妻综合| 欧美性猛交黑人性爽| 国产亚洲精品久久久久久毛片| 少妇熟女欧美另类| 丰满的人妻完整版| 亚洲无线观看免费| 99在线视频只有这里精品首页| 亚洲美女视频黄频| 麻豆一二三区av精品| 插阴视频在线观看视频| 国产精品福利在线免费观看| 一本精品99久久精品77| 校园人妻丝袜中文字幕| 99久久久亚洲精品蜜臀av| 直男gayav资源| 中国国产av一级| 一级黄色大片毛片| 国产69精品久久久久777片| 哪里可以看免费的av片| 国产伦理片在线播放av一区 | 国产成人精品一,二区 | 精品熟女少妇av免费看| 欧美极品一区二区三区四区| 国产精品爽爽va在线观看网站| 精品熟女少妇av免费看| 免费黄网站久久成人精品| 欧美性猛交╳xxx乱大交人| 我的老师免费观看完整版| 直男gayav资源| 校园人妻丝袜中文字幕| 亚洲激情五月婷婷啪啪| 亚洲人成网站高清观看| 女人十人毛片免费观看3o分钟| 神马国产精品三级电影在线观看| 又黄又爽又刺激的免费视频.| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品久久久久久婷婷小说 | 精品少妇黑人巨大在线播放 | 欧美日韩国产亚洲二区| 国产精品久久久久久精品电影小说 | 国产一区二区在线av高清观看|