Metabarcoding of fungal communities associated with bark beetles

Many species of fungi are closely allied with bark beetles, including many tree pathogens, but their species richness and patterns of distribution remain largely unknown. We established a protocol for metabarcoding of fungal communities directly from total genomic DNA extracted from individual beetles, showing that the ITS3/4 primer pair selectively amplifies the fungal ITS. Using three specimens of bark beetle from different species, we assess the fungal diversity associated with these specimens and the repeatability of these estimates in PCRs conducted with different primer tags. The combined replicates produced 727 fungal Operational Taxonomic Units (OTUs) for the specimen of Hylastes ater, 435 OTUs for Tomicus piniperda, and 294 OTUs for Trypodendron lineatum, while individual PCR reactions produced on average only 229, 54, and 31 OTUs for the three specimens, respectively. Yet, communities from PCR replicates were very similar in pairwise comparisons, in particular when considering species abundance, but differed greatly among the three beetle specimens. Different primer tags or the inclusion of amplicons in separate libraries did not impact the species composition. The ITS2 sequences were identified with the Lowest Common Ancestor approach and correspond to diverse lineages of fungi, including Ophiostomaceae and Leotiomycetes widely found to be tree pathogens. We conclude that Illumina MiSeq metabarcoding reliably captures fungal diversity associated with bark beetles, although numerous PCR replicates are recommended for an exhaustive sample. Direct PCR from beetle DNA extractions provides a rapid method for future surveys of fungal species diversity and their associations with bark beetles and environmental variables.     

High diversity and widespread occurrence of mitotic spore mats in ectomycorrhizal Pezizales

Fungal mitospores may function as dispersal units and/ or spermatia and thus play a role in distribution and/or mating of species that produce them. Mitospore production in ectomycorrhizal (EcM) Pezizales is rarely reported, but here we document mitospore production by a high diversity of EcM Pezizales on three continents, in both hemispheres. We sequenced the internal transcribed spacer (ITS) and partial large subunit (LSU) nuclear rDNA from 292 spore mats (visible mitospore clumps) collected in Argentina, Chile, China, Mexico and the USA between 2009 and 2012. We collated spore mat ITS sequences with 105 fruit body and 47 EcM root sequences to generate operational taxonomic units (OTUs). Phylogenetic inferences were made through analyses of both molecular data sets. A total of 48 OTUs from spore mats represented six independent EcM Pezizales lineages and included truffles and cup fungi. Three clades of seven OTUs have no known meiospore stage. Mitospores failed to germinate on sterile media, or form ectomycorrhizas on Quercus, Pinus and Populus seedlings, consistent with a hypothesized role of spermatia. The broad geographic range, high frequency and phylogenetic diversity of spore mats produced by EcM Pezizales suggests that a mitospore stage is important for many species in this group in terms of mating, reproduction and/or dispersal.     

Assessment of fungal diversity in deep-sea sediments by multiple primer approach

Increasing evidence of the fungal diversity in deep-sea sediments has come from amplification of environmental DNA with fungal specific or eukaryote primer sets. In order to assess the fungal diversity in deep-sea sediments of the Central Indian Basin (CIB) at ~5,000 m depth, we amplified sediment DNA with four different primer sets. These were fungal-specific primer pair ITS1F/ITS4 (internal transcribed spacers), universal 18S rDNA primers NS1/NS2, Euk18S-42F/Euk18S-1492R and Euk18S-555F/Euk18S-1269R. One environmental library was constructed with each of the primer pairs, and 48 clones were sequenced per library. These sequences resulted in 8 fungal Operational Taxonomic Units (OTUs) with ITS and 19 OTUs with 18S rDNA primer sets respectively by taking into account the 2% sequence divergence cut-off for species delineation. These OTUs belonged to 20 distinct fungal genera of the phyla Ascomycota and Basidiomycota. Seven sequences were found to be divergent by 79–97% from the known sequences of the existing database and may be novel. A majority of the sequences clustered with known sequences of the existing taxa. The phylogenetic affiliation of a few fungal sequences with known environmental sequences from marine and hypersaline habitat suggests their autochthonous nature or adaptation to marine habitat. The amplification of sequences belonging to Exobasidiomycetes and Cystobasidiomycetes from deep-sea is being reported for the first time in this study. Amplification of fungal sequences with eukaryotic as well as fungal specific primers indicates that among eukaryotes, fungi appear to be a dominant group in the sampling site of the CIB.     

