近年来,高通量测序技术[1,2,3]已成为分析土壤根际微生物多样性研究的方法之一[4,5]。巨菌草（Pennisetum giganteum）和绿洲一号（Arundo donax cv.Lvzhou No.1）均为禾本科（Gramineae）年生草本植物,国内（河南、内蒙古、西藏、陕西、青海、甘肃等地区）、国外（卢旺达、斐济、莱索托和南非等国）均有分布[6]。这2种菌草不仅可作为食用菌、药用菌的培养基[7],还可作为动物饲料[8,9]、生物能源物质[10,11]及具有生态治理作用[12,13,14,15,16]。随着人们生活质量的提高、生态环境的改善及畜牧业的发展,迫切需要开展提高菌草产量和质量等方面的研究。
植物与根际土壤微生物的关系十分密切,根际土壤微生物种类和数量对植物营养养分循环[17]过程具有重要意义。根际土壤微生物影响植物根系的发育、促进植物的生长及抑制或减轻植物病害的发生等起一定的作用[18,19,20],其多样性及其生态功能成为目前国内外的研究热点[21,22]。真菌是土壤微生物的主要组成部分之一,对植物的生长发育等方面有着重要的生物学功能,如防治植物病害[23,24]、促进植株生长[25]和增强植物抗逆性[26]等。国内外学者对植物根际土真菌群落多样性等方面进行了研究,常见的根际土真菌属有镰刀菌属（Fusarium）、青霉属（Penicillium）、曲霉属（Aspergillus）等[27];从白菜软腐病根部土壤分离出草茎点霉（Phoma herbarum F217-1）能拮抗细菌性软腐病[28]。
菌草作为具有广泛用途植物之一,其根际土壤真菌多样性的相关研究尚未见报道。为探讨同一生长条件下不同菌草根际土壤微生物群落多样性,以福建农林大学菌草基地的巨菌草和绿洲一号根际土壤为研究对象,采用Illumina Miseq高通量测序技术对真菌ITS3-ITS4区片段进行测序并进行生物信息学分析,从分子水平上揭示巨菌草和绿洲一号根际土壤真菌群落的多样性,以期为深入研究巨菌草和绿洲一号根际土壤真菌潜在的生态功能研究提供一定的理论参考。
1 材料与方法
1.1 材料
样品采自福建农林大学菌草试验基地根际土壤。采集时随机选取巨菌草和绿洲一号植株各5株,抖落根系外围土壤,收取附着于根系表面的土壤做为根际土,并同时采取其同一地块的巨菌草和绿洲一号共同对照非根际土壤。带回实验室立即放入-80℃冰箱内备用。
1.2 方法
1.2.1 样品的采集 本研究将采自福建农林大学菌草基地菌草根系表面0~5 cm范围内10~20 cm土层的土壤视为根际土壤,而非根际土壤样品为距离菌草根系表面10~15 cm范围内10~20 cm土层的土壤。随机选取3个样点,多点采集的根际土和非根际土分别充分混合后,置于无菌自封袋中。将土壤样本带回实验室,置于-80 ℃冰箱中保存,备用。
1.2.2 土壤总DNA的提取 土壤DNA提取采用试剂盒（Fast DNA Spin Kit for Soil）（MP Biomedicals,USA）,具体步骤参照提取说明书[29]。用1%的琼脂糖凝胶电泳（0.5×TAE）检测DNA的质量,质量合格的DNA于-20 ℃冰箱保存备用。
1.2.3 土壤真菌ITS序列的PCR扩增及高通量测序 以稀释后的土壤DNA为模板,采用真菌通用引物ITS3-2024F（5°-GCATCGATGAAGAACGCAGC-3°）和ITS4-2409R（5°-TCCTCCGCTTATTGATATGC-3°）PCR扩增真菌ITS3-ITS4区域。PCR总反应体系为30 μL（Phusion Master Mix（2×）15 μL、Primer（2 μmol/L）3 μL（6 μmol/L）、DNA模板10 μL（5~10 ng）、H2O 2 μL）。反应程序：98 ℃预变性1 min;30 个循环包括（98 ℃,10 s;50 ℃,30 s;72 ℃,30 s）;72 ℃,5 min,3次重复。将同一样品PCR 产物混合后经电泳检测、割胶纯化、浓度检测等操作后,将样品送到诺禾致源生物有限公司,使用New England Biolabs 公司的NEB Next® Ultra™ DNA Library Prep Kit for Illumina建库试剂盒进行文库的构建,构建好的文库经过Qubit 定量和文库检测,合格后,使用HiSeq进行上机测序。
1.3 数据处理
利用Uparse软件[35]对所有样品的全部 Effective Tags进行聚类,默认以97%的一致性（identity）将序列聚类成为OTUs（operational taxonomic units）,同时会选取OTUs的代表性序列进行物种注释,用QIIME软件中的blast方法[36] 与Unit数据库[37]进行物种注释分析,并分别在门（phylum）、纲（class）、属（genus）分类水平统计各样本的群落组成。使用MUSCLE[38]软件进行快速多序列比对,得到所有OTUs代表序列的系统发生关系。最后以样品中数据量最少的为标准对各样品的数据进行均一化处理,后续的Alpha多样性分析和Beta多样性分析都是基于均一化处理后的数据。
使用QIIME软件计算Observed-species、Chao1、Shannon、Simpson、ACE、Goods-coverage,使用R软件绘制稀释曲线、Rank abundance曲线;使用R软件绘制PCoA图。
2 结果与分析
2.1 测序数据处理与分析
质量通过Illumina MiSeq高通量测序,9份土壤样品质控后得到了 525 292条有效序列,经过过滤和去除嵌合体之后得到了511 816条优质序列,优质序列的长度集中在200~450 bp之间,以97%相似度划分,共得到7267个OTU。
2.2 Chao1指数稀释曲线和Observed species指数稀释曲线
图1 多样性指数稀释曲线图A：Observed species 指数稀释曲线图;B：Chao1 指数稀释曲线图。 Fig. 1 Rarefaction curves of Alpha diversity index A: Dilution curve of observed species index; B: Chao1 index dilution curve. |
表1 菌草根际土壤真菌群落的多样性指数Tab. 1 Diversity index of samples |
| 样品Samples | 文库覆盖率 Coverage/% | 香浓-威纳指数Shannon-Wiener | Simpson 指数 Simpson index | ACE指数 ACE index | Chao1指数 Chao1 index |
|---|---|---|---|---|---|
| FSLV | 99.73±0.0003a | 6.6507±0.061a | 0.0323±0.0029c | 934.41±45.161a | 927.12±34.567a |
| FSJU | 99.73±0.0003a | 5.1013±0.084b | 0.1163±0.0047b | 754.76±49.140ab | 779.78±54.024ab |
| FSCK | 99.63±0.0003a | 4.3937±0.094c | 0.1763±0.0253a | 739.43±64.780b | 735.00±68.438b |
注：同列数据后不同小写字母表示差异显著（P<0.05）。 | |
Note: Different lowercase letters in the same column indicate significant difference at 0.05 level. |
2.3 2种菌草土壤真菌群落结构变化
2.3.1 门水平群落结构分析 由图2可知,巨菌草根际土壤真菌群落在门类水平上主要类群为子囊菌门（Ascomycota）、担子菌亚门（Basidiomycota）、接合菌门（Zygomycota）、壶菌门（Chytridiomycota）,相对丰度分别为59.51%、13.39%、1.26%、0.04%;绿洲一号根际土壤真菌群落在门类水平上主要类群为子囊菌门（Ascomycota）、担子菌亚门（Basidiomycota）、接合菌门（Zygomycota）、壶菌门（Chytridiomycota）,相对丰度分别为64.37%、20.91%、0.83%、0.29%;对照非根际土壤真菌群落在门类水平上主要类群为子囊菌门（Ascomycota）、担子菌亚门（Basidiomycota）、接合菌门（Zygomycota）、壶菌门（Chytridiomycota）,相对丰度分别为77.24%、7.02%、6.18%、0.05%。综上所述,菌草根际土和非根际土在门水平上丰富度最高的是子囊菌门（Ascomycota）。
2.3.2 纲水平群落结构分析 由图3可知,巨菌草根际土壤真菌群落在纲类水平上主要类群散囊菌纲（Eurotiomycetes）、座囊菌纲（Dothideomycetes）、粪壳菌纲（Sordariomycetes）、伞菌纲（Agaricomycetes）、锤舌菌纲（Leotiomycetes）、银耳纲（Tremellomycetes）,其相对丰度分别为2.72 %、2.75%、19.80%、12.03%、0.76%、0.67%;绿洲一号根际土壤真菌群落在纲类水平上主要类群散囊菌纲（Eurotiomycetes）、座囊菌纲（Dothideomycetes）、粪壳菌纲（Sordariomycetes）、伞菌纲（Agaricomycetes）、锤舌菌纲（Leotiomycetes）、银耳纲（Tremellomycetes）,其相对丰度分别为9.51 %、15.51%、28.17%、16.97%、1.82%、2.61%;对照非根际土壤真菌群落在纲类水平上主要类群散囊菌纲（Eurotiomycetes）、座囊菌纲（Dothideomycetes）、粪壳菌纲（Sordariomycetes）、伞菌纲（Agaricomycetes）、锤舌菌纲（Leotiomycetes）、银耳纲（Tremellomycetes）,其相对丰度分别为55.44 %、4.67%、6.29%、2.10%、8.85%、3.95%。综上所述,菌草根际土和非根际土在纲水平上丰富度最高的是粪壳菌纲（Sordariomycetes）。
2.3.3 属水平群落结构分析 由图4可知,巨菌草根际土壤真菌群落在属水平上主要类群为篮状菌属（Talaromyces）、深黄伞形霉属（Umbelopsis）、枝孢菌属（Cladosporium）、暗双孢菌属（Cordana）、镰刀菌属（Fusarium）、隐球菌属（Cryptococcus）、淡紫紫霉属（Purpureocillium）、 赤霉菌属（Gibberella）、灵芝属（Ganoderma）,相对丰度分别为0.21%、0.62%、1.45%、1.21%、1.72%、0.38%、0.50%、0.10%、0.36%;绿洲一号根际土壤真菌群落在属水平上主要类群为篮状菌属（Talaromyces）、深黄伞形霉属（Umbelopsis）、枝孢菌属（Cladosporium）、暗双孢菌属（Cordana）、镰刀菌属（Fusarium）、隐球菌属（Cryptococcus）、淡紫紫霉属（Purpureocillium）、赤霉菌属（Gibberella）、灵芝属（Ganoderma）,相对丰度分别为0.77%、1.10%、5.58%、2.68%、3.74%、2.50%、3.05%、0.03%、0.96%;非根际土壤真菌群落在属水平上主要类群为篮状菌属（Talaromyces）、深黄伞形霉属（Umbelopsis）、枝孢菌属（Cladosporium）、暗双孢菌属（Cordana）、镰刀菌属（Fusarium）、隐球菌属（Cryptococcus）、淡紫紫霉属（Purpureocillium）、赤霉菌属（Gibberella）、灵芝属（Ganoderma）,相对丰度分别为27.86%、5.20%、1.92%、0.16%、0.83%、3.73%、0.29%、0.23%、0.17%。综上所述,巨菌草根际土壤真菌群落在属水平上丰富度最高的类群为镰刀菌属（Fusarium）,绿洲一号根际土壤真菌群落在属水平上丰富度最高的类群为枝孢菌属（Cladosporium）,而非根际土壤真菌群落在属水平上丰富度最高的类群为篮状菌属（Talaromyces）。
2.4 2种菌草真菌群落丰度和多样性分析
Shannon指数和Simpson指数用于反映土壤样品中微生物的多样性。由表1可知,与对照组相比,巨菌草和绿洲一号土壤真菌Shannon指数显著高于非根际土壤,Simpson指数显著低于非根际土壤;绿洲一号土壤样本的Shannon-Wiener指数高于巨菌草根际土壤,为6.6507,而Simpson指数低于巨菌草根际土壤,为0.0323。可见,巨菌草和绿洲一号根际土壤样本中土壤的真菌种群多样性高于非根际土壤。绿洲一号根际土壤样本的真菌群落多样性略高于巨菌草根际土壤样本。
ACE指数和Chao1指数反映群落物种丰富度。由表1可知,与对照组相比,巨菌草和绿洲一号根际土壤真菌ACE指数、Chao1指数均比非根际土壤高,差异显著;绿洲一号根际土壤真菌ACE指数、Chao1指数均比比巨菌草根际土高,差异显著。结果表明,巨菌草和绿洲一号根际土壤样本中土壤的真菌丰度高于非根际土,绿洲一号根际土壤真菌丰富度高于巨菌草根际土。
2.5 2种菌草真菌群落的相关性分析
2.6 基于OTUs的主坐标分析及样品间的NMDS分析
为了解不同菌草根际土壤微生物群落组成,对样本的土壤真菌群落进行主坐标分析（PCoA）。不同样本表现出分散或聚集的分布情况,结果如图6所示。图6A中绿洲一号与PC2轴相距均较近而与PC1相对较远,说明绿洲一号中真菌菌群结构受主成分PC2的影响大,PC2一种主成分对样品中菌群结构的影响达到总影响因子的24.8%;巨菌草与PC1和PC2轴相距均较近,说明巨菌草中真菌菌群结构受主成分PC1和PC2的影响均大,PC1和PC2两种主成分对样品中菌群结构的影响分别达到总影响因子的58.91%和24.8%;对照与PC2轴相距较远,说明FSCK样品中真菌菌群结构受主成分PC2的影响不大。通过主坐标分析可见,真菌群落PC1和PC2可分别解释59.81%和24.8%差异性。绿洲一号、巨菌草及对照没有聚在一起,说明微生物群落差异大。
无度量多维标定法（Non-Metric Multi-Dimensional Scaling, NMDS）统计是一种适用于生态学研究的排序方法。NMDS是非线性模型,其设计目的是为克服线性模型（包括PCA、 PCoA）的缺点,更好地反映生态学数据的非线性结构。应用NMDS分析,根据样本中包含的物种信息,以点的形式反映在多维空间上,而对不同样本间的差异程度,则是通过点与点间的距离体现,能够反映样本的组间和组内差异等。由（图6B）可知,3份样品中,代表绿洲一号（FSLV）的3个样品距离其他样品较远,表明该样品的真菌群落结构存在彼此之间差异较大。
3 讨论
巨菌草和绿洲一号是具有广泛作用的代表性菌草,在水土保持、荒漠化治理及提高落后国家人民的生活水平等方面具有重要的作用。植物根际土壤微生物参与土壤养分的转化和循环,可提高土壤肥力[39]。土壤微生物的种类和数量影响土壤养分的效能,进一步影响植物的生长。土壤健康的指标之一就是土壤微生物群落的组成及其多样性。
目前菌草根际土微生物多样性的研究还未见报导。菌草试验点选择在同一地点、管理一致的不同菌草,因此,巨菌草和绿洲一号根际真菌群落多样性的差异主要反映了不同菌草所产生的根际效应。高通量测序技术作为一种简单、快速、准确的免培养分子生物学技术,为认识土壤微生物多样性、群落结构组成及其生态功能提供有利手段[40]。本研究利用高通量测序技术分析了福州菌草基地的巨菌草和绿洲一号的根际土壤与非根际土壤真菌种群结构的多样性,为菌草根际微生物的研究奠定基础。菌草根际土壤样本的ACE和Shannon-iener指数值均大于非根际土壤样本,表明根际土壤真菌群落在丰富度和多样性上均高于其非根际土壤。绿洲一号根际土壤样本的ACE和Shannon- iener指数值大于巨菌草根际土壤样本,表明绿洲一号根际土壤真菌群落在丰富度和多样性上均高于巨菌草根际土壤。从门水平来看,3个样品的优势门为子囊菌门（Ascomycota）,与一些学者的研究结果一致[41,42,43,44]。其说明子囊菌门（Ascomyota）真菌适应性广泛,在不同的土壤条件下都有分布。在属水平上,篮状菌属（Talaromyces）、深黄伞形霉属（Umbelopsis）、枝孢菌属（Cladosporium）、暗双孢菌属（Cordana）、镰刀菌属（Fusarium）、隐球菌属（Cryptococcus）、淡紫紫霉属（Purpureoillium）、赤霉菌属（Gibberella）、灵芝属（Ganoderma）是所有样品土壤真菌Top10 属中的共有属。姜华等[45]研究表明,芦苇根际土壤真菌的优势属有链格孢属(Alternaria)、青霉属(Penicillium)、拟青霉属(Paecilomyces)、木霉属(Trichoderma)和链格孢属(Alternaria)。这说明不同植物根际土壤优势真菌在属水平上是不同的[45]。
本研究结果表明,巨菌草和绿洲一号根际土真菌群落多样性与丰富度均高于非根际土。绿洲一号根际土真菌群落多样性与丰富度略高于巨菌草根际土真菌群落多样性与丰富度。