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

    不同生境黑老虎根際和內(nèi)生真菌多樣性分析

    2023-08-20 05:13:34劉濤馬楠金吉芬彭志軍李興忠汪得凱
    廣西植物 2023年7期
    關(guān)鍵詞:根際內(nèi)生生境

    劉濤 馬楠 金吉芬 彭志軍 李興忠 汪得凱

    摘 要:? 為探究不同生長(zhǎng)條件下黑老虎根際和根部?jī)?nèi)生真菌群落組成和多樣性及其與土壤環(huán)境因子的相關(guān)性,該文應(yīng)用Illumina高通量測(cè)序方法對(duì)貴州3個(gè)不同生境下黑老虎根際和根部?jī)?nèi)生真菌進(jìn)行了研究。結(jié)果表明:(1)3種生境下,根際土壤真菌OTU數(shù)量(3 867)遠(yuǎn)多于根部?jī)?nèi)生真菌(801),其中根際土壤真菌共有的OTU為72個(gè),共注釋到5個(gè)門、49個(gè)屬,大多為子囊菌門;屬水平上被孢霉屬、外瓶柄霉屬、 柱孢屬占比較高;根部?jī)?nèi)生真菌共有的OTU為14個(gè),共注釋到2個(gè)門、11個(gè)屬,子囊菌門(13個(gè),占比92.9%)占絕對(duì)優(yōu)勢(shì),屬水平上被孢霉屬、外瓶柄霉屬、 柱孢屬和叢赤殼屬占比最高; 所有樣本中,共有的OTU僅為6個(gè),注釋到2個(gè)門、5個(gè)屬,子囊菌門(5個(gè),占比83.3%)為優(yōu)勢(shì)門;在屬水平上,占比最高的為外瓶柄霉屬(2個(gè),33.3%),其余分別為被孢霉屬、柱孢屬和叢赤殼屬 。Alpha多樣性分析表明,根際土壤的真菌群落多樣性和豐富度均顯著高于根部?jī)?nèi)生真菌,而野生生境的真菌多樣性高于栽培生境。(2)在門水平上,3個(gè)生境下主要內(nèi)生真菌類群均為子囊菌門(Ascomycota)和擔(dān)子菌門(Basidiomycota),占總菌群的88.28%;在屬水平上,不同生境條件下,根際和根部?jī)?nèi)生真菌群落結(jié)構(gòu)差異明顯;栽培生境下,根部?jī)?nèi)生真菌菌群具有一定的偏好性,而野生生境下,根部?jī)?nèi)生真菌菌群均勻度更高;FUNGuild真菌群落功能預(yù)測(cè)顯示,栽培生境下的病理-腐生營(yíng)養(yǎng)型(pathotroph-saprotroph) 在根部?jī)?nèi)生真菌中占比較高,而野生生境下的腐生營(yíng)養(yǎng)型(saprotroph)及共生營(yíng)養(yǎng)型(symbiotroph)占比較高。(3)土壤環(huán)境因子對(duì)根部?jī)?nèi)生真菌和根際真菌的影響方式不同,其中土壤總鉀(TK)和土壤總磷(TP)與黑老虎根部?jī)?nèi)生真菌香農(nóng)指數(shù)和辛普森指數(shù)顯著正相關(guān),而土壤有機(jī)質(zhì)(SOM)、總氮(TN)和速效氮(AN)與黑老虎根際土壤真菌Ace指數(shù)和Chao1指數(shù)顯著正相關(guān)。綜上表明,SOM、TN、AN是影響黑老虎根際土壤真菌群落的主要土壤環(huán)境因子。

    關(guān)鍵詞: 黑老虎, ITS, Illumina高通量測(cè)序, 真菌群落, 土壤環(huán)境因子, 多樣性

    中圖分類號(hào):? Q939

    文獻(xiàn)標(biāo)識(shí)碼:? A

    文章編號(hào):? 1000-3142(2023)07-1222-14

    收稿日期:? 2022-12-16

    基金項(xiàng)目:? 浙江省重點(diǎn)研發(fā)項(xiàng)目(2021C02043); 浙江理工大學(xué)科研啟動(dòng)基金(19042142-Y)。

    第一作者: 劉濤(1982-),碩士,副研究員,研究方向?yàn)楣麡湓耘嘤N,(E-mail)737689429@qq.com。

    通信作者:? 汪得凱,博士,研究員,研究方向?yàn)樗幱弥参镔Y源,(E-mail)kay77@163.com。

    Diversity of rhizosphere and root endophytic fungi in

    Kadsura coccinea under different habitats

    LIU Tao1, MA Nan2, JIN Jifen3, PENG Zhijun1, LI Xinzhong1, WANG Dekai2*

    ( 1. Citrus Fruit Research Institute in Guizhou Province, Guiyang 550100, China; 2. College of Life Sciences and Medicine, Zhejiang Sci-Tech University,

    Hangzhou 310018, China; 3. Institute of Crop Germplasm Resources, Guizhou Academy of Agricultural Sciences, Guiyang 550100, China )

    Abstract:? Kadsura coccinea is a perennial vine medicinal plant,? and its growth, development, and authentic medicinal qualities are affected by growing habitat, rhizosphere soil, and endophytic fungi. In order to explore the composition and diversity of fungal community in K. coccinea under different habitats, as well as their correlations with soil environmental factors, the rhizosphere and root endophytic fungi of K. coccinea from three different habitats in Guizhou were studied based on Illumina high-throughput sequencing method. The results were as follows: (1) Under three habitats, the number of OTU detected in rhizosphere soil (3 867) was much more than that in root endophyte (801). The common OTU of rhizosphere soil was 72, which were annotated to 5 phyla and 49 genera, most of which were ascomycetes. At the genus level, the proportions of Mortierella, Exophiala, and Cylindrocarpon were relatively high. The common OTU of root endophytic fungi was 14, which were annotated to 2 phylums and 11 genera. Ascomycetes (13, accounting for 92.9%) were overwhelmingly dominant. At the genus level, Mortierella, Exophiala, Cylindrocarpon and Nectria were dominant genra. There were only 6 OTU in total, annotated to 2 phylums and 5 genera, and the ascomycetes phylum (5, accounting for 83.3%) were the dominant phylum. At the genus level, Exophiala (2, 33.3%) accounted for the highest proportion? and the rest were Mortierella, Cylindrocarpon and Nectria, respectively. Alpha diversity analysis showed that the diversity and richness of fungal communities in rhizosphere soil were significantly higher than those in root endophytic fungi, while the diversity of fungi under wild habits was higher than that under cultivated habits. (2) At the phylum level, the main endophytic fungal groups under three habitats were Ascomycota and Basidiomycota, accounting for 88.28% of the total flora. At the genus level, the community structures of endophytic fungi in rhizosphere and root were significantly different under different habitats. The endophytic fungi flora in the root had a certain preference under the cultivation habitat, while the evenness under the wild habitat? was higher than that under cultivation habitat. The prediction functions of fungal community based on FUNGuild showed that pathotroph-saprotroph type in the root endophytic fungi was relatively high under the cultural habit, while the proportion of saprophytic type and symbiotroph type was relatively high under the wild habit. (3) Soil environmental factors had different effects on root endophytic fungi and rhizosphere fungi in K. coccinea. Total potassium (TK) and total phosphorus (TP) in soil were positively correlated with Shannon index and Simpson index of endophytic fungi in K. coccinea? root; whereas, soil organic matter (SOM), total nitrogen (TN) and available nitrogen (AN) were positively correlated with Ace index and Chao1 index of rhizosphere soil fungi. In conclusion, SOM, TN and AN are the main soil environmental factors affecting the rhizosphere soil fungal community of the K. coccinea.

    Key words: Kadsura coccinea, ITS, Illumina high-throughput sequencing, fungal community, soil environmental factors, diversity

    黑老虎 (Kadsura coccinea)隸屬五味子科(Schisandraceae)南五味子屬(Kadsura),其果實(shí)富含多種氨基酸、維生素C、粗蛋白質(zhì)、可溶性糖等營(yíng)養(yǎng)成分,是一種歷史悠久的民間藥材植物。近年來,黑老虎因豐富的營(yíng)養(yǎng)和獨(dú)特的風(fēng)味而深受消費(fèi)者的歡迎,在貴州黔東南、湖北恩施及湖南通道等南方地區(qū),作為一種新型特色水果,成為脫貧致富的林下產(chǎn)業(yè)(劉濤等,2009;高漸飛等,2022)。黑老虎在醫(yī)藥、保健食品及食品領(lǐng)域展現(xiàn)出廣闊的應(yīng)用前景(Sritalahareuthai et al., 2020;王麗軍等,2021)。

    植物宿主與內(nèi)生真菌、根際土壤真菌交互作用,形成復(fù)雜的調(diào)控網(wǎng)絡(luò),共同維持生態(tài)系統(tǒng)的動(dòng)態(tài)平衡(Jia et al., 2016;Lu et al., 2021)。植物內(nèi)生真菌生活在宿主體內(nèi),不會(huì)引起宿主明顯病癥,與宿主協(xié)同進(jìn)化,利用宿主營(yíng)養(yǎng)供自身需要,同時(shí)在植物生長(zhǎng)發(fā)育、增強(qiáng)逆境抵抗能力和生物系統(tǒng)演化等方面發(fā)揮重要作用(Huang et al., 2019)。此外,內(nèi)生真菌在與宿主的長(zhǎng)期相互作用中,通過影響和調(diào)節(jié)宿主次生代謝產(chǎn)物的合成及累積影響宿主活性物質(zhì)的分布和含量,或可以產(chǎn)生與宿主植物類似的次生代謝物,廣泛應(yīng)用于農(nóng)學(xué)、醫(yī)學(xué)、工業(yè)和生物技術(shù)等領(lǐng)域(Tiwari & Bae, 2022)。根際真菌一方面利用植物中的碳水化合物進(jìn)行生長(zhǎng)發(fā)育,在有機(jī)質(zhì)轉(zhuǎn)化、碳-氮循環(huán)和促進(jìn)土壤環(huán)境健康等方面發(fā)揮重要作用,另一方面其自身新陳代謝所產(chǎn)生的產(chǎn)物反過來可以促進(jìn)或抑制植物的生長(zhǎng)(Bonfante & Anca, 2009; Fierer, 2017)。而同為木蘭科的五味子內(nèi)生真菌研究有一些報(bào)道,五味子屬植物中含有豐富的內(nèi)生真菌,從不同來源、不同組織中分離的內(nèi)生菌具有抗氧化、拮抗植物病原菌等功能(潘爭(zhēng)艷等,2007;Gao et al., 2020;張譽(yù)薺等,2021)。此外,從五味子中陸續(xù)分離的一些內(nèi)生真菌能夠產(chǎn)生與宿主植物類似的木脂素及三萜類,或具有轉(zhuǎn)化宿主次生代謝物的功能(Wang et al., 2017;Qin et al., 2019;Qin et al., 2020)。劉濤等(2023)利用高通量測(cè)序技術(shù)對(duì)黑老虎不同組織內(nèi)生真菌進(jìn)行了分析,明確了不同組織內(nèi)生真菌群落的結(jié)構(gòu)和功能特征。植物內(nèi)生真菌和根際真菌群落組成和結(jié)構(gòu)受物種、地理環(huán)境、氣候和土壤特性等多種因素影響(Patten & Glick,2002; Chaparro et al., 2013; 王艷等, 2019)。目前,對(duì)黑老虎及南五味子屬的研究多集中于營(yíng)養(yǎng)成分和藥理學(xué)研究,而對(duì)黑老虎根際真菌和內(nèi)生真菌組成和多樣性研究尚不多見。

    本研究采用ITS高通量測(cè)序技術(shù)研究了栽培和野生生境下多年生成熟黑老虎根際真菌和根部?jī)?nèi)生真菌群落結(jié)構(gòu)、組成和多樣性,并開展了真菌多樣性和土壤理化因子的相關(guān)性分析,以期為黑老虎資源保護(hù)、土壤根際微生態(tài)環(huán)境研究及有益微生物的發(fā)掘和利用提供基礎(chǔ)。

    1 材料與方法

    1.1 采樣地點(diǎn)及樣品處理

    選取貴州省3個(gè)黑老虎代表性研究樣地(3個(gè)樣地)。貴陽(GY)樣地,貴陽市貴州果樹研究所黑老虎種質(zhì)資源圃(106°39′38.33″ E、26°29′59.00″ N),海拔1 051 m,為多年生人工管理種質(zhì)資源圃,土壤類型為中厚層紅黃壤;榕江(RJ)樣地,黔東南苗族侗族自治州榕江縣古州鎮(zhèn)(108°29′15.18″ E、25°55′52.43″ N),海拔361 m,為野生樣地,土壤類型為中厚層紅黃壤為主,富含有機(jī)質(zhì);福泉(FQ)樣地,黔南布依族苗族自治州福泉市金山街道(107°24′15.48″ E、26°44′36.62″ N),海拔928 m,為野生樣地,土壤類型以石灰石為主,富含P,福泉有“亞洲磷都”之美譽(yù)。3個(gè)樣地分別代表不同栽培生境、不同土壤類型和不同酸堿度,代表了貴州黑老虎生長(zhǎng)的主要生境。

    于果實(shí)采收季節(jié),分別在3個(gè)樣地選取生長(zhǎng)多年、藤蔓高度木質(zhì)化的健康植株,用鏟子去除地表植被和雜質(zhì),使用小鋤頭輕輕刨開土壤,尋找主根周邊的細(xì)根(直徑≤2 mm),挖取少量帶土的細(xì)根(Berhongaray et al., 2013)。每3株樣品合并為1個(gè)生物學(xué)重復(fù),每個(gè)樣本設(shè)置3個(gè)生物學(xué)重復(fù),將樣品收集至無菌采樣袋中并密封保存,迅速置于冷鏈箱內(nèi),于4 h內(nèi)帶回實(shí)驗(yàn)室處理。

    將黑老虎根部的土壤抖落,留存附著在根部無法抖落的為根際土(rhizosphere soil,Rh)(Edwards et al., 2015),將根系樣品轉(zhuǎn)移至無菌50 mL離心管中(含20 mL無菌10 mmol·L-1 PBS溶液),置于搖床 (蘇州培英實(shí)驗(yàn)設(shè)備有限公司,培英,型號(hào):THZ-C-1),120 r·min-1,室溫下振蕩20 min(Yu et al., 2013)。使用無菌鑷子挑除根系,剩余懸浮液高速離心 (6 000 g,4 ℃) 20 min,沉淀即為根際土樣品,液氮速凍后,置于-80 ℃冰箱備用。剩余的根系樣品(root,R)用自來水沖去表面多余的根際土,進(jìn)行表面消毒(75%乙醇浸泡2 min,5%次氯酸處理3 min,無菌水清洗3次)獲得表面消毒的根組織樣品,液氮速凍后,置于-80 ℃冰箱用于DNA 提取。

    1.2 土壤理化指標(biāo)測(cè)定

    土壤pH值測(cè)定:稱取3 g土壤樣品于50 mL燒杯中,加入7.5 mL去除CO2的蒸餾水(水土比2.5∶1),用水平振蕩器劇烈振蕩2 min,使土粒充分分散后,放置30 min后用pH計(jì)測(cè)定,1 h內(nèi)完成測(cè)定。每個(gè)樣品測(cè)定3次,計(jì)算平均值。土壤總氮(total nitrogen, TN)采用半微量凱氏定氮法,總磷(total phosphorus, TP)采用NaOH堿熔-鉬銻抗分光光度法,土壤總鉀(total potassium, TK)測(cè)定采用NaOH熔融法,土壤速效氮(avaliable nitrogen, AN)的測(cè)定采用堿解擴(kuò)散法,土壤速效磷(available phosphorus, AP)的測(cè)定采用分光光度法,土壤速效鉀(available potassium, AK)的測(cè)定采用乙酸銨浸提法,土壤有機(jī)質(zhì)(soil organic matter, SOM)的測(cè)定采用重鉻酸鉀容量法。每個(gè)樣品重復(fù)3次。

    1.3 DNA提取及高通量測(cè)序

    根系樣品DNA提取采用植物基因組DNA提取試劑盒(貨號(hào):DP305,天根生化科技有限公司,北京)。根際土壤樣品DNA提取采用土壤DNA提取試劑盒(貨號(hào):DP336,天根生化科技有限公司,北京)。ITS引物:fITS7(5′-GTGARTCATCGAA

    TCTTTG-3′)和ITS4(5′-TCCTCCGCTTATTGATATG

    C-3′)。反應(yīng)體系 (25 μL):2×Phanta Max master mix 12.5 μL,正、反向引物(1 μmol·L-1)各2.5 μL,DNA μL (50 ng),ddH2O 6.5 μL。PCR 反應(yīng)條件:94 ℃/1 min-(94 ℃/10 s-50 ℃/30 s-72 ℃/45 s) ×32個(gè)循環(huán)-72 ℃/10 min。PCR產(chǎn)物經(jīng)2%瓊脂糖凝膠電泳確證,用AMPure XT beads (Beckman Coulter Genomics, Danvers, MA, USA)純化,Qubit(Invitrogen, USA)定量。擴(kuò)增子文庫的大小和數(shù)量分別用Agilent 2100生物分析儀(Agilent,USA)和Illumina (Kapa Biosciences, Woburn, MA,USA)文庫定量試劑盒進(jìn)行評(píng)估。高通量測(cè)序委托杭州聯(lián)川生物技術(shù)有限公司完成,所用平臺(tái)為高通量二代測(cè)序儀 (Illumina,HiSeq 4000,USA)。測(cè)序原始數(shù)據(jù)已提交至NCBI數(shù)據(jù)庫(項(xiàng)目號(hào)為PRJNA855296,樣品編號(hào)為SAMN29487930-SAMN29487947)。

    1.4 數(shù)據(jù)處理及分析

    高通量測(cè)序獲得的下機(jī)數(shù)據(jù)根據(jù)Barcode拆分后,去除Barcode和接頭序列,并進(jìn)行質(zhì)控、過濾,得到Clean data。使用PEAR(v0.9.6)軟件對(duì)reads進(jìn)行拼接合并獲得原始Tags。用Fqtrim篩選Q20>20的序列。利用Vsearch(v2.3.4)去除嵌合體獲得最終有效數(shù)據(jù)用于進(jìn)一步分析。樣本的α多樣性指數(shù)采用mothur(http://www.mothur.org/wiki/ Schloss_SOP#Alpha_diversity)軟件計(jì)算。應(yīng)用Usearch 軟件(vesion 7.0 http://drive5.com/uparse/)對(duì)有效數(shù)據(jù)序列在97%的相似度下進(jìn)行OTU (Operational Taxonomic Units)聚類,采用Qiime 2進(jìn)行OTU物種注釋分析,獲得樣品OTU的分類學(xué)信息。韋恩圖、PCoA分析、RDA分析等利用聯(lián)川云平臺(tái)(https://www.omicstudio.cn/index)進(jìn)行。利用FUNGuild軟件對(duì)真菌進(jìn)行功能預(yù)測(cè)分析(Nguyen et al., 2016)。采用SPSS 26對(duì)土壤理化性質(zhì)進(jìn)行顯著性差異分析。

    2 結(jié)果與分析

    2.1 3個(gè)樣地黑老虎根際土壤理化特性

    由表1可知,3個(gè)不同生境下,黑老虎根際土壤的理化指標(biāo)存在顯著差異,野生生境土壤主要養(yǎng)分指標(biāo)高于栽培生境。3個(gè)生境中,貴陽栽培生境(GY,pH=4.66)和榕江野生生境(RJ,pH=4.63)下的根際土壤為酸性(pH<6.5),福泉野生生境下的根際土壤pH呈中性(FQ,pH=7.22),主要是由于福泉土壤類型以石灰土最多,土壤pH呈中性,福泉素有“亞洲磷都”之美譽(yù),土壤富含P,此外K含量也較多。RJ根際土壤的TN、AN和SOM含量在3個(gè)樣地中最高;FQ根際土壤的TP、TK、AP含量在3個(gè)樣地中最高;GY栽培生境下根際土壤TN、TP、TK、AN、AK、AP及SOM含量均顯著低于野生生境(RJ和FQ)。野生生境下的土壤營(yíng)養(yǎng)元素和有機(jī)質(zhì)含量高于栽培生境,并且生態(tài)系統(tǒng)更為平衡。原因可能為野生生境為茂密的樹林,存在較多腐殖質(zhì),有機(jī)質(zhì)含量更為豐富且含水量較高,土壤較為疏松,礦質(zhì)營(yíng)養(yǎng)元素不易流失;而人工苗圃在管理過程中施用肥料不多,腐殖質(zhì)和有機(jī)質(zhì)偏少且土壤存在一定程度的板結(jié),土壤中的N、P、K等元素隨雨水流失較為嚴(yán)重。

    2.2 黑老虎根際土壤真菌及根部?jī)?nèi)生真菌測(cè)序

    對(duì)Illumina高通量測(cè)序數(shù)據(jù)進(jìn)行雙端拼接、質(zhì)量控制、嵌合體過濾后,進(jìn)行高質(zhì)量數(shù)據(jù)統(tǒng)計(jì),共獲得有效序列1 463 619個(gè),樣品有效reads為76 878~84 980對(duì),每個(gè)樣品平均為81 312對(duì)reads,有效率在93.59%~98.70%之間,Q20全部在99%以上。18個(gè)樣品的OTU稀釋曲線隨測(cè)序數(shù)目增加逐漸趨于平緩,并且各樣品測(cè)序覆蓋度均在99.95%以上(圖1),表明測(cè)序數(shù)據(jù)量能較全面地反映測(cè)序樣品的真菌群落組成,測(cè)序質(zhì)量評(píng)價(jià)各項(xiàng)指標(biāo)均符合ITS測(cè)序要求,可以進(jìn)行后續(xù)分析和研究。

    以≥97%的相似性聚類得到4 328個(gè)OTU,3個(gè)樣地根際土壤真菌和根部?jī)?nèi)生真菌分別檢測(cè)到OTU數(shù)為3 867個(gè)和801個(gè),其中根際土壤中各地區(qū)OTU 數(shù)從多到少依次為RJ_Rh(1 358個(gè))> FQ_Rh(1 090個(gè))> GY_Rh(994個(gè)),共有的OTU為72個(gè)。共有的OTU共注釋到5個(gè)門,絕大多數(shù)為子囊菌門(50個(gè),占比68.4%);共注釋到49個(gè)屬,占比前幾位的有未分類屬(7個(gè),9.7%)、被孢霉屬(6個(gè),8.3%), 外瓶柄霉屬、 柱孢屬、叢赤殼屬和隱球菌屬均為3個(gè),各占比4.2%。根部?jī)?nèi)生真菌中,各地區(qū)OTU數(shù)從多到少依次為FQ_R(330)> GY_R(243)>RJ_R(153),共有的OTU為14個(gè)。14個(gè)共有的OTU共注釋到2個(gè)門,子囊菌門有13個(gè),占比92.9%;共注釋到11個(gè)屬,占比最高的為外瓶柄霉屬(3個(gè),21.4%),其余屬均為1個(gè),各占比7.1%。根際土壤真菌和內(nèi)生真菌所有樣品共有的OTU僅為6個(gè),注釋到2個(gè)門,子囊菌門為5個(gè),占比83.3%;注釋到5個(gè)屬,占比最高的為外瓶霉屬(2個(gè),33.3%),其余分別為被孢霉屬、 柱孢屬和叢赤殼屬(圖2)。圖2結(jié)果顯示3個(gè)樣地樣品中,真菌OTU數(shù)量均為根際土壤>根內(nèi),F(xiàn)Q野生樣地OTU多樣性較高。共有的OTU中,被孢霉屬、外瓶柄霉屬、 柱孢屬和叢赤殼屬最為常見。

    2.3 黑老虎根際土壤真菌和根部?jī)?nèi)生真菌多樣性

    Alpha多樣性分析表明,3種生境下,根際土壤的真菌群落多樣性和豐富度均顯著高于根部?jī)?nèi)生真菌。3個(gè)樣地根際土壤真菌中,香農(nóng)指數(shù)(Shannon index)和辛普森指數(shù)(Simpson index)無顯著差異,RJ_Rh中Chao1指數(shù)和Ace指數(shù)顯著高于GY_Rh和FQ_Rh。3個(gè)樣地根部?jī)?nèi)生真菌中,F(xiàn)Q_R香農(nóng)指數(shù)顯著高于GY_R和RJ_R,GY_R辛普森指數(shù)顯著低于RJ_R和FQ_R,而RJ_R的Chao1指數(shù)和Ace指數(shù)顯著低于FQ_R,與GY_R無顯著差異(表2)。這說明了3個(gè)生境根際土壤真菌多樣性和優(yōu)勢(shì)度無顯著差異,RJ_R的根際真菌群落豐度高于其余兩地,3個(gè)樣地根部?jī)?nèi)生真菌中, FQ_R的真菌群落多樣性顯著高于其他兩地,GY_R栽培生境的真菌優(yōu)勢(shì)度顯著低于兩個(gè)野生生境(RJ_R和FQ_R)。

    通過無權(quán)重距離矩陣(unweighted_unifrac_distance_matrix)算法對(duì)3個(gè)樣地樣品進(jìn)行主坐標(biāo)分析(PCoA)分析,結(jié)果顯示,主坐標(biāo)1和主坐標(biāo)2分別解釋了25.25%和13.76%的多樣性,合計(jì)解釋了39.01%的多樣性,3個(gè)生境下的根部?jī)?nèi)生真菌群落聚集于第二、第三象限,根際土壤真菌群落聚集于第一、第四象限,顯示3個(gè)不同生境中根際土壤真菌更為相似,而根部?jī)?nèi)生真菌群落更為相似(圖3)。結(jié)果表明,黑老虎根部?jī)?nèi)生真菌受地域影響較小,具有明顯的種屬特異性。另外,從PCoA結(jié)果可以看出,同一類樣本的3個(gè)生物重復(fù)均能很好地聚在一起,說明本研究所采用的樣本重復(fù)性好,數(shù)據(jù)可靠性高。

    2.4 黑老虎根際土壤和根部?jī)?nèi)生真菌組成和結(jié)構(gòu)

    3個(gè)生境下共18個(gè)黑老虎根際土壤和根樣品中,共檢測(cè)到了真菌12個(gè)門,依次為GY(10個(gè))、RJ(9個(gè))和FQ(9個(gè)),3個(gè)樣地相差不大。3個(gè)樣地6組樣品中,真菌優(yōu)勢(shì)類群均為子囊菌門(Ascomycota)和擔(dān)子菌門(Basidiomycota),其中,子囊菌門占比最高, 占比范圍48.02%~88.38%,在栽培生境(GY_R)的根內(nèi)(Y1)占比高達(dá)88.38%,而在栽培生境的根際土壤(GY_Rh)中占比最低(48.02%)。GY_Rh和FQ_Rh中,存在26.17%和13.68%的未分類真菌(圖4:A)。

    在屬水平上,共注釋到509個(gè)屬,3個(gè)生境中注釋到的真菌屬依次為RJ(375)> FQ(299)>GY(274)。3個(gè)生境下,優(yōu)勢(shì)屬的結(jié)構(gòu)和豐度存在明顯差異,根部?jī)?nèi)生真菌中,栽培生境(GY_R)下,外瓶霉屬(Exophiala,77.81%)優(yōu)勢(shì)明顯,其他屬占比均不超過4%,2個(gè)野生生境下,RJ_R的優(yōu)勢(shì)菌群為Cladophialophora(21.68%)、黃絲菌屬(Cephalotheca,19.46%)和Flagelloscypha(19.32%),F(xiàn)Q_R的優(yōu)勢(shì)菌群為蠟殼耳目未分類屬(19.7%)、外瓶霉屬(15.20%)和球囊菌門未分類屬(Glomeromycota_unclassified,12.12%)。3個(gè)樣地的根際土壤真菌,被孢霉屬(Mortierella)占比均最高,分別為8.57%(GY_Rh)、21.43%(RJ_Rh)和13.53%(FQ_Rh)。此外, GY_Rh優(yōu)勢(shì)菌群還有癬囊腔菌屬(Plectosphaerella,4.34%)和外瓶霉 屬(4.15%), RJ_Rh中優(yōu)勢(shì)菌群為新叢赤殼屬(Neonectria,7.15%),Y3中優(yōu)勢(shì)菌群為癬囊腔菌屬(5.97%)和柱孢屬(Cylindrocarpon,4.31%)(圖4:B)。結(jié)果顯示,栽培生境下,根部?jī)?nèi)生真菌菌群具有一定的偏好性,部分菌群具有優(yōu)越性,而野生生境下,菌群均勻度更高。

    以相對(duì)豐度Top20的屬制作的熱圖顯示,3個(gè)生境下,根部?jī)?nèi)生真菌和根際真菌群落均存在顯著差異,其中根際土壤真菌群落聚為一類,與PCoA聚類圖結(jié)果一致,表明黑老虎根際土壤真菌群落結(jié)構(gòu)受生境影響較小,黑老虎根部?jī)?nèi)生真菌群落具有一定的物種特異性(圖 5)。

    2.5 黑老虎根際土壤真菌及根部?jī)?nèi)生真菌群落功能預(yù)測(cè)

    利用FUNGuild軟件對(duì)黑老虎根際土壤及根部?jī)?nèi)生真菌進(jìn)行功能預(yù)測(cè),結(jié)果顯示真菌7個(gè)營(yíng)養(yǎng)類型在3個(gè)生境的根際土壤和根組織中分布不均,其中病理-腐生營(yíng)養(yǎng)型和腐生營(yíng)養(yǎng)型兩個(gè)類型為優(yōu)勢(shì)菌群。3個(gè)樣地的根際土壤樣本中,腐生-共生營(yíng)養(yǎng)型占比相當(dāng)高, 而在根部?jī)?nèi)生真菌中占比較低。與內(nèi)生真菌相比,根際真菌中不同營(yíng)養(yǎng)類型真菌菌群的分布更為均勻。根部?jī)?nèi)生真菌中,栽培生境下(GY_R),病理-腐生營(yíng)養(yǎng)型占比較高,RJ_R的腐生營(yíng)養(yǎng)型占比較高,F(xiàn)Q_R的共生營(yíng)養(yǎng)型占比較高(圖 6:A)。

    3個(gè)生境下6組樣本預(yù)測(cè)得到61個(gè)生態(tài)共位群(guild),不同生態(tài)共位群在3個(gè)生境下的根際土壤真菌分布更為相似,GY_R群落中動(dòng)物病原-真菌寄生蟲-未定義腐生營(yíng)養(yǎng)型占比較高,而RJ_R中未定義腐生營(yíng)養(yǎng)型占比較高,F(xiàn)Q_R中動(dòng)物病原-真菌寄生蟲-未定義腐生營(yíng)養(yǎng)型和叢枝菌根占比較高(圖6:B)。

    2.6 黑老虎根際土壤真菌群落結(jié)構(gòu)與土壤環(huán)境因子相關(guān)性

    利用土壤理化因子與根部?jī)?nèi)生真菌和根際真菌群落多樣性指數(shù)分別做Pearson相關(guān)性分析。

    根部?jī)?nèi)生真菌中,TK與根部?jī)?nèi)生真菌香農(nóng)指數(shù)和辛普森指數(shù)正相關(guān), TP與香農(nóng)指數(shù)正相關(guān)(P<0.05),但未達(dá)極顯著水平(圖7:A)。根際土壤真菌中,土壤SOM、AN和TN與根際土壤真菌Ace指數(shù)和Chao1指數(shù)呈極顯著(P<0.01)正相關(guān)(圖7:B)。這說明土壤環(huán)境因子對(duì)黑老虎內(nèi)生真菌和根際土壤真菌群落結(jié)構(gòu)的影響不同,TK和TP主要影響內(nèi)生真菌群落的多樣性,而SOM、TN和AN主要影響根際土壤真菌群落的豐度,對(duì)根際土壤真菌群落結(jié)構(gòu)的影響大于內(nèi)生真菌。

    為探究土壤理化因子與根際真菌群落結(jié)構(gòu)的相關(guān)性,通過冗余分析(redundancy analysis, RDA)分別研究了土壤理化因子與根際真菌與門水平和豐度排序前20的屬之間的相互關(guān)系(圖 8)。在門水平上,第一、第二排序軸累計(jì)解釋率分別為62.5%和22.5%, 累計(jì)解釋率為85.0%, 能較好地反映黑老虎土壤理化因子與根際真菌門水平的相互關(guān)系。SOM、AN、TN和TK與擔(dān)子菌門(Basidiomycota)和接合菌門(Zygomycota)呈現(xiàn)較強(qiáng)的正相關(guān)。在屬水平上,第一、第二排序軸累計(jì)解釋率分別為75.6%和15.0%,累計(jì)解釋率為90.6%,能較好地反映黑老虎根際真菌群落與土壤理化因子之間的相互關(guān)系。TN、AN、SOM和TK與新叢枝殼屬、黃絲菌屬、被孢霉屬等箭頭方向一致,呈現(xiàn)較強(qiáng)的正相關(guān),與真菌_未分類(Fungi_unclassified)、外瓶霉屬、肉座菌屬(Hypocrea)等箭頭方向相反,呈現(xiàn)明顯的負(fù)相關(guān)。綜上表明,SOM、TN、AN和TK對(duì)黑老虎根際土壤真菌群落豐度影響較大,是影響黑老虎根際土壤真菌結(jié)構(gòu)的主要因素。

    3 討論與結(jié)論

    本研究利用高通量測(cè)序技術(shù)研究了不同生境下黑老虎根內(nèi)生真菌及根際土壤真菌群落組成和多樣性,3個(gè)樣地根際土壤和根部?jī)?nèi)生真菌共有的OTU分別為72個(gè)和14個(gè),根際土壤真菌和內(nèi)生真菌所有樣品共有的OTU僅為6個(gè),各地特有的OTU有較大差異。這表明根際土壤真菌具有一定的地域?qū);院退拗鲗R恍?,與其他物種結(jié)果相類似(周婕等,2019)。3個(gè)樣地中,黑老虎根際土壤真菌的OTU數(shù)量和多樣性均顯著高于根部?jī)?nèi)生真菌,土壤中含有微生物所需的食物來源和生長(zhǎng)環(huán)境,是地球上微生物含量最為豐富的區(qū)域(Roesch et al.,2007),與重樓(王艷等,2019)、黃精(樊銳鋒等,2021)、川穹(張秫華等,2021)等研究結(jié)果一致。此外,3個(gè)樣地根際土壤真菌多樣性(香農(nóng)指數(shù)和辛普森指數(shù))無顯著差異,原因可能為3個(gè)樣地均在貴州省內(nèi),地理位置、大氣候和生態(tài)環(huán)境具有一定的相似性。但是,3個(gè)樣地根際土壤真菌的豐富度(Chao1指數(shù)和Ace指數(shù))存在顯著差異,其中RJ_Rh根際土壤真菌豐富度顯著高于GY_Rh和FQ_Rh,de Boe等(2005)研究表明植物根系及其分泌物、動(dòng)植物殘?bào)w和SOM是土壤真菌可利用的主要能源,榕江土壤豐富的有機(jī)質(zhì)為真菌群落的繁殖提供了豐富的養(yǎng)分,可能是榕江根際土壤真菌豐富度高于另外兩地的原因。3個(gè)樣地根部?jī)?nèi)生真菌中, FQ_R的真菌群落多樣性顯著高于其他兩地,GY_R栽培生境的真菌優(yōu)勢(shì)度顯著低于兩個(gè)野生生境(RJ_R和FQ_R),可能是因?yàn)橐吧呈苋斯じ深A(yù),栽培生境的植被、土壤結(jié)構(gòu)和營(yíng)養(yǎng)元素趨于不均衡,導(dǎo)致栽培生境下黑老虎根內(nèi)生真菌的多樣性偏低。

    本研究發(fā)現(xiàn)子囊菌門和擔(dān)子菌門是黑老虎根際和根部?jī)?nèi)生真菌的優(yōu)勢(shì)類群,子囊菌門大都為陸生,可分解土壤中動(dòng)植物腐爛的遺骸和有機(jī)質(zhì),釋放有機(jī)物供植物吸收利用,其含量的高低與土壤肥力有關(guān)。擔(dān)子菌門常與植物形成共生菌根,促進(jìn)植物的生長(zhǎng)和發(fā)育(Miguel et al., 2020)。本研究結(jié)果與五味子屬(江林春等,2020)及其他中藥材如重樓(王艷等,2019)、黃精(樊銳鋒等,2021)、川穹(張秫華等,2021)等一致。在屬水平上,根部?jī)?nèi)生真菌中, GY_R、RJ_R和FQ_R的優(yōu)勢(shì)菌屬分別為外瓶霉屬(77.81%)、Cladophialophora(21.68%)和蠟殼耳目未分類屬(19.7%)、外瓶霉屬(15.20%)。外瓶霉屬在栽培生境(GY)的根中豐度較高(77.81%)。其原因可能為外瓶霉屬有在土壤、水和空氣等不同生態(tài)系統(tǒng)中有著廣泛的分布,多為人、畜和水產(chǎn)魚類等的真菌病害,在植物中,有耐重金屬特性的報(bào)道(王璐,2011)。栽培生境下,長(zhǎng)期施用化肥農(nóng)藥,受到人類活動(dòng)干擾較多,外瓶霉屬上升為優(yōu)勢(shì)菌群。根際土壤真菌中,3個(gè)樣地的優(yōu)勢(shì)菌屬均含被孢霉屬,而其他優(yōu)勢(shì)菌群各不相同。被孢霉屬為土壤碳及養(yǎng)分轉(zhuǎn)化的關(guān)鍵成員,在植物生長(zhǎng)發(fā)育、代謝中均發(fā)揮關(guān)鍵作用(寧琪等,2022)。在屬水平上,對(duì)不同生境下根際土壤和根部?jī)?nèi)生真菌的α多樣性分析顯示,不同生境下黑老虎根際土壤和根部?jī)?nèi)生真菌群落的豐度和多樣性存在顯著差異。FUNGuild功能預(yù)測(cè)結(jié)果顯示,7個(gè)營(yíng)養(yǎng)類型在6組樣本中分布不均,栽培生境下(GY)病理-腐生營(yíng)養(yǎng)型及動(dòng)物病原-真菌寄生蟲-未定義腐生營(yíng)養(yǎng)型真菌在黑老虎根內(nèi)生菌中比例最高,推測(cè)可能是由于在栽培生境下,由于生態(tài)平衡遭到破壞,造成植株體內(nèi)病原菌和寄生蟲等侵害機(jī)會(huì)增加,病原菌潛伏或以條件性致病方式寄生于植物根內(nèi),是導(dǎo)致栽培生境下植物更容易發(fā)生病害的原因之一(Rodriguez et al., 2009)。王艷等(2019)發(fā)現(xiàn),野生生境下的重樓內(nèi)生真菌和土壤真菌菌群豐度和多樣性顯著高于園區(qū)環(huán)境,與本研究結(jié)果較為一致。

    植物根際土壤和根內(nèi)真菌在植物生長(zhǎng)發(fā)育、低于逆境和維持土壤微生態(tài)環(huán)境中發(fā)揮重要作用,其結(jié)構(gòu)和組成受土壤環(huán)境因子如濕度、pH值、N、P、K、SOM等理化因素的影響(Roesch et al., 2007)。本研究結(jié)果顯示,3個(gè)不同生境下,黑老虎根際土壤的理化指標(biāo)存在明顯差異,野生生境下的土壤營(yíng)養(yǎng)元素和有機(jī)質(zhì)含量高于栽培生境,說明栽培生境下, 要注意合理施用營(yíng)養(yǎng)元素和有機(jī)肥,確保植物生長(zhǎng)的需要和土壤生態(tài)環(huán)境平衡。同時(shí),本研究結(jié)果也為后續(xù)的人工苗圃管理提供了重要信息。土壤理化因子與根部?jī)?nèi)生真菌多樣性相關(guān)分析顯示,土壤TK和TP與根部?jī)?nèi)生真菌香農(nóng)指數(shù)相關(guān)顯著,表明TK和TP主要影響根部?jī)?nèi)生真菌的多樣性,對(duì)豐度和優(yōu)勢(shì)度影響不明顯。根際土壤真菌中,土壤SOM、AN和TN與根際土壤真菌Ace和Chao1指數(shù)極顯著相關(guān)。這表明土壤環(huán)境因子對(duì)根際土壤真菌群落結(jié)構(gòu)影響大于根部?jī)?nèi)生真菌且主要影響根部?jī)?nèi)生真菌的多樣性,而對(duì)根際土壤真菌的影響主要為豐度。李銘杰等(2021)發(fā)現(xiàn),土壤SOM、pH、AP、AK與北蒼術(shù)根區(qū)土壤叢枝菌根真菌Shannon指數(shù)、Simpson指數(shù)呈正相關(guān),其中SOM對(duì)叢枝菌根真菌群落多樣性影響最大。楊蕊毓等(2022)研究發(fā)現(xiàn),川麥冬根圍土壤叢枝菌根多樣性指數(shù)及群落組成與土壤理化因子TK、TP、TN相關(guān)性顯著。金輝等(2021)對(duì)紫花針茅根際土壤和內(nèi)生真菌研究分析發(fā)現(xiàn),土壤理化因子pH、AK、Fe、Ca 和Mg 與Ace指數(shù)、Chao1指數(shù)和Shannon指數(shù)等真菌多樣性指數(shù)顯著相關(guān)(P<0.05)或極顯著相關(guān)(P<0.01),與本研究結(jié)果具有部分相似性。N是植物生長(zhǎng)發(fā)育所需的大量元素,植物通過根系分泌物和物生物互作,形成植物-土壤的氮素化學(xué)循環(huán),因此N與土壤微生物群落和根內(nèi)生菌均發(fā)生交互作用(Hobbie et al., 2007)。有機(jī)質(zhì)含量影響土壤真菌的群落結(jié)構(gòu),主要由于SOM在微生物和酶的作用下發(fā)生分解,同時(shí)為土壤微生物提供豐富的養(yǎng)分來促進(jìn)了土壤真菌的生長(zhǎng) (隋躍宇等, 2009; Smith et al., 2014)。

    土壤TN、AN、SOM和TK是影響黑老虎根際土壤真菌結(jié)構(gòu)的主要因素。RDA分析結(jié)果顯示,被孢霉屬與TN、AN和SOM呈現(xiàn)較強(qiáng)的正相關(guān)。被孢霉屬是接合菌亞門接合菌綱毛霉目被孢霉科,能顯著提高土壤可溶性有機(jī)碳、AN和AP的含量(寧琪等,2022)。本研究結(jié)果表明,多施有機(jī)肥、氮肥和鉀肥,可以改善土壤微生物菌群結(jié)構(gòu),促進(jìn)黑老虎生長(zhǎng)發(fā)育。本研究為今后黑老虎根內(nèi)有益真菌的分離與鑒定提供了重要信息,同時(shí)也為黑老虎的引種馴化和生態(tài)栽培提供了理論依據(jù)。

    參考文獻(xiàn):

    BERHONGARAY G, JANSSENS IA, KING JS, et al., 2013. Fine root biomass and turnover of two fast-growing poplar genotypes in a short-rotation coppice culture [J]. Plant Soil, 373(1-2):269-283.

    BONFANTE P, ANCA IA, 2009. Plants, mycorrhizal fungi, and bacteria: a network of interactions [J]. Ann Rev Microbiol, 63: 363-383.

    CHAPARRO JM, BADRI DV, VIVANCO JM, 2013. Rhizosphere microbiome assemblage is affected by plant development [J]. ISME J, 8(4):? 790-803.

    DE BOER W, FOLMAN LB, SUMMERBELL RC, et al., 2005. Living in a fungal world:? impact of fungi on soil bacterial niche development [J]. FEMS Microbiol Rev, 29(4):? 795-811.

    EDWARDS J, JOHNSON C, SANTOS-MEDELLíN C, et al., 2015. Structure, variation, and assembly of the root-associated microbiomes of rice [J]. Proc Natl Acad Sci, 112(8): E911-920.

    FAN RF, WANG RF, DU YQ, et al., 2021. Community composition and ecological function of rhizosphere fungi and medicinal parts endophytic fungi in Polygonatum sibiricum [J]. Guihaia, 41(5): 799-807.? [樊銳鋒, 王若凡, 杜艷秋, 等, 2021. 黃精根際及藥用部位內(nèi)生真菌群落組成和生態(tài)功能分析 [J]. 廣西植物, 41(5): 799-807.]

    FIERERN, 2017. Embracing the unknown:? disentangling the complexities of the soil microbiome [J]. Nat Rev Microbiol, 15(10): 579-590.

    GAO JF, ZHOU W, LIU N, 2022. Differential analysis of phenolic acids from different parts of Kadsura coccinea? [J]. J Chin Med Mat, 45(2): 351-358.? [高漸飛, 周瑋, 劉妮, 2022. 黑老虎不同部位酚酸類成分差異分析 [J]. 中藥材, 45(2): 351-358.]

    GAO Y, ZHOU J, RUAN HL, 2020. Trichothecenes from an endophytic fungus Alternaria sp. sb23 [J]. Planta Med, 86(13/14): 976-982.

    HOBBIE SE, OGDAHL M, CHOROVER J, et al., 2007. Tree species effects on soil organic matter dynamics:? the role of soil cation composition? [J]. Ecosystems, 10(6):? 999-1018.

    HUANG SZ, DUAN LP, WANG H, et al., 2019. Two new AChE inhibitors isolated from Li folk herb Heilaohu “Kadsura coccinea” stems? [J]. Molecules, 24(19): 3628.

    JIA M, CHEN L, XIN HL, et al., 2016. A friendly relationship between endophytic fungi and medicinal plants:? a systematic review [J]. Front Microbiol, 7: 906.

    JIANG LC, WEI JM, BAI SS, et al., 2020. Characteristics of rhizosphere soil fungi community of Schisandra sphenanthera Rehd. Et Wils [J]. J Yanan Univ (Nat Sci Ed), 39(4): 15-20.? [江林春, 韋潔敏, 白尚尚, 等, 2020. 華中五味子根際土壤真菌群落特征研究 [J]. 延安大學(xué)學(xué)報(bào)(自然科學(xué)版), 39(4): 15-20.]

    JIN H, YANG XY, LIU HY, et al., 2021. The structure of rhizosphere soil and endophytic fungal communities associated with Stipa purpurea and their correlation with soil environmental factors [J]. Acta Microbiol Sin, 61(11): 3520-3541.? [金輝, 楊曉燕, 柳皓月, 等, 2021. 紫花針茅根際和體內(nèi)真菌群落結(jié)構(gòu)及與土壤環(huán)境因子的相關(guān)性研究 [J]. 微生物學(xué)報(bào),? 61(11): 3520-3541.]

    LI MJ, ZHOU ZJ, XING LJ, et al., 2021. Correlation analysis between diversity of AM fungi and nutrient in root-zone soil of Atractylodes chinensis [J]. Chin J Biol Control, 37(6): 1288-1297.? [李銘杰, 周志杰, 邢禮軍, 等, 2021. 北蒼術(shù)根區(qū)土壤中AM真菌多樣性及其與土壤養(yǎng)分相關(guān)性分析 [J]. 中國(guó)生物防治學(xué)報(bào), 37(6): 1288-1297.]

    LIU T, MA N, JIN JF, et al., 2023. Composition structure and ecological function analyses of endophytic and rhizosphere soil fungi in Kadsura coccinea [J]. Guihaia, 43(5): 869-879.? [劉濤, 馬楠, 金吉芬, 等, 2023. 黑老虎內(nèi)生真菌及根際土壤真菌的群落結(jié)構(gòu)與生態(tài)功能分析? [J]. 廣西植物, 43(5): 869-879.]

    LIU T, WEI Q, WU XB, et al., 2009. Research status and development and utilization prospect of wild fruit species Kadsura coccinea (Lem.) A. C. Smith [J]. J Anhui Agric Sci, 37(23): 10971-10972.? [劉濤, 韋茜, 吳小波, 等, 2009. 野生水果黑老虎的研究現(xiàn)狀及開發(fā)利用前景 [J]. 安徽農(nóng)業(yè)科學(xué), 37(23): 10971-10972.]

    LU HY, WEI TY, LOU HH, et al., 2021. A critical review on communication mechanism within plant-endophytic fungi interactions to cope with biotic and abiotic stresses [J]. J Fung (Basel), 7(9): 719.

    MIGUEL RM, JOSE GMV, MEIKE P, et al., 2020. Diversity of fungi in soils with different degrees of degradation in Germany and Panama? [J]. Mycobiology, 48(1): 20-28.

    NGUYEN NH, SONG Z, BATES ST, et al., 2016. FUNGuild:? An open annotation tool for parsing fungal community datasets by ecological guild [J]. Fung Ecol, 20:? 241-248.

    NING Q, CHEN L, LI F, et al., 2022. Effects of Mortierella on nutrient availability and straw decomposition in soil [J]. Acta Pedol Sin, 59(1): 206-217.? [寧琪, 陳林, 李芳, 等, 2022. 被孢霉對(duì)土壤養(yǎng)分有效性和秸稈降解的影響 [J]. 土壤學(xué)報(bào), 59(1):? 206-217.]

    PAN ZY, FU JF, ZHOU RJ, et al., 2007. Preliminary study of soil fungi in Schisandra chinensis garden and screening antagonism strains [J]. J Jilin Agric Univ, 29(6): 636-639.? [潘爭(zhēng)艷, 傅俊范, 周如軍, 等, 2007. 五味子園根際真菌多樣性初探及拮抗菌株篩選 [J]. 吉林農(nóng)業(yè)大學(xué)學(xué)報(bào), 29(6): 636-639.]

    PATTENCL, GLICK BR, 2002. Role of Pseudomonas putida indoleacetic acid in development of the host plant root system [J]. Appl Environ Microb, 68(8): 3795-3795.

    QIN D, SHEN WY, WANG JQ, et al., 2019. Enhanced production of unusual triterpenoids from Kadsura angustifolia fermented by a symbiont endophytic fungus, Penicillium sp. SWUKD4.1850 [J]. Phytochemistry, 158: 56-66.

    QIN D, SHEN WY, GAO TC, et al., 2020. Kadanguslactones A-E, further oxygenated terpenoids from Kadsura angustifolia fermented by a symbiotic endophytic fungus, Penicillium ochrochloron SWUKD4.1850 [J]. Phytochemistry, 174: 112335.

    RODRIGUEZ RJ, WHITE JF, ARNOLD AE, et al., 2009. Fungal endophytes:? diversity and functional roles [J]. New Phytol, 182(2): 314-330.

    ROESCH LF, FULTHORPE RR, RIVA A, et al., 2007. Triplett EW. Pyrosequencing enumerates and contrasts soil microbial diversity [J]. ISME J, 1(4):? 283-290.

    SMITH AP, MARíN-SPIOTTA E, DE GRAAFF MA, et al., 2014. Microbial community structure varies across soil organic matter aggregate pools during tropical land cover change [J]. Soil Biol Biochem, 77:? 292-303.

    SRITALAHAREUTHAI V, TEMVIRIYANUKUL P, ON-NOM N, et al., 2020. Phenolic profiles, antioxidant, and inhibitory activities of Kadsura heteroclita (Roxb.) Craib and Kadsura coccinea (Lem.) A.C. Sm [J]. Foods, 9(9): 1222.

    SUI YY, JIAO XG, GAO CS, et al., 2009. The relationship among organic matter content and soil microbial biomass and soil enzyme activities [J]. Chin J Soil Sci, 40(5):? 1036-1039. [隋躍宇, 焦曉光, 高崇生, 等, 2009. 土壤有機(jī)質(zhì)含量與土壤微生物量及土壤酶活性關(guān)系的研究 [J]. 土壤通報(bào), 40(5):? 1036-1039.]

    TIWARI P, BAE H, 2022. Endophytic fungi:? key insights, emerging prospects, and challenges in natural product drug discovery [J]. Microorganisms, 10(2): 360.

    WANG L, 2011. Resistance screening on heavy metal and taxonomic studies of strains belong to Phialophora and Exophiala? [D]. Kumming: Yunnan University.? [王璐,? 2011. 瓶霉屬(Phialophora)和外瓶霉屬(Exophiala)菌株重金屬耐性篩選及分類學(xué)研究 [D]. 昆明: 云南大學(xué).]

    WANG L, QIN D, ZANG K, et al., 2017. Metabolites from the co-culture of nigranoic acid and Umbelopsis dimorpha SWUKD3.1410, an endophytic fungus from Kadsura angustifolia [J]. Nat Prod Res, 31(12): 1414-1421.

    WANG LJ,LIAO SQ, LIANG J, et al., 2021. Analysis and evaluation of nutritional components of Kadsura coccinea seeds [J]. Chin Oils Fats, 46(12): 112-117.? [王麗軍, 廖蘇奇, 梁潔, 等, 2021. 黑老虎種子的營(yíng)養(yǎng)成分分析及評(píng)價(jià) [J]. 中國(guó)油脂, 46(12): 112-117.]

    WANG Y, CHANG F, CHENG HY, et al., 2019. Fungal diversity and differential analysis of fungi community composition in Paris polyphylla var. chinensis rhizosphere and endophytic [J]. Chin Trad Herb Drugs, 50(5): 1232-1237.? [王艷, 常帆, 程虎印, 等, 2019. 重樓根際及藥用部位內(nèi)生真菌多樣性與群落結(jié)構(gòu)差異分析 [J]. 中草藥, 50(5): 1232-1237.]

    YANG RY, DENG QZ, TIAN LP, et al., 2022. Diversity of arbuscular mycorrhizal fungi in the rhizosphere of Ophiopogon japonicus in different habitats [J]. Acta Bot Boreal-Occident Sin, 42(1): 145-153.? [楊蕊毓, 鄧錡璋, 田麗平, 等, 2022. 不同生境下川麥冬根圍土壤叢枝菌根真菌多樣性 [J]. 西北植物學(xué)報(bào),? 42(1): 145-153.]

    YU J, ZHOU XF, YANG SJ, et al., 2013. Design and application of specific 16S rDNA-targeted primers for assessing endophytic diversity in Dendrobium officinale using nested PCR-DGGE [J]. Appl Microbiol Biotechnol, 97(22): 9825-9836.

    ZHANG SH, FANG Q, JIA HM, et al., 2021. Difference analysis of fungal community among bulk soil, rhizosphere, and rhizomes of Ligusticum chuanxiong Hort [J]. Biotechnol Bull, 37(4): 56-69.? [張秫華, 方千, 賈紅梅, 等, 2021. 川芎非根際、根際及根莖內(nèi)生真菌群落差異分析 [J]. 生物技術(shù)通報(bào), 37(4): 56-69.]

    ZHANG YQ, XU YH, WEN XS, et al., 2021. Response surface methodology for optimizing fermentation process of the endophytic fungus Purpureocillium lilacinum WG9 preventing Schisandra chinensis brown spot disease [J]. Mycosystema, 40(6): 1561-1574.? [張譽(yù)薺, 許永華, 文湘穗, 等, 2021. 響應(yīng)面法優(yōu)化北五味子褐斑病內(nèi)生生防真菌淡紫擬青霉WG9發(fā)酵工藝及發(fā)酵產(chǎn)物穩(wěn)定性研究 [J]. 菌物學(xué)報(bào), 40(6): 1561-1574.]

    ZHOU J, MIAO YF, FANG K, et al., 2019. Diversity of the endophytic and rhizosphere soil fungi of Ageratina adenophora [J]. Ecol Sci, 38(5):? 1-7.? [周婕, 苗一方, 方楷, 等, 2019. 紫莖澤蘭內(nèi)生真菌及其根際土壤真菌的多樣性研究 [J]. 生態(tài)科學(xué), 38(5): 1-7.]

    (責(zé)任編輯 李 莉 王登惠)

    猜你喜歡
    根際內(nèi)生生境
    根際微生物對(duì)植物與土壤交互調(diào)控的研究進(jìn)展
    棗樹適應(yīng)干旱生境研究進(jìn)展
    落葉果樹(2021年6期)2021-02-12 01:29:10
    植物內(nèi)生菌在植物病害中的生物防治
    內(nèi)生微生物和其在作物管理中的潛在應(yīng)用
    “黨建+”激活鄉(xiāng)村發(fā)展內(nèi)生動(dòng)力
    授人以漁 激活脫貧內(nèi)生動(dòng)力
    黃花蒿葉水提物對(duì)三七根際尖孢鐮刀菌生長(zhǎng)的抑制作用
    促植物生長(zhǎng)根際細(xì)菌HG28-5對(duì)黃瓜苗期生長(zhǎng)及根際土壤微生態(tài)的影響
    金秀瑤族長(zhǎng)鼓舞傳承的文化生境
    體育科技(2016年2期)2016-02-28 17:06:07
    不同生境中入侵植物喜旱蓮子草與本地種接骨草的光合特性
    午夜视频国产福利| 我的老师免费观看完整版| 欧美xxⅹ黑人| 欧美亚洲 丝袜 人妻 在线| 国产有黄有色有爽视频| 成人毛片a级毛片在线播放| 国产乱人视频| 久久精品久久久久久噜噜老黄| av免费在线看不卡| 2018国产大陆天天弄谢| 久久青草综合色| 黄色怎么调成土黄色| 久久精品久久精品一区二区三区| 2021少妇久久久久久久久久久| 久久久久久久国产电影| 蜜桃在线观看..| 国产伦在线观看视频一区| 你懂的网址亚洲精品在线观看| 日韩av在线免费看完整版不卡| 亚洲真实伦在线观看| 亚洲中文av在线| 亚洲第一区二区三区不卡| 18禁动态无遮挡网站| 精品久久久噜噜| 秋霞伦理黄片| 国产爽快片一区二区三区| 久久久a久久爽久久v久久| 亚洲第一av免费看| 色网站视频免费| 中文字幕av成人在线电影| 日韩中字成人| 国产精品一区二区在线不卡| 观看美女的网站| 亚洲欧美日韩卡通动漫| 日日摸夜夜添夜夜爱| 国产老妇伦熟女老妇高清| 秋霞伦理黄片| av福利片在线观看| 亚洲美女视频黄频| 国产高清国产精品国产三级 | 丰满迷人的少妇在线观看| 一区二区三区四区激情视频| av卡一久久| 亚洲四区av| 国产乱来视频区| 国产深夜福利视频在线观看| 小蜜桃在线观看免费完整版高清| 亚洲中文av在线| 欧美精品人与动牲交sv欧美| 国产精品av视频在线免费观看| 少妇人妻久久综合中文| 在线观看美女被高潮喷水网站| 好男人视频免费观看在线| 欧美xxxx性猛交bbbb| 成年av动漫网址| 日本欧美视频一区| 大又大粗又爽又黄少妇毛片口| 国产成人精品福利久久| 大话2 男鬼变身卡| 在线观看av片永久免费下载| 一区在线观看完整版| 国产伦精品一区二区三区视频9| 亚洲内射少妇av| 久久6这里有精品| 久久久久精品性色| 最近最新中文字幕大全电影3| 国产成人午夜福利电影在线观看| 亚洲欧美一区二区三区黑人 | av专区在线播放| 亚洲国产精品999| 成人免费观看视频高清| 免费av不卡在线播放| 国产免费又黄又爽又色| 在线播放无遮挡| 日韩 亚洲 欧美在线| 欧美zozozo另类| 国产男人的电影天堂91| 91aial.com中文字幕在线观看| av在线观看视频网站免费| 男人舔奶头视频| 亚洲激情五月婷婷啪啪| 成人毛片60女人毛片免费| av.在线天堂| 欧美最新免费一区二区三区| 国产黄片美女视频| 国产精品一区二区在线观看99| 国产白丝娇喘喷水9色精品| av福利片在线观看| 免费人妻精品一区二区三区视频| 国产精品爽爽va在线观看网站| 午夜免费鲁丝| 天美传媒精品一区二区| 偷拍熟女少妇极品色| 伦理电影免费视频| av一本久久久久| 99热国产这里只有精品6| 久久青草综合色| 97在线视频观看| 亚洲伊人久久精品综合| 99久久中文字幕三级久久日本| 九色成人免费人妻av| 亚洲欧美精品专区久久| 全区人妻精品视频| 国产午夜精品一二区理论片| 99久久精品一区二区三区| 亚洲国产毛片av蜜桃av| 国产亚洲一区二区精品| 免费黄网站久久成人精品| 久热久热在线精品观看| 最后的刺客免费高清国语| 欧美性感艳星| 黑人高潮一二区| 亚洲国产欧美人成| 国产亚洲午夜精品一区二区久久| 91aial.com中文字幕在线观看| 中文字幕av成人在线电影| 91久久精品电影网| 亚洲av欧美aⅴ国产| 国产女主播在线喷水免费视频网站| 国产免费一级a男人的天堂| 久久久久久九九精品二区国产| 国产高潮美女av| 永久网站在线| 亚洲av日韩在线播放| 精品熟女少妇av免费看| 性高湖久久久久久久久免费观看| 久久热精品热| 国产精品人妻久久久影院| 中文乱码字字幕精品一区二区三区| 内地一区二区视频在线| 欧美xxxx性猛交bbbb| 国产精品99久久99久久久不卡 | 亚洲无线观看免费| 中文字幕av成人在线电影| 91久久精品国产一区二区三区| 麻豆乱淫一区二区| av天堂中文字幕网| 又粗又硬又长又爽又黄的视频| 中文乱码字字幕精品一区二区三区| 久久精品夜色国产| 一个人看的www免费观看视频| 日本av免费视频播放| 精品久久久久久久久av| 三级国产精品片| 亚洲av成人精品一区久久| av天堂中文字幕网| 各种免费的搞黄视频| 国产成人a∨麻豆精品| 国产69精品久久久久777片| 欧美+日韩+精品| 蜜臀久久99精品久久宅男| 大码成人一级视频| 成年人午夜在线观看视频| 97在线人人人人妻| 国产成人精品福利久久| 中文字幕av成人在线电影| 久久人人爽av亚洲精品天堂 | 成人高潮视频无遮挡免费网站| 国产爱豆传媒在线观看| 亚洲美女视频黄频| 精品人妻熟女av久视频| 人人妻人人看人人澡| 三级国产精品欧美在线观看| 国产精品欧美亚洲77777| 插阴视频在线观看视频| 热99国产精品久久久久久7| 中文字幕免费在线视频6| 国产探花极品一区二区| 欧美日韩视频精品一区| 色视频在线一区二区三区| 国产黄色视频一区二区在线观看| 日日啪夜夜爽| 国产精品久久久久久久电影| 99热网站在线观看| av卡一久久| 99热这里只有精品一区| 久久久久久久久久久免费av| 最近最新中文字幕免费大全7| av在线播放精品| 三级国产精品片| 日韩欧美一区视频在线观看 | 日本黄大片高清| 搡女人真爽免费视频火全软件| 国精品久久久久久国模美| 国产中年淑女户外野战色| 免费黄频网站在线观看国产| 亚洲精品自拍成人| 国产精品一区二区在线不卡| 九九爱精品视频在线观看| 日本黄大片高清| 青春草亚洲视频在线观看| 欧美极品一区二区三区四区| 熟女电影av网| av又黄又爽大尺度在线免费看| 亚洲精品久久午夜乱码| 美女视频免费永久观看网站| 超碰av人人做人人爽久久| 国语对白做爰xxxⅹ性视频网站| 少妇人妻精品综合一区二区| 建设人人有责人人尽责人人享有的 | 热re99久久精品国产66热6| 国产乱人偷精品视频| 人妻 亚洲 视频| 麻豆成人av视频| 欧美变态另类bdsm刘玥| 国产精品99久久久久久久久| 女性生殖器流出的白浆| 老司机影院毛片| 高清午夜精品一区二区三区| 成人亚洲欧美一区二区av| 国产探花极品一区二区| 国产欧美亚洲国产| 插逼视频在线观看| av网站免费在线观看视频| 亚洲精品视频女| 五月伊人婷婷丁香| 亚洲av.av天堂| 2018国产大陆天天弄谢| 激情 狠狠 欧美| 亚洲,一卡二卡三卡| 久久久久精品性色| 特大巨黑吊av在线直播| 日韩视频在线欧美| 亚洲色图av天堂| 亚洲综合色惰| 一本一本综合久久| 欧美变态另类bdsm刘玥| 最黄视频免费看| 国产精品三级大全| 尤物成人国产欧美一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 各种免费的搞黄视频| 在线 av 中文字幕| 久久亚洲国产成人精品v| 蜜臀久久99精品久久宅男| 国国产精品蜜臀av免费| 亚洲精品中文字幕在线视频 | 少妇人妻精品综合一区二区| 亚洲国产精品999| 日本av免费视频播放| 国精品久久久久久国模美| 亚洲精品乱码久久久v下载方式| 精品国产一区二区三区久久久樱花 | av播播在线观看一区| 偷拍熟女少妇极品色| 久久久久精品性色| 80岁老熟妇乱子伦牲交| 一级毛片 在线播放| 男人和女人高潮做爰伦理| 国产精品伦人一区二区| 国产av码专区亚洲av| 99久久中文字幕三级久久日本| 99久国产av精品国产电影| 亚洲av中文字字幕乱码综合| 国产黄频视频在线观看| 天天躁日日操中文字幕| 一级毛片aaaaaa免费看小| av国产免费在线观看| 简卡轻食公司| 久久人妻熟女aⅴ| 97在线视频观看| 一级黄片播放器| 国产伦理片在线播放av一区| 少妇人妻 视频| 丰满少妇做爰视频| 人妻 亚洲 视频| 深夜a级毛片| 亚洲成人手机| 国产免费又黄又爽又色| 日本与韩国留学比较| 亚洲天堂av无毛| 久久韩国三级中文字幕| 久久这里有精品视频免费| 国产精品成人在线| 久久6这里有精品| 成人午夜精彩视频在线观看| 国产 一区精品| 国产免费一区二区三区四区乱码| 国产在线男女| 国产一区二区三区av在线| 日本欧美国产在线视频| 久热这里只有精品99| 91aial.com中文字幕在线观看| 天天躁日日操中文字幕| 亚洲人成网站在线观看播放| 久久精品国产自在天天线| 日本一二三区视频观看| 人人妻人人爽人人添夜夜欢视频 | 精品人妻偷拍中文字幕| 日韩成人av中文字幕在线观看| 久久精品国产鲁丝片午夜精品| 国产av码专区亚洲av| 夫妻午夜视频| 欧美极品一区二区三区四区| 国产黄色免费在线视频| 卡戴珊不雅视频在线播放| 精品亚洲成a人片在线观看 | 日韩强制内射视频| 国产一区有黄有色的免费视频| 国产爱豆传媒在线观看| 亚洲综合精品二区| 麻豆成人av视频| 精品久久国产蜜桃| 亚洲av不卡在线观看| 身体一侧抽搐| 精品一区在线观看国产| 极品少妇高潮喷水抽搐| 成人黄色视频免费在线看| 18禁裸乳无遮挡免费网站照片| 成人综合一区亚洲| 国产伦精品一区二区三区视频9| 男人和女人高潮做爰伦理| 免费人成在线观看视频色| 亚洲综合色惰| 在线免费观看不下载黄p国产| 国产在线免费精品| 亚洲内射少妇av| 日本黄色片子视频| 中文字幕av成人在线电影| 欧美zozozo另类| 国产色婷婷99| 欧美日韩精品成人综合77777| 久久久久精品性色| 亚洲欧美日韩东京热| 亚州av有码| 国产乱来视频区| 夫妻性生交免费视频一级片| 大片电影免费在线观看免费| 成年免费大片在线观看| 制服丝袜香蕉在线| videossex国产| 成人综合一区亚洲| 又大又黄又爽视频免费| 春色校园在线视频观看| 久久久久久久久久久免费av| 一级毛片电影观看| av线在线观看网站| 国产爱豆传媒在线观看| videossex国产| 小蜜桃在线观看免费完整版高清| 亚洲伊人久久精品综合| 久久精品久久久久久久性| 精品少妇黑人巨大在线播放| 久久久久久伊人网av| 最近的中文字幕免费完整| 成人高潮视频无遮挡免费网站| 777米奇影视久久| 国产精品蜜桃在线观看| 国产成人精品久久久久久| 一级二级三级毛片免费看| 国产精品伦人一区二区| 国产永久视频网站| 嫩草影院入口| 成人一区二区视频在线观看| 久久精品国产自在天天线| 国内精品宾馆在线| 大又大粗又爽又黄少妇毛片口| 成人国产av品久久久| 又黄又爽又刺激的免费视频.| 国产爱豆传媒在线观看| 欧美成人a在线观看| 久久久久久久久久成人| 日韩三级伦理在线观看| 久久久久久人妻| 亚洲欧美一区二区三区黑人 | 国产精品国产av在线观看| 久久久久精品久久久久真实原创| 国产精品久久久久成人av| 在线播放无遮挡| 亚洲一级一片aⅴ在线观看| 国产亚洲一区二区精品| 欧美日韩一区二区视频在线观看视频在线| 天堂8中文在线网| 久久精品国产a三级三级三级| 成人毛片a级毛片在线播放| 精品国产三级普通话版| 水蜜桃什么品种好| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲5aaaaa淫片| 日韩免费高清中文字幕av| 亚洲最大成人中文| 免费久久久久久久精品成人欧美视频 | 日韩欧美 国产精品| 亚洲va在线va天堂va国产| 欧美日韩综合久久久久久| 久久国产亚洲av麻豆专区| 免费看不卡的av| 亚州av有码| 夜夜爽夜夜爽视频| 久久人人爽人人爽人人片va| 国产精品99久久99久久久不卡 | 九九爱精品视频在线观看| 国产一区二区在线观看日韩| 国产v大片淫在线免费观看| 寂寞人妻少妇视频99o| 高清欧美精品videossex| 亚洲国产高清在线一区二区三| 成人无遮挡网站| 精品酒店卫生间| 久久人人爽av亚洲精品天堂 | 国语对白做爰xxxⅹ性视频网站| 岛国毛片在线播放| 久久综合国产亚洲精品| 久久鲁丝午夜福利片| 中文精品一卡2卡3卡4更新| 男人狂女人下面高潮的视频| 深夜a级毛片| 国产精品久久久久久精品古装| 精品少妇黑人巨大在线播放| 高清av免费在线| av播播在线观看一区| 波野结衣二区三区在线| a级毛片免费高清观看在线播放| 另类亚洲欧美激情| 久久久久久久精品精品| 在线播放无遮挡| 国产精品国产三级专区第一集| 色视频在线一区二区三区| 亚洲av免费高清在线观看| 国产爱豆传媒在线观看| 五月天丁香电影| 极品教师在线视频| 女人久久www免费人成看片| 人体艺术视频欧美日本| 只有这里有精品99| 麻豆精品久久久久久蜜桃| 日韩在线高清观看一区二区三区| 成人免费观看视频高清| 黄色欧美视频在线观看| 成人特级av手机在线观看| 亚洲欧美一区二区三区国产| 亚洲人成网站在线观看播放| 自拍偷自拍亚洲精品老妇| 观看美女的网站| 免费黄色在线免费观看| 亚洲欧美成人精品一区二区| 国产 一区精品| 国产中年淑女户外野战色| 亚洲第一区二区三区不卡| 久久久久精品久久久久真实原创| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩亚洲高清精品| 国内少妇人妻偷人精品xxx网站| 日韩强制内射视频| 成人二区视频| 少妇高潮的动态图| 色吧在线观看| 国产精品一区www在线观看| 内射极品少妇av片p| 搡女人真爽免费视频火全软件| 国产精品久久久久成人av| 精品酒店卫生间| 国产精品欧美亚洲77777| av在线app专区| 亚洲精品,欧美精品| 男人爽女人下面视频在线观看| 美女脱内裤让男人舔精品视频| 啦啦啦视频在线资源免费观看| 国产欧美亚洲国产| 丰满迷人的少妇在线观看| 国产av一区二区精品久久 | 日韩大片免费观看网站| 国产乱人偷精品视频| 日韩中文字幕视频在线看片 | 国产黄色免费在线视频| 日本欧美国产在线视频| 交换朋友夫妻互换小说| 能在线免费看毛片的网站| 亚洲av日韩在线播放| 夜夜骑夜夜射夜夜干| 久久久久久久久久人人人人人人| 综合色丁香网| 国产午夜精品一二区理论片| 亚洲国产欧美人成| 又爽又黄a免费视频| 亚洲精品乱久久久久久| av在线app专区| 亚洲欧美成人综合另类久久久| 欧美成人一区二区免费高清观看| 国产又色又爽无遮挡免| 黄片wwwwww| 超碰97精品在线观看| 国产毛片在线视频| 精品一区在线观看国产| 免费在线观看成人毛片| 大香蕉久久网| 十分钟在线观看高清视频www | 日日啪夜夜爽| 26uuu在线亚洲综合色| 99热这里只有是精品在线观看| 2022亚洲国产成人精品| 久久精品国产鲁丝片午夜精品| 80岁老熟妇乱子伦牲交| 国产精品免费大片| 久久久久久久国产电影| 久久久久精品久久久久真实原创| 久久6这里有精品| 夜夜看夜夜爽夜夜摸| 久久国内精品自在自线图片| 亚洲电影在线观看av| 啦啦啦中文免费视频观看日本| 18+在线观看网站| 国产精品三级大全| 99热这里只有精品一区| 国产一区二区三区av在线| 日本黄色片子视频| 日本与韩国留学比较| 欧美xxxx性猛交bbbb| av福利片在线观看| 久久99热这里只频精品6学生| 久久久欧美国产精品| 久久久久视频综合| 亚洲国产欧美人成| 啦啦啦啦在线视频资源| 国产av精品麻豆| 青春草亚洲视频在线观看| av在线app专区| 少妇的逼水好多| www.色视频.com| 国产精品伦人一区二区| 波野结衣二区三区在线| freevideosex欧美| 最后的刺客免费高清国语| 欧美+日韩+精品| 久久人人爽人人爽人人片va| 亚洲成人一二三区av| 91aial.com中文字幕在线观看| 99久国产av精品国产电影| 国产在线免费精品| 免费观看av网站的网址| 国产亚洲精品久久久com| www.av在线官网国产| 91久久精品国产一区二区三区| 成年av动漫网址| 国产永久视频网站| 亚洲va在线va天堂va国产| 国产成人精品久久久久久| 久久久国产一区二区| 啦啦啦啦在线视频资源| 国产极品天堂在线| 久久精品国产鲁丝片午夜精品| 少妇人妻久久综合中文| 人妻系列 视频| 中文字幕av成人在线电影| 国产美女午夜福利| 免费高清在线观看视频在线观看| a级毛色黄片| 黄色日韩在线| 哪个播放器可以免费观看大片| 一区二区三区乱码不卡18| 亚洲,欧美,日韩| 在线观看国产h片| 日韩视频在线欧美| 99热国产这里只有精品6| 免费黄频网站在线观看国产| 日韩三级伦理在线观看| av在线老鸭窝| 日日撸夜夜添| 一个人免费看片子| av又黄又爽大尺度在线免费看| 高清毛片免费看| 亚洲真实伦在线观看| 六月丁香七月| 亚洲成人av在线免费| 插阴视频在线观看视频| 三级国产精品欧美在线观看| 丰满少妇做爰视频| 免费人妻精品一区二区三区视频| 丰满少妇做爰视频| 我的老师免费观看完整版| 欧美xxxx黑人xx丫x性爽| 一级毛片黄色毛片免费观看视频| 亚洲婷婷狠狠爱综合网| 日日啪夜夜撸| 日本欧美国产在线视频| 精品人妻偷拍中文字幕| 在线免费十八禁| 91午夜精品亚洲一区二区三区| 国产亚洲一区二区精品| 肉色欧美久久久久久久蜜桃| 久久久色成人| 精品人妻偷拍中文字幕| 亚洲精品成人av观看孕妇| 国产伦精品一区二区三区视频9| 精品人妻一区二区三区麻豆| 大码成人一级视频| 精品亚洲成国产av| 精品久久久久久久久av| 日日摸夜夜添夜夜添av毛片| 伦理电影免费视频| 亚洲精华国产精华液的使用体验| 欧美xxxx黑人xx丫x性爽| 九九在线视频观看精品| 久久久色成人| av天堂中文字幕网| 亚洲av免费高清在线观看| 色综合色国产| 午夜老司机福利剧场| 两个人的视频大全免费| 国产乱来视频区| 国产午夜精品久久久久久一区二区三区| 亚洲不卡免费看| 国产乱人视频| 日本黄色片子视频| 九色成人免费人妻av| 在线观看人妻少妇| 国产亚洲欧美精品永久| 国产色婷婷99| 国产高清国产精品国产三级 | 久久久精品免费免费高清| 一二三四中文在线观看免费高清| 永久网站在线| 国产淫片久久久久久久久| 亚洲精品日韩在线中文字幕| 美女高潮的动态|