青藏高原食草动物粪栖真菌的多样性

分别利用真菌通用引物和厌氧真菌特异引物构建了西藏地区3种反刍动物和1种单胃动物共8份新鲜粪样的ITS克隆文库,以通过系统发育分析解析其中好氧真菌与厌氧真菌的多样性。通用引物ITS文库测得324条真菌序列,分别属于子囊菌门Ascomycota 3个目、担子菌门Basidiomycota 2个目、接合菌门Zygomycota的1个目和严格厌氧的新丽鞭毛菌门Neocallimastigomycota,共24个OTUs。其中,子囊菌门相对丰度最高,占80.6%;新丽鞭毛菌门相对丰度最低,仅占0.6%。大部分OTUs与已知真菌属、种关系较远。厌氧真菌特异引物文库测得661条序列,全部属于新丽鞭毛菌门,包括所有已知的厌氧霉属Anaeromyces、盲肠菌属Caecomyces、肠霉属Cyllamyces、新丽鞭毛属Neocallimastix、奥式霉属Orpinomyces、胃梨囊霉属Piromyces 6个属和3个未培养的属级类群(NG9、NG10、NG11),共29个OTUs。其中3个已知的单中心属存在于所有反刍动物样品中,并以Piromyces相对丰度最高(37.4%)。单胃动物马粪样中全部为NG9类群。NG9是本研究新发现的属级类群,研究中同时揭示有多个未培养种和潜在的新种。研究结果证明青藏高原反刍动物粪栖真菌多样性较高,并存在丰富的未培养种和潜在的新属及新种。     

Fungal Community Analysis by Large-Scale Sequencing of Environmental Samples

Fungi are an important and diverse component of soil communities, but these communities have proven difficult to study in conventional biotic surveys. We evaluated soil fungal diversity at two sites in a temperate forest using direct isolation of small-subunit and internal transcribed spacer (ITS) rRNA genes by PCR and high-throughput sequencing of cloned fragments. We identified 412 sequence types from 863 fungal ITS sequences, as well as 112 ITS sequences from other eukaryotic microorganisms. Equal proportions of Basidiomycota and Ascomycota sequences were present in both the ITS and small-subunit libraries, while members of other fungal phyla were recovered at much lower frequencies. Many sequences closely matched sequences from mycorrhizal, plant-pathogenic, and saprophytic fungi. Compositional differences were observed among samples from different soil depths, with mycorrhizal species predominating deeper in the soil profile and saprophytic species predominating in the litter layer. Richness was consistently lowest in the deepest soil horizon samples. Comparable levels of fungal richness have been observed following traditional specimen-based collecting and culturing surveys, but only after much more extensive sampling. The high rate at which new sequence types were recovered even after sampling 863 fungal ITS sequences and the dominance of fungi in our libraries relative to other eukaryotes suggest that the abundance and diversity of fungi in forest soils may be much higher than previously hypothesized.     

Fungal palaeodiversity revealed using high‐throughput metabarcoding of ancient DNA from arctic permafrost

The taxonomic and ecological diversity of ancient fungal communities was assessed by combining next generation sequencing and metabarcoding of DNA preserved in permafrost. Twenty‐six sediment samples dated 16 000–32 000 radiocarbon years old from two localities in Siberia were analysed for fungal ITS. We detected 75 fungal OTUs from 21 orders representing three phyla, although rarefaction analyses suggested that the full diversity was not recovered despite generating an average of 6677 ± 3811 (mean ± SD) sequences per sample and that preservation bias likely has considerable effect on the recovered DNA. Most OTUs (75.4%) represented ascomycetes. Due to insufficient sequencing depth, DNA degradation and putative preservation biases in our samples, the recovered taxa probably do not represent the complete historic fungal community, and it is difficult to determine whether the fungal communities varied geographically or experienced a composition shift within the period of 16 000–32 000 bp . However, annotation of OTUs to functional ecological groups provided a wealth of information on the historic communities. About one‐third of the OTUs are presumed plant‐associates (pathogens, saprotrophs and endophytes) typical of graminoid‐ and forb‐rich habitats. We also detected putative insect pathogens, coprophiles and keratinophiles likely associated with ancient insect and herbivore faunas. The detection of putative insect pathogens, mycoparasites, aquatic fungi and endophytes broadens our previous knowledge of the diversity of fungi present in Beringian palaeoecosystems. A large group of putatively psychrophilic/psychrotolerant fungi was also detected, most likely representing a modern, metabolically active fungal community.     

