Exploring the accuracy of amplicon‐based internal transcribed spacer markers for a fungal community

With the continual improvement in high‐throughput sequencing technology and constant updates to fungal reference databases, the use of amplicon‐based DNA markers as a tool to reveal fungal diversity and composition in various ecosystems has become feasible. However, both primer selection and the experimental procedure require meticulous verification. Here, we computationally and experimentally evaluated the accuracy and specificity of three widely used or newly designed internal transcribed spacer (ITS) primer sets (ITS1F/ITS2, gITS7/ITS4 and 5.8S‐Fun/ITS4‐Fun). In silico evaluation revealed that primer coverage varied at different taxonomic levels due to differences in degeneracy and the location of primer sets. Using even and staggered mock community standards, we identified different proportions of chimeric and mismatch reads generated by different primer sets, as well as great variation in species abundances, suggesting that primer selection would affect the results of amplicon‐based metabarcoding studies. Choosing proofreading and high‐fidelity polymerase (KAPA HiFi) could significantly reduce the percentage of chimeric and mismatch sequences, further reducing inflation of operational taxonomic units. Moreover, for two types of environmental fungal communities, plant endophytic and soil fungi, it was demonstrated that the three primer sets could not reach a consensus on fungal community composition or diversity, and that primer selection, not experimental treatment, determines observed soil fungal community diversity and composition. Future DNA marker surveys should pay greater attention to potential primer effects and improve the experimental scheme to increase credibility and accuracy.     

Communities of wood-inhabiting fungi in dead pine logs along a geographical gradient in Japan

Fungi are the main agents of coarse woody debris decomposition in forest ecosystems. We examined the associations of environmental variables with fungal community structures in dead pine logs at 12 geographically distant sites using amplicon pyrosequencing of fungal ITS rDNA. A total of 575 operational taxonomic units (OTUs) were identified based on clustering at 97% similarity. Among the known fungal ecological groups, saprotrophic fungi generally showed highest frequency of occurrence and were positively associated with mean annual temperature (MAT) and log diameter. Wood decay fungi with unknown decay type were positively associated with pine wilt disease and negatively associated with log diameter. Ordination analysis of the 42 most prevalent OTUs showed that MAT and annual precipitation significantly explained the observed fungal community structure. These results suggested that climate conditions and site history differentially effect structure fungal communities in pine logs among different ecological groups.     

Employing 454 amplicon pyrosequencing to reveal intragenomic divergence in the internal transcribed spacer rDNA region in fungi

The rDNA internal transcribed spacer (ITS) region has been accepted as a DNA barcoding marker for fungi and is widely used in phylogenetic studies; however, intragenomic ITS variability has been observed in a broad range of taxa, including prokaryotes, plants, animals, and fungi, and this variability has the potential to inflate species richness estimates in molecular investigations of environmental samples. In this study 454 amplicon pyrosequencing of the ITS1 region was applied to 99 phylogenetically diverse axenic single‐spore cultures of fungi (Dikarya: Ascomycota and Basidiomycota) to investigate levels of intragenomic variation. Three species (one Basidiomycota and two Ascomycota), in addition to a positive control species known to contain ITS paralogs, displayed levels of molecular variation indicative of intragenomic variation; taxon inflation due to presumed intragenomic variation was ≈9%. Intragenomic variability in the ITS region appears to be widespread but relatively rare in fungi (≈3–5% of species investigated in this study), suggesting this problem may have minor impacts on species richness estimates relative to PCR and/or pyrosequencing errors. Our results indicate that 454 amplicon pyrosequencing represents a powerful tool for investigating levels of ITS intragenomic variability across taxa, which may be valuable for better understanding the fundamental mechanisms underlying concerted evolution of repetitive DNA regions.     

ITS2 metabarcoding analysis complements lichen mycobiome diversity data

