• <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)展
    不卡av一区二区三区| 一本—道久久a久久精品蜜桃钙片| 亚洲激情五月婷婷啪啪| 亚洲熟女精品中文字幕| 成人国产麻豆网| 亚洲欧美成人综合另类久久久| 桃花免费在线播放| 色综合欧美亚洲国产小说| www.熟女人妻精品国产| 高清不卡的av网站| 久久国产亚洲av麻豆专区| 丰满少妇做爰视频| 免费日韩欧美在线观看| 午夜激情av网站| 老汉色∧v一级毛片| 中文字幕制服av| 日本91视频免费播放| 秋霞伦理黄片| 亚洲久久久国产精品| 亚洲成人av在线免费| 欧美精品人与动牲交sv欧美| 操出白浆在线播放| 精品一区二区免费观看| 国产一区二区 视频在线| 日日摸夜夜添夜夜爱| 我的亚洲天堂| 啦啦啦在线观看免费高清www| 精品久久久久久电影网| 少妇人妻 视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩视频精品一区| 国产精品一国产av| 久久性视频一级片| 久久精品亚洲av国产电影网| 麻豆av在线久日| 日韩av不卡免费在线播放| 久久99精品国语久久久| 日本午夜av视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧洲国产日韩| 亚洲av中文av极速乱| 国产精品 国内视频| 午夜日本视频在线| 自拍欧美九色日韩亚洲蝌蚪91| 国语对白做爰xxxⅹ性视频网站| 熟妇人妻不卡中文字幕| 九草在线视频观看| 街头女战士在线观看网站| 青春草亚洲视频在线观看| 久久久久精品性色| 亚洲欧美一区二区三区久久| 九色亚洲精品在线播放| 黄色怎么调成土黄色| 国产97色在线日韩免费| 亚洲伊人久久精品综合| 69精品国产乱码久久久| 日韩不卡一区二区三区视频在线| 色婷婷av一区二区三区视频| www.精华液| 女性生殖器流出的白浆| 美女高潮到喷水免费观看| av有码第一页| 欧美最新免费一区二区三区| 日韩av在线免费看完整版不卡| 日韩伦理黄色片| 日韩一卡2卡3卡4卡2021年| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产av新网站| 国产视频首页在线观看| 精品少妇内射三级| 国产一级毛片在线| 两个人免费观看高清视频| 亚洲欧美精品自产自拍| 国产精品无大码| 一边亲一边摸免费视频| 亚洲自偷自拍图片 自拍| 18禁观看日本| 国产一区二区三区综合在线观看| 两性夫妻黄色片| 少妇的丰满在线观看| 青春草亚洲视频在线观看| 男人爽女人下面视频在线观看| 老司机靠b影院| h视频一区二区三区| 久久人妻熟女aⅴ| 日韩制服骚丝袜av| 色婷婷av一区二区三区视频| 精品少妇内射三级| 天天躁夜夜躁狠狠躁躁| 亚洲精品一二三| 国产成人精品久久久久久| 在线观看免费午夜福利视频| 国产高清国产精品国产三级| 国产免费一区二区三区四区乱码| 如何舔出高潮| 美女中出高潮动态图| 国产高清不卡午夜福利| 99久久精品国产亚洲精品| 国产99久久九九免费精品| 建设人人有责人人尽责人人享有的| 欧美黑人精品巨大| 母亲3免费完整高清在线观看| 亚洲精品国产一区二区精华液| av不卡在线播放| 免费黄网站久久成人精品| 久久精品国产亚洲av涩爱| 午夜av观看不卡| 精品一区二区三区av网在线观看 | 一级毛片电影观看| 欧美人与性动交α欧美软件| 女人高潮潮喷娇喘18禁视频| 国产在线免费精品| 亚洲精品,欧美精品| 免费av中文字幕在线| 9热在线视频观看99| 无遮挡黄片免费观看| 最近中文字幕2019免费版| 又大又黄又爽视频免费| 中国国产av一级| 欧美日韩精品网址| 天天影视国产精品| 国产精品一区二区在线观看99| 日本av免费视频播放| 国产日韩一区二区三区精品不卡| av国产久精品久网站免费入址| 九九爱精品视频在线观看| 美女高潮到喷水免费观看| 亚洲欧洲国产日韩| 亚洲色图综合在线观看| 性色av一级| 尾随美女入室| 国产一区二区三区av在线| 亚洲,一卡二卡三卡| 丰满饥渴人妻一区二区三| 赤兔流量卡办理| 日韩大码丰满熟妇| 91老司机精品| 亚洲五月色婷婷综合| 一区二区日韩欧美中文字幕| av网站免费在线观看视频| 国产男女内射视频| 91精品伊人久久大香线蕉| 精品国产超薄肉色丝袜足j| 午夜福利乱码中文字幕| av网站免费在线观看视频| 亚洲精品国产一区二区精华液| 咕卡用的链子| 国产精品蜜桃在线观看| 精品人妻熟女毛片av久久网站| 久久久亚洲精品成人影院| 久久久国产一区二区| 亚洲自偷自拍图片 自拍| 久久久国产欧美日韩av| 天天添夜夜摸| 日韩大片免费观看网站| 少妇猛男粗大的猛烈进出视频| 免费黄网站久久成人精品| 99re6热这里在线精品视频| 婷婷色麻豆天堂久久| xxx大片免费视频| 一本大道久久a久久精品| 嫩草影视91久久| 亚洲中文av在线| 久久久久久久久久久久大奶| 少妇人妻久久综合中文| 精品一区在线观看国产| 国产1区2区3区精品| 久久久久国产精品人妻一区二区| 99国产精品免费福利视频| 国产麻豆69| av女优亚洲男人天堂| 91精品三级在线观看| 亚洲免费av在线视频| 欧美精品亚洲一区二区| 亚洲人成77777在线视频| 亚洲av福利一区| 久久精品久久久久久久性| 久久国产精品大桥未久av| 国产亚洲欧美精品永久| 国产精品久久久久久精品古装| 午夜福利在线免费观看网站| 国产激情久久老熟女| 97精品久久久久久久久久精品| 超碰97精品在线观看| 人人澡人人妻人| 在线 av 中文字幕| 高清在线视频一区二区三区| 高清不卡的av网站| 亚洲综合色网址| 国产精品 欧美亚洲| 两个人看的免费小视频| 咕卡用的链子| 三上悠亚av全集在线观看| 国产成人午夜福利电影在线观看| 丝袜美腿诱惑在线| 桃花免费在线播放| 建设人人有责人人尽责人人享有的| 亚洲人成77777在线视频| 一级毛片电影观看| 亚洲av成人不卡在线观看播放网 | 午夜福利乱码中文字幕| av片东京热男人的天堂| 蜜桃国产av成人99| av.在线天堂| 在线观看国产h片| 九九爱精品视频在线观看| 美女大奶头黄色视频| 卡戴珊不雅视频在线播放| 男女无遮挡免费网站观看| 观看av在线不卡| 2021少妇久久久久久久久久久| 国产一卡二卡三卡精品 | 亚洲av综合色区一区| 性高湖久久久久久久久免费观看| 国产成人91sexporn| 色视频在线一区二区三区| 深夜精品福利| 亚洲精品视频女| 两个人免费观看高清视频| a级片在线免费高清观看视频| 纵有疾风起免费观看全集完整版| av免费观看日本| 久久精品久久精品一区二区三区| 亚洲国产欧美日韩在线播放| 免费看不卡的av| 伦理电影大哥的女人| 亚洲人成电影观看| 亚洲一级一片aⅴ在线观看| 国产av国产精品国产| 韩国高清视频一区二区三区| 一级爰片在线观看| 久久性视频一级片| 久久久久国产精品人妻一区二区| 制服丝袜香蕉在线| 国产免费一区二区三区四区乱码| 亚洲国产日韩一区二区| 亚洲av电影在线观看一区二区三区| 涩涩av久久男人的天堂| 亚洲熟女精品中文字幕| 亚洲欧美中文字幕日韩二区| 免费高清在线观看日韩| 亚洲精品国产一区二区精华液| 午夜激情av网站| tube8黄色片| 日本爱情动作片www.在线观看| 久久精品国产综合久久久| 精品亚洲成a人片在线观看| 久久人妻熟女aⅴ| 亚洲综合色网址| 色婷婷久久久亚洲欧美| 男女无遮挡免费网站观看| 欧美日韩综合久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 超碰成人久久| 青春草视频在线免费观看| 欧美在线一区亚洲| 美女扒开内裤让男人捅视频| a 毛片基地| 国产伦理片在线播放av一区| 成人黄色视频免费在线看| 欧美久久黑人一区二区| a 毛片基地| 青春草视频在线免费观看| 亚洲国产精品一区二区三区在线| 美女扒开内裤让男人捅视频| 国产日韩欧美视频二区| 亚洲精品一二三| 一级毛片 在线播放| 国产无遮挡羞羞视频在线观看| 永久免费av网站大全| 大片免费播放器 马上看| 人妻一区二区av| 国产精品免费大片| 国产乱人偷精品视频| 亚洲天堂av无毛| 狂野欧美激情性bbbbbb| 亚洲成人av在线免费| 精品少妇久久久久久888优播| 久久精品熟女亚洲av麻豆精品| 性少妇av在线| 国产成人欧美| 日韩人妻精品一区2区三区| 五月天丁香电影| 夫妻性生交免费视频一级片| 亚洲,一卡二卡三卡| 亚洲精华国产精华液的使用体验| 国产精品免费大片| 日韩电影二区| 国产97色在线日韩免费| av在线app专区| 日本色播在线视频| 丰满饥渴人妻一区二区三| 美女大奶头黄色视频| 无限看片的www在线观看| 欧美日本中文国产一区发布| 91老司机精品| 国产精品久久久人人做人人爽| 国产精品人妻久久久影院| 天天躁夜夜躁狠狠躁躁| 波多野结衣一区麻豆| 亚洲精品成人av观看孕妇| 男女下面插进去视频免费观看| 一区二区三区乱码不卡18| 不卡视频在线观看欧美| 制服人妻中文乱码| 另类亚洲欧美激情| 亚洲精品国产色婷婷电影| 一级爰片在线观看| 久久久久国产一级毛片高清牌| 男人舔女人的私密视频| 视频区图区小说| 黄色 视频免费看| 亚洲精品视频女| 曰老女人黄片| 亚洲七黄色美女视频| 欧美日韩国产mv在线观看视频| 国产探花极品一区二区| 久久精品久久精品一区二区三区| 久久久久精品人妻al黑| av国产精品久久久久影院| 亚洲欧美成人精品一区二区| 久久影院123| 国产av一区二区精品久久| 中文精品一卡2卡3卡4更新| 精品人妻一区二区三区麻豆| 一级毛片 在线播放| 久久精品国产a三级三级三级| 美女主播在线视频| 999久久久国产精品视频| 国产一区二区三区av在线| 麻豆乱淫一区二区| 欧美亚洲 丝袜 人妻 在线| 老司机在亚洲福利影院| 母亲3免费完整高清在线观看| 在线 av 中文字幕| 伊人亚洲综合成人网| 国产片特级美女逼逼视频| 亚洲欧美精品自产自拍| 一边摸一边抽搐一进一出视频| 欧美日韩一区二区视频在线观看视频在线| 高清欧美精品videossex| 色网站视频免费| 亚洲国产欧美网| 国产在线免费精品| 99久久综合免费| videosex国产| 毛片一级片免费看久久久久| 黑丝袜美女国产一区| 嫩草影院入口| 国产一区二区激情短视频 | 69精品国产乱码久久久| 三上悠亚av全集在线观看| 别揉我奶头~嗯~啊~动态视频 | 国产精品免费视频内射| 丝袜人妻中文字幕| 建设人人有责人人尽责人人享有的| 18在线观看网站| 国产 精品1| 老鸭窝网址在线观看| 久久久久网色| 亚洲欧美精品自产自拍| 久久久久精品性色| 婷婷色麻豆天堂久久| 亚洲精品成人av观看孕妇| 亚洲精品日韩在线中文字幕| 一边摸一边做爽爽视频免费| 国产精品一国产av| 999久久久国产精品视频| bbb黄色大片| 在线亚洲精品国产二区图片欧美| 国产精品秋霞免费鲁丝片| 免费黄网站久久成人精品| 成人国产麻豆网| 日韩制服丝袜自拍偷拍| 操美女的视频在线观看| 国产日韩欧美在线精品| 午夜激情久久久久久久| 午夜福利乱码中文字幕| 成人18禁高潮啪啪吃奶动态图| 亚洲第一青青草原| 午夜久久久在线观看| avwww免费| 亚洲欧洲日产国产| 国产又色又爽无遮挡免| 校园人妻丝袜中文字幕| 看十八女毛片水多多多| 最近手机中文字幕大全| 男女免费视频国产| 最近的中文字幕免费完整| 这个男人来自地球电影免费观看 | 精品一区二区免费观看| 午夜影院在线不卡| 丰满乱子伦码专区| 91精品三级在线观看| 久久精品久久精品一区二区三区| 波野结衣二区三区在线| av电影中文网址| 天堂中文最新版在线下载| 不卡av一区二区三区| 成人毛片60女人毛片免费| 黄色视频在线播放观看不卡| 国产精品偷伦视频观看了| 国产一级毛片在线| 亚洲伊人久久精品综合| 精品一区二区三卡| 国产高清国产精品国产三级| 亚洲精品美女久久av网站| 午夜福利视频精品| 黄频高清免费视频| 最近的中文字幕免费完整| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产区一区二| 成人国产麻豆网| 涩涩av久久男人的天堂| 国产男女内射视频| 一本一本久久a久久精品综合妖精| 久久久久久久久免费视频了| 欧美激情极品国产一区二区三区| 人妻人人澡人人爽人人| 国产成人精品久久二区二区91 | 黄频高清免费视频| 97精品久久久久久久久久精品| 欧美日韩成人在线一区二区| 悠悠久久av| 久久久久精品久久久久真实原创| 天天躁狠狠躁夜夜躁狠狠躁| 丰满乱子伦码专区| 日韩,欧美,国产一区二区三区| 亚洲成人手机| 狠狠婷婷综合久久久久久88av| 美女国产高潮福利片在线看| 91国产中文字幕| 亚洲精品国产av蜜桃| 日日撸夜夜添| 日本欧美视频一区| 国产片内射在线| 久久精品国产亚洲av涩爱| 中文字幕高清在线视频| 亚洲av日韩在线播放| 久久久久国产精品人妻一区二区| 一区在线观看完整版| 亚洲伊人色综图| xxx大片免费视频| 另类亚洲欧美激情| 啦啦啦视频在线资源免费观看| 波野结衣二区三区在线| 曰老女人黄片| 黄片无遮挡物在线观看| 亚洲国产精品国产精品| 熟女少妇亚洲综合色aaa.| 国产不卡av网站在线观看| 久久精品久久精品一区二区三区| 精品国产乱码久久久久久男人| 老司机在亚洲福利影院| 亚洲熟女精品中文字幕| 丁香六月天网| 男女高潮啪啪啪动态图| 国产在线免费精品| 在现免费观看毛片| 青草久久国产| 久久久久久人妻| 中国三级夫妇交换| 亚洲精华国产精华液的使用体验| 日日摸夜夜添夜夜爱| 狠狠婷婷综合久久久久久88av| 国产黄色视频一区二区在线观看| 日韩精品有码人妻一区| 国产日韩欧美视频二区| 午夜老司机福利片| 天天影视国产精品| 久久久久久久久免费视频了| 99精国产麻豆久久婷婷| 久久久久国产精品人妻一区二区| 操出白浆在线播放| 新久久久久国产一级毛片| 高清av免费在线| 国产精品二区激情视频| av电影中文网址| 亚洲激情五月婷婷啪啪| 久久久久久人妻| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精华国产精华液的使用体验| 中文字幕高清在线视频| 欧美精品人与动牲交sv欧美| 女的被弄到高潮叫床怎么办| 十分钟在线观看高清视频www| 考比视频在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲七黄色美女视频| 黄网站色视频无遮挡免费观看| 午夜老司机福利片| 日韩欧美一区视频在线观看| 日日爽夜夜爽网站| 天堂中文最新版在线下载| 亚洲熟女毛片儿| 国产极品粉嫩免费观看在线| 亚洲国产看品久久| 天美传媒精品一区二区| 精品免费久久久久久久清纯 | 中文字幕人妻丝袜一区二区 | 色精品久久人妻99蜜桃| 狠狠精品人妻久久久久久综合| 亚洲国产欧美在线一区| 99国产精品免费福利视频| 美国免费a级毛片| 久久免费观看电影| 自线自在国产av| 成人国产av品久久久| 欧美中文综合在线视频| 成人午夜精彩视频在线观看| 久久韩国三级中文字幕| 电影成人av| 天天躁夜夜躁狠狠躁躁| 国产成人啪精品午夜网站| 在线天堂中文资源库| 在现免费观看毛片| av免费观看日本| 日本欧美视频一区| 久久青草综合色| 国产成人精品久久久久久| 婷婷色av中文字幕| 亚洲av福利一区| 日本黄色日本黄色录像| 久久久久网色| 免费在线观看黄色视频的| 日韩电影二区| 卡戴珊不雅视频在线播放| av网站免费在线观看视频| 在线亚洲精品国产二区图片欧美| 国产男女内射视频| 亚洲av日韩精品久久久久久密 | 另类精品久久| 国产精品人妻久久久影院| 亚洲,欧美,日韩| 一级毛片我不卡| av在线观看视频网站免费| 九九爱精品视频在线观看| 夫妻午夜视频| 色94色欧美一区二区| 久热爱精品视频在线9| 精品一区二区三区四区五区乱码 | 久久久国产精品麻豆| 九色亚洲精品在线播放| 精品一区二区三区四区五区乱码 | 波野结衣二区三区在线| 久久人人爽人人片av| 精品一区二区三区四区五区乱码 | 国产精品一区二区在线观看99| 久热爱精品视频在线9| 精品久久久久久电影网| a级片在线免费高清观看视频| av电影中文网址| 两性夫妻黄色片| 日本欧美视频一区| 黄色毛片三级朝国网站| 色精品久久人妻99蜜桃| 精品人妻一区二区三区麻豆| 黄片播放在线免费| 国产福利在线免费观看视频| 欧美成人精品欧美一级黄| 国产激情久久老熟女| 久久久久久久久久久免费av| 亚洲欧美日韩另类电影网站| 成人国产av品久久久| xxxhd国产人妻xxx| 国产免费福利视频在线观看| 国产欧美日韩一区二区三区在线| 多毛熟女@视频| 国产男女超爽视频在线观看| 一边摸一边做爽爽视频免费| 亚洲成国产人片在线观看| 亚洲精品美女久久久久99蜜臀 | 欧美变态另类bdsm刘玥| 咕卡用的链子| 国产精品无大码| 美女中出高潮动态图| 国产成人a∨麻豆精品| 汤姆久久久久久久影院中文字幕| 欧美av亚洲av综合av国产av | 高清黄色对白视频在线免费看| 午夜精品国产一区二区电影| 国产视频首页在线观看| 中文字幕人妻丝袜制服| 丝袜在线中文字幕| 黑人猛操日本美女一级片| 99国产综合亚洲精品| 亚洲精品中文字幕在线视频| 又粗又硬又长又爽又黄的视频| 成人免费观看视频高清| 黄网站色视频无遮挡免费观看| 亚洲国产欧美在线一区| videosex国产| 激情视频va一区二区三区| 久久久国产一区二区| 亚洲av电影在线进入| 久久久久久免费高清国产稀缺| 日韩av免费高清视频| 波野结衣二区三区在线| 国产无遮挡羞羞视频在线观看| 免费久久久久久久精品成人欧美视频| 欧美最新免费一区二区三区| 亚洲av综合色区一区| 多毛熟女@视频| 你懂的网址亚洲精品在线观看| 国产一级毛片在线| av一本久久久久| 欧美激情高清一区二区三区 | 高清av免费在线| 国产男女超爽视频在线观看| 国产精品久久久久久精品古装| www.精华液| 欧美激情高清一区二区三区 | 熟女少妇亚洲综合色aaa.| 黑丝袜美女国产一区| 欧美成人精品欧美一级黄|