不同扩增引物对高通量测序分析盐角草内生真菌多样性的影响

【背景】高通量测序分析作为深入了解环境微生物群落组成的重要方法,已成为植物内生真菌多样性研究的有效手段,然而由于引物的扩增差异,采用不同引物可对实验结果分析造成影响。同时,盐角草作为世界上最耐盐的植物之一,存在着多种功能性的内生真菌,而较为全面介绍其内生真菌组成和多样性的报道鲜见。【目的】为了揭示盐角草内生真菌的多样性,解析不同扩增引物对内生菌多样性分析的影响。【方法】分别采用真菌高通量测序常用引物对ITS1-5F、ITS1-1F、ITS2对采自乌鲁木齐达坂城盐湖的盐角草内生真菌进行扩增,开展其内生真菌OTU的分析。【结果】通过不同引物对扩增并测序共获得102个盐角草内生真菌OTU,涉及真菌界8个门和未分类菌群,其中子囊菌门(Ascomycota)占绝对优势,其次为担子菌门(Basidiomycota);在属层次上,盐角草内生真菌共涉及64个属及20个未分类属,其中Alternaria、Cladosporium、Podospora等3个属为盐角草内生真菌优势菌群。对不同引物对扩增测序结果分析表明,不同引物对扩增对分析内生真菌OTU数量和种类具有明显的影响,在全部所得的102个OTU中,ITS1-5F引物对获得44个OTU、ITS1-1F引物对获得55个OTU、ITS2引物对获得25个OTU,但以上3对引物扩增均检测到的OTU数仅为5个。物种组成和多样性分析表明,内生真菌多样性分析中采用以ITS1-1F为主,ITS1-5F为辅的分析策略,可较为全面地展现内生真菌的多样性。【结论】盐角草存在较为丰富的内生真菌资源,不同扩增引物对高通量分析盐角草内生真菌组成和分布具有明显的影响。     

Fungal community analysis by large-scale sequencing of environmental samples

Fungi are an important and diverse component of soil communities, but these communities have proven difficult to study in conventional biotic surveys. We evaluated soil fungal diversity at two sites in a temperate forest using direct isolation of small-subunit and internal transcribed spacer (ITS) rRNA genes by PCR and high-throughput sequencing of cloned fragments. We identified 412 sequence types from 863 fungal ITS sequences, as well as 112 ITS sequences from other eukaryotic microorganisms. Equal proportions of Basidiomycota and Ascomycota sequences were present in both the ITS and small-subunit libraries, while members of other fungal phyla were recovered at much lower frequencies. Many sequences closely matched sequences from mycorrhizal, plant-pathogenic, and saprophytic fungi. Compositional differences were observed among samples from different soil depths, with mycorrhizal species predominating deeper in the soil profile and saprophytic species predominating in the litter layer. Richness was consistently lowest in the deepest soil horizon samples. Comparable levels of fungal richness have been observed following traditional specimen-based collecting and culturing surveys, but only after much more extensive sampling. The high rate at which new sequence types were recovered even after sampling 863 fungal ITS sequences and the dominance of fungi in our libraries relative to other eukaryotes suggest that the abundance and diversity of fungi in forest soils may be much higher than previously hypothesized. All sequences were deposited in GenBank, with accession numbers AY 969316 to AY 970290 for the ITS sequences and AY 969135 to AY 969315 for the SSU sequences.     

Diversity and host recurrence of fungi associated with the bleached leaf litter in a subtropical forest