Lichen thalli harbor complex fungal communities (mycobiomes) of species with divergent trophic and ecological strategies. The complexity and diversity of lichen mycobiomes are still largely unknown, despite surveys combining culture-based methods and high-throughput sequencing (HTS). The results of such surveys are strongly influenced by the barcode locus chosen, its sensitivity in discriminating taxa, and the depth to which public sequence repositories cover the phylogenetic spectrum of fungi. Here, we use HTS of the internal transcribed spacer 2 (ITS2) to assess the taxonomic composition and diversity of a well-characterized, alpine rock lichen community that includes thalli symptomatically infected by lichenicolous fungi as well as asymptomatic thalli. Taxa belonging to the order Chaetothyriales are the major components of the observed lichen mycobiomes. We predict sequences representative of lichenicolous fungi characterized morphologically and assess their asymptomatic presence in lichen thalli. We demonstrated the limitations of metabarcoding in fungi and show how the estimation of species diversity widely differs when ITS1 or ITS2 are used as barcode, and particularly biases the detection of Basidiomycota. The complementary analysis of both ITS1 and ITS2 loci is therefore required to reliably estimate the diversity of lichen mycobiomes.     

Comparison of ITS1 and ITS2 rDNA in 454 sequencing of hyperdiverse fungal communities

Next-generation amplicon sequencing is a powerful tool in ecological studies of fungi. Technological development suggests that short fragment high-throughput techniques, e.g. Illumina, will gain importance in fungal community analyses. Thus there is a need for short (250 bp) and informative molecular identifiers. Here we compared ITS1 vs. ITS2 rDNA using empirical data from a study of hyperdiverse leaf-associated fungal communities. Our results suggest that ITS2 may be more variable and recovers more of the molecular diversity. We confirm an earlier in silico study showing that ITS1 and ITS2 yielded somewhat different taxonomic community compositions when blasted against public databases. However, we demonstrate that both ITS1 and ITS2 reveal similar patterns in community structure when analyzed in a community ecology context.     

Clustering of fungal community internal transcribed spacer sequence data obscures taxonomic diversity

Next‐generation DNA sequencing has enabled a rapid expansion in the size of molecular fungal ecology studies employing the nuclear internal transcribed spacer (ITS) region. Many sequence‐processing pipelines and protocols require sequence clustering to generate operational taxonomic units (OTUs) based on sequence similarity as a step to reduce total data quantity and complexity prior to taxonomic assignment. However, the consequences of ITS sequence clustering in regard to sample taxonomic coverage have not been carefully examined. Here we demonstrate that typically used clustering thresholds for fungal ITS sequences result in statistically significant losses in taxonomic coverage. Analyses using environmentally derived fungal sequences indicated an average of 3.1% of species went undetected (P < 0.05) if the sequences were denoised and clustered at a 97% threshold prior to taxonomic assignment. Additionally, an in silico analysis using a reference fungal ITS database suggested that approximately 25% of species went undetected if the sequences were clustered prior to taxonomic assignment. Finally, analysis of sequences derived from pure‐cultured fungal isolates of known identity indicated sequence denoising and clustering were not critical in improving identification accuracy.     

Patterns of Symbiodinium (Dinophyceae) diversity and assemblages among diverse hosts and the coral reef environment of Lizard Island,Australia

Large‐scale environmental disturbances may impact both partners in coral host–Symbiodinium systems. Elucidation of the assembly patterns in such complex and interdependent communities may enable better prediction of environmental impacts across coral reef ecosystems. In this study, we investigated how the community composition and diversity of dinoflagellate symbionts in the genus Symbiodinium were distributed among 12 host species from six taxonomic orders (Actinaria, Alcyonacea, Miliolida, Porifera, Rhizostoma, Scleractinia) and in the reef water and sediments at Lizard Island, Great Barrier Reef before the 3rd Global Coral Bleaching Event. 454 pyrosequencing of the ITS2 region of Symbiodinium yielded 83 operational taxonomic units (OTUs) at a 97% similarity cut‐off. Approximately half of the Symbiodinium OTUs from reef water or sediments were also present in symbio. OTUs belonged to six clades (A‐D, F‐G), but community structure was uneven. The two most abundant OTUs (100% matches to types C1 and A3) comprised 91% of reads and OTU C1 was shared by all species. However, sequence‐based analysis of these dominant OTUs revealed host species specificity, suggesting that genetic similarity cut‐offs of Symbiodinium ITS2 data sets need careful evaluation. Of the less abundant OTUs, roughly half occurred at only one site or in one species and the background Symbiodinium communities were distinct between individual samples. We conclude that sampling multiple host taxa with differing life history traits will be critical to fully understand the symbiont diversity of a given system and to predict coral ecosystem responses to environmental change and disturbance considering the differential stress response of the taxa within.     

