BackgroundDuring the last 15 years the internal transcribed spacer (ITS) of nuclear DNA has been used as a target for analyzing fungal diversity in environmental samples, and has recently been selected as the standard marker for fungal DNA barcoding. In this study we explored the potential amplification biases that various commonly utilized ITS primers might introduce during amplification of different parts of the ITS region in samples containing mixed templates ('environmental barcoding'). We performed in silico PCR analyses with commonly used primer combinations using various ITS datasets obtained from public databases as templates.ResultsSome of the ITS primers, such as ITS1-F, were hampered with a high proportion of mismatches relative to the target sequences, and most of them appeared to introduce taxonomic biases during PCR. Some primers, e.g. ITS1-F, ITS1 and ITS5, were biased towards amplification of basidiomycetes, whereas others, e.g. ITS2, ITS3 and ITS4, were biased towards ascomycetes. The assumed basidiomycete-specific primer ITS4-B only amplified a minor proportion of basidiomycete ITS sequences, even under relaxed PCR conditions. Due to systematic length differences in the ITS2 region as well as the entire ITS, we found that ascomycetes will more easily amplify than basidiomycetes using these regions as targets. This bias can be avoided by using primers amplifying ITS1 only, but this would imply preferential amplification of 'non-dikarya' fungi.ConclusionsWe conclude that ITS primers have to be selected carefully, especially when used for high-throughput sequencing of environmental samples. We suggest that different primer combinations or different parts of the ITS region should be analyzed in parallel, or that alternative ITS primers should be searched for.

Novel high-throughput sequencing methods outperform earlier approaches in terms of resolution and magnitude. They enable identification and relative quantification of community members and offer new insights into fungal community ecology. These methods are currently taking over as the primary tool to assess fungal communities of plant-associated endophytes, pathogens, and mycorrhizal symbionts, as well as free-living saprotrophs.Taking advantage of the collective experience of six research groups, we here review the different stages involved in fungal community analysis, from field sampling via laboratory procedures to bioinformatics and data interpretation. We discuss potential pitfalls, alternatives, and solutions.Highlighted topics are challenges involved in: obtaining representative DNA/RNA samples and replicates that encompass the targeted variation in community composition, selection of marker regions and primers, options for amplification and multiplexing, handling of sequencing errors, and taxonomic identification.Without awareness of methodological biases, limitations of markers, and bioinformatics challenges, large-scale sequencing projects risk yielding artificial results and misleading conclusions.

TTh he e p pl la an nt t c ce el ll l w wa al ll l d de ec co om mp po os si in ng g m ma ac ch hi in ne er ry y u un nd de er rl li ie es s t th he e f fu un nc ct ti io on na al l d di iv ve er rs si it ty y o of f f fo or re es st t f fu un ng gi i T Th he e p pl la an nt t c ce el ll l w wa al ll l d de ec co om mp po os si in ng g m ma ac ch hi in ne er ry y u un nd de er rl li ie es s t th he e f fu un nc ct ti io on na al l d di iv ve er rs si it ty y o of f f fo or re es st t f fu un ng gi i

No abstract

The cryptic lifestyle of most fungi necessitates molecular identification of the guild in environmental studies.Over the past decades, rapid development and affordability of molecular tools have tremendously improved insights of the fungal diversity in all ecosystems and habitats. Yet, in spite of the progress of molecular methods, knowledge about functional properties of the fungal taxa is vague and interpretation of environmental studies in an ecologically meaningful manner remains challenging. In order to facilitate functional assignments and ecological interpretation of environmental studies we introduce a user friendly traits and character database FungalTraits operating at genus and species hypothesis levels. Combining the information from previous efforts such as FUNGuild and Fun Fun together with involvement of expert knowledge, we reannotated 10210 and 151 fungal and Stramenopila genera, respectively. This resulted in a stand-alone spreadsheet dataset covering 17 lifestyle related traits of fungal and Stramenopila genera, designed for rapid functional assignments of environmental studies. In order to assign the trait states to fungal species hypotheses, the scientific community of experts manually categorised and assigned available trait information to 697413 fungal ITS sequences. On the basis of those sequences we were able to summarise trait and host information into 92623 fungal species hypotheses at 1% dissimilarity threshold.

The nuclear ribosomal Internal Transcribed Spacer ITS region is widely used as a DNA metabarcoding marker to characterize the diversity and composition of fungal communities. In amplicon pyrosequencing studies of fungal diversity, one of the spacers ITS1 or ITS2 of the ITS region is normally used. In this methodological study we evaluate the usability of ITS1 vs. ITS2 as a DNA metabarcoding marker for fungi. We analyse three data sets: two comprising ITS1 and ITS2 sequences of known taxonomic affiliations and a third comprising ITS1 and ITS2 environmental amplicon pyrosequencing data. Clustering analyses of sequences with known taxonomy using the bioinformatics pipeline ClustEx revealed that a 97% similarity cut-off represent a reasonable threshold for estimating the number of known species in the data sets for both ITS1 and ITS2. However, no single threshold value worked well for all fungi at the same time within the curated UNITE database, and we found that the Operational Taxonomic Unit (OTU) concept is not easily translated into the level of species because many species are distributed over several clusters. Clustering analyses of the 134 692 ITS1 and ITS2 pyrosequences using a 97% similarity cut-off revealed a high similarity between the two data sets when it comes to taxonomic coverage. Although some groups are under- or unrepresented in the two data sets due to, e.g. primer mismatches, our results indicate that ITS1 and ITS2 to a large extent yield similar results when used as DNA metabarcodes for fungi.

Metabarcoding approaches use total and typically degraded DNA from environmental samples to analyse biotic assemblages and can potentially be carried out for any kinds of organisms in an ecosystem. These analyses rely on specific markers, here called metabarcodes, which should be optimized for taxonomic resolution, minimal bias in amplification of the target organism group and short sequence length. Using bioinformatic tools, we developed metabarcodes for several groups of organisms: fungi, bryophytes, enchytraeids, beetles and birds. The ability of these metabarcodes to amplify the target groups was systematically evaluated by (i) in silico PCRs using all standard sequences in the EMBL public database as templates, (ii) in vitro PCRs of DNA extracts from surface soil samples from a site in Varanger, northern Norway and (iii) in vitro PCRs of DNA extracts from permanently frozen sediment samples of late‐Pleistocene age (∼16 000–50 000 years bp) from two Siberian sites, Duvanny Yar and Main River. Comparison of the results from the in silico PCR with those obtained in vitro showed that the in silico approach offered a reliable estimate of the suitability of a marker. All target groups were detected in the environmental DNA, but we found large variation in the level of detection among the groups and between modern and ancient samples. Success rates for the Pleistocene samples were highest for fungal DNA, whereas bryophyte, beetle and bird sequences could also be retrieved, but to a much lesser degree. The metabarcoding approach has considerable potential for biodiversity screening of modern samples and also as a palaeoecological tool.

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.