Fungi that selectively remove lignin cause extensive bleaching of leaf litter, which results in the acceleration of litter decomposition. The taxonomic diversity and host recurrence of saprotrophic fungi associated with bleached leaf litter were investigated in a subtropical forest in southern Japan. A total of 211 fungal isolates and sporocarps were obtained from bleached portions of leaf litter of 33 plant species in 18 families and were classified into 83 operational taxonomic units (OTUs) at the 97% similarity level of the ITS rDNA sequence. These fungal OTUs belonged to Rhytismataceae and Xylariaceae in the Ascomycota, and Marasmiaceae, Mycenaceae, Lachnocladiaceae, and Physalacriaceae in the Basidiomycota. OTUs in Rhytismataceae and Marasmiaceae showed a significantly higher degree of recurrence for plant species than simulated networks where partners were associated randomly. In contrast, OTUs in Xylariaceae and Mycenaceae showed no significant recurrence for plant species. Such differing degrees of recurrence for plant species implied different roles of fungal families in leaf litter decomposition.     

南海长茎葡萄蕨藻共附生真菌的多样性

为了揭示南海长茎葡萄蕨藻Caulerpa lentillifera共附生真菌的群落结构,并分析自来水清洗对真菌群落多样性的影响,本研究以收集于深圳南海各月份长茎葡萄蕨藻为材料,利用Illumina MiSeq 2500测序平台进行真菌ITS1扩增子测序,分析自来水清洗前后真菌类群组成和多样性的差异。研究发现,从所有18个样品中总共获得914个OTUs,检测出真菌3门12纲32目50科74属106种,其中优势属为曲霉属Aspergillus、单胞瓶霉属Phialemonium、枝孢属Cladosporium、链座菌属Catenulostroma和茎点霉属Phoma;不同月份真菌Alpha多样性指数显示南海长茎葡萄蕨藻共附生真菌丰度和多样性均随着温度降低而降低;自来水清洗前后南海长茎葡萄蕨藻共附生真菌的多样性及群落结构存在一定差异,清洗后真菌OTUs和Alpha多样性指数均有所降低,其中Chao1和Ace指数为显著性下降,表明自来水清洗对真菌群落组成丰度及多样性影响显著,这对长茎葡萄蕨藻的食品安全具有积极意义。     

Molecular characterization of airborne fungal spores in boreal forests of contrasting human disturbance

In this study we present a new approach to characterize fungal diversity with DNA sequencing of mycelium grown from trapped airborne spores. Fungal spores were extracted systematically from air in three boreal forest sites (clear-cut, young and old-growth forests) using an air sampling device. Internal transcribed spacer (ITS) sequences from the nuclear ribosomal DNA (nrDNA) were generated, and the sequences most likely taxon affinities were established through DNA homology searches. Phylogenetic analyses were used to classify similar sequences into operational taxonomic units (OTUs). The analyses indicated that a total of 84 different OTUs had been sampled, 24 basidiomycetes and 60 ascomycetes. OTUs belonging to the ascomycete orders Helotiales and Pleosporales were most frequent (31 and 18 respectively). A total of 54, 29 and 33 OTUs were sampled, respectively, in the old-growth, young and clear-cut forest sites. Although heavy generalization should be avoided due to few replicates, the results could indicate that old-growth boreal forests have significantly higher airborne fungal species richness than recently managed forests. The study shows that the spore-trapping approach has a great potential for targeting and studying anonymous fungi.     

Ectomycorrhizal fungi diversity in a white sand forest in western Amazonia

The genera Dicymbe and Aldina (Fabaceae) host ectomycorrhizal fungi (EcM) and are common in white sand forests (WSFs), a highly specialized habitat with a high level of plant endemism compared with terra-firme forests. In this study, we visited four times a 1-ha permanent plot established in a small patch of a WSF in the south of Colombia Amazonia. Forty-eight species of EcM fungi were recovered from sporocarps and 15 ITS species-level were detected from root tips. Seventeen species were new reports to Colombia and seven corresponded to undescribed species. These results confirm that this WSF supports a significant EcM fungal diversity. Most of the species found in this study have been previously reported to be associated with other legume and/or dipterocarp species from geographically distant forests. The long-distance occurrence combined with low host specificity, suggest the possibility of gene flow between geographically distant populations of EcM fungi in neotropical lowland rainforests.     

三种杓兰根相关真菌多样性和生态功能

