• <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
    不同生境中入侵植物喜旱蓮子草與本地種接骨草的光合特性
    无遮挡黄片免费观看| 精品不卡国产一区二区三区| 99久国产av精品| 美女内射精品一级片tv| 国产精品久久久久久精品电影| 国产精品无大码| 久久久成人免费电影| 天堂√8在线中文| 91在线精品国自产拍蜜月| 欧美日韩在线观看h| 99久久精品热视频| 露出奶头的视频| 国产精品一及| 最新中文字幕久久久久| 亚洲欧美日韩无卡精品| 久久精品夜夜夜夜夜久久蜜豆| www.色视频.com| 熟女电影av网| 国产精品乱码一区二三区的特点| 人妻丰满熟妇av一区二区三区| 精品日产1卡2卡| 国产成人aa在线观看| 亚洲久久久久久中文字幕| 免费人成在线观看视频色| 国产精品不卡视频一区二区| 日韩欧美国产在线观看| 免费一级毛片在线播放高清视频| 亚洲最大成人中文| 一个人看视频在线观看www免费| 高清午夜精品一区二区三区 | 亚洲欧美日韩东京热| 国产精品久久久久久久电影| ponron亚洲| 成年女人看的毛片在线观看| 国产美女午夜福利| 国产美女午夜福利| 欧美3d第一页| 国产乱人偷精品视频| 国产精品久久久久久久久免| 国内久久婷婷六月综合欲色啪| 成人二区视频| 中文字幕人妻熟人妻熟丝袜美| 看免费成人av毛片| 亚洲精品亚洲一区二区| 精品午夜福利在线看| 国模一区二区三区四区视频| 中出人妻视频一区二区| 国产精品电影一区二区三区| 国产精品亚洲一级av第二区| 一进一出抽搐gif免费好疼| 寂寞人妻少妇视频99o| 美女大奶头视频| 床上黄色一级片| 毛片一级片免费看久久久久| 成人av在线播放网站| 亚洲人成网站在线播| 一级毛片aaaaaa免费看小| 亚洲人成网站在线观看播放| 看片在线看免费视频| 黄色日韩在线| 国产亚洲精品久久久久久毛片| 国产熟女欧美一区二区| 黄色配什么色好看| 国产精品嫩草影院av在线观看| 色在线成人网| 日产精品乱码卡一卡2卡三| www日本黄色视频网| 欧美bdsm另类| 久久久久久伊人网av| 日韩欧美免费精品| 午夜免费激情av| 99热6这里只有精品| 在线免费观看不下载黄p国产| 色视频www国产| 国产色爽女视频免费观看| 国产亚洲精品久久久com| 精品久久久久久成人av| 最近最新中文字幕大全电影3| 淫妇啪啪啪对白视频| 91午夜精品亚洲一区二区三区| 国产高清不卡午夜福利| 一区二区三区高清视频在线| 男女下面进入的视频免费午夜| 亚洲欧美日韩无卡精品| 热99在线观看视频| 国产精品久久久久久精品电影| 97在线视频观看| 女生性感内裤真人,穿戴方法视频| 天美传媒精品一区二区| 此物有八面人人有两片| 国产成人精品久久久久久| 久久热精品热| 热99在线观看视频| 美女免费视频网站| 五月伊人婷婷丁香| 一a级毛片在线观看| 免费高清视频大片| 日本爱情动作片www.在线观看 | 综合色丁香网| 1024手机看黄色片| 亚洲第一电影网av| 中文资源天堂在线| 国产精品综合久久久久久久免费| 此物有八面人人有两片| 欧美高清性xxxxhd video| 真实男女啪啪啪动态图| 中文字幕久久专区| 99热精品在线国产| 免费无遮挡裸体视频| 看非洲黑人一级黄片| 国产精品精品国产色婷婷| 久久久久性生活片| 看黄色毛片网站| 悠悠久久av| 少妇的逼好多水| 一本久久中文字幕| 中文字幕人妻熟人妻熟丝袜美| 国产高清不卡午夜福利| 久久精品夜夜夜夜夜久久蜜豆| 69av精品久久久久久| 国产亚洲欧美98| 欧美+日韩+精品| 人人妻人人澡欧美一区二区| 听说在线观看完整版免费高清| 在线看三级毛片| 看片在线看免费视频| 精品国产三级普通话版| 成熟少妇高潮喷水视频| 国产精华一区二区三区| 国产色爽女视频免费观看| 在线国产一区二区在线| 日产精品乱码卡一卡2卡三| 久久久精品欧美日韩精品| 能在线免费观看的黄片| 免费av观看视频| 自拍偷自拍亚洲精品老妇| 日本免费a在线| 欧美丝袜亚洲另类| 亚洲欧美成人精品一区二区| 亚洲精品亚洲一区二区| 亚洲最大成人av| 国产日本99.免费观看| 中文资源天堂在线| 美女大奶头视频| 两个人视频免费观看高清| 九九爱精品视频在线观看| 亚洲丝袜综合中文字幕| 男女啪啪激烈高潮av片| or卡值多少钱| 非洲黑人性xxxx精品又粗又长| 日日摸夜夜添夜夜添小说| 日韩精品青青久久久久久| 免费观看人在逋| 99久国产av精品| 日日摸夜夜添夜夜添av毛片| 又黄又爽又免费观看的视频| 亚洲精品日韩在线中文字幕 | 99久久九九国产精品国产免费| 麻豆成人午夜福利视频| 97超视频在线观看视频| 亚洲国产精品sss在线观看| 男女边吃奶边做爰视频| 激情 狠狠 欧美| 午夜福利18| 免费电影在线观看免费观看| 香蕉av资源在线| 高清毛片免费观看视频网站| 国产欧美日韩精品一区二区| 嫩草影视91久久| 深夜精品福利| 免费观看的影片在线观看| 一级毛片aaaaaa免费看小| ponron亚洲| 久久久久久久午夜电影| 99热6这里只有精品| 成人永久免费在线观看视频| 国产一区二区在线观看日韩| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品乱码久久久v下载方式| 成年版毛片免费区| 欧美又色又爽又黄视频| 日本精品一区二区三区蜜桃| 两个人视频免费观看高清| 亚洲图色成人| 熟女电影av网| 精品久久久久久久久av| 三级毛片av免费| 91在线观看av| 午夜影院日韩av| 色综合站精品国产| 国产片特级美女逼逼视频| 久久久a久久爽久久v久久| 黄色一级大片看看| av福利片在线观看| 国产成人a区在线观看| 国产精品电影一区二区三区| 国产伦一二天堂av在线观看| 亚洲人成网站在线播放欧美日韩| 97超视频在线观看视频| 久99久视频精品免费| 天堂网av新在线| 不卡一级毛片| 91麻豆精品激情在线观看国产| 国产成人影院久久av| 亚洲精品456在线播放app| 国产精品电影一区二区三区| 在线观看66精品国产| 国产成人91sexporn| 男插女下体视频免费在线播放| 日本在线视频免费播放| 黄色日韩在线| 人人妻,人人澡人人爽秒播| 床上黄色一级片| 国产精华一区二区三区| 免费观看精品视频网站| 日本成人三级电影网站| 国产成人91sexporn| 色播亚洲综合网| 亚洲欧美日韩卡通动漫| 在线a可以看的网站| 国产av在哪里看| 91久久精品国产一区二区三区| 国内精品宾馆在线| 老女人水多毛片| 日韩在线高清观看一区二区三区| 国产av在哪里看| 久久午夜福利片| av在线老鸭窝| 亚洲av五月六月丁香网| 男女做爰动态图高潮gif福利片| 久久午夜亚洲精品久久| 在线播放无遮挡| 成人三级黄色视频| 日本免费a在线| 美女高潮的动态| 亚洲精品一卡2卡三卡4卡5卡| 日韩欧美三级三区| 国产日本99.免费观看| 免费一级毛片在线播放高清视频| 亚洲自拍偷在线| 99国产精品一区二区蜜桃av| 老师上课跳d突然被开到最大视频| 国产高潮美女av| 久久久精品欧美日韩精品| 免费看日本二区| 久久亚洲国产成人精品v| 久久久久久九九精品二区国产| 高清日韩中文字幕在线| 春色校园在线视频观看| 中出人妻视频一区二区| 小蜜桃在线观看免费完整版高清| av中文乱码字幕在线| 国产一区二区三区在线臀色熟女| 久久久久国内视频| 日韩欧美国产在线观看| 不卡一级毛片| 麻豆国产97在线/欧美| 欧美xxxx黑人xx丫x性爽| 久久久国产成人精品二区| 人人妻人人看人人澡| 又粗又爽又猛毛片免费看| 成年女人永久免费观看视频| 午夜激情欧美在线| 桃色一区二区三区在线观看| 亚洲av中文字字幕乱码综合| 少妇熟女aⅴ在线视频| 三级经典国产精品| 国产aⅴ精品一区二区三区波| 成人亚洲欧美一区二区av| 搡女人真爽免费视频火全软件 | 亚洲七黄色美女视频| 久久精品夜色国产| 国产午夜精品论理片| 身体一侧抽搐| 一个人免费在线观看电影| 啦啦啦啦在线视频资源| 欧美潮喷喷水| 村上凉子中文字幕在线| 日韩欧美免费精品| 又爽又黄a免费视频| 男女啪啪激烈高潮av片| 久久久久久国产a免费观看| 麻豆乱淫一区二区| 成人一区二区视频在线观看| 亚洲专区国产一区二区| 久久久成人免费电影| 最近中文字幕高清免费大全6| 亚洲欧美清纯卡通| 熟妇人妻久久中文字幕3abv| 久久精品国产清高在天天线| 亚洲欧美日韩卡通动漫| 人妻久久中文字幕网| 一级黄色大片毛片| 精品人妻视频免费看| 国产精品三级大全| 亚洲国产精品成人久久小说 | 麻豆国产av国片精品| 国产精品嫩草影院av在线观看| 国产男靠女视频免费网站| 国产乱人视频| 91av网一区二区| 欧美性感艳星| 免费黄网站久久成人精品| 国产女主播在线喷水免费视频网站 | 有码 亚洲区| 99热6这里只有精品| 青春草视频在线免费观看| 亚洲av五月六月丁香网| 黑人高潮一二区| 夜夜爽天天搞| 黄色欧美视频在线观看| 搡老岳熟女国产| 亚洲精品久久国产高清桃花| 成人欧美大片| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美成人综合另类久久久 | 看黄色毛片网站| 日本-黄色视频高清免费观看| 精品人妻偷拍中文字幕| 亚洲18禁久久av| 美女cb高潮喷水在线观看| 国产色爽女视频免费观看| 色在线成人网| 免费在线观看成人毛片| 久久热精品热| 能在线免费观看的黄片| 不卡一级毛片| 男女之事视频高清在线观看| 日本熟妇午夜| 欧美一区二区亚洲| 91狼人影院| 欧美性猛交黑人性爽| 国产综合懂色| 美女大奶头视频| 韩国av在线不卡| 日本免费一区二区三区高清不卡| 免费观看人在逋| 精品久久久久久久久久久久久| 免费无遮挡裸体视频| 丰满的人妻完整版| 亚洲欧美成人综合另类久久久 | 女人十人毛片免费观看3o分钟| 男人狂女人下面高潮的视频| 97超碰精品成人国产| 国产男靠女视频免费网站| 俄罗斯特黄特色一大片| 亚洲图色成人| 又粗又爽又猛毛片免费看| 俄罗斯特黄特色一大片| 欧美极品一区二区三区四区| 夜夜看夜夜爽夜夜摸| 嫩草影院精品99| 欧美3d第一页| 国语自产精品视频在线第100页| 美女大奶头视频| 亚洲欧美日韩卡通动漫| 亚洲不卡免费看| 成人亚洲欧美一区二区av| 精品国产三级普通话版| 色av中文字幕| 亚洲精品成人久久久久久| 中文字幕精品亚洲无线码一区| 免费观看在线日韩| 午夜老司机福利剧场| 欧美日韩综合久久久久久| 日日摸夜夜添夜夜爱| 欧美日韩国产亚洲二区| 乱人视频在线观看| 波多野结衣巨乳人妻| 亚洲熟妇中文字幕五十中出| 欧美一级a爱片免费观看看| 国产精品人妻久久久影院| 在线播放国产精品三级| 久久婷婷人人爽人人干人人爱| 亚洲国产精品久久男人天堂| 搡女人真爽免费视频火全软件 | 国产成人精品久久久久久| 最近手机中文字幕大全| 你懂的网址亚洲精品在线观看 | 一级毛片电影观看 | 久久九九热精品免费| 国产亚洲精品久久久久久毛片| 欧美激情在线99| 午夜视频国产福利| 偷拍熟女少妇极品色| 中文字幕av在线有码专区| 99热网站在线观看| 又黄又爽又刺激的免费视频.| 国产成人91sexporn| 亚洲第一电影网av| 亚洲自偷自拍三级| 成人欧美大片| 99九九线精品视频在线观看视频| 国产伦一二天堂av在线观看| 亚洲在线观看片| 国产一区二区在线观看日韩| 欧美极品一区二区三区四区| 日韩欧美精品v在线| 成年女人毛片免费观看观看9| 深夜a级毛片| 性欧美人与动物交配| 日本五十路高清| 最近在线观看免费完整版| 身体一侧抽搐| 一进一出抽搐动态| 蜜臀久久99精品久久宅男| 国产精品一二三区在线看| 男人舔奶头视频| 国产中年淑女户外野战色| 人妻夜夜爽99麻豆av| 免费观看的影片在线观看| 插逼视频在线观看| 最近手机中文字幕大全| 我的老师免费观看完整版| 国产精华一区二区三区| 亚洲经典国产精华液单| 日本欧美国产在线视频| 色哟哟哟哟哟哟| 国产片特级美女逼逼视频| 久久久久久九九精品二区国产| 日本一二三区视频观看| 一区二区三区高清视频在线| 午夜激情福利司机影院| 国产男靠女视频免费网站| 一级av片app| 亚洲精品日韩av片在线观看| 国产精品伦人一区二区| 国产午夜精品久久久久久一区二区三区 | 精品久久久久久久久久免费视频| 变态另类丝袜制服| 国产亚洲精品综合一区在线观看| 国产精品女同一区二区软件| 国产精品一区二区三区四区免费观看 | 一级黄片播放器| 日韩,欧美,国产一区二区三区 | 岛国在线免费视频观看| 免费高清视频大片| 久久久国产成人精品二区| 成年女人看的毛片在线观看| 中文字幕熟女人妻在线| 99热这里只有精品一区| 国产三级中文精品| 日韩高清综合在线| 91狼人影院| 一进一出好大好爽视频| 成年女人毛片免费观看观看9| 欧美日韩一区二区视频在线观看视频在线 | 干丝袜人妻中文字幕| 91av网一区二区| 色尼玛亚洲综合影院| 香蕉av资源在线| 成人美女网站在线观看视频| 九九热线精品视视频播放| 中国美白少妇内射xxxbb| 免费观看人在逋| 一个人看的www免费观看视频| 日韩高清综合在线| 男女啪啪激烈高潮av片| 久久草成人影院| 国产伦精品一区二区三区视频9| 老师上课跳d突然被开到最大视频| 黑人高潮一二区| 日日摸夜夜添夜夜添小说| 最近2019中文字幕mv第一页| ponron亚洲| 在线免费观看不下载黄p国产| 亚洲自拍偷在线| 亚洲不卡免费看| 国产不卡一卡二| 十八禁国产超污无遮挡网站| 亚洲精品粉嫩美女一区| 亚洲,欧美,日韩| 久久韩国三级中文字幕| 日本-黄色视频高清免费观看| 美女 人体艺术 gogo| 国产 一区精品| 亚洲精品乱码久久久v下载方式| 精品乱码久久久久久99久播| 午夜免费男女啪啪视频观看 | av专区在线播放| 五月玫瑰六月丁香| 欧美成人精品欧美一级黄| 中文字幕免费在线视频6| 国产精品无大码| 韩国av在线不卡| 国产乱人视频| 国产黄色视频一区二区在线观看 | 久久6这里有精品| 国产v大片淫在线免费观看| 国产午夜精品久久久久久一区二区三区 | 久久天躁狠狠躁夜夜2o2o| 午夜日韩欧美国产| 久久久久久伊人网av| 男女之事视频高清在线观看| 蜜桃久久精品国产亚洲av| 亚洲美女黄片视频| 亚洲人与动物交配视频| 国产在线男女| 日韩中字成人| 一级黄片播放器| 国内精品美女久久久久久| 日韩强制内射视频| 美女高潮的动态| 两个人的视频大全免费| 欧美日韩一区二区视频在线观看视频在线 | 少妇熟女欧美另类| 久久鲁丝午夜福利片| 免费高清视频大片| 成人午夜高清在线视频| 99热精品在线国产| 神马国产精品三级电影在线观看| 综合色丁香网| 18+在线观看网站| 亚洲中文字幕日韩| 免费人成视频x8x8入口观看| 精品久久久久久久人妻蜜臀av| 我的女老师完整版在线观看| 久久久国产成人免费| 插逼视频在线观看| 噜噜噜噜噜久久久久久91| 天堂网av新在线| av免费在线看不卡| 国内久久婷婷六月综合欲色啪| 久久精品夜夜夜夜夜久久蜜豆| 麻豆av噜噜一区二区三区| 国内少妇人妻偷人精品xxx网站| 亚洲高清免费不卡视频| 亚洲,欧美,日韩| 中文字幕熟女人妻在线| 日韩欧美在线乱码| 日韩欧美精品v在线| 麻豆成人午夜福利视频| 欧美激情国产日韩精品一区| 床上黄色一级片| 天堂影院成人在线观看| 日韩欧美三级三区| 国产精品一区www在线观看| 久久久久久九九精品二区国产| 偷拍熟女少妇极品色| 国产精品一二三区在线看| АⅤ资源中文在线天堂| 亚洲国产精品国产精品| 国产高清激情床上av| 女的被弄到高潮叫床怎么办| av在线蜜桃| 久久久久九九精品影院| 亚洲av二区三区四区| 夜夜爽天天搞| 看十八女毛片水多多多| 99国产精品一区二区蜜桃av| 成人亚洲欧美一区二区av| 性色avwww在线观看| 亚洲中文字幕一区二区三区有码在线看| 久久综合国产亚洲精品| 亚洲成av人片在线播放无| 国产一区亚洲一区在线观看| 十八禁国产超污无遮挡网站| 亚洲自偷自拍三级| 国产又黄又爽又无遮挡在线| 国产av在哪里看| 69av精品久久久久久| 人妻制服诱惑在线中文字幕| 国产精品人妻久久久影院| 国产在视频线在精品| 久久国产乱子免费精品| 精品久久久久久成人av| 2021天堂中文幕一二区在线观| 一级av片app| 久久精品国产自在天天线| 久久国内精品自在自线图片| 亚洲精品一卡2卡三卡4卡5卡| 欧美xxxx黑人xx丫x性爽| 亚洲国产日韩欧美精品在线观看| 中国国产av一级| 欧美+日韩+精品| 3wmmmm亚洲av在线观看| 午夜老司机福利剧场| 成人特级黄色片久久久久久久| 美女内射精品一级片tv| 尾随美女入室| 我要搜黄色片| 少妇熟女aⅴ在线视频| 亚洲av熟女| 少妇高潮的动态图| 夜夜爽天天搞| 全区人妻精品视频| 欧美性感艳星| 少妇高潮的动态图| 高清毛片免费观看视频网站| 日本-黄色视频高清免费观看| 99久国产av精品国产电影| 国产精品综合久久久久久久免费| 欧美丝袜亚洲另类| 国产男靠女视频免费网站| 九九爱精品视频在线观看| 热99在线观看视频| 成人综合一区亚洲| 中文字幕熟女人妻在线| 伦理电影大哥的女人| 日韩高清综合在线| 成人特级黄色片久久久久久久| 国产麻豆成人av免费视频| 无遮挡黄片免费观看| 精品99又大又爽又粗少妇毛片| 婷婷六月久久综合丁香| 亚洲天堂国产精品一区在线| 日韩三级伦理在线观看| 看黄色毛片网站| 一个人免费在线观看电影| 又黄又爽又免费观看的视频| 国产精品爽爽va在线观看网站| av专区在线播放| 美女xxoo啪啪120秒动态图| 国产亚洲精品av在线| 精品免费久久久久久久清纯|