Intragenomic variation in the ITS rDNA region obscures phylogenetic relationships and inflates estimates of operational taxonomic units in genus Laetiporus

Regions of rDNA are commonly used to infer phylogenetic relationships among fungal species and as DNA barcodes for identification. These regions occur in large tandem arrays, and concerted evolution is believed to reduce intragenomic variation among copies within these arrays, although some variation still might exist. Phylogenetic studies typically use consensus sequencing, which effectively conceals most intragenomic variation, but cloned sequences containing intragenomic variation are becoming prevalent in DNA databases. To understand effects of using cloned rDNA sequences in phylogenetic analyses we amplified and cloned the ITS region from pure cultures of six Laetiporus species and one Wolfiporia species (Basidiomycota, Polyporales). An average of 66 clones were selected randomly and sequenced from 21 cultures, producing a total of 1399 interpretable sequences. Significant variation (≥ 5% variation in sequence similarity) was observed among ITS copies within six cultures from three species clades (L. cincinnatus, L. sp. clade J, and Wolfiporia dilatohypha) and phylogenetic analyses with the cloned sequences produced different trees relative to analyses with consensus sequences. Cloned sequences from L. cincinnatus fell into more than one species clade and numerous cloned L. cincinnatus sequences fell into entirely new clades, which if analyzed on their own most likely would be recognized as "undescribed" or "novel" taxa. The use of a 95% cut off for defining operational taxonomic units (OTUs) produced seven Laetiporus OTUs with consensus ITS sequences and 20 OTUs with cloned ITS sequences. The use of cloned rDNA sequences might be problematic in fungal phylogenetic analyses, as well as in fungal bar-coding initiatives and efforts to detect fungal pathogens in environmental samples.     

An open source chimera checker for the fungal ITS region