【背景】除了菌根真菌(Orchid mycorrhizal fungi,OrMF)外,兰科植物根中还有其它内生真菌,称为根相关真菌(Root-associated fungi,RAF)。【目的】采用分离培养的方法获得同一栖息地针叶林和灌木林两种不同生境西藏杓兰、黄花杓兰和无苞杓兰的RAF菌株,研究其真菌谱系、多样性和生态功能结构。【方法】从杓兰根碎屑中分离RAF,通过总DNA提取、PCR扩增及测序得到ITS(Internaltranscribedspacer)序列;进行系统发育和多样性分析,并通过NCBI数据库比对得到相似性最高序列的注释信息来分析RAF生态学特性。【结果】共分离得到278株RAF,25种OTU类型,包括23个子囊菌门OTU,2个毛霉菌门OTU。RAF物种丰富度分析发现西藏杓兰的较黄花杓兰高,不同生境没有显著差异;不同杓兰物种较不同生境的RAF群落分化程度高。生态功能分析显示25个OTU包括共生型、腐生型和致病型3种营养型,以及外生菌根菌群、植物病原菌群、内生真菌群、动物病原菌群、真菌寄生菌群、杜鹃花类菌根群、未定义的腐生菌群和不确定型8种共位群。【结论】阐明不同生境采集的不同杓兰中RAF的分布特点和生态功能,为未来研究RAF与杓兰属植物的共生关系奠定基础。     

基于rDNA ITS序列的油茶叶片内生真菌多样性研究(英文)

油茶是我国热带和亚热带地区广泛栽培的一类重要的经济树种,其种子可生产优质食用油(茶油)。为了比较全面地了解油茶的内生真菌多样性,采用基于rDNA ITS的免培养法从油茶叶片中提取了总DNA,再从总DNA样本中直接扩增了内生真菌核糖体RNA基因内转录间隔区序列(ITS),进而构建了ITS克隆文库。共计从50个克隆子中获得了22种不同的克隆序列,其中3个为嵌合体,7个归属为植物,其余克隆序列根据序列相似性和系统发育分析归为12个不同的分类操作单元(OTUs),全部为子囊菌,分属于4纲4目,其中Aspergillus属2种,Mycosphaerellaceae科1种,Sordariales目4种,Helotiales目5种,Aspergillus属真菌是优势菌。结果表明油茶叶片内生真菌的种类分布较广。     

基于DNA复合条形码技术的蝗虫肠道共生真菌多样性研究

利用DNA复合条形码技术,研究了11个样本的蝗虫肠道共生真菌的多样性。结果显示:ITS在所研究的物种中鉴定了5门16纲29目40属2786 OTU真菌。肠道真菌群落组成分析结果显示:所有物种肠道真菌类群中含量最高的是木耳菌目和银耳目,其中斑翅蝗科的真菌类群多样性相对最高,斑腿蝗科的真菌类群多样性相对最低,表明各蝗虫肠道之间存在着明显的菌群多样性变化。α多样性分析结果显示:斑翅蝗科的共生真菌群落丰富度和多样性最高,斑腿蝗科的则最低。β多样性分析结果显示:(1)同科的各个种的肠道真菌群落结构差异性较小,不同科的种的肠道真菌群落结构差异性较大;(2)剑角蝗科的肠道真菌群落结构与其他物种的相似性均相对较低,而且在两个不同取样地得到的中华剑角蝗的真菌群落结构相似性也相对较低。聚类分析结果显示:(1)同科的蝗虫肠道真菌首先聚到一起,且群落相似性也相对较高;(2)布勒掷孢酵母属、内疣衣属和外瓶霉属3个属在蝗虫肠道真菌中是优势菌属。     

北京地区手参内生真菌的区系组成分析

