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

    大別山五針?biāo)筛H微生物和內(nèi)生菌群落特征及功能多樣性

    2025-02-16 00:00:00項小燕阮思睿趙菲菲徐偉芳
    廣西植物 2025年1期
    關(guān)鍵詞:高通量測序多樣性

    摘 要: "為探明大別山五針?biāo)筛H微生物和內(nèi)生菌群落特征及功能多樣性,該文以大別山五針?biāo)筛H土壤和根部組織為研究對象,利用高通量測序技術(shù)分析根際微生物和根部內(nèi)生菌群落特征,并利用PICRUSt和FUNGuild軟件分別對細(xì)菌群落和真菌群落進行功能預(yù)測。結(jié)果表明:(1)根際微生物總操作分類單元(OTUs)數(shù)目和α-多樣性指數(shù)均明顯大于根部內(nèi)生菌。(2)根際微生物主要由厚壁菌門(如李斯特氏菌屬)、變形菌門(如醋桿菌屬)、子囊菌門(如鐮刀菌屬和毛殼屬)等菌群構(gòu)成,而根部內(nèi)生菌則以厚壁菌門(如乳桿菌屬和芽孢菌屬)、子囊菌門(如枝孢菌屬)、擔(dān)子菌門(如須腹菌屬和針菌屬)菌群為主。(3)功能預(yù)測結(jié)果顯示,根際細(xì)菌群落在細(xì)胞運動性、環(huán)境適應(yīng)性、氨基酸代謝等功能方面較為突出,而根部內(nèi)生細(xì)菌功能主要體現(xiàn)在酶家族、消化系統(tǒng)、能量代謝等方面;根際真菌群落和根部內(nèi)生真菌均含有共生營養(yǎng)型、腐生營養(yǎng)型、病理營養(yǎng)型,以及多種交叉營養(yǎng)型的真菌類群,并且外生菌根更多富集在根部組織中。該研究認(rèn)為,大別山五針?biāo)筛H微生物與根部內(nèi)生菌在菌群豐富度、群落結(jié)構(gòu)組成及功能多樣性上均存在差異。該研究結(jié)果為認(rèn)識大別山五針?biāo)筛H微生物和內(nèi)生菌的結(jié)構(gòu)與功能提供了依據(jù),也為后期開發(fā)促進大別山五針?biāo)缮L、防治其病害的菌劑提供了一定的理論基礎(chǔ)。

    關(guān)鍵詞: 大別山五針?biāo)桑?根際微生物, 內(nèi)生菌, 高通量測序, 多樣性

    中圖分類號: "Q948

    文獻標(biāo)識碼: "A

    文章編號: "1000-3142(2025)01-0031-13

    基金項目: "安徽省自然科學(xué)基金(2108085QH374); 安慶師范大學(xué)科研項目(100001195); 皖西南生物多樣性研究與生態(tài)保護安徽省重點實驗室開放基金(Wxn202303); 中藥研究與開發(fā)安徽省重點實驗室開放基金(AKLPDCM202310)。

    第一作者: 項小燕(1981—),博士,教授,主要從事瀕危物種資源保護與利用研究,(E-mail)xiaoyanxiang@aqnu.edu.cn。

    *通信作者: "徐偉芳,博士,副教授,研究方向為環(huán)境微生物多樣性分析,(E-mail)xuweifang@ahtcm.edu.cn。

    Community characteristics and functional diversity

    of the rhizosphere microorganism and endophyte

    from Pinus dabeshanensis

    XIANG Xiaoyan"RUAN Sirui"ZHANG Ying2, ZHAO Feifei2, XU Weifang2*

    ( 1. Key Laboratory of Biodiversity and Ecology Conservation in Southwest Anhui, College of Life Sciences, Anqing Normal University,

    Anqing 24613 Anhui, China; 2. Anhui Province Key Laboratory of Research & Development of Chinese Medicine,

    College of Pharmacy, Anhui University of Traditional Chinese Medicine, Hefei 230012, China )

    Abstract: "Pinus dabeshanensis is an endangered plant. In order to explore the community characteristics and functional diversity of plant-associated microbes, the compositions of rhizosphere microbiota and root endophyte from P. dabeshanensis were analyzed based on the high-throughput sequencing technology. And the functions of bacterial and fungal communities were predicted by using PICRUSt and FUNGuild software, respectively. The results were as follows: (1) The rhizosphere soil samples exhibited higher operational taxonomic units (OTUs) and α-diversity indices than those of the root endophyte. (2) Rhizosphere microbiota were mainly composed of microflora such as Firmicutes (genus: Listeria), Proteobacteria (genus: Acetobacter), Ascomycetes (genus: Fusarium and Chaetomium) and so on, while root endophyte were mainly composed of Firmicutes (genus: Lactobacillus and Bacillus), Ascomycetes (genus: Cladosporium), and Basidiomycetes (genus: Pleurotus and Chaetomium). (3) The functional annotation results showed that the rhizosphere microbial community was prominent in cell mobility, environmental adaptability, and amino acid metabolism, etc., while the root endophyte community was prominent in enzyme family, digestive system, and energy metabolism, etc. Both rhizosphere fungal community and root endophytes contained symbiotic, saprophytic, pathophytic and cross-trophic fungal groups, and ectomycorrhiza were more abundant in root tissue. Collectively, there were differences in flora richness, community structure and functional diversity between rhizosphere microorganism and root endophyte of P. dabeshanensis. The study contributes to knowledge of the structure and function of rhizosphere microorganism and endophyte of P. dabeshanensis, and "provides a theoretical basis for the later development of potential microbial agents for use in promoting plant growth and controlling plant diseases.

    Key words: Pinus dabeshanensis, rhizosphere microorganism, endophyte, high-throughput sequencing, diversity

    大別山五針?biāo)桑≒inus dabeshanensis)是中國特有物種,主要分布于安徽省的岳西縣和金寨縣,以及湖北省的英山縣(吳甘霖, 2016)。目前,該物種僅存200余株且處于逐漸衰退趨勢,被列為國家一級珍稀瀕危保護植物(王雷宏等,2014)。有關(guān)大別山五針?biāo)傻难芯恐饕性谖锓N競爭(王雷宏等,2014;項小燕等,2015)、種群結(jié)構(gòu)(項小燕等, 2016a)、種子萌發(fā)特性(韓建偉等,2014)、花粉活力(項小燕等, 2016b)等方面,而有關(guān)其微生物菌群的研究鮮有報道。

    植物體在自然界的實際生存狀態(tài)被認(rèn)為是微生物和植物相互依存的狀態(tài)(石晶盈等, 2006)。植物內(nèi)生菌是指廣泛存在于健康植物組織中,與宿主植物共生且不引起宿主植物病害的微生物(黃雪珍和趙龍飛, 2023)。有研究推測,植物內(nèi)生菌來自土壤,土壤微生物通過侵入宿主植物,并與宿主植物共同進化逐漸成為植物內(nèi)生菌(Compant et al., 2010; Afzal et al., 2019) 。內(nèi)生菌種群組成和數(shù)量可能與種植地土壤微生物及植物根際微生物的種群結(jié)構(gòu)和數(shù)量有關(guān)(趙娟等, 2015)。也有研究表明,植物內(nèi)生菌與根際微生物種類多樣性存在一定的相關(guān)性,一些內(nèi)生細(xì)菌源于土壤(Mcinroy amp; Kloepper, 1994)。近年來,科技工作者在桑樹(Xie et al., 2021)、三七(Tan et al., 2017)、大豆(Pérez-Jaramillo et al., 2019)等植物根際-土壤-微生物之間的相互作用方面做了大量的探索工作。大別山五針?biāo)墒蔷窒薹植荚谖覈髣e山區(qū)域的五針?biāo)山M唯一的一種植物,是松類樹種地理分布的特殊現(xiàn)象。然而,人們對這種特定生境下的大別山五針?biāo)筛滴⑸锶郝涮卣魅狈ο到y(tǒng)的認(rèn)知,根系微生物群落結(jié)構(gòu)與功能之間如何聯(lián)系尚不清楚。隨著分子生物學(xué)和測序技術(shù)的不斷發(fā)展,基于高通量測序技術(shù)研究環(huán)境微生物多樣性特征的報道日益增多(秦宇蒙等, 2021)。探明大別山五針?biāo)筛迪嚓P(guān)微生物的種群分布特征,并預(yù)測其功能,對解析大別山五針?biāo)傻纳鷳B(tài)生存機制具有重要意義。

    Reinhold-Hurek等(2015)研究表明,根際微生物與植物內(nèi)生菌是促進宿主植物生長、增強宿主植物抗逆性的兩大類群。有益根際微生物通過多種機制促進植物生長,并維持植物健康,包括促進植物營養(yǎng)吸收、增強植物抗逆性、調(diào)節(jié)植物生理,并抵御土傳病害等(Gopal amp; Gupta, 2016)。植物根系會向根際分泌大量化合物,主要包含水、離子、酶、氧氣、黏液等初級代謝和次級代謝產(chǎn)物。其中,有許多被認(rèn)為可以介導(dǎo)根際中植物-微生物間的相互作用,對植物的生長和適應(yīng)性起著至關(guān)重要的作用。此外,內(nèi)生菌不僅通過分泌一些代謝產(chǎn)物如生長素、細(xì)胞分裂素等,促進植物細(xì)胞分裂、根系發(fā)育、幼苗生長,而且還通過促進植物對土壤中的營養(yǎng)物質(zhì)吸收而對宿主的生長發(fā)育產(chǎn)生積極影響(劉麗輝等, 2020)。盡管越來越多的研究表明根系相關(guān)微生物為其宿主植物提供了許多有益的功能,但是大別山五針?biāo)筛H土壤菌群結(jié)構(gòu)和根部內(nèi)生菌群落分布特征如何?有哪些微生物類群可能成為促進大別山五針?biāo)缮L發(fā)育的重要資源?這是本研究旨在解決的兩大科學(xué)問題。

    因此,本研究采集大別山五針?biāo)傻母H土壤和根部組織樣品,利用高通量測序技術(shù),對根際土壤和根內(nèi)部中的細(xì)菌和真菌群落組成進行系統(tǒng)鑒定,分析定殖在大別山五針?biāo)筛抵械膬?yōu)勢微生物類群,并對其生物功能進行探討,以期揭示大別山五針?biāo)筛肯嚓P(guān)微生物的群落分布特征,為解析大別山五針?biāo)傻沫h(huán)境適應(yīng)性以及人工復(fù)壯大別山五針?biāo)煞N群提供重要理論支撐。

    1 材料與方法

    1.1 樣品采集

    2020年11月上旬,于大別山五針?biāo)勺匀环N群(115°24′ E、30°44′ N)分布地安徽省岳西縣大王溝海拔1 050~1 100 m范圍內(nèi),采集陰坡健康植株樣品。采集樣品時,隨機選取3個3 m × 3 m的樣方,在每個樣方中選取1株健康植株。去除土壤表面落葉和雜草后,在植株根部一側(cè)小心挖掘,待根系露出后,挖掘部分健康根系(長度約為10 cm),采用“抖根法”(王香生等,2022),收集大別山五針?biāo)筛H東南西北4個不同方向土壤混合樣本,編號為CT,作為根際土壤DNA樣本;取同一植株東南西北4個不同方向的混合根樣為根部DNA樣本,編號為CR。待所有樣本收集完畢,立即放入冰盒帶回實驗室。

    1.2 基因組DNA的提取

    對根際土壤樣品,依據(jù)PowerSoil"DNA Isolation kit(MoBio, U.S.)試劑盒說明書中的方法提取總DNA。根部樣品經(jīng)75%酒精表面消毒后,利用改良的十六烷基三甲基溴化銨(CTAB)方法(宮強等,2005)提取總DNA。利用瓊脂糖凝膠電泳檢測DNA的純度,取適量的樣本DNA于離心管中,使用無菌水稀釋樣本至1 ng·μL-1。

    1.3 PCR擴增

    以稀釋后的基因組DNA為模板,根據(jù)測序區(qū)域的選擇,使用帶Barcode的特異引物和New England Biolabs公司的Phusion"High-Fidelity PCR Master Mix with GCBuffer以及高效高保真酶進行PCR擴增。其中,細(xì)菌采用16S rDNA V4區(qū)基因(515F: 5′-GTGCCAGCMGCCGCGGTAA-3′和806R: 5′-GGACTACHVGGGTWTCTAAT-3′),真菌采用ITS1區(qū)基因(ITS5-1737F: 5′-GGAAGTAAAAGTCGT

    AACAAGG-3′和ITS2-2043R: 5′-GCTGCGTTCTTCAT

    CGATGC-3′)進行擴增。PCR反應(yīng)體系:Phusion Master Mix(2×)15 μL,Primer(2 μmol·-1)3 μL,gDNA(1 ng·μL-1)1 μL,H2O 11 μL。反應(yīng)程序:98 ℃預(yù)變性1 min,30個循環(huán)(98 ℃, 10 s;50 ℃, 30 s;72 ℃,30 s);72 ℃,5 min。采用Bio-rad T100梯度PCR儀,PCR產(chǎn)物使用2%濃度的瓊脂糖凝膠電泳檢測,根據(jù)PCR產(chǎn)物濃度進行等濃度混樣,充分混勻后使用1×TAE濃度2%的瓊脂糖膠電泳純化PCR產(chǎn)物。其中,產(chǎn)物純化試劑盒使用Thermo Scientific公司GeneJET膠回收試劑盒。選擇主帶大小在400~450 bp之間的序列,割膠回收目標(biāo)條帶。

    1.4 文庫構(gòu)建和上機測序

    使用TruSeq"DNAPCR-Free Sample Preparation Kit建庫試劑盒進行文庫構(gòu)建,構(gòu)建好的文庫經(jīng)過Qubit和Q-PCR定量,文庫合格后,使用Nova Seq6000進行上機測序。所得樣品DNA委托北京諾禾致源生物信息科技有限公司完成測序。

    1.5 測序數(shù)據(jù)處理分析

    根據(jù)Barcode序列和PCR擴增引物序列,從下機數(shù)據(jù)中拆分出各樣本數(shù)據(jù),截去Barcode和引物序列后,使用FLAS軟件進行拼接,得到原始序列數(shù)據(jù)。經(jīng)過嚴(yán)格的過濾處理(Bokulich et al., 2013)后,參照Qiime(Caporaso et al., 2010)的序列質(zhì)量控制流程,通過與物種注釋數(shù)據(jù)庫進行比對(Rognes et al., 2016),檢測且去除嵌合體序列,得到最終的有效數(shù)據(jù)。利用Uparse軟件對所有樣本的全部有效數(shù)據(jù)進行聚類,以97%的一致性將序列聚類成為操作分類單元(operational taxonomic unit, OTU),同時選取出現(xiàn)頻數(shù)最高的序列作為OTUs的代表序列。通過對OTUs進行Venn圖分析和α-多樣性計算,得到不同樣本共有和特有OTUs信息以及物種豐富度和均勻度信息。利用獨立樣本t檢驗分析不同樣本的差異顯著性。此外,還對OTUs序列進行物種注釋,細(xì)菌用Mothur方法與SILVA138(Edgar, 2013)的SSUrRNA數(shù)據(jù)庫(Wang et al., 2007)進行物種注釋分析(設(shè)定閾值為0.8~1);真菌用Qiime軟件(Version 1.9.1)中的Blast方法(Altschul et al., 1990)與Unit(v8.2)數(shù)據(jù)庫(Kljalg et al., 2013)進行物種注釋分析,得到對應(yīng)的物種信息和基于物種的豐度分布情況。通過R軟件(Version 3.0.3)繪制PCoA圖,探究不同樣本群落結(jié)構(gòu)的差異。物種組成分析能夠反映樣品在分類學(xué)水平的群落結(jié)構(gòu),本研究在門水平上選取了豐度排名前5的物種,在屬水平上選取豐度前20的物種進行詳細(xì)的物種群落結(jié)構(gòu)組成分析。此外,本研究還利用PICRUSt軟件(KEGG數(shù)據(jù)庫)和FUNGuild軟件對大別山五針?biāo)筛H微生物總豐度前35的物種進行了功能預(yù)測。

    2 結(jié)果與分析

    2.1 測序概況分析

    經(jīng)細(xì)菌16S rRNA基因測序分析,細(xì)菌共獲得380 377條有效序列,其中3份根際土壤樣品的平均有效序列為62 202條,平均序列長度為416 bp;3份根部組織樣品的平均有效序列為64 591條,平均序列長度為409 bp。經(jīng)真菌ITS測序分析,真菌共獲得389 661條有效序列,其中3份根際土壤樣品的平均有效序列為66 142條,平均序列長度為261 bp;3份根部組織樣品的平均有效序列為63 745條,平均序列長度為268 bp。獨立樣本t檢驗分析結(jié)果顯示,根際土壤樣品和根部組織樣品有效序列數(shù)差異不顯著(P>0.05)。

    由圖1可知,細(xì)菌和真菌的稀釋曲線上升到一定高度后趨于平穩(wěn),說明測序的深度已足夠代表樣品中大多數(shù)細(xì)菌和真菌的種類,并且樣本測序數(shù)量也已足夠,可用于后續(xù)生物信息學(xué)分析。

    2.2 大別山五針?biāo)筛H微生物和內(nèi)生菌OTU分析

    依據(jù)97%相似度對所得序列進行OTU聚類,從根際土壤樣品中共獲得1 843個細(xì)菌OTU數(shù)目和821個真菌OTU數(shù)目,從根部組織中獲得1 073個內(nèi)生細(xì)菌OTU數(shù)目和750個真菌OTU數(shù)目。因此,無論是根際菌群還是根部內(nèi)生菌,其細(xì)菌OTU數(shù)目均大于真菌OTU數(shù)目。同時,我們還發(fā)現(xiàn),根際土壤OTU數(shù)目高于根部內(nèi)生菌OTU數(shù)目:(1)對于細(xì)菌來說,CT的總OTU數(shù)目(1 843)和特有的OTU數(shù)目(1 492)均大于CR的總OTU數(shù)目(1 073)和特有的OTU數(shù)目(722)(圖2: A);(2)對于真菌來說,CT的總OTU數(shù)目(821)和特有的OTU數(shù)目(738)亦均大于CR的總OTU數(shù)目(750)和特有的OTU數(shù)目(667)(圖2: B)。

    2.3 大別山五針?biāo)筛H菌和內(nèi)生菌α-多樣性分析

    由綜合分析細(xì)菌、真菌的α-多樣性指數(shù)可知,根際土壤樣品中的微生物在細(xì)菌和真菌方面的多樣性和物種豐富度均高于根部組織樣品,但差異不顯著。根際土壤樣品細(xì)菌的Observed-species指數(shù)、Chao指數(shù)、ACE指數(shù)分別為1 964.67、2 071.78、2 082.44,均大于根部組織樣品(1 068.67、1 209.03、1 230.83),但無顯著性差異(P>0.05);而Shannon指數(shù)為7.80,大于根部組織樣品(3.20),并且有顯著性差異(P<0.05),根際土壤樣品真菌的Observed-species指數(shù)、Chao指數(shù)、ACE指數(shù)和Shannon指數(shù)分別為381.67、406.95、406.49和4.82,亦均大于根部組織樣品的310.67、329.04、329.75和3.72,并且均無顯著性差異(P>0.05)(表1)。這說明大別山五針?biāo)筛H土壤真菌的物種豐富度和多樣性高于根部內(nèi)生真菌,但兩者無顯著性差異。

    2.4 大別山五針?biāo)筛H微生物和內(nèi)生菌β-多樣性分析

    PCoA分析(主坐標(biāo)分析),即通過樣品間的空間距離來反映樣品組間群落結(jié)構(gòu)的差異性。其中,組間的空間距離越遠(yuǎn),說明樣品間群落結(jié)構(gòu)差異越大(高嵩等,2021)。本文基于Unweichted Unitrac 距離進行PCoA分析,結(jié)果顯示無論是細(xì)菌還是真菌,大別山五針?biāo)筛H土壤微生物菌群和根部內(nèi)生菌菌群坐標(biāo)點的空間距離都較遠(yuǎn),說明兩組樣品間群落結(jié)構(gòu)存在差異。由圖3可知,根際土壤樣品與根部組織樣品同組間的距離較近,表明組間群落結(jié)構(gòu)組成相似,而根際土壤樣品與根部組織樣品不同組間距離較遠(yuǎn),表明兩組樣品的群落結(jié)構(gòu)相似性較低。

    2.5 大別山五針?biāo)筛H微生物和內(nèi)生菌群落結(jié)構(gòu)組成分析

    2.5.1 細(xì)菌 從大別山五針?biāo)筛H土壤樣品中獲得的細(xì)菌分布于174綱88門342目437科671屬,根部組織樣品中獲得的細(xì)菌分布于91綱40門195目255科318屬。在門分類水平上,厚壁菌門(Firmicutes)、unidentified-Bacteria和變形菌門(Proteobacteria)為兩組樣品的共有菌門。其中,前兩者在根際土壤樣品中相對豐度高于根部組織樣品(厚壁菌門: 31.31%gt;5.48%;unidentified-Bacteria: 17.70%gt;2.35%),而后者在根際土壤樣品中相對豐度低于根部組織樣品(變形菌門: 9.50%lt;23.23%)。藍(lán)藻菌門(Cyanobacteria)作為根部組織樣品特有的菌門,其相對豐度為61.82% (圖4)。在屬分類水平上,根際土壤樣品乳桿菌屬(Lactobacillus,0.73%)與根部組織樣品(2.60%)有顯著差異(P<0.05),并且相對豐度低于根部組織樣品。根際土壤樣品的優(yōu)勢菌屬為李斯特氏菌屬(Planococcus, 14.46%),其次為unidentified-Actinomarinales(5.71%)、芽孢桿菌屬(Bacillus, 3.90%)和擬桿菌屬(Bacteroides, 2.85%);而根部組織樣品的優(yōu)勢菌屬為unidentified-Chloroplast(61.46%),其次為unidentified-Mitochondria(14.10%)、乳桿菌屬(Lactobacillus,2.60%)和孿生菌屬(Macrococcus,2.44%)(表2)。

    2.5.2真菌 從大別山五針?biāo)筛H土壤樣品中獲得的真菌分布于11門187屬,從根部組織樣品中獲得的真菌分布于9門262屬。在門分類水平上,子囊菌門(Ascomycota)、球囊菌門(Glomeromycota)、壺菌門(Chytridiomycota)、羅茲菌門(Rozellomycota)和擔(dān)子菌門(Basidiomycota)為兩組樣品的共有菌門。其中,前四者在根際土壤樣品中相對豐度均高于根部組織樣品(子囊菌門: 50.56%>47.74%; 球囊菌門: 6.75%>0.18%;壺菌門: 5.99%>0.07%; 羅茲菌門: 5.48%>0.08%);而擔(dān)子菌門相對豐度明顯低于根部組織樣品(擔(dān)子菌門: 1.40%<40.69%)(圖4)。在屬分類水平上,籃狀菌屬(Talaromyces)、Aspergil和鐮刀菌屬(Fusarium)為兩組樣品所共有,相對豐度不同,并且不存在顯著性差異(P>0.05)。在屬分類水平上,根際土壤樣品中的優(yōu)勢菌屬為Melanocarpus和Pithoascus,其相對豐度分別為14.92%和3.22%;而根部組織樣品中的優(yōu)勢菌屬為unidentified_Helotiales_sp、須腹菌屬(Rhizopogon)和Thozetella,其相對豐度分別為25.67%、25.37%和10.00% (表2)。

    2.6 大別山五針?biāo)筛H微生物和內(nèi)生菌的功能預(yù)測

    本研究利用PICRUSt軟件分析大別山五針?biāo)傻母H細(xì)菌和根部內(nèi)生細(xì)菌功能預(yù)測,并比較其代謝功能通路差異。結(jié)果如圖5所示,除了不可分類通路(Unclassified, None)以外,兩組樣品在一級功能層上共有6個代謝通路,分別為新陳代謝(Metabolism)、遺傳信息處理(Genetic_Information_Processing)、環(huán)境信息處理(Environmental_Information_Processing)、細(xì)胞進程(Cellular_Processes)、人類疾?。℉uman_Diseases)和有機系統(tǒng)(Organismal_Systems)。其中,新陳代謝通路由兩組樣品的主要功能組成,在根際細(xì)菌的占比為50.25%,在根部內(nèi)生細(xì)菌的占比為49.66%。同時,對大別山五針?biāo)筛H細(xì)菌和根部內(nèi)生細(xì)菌基因二級功能層進行預(yù)測分析,結(jié)果發(fā)現(xiàn)兩組樣品的二級功能層主要由膜運輸(Membrane_Transport)、氨基酸代謝(Amino_Acid_Metabolism)、碳水化合物代謝(Carbohydrate_Metabolism)、能量代謝(Energy_Metabolism)和輔助因子和維生素的代謝(Metabolism_of_Cofactors_and_Vitamins)等41個子功能組成,并且功能豐富。對二級代謝通路預(yù)測基因相對豐度較高的子功能進行分析,結(jié)果表明這些子功能基因的相對豐度在兩組樣品中存在差異。其中,根際細(xì)菌在膜運輸、氨基酸代謝、碳水化合物代謝方面特征不明顯,脂質(zhì)代謝和異生物質(zhì)生物降解代謝方面高于根部內(nèi)生細(xì)菌。在能量代謝、復(fù)制和修復(fù)、輔因子和維生素的代謝和翻譯等子功能上,根際細(xì)菌樣品低于根部內(nèi)生細(xì)菌。細(xì)胞運動性、環(huán)境適應(yīng)性等功能主要富集于根際細(xì)菌群落,而根部內(nèi)生細(xì)菌功能主要體現(xiàn)在酶家族、消化系統(tǒng)等方面。

    基于FUNGuild對大別山五針?biāo)筛H真菌和根部內(nèi)生真菌的功能進行預(yù)測,并按照真菌營養(yǎng)模式進行分類。結(jié)果顯示CR和CT樣品中均富集到共生營養(yǎng)型、腐生營養(yǎng)型、病理營養(yǎng)型及多種交叉營養(yǎng)型的真菌類群,并且以前兩者為主要響應(yīng)功能群,具體包括未知腐生營養(yǎng)模式(Undefined_Saprotroph)、外生菌根(Ectomycorrhizal)、植物病原體土壤腐生物木材腐生物(Plant_Pathogen-Soil_Saprotroph-Wood_Saprotroph)、動物病原體內(nèi)生菌真菌寄生物植物病原體木材腐生物(Animal_Pathogen-Endophyte-Fungal_Parasite-Plant_Pathogen-Wood_Saprotroph)、動物病原體內(nèi)生菌植物病原體木材腐生物(Animal_Pathogen-Endophyte-Plant_Pathogen-Wood_Saprotroph)等(圖6: A)。其中,進一步對被分類的真菌營養(yǎng)模式進行差異性分析后發(fā)現(xiàn),僅外生菌根、動物病原體內(nèi)生菌真菌寄生物植物病原體木材腐生物、動物病原體內(nèi)生菌植物病原體木材腐生物3種營養(yǎng)模式的真菌在兩者樣品中存在顯著性差異(Plt;0.05)。其中,外生菌根在CR樣品中所占比例為25.99%,顯著大于其在CT樣品中所占豐度比例(0.44%) (P=0.048)(圖6: B)。

    3 討論與結(jié)論

    3.1 根際微生物與根部內(nèi)生菌群落組成分析

    Observed-species指數(shù)、Chao指數(shù)和ACE指數(shù)用于反映微生物的物種豐富度,而Shannon指數(shù)則用于反映物種的多樣性(武亞婷等, 2019)。本研究中,根際微生物的總OTUs數(shù)目、Observed-species指數(shù)、Chao指數(shù)和ACE指數(shù)均明顯大于根部內(nèi)生菌群。大量研究結(jié)果也證實,根際微生物的群落豐富度與多樣性均高于根部內(nèi)生菌。有關(guān)韭菜、黔中金蕎麥不同生態(tài)位群落組成的研究發(fā)現(xiàn),韭菜和金蕎麥根際土壤微生物的OTUs數(shù)目、Chao指數(shù)、Shannon指數(shù)、Simpson even指數(shù)均大于根部內(nèi)生菌(付靜等, 2018;張濤等, 2022;孫楠等, 2023)。肖健等(2024)的研究表明,甘蔗根際土壤樣品中真菌的Chao指數(shù)、Shannon指數(shù)均大于根系內(nèi)生真菌。

    本研究對于細(xì)菌來說, 厚壁菌門和變形菌門為兩個樣品的共有菌門,而與根部組織樣品相比,根際土壤樣品中還含有少量的脫硫桿菌門和擬桿菌門。其中,乳桿菌屬為兩組樣品的共有菌屬,根際土壤樣品主要菌屬為李斯特氏菌屬、芽孢桿菌屬、醋桿菌屬和擬桿菌屬,根部組織樣品主要菌屬為伯克氏菌屬、不動桿菌屬、腸桿菌屬和慢生根瘤菌屬。對于真菌來說,子囊菌門和擔(dān)子菌門為兩個樣品所共有,而與根部組織樣品相比,根際土壤樣品中還含有少量的壺菌門和羅茲菌門。其中,兩組樣品的共有菌屬為鐮刀菌屬和籃狀菌屬,根際土壤樣品的主要菌屬為毛殼屬和青霉屬,根部組織樣品的主要菌屬為須腹菌屬和針菌屬。

    張秋玉等(2022)研究發(fā)現(xiàn)莼菜內(nèi)生細(xì)菌的優(yōu)勢菌群為變形菌門、厚壁菌門、放線菌門和擬桿菌門,安超等(2021)研究的野生竹根七根際土的細(xì)菌主要分布在變形菌門和放線菌門,以及迪拉熱·海米提等(2021)在鹽穗木的根際土壤細(xì)菌中發(fā)現(xiàn)優(yōu)勢菌門為放線菌門和變形菌門,而內(nèi)生真菌的優(yōu)勢門多為子囊菌門和擔(dān)子菌門,這些結(jié)果與本研究的相似。這說明不同植物根際微生物或內(nèi)生菌既存在差異,又具有一定的共性,尤其是變形菌門。Tang等(2020)研究表明,變形菌廣泛存在于土壤環(huán)境中,適應(yīng)能力強,在生理、形態(tài)和代謝等方面呈現(xiàn)出豐富的多樣性,對土壤中的碳氮循環(huán)及生態(tài)系統(tǒng)具有重要的影響。石水琴等(2022)對大別山五針?biāo)筛H細(xì)菌的研究結(jié)果與本研究的有所不同,前者主要種類為變形菌門、放線菌門和酸桿菌門,推測其原因可能是生境不同導(dǎo)致的差異。前者取樣地點為大別山五針?biāo)扇斯ち?,而本研究為天然林,植物的生境存在較大差別。

    3.2 根際微生物與根部內(nèi)生菌功能預(yù)測分析

    據(jù)前人研究發(fā)現(xiàn),根部微生物對植物的生長有很大的影響,如慢生根瘤菌屬、腸桿菌屬、伯克霍爾德氏菌屬、擬桿菌屬、不動桿菌屬和芽孢桿菌屬等都被證實可促進植物生長(張典利等,2018)。朱滕滕等(2023)研究發(fā)現(xiàn)芽孢桿菌屬、不動桿菌屬的菌株對硒砂瓜幼苗的多種生長指標(biāo)均有促進作用。江美彥等(2022)研究發(fā)現(xiàn)分離得到的芽孢桿菌屬的菌株對白芷的生長具有較強的促進作用。本研究在兩組樣品的PICRUSt功能預(yù)測中,氨基酸代謝、碳水化合物代謝及輔助因子和維生素代謝這3種子功能代謝占比較大。氨基酸是構(gòu)成蛋白質(zhì)的基本單位,游離的氨基酸可以通過調(diào)節(jié)滲透壓幫助植物抵御環(huán)境的脅迫(燕輝等,2012)。碳水化合物是植物生長的重要營養(yǎng)物質(zhì),林夏珍等(2021)研究表明,非結(jié)構(gòu)性碳水化合物在一定程度上可以增加植物的抗逆性。維生素是普遍存在于植物中的重要營養(yǎng)素,安華明等(2004)研究發(fā)現(xiàn)植物體內(nèi)合成的維生素C在抗氧化和自由基清除、光合作用和光保護、細(xì)胞生長和分裂以及一些重要次生代謝物和乙烯的合成等方面具有非常重要的生理功能。由此推測,大別山五針?biāo)筛H土壤樣品與根部組織樣品中不同細(xì)菌組成不同功能的基因,從而影響到氨基酸、碳水化合物及輔助因子和維生素等方面的代謝,對促進大別山五針?biāo)傻纳L,提高其對環(huán)境的抵抗性等方面發(fā)揮了重要作用。

    陳羽彤等(2024)研究顯示,土壤真菌群落的變化可以增加生態(tài)系統(tǒng)穩(wěn)定性,與宿主植物的生長發(fā)育密切相關(guān)。本研究發(fā)現(xiàn),大別山五針?biāo)筛恐邪采鸂I養(yǎng)型、腐生營養(yǎng)型、病理營養(yǎng)型及多種交叉營養(yǎng)型等四大真菌類群,并且以共生營養(yǎng)型和腐生營養(yǎng)型為主要響應(yīng)功能群,這與馬源等(2024)的研究結(jié)果一致。一般來說,腐生營養(yǎng)型真菌大多歸屬為子囊菌門,在分解土壤有機質(zhì)和營養(yǎng)物質(zhì)循環(huán)方面起重要的作用。菌根真菌是共生營養(yǎng)型的重要代表,其與植物存在共生關(guān)系,菌根真菌與寄主植物相互作用有助于植物吸收養(yǎng)分,菌根真菌也可以從寄主植物根中提取脂質(zhì)和碳水化合物。外生菌根真菌是森林生態(tài)系統(tǒng)中與樹木共生的重要成員,它與植物跨界共生具有復(fù)雜的分子互作過程(俞嘉瑞和袁海生,2023)。本研究中,大別山五針?biāo)筛績?nèi)生真菌中的外生菌根類群顯著高于土壤根際真菌類群,說明外生菌根在促進大別山五針?biāo)缮L發(fā)育方面起到十分重要的作用。有關(guān)外生菌根的生理活動對其宿主具體的影響機制有待后續(xù)進一步探索。另外,本研究仍有33.88%~68.31%的根部真菌功能未被鑒定,說明真菌群落功能的復(fù)雜性仍需深入研究。

    本研究采用高通量測序分析了大別山五針?biāo)筛课⑸锏亩鄻有院腿郝浣Y(jié)構(gòu),并預(yù)測了其功能。本研究結(jié)果顯示,根際土壤樣品微生物的多樣性和物種豐富度在細(xì)菌和真菌方面均高于根部組織樣品,并且兩組樣品在細(xì)菌和真菌方面的群落結(jié)構(gòu)組成差異較大。本研究功能注釋結(jié)果暗示,特定的細(xì)菌(不動桿菌屬和芽孢桿菌屬等)和真菌(外生菌根真菌等)可能在大別山五針?biāo)傻纳L與環(huán)境適應(yīng)中發(fā)揮了重要作用。然而,本研究僅基于PICRUSt和FUNGuild對大別山五針?biāo)傻募?xì)菌和真菌功能進行了初步探索,各類群(尤其是關(guān)鍵類群)的確切功能尚不明晰,后續(xù)將利用宏基因組測序技術(shù)深入探索大別山五針?biāo)傻娜郝浣Y(jié)構(gòu)組成及其相關(guān)的功能注釋。此外,后續(xù)研究還將在本研究基礎(chǔ)上通過改良培養(yǎng)基類型、模擬大別山自然生態(tài)條件,盡可能富集到多樣的可培養(yǎng)微生物分離株,并從相關(guān)的菌種庫中分離篩選具有抗逆促生功效的有益菌株,為促進大別山五針?biāo)傻纳L發(fā)育提供參考依據(jù)。

    參考文獻:

    AFZAL I, SHINWARI ZK, SIKANDAR S, et al., 2019. Plant beneficial endophytic bacteria: Mechanisms, diversity, host range and genetic determinants "[J]. Microbiological Research, 221: 36-49.

    ALTSCHUL SF, GISH W, MILLER W, et al., 1990. Basic local alignment search tool "[J]. Journal Molecular Biology, 215(3): 403-410.

    AN C, MA SJ, XUE WJ, et al., 2021. Analysis of endophytic and rhizosphere bacterial community characteristics of Tupistra chinensis Baker based on high-throughput sequencing technology "[J]. Genomics and Applied Biology, 40(5): 2076-2082. "[安超, 馬賽箭, 薛文嬌, 等, 2021. 基于高通量測序的野生竹根七根際及內(nèi)生細(xì)菌群落結(jié)構(gòu)特征分析 [J]. 基因組學(xué)與應(yīng)用生物學(xué), 40(5): 2076-2082.]

    AN HM, CHEN LG, FAN WG, et al., 2004. Advances in research on function, biosynthesis and metabolism of ascorbic acid in higher plants "[J]. Chinese Bulletin of Botany, 21(5): 608-617. "[安華明, 陳力耕, 樊衛(wèi)國, 等, 2004. 高等植物中維生素C的功能、合成及代謝研究進展 [J]. 植物學(xué)通報, 21(5): 608-617.]

    BOKULICH NA, SUBRAMANIAN S, FAITH JJ, et al., 2013. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing "[J]. Nature Methods, 10(1): 57-59.

    CAPORASO JG, KUCZYNKI J, STOMBAUGH J, et al., 2010. QIIME allows analysis of high-throughput community sequencing data "[J]. Nature Methods, 7(5): 335-336.

    CHEN YT, ZHANG LM, WANG SQ, et al., 2024. Soil fungal community structure and functional group in response to different soil organic carbon inputs in the temperate forest during the freeze-thaw season "[J]. Acta Ecologica Sinica, 44(6): 2244-2255. "[陳羽彤, 張利敏,王思琪, 等, 2024. 凍融季溫帶森林土壤真菌群落結(jié)構(gòu)及功能類群對不同土壤有機碳輸入的響應(yīng) [J]. 生態(tài)學(xué)報, 44(6): 2244-2255. ]

    COMPANT S, CLMENT C, SESSITSCH A, 2010. Plant growth-promoting bacteria in the rhizo-"and endosphere of plants: Their role, colonization, mechanisms involved and prospects for utilization "[J]. Soil Biology and "Biochemistry, 42: 669-678.

    DILARE HMT, FAN YH, WANG WN, et al., 2021. Analysis of microbial communities in leaves and rhizosphere soil of Halostachys capsica by high-throughput sequencing "[J]. Xinjiang Agricultural Sciences, 58(4): 731-740. "[迪拉熱·海米提, 樊永紅, 王偉楠, 等. 2021. 鹽穗木葉片及根際土壤微生物群落高通量分析 [J]. 新疆農(nóng)業(yè)科學(xué), 58(4): 731-740.]

    EDGAR RC, 2013. UPARSE: highly accurate OTU sequences from microbial amplicon reads "[J]. Nature Methods, 10(10): 996-998.

    FU J, LIU YF, MEI M, et al., 2018. Comparing analysis of diversity of Chinese leek rhizospheritic and endophytic fungal communities "[J]. Acta Agriculturae Boreali-Sinica, 33(z1):246-252. "[付靜, 劉堰鳳, 梅眉, 等, 2018. 韭菜根際真菌和內(nèi)生真菌種群多樣性分析 [J]. 華北農(nóng)學(xué)報, 33(增刊): 246-252.]

    GAO S, SUN WS, YU CL, et al., 2021. Analysis of fungal community structure in gentian rhizosphere soil based on high-throughput sequencing "[J]. Jiangsu Agricultural Sciences, 49(12): 190-195. "[高嵩, 孫文松, 于春雷, 等, 2021. 基于高通量測序的龍膽草根際土壤真菌群落結(jié)構(gòu)分析 [J]. 江蘇農(nóng)業(yè)科學(xué), 49(12): 190-195.]

    GONG Q, GUAN DM, WANG YB, et al., 2005. Meth extracted Escherichia coli genomic DNA using the fast ways "[J]. Marine Environmental Science, 24(4): 63-66. "[宮強, 關(guān)道明, 王耀兵, 等, 2005. 大腸桿菌總DNA快速提取方法的比較研究 [J]. 海洋環(huán)境科學(xué), 24(4): 63-66.]

    GOPAL M, GUPTA A, 2016. Microbiome selection could spur next-generation plant breeding strategies "[J]. Frontiers Microbiology, 7: 1971.

    HAN JW, ZHANG ZY, WANG EM, et al., 2014. Pinus dabeshanensis seed characteristics and the promotion of seed germination "[J]. Chinese Agricultural Science Bulletin, 30(1): 5-10. "[韓建偉, 張智勇, 王恩茂, 等, 2014. 大別山五針?biāo)煞N子特性及促進種子萌發(fā)的研究 [J]. 中國農(nóng)學(xué)通報, 30(1): 5-10.]

    HUANG XZ, ZHAO LF, 2023. Mechanism of endophytes of medicinal plants in promoting the growth of host plants "[J]. Microbiology China, 50(4):1653-1665. "[黃雪珍, 趙龍飛, 2023. 藥用植物內(nèi)生菌對宿主植物促生作用機制研究進展 [J]. 微生物學(xué)通報, 50(4):1653-1665.]

    JIANG MY, ZHOU Y, LIU RL, et al., 2022. Screening and plant growth promoting of grow-promoting bacteria in rhizosphere bacteria of Angelica dahurica var. formosana "[J]. Biotechnology Bulletin, 38(8): 167-178. "[江美彥, 周楊, 劉仁浪, 等, 2022. 白芷根際促生菌的篩選及其促生效果研究 [J]. 生物技術(shù)通報, 38(8): 167-178.]

    KLJALG U, NILSSON RH, ABARENKOV K, et al., 2013. Towards a unified paradigm for sequence-based identification of fungi "[J]. Molecular Ecology, 22(21): 5271-5277.

    LIN XZ, LIN L, DONG TT, et al., 2021. Effects of non-structural carbohydrate and nitrogen allocation on the ability of Populus deltoides and P. cathayana to resist soil salinity stress "[J]. Chinese Journal of "Plant Ecology, 45(9): 961-971. [林夏珍, 林劉, 董婷婷, 等, 2021. 非結(jié)構(gòu)性碳水化合物與氮分配對美洲黑楊和青楊耐鹽能力的影響 [J]. 植物生態(tài)學(xué)報, 45(9): 961-971.]

    LIU LH, JIANG HM, QU YC, et al., 2020. Identification and growth promotion of endophytic bacteria isolated from Oryza meridionalis "[J]. Chinese Journal of Applied Environmental Biology, 26(5): 1051-1058. "[劉麗輝, 蔣慧敏, 區(qū)宇程, 等, 2020. 南方野生稻內(nèi)生細(xì)菌的分離鑒定及促生作用 [J]. 應(yīng)用與環(huán)境生物學(xué)報, 26(5): 1051-1058.]

    MA Y, WANG XL, MA YS, et al., 2024. Effects of the degree of alpine meadow degradation on the rhizosphere soil fungal community and the ecological network of dominant species "[J]. Acta Prataculturae Sinica, 33(2): 125-137. "[馬源, 王曉麗, 馬玉壽, 等, 2024. 高寒草甸退化程度對優(yōu)勢物種根際土壤真菌群落和生態(tài)網(wǎng)絡(luò)的影響 [J]. 草業(yè)學(xué)報, 33(2): 125-137.]

    MCINROY JA, KLOEPPER WJ, 1994. Studies on indigenous endophytic bacteria of sweet corn and cotton "[J]. Molecular Ecology of "Rhizosphere Microorganisms.https://doi.org/10.1002/9783527615810.ch2.

    PREZ-JARAMILLO JE,DE HOLLANDER M, RAMREZ CA, et al., 2019. Deciphering rhizosphere microbiome assembly of wild and modern common bean (Phaseolus vulgaris) in native and agricultural soils from Colombia "[J]. Microbiome, 7: 1-16.

    QIN YM, ZHOU XL, GUAN QL, et al., 2021. Fungal diversity in natural tomato fermentation as shown by high-throughput sequencing "[J]. Fujian Journal of "Agricultural Sciences, 36(9): 1110-1118. "[秦宇蒙, 周笑犁, 管慶林, 等. 2021. 基于高通量測序分析番茄自然發(fā)酵過程中的真菌多樣性 [J]. 福建農(nóng)業(yè)學(xué)報, 36(9): 1110-1118.]

    REINHOLD-HUREK B, BNGER W, BURBANO CS, et al., 2015. Roots shaping their microbiome: global hotspots for microbial activity "[J]. Annual Review of "Phytopathology, 53(1): 403-424.

    ROGNES T, FLOURI T, NICHOLS B, et al., 2016. VSEARCH: a versatile open source tool for metagenomics "[J]. PeerJ, 4: e2584.

    SHI JY, CHEN WX, LIU AY, 2006. Advances in the study of endophytes and their effects on control of plant diseases "[J]. Acta Ecologica Sinica, 26(7): 2395-2401. "[石晶盈, 陳維信, 劉愛媛, 2006. 植物內(nèi)生菌及其防治植物病害的研究進展 [J]. 生態(tài)學(xué)報, 26(7): 2395-2401.]

    SHI SQ, QIN HG, ZHANG JJ, et al., 2022. Characteristics and function analysis of rhizosphere bacterial community of endangered plant Pinus dabeshanensis "[J]. Chinese Bulletin of "Botany, 57(4): 457-467. "[石水琴, 秦華光, 張靜靜, 等, 2022. 瀕危植物大別山五針?biāo)筛H細(xì)菌群落特征與功能分析 [J]. 植物學(xué)報, 57(4): 457-467.]

    SUN N, LIANG FF, CHEN JH, et al., 2023. Effects of different niches of Chinese chives on the composition of microbial community and functional characters [J]. Journal of "China Agricultural University, 28(3): 61-73. "[孫楠, 梁方方, 陳建華, 等, 2023. 韭菜不同生態(tài)位微生物群落組成及功能特征 [J]. 中國農(nóng)業(yè)大學(xué)學(xué)報, 28(3): 61-73.]

    TAN Y, CUI Y, LI H, et al., 2017. Rhizospheric soil and root endogenous fungal diversity and composition in response to continuous Panax notoginseng cropping practices [J]. Microbiological Research, 194: 10-19.

    TANG HM, XIAO XP, LIC, et al., 2020. Microbial carbon source utilization in rice rhizosphere and nonrhizosphere soils with short-term manure N input rate in paddy field [J]. Scientific Reports, 10:6487.

    WANG LH, HUANG QF, PU FG, et al., 2014. Interspecific relationship between Pinus dabeshanensis and dominant species in the community at Tianma National Nature Reserve "[J]. Resources and Environment in the Yangtze Basin, 23(7): 918-922. "[王雷宏, 黃慶豐, 蒲發(fā)光, 等, 2014. 天馬自然保護區(qū)大別山五針?biāo)膳c群落中優(yōu)勢種的種間關(guān)系 [J]. 長江流域資源與環(huán)境, 23(7): 918-922.]

    WANG Q, GARRITY GM, TIEDJE JM, et al., 2007. Nave bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy "[J]. Applied and Environmental Microbiology, 73(16): 5261-5267.

    WANG XS, LIAN YH, GUO H, et al., 2022. Effects of wheat/safflower intercropping on rhizosphere microbial community function and structure "[J]. Chinese Journal of Eco-Agriculture, 31(4): 516-529. "[王香生, 連延浩, 郭輝, 等, 2022. 小麥紅花間作系統(tǒng)下根際微生物群落結(jié)構(gòu)及功能分析 [J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(中英文), 31(4): 516-529.]

    WU GL, XIANG XY, DUAN RY, et al., 2016. Research progress and conservation measures of Pinus dabeshanensis "[J]. Journal of "Anqing Teachers College (Natural Science Edition), 22(4): 97-99. "[吳甘霖, 項小燕, 段仁燕, 等, 2016. 大別山五針?biāo)傻难芯窟M展及保護對策 [J]. 安慶師范學(xué)院學(xué)報(自然科學(xué)版), 22(4): 97-99.]

    WU YT, DU MY, HE HH, et al., 2019. Microbial diversity analysis of natural fermented Chili Sauce from different regions in Xinjiang by high-throughput sequencing "[J]. Food and Fermentation Industries, 45(21): 221-228. "[武亞婷, 杜木英, 何歡歡, 等, 2019. 基于高通量測序技術(shù)分析新疆不同地區(qū)自然發(fā)酵辣椒醬微生物群落多樣性 [J]. 食品與發(fā)酵工業(yè), 45(21): 221-228.]

    XIANG XY, WU GL, DUAN RY, et al., 2015. Intraspecific and interspecific competition of Pinus dabeshanesis "[J]. Acta Ecologica Sinica, 35(2): 389-395. "[項小燕, 吳甘霖, 段仁燕, 等, 2015. 大別山五針?biāo)煞N內(nèi)和種間競爭強度 [J]. 生態(tài)學(xué)報, 35(2): 389-395.]

    XIANG XY, WU GL, DUAN RY, et al., 2016a. Studies on population structure and dynamics of Pinus dabeshanensis "[J]. Resources and Environment in the Yangtze Basin, 25(1): 55-62. "[項小燕, 吳甘霖, 段仁燕, 等, 2016a. 大別山五針?biāo)煞N群結(jié)構(gòu)及動態(tài)研究 [J]. 長江流域資源與環(huán)境, 25(1): 55-62.]

    XIANG XY, WU GL, DUAN RY, et al., 2016b. Change in pollen viability of Pinus dabeshanensis at different storage temperatures "[J]. Journal of Plant Resources Environment, 25(2): 114-116. "[項小燕, 吳甘霖, 段仁燕, 等. 2016b. 不同貯藏溫度下大別山五針?biāo)苫ǚ刍盍Φ淖兓?[J]. 植物資源與環(huán)境學(xué)報, 25(2): 114-116.]

    XIAO J, YANG SD, TIAN L, et al., 2024. Effects of long-term applications with nitrogen fertilizers on soil and endophytic fungal compositions in rhizospheres and roots of sugarcanes "[J]. Chinese Journal of Tropical Crops, 45(5): 1040-1051. "[肖健, 楊尚東, 闐梁, 等, 2024. 長期施氮對甘蔗根際土壤真菌和根系內(nèi)生真菌群落結(jié)構(gòu)的影響 [J]. 熱帶作物學(xué)報, 45(5): 1040-1051.]

    XIE J, XU W, ZHANG M, et al., 2021. The impact of the endophytic bacterial community on mulberry tree growth in the Three Gorges Reservoir ecosystem, China [J]. Environmental Microbiology, 23(4): 1858-1875.

    YAN H, PENG XB, XUE JJ, 2012. Effects of NaCl stress on leaf photosynthesis characteristics and free amino acid metabolism of Heyedysarum scoparium "[J]. Chinese Journal of "Applied Ecology, 23(7): 1790-1796. "[燕輝, 彭曉邦, 薛建杰, 2012. NaCl脅迫對花棒葉片光合特性及游離氨基酸代謝的影響 [J]. 應(yīng)用生態(tài)學(xué)報, 23(7): 1790-1796.]

    YU JR, YUAN HS, 2023. Research progress on symbiotic interaction and host selection mechanisms of ectomycorrhizal fungi "[J]. Mycosystema, 42(1): 86-100. [俞嘉瑞, 袁海生, 2023. 外生菌根真菌的共生互作和宿主選擇機制研究進展 [J]. 菌物學(xué)報, 42(1): 86-100.]

    ZHANG QY, CHEN J, ZHOU XM, et al., 2022. Diversity of entophytic bacteria in different tissues of Brasenia schreberi "[J]. Journal of South-Central Minzu University (Natural Science Edition), 41(3): 285-291. "[張秋玉, 陳晶, 周雪妹, 等, 2022. 莼菜不同組織內(nèi)生細(xì)菌多樣性的比較分析 [J]. 中南民族大學(xué)學(xué)報(自然科學(xué)版), 41(3): 285-291.]

    ZHANG T, ZHOU SX, TANG YJ, et al., 2022. Diversity analysis of Fagopyrum dibotrys Qian Zhong endophytes and rhizosphere soil microorganism by high-throughput sequencing "[J]. Shandong Agricultural Science, 54(11): 70-75. "[張濤, 周思旋, 唐遠(yuǎn)江, 等, 2022. 基于高通量測序分析黔中金蕎麥內(nèi)生菌及其根際土壤微生物多樣性 [J]. 山東農(nóng)業(yè)科學(xué), 54(11): 70-75.]

    ZHAO J, LIU T, PAN L, et al., 2015. Isolation and identification of root endophytic and rhizosphere bacteria of rice landraces in Yuanyang Terrace, China "[J]. Chinese Journal of Applied Ecology, 26(12): 3737-3745. "[趙娟, 劉濤, 潘磊, 等, 2015. 元陽梯田地方水稻品種根部內(nèi)生菌及根際微生物的分離與鑒定 [J]. 應(yīng)用生態(tài)學(xué)報, 26(12): 3737-3745.]

    ZHANG DL, MENG Z, QI WZ, et al., 2018. The research and application status of plant growth promoting rhizobacteria "[J]. World Pesticides, 40(6): 37-43. "[張典利, 孟臻, 齊文哲, 等, 2018. 植物根際促生菌的研究與應(yīng)用現(xiàn)狀 [J]. 世界農(nóng)藥, 40(6): 37-43.]

    ZHU TT, HE XL, MIAN XJ, et al., 2023. Growth-promoting bacteria screening from Selenium watermelon rhizosphere and their growth promoting effects "[J]. Agricultural Research in the Arid Areas, 41(2): 221-229. [朱滕滕, 何學(xué)蓮, 勉小娟, 等, 2023. 硒砂瓜根際促生菌篩選鑒定及其促生效應(yīng)研究 [J]. 干旱地區(qū)農(nóng)業(yè)研究, 41(2): 221-229.]

    LOU J, LIU Y, YAN L, 2014. Review of high-throughput sequencing techniques in studies of soil microbial diversity "[J]. Chinese Agricultural Science Bulletin, 30(15): 256-260. "[樓駿, 柳勇, 李延, 2014. 高通量測序技術(shù)在土壤微生物多樣性研究中的研究進展 [J]. 中國農(nóng)學(xué)通報, 30(15): 256-260.]

    (責(zé)任編輯 蔣巧媛 王登惠)

    猜你喜歡
    高通量測序多樣性
    環(huán)狀RNA在疾病發(fā)生中的作用
    川明參輪作對煙地土壤微生物群落結(jié)構(gòu)的影響
    多穗柯轉(zhuǎn)錄組分析及黃酮類化合物合成相關(guān)基因的挖掘
    人參根際真菌群落多樣性及組成的變化
    LncRNAs作為miRNA的靶模擬物調(diào)節(jié)miRNA
    由古典戲曲看“代言體”在中國的前世今生
    戲劇之家(2016年22期)2016-11-30 15:13:39
    淺談新時期群文輔導(dǎo)工作的特征
    新時期群文輔導(dǎo)工作的特征
    海洋微生物次生代謝的生物合成機制
    科技資訊(2016年19期)2016-11-15 10:39:12
    舞蹈表演的表現(xiàn)形式多樣性研究
    人間(2016年27期)2016-11-11 16:27:23
    亚洲 国产 在线| 91麻豆av在线| 国产欧美日韩综合在线一区二区| 51午夜福利影视在线观看| 老司机影院毛片| 老汉色av国产亚洲站长工具| 啦啦啦 在线观看视频| 亚洲熟妇熟女久久| 久久av网站| kizo精华| 一个人免费看片子| 丁香欧美五月| h视频一区二区三区| 久久精品国产亚洲av香蕉五月 | 90打野战视频偷拍视频| 1024视频免费在线观看| 最新的欧美精品一区二区| 亚洲中文av在线| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美在线一区| 曰老女人黄片| 免费少妇av软件| 久久中文字幕人妻熟女| 日韩一卡2卡3卡4卡2021年| 50天的宝宝边吃奶边哭怎么回事| 飞空精品影院首页| 午夜福利免费观看在线| 成人国语在线视频| 欧美激情久久久久久爽电影 | 大陆偷拍与自拍| 欧美日韩一级在线毛片| 精品卡一卡二卡四卡免费| 日韩免费av在线播放| 99re6热这里在线精品视频| 丝袜人妻中文字幕| 亚洲情色 制服丝袜| 曰老女人黄片| 久久天堂一区二区三区四区| 1024视频免费在线观看| 正在播放国产对白刺激| 少妇 在线观看| 无遮挡黄片免费观看| 欧美变态另类bdsm刘玥| 国产高清激情床上av| 国产成人免费无遮挡视频| 欧美 亚洲 国产 日韩一| 在线 av 中文字幕| 天堂俺去俺来也www色官网| 男人舔女人的私密视频| 18禁观看日本| 丰满饥渴人妻一区二区三| 麻豆成人av在线观看| 午夜视频精品福利| 少妇精品久久久久久久| 激情视频va一区二区三区| 久久久精品区二区三区| 少妇被粗大的猛进出69影院| 日韩欧美三级三区| 免费女性裸体啪啪无遮挡网站| 一级,二级,三级黄色视频| 亚洲熟女精品中文字幕| 男女下面插进去视频免费观看| 亚洲一码二码三码区别大吗| a级毛片在线看网站| 国产欧美日韩综合在线一区二区| 欧美一级毛片孕妇| www.熟女人妻精品国产| 丝袜美足系列| 欧美精品亚洲一区二区| 中文字幕制服av| 精品一区二区三卡| 91精品三级在线观看| 黄频高清免费视频| 男女边摸边吃奶| 国产精品美女特级片免费视频播放器 | 国产在线一区二区三区精| 欧美激情 高清一区二区三区| 嫁个100分男人电影在线观看| 搡老乐熟女国产| 午夜两性在线视频| 国产三级黄色录像| 国产成人一区二区三区免费视频网站| 不卡av一区二区三区| 国产xxxxx性猛交| 国产成人精品久久二区二区91| 欧美日韩亚洲国产一区二区在线观看 | 高清在线国产一区| 啦啦啦在线免费观看视频4| 精品国产超薄肉色丝袜足j| 国产一卡二卡三卡精品| bbb黄色大片| 熟女少妇亚洲综合色aaa.| 在线观看免费高清a一片| 亚洲专区国产一区二区| 国产在线一区二区三区精| 亚洲九九香蕉| 狂野欧美激情性xxxx| 在线 av 中文字幕| 久久久久精品人妻al黑| 一本色道久久久久久精品综合| 免费在线观看视频国产中文字幕亚洲| 80岁老熟妇乱子伦牲交| 桃红色精品国产亚洲av| 亚洲国产欧美网| 精品亚洲成a人片在线观看| 国产av精品麻豆| 可以免费在线观看a视频的电影网站| 不卡av一区二区三区| 在线观看免费视频网站a站| 一级,二级,三级黄色视频| 亚洲成国产人片在线观看| 女性被躁到高潮视频| av欧美777| 免费看十八禁软件| 欧美大码av| 国产无遮挡羞羞视频在线观看| 日本黄色日本黄色录像| 可以免费在线观看a视频的电影网站| 欧美性长视频在线观看| 美女高潮喷水抽搐中文字幕| 国产一区二区在线观看av| 一本久久精品| 久久毛片免费看一区二区三区| 男女边摸边吃奶| 国产成人av激情在线播放| 在线亚洲精品国产二区图片欧美| 国产有黄有色有爽视频| 丰满人妻熟妇乱又伦精品不卡| 天堂俺去俺来也www色官网| 中文字幕人妻丝袜一区二区| 亚洲熟女精品中文字幕| 精品国内亚洲2022精品成人 | 欧美+亚洲+日韩+国产| 国产精品av久久久久免费| 成年女人毛片免费观看观看9 | 国产无遮挡羞羞视频在线观看| 国产在线免费精品| 老司机福利观看| av免费在线观看网站| 国产伦理片在线播放av一区| 亚洲国产av新网站| 免费日韩欧美在线观看| 男女床上黄色一级片免费看| 国产97色在线日韩免费| 建设人人有责人人尽责人人享有的| 国产一区二区三区综合在线观看| 老司机午夜福利在线观看视频 | 国产有黄有色有爽视频| 精品第一国产精品| 一个人免费看片子| 久久精品亚洲av国产电影网| 又黄又粗又硬又大视频| 亚洲人成电影免费在线| 在线av久久热| 国产老妇伦熟女老妇高清| 久久久精品免费免费高清| 日韩大码丰满熟妇| 91精品国产国语对白视频| 亚洲美女黄片视频| 亚洲精华国产精华精| 久久影院123| 纯流量卡能插随身wifi吗| 少妇裸体淫交视频免费看高清 | 亚洲国产精品一区二区三区在线| 国产高清国产精品国产三级| 免费看a级黄色片| 欧美精品人与动牲交sv欧美| 热re99久久国产66热| 国产精品久久久av美女十八| 日本av免费视频播放| 王馨瑶露胸无遮挡在线观看| 久久久精品免费免费高清| 麻豆国产av国片精品| 女人爽到高潮嗷嗷叫在线视频| 日韩成人在线观看一区二区三区| 久久精品亚洲av国产电影网| 免费日韩欧美在线观看| 国产一区二区激情短视频| 国产精品亚洲av一区麻豆| 在线观看66精品国产| 亚洲精品美女久久久久99蜜臀| 亚洲人成伊人成综合网2020| 午夜两性在线视频| 亚洲av第一区精品v没综合| 久久午夜综合久久蜜桃| 少妇精品久久久久久久| 最黄视频免费看| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利影视在线免费观看| 欧美黑人欧美精品刺激| 国产男女内射视频| 黑人巨大精品欧美一区二区蜜桃| 99riav亚洲国产免费| 亚洲av片天天在线观看| 黄片小视频在线播放| 性少妇av在线| 亚洲精品国产区一区二| 天天操日日干夜夜撸| 黄色a级毛片大全视频| 国产日韩欧美亚洲二区| 免费看十八禁软件| 99国产精品一区二区三区| 国产精品免费一区二区三区在线 | 日韩大片免费观看网站| 国产精品1区2区在线观看. | 欧美日韩亚洲综合一区二区三区_| 国产男女内射视频| www.自偷自拍.com| 久久精品亚洲av国产电影网| 下体分泌物呈黄色| 成在线人永久免费视频| 精品国内亚洲2022精品成人 | 国产又色又爽无遮挡免费看| 极品人妻少妇av视频| 侵犯人妻中文字幕一二三四区| 性高湖久久久久久久久免费观看| 777米奇影视久久| 亚洲精品成人av观看孕妇| 成人亚洲精品一区在线观看| 欧美久久黑人一区二区| 欧美在线黄色| 国产精品美女特级片免费视频播放器 | 巨乳人妻的诱惑在线观看| 这个男人来自地球电影免费观看| 亚洲av片天天在线观看| 在线观看免费视频日本深夜| 黄频高清免费视频| 亚洲国产欧美网| 蜜桃在线观看..| 日本vs欧美在线观看视频| 一级黄色大片毛片| 日韩视频一区二区在线观看| 天堂8中文在线网| 免费黄频网站在线观看国产| 国产欧美日韩综合在线一区二区| 成人特级黄色片久久久久久久 | 精品欧美一区二区三区在线| 欧美日韩av久久| 久久 成人 亚洲| 一本大道久久a久久精品| 伊人久久大香线蕉亚洲五| 国产不卡av网站在线观看| 大陆偷拍与自拍| 一区二区av电影网| 啦啦啦在线免费观看视频4| 亚洲欧美激情在线| 一区二区日韩欧美中文字幕| 99精国产麻豆久久婷婷| 成人黄色视频免费在线看| 国产精品一区二区免费欧美| 黑人巨大精品欧美一区二区mp4| 亚洲国产看品久久| 男女免费视频国产| 午夜激情久久久久久久| 成年版毛片免费区| 18禁国产床啪视频网站| 一边摸一边抽搐一进一小说 | 国产精品免费大片| 女同久久另类99精品国产91| 国产精品免费一区二区三区在线 | 啦啦啦 在线观看视频| 一区在线观看完整版| 久久午夜综合久久蜜桃| 啦啦啦免费观看视频1| 69精品国产乱码久久久| 麻豆乱淫一区二区| 欧美精品一区二区大全| 中文字幕精品免费在线观看视频| 人人妻人人澡人人爽人人夜夜| 在线播放国产精品三级| 欧美在线黄色| 国产69精品久久久久777片 | 国产淫片久久久久久久久 | 国产精品女同一区二区软件 | 欧美性猛交黑人性爽| 国产成人aa在线观看| 国产亚洲精品综合一区在线观看| 国产高清视频在线观看网站| 国产精品亚洲av一区麻豆| 色哟哟哟哟哟哟| 精品国产乱码久久久久久男人| 啦啦啦观看免费观看视频高清| 草草在线视频免费看| 午夜激情欧美在线| 大型黄色视频在线免费观看| 99精品久久久久人妻精品| 在线a可以看的网站| 国产精品精品国产色婷婷| 精品一区二区三区av网在线观看| 1024香蕉在线观看| x7x7x7水蜜桃| 午夜福利高清视频| 毛片女人毛片| 欧美三级亚洲精品| 俄罗斯特黄特色一大片| 日本三级黄在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精品日韩av在线免费观看| 在线播放国产精品三级| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 熟妇人妻久久中文字幕3abv| 久久中文字幕人妻熟女| 亚洲国产欧美网| 亚洲熟妇熟女久久| 精品国产乱子伦一区二区三区| 久久精品国产亚洲av香蕉五月| 欧美丝袜亚洲另类 | 很黄的视频免费| 激情在线观看视频在线高清| www.999成人在线观看| 亚洲av电影在线进入| av片东京热男人的天堂| 色综合欧美亚洲国产小说| or卡值多少钱| 久久草成人影院| 国产免费男女视频| 色综合亚洲欧美另类图片| 午夜福利在线观看吧| 黑人欧美特级aaaaaa片| 人妻夜夜爽99麻豆av| 久久精品国产清高在天天线| 麻豆国产av国片精品| 人妻丰满熟妇av一区二区三区| 日韩欧美国产一区二区入口| 一个人免费在线观看的高清视频| 99热这里只有精品一区 | 人妻夜夜爽99麻豆av| 国产三级黄色录像| av女优亚洲男人天堂 | 亚洲精品美女久久av网站| 欧美3d第一页| 日韩精品中文字幕看吧| 长腿黑丝高跟| 国产精品久久电影中文字幕| 老熟妇仑乱视频hdxx| 午夜福利在线观看吧| 99久久精品国产亚洲精品| 日日摸夜夜添夜夜添小说| 91av网一区二区| a级毛片a级免费在线| 少妇的丰满在线观看| 无限看片的www在线观看| 长腿黑丝高跟| 两个人视频免费观看高清| 非洲黑人性xxxx精品又粗又长| 精品不卡国产一区二区三区| 91av网站免费观看| 中文在线观看免费www的网站| 亚洲av成人不卡在线观看播放网| 岛国在线免费视频观看| av黄色大香蕉| 啦啦啦免费观看视频1| 欧美日韩福利视频一区二区| 国产毛片a区久久久久| 偷拍熟女少妇极品色| 亚洲av成人精品一区久久| 日韩欧美精品v在线| 午夜免费激情av| 可以在线观看毛片的网站| 热99在线观看视频| 欧美国产日韩亚洲一区| 婷婷六月久久综合丁香| 亚洲第一电影网av| 日韩欧美一区二区三区在线观看| 亚洲第一电影网av| 91麻豆精品激情在线观看国产| 成人18禁在线播放| 国产精品久久久人人做人人爽| 人妻久久中文字幕网| 午夜影院日韩av| 国产97色在线日韩免费| 亚洲成人精品中文字幕电影| 亚洲成a人片在线一区二区| 精品人妻1区二区| 母亲3免费完整高清在线观看| 国产精品一区二区免费欧美| 脱女人内裤的视频| 亚洲 欧美一区二区三区| 亚洲国产精品999在线| 国产亚洲精品久久久com| 波多野结衣巨乳人妻| 午夜福利欧美成人| 香蕉国产在线看| 久久精品亚洲精品国产色婷小说| 国产蜜桃级精品一区二区三区| 女警被强在线播放| 免费无遮挡裸体视频| 国产一区二区激情短视频| 欧美日韩亚洲国产一区二区在线观看| 两个人的视频大全免费| 天天添夜夜摸| 欧美一级毛片孕妇| 午夜福利高清视频| 99久久无色码亚洲精品果冻| 小蜜桃在线观看免费完整版高清| 国产高潮美女av| 两个人视频免费观看高清| 怎么达到女性高潮| 无遮挡黄片免费观看| 欧美日韩福利视频一区二区| 国产乱人伦免费视频| 亚洲无线在线观看| 日本一本二区三区精品| 丁香欧美五月| 麻豆av在线久日| 国产熟女xx| 精品国产三级普通话版| 国产日本99.免费观看| 亚洲精品一卡2卡三卡4卡5卡| 成人特级av手机在线观看| 观看免费一级毛片| 欧美又色又爽又黄视频| 高清毛片免费观看视频网站| 99久久99久久久精品蜜桃| 国产av不卡久久| 校园春色视频在线观看| 国产黄片美女视频| 黄色视频,在线免费观看| 色综合欧美亚洲国产小说| 日韩人妻高清精品专区| 美女cb高潮喷水在线观看 | 丁香欧美五月| 毛片女人毛片| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人澡欧美一区二区| 成年女人永久免费观看视频| 日韩欧美在线二视频| 亚洲精品在线美女| 亚洲性夜色夜夜综合| 免费观看精品视频网站| 亚洲片人在线观看| 免费看日本二区| 亚洲中文字幕一区二区三区有码在线看 | 两个人视频免费观看高清| 91老司机精品| 亚洲av成人精品一区久久| 成人三级黄色视频| 精品电影一区二区在线| 国产免费男女视频| 婷婷亚洲欧美| 婷婷六月久久综合丁香| 99久国产av精品| 久久香蕉精品热| 亚洲性夜色夜夜综合| 国内精品久久久久精免费| 国产高潮美女av| 亚洲aⅴ乱码一区二区在线播放| 成人一区二区视频在线观看| 一区二区三区激情视频| 一二三四社区在线视频社区8| 久久久国产成人精品二区| 两个人的视频大全免费| 国产精华一区二区三区| 香蕉久久夜色| 淫秽高清视频在线观看| 国产av麻豆久久久久久久| 国产精品香港三级国产av潘金莲| 夜夜夜夜夜久久久久| 精品国产亚洲在线| 制服丝袜大香蕉在线| 国产熟女xx| 97超视频在线观看视频| 午夜福利成人在线免费观看| 成人三级黄色视频| 99久久国产精品久久久| 熟女电影av网| 在线观看午夜福利视频| 欧美乱码精品一区二区三区| or卡值多少钱| 一本精品99久久精品77| 国产蜜桃级精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 老汉色∧v一级毛片| 午夜福利视频1000在线观看| 欧美又色又爽又黄视频| 很黄的视频免费| 精华霜和精华液先用哪个| 制服人妻中文乱码| 中文字幕熟女人妻在线| 亚洲av成人一区二区三| 嫩草影院入口| 99视频精品全部免费 在线 | 日本黄大片高清| 禁无遮挡网站| 18禁观看日本| 精品99又大又爽又粗少妇毛片 | 亚洲无线观看免费| 淫秽高清视频在线观看| 在线观看舔阴道视频| 99国产精品一区二区蜜桃av| 女同久久另类99精品国产91| 热99re8久久精品国产| 亚洲国产看品久久| 性欧美人与动物交配| 香蕉av资源在线| 国产人伦9x9x在线观看| 精品久久久久久久毛片微露脸| 亚洲欧美日韩东京热| 青草久久国产| 欧美一区二区精品小视频在线| 亚洲精品456在线播放app | 人妻丰满熟妇av一区二区三区| 变态另类丝袜制服| 两性午夜刺激爽爽歪歪视频在线观看| 国产成年人精品一区二区| 欧美成人免费av一区二区三区| 黑人欧美特级aaaaaa片| 在线观看免费视频日本深夜| 国内毛片毛片毛片毛片毛片| 国产成+人综合+亚洲专区| 最近视频中文字幕2019在线8| 亚洲在线自拍视频| 美女大奶头视频| 国产日本99.免费观看| 欧美日韩黄片免| 一本一本综合久久| 亚洲国产高清在线一区二区三| 久久久久九九精品影院| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品国产精品久久久不卡| 淫妇啪啪啪对白视频| 一进一出抽搐gif免费好疼| 少妇人妻一区二区三区视频| 一本久久中文字幕| 国产乱人伦免费视频| 18禁美女被吸乳视频| 手机成人av网站| 麻豆久久精品国产亚洲av| 最近最新免费中文字幕在线| 欧美激情久久久久久爽电影| 欧美在线一区亚洲| 精品福利观看| 欧美xxxx黑人xx丫x性爽| 成人特级黄色片久久久久久久| 一本久久中文字幕| 国产午夜福利久久久久久| 久久性视频一级片| 亚洲一区二区三区色噜噜| 青草久久国产| 欧美日本亚洲视频在线播放| 两个人看的免费小视频| 久久中文字幕人妻熟女| 午夜福利在线观看免费完整高清在 | 男人舔女人的私密视频| 成人亚洲精品av一区二区| 国产午夜精品久久久久久| 欧美日韩精品网址| 日韩有码中文字幕| 黄色女人牲交| 色综合婷婷激情| 国产精品久久视频播放| 女人高潮潮喷娇喘18禁视频| 制服人妻中文乱码| 黄频高清免费视频| 三级国产精品欧美在线观看 | 国产高清视频在线观看网站| 国产精品电影一区二区三区| 夜夜爽天天搞| 久久国产精品影院| 午夜福利在线观看吧| 亚洲中文字幕一区二区三区有码在线看 | 久久久久性生活片| 91av网一区二区| 亚洲在线自拍视频| 99精品欧美一区二区三区四区| 中文字幕av在线有码专区| 免费大片18禁| 国产熟女xx| а√天堂www在线а√下载| 不卡av一区二区三区| 嫩草影院入口| 国产精品一区二区精品视频观看| 亚洲欧美日韩卡通动漫| 日本成人三级电影网站| 久久性视频一级片| 亚洲 国产 在线| 免费在线观看亚洲国产| 亚洲自拍偷在线| 久久久国产成人精品二区| 久久这里只有精品中国| 久久伊人香网站| 国产欧美日韩精品亚洲av| 久久久久久久久中文| 精品免费久久久久久久清纯| 免费观看精品视频网站| 男女床上黄色一级片免费看| 国产欧美日韩精品一区二区| 叶爱在线成人免费视频播放| 啦啦啦观看免费观看视频高清| 99在线人妻在线中文字幕| 99国产精品99久久久久| 国产精品女同一区二区软件 | 村上凉子中文字幕在线| 成人18禁在线播放| 国产精品一区二区三区四区久久| 国产人伦9x9x在线观看| 老汉色av国产亚洲站长工具| 久久久久免费精品人妻一区二区| 一夜夜www| 久久国产精品影院| 99在线视频只有这里精品首页| 男女午夜视频在线观看| 欧美丝袜亚洲另类 | 亚洲av熟女| 老熟妇乱子伦视频在线观看| 真人一进一出gif抽搐免费| 999精品在线视频| 欧美丝袜亚洲另类 | tocl精华| 女人高潮潮喷娇喘18禁视频| 国产精品亚洲av一区麻豆| 国产免费av片在线观看野外av| 精品国产乱子伦一区二区三区| 国产精品av久久久久免费|