The internal transcribed spacer (ITS) region of the nuclear ribosomal repeat unit holds a central position in the pursuit of the taxonomic affiliation of fungi recovered through environmental sampling. Newly generated fungal ITS sequences are typically compared against the International Nucleotide Sequence Databases for a species or genus name using the sequence similarity software suite blast . Such searches are not without complications however, and one of them is the presence of chimeric entries among the query or reference sequences. Chimeras are artificial sequences, generated unintentionally during the polymerase chain reaction step, that feature sequence data from two (or possibly more) distinct species. Available software solutions for chimera control do not readily target the fungal ITS region, but the present study introduces a blast -based open source software package (available at http://www.emerencia.org/chimerachecker.html ) to examine newly generated fungal ITS sequences for the presence of potentially chimeric elements in batch mode. We used the software package on a random set of 12 300 environmental fungal ITS sequences in the public sequence databases and found 1.5% of the entries to be chimeric at the ordinal level after manual verification of the results. The proportion of chimeras in the sequence databases can be hypothesized to increase as emerging sequencing technologies drawing from pooled DNA samples are becoming important tools in molecular ecology research.     

Multiple markers pyrosequencing reveals highly diverse and host-specific fungal communities on the mangrove trees Avicennia marina and Rhizophora stylosa

Fungi are important actors in ecological processes and trophic webs in mangroves. Although saprophytic fungi occurring in the intertidal part of mangrove have been well studied, little is known about the diversity and structure of the fungal communities in this ecosystem or about the importance of functional groups like pathogens and mutualists. Using tag-encoded 454 pyrosequencing of the ITS1, ITS2, nu-ssu-V5 and nu-ssu-V7 regions, we studied and compared the fungal communities found on the marine and aerial parts of Avicennia marina and Rhizophora stylosa trees in a mangrove in New Caledonia. A total of 209,544 reads were analysed, corresponding to several thousand molecular operational taxonomic units (OTU). There is a marked zonation in the species distribution, with most of the OTU being found specifically in one of the microhabitat studied. Ascomycetes are the dominant phylum (82%), Basidiomycetes are very rare (3%), and 15% of the sequences correspond to unknown taxa. Our results indicate that host specificity is a key factor in the distribution of the highly diverse fungal communities, in both the aerial and intertidal parts of the trees. This study also validates the usefulness of multiple markers in tag-encoded pyrosequencing to consolidate and refine the assessment of the taxonomic diversity.     

ITS1 metabarcoding highlights low specificity of lichen mycobiomes at a local scale

As self‐supporting and long‐living symbiotic structures, lichens provide a habitat for many other organisms beside the traditionally considered lichen symbionts—the myco‐ and the photobionts. The lichen‐inhabiting fungi either develop diagnostic phenotypes or occur asymptomatically. Because the degree of specificity towards the lichen host is poorly known, we studied the diversity of these fungi among neighbouring lichens on rocks in an alpine habitat. Using a sequencing metabarcoding approach, we show that lichen mycobiomes clearly reflect the overlap of multiple ecological sets of taxa, which differ in their trophic association with lichen thalli. The lack of specificity to the lichen mycobiome is further supported by the lack of community structure observed using clustering and ordination methods. The communities encountered across samples largely result from the subsampling of a shared species pool, in which we identify three major ecological components: (i) a generalist environmental pool, (ii) a lichenicolous/endolichenic pool and (iii) a pool of transient species. These taxa majorly belong to the fungal classes Dothideomycetes, Eurotiomycetes and Tremellomycetes with close relatives in adjacent ecological niches. We found no significant evidence that the phenotypically recognized lichenicolous fungi influence the occurrence of the other asymptomatic fungi in the host thalli. We claim that lichens work as suboptimal habitats or as a complex spore and mycelium bank, which modulate and allow the regeneration of local fungal communities. By performing an approach that minimizes ambiguities in the taxonomic assignments of fungi, we present how lichen mycobiomes are also suitable targets for improving bioinformatic analyses of fungal metabarcoding.     

Does peptide-nucleic acid (PNA) clamping of host plant DNA benefit ITS1 amplicon-based characterization of the fungal endophyte community?

Fungal endophyte community amplicon sequencing can lose a significant number of informative reads due to host-plant co-amplification. Blocking of plant-specific sequences with peptide nucleic acid (PNA) clamps has been shown to improve metrics of detected microbial diversity in studies targeting 16S and 18S regions of rRNA genes. However, PNA clamping has not been applied to the plant ITS region of rRNA gene – a widely accepted fungal marker. By applying PNA clamping technique to ITS amplicon sequencing of the endophytic fungal community of elderberry this study shows that PNA clamping significantly reduces host-plant co-amplification with the universal ITS1/ITS4 primer set. However, PNA clamping in combination with the discriminatory ITS1F/ITS2 primer set did not improve the metrics of fungal endophyte community ITS amplicon Illumina sequencing. This study shows that PNA clamping does not add practical benefit to taxonomic profiling of plant-associated fungal communities if the primers are already specific enough to exclude amplification of host DNA.     

Fungal sporocarps house diverse and host-specific communities of fungicolous fungi

Sporocarps (fruit bodies) are the sexual reproductive stage in the life cycle of many fungi. They are highly nutritious and consequently vulnerable to grazing by birds and small mammals, and invertebrates, and can be infected by microbial and fungal parasites and pathogens. The complexity of communities thriving inside sporocarps is largely unknown. In this study, we revealed the diversity, taxonomic composition and host preference of fungicolous fungi (i.e., fungi that feed on other fungi) in sporocarps. We carried out DNA metabarcoding of the ITS2 region from 176 sporocarps of 11 wood-decay fungal host species, all collected within a forest in northeast Finland. We assessed the influence of sporocarp traits, such as lifespan, morphology and size, on the fungicolous fungal community. The level of colonisation by fungicolous fungi, measured as the proportion of non-host ITS2 reads, varied between 2.8–39.8% across the 11 host species and was largely dominated by Ascomycota. Host species was the major determinant of the community composition and diversity of fungicolous fungi, suggesting that host adaptation is important for many fungicolous fungi. Furthermore, the alpha diversity was consistently higher in short-lived and resupinate sporocarps compared to long-lived and pileate ones, perhaps due to a more hostile environment for fungal growth in the latter too. The fungicolous fungi represented numerous lineages in the fungal tree of life, among which a significant portion was poorly represented with reference sequences in databases.Subject terms: Microbial ecology, Molecular evolution     

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.     

High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples

The kingdom Fungi is estimated to include 1.5 million or more species, playing key roles as decomposers, mutualists, and parasites in every biome on the earth. To comprehensively understand the diversity and ecology of this huge kingdom, DNA barcoding targeting the internal transcribed spacer (ITS) region of the nuclear ribosomal repeat has been regarded as a prerequisite procedure. By extensively surveying ITS sequences in public databases, we designed new ITS primers with improved coverage across diverse taxonomic groups of fungi compared to existing primers. An in silico analysis based on public sequence databases indicated that the newly designed primers matched 99% of ascomycete and basidiomycete ITS taxa (species, subspecies or varieties), causing little taxonomic bias toward either fungal group. Two of the newly designed primers could inhibit the amplification of plant sequences and would enable the selective investigation of fungal communities in mycorrhizal associations, soil, and other types of environmental samples. Optimal PCR conditions for the primers were explored in an in vitro investigation. The new primers developed in this study will provide a basis for ecological studies on the diversity and community structures of fungi in the era of massive DNA sequencing.     

Optimal sequence similarity thresholds for clustering of molecular operational taxonomic units in DNA metabarcoding studies

Clustering approaches are pivotal to handle the many sequence variants obtained in DNA metabarcoding data sets, and therefore they have become a key step of metabarcoding analysis pipelines. Clustering often relies on a sequence similarity threshold to gather sequences into molecular operational taxonomic units (MOTUs), each of which ideally represents a homogeneous taxonomic entity (e.g., a species or a genus). However, the choice of the clustering threshold is rarely justified, and its impact on MOTU over-splitting or over-merging even less tested. Here, we evaluated clustering threshold values for several metabarcoding markers under different criteria: limitation of MOTU over-merging, limitation of MOTU over-splitting, and trade-off between over-merging and over-splitting. We extracted sequences from a public database for nine markers, ranging from generalist markers targeting Bacteria or Eukaryota, to more specific markers targeting a class or a subclass (e.g., Insecta, Oligochaeta). Based on the distributions of pairwise sequence similarities within species and within genera, and on the rates of over-splitting and over-merging across different clustering thresholds, we were able to propose threshold values minimizing the risk of over-splitting, that of over-merging, or offering a trade-off between the two risks. For generalist markers, high similarity thresholds (0.96–0.99) are generally appropriate, while more specific markers require lower values (0.85–0.96). These results do not support the use of a fixed clustering threshold. Instead, we advocate careful examination of the most appropriate threshold based on the research objectives, the potential costs of over-splitting and over-merging, and the features of the studied markers.     

Nutrient enrichment increased species richness of leaf litter fungal assemblages in a tropical forest

Microbial communities play a major role in terrestrial ecosystem functioning, but the determinates of their diversity and functional interactions are not well known. In this study, we explored leaf litter fungal diversity in a diverse Panama lowland tropical forest in which a replicated factorial N, P, K and micronutrient fertilization experiment of 40 × 40 m plots had been ongoing for nine years. We extracted DNA from leaf litter samples and used fungal‐specific amplification and a 454 pyrosequencing approach to sequence two loci, the nuclear ribosomal internal transcribed spacer (ITS) region and the nuclear ribosomal large subunit (LSU) D1 region. Using a 95% sequence similarity threshold for ITS1 spacer recovered a total of 2523 OTUs, and the number of unique ITS1 OTUs per 0.5–1.0 g leaf litter sample ranged from 55 to 177. Ascomycota were the dominant phylum among the leaf litter fungi (71% of the OTUs), followed by Basidiomycota (26% of the OTUs). In contrast to our expectations based on temperate ecosystems, long‐term addition of nutrients increased, rather than decreased, species richness relative to controls. Effect of individual nutrients was more subtle and seen primarily as changes in community compositions especially at lower taxonomic levels, rather than as significant changes in species richness. For example, plots receiving P tended to show a greater similarity in community composition compared to the other nutrient treatments, the +PK, +NK and +NPK plots appeared to be more dominated by the Nectriaceae than other treatments, and indicator species for particular nutrient combinations were identified.     

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.     

Algal and Fungal Diversity in Antarctic Lichens

The composition of lichen ecosystems except mycobiont and photobiont has not been evaluated intensively. In addition, recent studies to identify algal genotypes have raised questions about the specific relationship between mycobiont and photobiont. In the current study, we analyzed algal and fungal community structures in lichen species from King George Island, Antarctica, by pyrosequencing of eukaryotic large subunit (LSU) and algal internal transcribed spacer (ITS) domains of the nuclear rRNA gene. The sequencing results of LSU and ITS regions indicated that each lichen thallus contained diverse algal species. The major algal operational taxonomic unit (OTU) defined at a 99% similarity cutoff of LSU sequences accounted for 78.7–100% of the total algal community in each sample. In several cases, the major OTUs defined by LSU sequences were represented by two closely related OTUs defined by 98% sequence similarity of ITS domain. The results of LSU sequences indicated that lichen‐associated fungi belonged to the Arthoniomycetes, Eurotiomycetes, Lecanoromycetes, Leotiomycetes, and Sordariomycetes of the Ascomycota, and Tremellomycetes and Cystobasidiomycetes of the Basidiomycota. The composition of major photobiont species and lichen‐associated fungal community were mostly related to the mycobiont species. The contribution of growth forms or substrates on composition of photobiont and lichen‐associated fungi was not evident.     

ITS all right mama: investigating the formation of chimeric sequences in the ITS2 region by DNA metabarcoding analyses of fungal mock communities of different complexities

The formation of chimeric sequences can create significant methodological bias in PCR‐based DNA metabarcoding analyses. During mixed‐template amplification of barcoding regions, chimera formation is frequent and well documented. However, profiling of fungal communities typically uses the more variable rDNA region ITS. Due to a larger research community, tools for chimera detection have been developed mainly for the 16S/18S markers. However, these tools are widely applied to the ITS region without verification of their performance. We examined the rate of chimera formation during amplification and 454 sequencing of the ITS2 region from fungal mock communities of different complexities. We evaluated the chimera detecting ability of two common chimera‐checking algorithms: perseus and uchime . Large proportions of the chimeras reported were false positives. No false negatives were found in the data set. Verified chimeras accounted for only 0.2% of the total ITS2 reads, which is considerably less than what is typically reported in 16S and 18S metabarcoding analyses. Verified chimeric ‘parent sequences’ had significantly higher per cent identity to one another than to random members of the mock communities. Community complexity increased the rate of chimera formation. GC content was higher around the verified chimeric break points, potentially facilitating chimera formation through base pair mismatching in the neighbouring regions of high similarity in the chimeric region. We conclude that the hypervariable nature of the ITS region seems to buffer the rate of chimera formation in comparison with other, less variable barcoding regions, due to shorter regions of high sequence similarity.