手参Gymnadenia conopsea是一种地生兰科植物,也是中国传统中药材之一。研究以ITS1为目标序列,应用 Illumina MiSeq高通量测序技术对北京地区松山和百花山的手参内生真菌区系组成进行测定分析。结果表明,从手参的根中获得有效序列50 420条,209个可操作分类单元（OTUs）,分属于3门9纲28目54科,包括盘菌纲Pezizomycetes（3个OTU）、伞菌纲Agaricomycetes（21个OTU）、锤舌菌纲Leotiomycetes（18个OTU）、粪壳菌纲Sordariomycetes（51个OTU）、接合菌纲Zygomycetes（18个OTU）、银耳纲Tremellomycetes（19个OTU）、散囊菌纲Eurotiomycetes（22个OTU）、酵母纲Saccharomycetes（8个OTU）、座囊菌纲Dothideomycetes（14个OTU）。其中,百花山手参根中内生真菌优势类群是鸡油菌目Cantharellales（38.62%）和柔膜菌目Helotiales（24.64%）真菌;松山手参的优势内生菌类群是盘菌目Pezizales（66.19%）和鸡油菌目Cantharellales（29.75%）真菌。这些结果利于系统了解手参根部内生真菌区系组成以及优势内生真菌类群,为今后开展相关内生真菌在手参人工繁育中的应用奠定了基础。     

长茎羊耳蒜菌根真菌类群的研究

研究了广西雅长自然保护区和云南西双版纳自然保护区共3个产地的兰科植物羊耳蒜属长茎羊耳蒜Liparis viridiflora的菌根真菌类群区系组成.根内菌根真菌的核糖体基因内转录间隔区序列(rDNA-ITS)采用PCR技术扩增,克隆,测序并构建系统发育树.结果表明,长茎羊耳蒜根内所检测到的真菌大部分为胶膜菌科Tulasnellaceae真菌;根据序列相似性和系统发育分析,所有真菌可归为12个可操作分类单元(OTU),其中胶膜菌科有7个OTUs,达到总数的90.6%,为优势类群.菌根真菌多样性及区系组成在3个不同产地样本之间存在一定的差异;菌根真菌可能和兰科植物的生境适应性存在一定的相关性.     

Assessment of root-associated fungal communities colonizing two species of tropical grasses reveals incongruence to fungal communities of North American native grasses

This study characterized the root-associated fungal (RAF) communities inhabiting Eustachys petraea and Panicum maximum, two tropical and weedy grass species on an island off the coast of Honduras, Central America. Media-based analyses revealed that fungi colonizing these grasses exhibited similar morphotypes (albeit in different proportions). Conversely, molecular-based analyses suggested that the community of fungal OTUs were dissimilar, with several OTUs commonly present on only one plant and no AMF sequences present in either plant species. When compared to various datasets of RAF obtained from the North American mainland (with similar methods), the Honduran grasses harbored a different set of fungal OTUs. Interestingly, some cosmopolitan taxonomic clades commonly encountered on the mainland were also isolated in Honduras, though based on ITS rDNA the fungi were only related distantly to those on the mainland. Additional work is necessary to determine how the RAF communities change over geographic distance and whether invasive grasses are constitutionally set up to accept relationships with a different suite of RAF species.     

Changes in the root-associated fungal communities along a primary succession gradient analysed by 454 pyrosequencing

We investigated changes in the root-associated fungal communities associated with the ectomycorrhizal herb Bistorta vivipara along a primary succession gradient using 454 amplicon sequencing. Our main objective was to assess the degree of variation in fungal richness and community composition as vegetation cover increases along the chronosequence. Sixty root systems of B. vivipara were sampled in vegetation zones delimited by dated moraines in front of a retreating glacier in Norway. We extracted DNA from rinsed root systems, amplified the ITS1 region using fungal-specific primers and analysed the amplicons using 454 sequencing. Between 437 and 5063 sequences were obtained from each root system. Clustering analyses using a 98.5% sequence similarity cut-off yielded a total of 470 operational taxonomic units (OTUs), excluding singletons. Between eight and 41 fungal OTUs were detected within each root system. Already in the first stage of succession, a high fungal diversity was present in the B. vivipara root systems. Total number of OTUs increased significantly along the gradient towards climax vegetation, but the average number of OTUs per root system stayed unchanged. There was a high patchiness in distribution of fungal OTUs across root systems, indicating that stochastic processes to a large extent structure the fungal communities. However, time since deglaciation had impact on the fungal community structure, as a systematic shift in the community composition was observed along the chronosequence. Ectomycorrhizal basidiomycetes were the dominant fungi in the roots of B. vivipara, when it comes to both number of OTUs and number of sequences.