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

    基于高通量測(cè)序的明日葉不同組織內(nèi)生菌群落結(jié)構(gòu)分析

    2023-11-02 16:01:07劉雯雯張旭黃奇胡連清陳露周萬(wàn)海魏琴馮瑞章趙鑫
    關(guān)鍵詞:內(nèi)生真菌多樣性

    劉雯雯 張旭 黃奇 胡連清 陳露 周萬(wàn)海 魏琴 馮瑞章 趙鑫

    DOI:10.3969/j.issn.2095-1191.2023.06.004

    摘要:【目的】分析明日葉不同組織內(nèi)生菌的群落結(jié)構(gòu)和多樣性,為探究明日葉和內(nèi)生菌互作機(jī)制及功能菌株篩選提供參考?!痉椒ā恳运拇ㄒ速e栽培的藥食兼用型明日葉為試材,采用Illumina MiSeq高通量測(cè)序技術(shù)確定根、莖、葉和果實(shí)組織內(nèi)生真菌和內(nèi)生細(xì)菌群落多樣性?!窘Y(jié)果】根據(jù)97%的序列相似性水平,將內(nèi)生真菌和細(xì)菌的reads分別歸類為1023和1871個(gè)操作分類單元(OTU)。韋恩圖和熱圖顯示明日葉各組織內(nèi)生真菌和內(nèi)生細(xì)菌群落結(jié)構(gòu)差異性較高,種群豐度存在一定的差異。α多樣性分析發(fā)現(xiàn)明日葉的葉內(nèi)生真菌種類數(shù)最多,根內(nèi)生真菌群落多樣性最低,果實(shí)內(nèi)生真菌群落豐富度最低;而根含有內(nèi)生細(xì)菌種類數(shù)最多,莖內(nèi)生細(xì)菌群落多樣性最低,果實(shí)內(nèi)生細(xì)菌群落豐富度最低。β多樣性分析表明,組織特征對(duì)明日葉內(nèi)生真菌或內(nèi)生細(xì)菌群落結(jié)構(gòu)有顯著影響。明日葉4個(gè)不同組織中的優(yōu)勢(shì)真菌門和細(xì)菌門分別為子囊菌門(Ascomycota)和變形菌門(Proteobacteria)。根、莖、葉和果實(shí)中內(nèi)生真菌優(yōu)勢(shì)屬分別為Serendipita、Paramycosphaerella、多臂菌屬(Trichomerium)和Apiotrichum,內(nèi)生細(xì)菌優(yōu)勢(shì)屬分別為羅爾斯通菌屬(Ralstonia)、羅爾斯通菌屬、甲基桿菌屬(Methylobacterium)和弧菌屬(Vibrio)。線性判別分析(LEfSe)表明,除根外,明日葉的莖、葉和果實(shí)中內(nèi)生真菌重要差異指示種分別是Paramycosphaerella、蒼白尾孢屬(Pallidocercospora)和橫斷孢屬(Strelitziana);莖、葉和果實(shí)中內(nèi)生細(xì)菌重要差異指示種分別是羅爾斯通菌屬、甲基桿菌屬和亮桿菌屬(Leucobacter)。經(jīng)FUNGuild和COG功能預(yù)測(cè)發(fā)現(xiàn),不同組織內(nèi)生菌群功能豐富,除具有固氮、促生和抑菌等功能的菌群外,還存在較大比例未知功能菌群?!窘Y(jié)論】明日葉中蘊(yùn)含豐富多樣的內(nèi)生菌群,在不同組織中的群落結(jié)構(gòu)、多樣性、優(yōu)勢(shì)菌屬及差異指示種存在保守性與差異性,組織特征差異是影響內(nèi)生菌群的重要因素;不同組織內(nèi)生菌群生態(tài)功能亦差異較大,且大量為未鑒定出種屬和未定義的功能菌群。

    關(guān)鍵詞:明日葉;內(nèi)生真菌;內(nèi)生細(xì)菌;多樣性;功能預(yù)測(cè)

    中圖分類號(hào):S567.239? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼:A 文章編號(hào):2095-1191(2023)06-1622-11

    Community structure analysis of the endophyte in different

    tissues of Angelica keiskei (Miq.) Koidz. based on

    high-throughput sequencing

    LIU Wen-wen1,2, ZHANG Xu3, HUANG Qi2, HU Lian-qing2, CHEN Lu2,

    ZHOU Wan-hai1,2*, WEI Qin1, FENG Rui-zhang1*, ZHAO Xin1

    (1Sichuan Cinnamomum longepaniculatum Engineering Technology Research Center, Yibin University,Yibin,

    Sichuan? 644007, China; 2Faculty of Agriculture,F(xiàn)orestry and Food Engineering,Yibin University,Yibin,

    Sichuan? 644007,China; 3School of Life Science, Hainan University,Haikou, Hainan? 570228,China)

    Abstract:【Objective】The community structure and diversity of endophyte in different tissues of Angelica keiskei (Miq.)Koidz. was analyzed, to provide a theoretical basis to explore the interaction mechanism and functional strain screening of A. keiskei and endophyte. 【Method】The medicinal and edible plant A. keiskei, which was cultivated in Yibin, Sichuan, was taken as the object. The diversity of endophytic fungi and bacteria in the root, stem, leave and fruit tissues was examined by high-throughput sequencing using Illumina MiSeq. 【Result】According to the 97% sequence similarity, endophytic fungal and endophytic bacterial reads could be indentified into 1023 and 1871 operational taxonomic units(OTU), respectively. A Venn diagram and a heatmap showed that there were fairly high community structure difference between endophytic fungi or endophytic bacteria in the four different parts of tissues of A. keiskei, and the abundances of populations differed. According to α diversity analysis, the variety of endophytic fungi in the leaves of A. keiskei was the most, the diversity of endophytic fungal community in the roots was the lowest, and the abundance of endophytic fungal community in the fruits was the lowest. The variety of endophytic bacteria in the roots of A. keiskei was the most, the diversity of endophytic bacterial community in the stems was the lowest, and the abundance of endophytic bacterial community in the fruits was the lowest. β diversity analysis showed that the characteristics of tissues had a significant effect on the community structure of the endophytic fungi or endophytic bacteria in A. keiskei. Ascomycota and Proteobacteria were the dominant fungal phylum and bacterial phylum in the four different parts of tissues of A. keiskei respectively. The dominant genera of endophytic fungi in the roots, stems, leaves and fruits were Serendipita, Paramycosphaerella, Trichomerium and Apiotrichum respectively, and the dominant genera of endophytic bacteria were Ralstonia, Methylobacterium and Vibrio respectively. Linear discriminant analysis showed that, the most important biomarkers of the endophytic fungi in the stems, leaves and fruits were Paramycosphaerella, Pallidocercospora and Strelitziana respectively. The most important biomarkers of endophytic bacteria in stems, leaves and fruits were Ralstonia, Methylobacterium and Leucobacter respectively. FUNGuild and COG functional predictions showed that the function of endophyte in different tissues was abundant, and there was a large proportion of unknown functional endophyte in addition to those that could fix nitrogen, promote growth and inhibit bacteria. 【Conclusion】There are abundant and diverse endophytes in A. keiskei. The community structure, diversity, dominant genera and biomarkers in different tissues are conservative and different, and the difference of tissue characteristics is an important factor affecting the endophyte. The ecological functions of endophyte in different tissues are also quite different, and a large number of species and functional endophyte are not identified.

    Key words: Angelica keiskei(Miq.) Koidz.; endophytic fungi; endophytic bacteria; diversity; function prediction

    Foundation items:Sichuan National Natural Science Foundation (2022NSFSC0176);Application Project of Sichuan Science and Technology Department(2018JY0495);High-level Innovation and Entrepreneurship Team Project of Yibin Science and Technology Bureau(2021YGC03)

    0 引言

    【研究意義】明日葉[Angelica keiskei (Miq.) Koidz.]屬傘形科當(dāng)歸屬草本植物,是一種經(jīng)濟(jì)價(jià)值極高的藥食兩用植物(Gong et al.,2020)。藥用植物的次生產(chǎn)物積累、藥用品質(zhì)和抗病等與組織內(nèi)微生物有復(fù)雜的作用關(guān)系,內(nèi)生菌能對(duì)宿主植物物質(zhì)循環(huán)和能量轉(zhuǎn)換過(guò)程中的信號(hào)識(shí)別、信號(hào)轉(zhuǎn)導(dǎo)和基因表達(dá)等產(chǎn)生顯著影響,從而促進(jìn)植物對(duì)營(yíng)養(yǎng)成分的吸收和代謝(蔡媛等,2021),部分內(nèi)生菌還可代謝產(chǎn)生與宿主相同或相似作用的天然活性成分,此特征已成為挖掘新型藥物的重要來(lái)源(郭文秀等,2022)。因此,開(kāi)展特色藥食兩用植物明日葉不同組織的內(nèi)生菌群落結(jié)構(gòu)及多樣性研究,有助于揭示明日葉不同組織內(nèi)生菌群落結(jié)構(gòu)特征,且對(duì)挖掘可提高其有效成分含量的內(nèi)生菌資源及開(kāi)發(fā)促進(jìn)明日葉生長(zhǎng)、生物轉(zhuǎn)化和品質(zhì)提升的方法具有重要意義?!厩叭搜芯窟M(jìn)展】目前已對(duì)蛇足石衫、黃精和白蘞等多種藥用植物內(nèi)生菌群落結(jié)構(gòu)、多樣性和組織分布特征等開(kāi)展了相關(guān)研究(Fan et al.,2020;李鑒濱等,2022),亦對(duì)內(nèi)生菌響應(yīng)宿主植物內(nèi)環(huán)境因子及調(diào)節(jié)宿主植物生長(zhǎng)、活性次生代謝物積累和分布進(jìn)行了分析(黃雪珍和趙龍飛,2023),發(fā)現(xiàn)內(nèi)生菌分布于植物器官或組織中,每個(gè)植物器官均代表內(nèi)生細(xì)菌群落特有的生態(tài)位。Dong等(2018)對(duì)中草藥三七不同部位的內(nèi)生細(xì)菌進(jìn)行多樣性分析,得出根部的慢生型大豆根瘤菌、鞘氨醇單胞菌和苯基桿菌等內(nèi)生細(xì)菌數(shù)量均高于莖部和葉部。許國(guó)琪等(2021)發(fā)現(xiàn)青藏高原冰川棘豆不同組織間形成了保守性內(nèi)生細(xì)菌和特異性內(nèi)生真菌共存的群落結(jié)構(gòu),為其不同組織適應(yīng)極端環(huán)境提供幫助。陳美琪等(2023)發(fā)現(xiàn)對(duì)虎耳草不同組織內(nèi)生真菌的多樣性指數(shù)均高于多數(shù)已報(bào)道的其他藥用植物,且根部?jī)?nèi)生真菌的多樣性明顯高于莖部和葉部。張麗芳等(2023)發(fā)現(xiàn)火棘根中內(nèi)生細(xì)菌物種最豐富,參與各種代謝調(diào)控的細(xì)菌豐度最高。關(guān)于明日葉的研究主要集中在栽培技術(shù)、活性物質(zhì)提取與鑒定及中藥藥性探析等方面。關(guān)于內(nèi)生菌,本課題組前期采用純培養(yǎng)法篩選明日葉不同器官可培養(yǎng)的內(nèi)生真菌和產(chǎn)黃酮菌株,淺析了明日葉不同組織內(nèi)生真菌群落結(jié)構(gòu)和多樣性,但分離的菌株數(shù)量和類型較少,多樣性較低(劉雯雯等,2022)。大量研究表明傳統(tǒng)的平板分離法并不能真正反映內(nèi)生菌群落特征,高通量測(cè)序技術(shù)能提供更加可靠和深入的分析(陳美琪等,2023)?!颈狙芯壳腥朦c(diǎn)】有關(guān)明日葉內(nèi)生細(xì)菌的群落結(jié)構(gòu)及組織差異尚缺乏系統(tǒng)的認(rèn)識(shí)?!緮M解決的關(guān)鍵問(wèn)題】收集四川宜賓栽培的藥食兼用型明日葉的根、莖、葉和果實(shí)組織,基于高通量測(cè)序技術(shù)分析明日葉不同組織中內(nèi)生真菌和內(nèi)生細(xì)菌群落組成和多樣性,并比較內(nèi)生真菌和內(nèi)生細(xì)菌群落結(jié)構(gòu)差異,探究明日葉中特有的微生物資源及其功能,為探尋功能微生物資源和明日葉的功能開(kāi)發(fā)應(yīng)用,挖掘能合成或促進(jìn)植物活性產(chǎn)物的菌株資源提供理論依據(jù)和技術(shù)支撐。

    1 材料與方法

    1. 1 試驗(yàn)區(qū)概況

    明日葉采樣點(diǎn)位于四川省宜賓市宜賓學(xué)院試驗(yàn)農(nóng)場(chǎng)(東經(jīng)104°59′76″,北緯28°47′77″)。屬中亞熱帶季風(fēng)濕潤(rùn)氣候,氣候溫和,熱量豐足,雨量充沛,光照適宜,冬暖春早,無(wú)霜期為347 d,年平均溫度17.5 ℃,年均降水量1070.4 mm。

    1. 2 樣品采集和處理

    在采集地選擇不同位點(diǎn)的5株3年生且長(zhǎng)勢(shì)良好的明日葉植株,將整棵植株挖出,盡量保持根系完整,立即放入內(nèi)置冰袋的冰盒內(nèi)運(yùn)送至實(shí)驗(yàn)室,用大量無(wú)菌水對(duì)明日葉植株表面和根系附著泥土進(jìn)行沖洗后,用無(wú)菌剪刀將不同位品按根、莖、葉和果實(shí)組織分開(kāi)。將不同位點(diǎn)每株明日葉不同組織的混合樣品作為一個(gè)生物學(xué)重復(fù),共設(shè)3個(gè)重復(fù)。稱取50 g的各組織樣品于超凈工作臺(tái)中無(wú)菌水沖洗多次,再依次用75%乙醇(40 s)、3% NaClO(3 min)、75%乙醇(40 s)無(wú)菌化處理,最后無(wú)菌水洗滌3次,無(wú)菌濾紙吸干表面水分,置于-80 ℃超低溫冰箱中保存,用于明日葉內(nèi)生菌群落結(jié)構(gòu)分析。

    1. 3 明日葉微生物DNA提取與測(cè)序

    明日葉不同組織樣本的總基因組DNA使用FastDNA?SPIN Kit for Soil試劑盒進(jìn)行提取,借助1%瓊脂糖凝膠電泳和NanoDrop 2000超微量分光光度計(jì)分別測(cè)定DNA提取質(zhì)量和DNA濃度、純度。PCR使用引物ITS1F(5'-CTTGGTCATTTAGAGGA AGTAA-3')和ITS2R(5'-GCTGCGTTCTTCATCGAT GC-3')對(duì)內(nèi)生真菌ITS1區(qū)進(jìn)行擴(kuò)增;使用引物799F(5'-AACMGGATTAGATACCCKG-3')和1193R(5'-AC GTCATCCCCACCTTCC-3')對(duì)內(nèi)生細(xì)菌16S rDNA的V5~V7區(qū)進(jìn)行擴(kuò)增。利用凝膠回收試劑盒(Gel Extraction Kit)對(duì)擴(kuò)增目標(biāo)片段切膠回收,2%瓊脂糖凝膠電泳和Quantus?熒光儀對(duì)產(chǎn)物檢測(cè)并定量,使用NEXTflexTM DNA-Seq Kit建庫(kù)。最后利用Illumina公司Miseq PE300平臺(tái)進(jìn)行測(cè)序。上述工作均委托上海美吉生物醫(yī)藥科技有限公司完成。

    1. 4 測(cè)序數(shù)據(jù)處理

    將測(cè)序獲得的數(shù)據(jù)過(guò)濾去除引物和接頭序列、非特異性擴(kuò)增序列和嵌合體,得到有效序列用于進(jìn)一步分析。按照97%相似水平進(jìn)行操作分類單元(Operation taxonomic unit,OTU)聚類,利用R軟件繪制韋恩圖后展開(kāi)OTU聚類分析。α多樣性指數(shù)包括Shannon指數(shù)(反映內(nèi)生菌群落多樣性)和Chao1指數(shù)(反映內(nèi)生菌群落豐富度)。通過(guò)對(duì)OTU數(shù)據(jù)分析,利用柱狀圖統(tǒng)計(jì)樣本中內(nèi)生真菌和內(nèi)生細(xì)菌在門和屬水平上群落結(jié)構(gòu)的分類學(xué)組成,確定優(yōu)勢(shì)菌門和優(yōu)勢(shì)菌屬?;谥饕狾TU豐度信息繪制熱圖,反映不同樣本間物種分布的相似性和差異性。在ANOSIM組間差異檢驗(yàn)下進(jìn)行NMDS分析,用于判斷分組是否有意義。根據(jù)分類學(xué)組成對(duì)樣本進(jìn)行線性判別分析(LEfSe),找出對(duì)樣本劃分產(chǎn)生顯著性差異影響的群落或物種。使用FUNGuild和EggNOG數(shù)據(jù)庫(kù)分別預(yù)測(cè)內(nèi)生真菌和內(nèi)生細(xì)菌的功能,統(tǒng)計(jì)內(nèi)生真菌Trophic Mode(營(yíng)養(yǎng)方式)和Guild(功能分組)豐度結(jié)果,比對(duì)內(nèi)生細(xì)菌OTU的COG家族信息。

    1. 5 統(tǒng)計(jì)分析

    試驗(yàn)數(shù)據(jù)采用SPSS 22.0和Excel 2007進(jìn)行統(tǒng)計(jì)分析。利用上海美吉生物醫(yī)藥科技有限公司的I-sanger云數(shù)據(jù)分析平臺(tái)進(jìn)行在線數(shù)據(jù)分析。

    2 結(jié)果與分析

    2. 1 明日葉植物內(nèi)生菌OTU數(shù)量

    明日葉內(nèi)生真菌OTUs總數(shù)為1023個(gè),其中根、莖、葉和果實(shí)分別含有210、210、426和177個(gè),4個(gè)組織共有內(nèi)生真菌OTUs 29個(gè);葉部特有OTUs數(shù)量最多,為249個(gè),莖部次之,果實(shí)特有OTUs數(shù)量最少,僅54個(gè)(圖1-A)。明日葉各組織內(nèi)生細(xì)菌共獲得1871個(gè)OTUs,其中根、莖、葉和果實(shí)分別含有1114、243、298和216個(gè),4個(gè)組織共有內(nèi)生細(xì)菌OTUs 55個(gè);根特有OTUs數(shù)量最多,為884個(gè),其次是葉116個(gè),莖最少,僅41個(gè)(圖1-B)。

    2. 2 明日葉植物內(nèi)生菌α和β多樣性

    2. 2. 1 內(nèi)生真菌α多樣性 測(cè)序共獲得247037條內(nèi)生真菌有效序列,平均測(cè)序覆蓋率達(dá)99.962%。由Shannon指數(shù)箱式圖(圖2-A)可知,明日葉的根與莖、根與葉差異極顯著(P<0.01,下同),其他組織間無(wú)顯著差異(P>0.05,下同),Shannon指數(shù)大小依次是葉>莖>果實(shí)>根;Chao1指數(shù)箱式圖(圖2-B)顯示葉與其他3個(gè)組織間均差異極顯著,Chao1指數(shù)大小依次為葉>莖>根>果實(shí)。說(shuō)明葉部?jī)?nèi)生真菌群落多樣性和豐富度均最好。

    2. 2. 2 內(nèi)生細(xì)菌α多樣性 測(cè)序共獲得87178條內(nèi)生細(xì)菌有效序列,平均測(cè)序覆蓋率為99.808%。結(jié)合Shannon指數(shù)和Chao1指數(shù)箱式圖(圖3)可知,明日葉根的Shannon指數(shù)與其他組織間均差異極顯著,而其他組織間無(wú)顯著性差異;果實(shí)的Chaol指數(shù)與根和葉間差異顯著(P<0.05,下同)。因此,在Shannon指數(shù)上各組織排序依次是根>葉>果實(shí)>莖,Chao1指數(shù)大小依次為根>葉>莖>果實(shí)。說(shuō)明根部?jī)?nèi)生細(xì)菌的群落多樣性和豐富度均最好,葉部其次。

    2. 2. 3 內(nèi)生菌β多樣性 基于ANOSIM組間差異檢驗(yàn)的NMDS分析表明,組織間內(nèi)生真菌OTU(R=0.4815,P=0.007)和內(nèi)生細(xì)菌OTU(R=0.7407,P=0.002)均存在極顯著差異(圖4)。Stress主要用于檢驗(yàn)NMDS分析結(jié)果的優(yōu)劣,通常認(rèn)為Stress<0.1時(shí),不同樣本具有一個(gè)好的排序,內(nèi)生真菌NMDS分析Stress=0.090,內(nèi)生細(xì)菌Stress=0.073,因此在內(nèi)生真菌和內(nèi)生細(xì)菌菌群結(jié)構(gòu)的差異分析中,分組是一個(gè)好的排序。

    2. 3 明日葉植物內(nèi)生菌組成分析

    在門、屬分類水平上對(duì)明日葉不同組織內(nèi)生菌進(jìn)行群落組成和相對(duì)豐度分析,由于許多物種占比小,未有代表意義,因此將相對(duì)豐度<1%的物種歸于其他類(others)。

    2. 3. 1 內(nèi)生真菌組成 700個(gè)內(nèi)生真菌OTUs比對(duì)到7個(gè)門和228個(gè)屬。在門水平上,根部菌門數(shù)最多,為6個(gè)(表1)。由圖5可知,子囊菌門(Ascomycota)在根、莖、葉和果實(shí)組織中均占主要優(yōu)勢(shì),其中在葉部占比最高;而其他類門(p_others)和擔(dān)子菌門(Basidiomycota)在果實(shí)占比均高于根、莖和葉,未知真菌門(p_unclassified Fungi)在莖中相對(duì)占比最多。在屬分類水平上,根有81個(gè)菌屬,未知真菌屬(g_unclassified Fungi)、g_unclassified Dictyosporiaceae、Serendipita、Campylospora、Codinaea、子囊菌門未分類菌屬(g_unclassified Ascomycota)和g_unclassified Tremellodendropsidales 等7個(gè)菌屬的相對(duì)豐度占比大于3.00%,其他類屬(g_others)占比5.51%;莖有93個(gè)菌屬,未知真菌屬、Paramycosphaerella、g_unclassified Trichomeriaceae、多臂菌屬(Trichomerium)、子囊菌門未分類菌屬、橫斷孢屬(Strelitziana)和枝孢屬(Cladosporium)等7個(gè)菌屬的相對(duì)豐度占比大于3.00%,其他類屬占比11.79%;葉有156個(gè)菌屬,g_unclassified Trichomeriaceae、多臂菌屬、橫斷孢屬、枝孢屬和未知真菌屬等5個(gè)菌屬的相對(duì)豐度占比大于3.00%,其他類屬占比11.99%;果實(shí)有81個(gè)菌屬,Apiotrichum、g_unclassified Trichomeriaceae、橫斷孢屬、多臂菌屬、鏈格孢屬(Alternaria)、枝孢屬和未知真菌屬等7個(gè)菌屬的相對(duì)豐度占比大于3.00%,其他類屬占比9.78%。4個(gè)組織共有內(nèi)生真菌屬32個(gè),占總屬類的14.03%(圖6)。

    2. 3. 2 內(nèi)生細(xì)菌組成 1393個(gè)內(nèi)生細(xì)菌OTUs比對(duì)到30個(gè)門和550個(gè)屬(表2)。由圖7可知,在門分類水平上,變形菌門(Proteobacteria)在根、莖、葉和果實(shí)中占主要優(yōu)勢(shì)。擬桿菌門(Bacteroidetes)、衣原體門(Chlamydiae)、綠彎菌門(Chloroflexi)、Dependentiae及其他類菌門在根中相對(duì)占比均最高;放線菌門(Actinobacteria)、異常球菌—棲熱菌門(Deinococcus-Thermus)在葉中相對(duì)占比最多;厚壁菌門(Firmicutes)在果實(shí)的相對(duì)占比均高于其他組織;而酸桿菌門(Acidobacteria)在莖中相對(duì)占比最高。在屬分類水平上,根有454個(gè)菌屬,僅羅爾斯通菌屬(Ralstonia)和短根瘤菌屬(Bradyrhizobium)2個(gè)菌屬相對(duì)豐度占比大于3.00%,其他類屬占比52.66%;莖有172個(gè)菌屬,羅爾斯通菌屬、弧菌屬(Vibrio)、紅球菌屬(Rhodococcus)、假單胞菌屬(Pseudomonas)、甲基桿菌屬(Methylobacterium)、伯克氏科未分類菌屬(g_unclassified Burkholderiaceae)、Aureimonas和鞘脂單胞菌屬(Sphingomonas)等8個(gè)菌屬的相對(duì)豐度占比大于3.00%,其他類屬占比18.54%;葉有169個(gè)菌屬,甲基桿菌屬、鞘脂單胞菌屬、異常球菌屬(Deinococcus)、胺桿菌屬(Amnibacterium)、紅球菌屬和馬賽菌屬(Massilia)6個(gè)菌屬相對(duì)豐度占比大于3.00%,其他類屬占比10.50%;果實(shí)有147個(gè)菌屬,弧菌屬、Ruegeria、甲基桿菌屬、羅爾斯通菌屬、亮桿菌屬(Leucobacter)和紅球菌屬等6個(gè)菌屬相對(duì)豐度占比大于3.00%,其他類屬占比11.81%。4個(gè)組織共有內(nèi)生細(xì)菌屬58個(gè),占總屬類的10.55%(圖8)。

    2. 3. 3 內(nèi)生菌群落結(jié)構(gòu)熱圖分析 對(duì)明日葉不同組織內(nèi)生真菌和內(nèi)生細(xì)菌測(cè)序所得的OTU結(jié)果在屬水平上進(jìn)行分類,并利用相對(duì)豐度信息繪制熱圖,選取在屬水平上前30的物種進(jìn)行聚類分析(圖9),在內(nèi)生真菌方面,葉和果實(shí)聚為一類,而對(duì)于內(nèi)生細(xì)菌,莖和果實(shí)聚為一類。根與其他3個(gè)組織的菌群物種組成和物種豐度差距均最大。

    2. 4 明日葉不同組織內(nèi)生菌群落差異性

    采用LEfSe分析方法進(jìn)行明日葉不同組織內(nèi)生菌差異顯著分析,尋找具有差異的生物標(biāo)識(shí)(Biomarker)。設(shè)置LDA閾值為4,從圖10-A可看出,除根外,Paramycosphaerella是莖中起重要作用的內(nèi)生真菌類群;散囊菌綱(Eurotiomycetes)、刺盾炱目(Chaetothyriales)、Trichomeriaceae、蒼白尾孢屬(Pallidocercospora)、多臂菌屬、擲孢酵母屬、莖點(diǎn)霉屬、亞隔孢殼科未分類菌屬(g_unclassified Didymellaceae)、煤炱目未分類菌科(f_unclassified Capnodiales)和煤炱目未分類菌屬(g_unclassified Capnodiales)是葉中起重要作用的內(nèi)生真菌類群;橫斷孢屬和刺盾炱目的fam_Incertae_sedis科是果中起重要作用的內(nèi)生真菌類群。從圖10-B可知,根中起重要作用的內(nèi)生細(xì)菌類群是擬桿菌綱(Bacteroidia)、擬桿菌門、黃單胞菌目(Xanthomonadales)、黃單胞菌科(Xanthomonadaceae)、生絲微菌科(Hyphomicrobiaceae)和Gaiellaceae;莖中起重要作用的內(nèi)生真菌類群是γ-變形菌綱(Gammaproteobacteria)、羅爾斯通菌屬、假單胞菌目(Pseudomonadales)、假單胞菌科(Pseudomonadaceae)、假單胞菌屬和黃桿菌目(Flavobacteriales);葉中起重要作用的內(nèi)生細(xì)菌類群是甲基桿菌屬、鞘脂單胞菌目(Sphingomonadales)、鞘脂單胞菌科(Sphingomonadaceae)和鞘脂單胞菌屬;果實(shí)中起重要作用的內(nèi)生細(xì)菌類群是微球菌目(Micrococcales)和亮桿菌屬。

    2. 5 明日葉內(nèi)生菌菌群功能組成

    2. 5. 1 內(nèi)生真菌 通過(guò)FUNGuild數(shù)據(jù)庫(kù)內(nèi)生真菌群落功能預(yù)測(cè)(圖11)可知,從營(yíng)養(yǎng)類型看,腐生型(57.67%)、未知(29.79%)和共生型(10.42%)在根組織中占比較大;莖中未知型(35.36%)、腐生—病理寄生型(30.89%)和共生型(13.01%)占優(yōu)勢(shì);葉中共生型(52.33%)、病理寄生型(14.49%)和腐生—共生—病理寄生型(9.17%)占比較高;果實(shí)中共生型(28.04%)、腐生型(26.19%)、病理寄生型(16.96%)和腐生—共生—病理寄生型(12.21%)具有優(yōu)勢(shì)。根據(jù)Guild分類,又分為植物病原菌、附生真菌、植物內(nèi)生菌、土壤腐生菌、蘭花類菌根真菌、植物腐生菌—木腐菌等17種。除未知腐生菌和未定義功能外,根中菌根真菌(15.05%),莖中植物病原菌—未定義腐生菌(24.14%)和附生真菌(12.55%),葉中植物病原菌(14.48)、附生真菌(32.99%)和內(nèi)生菌(19.27%),以及果實(shí)中植物病原菌(16.36%)、附生真菌(16.78%)和土壤腐生菌(21.55%)功能菌群在各組織Guild分類中的占比均超過(guò)10.00%。

    2. 5. 2 內(nèi)生細(xì)菌 通過(guò)比對(duì)EggNOG數(shù)據(jù)庫(kù)得到24種COG功能,包括RNA加工和修飾(A),染色質(zhì)結(jié)構(gòu)和動(dòng)態(tài)變化(B),能量產(chǎn)生與轉(zhuǎn)換(C),細(xì)胞循環(huán)調(diào)控、細(xì)胞分裂、染色體分裂(D),氨基酸轉(zhuǎn)運(yùn)和代謝(E),核苷酸轉(zhuǎn)運(yùn)與代謝(F),碳水化合物的運(yùn)輸和代謝(G),輔酶轉(zhuǎn)運(yùn)與代謝(H),脂質(zhì)轉(zhuǎn)運(yùn)和代謝(I),翻譯、核糖體結(jié)構(gòu)和生物合成(J),翻譯(K),復(fù)制、重建和修復(fù)(L),細(xì)胞壁/細(xì)胞膜/膜結(jié)構(gòu)的生物合成(M),細(xì)胞活力(N),翻譯后修飾、蛋白質(zhì)周轉(zhuǎn)、伴侶(O),無(wú)機(jī)離子轉(zhuǎn)運(yùn)與代謝(P),次生代謝產(chǎn)物的合成、轉(zhuǎn)運(yùn)和代謝(Q),一般功能預(yù)測(cè)(R),未知功能基因家族(S),信號(hào)轉(zhuǎn)導(dǎo)機(jī)制(T)、細(xì)胞內(nèi)轉(zhuǎn)運(yùn)、分泌和液泡運(yùn)輸(U)、防御機(jī)制(V)、細(xì)胞外結(jié)構(gòu)(W),以及細(xì)胞骨架(Z),其豐度隨組織的不同而異。明日葉不同組織內(nèi)生細(xì)菌COG功能組成分析顯示,根中排名前10的COG功能分布是S>E>R>C>K>M>T>P>G>L,莖中S>E>R>C>K>P>T>M>G>J,葉中S>E>R>C>K>P>G>M>T>I,果實(shí)中S>E>R>K>C>G>P>M>J>T。在明日葉根(2.58%)、莖(2.59%)、葉(2.94%)、果實(shí)(2.29%)組織中Q相關(guān)基因均有表達(dá)。

    3 討論

    3. 1 不同組織中內(nèi)生菌多樣性和豐富度的差異分析

    相關(guān)研究認(rèn)為,宿主植物群落多樣性受植物種類及同種植物不同組織的結(jié)構(gòu)、功能和特有化學(xué)成分等因素影響(Deng et al.,2019)。通過(guò)對(duì)冰川棘豆內(nèi)生菌多樣性研究可知,內(nèi)生真菌和內(nèi)生細(xì)菌群落結(jié)構(gòu)和多樣性存在差異(許國(guó)琪等,2021);對(duì)蛇足石杉內(nèi)生真菌的研究(Fan et al.,2020)發(fā)現(xiàn),其葉部具有相對(duì)莖和根組織更為豐富的內(nèi)生菌群落組成。本研究中,明日葉4個(gè)不同組織均定殖豐富的內(nèi)生真菌和內(nèi)生細(xì)菌,但其分布存在器官組織差異性,葉中也表現(xiàn)出較高的內(nèi)生真菌群落多樣性和豐富度,根部最低,而根中內(nèi)生細(xì)菌群落多樣性和豐富度最高。究其原因,可能葉部除來(lái)源于莖部真菌外,較大的表面積利于來(lái)自空氣、水分或其他生物的真菌附著并進(jìn)入植物體,植物光合作用產(chǎn)生的豐富營(yíng)養(yǎng)物質(zhì)也為真菌共生創(chuàng)造優(yōu)越條件。但也有研究發(fā)現(xiàn)根中具有最高的內(nèi)生真菌豐富度和多樣性(Wang et al., 2022),可能與宿主種類、生長(zhǎng)環(huán)境不同有關(guān),還應(yīng)做進(jìn)一步研究。有報(bào)道得出中藥石蓮根部?jī)?nèi)生細(xì)菌的多樣性和豐富度顯著高于葉和莖(Huo et al.,2020),可能由于內(nèi)生細(xì)菌常以根際土壤微生物為起點(diǎn),經(jīng)根部表面和內(nèi)部開(kāi)始定殖,隨著蒸騰作用沿著維管束向莖部和葉部遷移,且微生物對(duì)根部浸染程度強(qiáng)于地面部分,從而導(dǎo)致根部?jī)?nèi)生細(xì)菌群落多樣性高于其他組織。

    3. 2 不同組織中內(nèi)生菌的群落組成分析

    通過(guò)聚類熱圖分析,明日葉的根、莖、葉和果實(shí)微生物菌群組成差異較大,但莖、葉和果實(shí)組織內(nèi)生真菌和內(nèi)生細(xì)菌群落結(jié)構(gòu)相近,由于三者發(fā)育來(lái)源均屬于地面組織,可能是三者內(nèi)生菌結(jié)構(gòu)更為類似的原因,說(shuō)明明日葉組織是影響其內(nèi)生菌種群結(jié)構(gòu)與多樣性變化因素之一。

    白蘞、蛇足石杉等藥用植物的內(nèi)生菌中子囊菌門為優(yōu)勢(shì)類群(Fan et al.,2020;李鑒濱等,2022),子囊菌門是真菌中種類最多的類群,占真菌總數(shù)的40%,可能與奠基者效應(yīng)及其遺傳特性有關(guān)(張淼等,2020),且在根際土壤、藏茶、葡萄表皮等中均占比較高。本研究中子囊菌門亦為不同組織內(nèi)生真菌優(yōu)勢(shì)菌門,該門類中大多數(shù)為腐生菌,具有較強(qiáng)的物質(zhì)分解和合成能力,在養(yǎng)分循環(huán)中起著重要作用(張愛(ài)梅等, 2021),是最易與植物互惠共生的門類,因此該菌門對(duì)明日葉不同組織的養(yǎng)分循環(huán)和代謝具有意義。高通量測(cè)序的明日葉不同組織優(yōu)勢(shì)菌屬與純培養(yǎng)內(nèi)生真菌結(jié)果有較大差異(劉雯雯等,2022),說(shuō)明傳統(tǒng)培養(yǎng)法解析植物內(nèi)生真菌群落多樣性的深度與高通量技術(shù)有較大差別。此外,莖、葉和果實(shí)共有的多臂菌屬具有促進(jìn)宿主植物產(chǎn)生藥用價(jià)值物質(zhì)的作用(Nassal et al.,2018)。僅存在于根中的Serendipita能產(chǎn)生促進(jìn)宿主互作的效應(yīng)分子(Mahdi et al.,2022)。Campylospora、橫斷孢屬和Apiotrichum等菌屬可合成有用的生物活性物質(zhì)或?qū)闹髦参锞哂兄匾纳鷳B(tài)效應(yīng)。本研究中仍有一些OTU序列在不同的分類水平上未被注釋分類,表明明日葉體內(nèi)還蘊(yùn)含有大量未知的內(nèi)生真菌,有待深入研究。

    明日葉4個(gè)組織中優(yōu)勢(shì)細(xì)菌門類為變形菌門。變形菌門是植物內(nèi)部存在的最廣泛的細(xì)菌門類,物種和遺傳多樣性豐富,具有生物固氮、病蟲害防治和土壤修復(fù)等輔助作用,莼菜嫩芽、莖和葉(張秋玉等, 2022)及甜櫻桃(趙柏霞和閆建芳,2022)不同組織內(nèi)生細(xì)菌分析均發(fā)現(xiàn),優(yōu)勢(shì)細(xì)菌門類為變形菌門,且變形菌門豐度越高,組織內(nèi)碳氮養(yǎng)分及有機(jī)質(zhì)含量就越多。放線菌門、厚壁菌門和擬桿菌門是促進(jìn)生物進(jìn)化、挖掘農(nóng)用抗生素菌群的重要種質(zhì)資源來(lái)源,其在明日葉中占比也較高。此外,本研究利用高通量測(cè)序結(jié)果注釋到的一些具有不同功能的優(yōu)勢(shì)菌或核心種,如假單胞菌屬的某些菌株可分泌生長(zhǎng)素和抗生素,作為生物防治菌(王紅瑩等, 2021)。羅爾斯通菌屬具有碳固定和能量利用的代謝特點(diǎn)。馬賽菌屬?gòu)V泛存在于水體、植物根際和葉際等環(huán)境中,具有分泌生長(zhǎng)素和酶、增強(qiáng)植物抗逆性、耐受重金屬等功能(陳翔,2021)。分布于根中的短根瘤菌屬具有固氮作用;葉中甲基桿菌屬和鞘脂單胞菌屬具有促生和抑制病原體的作用,表明組織特征是影響明日葉內(nèi)生菌群落結(jié)構(gòu)的重要因素,明日葉不同組織內(nèi)生菌的種類、豐度和功能有較大差異。同時(shí),明日葉組織中還存在大量未被注釋的微生物類群,可為挖掘促生、抗逆及醫(yī)藥等有重要應(yīng)用價(jià)值的新種提供基礎(chǔ)。

    3. 3 不同組織內(nèi)生菌功能預(yù)測(cè)分析

    內(nèi)生真菌與植物宿主的演化選擇過(guò)程中,為應(yīng)對(duì)植物不同時(shí)期、不同組織的生理狀態(tài),具備共生、腐生和寄生生存策略,并發(fā)揮自生重要的生態(tài)作用,而內(nèi)生細(xì)菌通過(guò)與宿主植物間基因的橫向轉(zhuǎn)移,參與宿主關(guān)鍵的生命活動(dòng)(孫紅敏等,2016)。本研究利用FUNGuild對(duì)內(nèi)生真菌群落營(yíng)養(yǎng)模式和功能分析發(fā)現(xiàn),不同組織中未知營(yíng)養(yǎng)型和功能型均占比較大,說(shuō)明FUNGuild數(shù)據(jù)庫(kù)解析不全面,也反映出明日葉存在大量待挖掘的內(nèi)生真菌資源。根中腐生型功能群較高,可能因?yàn)楦佑|土壤,根際土壤中腐生真菌種類多,腐生真菌進(jìn)入根中成為內(nèi)生真菌的概率更高,同時(shí)腐生真菌可產(chǎn)生氧化酶和水解酶,利于碳水化合物的分級(jí),與養(yǎng)分循環(huán)密切相關(guān)(寧琪等,2022)。莖中病理腐生型功能群較多,葉和果實(shí)中共生型功能群占比較高,可能由于植物內(nèi)生真菌在不同組織中發(fā)揮不同功能,導(dǎo)致不同組織優(yōu)勢(shì)功能類群分布不同。針對(duì)內(nèi)生細(xì)菌功能,氨基酸運(yùn)輸與代謝和能量產(chǎn)生與轉(zhuǎn)換的相關(guān)基因在各組織中均占比較高,為各組織提供基礎(chǔ)保障;葉組織無(wú)機(jī)離子運(yùn)輸與代謝相關(guān)基因比其他組織高,可能與葉活動(dòng)功能代謝快有關(guān)。次生代謝產(chǎn)物是內(nèi)生菌發(fā)揮功能作用的物質(zhì)基礎(chǔ)(李盼盼等,2018),明日葉不同組織中均存在一定比例的次生代謝產(chǎn)物生物合成、分解代謝相關(guān)基因,為微生物參與宿主植物代謝提供可能。

    4 結(jié)論

    明日葉中蘊(yùn)含豐富多樣的內(nèi)生菌群,在不同組織中的群落結(jié)構(gòu)、多樣性、優(yōu)勢(shì)菌屬及差異指示種存在保守性與差異性,組織特征差異是影響內(nèi)生菌群的重要因素;不同組織內(nèi)生菌群生態(tài)功能亦差異較大,且大量為未鑒定出種屬和未定義的功能菌群。

    參考文獻(xiàn):

    蔡媛,劉浩,孔文平,鐘燦,謝景,王勇慶,黃建華,張水寒. 2021. 多花黃精內(nèi)生菌群落結(jié)構(gòu)多樣性及其與有效成分含量相關(guān)性研究[J]. 中草藥,52(13):4023-4030. [Cai Y,Liu H,Kong W P,Zhong C,Xie J,Wang Y Q,Huang J H,Zhang S H. 2021. Study on community and diversity of endophyte and its correlation with active ingredient content in Polygonatum cyrtonema[J]. Chinese Traditional and Herbal Drugs,52(13):4023-4030.] doi:10.7501/j.issn.0253-2670.2021.13.025.

    陳美琪,王興,黎俊彥,譚綺雯,李鑒濱,姚華雄,張玉鑫,謝思韻,鄧祖軍. 2023. 基于高通量測(cè)序和組織分離法的虎耳草內(nèi)生真菌多樣性分析及其抗氧化活性研究[J]. 中草藥,54(6):1924-1934. [Chen M Q,Wang X,Li J Y,Tan Q W,Li J B,Yao H X,Zhang Y X,Xie S Y,Deng Z J. 2023. Study of endophytic fungi diversity of Saxifraga stolonifera based on high-throughput sequencing method and tissue isolation method and their antioxidant activities[J]. Chinese Traditional and Herbal Drugs,54(6):1924-1934.] doi:10.7501/j.issn.0253-2670.2023.06.025.

    陳翔. 2021. 人工濕地污水處理系統(tǒng)馬賽菌(Massilia)組成、分布特征及影響因素研究[D]. 青島:青島理工大學(xué). [Chen X. 2021. Study on the composition distribution characteristics and influencing factors of Massilia in constructed wetland wastewater treatment system[D]. Qingdao:Qingdao Technological University.]

    郭文秀,巨鳳,沈芋蓉,譚璐,李新愛(ài),郭大樂(lè),鄧赟. 2022. 川芎內(nèi)生菌Pseudeurotium ovale代謝產(chǎn)物研究[J]. 天然產(chǎn)物研究與開(kāi)發(fā),34(2):239-244. [Guo W X,Ju F,Shen Y R,Tan L,Li X A,Guo D L,Deng Y. 2022. Study on metabolites of Pseudeurotium ovale isolated from Ligusticum chuanxiong[J]. Natural Product Research and Deve lopment,34(2):239-244.] doi:10.16333/j.1001-6880. 2022.2.009.

    黃雪珍,趙龍飛. 2023. 藥用植物內(nèi)生菌對(duì)宿主植物促生作用機(jī)制研究進(jìn)展[J]. 微生物學(xué)通報(bào),50(4):1653-1665. [Huang X Z,Zhao L F. 2023. Mechanism of endophytes of medicinal plants in promoting the growth of host plants[J]. Microbiology China,50(4):1653-1665.] doi:10.13344/j.microbiol.china.221120.

    李鑒濱,張?jiān)龇?,譚綺雯,陳美琪,姚華雄,鄧祖軍. 2022. 基于高通量解析和可培養(yǎng)技術(shù)的白蘞內(nèi)生真菌類群分析及分離菌株的抗痤瘡致病菌活性[J/OL]. 微生物學(xué)通報(bào). 2022-11-24. https://kns.cnki.net/kcms/detail/11.1996.q.20221122.1551.003.html. [Li J B,Zhang Z F,Tan Q W,Chen M Q,Yao H X,Deng Z J. 2022. Taxonomic analysis of fungal endophytes from Ampelopsis japonica (Thunb.) Makino by high-throughput sequencing and culture-dependent methods and the activity of the isolated strains against acne-associated pathogens[J/OL]. Microbiology China. 2022-11-24. https://kns.cnki.net/kcms/detail/11.1996.q.20221122.1551.003.html.]

    李盼盼,袁曉龍,李金海,申國(guó)明,高林,李青誠(chéng),高加明,張鵬. 2018. 湖北煙草內(nèi)生真菌生物多樣性和種群結(jié)構(gòu)分析[J]. 微生物學(xué)報(bào),58(10):1853-1863. [Li P P,Yuan X L,Li J H,Shen G M,Gao L,Li Q C,Gao J M,Zhang P. 2018.? Biodiversity and community structure of endophy tic fungi isolated from Nicotiana tabacum L.[J]. Acta Microbiologica Sinica,58(10):1853-1863.] doi:10.13343/j.cnki.wsxb.20180125.

    劉雯雯,馮瑞章,魏琴,周萬(wàn)海,王虹,杜永華. 2022. 明日葉可培養(yǎng)內(nèi)生真菌分離鑒定及功能性篩選[J]. 食品與發(fā)酵工業(yè),48(15):137-144. [Liu W W,F(xiàn)eng R Z,Wei Q,Zhou W H,Wang H,Du Y H. 2022. Isolation,identification and functional screening of culturable endophytic fungi from Angelica keiskei[J]. Food and Fermentation Industries,48(15):137-144.] doi:10.13995/j.cnki.11-1802/ts.029994.

    寧琪,陳林,李芳,張叢志,馬東豪,蔡澤江,張佳寶. 2022. 被孢霉對(duì)土壤養(yǎng)分有效性和秸稈降解的影響[J]. 土壤學(xué)報(bào),59(1):206-217. [Ning Q,Chen L,Li F,Zhang C Z,Ma D H,Cai Z J,Zhang J B. 2022. Effects of Mortierella on nutrient availability and straw decomposition in soil[J]. Acta Pedologica Sinica,59(1):206-217.] doi:10. 11766/trxb202006020213.

    孫紅敏,魏玉珍,方曉梅,余利巖,張玉琴. 2016. 蛇足石杉內(nèi)生細(xì)菌多樣性[J]. 微生物學(xué)報(bào),56(4):614-628. [Sun H M,Wei Y Z,F(xiàn)ang X M,Yu L Y,Zhang Y Q. 2016. Diversity of endophytic bacteria isolated from Huperzia serrata[J]. Acta Microbiologica Sinica,56(4):614-628.] doi:10.13343/j.cnki.wsxb.20150275.

    王紅瑩,李春燕,宋發(fā)軍,林愛(ài)華,孟艷艷. 2021. 一株白及內(nèi)生細(xì)菌的分離鑒定及抑菌機(jī)理初步研究[J]. 中南民族大學(xué)學(xué)報(bào)(自然科學(xué)版),40(2):246-251. [Wang H Y,Li C Y,Song F J,Lin A H,Meng Y Y. 2021. Isolation and identification of an endophytic bacterium from Bletilla striata and preliminary study of antibacterial mechanism[J]. Journal of South-central Minzu University(Na tural Science Edition),40(2):246-251.] doi:10.12130/znmdzk.20210305.

    許國(guó)琪,劉怡萱,曹鵬熙,姬亞麗,李敬科,李小燕,劉星. 2021. 基于Illumina MiSeq測(cè)序技術(shù)對(duì)冰川棘豆(Oxytropis glacialis)不同組織內(nèi)生菌多樣性的研究[J]. 生態(tài)學(xué)報(bào),41(12):4993-5003. [Xu G Q,Liu Y X,Cao P X,Ji Y L,Li J K,Li X Y,Liu X. 2021. Endophytes diversity of Oxytropis glacialis at different tissues based on Illumina MiSeq sequencing[J]. Acta Ecologica Sinica,41(12):4993-5003.] doi:10.5846/stxb202005081148.

    張愛(ài)梅,殷一然,孔維寶,朱學(xué)泰,楊穎麗. 2021. 西藏沙棘5種不同組織內(nèi)生細(xì)菌多樣性[J]. 生物多樣性,29(9):1236-1244. [Zhang A M,Yin Y R,Kong W B,Zhu X T,Yang Y L. 2021. Diversity of endophytic bacteria in five types of tissues of Hippophae tibetana[J]. Biodiversity Science,29(9):1236-1244.] doi:10.17520/biods.2021034.

    張麗芳,胡海林,桂騰茸,高秀,徐晴芳,朱鈴,徐啟賀,蔡建. 2023. 火棘不同組織內(nèi)生細(xì)菌群落多樣性[J/OL]. 廣西植物. 2023-04-13. https://kns.cnki.net/kcms/detail/45.1134. q.20230411.1651.014.html. [Zhang L F,Hu H L,Gui T R,Gao X,Xu Q F,Zhu L,Xu Q H,Cai J. 2023. Diversity of endophytic bacterial community in different Pyracantha fortuneana tissues[J/OL]. Guihaia. 2023-04-13. https://kns.cnki.net/kcms/detail/45.1134.q.20230411. 1651.014.html.]

    張淼,陳裕鳳,陳龍,黃飄玲,韋露玲. 2020. 不同地區(qū)藥用植物兩面針根際土壤真菌種群多樣性差異分析[J]. 生物技術(shù)通報(bào),36(9):167-179. [Zhang M,Chen Y F,Chen L,Huang P L,Wei L L. 2020. Difference analysis of the community diversity of fungi in the rhizosphere soil of Zanthoxylum nitidum(Roxb.)DC in different regions[J]. Biotechnology Bulletin,36(9):167-179.] doi:10.13560/j.cnki.biotech.bull.1985.2020-0667.

    張秋玉,陳晶,周雪妹,龔漢雨,余光輝. 2022. 莼菜不同組織內(nèi)生細(xì)菌多樣性的比較分析[J]. 中南民族大學(xué)學(xué)報(bào)(自然科學(xué)版),41(3):285-291. [Zhang Q Y,Chen J,Zhou X M,Gong H Y,Yu G H. 2022. Diversity of entophytic bacteria in different tissues of Brasenia schreberi[J]. Journal of South-central Minzu University(Natural Scien ce Edition),41(3):285-291.] doi:10.12130/znmdzk. 20220305.

    趙柏霞,閆建芳. 2022. 高通量技術(shù)分析‘砂蜜豆’甜櫻桃不同組織內(nèi)生細(xì)菌多樣性[J]. 中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào),25(3):66-77. [Zhao B X,Yan J F. 2022. Diversity analysis of endophytic bacterial community in different tissues of ‘Summit’ sweet cherry using high-throughput sequencing[J]. Journal of Agricultural Science and Technology,25(3):66-77.] doi:10.13304/j.nykjdb.2021.0647.

    Deng Z S,Liu X D,Zhang B C,Jiao S,Qi X Y,Sun Z H,He X L,Liu Y Z,Li J,Chen K K,Lin Z X,Jiang Y Y. 2019. The root endophytic fungi community structure of Pennisetum sinese from four representative provinces in China[J]. Microorganisms,7(9):332. doi:10.3390/microorganisms7090332.

    Dong L L,Cheng R Y,Xiao L N,Wei F G,Wei G F,Xu J,Wang Y,Guo X T,Chen Z J,Chen S L. 2018. Diversity and composition of bacterial endophytes among plant parts of Panax notoginseng[J]. Chinese Medicine,13(1):41. doi:10.1186/s13020-018-0198-5.

    Fan S P,Miao L Y,Li H D,Lin A H,Song F J,Zhang P. 2020. Illumina-based analysis yields new insights into the diversity and composition of endophytic fungi in cultivated Huperzia serrata[J]. PLoS One,15(11):e0242258. doi:10.1371/journal.pone.0242258.

    Gong S X,Xu B,Gu X Z,Li W H,Yu Y,Zhang Wei,Wang Z W. 2020. Study on the effects of sugar alcohols and Angelica keiskei flour on cookie quality,antioxidant,and nutrition[J]. Cereal Chemistry,97(3):714-722. doi:10. 1002/cche.10287.

    Huo X W,Wang Y,Zhang D W,Gao T,Liu M M. 2020. Characteristics and diversity of endophytic bacteria in endangered Chinese herb Glehnia littoralis based on illumina sequencing[J]. Polish Journal of Microbiology,69(3):283-291. doi:10.33073/pjm-2020-031.

    Mahdi L K,Miyauchi S,Uhlmann C,Garrido-Oter R,Langen G,Wawra S,Niu Y L,Guan R,Robertson-Albertyn S,Bulgarelli D,Parker J E,Zuccaro A. 2022. The fungal root endophyte Serendipita vermifera displays inter-kingdom synergistic beneficial effects with the microbiota in Arabidopsis thaliana and barley[J]. ISME Journal,16(3):876-889. doi:10.1038/s41396-021-01138-y.

    Nassal D,Spohn M,Eltlbany N,Jacquiod S,Smalla K,Marhan S,Kandeler E. 2018. Effects of phosphorus-mobilizing bacteria on tomato growth and soil microbial activity[J]. Plant and Soil,427(1-2):17-37. doi:10.1007/s11104-017-3528-y.

    Wang R,Zhang Q,Ju M,Yan S,Zhang Q,Gu P. 2022. The endophytic fungi diversity,community structure,and ecological function prediction of Sophora alopecuroides in Ningxia,China[J]. Microorganisms,10(11):2099. doi:10.3390/microorganisms10112099.

    (責(zé)任編輯 鄧慧靈)

    收稿日期:2023-04-15

    基金項(xiàng)目:四川省自然科學(xué)基金項(xiàng)目(2022NSFSC0176);四川省科學(xué)技術(shù)廳應(yīng)用基礎(chǔ)項(xiàng)目(2018JY0495);宜賓市科技局高端創(chuàng)新創(chuàng)業(yè)團(tuán)隊(duì)項(xiàng)目(2021YGC03)

    通訊作者:周萬(wàn)海(1979-),https://orcid.org/0000-0002-8013-0029,博士,副教授,主要從事農(nóng)業(yè)資源開(kāi)發(fā)利用研究工作,E-mail:wanhaizhou@126.com;馮瑞章(1978-),https://orcid.org/0000-0002-1601-5843,博士,教授,主要從事微生物資源開(kāi)發(fā)利用研究工作,E-mail:ruizhangfeng@126.com

    第一作者:劉雯雯(1991-),https://orcid.org/0000-0002-8659-430X,主要從事微生物學(xué)研究工作,E-mail:965501606@qq.com

    猜你喜歡
    內(nèi)生真菌多樣性
    黃草烏根部?jī)?nèi)生真菌的分離與鑒定
    內(nèi)生真菌在中藥現(xiàn)代化中的作用
    由古典戲曲看“代言體”在中國(guó)的前世今生
    戲劇之家(2016年22期)2016-11-30 15:13:39
    淺談新時(shí)期群文輔導(dǎ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)查分析
    4株石杉?jí)A甲生產(chǎn)菌原生質(zhì)體制備與再生條件的研究
    茶樹(shù)內(nèi)生真菌研究進(jìn)展
    七月丁香在线播放| 丝瓜视频免费看黄片| 午夜久久久在线观看| 2022亚洲国产成人精品| 国产欧美亚洲国产| 亚洲精品一二三| 一级爰片在线观看| 亚洲国产欧美日韩在线播放| 精品国产国语对白av| 国产精品一区二区在线不卡| 亚洲国产欧美在线一区| av在线观看视频网站免费| 免费观看av网站的网址| 成年人免费黄色播放视频| 欧美精品高潮呻吟av久久| 女的被弄到高潮叫床怎么办| 精品久久久精品久久久| 国产精品三级大全| 日日撸夜夜添| 18在线观看网站| 欧美最新免费一区二区三区| 尾随美女入室| 精品一品国产午夜福利视频| 精品亚洲成国产av| 日本黄色片子视频| 一本久久精品| 欧美另类一区| 国产一区二区在线观看av| 大香蕉久久成人网| 99热国产这里只有精品6| .国产精品久久| 国产白丝娇喘喷水9色精品| 伦精品一区二区三区| kizo精华| 在线观看免费视频网站a站| 18禁动态无遮挡网站| 99热6这里只有精品| 卡戴珊不雅视频在线播放| 青春草亚洲视频在线观看| 黑人高潮一二区| 十分钟在线观看高清视频www| av免费观看日本| 欧美人与性动交α欧美精品济南到 | 国产一级毛片在线| 国产精品久久久久成人av| 国产高清有码在线观看视频| 免费大片18禁| av一本久久久久| 97超视频在线观看视频| 亚洲综合色惰| 国产熟女午夜一区二区三区 | 女性被躁到高潮视频| 大香蕉97超碰在线| 一区二区三区乱码不卡18| 视频中文字幕在线观看| 日韩成人av中文字幕在线观看| 欧美性感艳星| 国产亚洲午夜精品一区二区久久| av免费在线看不卡| 黑人猛操日本美女一级片| 免费av中文字幕在线| 亚洲精品色激情综合| 中文精品一卡2卡3卡4更新| 中文精品一卡2卡3卡4更新| 亚洲综合色网址| 一级a做视频免费观看| 黑人高潮一二区| 国产国拍精品亚洲av在线观看| 日本wwww免费看| 99热国产这里只有精品6| 亚洲怡红院男人天堂| 成人国产麻豆网| 26uuu在线亚洲综合色| 免费播放大片免费观看视频在线观看| 日韩不卡一区二区三区视频在线| 搡女人真爽免费视频火全软件| av福利片在线| 国产一区二区三区av在线| 日本色播在线视频| 一区二区三区精品91| 在线观看免费日韩欧美大片 | 欧美精品亚洲一区二区| 亚洲欧美一区二区三区国产| 老司机亚洲免费影院| 久久久久网色| 久久久久网色| 免费黄色在线免费观看| 亚洲人与动物交配视频| 亚洲欧美一区二区三区国产| 热re99久久精品国产66热6| 婷婷色av中文字幕| 中文字幕制服av| 久久久久久久亚洲中文字幕| 欧美97在线视频| 久久久欧美国产精品| 大香蕉久久网| 国产男女超爽视频在线观看| 全区人妻精品视频| 国产女主播在线喷水免费视频网站| 边亲边吃奶的免费视频| 久久久久久久久久久久大奶| 亚洲av欧美aⅴ国产| 日产精品乱码卡一卡2卡三| 免费高清在线观看视频在线观看| 国产精品国产三级专区第一集| 国产午夜精品久久久久久一区二区三区| 日韩人妻高清精品专区| 国产探花极品一区二区| 男男h啪啪无遮挡| 亚洲第一av免费看| 精品一区二区三卡| 国产午夜精品一二区理论片| 国产精品熟女久久久久浪| videossex国产| av电影中文网址| 午夜激情福利司机影院| 欧美亚洲日本最大视频资源| 自线自在国产av| 精品亚洲成a人片在线观看| 少妇精品久久久久久久| 亚洲国产精品999| 最近中文字幕2019免费版| 久久久久久久久久久丰满| 麻豆精品久久久久久蜜桃| 亚洲欧美一区二区三区国产| 亚洲精品成人av观看孕妇| 亚洲成人一二三区av| 国产亚洲av片在线观看秒播厂| 一级毛片aaaaaa免费看小| 国产精品人妻久久久久久| 十八禁网站网址无遮挡| 国产av精品麻豆| 男女边摸边吃奶| 亚洲av成人精品一二三区| 各种免费的搞黄视频| 黄色视频在线播放观看不卡| 97精品久久久久久久久久精品| 免费看av在线观看网站| 亚洲美女搞黄在线观看| 久久久久久久亚洲中文字幕| av有码第一页| 国产永久视频网站| 免费人妻精品一区二区三区视频| 国产精品久久久久久久电影| 国产精品久久久久久久电影| 亚洲欧洲精品一区二区精品久久久 | 麻豆成人av视频| 国产精品国产三级国产av玫瑰| 王馨瑶露胸无遮挡在线观看| 97在线视频观看| 一本久久精品| 午夜影院在线不卡| 高清午夜精品一区二区三区| 一边亲一边摸免费视频| 国产精品一区www在线观看| 国产午夜精品一二区理论片| 伊人亚洲综合成人网| 亚洲av.av天堂| 国精品久久久久久国模美| 亚洲国产精品国产精品| 母亲3免费完整高清在线观看 | 男的添女的下面高潮视频| 日韩三级伦理在线观看| 曰老女人黄片| 国产精品国产av在线观看| 国产精品 国内视频| 一级毛片电影观看| 久久久久久久国产电影| 欧美精品人与动牲交sv欧美| 天美传媒精品一区二区| 在线免费观看不下载黄p国产| 国产亚洲av片在线观看秒播厂| 制服丝袜香蕉在线| 国产综合精华液| 久久久久久久久久久丰满| 夜夜看夜夜爽夜夜摸| 2018国产大陆天天弄谢| 香蕉精品网在线| 岛国毛片在线播放| 亚洲内射少妇av| 99热6这里只有精品| 欧美+日韩+精品| 九九久久精品国产亚洲av麻豆| 亚洲精品国产av成人精品| 国产精品国产三级国产专区5o| 另类精品久久| 女性被躁到高潮视频| 精品国产国语对白av| 国产乱人偷精品视频| 在线观看人妻少妇| 国产av精品麻豆| 一本久久精品| 我要看黄色一级片免费的| 欧美日韩av久久| 女的被弄到高潮叫床怎么办| 亚洲欧美色中文字幕在线| 一区二区日韩欧美中文字幕 | 日韩成人伦理影院| 色婷婷av一区二区三区视频| 国产日韩欧美视频二区| 久久久久国产网址| 国产亚洲欧美精品永久| 精品99又大又爽又粗少妇毛片| 大陆偷拍与自拍| 午夜久久久在线观看| 久久精品久久久久久久性| 高清视频免费观看一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 久久人人爽人人片av| 国产欧美亚洲国产| 91久久精品国产一区二区三区| 欧美 日韩 精品 国产| 亚洲成色77777| 日本猛色少妇xxxxx猛交久久| 国产亚洲一区二区精品| 满18在线观看网站| 一区二区三区四区激情视频| 黄色欧美视频在线观看| 少妇人妻 视频| 国产老妇伦熟女老妇高清| av专区在线播放| 日本午夜av视频| 亚洲国产最新在线播放| 色5月婷婷丁香| 日日撸夜夜添| 亚洲国产欧美日韩在线播放| 大话2 男鬼变身卡| 精品少妇内射三级| 午夜福利,免费看| 搡女人真爽免费视频火全软件| 亚洲成人手机| 嘟嘟电影网在线观看| 大码成人一级视频| 最黄视频免费看| 久久久国产一区二区| 免费大片18禁| 丰满饥渴人妻一区二区三| 午夜影院在线不卡| 精品人妻熟女毛片av久久网站| 国产欧美亚洲国产| 你懂的网址亚洲精品在线观看| 国产极品天堂在线| 国产精品人妻久久久影院| 日韩成人av中文字幕在线观看| 久久久久网色| 欧美 日韩 精品 国产| 国产精品不卡视频一区二区| 国产免费福利视频在线观看| 国产精品一国产av| 国产男女内射视频| 欧美 亚洲 国产 日韩一| 国产爽快片一区二区三区| 国产探花极品一区二区| 男女边摸边吃奶| 国精品久久久久久国模美| 大又大粗又爽又黄少妇毛片口| 日本午夜av视频| 亚洲国产精品成人久久小说| 老司机影院毛片| 亚洲av电影在线观看一区二区三区| av线在线观看网站| 男人操女人黄网站| 熟女人妻精品中文字幕| 国产免费福利视频在线观看| 午夜91福利影院| 免费看不卡的av| 亚洲av男天堂| 欧美 日韩 精品 国产| 久久热精品热| 春色校园在线视频观看| 久久久国产精品麻豆| 欧美另类一区| 男人添女人高潮全过程视频| 99热网站在线观看| 国产精品蜜桃在线观看| 大香蕉久久成人网| 丝袜脚勾引网站| 波野结衣二区三区在线| 观看美女的网站| 精品视频人人做人人爽| 国产成人免费观看mmmm| 亚洲人成网站在线播| 久久热精品热| 欧美日韩av久久| 波野结衣二区三区在线| 亚洲精品av麻豆狂野| 考比视频在线观看| 国产一区亚洲一区在线观看| 国产在线一区二区三区精| 伊人亚洲综合成人网| 亚洲精品国产av成人精品| 免费大片黄手机在线观看| 有码 亚洲区| 十八禁高潮呻吟视频| 日韩一本色道免费dvd| 欧美日韩国产mv在线观看视频| 久久久精品免费免费高清| 一级毛片 在线播放| 欧美丝袜亚洲另类| 亚洲精品亚洲一区二区| 51国产日韩欧美| .国产精品久久| av在线播放精品| 成人午夜精彩视频在线观看| 一级片'在线观看视频| 蜜桃久久精品国产亚洲av| 国产精品国产三级专区第一集| 99九九在线精品视频| 人人妻人人爽人人添夜夜欢视频| 免费久久久久久久精品成人欧美视频 | 人人澡人人妻人| 丝袜喷水一区| 亚洲国产av影院在线观看| 秋霞在线观看毛片| 男女啪啪激烈高潮av片| 国产片特级美女逼逼视频| 国产精品.久久久| 亚洲av电影在线观看一区二区三区| 日韩视频在线欧美| 大码成人一级视频| 亚洲av成人精品一区久久| 国产欧美日韩综合在线一区二区| 嘟嘟电影网在线观看| 最新的欧美精品一区二区| 在线 av 中文字幕| 男女免费视频国产| 国内精品宾馆在线| 在线天堂最新版资源| 王馨瑶露胸无遮挡在线观看| 成人亚洲精品一区在线观看| 一二三四中文在线观看免费高清| 91精品三级在线观看| a级毛片免费高清观看在线播放| 午夜福利视频在线观看免费| 女人精品久久久久毛片| 亚洲图色成人| 亚洲怡红院男人天堂| 香蕉精品网在线| 国产午夜精品久久久久久一区二区三区| 国产成人午夜福利电影在线观看| 色哟哟·www| 久久久国产欧美日韩av| 国产欧美另类精品又又久久亚洲欧美| 一个人免费看片子| 最近最新中文字幕免费大全7| 久久久久久久大尺度免费视频| 91国产中文字幕| 日韩 亚洲 欧美在线| 免费观看在线日韩| 99久久精品国产国产毛片| 久久国产精品大桥未久av| 大香蕉久久网| 国产综合精华液| 欧美老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图 | 日韩不卡一区二区三区视频在线| 在线 av 中文字幕| 两个人免费观看高清视频| 99热全是精品| 精品一品国产午夜福利视频| 欧美日韩成人在线一区二区| av国产久精品久网站免费入址| 久久久欧美国产精品| 久久久国产一区二区| 亚洲不卡免费看| 欧美激情极品国产一区二区三区 | 欧美日韩视频精品一区| 午夜福利视频精品| 免费黄频网站在线观看国产| 成人毛片a级毛片在线播放| 亚洲一级一片aⅴ在线观看| 欧美xxxx性猛交bbbb| 五月开心婷婷网| 国产欧美日韩综合在线一区二区| 男男h啪啪无遮挡| 如何舔出高潮| 亚洲欧美一区二区三区黑人 | 国产在线免费精品| 日韩中字成人| 91精品伊人久久大香线蕉| 新久久久久国产一级毛片| 九九在线视频观看精品| 人人妻人人澡人人看| 人妻人人澡人人爽人人| 国产成人一区二区在线| 久久精品国产自在天天线| 美女福利国产在线| 一区二区日韩欧美中文字幕 | 精品人妻偷拍中文字幕| 在线观看三级黄色| 国产精品久久久久久精品电影小说| 久热这里只有精品99| 久久久久久久久大av| 蜜桃国产av成人99| 国产免费福利视频在线观看| 黑人高潮一二区| 亚洲国产欧美日韩在线播放| 欧美+日韩+精品| 亚洲av福利一区| 亚洲精品一区蜜桃| 日韩一本色道免费dvd| 亚洲综合色网址| 久久精品国产a三级三级三级| 天天躁夜夜躁狠狠久久av| 一区二区三区四区激情视频| 少妇被粗大的猛进出69影院 | 22中文网久久字幕| 高清欧美精品videossex| 欧美激情国产日韩精品一区| 夫妻午夜视频| av在线app专区| 视频区图区小说| 男人爽女人下面视频在线观看| 一级毛片黄色毛片免费观看视频| 亚洲国产精品一区三区| 成年av动漫网址| 在线观看免费日韩欧美大片 | 汤姆久久久久久久影院中文字幕| 老司机影院成人| 国产成人精品无人区| videossex国产| 色5月婷婷丁香| av一本久久久久| 国产色婷婷99| www.色视频.com| 国产一级毛片在线| 2022亚洲国产成人精品| 午夜激情av网站| 国产免费福利视频在线观看| 精品久久久久久久久av| 国产无遮挡羞羞视频在线观看| 22中文网久久字幕| 国产精品熟女久久久久浪| 亚洲久久久国产精品| 国产亚洲午夜精品一区二区久久| 午夜影院在线不卡| 亚洲精品色激情综合| 免费黄网站久久成人精品| av播播在线观看一区| 亚洲成人手机| 国产精品99久久久久久久久| 精品久久国产蜜桃| 久久久久久久久久久丰满| av福利片在线| 秋霞在线观看毛片| 国产精品.久久久| 精品少妇内射三级| 成人免费观看视频高清| 日本爱情动作片www.在线观看| 日日摸夜夜添夜夜添av毛片| 高清av免费在线| 少妇精品久久久久久久| 日日撸夜夜添| 丝袜喷水一区| 亚洲国产av新网站| 亚洲精品456在线播放app| 日韩av在线免费看完整版不卡| 插阴视频在线观看视频| 91精品一卡2卡3卡4卡| 亚洲综合色网址| 老女人水多毛片| 亚洲伊人久久精品综合| 少妇熟女欧美另类| 亚洲精品aⅴ在线观看| 中文字幕精品免费在线观看视频 | 日本黄色日本黄色录像| 爱豆传媒免费全集在线观看| 涩涩av久久男人的天堂| 国产熟女午夜一区二区三区 | 全区人妻精品视频| 啦啦啦中文免费视频观看日本| 久久亚洲国产成人精品v| 91久久精品国产一区二区三区| 久久婷婷青草| 成人国产麻豆网| 成人二区视频| 18在线观看网站| 亚洲精品成人av观看孕妇| 美女福利国产在线| 久久久久网色| 亚洲av二区三区四区| 男女边吃奶边做爰视频| 成年女人在线观看亚洲视频| 啦啦啦在线观看免费高清www| 日日啪夜夜爽| 国产成人91sexporn| 韩国av在线不卡| 日本-黄色视频高清免费观看| 久久久精品94久久精品| 亚洲av成人精品一区久久| 亚洲国产精品专区欧美| 免费黄网站久久成人精品| 国产国语露脸激情在线看| 亚洲精品国产av蜜桃| 国产欧美日韩综合在线一区二区| 国产深夜福利视频在线观看| 桃花免费在线播放| 丁香六月天网| 桃花免费在线播放| 亚洲欧美色中文字幕在线| 大码成人一级视频| 日韩精品免费视频一区二区三区 | 亚洲av中文av极速乱| 日韩 亚洲 欧美在线| 一级,二级,三级黄色视频| 少妇人妻精品综合一区二区| 欧美精品人与动牲交sv欧美| 老司机影院成人| 成人毛片60女人毛片免费| 男人爽女人下面视频在线观看| 午夜日本视频在线| 精品人妻熟女av久视频| 亚洲精品乱码久久久久久按摩| 99热全是精品| 中文欧美无线码| 国产一级毛片在线| 亚洲色图 男人天堂 中文字幕 | 一本大道久久a久久精品| 老司机影院毛片| 亚洲精品国产av成人精品| 亚洲第一区二区三区不卡| 国产亚洲精品久久久com| 夜夜看夜夜爽夜夜摸| 国内精品宾馆在线| 在线看a的网站| 国产片特级美女逼逼视频| 欧美精品一区二区大全| 水蜜桃什么品种好| 只有这里有精品99| 欧美丝袜亚洲另类| 69精品国产乱码久久久| 欧美丝袜亚洲另类| 国产一区二区在线观看日韩| 97在线人人人人妻| av免费在线看不卡| 涩涩av久久男人的天堂| 大码成人一级视频| videos熟女内射| 一本色道久久久久久精品综合| 免费少妇av软件| 亚洲精品国产色婷婷电影| 国产精品人妻久久久久久| 精品少妇久久久久久888优播| 亚洲国产精品999| 欧美日韩视频高清一区二区三区二| 日日摸夜夜添夜夜添av毛片| xxxhd国产人妻xxx| 五月开心婷婷网| 日日摸夜夜添夜夜爱| 国产高清不卡午夜福利| 99九九在线精品视频| 精品熟女少妇av免费看| 亚洲国产精品成人久久小说| 成人毛片a级毛片在线播放| a 毛片基地| 91久久精品电影网| 亚洲欧美一区二区三区国产| 制服丝袜香蕉在线| 国产片特级美女逼逼视频| 成年美女黄网站色视频大全免费 | 久久久久久久久久成人| 久久午夜综合久久蜜桃| 国产色婷婷99| 亚洲av福利一区| 久久久国产欧美日韩av| 18禁裸乳无遮挡动漫免费视频| 蜜桃久久精品国产亚洲av| 久久久久久伊人网av| 国产亚洲av片在线观看秒播厂| 国产精品国产三级国产av玫瑰| 国产日韩欧美在线精品| 久久久久久久久大av| 18禁观看日本| 女人精品久久久久毛片| 成年女人在线观看亚洲视频| 国产视频首页在线观看| 亚洲在久久综合| 久久 成人 亚洲| 熟女电影av网| 成年人午夜在线观看视频| 日韩一本色道免费dvd| 亚洲av不卡在线观看| 99热网站在线观看| 国产深夜福利视频在线观看| av免费在线看不卡| 久久精品熟女亚洲av麻豆精品| 国产视频首页在线观看| kizo精华| av天堂久久9| 色94色欧美一区二区| 精品久久久精品久久久| 91午夜精品亚洲一区二区三区| 久久久精品区二区三区| 久久精品久久久久久噜噜老黄| 丝瓜视频免费看黄片| 黑人高潮一二区| 精品人妻在线不人妻| 成人18禁高潮啪啪吃奶动态图 | 亚洲国产精品一区二区三区在线| 免费日韩欧美在线观看| a级毛片在线看网站| 少妇人妻精品综合一区二区| 欧美bdsm另类| 蜜桃久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 尾随美女入室| 好男人视频免费观看在线| 高清黄色对白视频在线免费看| 国产精品99久久99久久久不卡 | 欧美日韩国产mv在线观看视频| 亚洲欧美一区二区三区黑人 | 亚洲综合色网址| a级片在线免费高清观看视频| 色视频在线一区二区三区| 国产在线免费精品|