Approaching the taxonomic affiliation of unidentified sequences in public databases--an example from the mycorrhizal fungi

Lentinula madagasikarensis sp. nov., a relative of shiitake mushrooms from Madagascar

We describe the first species of Lentinula from Africa, Lentinula madagasikarensis sp. nov. The new taxon, which was collected from central Madagascar, is strikingly similar to L. edodes, the shiitake mushroom. A BLAST search using ITS sequences from L. madagasikarensis as the query retrieves a mix of Lentinula, Gymnopus, Marasmiellus, and other members of Omphalotaceae as the top hits. A 28S phylogeny of the Omphalotaceae confirms placement of L. madagasikarensis within Lentinula. An ITS phylogeny places L. madagasikarensis as the sister group of L. aciculospora, which is a neotropical species. Lentinula madagasikarensis is characterized by robust basidiomata with vinaceous pilei, prominent floccose scales near the pileus margin, florets of sphaeropedunculate cheilocystidia, and subcylindrical basidiospores. This report constitutes a 4 000-mile, trans-oceanic range extension for Lentinula.

Species diversity, phylogeny, endemism and geography of the truffle genus Tuber in China based on morphological and molecular data

The genus Tuber (Tuberaceae, Pezizales) is an important fungal group of Ascomycota both economically and ecologically. However, the species diversity, phylogenetic relationships, and geographic distribution of Tuber species in China remains poorly understood, primarily because descriptions of many new species relied heavily on morphological features with molecular data either not sought or ignored. The misapplication of European and North American names further added to confusion regarding the taxonomy of Chinese Tuber species. In this study, we examined more than 1000 specimens from China, and performed a comprehensive phylogenetic analysis for Chinese Tuber species using ITS sequences and multilocus sequence data. To infer the phylogeny of Chinese Tuber spp., 11 molecular datasets were assembled, including a concatenated internal transcribed spacers of the nuc rDNA (ITS), nuc rDNA 28S subunit (LSU), translation elongation factor 1-alpha (tef1-α), and RNA polymerase II subunit (rpb2) dataset as well as 10 ITS datasets (totally including 1435 sequences from 828 collections with 597 newly generated sequences, and 168 sequences from the types of 63 species). Our phylogenetic tree based on a concatenated multilocus dataset revealed that all Chinese Tuber species nested in nine phylogenetic clades (phylogroups), including Aestivum, Excavatum, Latisporum, Macrosporum, Maculatum, Melanosporum, Puberulum, Rufum and Turmericum. Of these, five phylogroups (Macrosporum, Maculatum, Melanosporum, Puberulum and Rufum) are shared across the continents of Asia, Europe and North America; two phylogroups (Aestivum and Excavatum) are shared by Europe and Asia; and the phylogroups Turmericum and Latisporum are endemic only to Asia. Phylogenetic trees based on 10 ITS datasets confirmed the presence of at least 82 phylogenetic species in China. Of these, 53 are identified as known species, including three new records for China, and 25 species are identified as new to science. Of the new species, nine are described and illustrated in this paper, and the others remain unnamed due to the paucity or absence of ascomatal materials. Accordingly, the confirmed, excluded and doubtful Tuber species in China are discussed. Tuber species showed high endemism. Of the 82 phylogenetic species found in China, 68 species occur only in China, six species are also found in other regions in Asia, and only eight species (T. anniae, T. excelsum-reticulatum, T. formosanum, T. maculatum, T. wenchuanense, Tuber sp. CHN-3, Tuber sp. CHN-10 and Tuber sp. CHN-11) are shared with other continents. Most Tuber species have a small and limited distribution in China, but a few, such as T. formosanum and T. parvomurphium, are widely distributed across China. Some phylogenetically closely related species, such as T. liaotongense and T. subglobosum, as well as T. xuanhuaense and T. lijiangense, show a pattern of allopatric distribution.

Long- and short-read metabarcoding technologies reveal similar spatiotemporal structures in fungal communities

Fungi form diverse communities and play essential roles in many terrestrial ecosystems, yet there are methodological challenges in taxonomic and phylogenetic placement of fungi from environmental sequences. To address such challenges, we investigated spatiotemporal structure of a fungal community using soil metabarcoding with four different sequencing strategies: short-amplicon sequencing of the ITS2 region (300-400 bp) with Illumina MiSeq, Ion Torrent Ion S5 and PacBio RS II, all from the same PCR library, as well as long-amplicon sequencing of the full ITS and partial LSU regions (1200-1600 bp) with PacBio RS II. Resulting community structure and diversity depended more on statistical method than sequencing technology. The use of long-amplicon sequencing enables construction of a phylogenetic tree from metabarcoding reads, which facilitates taxonomic identification of sequences. However, long reads present issues for denoising algorithms in diverse communities. We present a solution that splits the reads into shorter homologous regions prior to denoising, and then reconstructs the full denoised reads. In the choice between short and long amplicons, we suggest a hybrid approach using short amplicons for sampling breadth and depth, and long amplicons to characterize the local species pool for improved identification and phylogenetic analyses.

Glacier retreat in the High Arctic: opportunity or threat for ectomycorrhizal diversity?

Climate change causes Arctic glaciers to retreat faster, exposing new areas for colonization. Several pioneer plants likely to colonize recent deglaciated, nutrient-poor areas depend on fungal partners for successful establishment. Little is known about general patterns or characteristics of facilitating fungal pioneers and how they vary with regional climate in the Arctic. The High Arctic Archipelago Svalbard represents an excellent study system to address these questions, as glaciers cover ∼60% of the land surface and recent estimations suggest at least 7% reduction of glacier area since 1960s. Roots of two ectomycorrhizal (ECM) plants (Salix polaris and Bistorta vivipara) were sampled in eight glacier forelands. Associated ECM fungi were assessed using DNA metabarcoding. About 25% of the diversity was unknown at family level, indicating presence of undescribed species. Seven genera dominated based on richness and abundance, but their relative importance varied with local factors. The genus Geopora showed surprisingly high richness and abundance, particularly in dry, nutrient-poor forelands. Such forelands will diminish along with increasing temperature and precipitation, and faster succession. Our results support a taxonomical shift in pioneer ECM diversity with climate change, and we are likely to lose unknown fungal diversity, without knowing their identity or ecological importance.

Mycobiome diversity: high-throughput sequencing and identification of fungi

Fungi are major ecological players in both terrestrial and aquatic environments by cycling organic matter and channelling nutrients across trophic levels. High-throughput sequencing (HTS) studies of fungal communities are redrawing the map of the fungal kingdom by hinting at its enormous - and largely uncharted - taxonomic and functional diversity. However, HTS approaches come with a range of pitfalls and potential biases, cautioning against unwary application and interpretation of HTS technologies and results. In this Review, we provide an overview and practical recommendations for aspects of HTS studies ranging from sampling and laboratory practices to data processing and analysis. We also discuss upcoming trends and techniques in the field and summarize recent and noteworthy results from HTS studies targeting fungal communities and guilds. Our Review highlights the need for reproducibility and public data availability in the study of fungal communities. If the associated challenges and conceptual barriers are overcome, HTS offers immense possibilities in mycology and elsewhere.

Over 2.5 million COI sequences in GenBank and growing

The increasing popularity of cytochrome c oxidase subunit 1 (COI) DNA metabarcoding warrants a careful look at the underlying reference databases used to make high-throughput taxonomic assignments. The objectives of this study are to document trends and assess the future usability of COI records for metabarcode identification. The number of COI records deposited to the NCBI nucleotide database has increased by a geometric average of 51% per year, from 8,137 records deposited in 2003 to a cumulative total of ~ 2.5 million by the end of 2017. About half of these records are fully identified to the species rank, 92% are at least 500 bp in length, 74% have a country annotation, and 51% have latitude-longitude annotations. To ensure the future usability of COI records in GenBank we suggest: 1) Improving the geographic representation of COI records, 2) Improving the cross-referencing of COI records in the Barcode of Life Data System and GenBank to facilitate consolidation and incorporation into existing bioinformatic pipelines, 3) Adherence to the minimum information about a marker gene sequence guidelines, and 4) Integrating metabarcodes from eDNA and mixed community studies with existing reference sequences. The growth of COI reference records over the past 15 years has been substantial and is likely to be a resource across many fields for years to come.

Formal description of sequence-based voucherless Fungi: promises and pitfalls, and how to resolve them

There is urgent need for a formal nomenclature of sequence-based, voucherless Fungi, given that environmental sequencing has accumulated more than one billion fungal ITS reads in the Sequence Read Archive, about 1,000 times as many as fungal ITS sequences in GenBank. These unnamed Fungi could help to bridge the gap between 115,000 to 140,000 currently accepted and 2.2 to 3.8 million predicted species, a gap that cannot realistically be filled using specimen or culture-based inventories. The Code never aimed at placing restrictions on the nature of characters chosen for taxonomy, and the requirement for physical types is now becoming a constraint on the advancement of science. We elaborate on the promises and pitfalls of sequence-based nomenclature and provide potential solutions to major concerns of the mycological community. Types of sequence-based taxa, which by default lack a physical specimen or culture, could be designated in four alternative ways: (1) the underlying sample ('bag' type), (2) the DNA extract, (3) fluorescent in situ hybridization (FISH), or (4) the type sequence itself. Only (4) would require changes to the Code and the latter would be the most straightforward approach, complying with three of the five principal functions of types better than physical specimens. A fifth way, representation of the sequence in an illustration, has been ruled as unacceptable in the Code. Potential flaws in sequence data are analogous to flaws in physical types, and artifacts are manageable if a stringent analytical approach is applied. Conceptual errors such as homoplasy, intragenomic variation, gene duplication, hybridization, and horizontal gene transfer, apply to all molecular approaches and cannot be used as a specific argument against sequence-based nomenclature. The potential impact of these phenomena is manageable, as phylogenetic species delimitation has worked satisfactorily in Fungi. The most serious shortcoming of sequence-based nomenclature is the likelihood of parallel classifications, either by describing taxa that already have names based on physical types, or by using different markers to delimit species within the same lineage. The probability of inadvertently establishing sequence-based species that have names available is between 20.4 % and 1.5 % depending on the number of globally predicted fungal species. This compares favourably to a historical error rate of about 30 % based on physical types, and this rate could be reduced to practically zero by adding specific provisions to this approach in the Code. To avoid parallel classifications based on different markers, sequence-based nomenclature should be limited to a single marker, preferably the fungal ITS barcoding marker; this is possible since sequence-based nomenclature does not aim at accurate species delimitation but at naming lineages to generate a reference database, independent of whether these lineages represent species, closely related species complexes, or infraspecies. We argue that clustering methods are inappropriate for sequence-based nomenclature; this approach must instead use phylogenetic methods based on multiple alignments, combined with quantitative species recognition methods. We outline strategies to obtain higher-level phylogenies for ITS-based, voucherless species, including phylogenetic binning, 'hijacking' species delimitation methods, and temporal banding. We conclude that voucherless, sequence-based nomenclature is not a threat to specimen and culture-based fungal taxonomy, but a complementary approach capable of substantially closing the gap between known and predicted fungal diversity, an approach that requires careful work and high skill levels.

Scaling up: A guide to high-throughput genomic approaches for biodiversity analysis

The purpose of this review is to present the most common and emerging DNA-based methods used to generate data for biodiversity and biomonitoring studies. As environmental assessment and monitoring programmes may require biodiversity information at multiple levels, we pay particular attention to the DNA metabarcoding method and discuss a number of bioinformatic tools and considerations for producing DNA-based indicators using operational taxonomic units (OTUs), taxa at a variety of ranks and community composition. By developing the capacity to harness the advantages provided by the newest technologies, investigators can "scale up" by increasing the number of samples and replicates processed, the frequency of sampling over time and space, and even the depth of sampling such as by sequencing more reads per sample or more markers per sample. The ability to scale up is made possible by the reduced hands-on time and cost per sample provided by the newest kits, platforms and software tools. Results gleaned from broad-scale monitoring will provide opportunities to address key scientific questions linked to biodiversity and its dynamics across time and space as well as being more relevant for policymakers, enabling science-based decision-making, and provide a greater socio-economic impact. As genomic approaches are continually evolving, we provide this guide to methods used in biodiversity genomics.

New light on names and naming of dark taxa

A growing proportion of fungal species and lineages are known only from sequence data and cannot be linked to any physical specimen or resolved taxonomic name. Such fungi are often referred to as “dark taxa” or “dark matter fungi”. As they lack a taxonomic identity in the form of a name, they are regularly ignored in many important contexts, for example in legalisation and species counts. It is therefore very urgent to find a system to also deal with these fungi. Here, issues relating to the taxonomy and nomenclature of dark taxa are discussed and a number of questions that the mycological community needs to consider before deciding on what system/s to implement are highlighted.

Ribosomal DNA and Plastid Markers Used to Sample Fungal and Plant Communities from Wetland Soils Reveals Complementary Biotas

Though the use of metagenomic methods to sample below-ground fungal communities is common, the use of similar methods to sample plants from their underground structures is not. In this study we use high throughput sequencing of the ribulose-bisphosphate carboxylase large subunit (rbcL) plastid marker to study the plant community as well as the internal transcribed spacer and large subunit ribosomal DNA (rDNA) markers to investigate the fungal community from two wetland sites. Observed community richness and composition varied by marker. The two rDNA markers detected complementary sets of fungal taxa and total fungal composition clustered according to primer rather than by site. The composition of the most abundant plants, however, clustered according to sites as expected. We suggest that future studies consider using multiple genetic markers, ideally generated from different primer sets, to detect a more taxonomically diverse suite of taxa compared with what can be detected by any single marker alone. Conclusions drawn from the presence of even the most frequently observed taxa should be made with caution without corroborating lines of evidence.

Ectomycorrhizal fungi mediate indirect effects of a bark beetle outbreak on secondary chemistry and establishment of pine seedlings

Dendroctonus ponderosae has killed millions of Pinus contorta in western North America with subsequent effects on stand conditions, including changes in light intensity, needle deposition, and the composition of fungal community mutualists, namely ectomycorrhizal fungi. It is unknown whether these changes in stand conditions will have cascading consequences for the next generation of pine seedlings. To test for transgenerational cascades on pine seedlings, we tested the effects of fungal inoculum origin (beetle-killed or undisturbed stands), light intensity and litter (origin and presence) on seedling secondary chemistry and growth in a glasshouse. We also tracked survival of seedlings over two growing seasons in the same stands from which fungi and litter were collected. Fungal communities differed by inoculum origin. Seedlings grown with fungi collected from beetle-killed stands had lower monoterpene concentrations and fewer monoterpene compounds present compared with seedlings grown with fungi collected from undisturbed stands. Litter affected neither monoterpenes nor seedling growth. Seedling survival in the field was lower in beetle-killed than in undisturbed stands. We demonstrate that stand mortality caused by prior beetle attacks of mature pines have cascading effects on seedling secondary chemistry, growth and survival, probably mediated through effects on below-ground mutualisms.

Multi-locus phylogeny of lethal amanitas: implications for species diversity and historical biogeography

BACKGROUND

Lethal amanitas (Amanita section Phalloideae) are a group of wild, fatal mushrooms causing many poisoning cases worldwide. However, the diversity and evolutionary history of these lethal mushrooms remain poorly known due to the limited sampling and insufficient gene fragments employed for phylogenetic analyses. In this study, five gene loci (nrLSU, ITS, rpb2, ef1-α and β-tubulin) with a widely geographic sampling from East and South Asia, Europe, North and Central America, South Africa and Australia were analysed with maximum-likelihood, maximum-parsimony and Bayesian inference methods. Biochemical analyses were also conducted with intention to detect amatoxins and phalloidin in 14 representative samples.

RESULT

Lethal amanitas were robustly supported to be a monophyletic group after excluding five species that were provisionally defined as lethal amanitas based on morphological studies. In lethal amanitas, 28 phylogenetic species were recognised by integrating molecular phylogenetic analyses with morphological studies, and 14 of them represented putatively new species. The biochemical analyses indicated a single origin of cyclic peptide toxins (amatoxins and phalloidin) within Amanita and suggested that this kind of toxins seemed to be a synapomorphy of lethal amanitas. Molecular dating through BEAST and biogeographic analyses with LAGRANGE and RASP indicated that lethal amanitas most likely originated in the Palaeotropics with the present crown group dated around 64.92 Mya in the early Paleocene, and the East Asia-eastern North America or Eurasia-North America-Central America disjunct distribution patterns were primarily established during the middle Oligocene to Miocene.

CONCLUSION

The cryptic diversity found in this study indicates that the species diversity of lethal amanitas is strongly underestimated under the current taxonomy. The intercontinental sister species or sister groups relationships among East Asia and eastern North America or Eurasia-North America-Central America within lethal amanitas are best explained by the diversification model of Palaeotropical origin, dispersal via the Bering Land Bridge, followed by regional vicariance speciation resulting from climate change during the middle Oligocene to the present. These findings indicate the importance of both dispersal and vicariance in shaping the intercontinental distributions of these ectomycorrhizal fungi.

BLASTGrabber: a bioinformatic tool for visualization, analysis and sequence selection of massive BLAST data

Background

Advances in sequencing efficiency have vastly increased the sizes of biological sequence databases, including many thousands of genome-sequenced species. The BLAST algorithm remains the main search engine for retrieving sequence information, and must consequently handle data on an unprecedented scale. This has been possible due to high-performance computers and parallel processing. However, the raw BLAST output from contemporary searches involving thousands of queries becomes ill-suited for direct human processing. Few programs attempt to directly visualize and interpret BLAST output; those that do often provide a mere basic structuring of BLAST data.

Results

Here we present a bioinformatics application named BLASTGrabber suitable for high-throughput sequencing analysis. BLASTGrabber, being implemented as a Java application, is OS-independent and includes a user friendly graphical user interface. Text or XML-formatted BLAST output files can be directly imported, displayed and categorized based on BLAST statistics. Query names and FASTA headers can be analysed by text-mining. In addition to visualizing sequence alignments, BLAST data can be ordered as an interactive taxonomy tree. All modes of analysis support selection, export and storage of data. A Java interface-based plugin structure facilitates the addition of customized third party functionality.

Conclusion

The BLASTGrabber application introduces new ways of visualizing and analysing massive BLAST output data by integrating taxonomy identification, text mining capabilities and generic multi-dimensional rendering of BLAST hits. The program aims at a non-expert audience in terms of computer skills; the combination of new functionalities makes the program flexible and useful for a broad range of operations.

Tissue storage and primer selection influence pyrosequencing-based inferences of diversity and community composition of endolichenic and endophytic fungi

Next-generation sequencing technologies have provided unprecedented insights into fungal diversity and ecology. However, intrinsic biases and insufficient quality control in next-generation methods can lead to difficult-to-detect errors in estimating fungal community richness, distributions and composition. The aim of this study was to examine how tissue storage prior to DNA extraction, primer design and various quality-control approaches commonly used in 454 amplicon pyrosequencing might influence ecological inferences in studies of endophytic and endolichenic fungi. We first contrast 454 data sets generated contemporaneously from subsets of the same plant and lichen tissues that were stored in CTAB buffer, dried in silica gel or freshly frozen prior to DNA extraction. We show that storage in silica gel markedly limits the recovery of sequence data and yields a small fraction of the diversity observed by the other two methods. Using lichen mycobiont sequences as internal positive controls, we next show that despite careful filtering of raw reads and utilization of current best-practice OTU clustering methods, homopolymer errors in sequences representing rare taxa artificially increased estimates of richness c. 15-fold in a model data set. Third, we show that inferences regarding endolichenic diversity can be improved using a novel primer that reduces amplification of the mycobiont. Together, our results provide a rationale for selecting tissue treatment regimes prior to DNA extraction, demonstrate the efficacy of reducing mycobiont amplification in studies of the fungal microbiomes of lichen thalli and highlight the difficulties in differentiating true information about fungal biodiversity from methodological artefacts.

Rapid and accurate taxonomic classification of insect (class Insecta) cytochrome c oxidase subunit 1 ( COI ) DNA barcode sequences using a naïve Bayesian classifier

Current methods to identify unknown insect (class Insecta) cytochrome c oxidase (COI barcode) sequences often rely on thresholds of distances that can be difficult to define, sequence similarity cut-offs, or monophyly. Some of the most commonly used metagenomic classification methods do not provide a measure of confidence for the taxonomic assignments they provide. The aim of this study was to use a naïve Bayesian classifier (Wang et al.Applied and Environmental Microbiology, 2007; 73: 5261) to automate taxonomic assignments for large batches of insect COI sequences such as data obtained from high-throughput environmental sequencing. This method provides rank-flexible taxonomic assignments with an associated bootstrap support value, and it is faster than the blast-based methods commonly used in environmental sequence surveys. We have developed and rigorously tested the performance of three different training sets using leave-one-out cross-validation, two field data sets, and targeted testing of Lepidoptera, Diptera and Mantodea sequences obtained from the Barcode of Life Data system. We found that type I error rates, incorrect taxonomic assignments with a high bootstrap support, were already relatively low but could be lowered further by ensuring that all query taxa are actually present in the reference database. Choosing bootstrap support cut-offs according to query length and summarizing taxonomic assignments to more inclusive ranks can also help to reduce error while retaining the maximum number of assignments. Additionally, we highlight gaps in the taxonomic and geographic representation of insects in public sequence databases that will require further work by taxonomists to improve the quality of assignments generated using any method.

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.

Fungal high-throughput taxonomic identification tool for use with next-generation sequencing (FHiTINGS)

Improvements in DNA sequencing technology provide unprecedented opportunities to explore fungal diversity, but also present challenges in data analysis due to the large number of sequences generated. Here, we describe an open source software program "FHiTINGS" that utilizes the output of a BLASTn (blastall) search to rapidly identify, classify, and parse internal transcribed spacer (ITS) DNA sequences produced in fungal ecology studies that utilize next-generation DNA sequencing. This tool was designed for use with 454 pyrosequencing and is also appropriate for use with any sequencing platform that allows for BLAST searches against the indicated ITS database. For each sequence, FHiTINGS uses the lowest common ancestor method (LCA) to produce a single identification from BLAST output results, and then assigns taxonomic ranks from species through kingdom when possible for each sequence based on the Index Fungorum database. The program then sums and sorts this data into tabular form to enable rapid analysis of the sample, including α-diversity measures or richness. In silico testing demonstrates the time required to analyze and classify 1000 sequences is reduced from over 2 h by manual sorting to <1 min of computational time when using FHiTINGS, and that the classification output from the software is consistent with that derived from manual sorting of the data.

Filling gaps in biodiversity knowledge for macrofungi: contributions and assessment of an herbarium collection DNA barcode sequencing project

Despite recent advances spearheaded by molecular approaches and novel technologies, species description and DNA sequence information are significantly lagging for fungi compared to many other groups of organisms. Large scale sequencing of vouchered herbarium material can aid in closing this gap. Here, we describe an effort to obtain broad ITS sequence coverage of the approximately 6000 macrofungal-species-rich herbarium of the Museum of Natural History in Venice, Italy. Our goals were to investigate issues related to large sequencing projects, develop heuristic methods for assessing the overall performance of such a project, and evaluate the prospects of such efforts to reduce the current gap in fungal biodiversity knowledge. The effort generated 1107 sequences submitted to GenBank, including 416 previously unrepresented taxa and 398 sequences exhibiting a best BLAST match to an unidentified environmental sequence. Specimen age and taxon affected sequencing success, and subsequent work on failed specimens showed that an ITS1 mini-barcode greatly increased sequencing success without greatly reducing the discriminating power of the barcode. Similarity comparisons and nonmetric multidimensional scaling ordinations based on pairwise distance matrices proved to be useful heuristic tools for validating the overall accuracy of specimen identifications, flagging potential misidentifications, and identifying taxa in need of additional species-level revision. Comparison of within- and among-species nucleotide variation showed a strong increase in species discriminating power at 1-2% dissimilarity, and identified potential barcoding issues (same sequence for different species and vice-versa). All sequences are linked to a vouchered specimen, and results from this study have already prompted revisions of species-sequence assignments in several taxa.

Multiple origins of sequestrate basidiomes within Entoloma inferred from molecular phylogenetic analyses

The genus Entoloma comprises diverse trophic modes and basidiome morphologies. Although Entoloma includes some sequestrate species, their origins are not clearly understood in relation to phylogenetic position and trophic status. In this study, we collected 34 sequestrate Entoloma specimens in Japan over a 9-y period. Their identities and phylogenetic positions were determined by molecular analyses using three nuclear loci [internal transcribed spacer (ITS) and large subunit (LSU) regions of rDNA and RNA polymerase II second LSU (rpb2) gene]. Based on species delimitation of 97 % sequence matches in the ITS region, which is a suitable region for species-level identification of higher fungi, we identified four sequestrate Entoloma species. Molecular phylogenetic analyses that included all related sequences in the International Nucleotide Sequence Database revealed that the four sequestrate Entoloma species belonged to two major phylogroups. One of the phylogroups was Inocephalus-Cyanula, which is composed only of saprotrophic species. Three of the Japanese sequestrate species, as well as three previously known sequestrate species from other regions, fell into at least two independent clades in this phylogroup, indicating multiple origins of sequestrate forms within this saprotrophic lineage. Another phylogroup, Rhodopolioid, was also shown to include a sequestrate species for the first time. Because the Rhodopolioid phylogroup is composed exclusively of mycorrhizal species (ectomycorrhizal and tuberculate mycorrhizal species), the sequestrate form may also have evolved from a mycorrhizal ancestor. Our results suggest that sequestrate basidiomes have evolved multiple times, irrespective of their trophic status in Entoloma. Finally, based on molecular and morphological characteristics, here we describe two new sequestrate Entoloma species, i.e., Entoloma prismaticum Sasaki, Kinoshita et Nara, sp. nov. and Entoloma hypogaeum Sasaki, Kinoshita et Nara, sp. nov.

How well do ITS rDNA sequences differentiate species of true morels (Morchella)?

Arguably more mycophiles hunt true morels (Morchella) during their brief fruiting season each spring in the northern hemisphere than any other wild edible fungus. Concerns about overharvesting by individual collectors and commercial enterprises make it essential that science-based management practices and conservation policies are developed to ensure the sustainability of commercial harvests and to protect and preserve morel species diversity. Therefore, the primary objectives of the present study were to: (i) investigate the utility of the ITS rDNA locus for identifying Morchella species, using phylogenetic species previously inferred from multilocus DNA sequence data as a reference; and (ii) clarify insufficiently identified sequences and determine whether the named sequences in GenBank were identified correctly. To this end, we generated 553 Morchella ITS rDNA sequences and downloaded 312 additional ones generated by other researchers from GenBank using emerencia and analyzed them phylogenetically. Three major findings emerged: (i) ITS rDNA sequences were useful in identifying 48/62 (77.4%) of the known phylospecies; however, they failed to identify 12 of the 22 species within the species-rich Elata Subclade and two closely related species in the Esculenta Clade; (ii) at least 66% of the named Morchella sequences in GenBank are misidentified; and (iii) ITS rDNA sequences of up to six putatively novel Morchella species were represented in GenBank. Recognizing the need for a dedicated Web-accessible reference database to facilitate the rapid identification of known and novel species, we constructed Morchella MLST (http://www.cbs.knaw.nl/morchella/), which can be queried with ITS rDNA sequences and those of the four other genes used in our prior multilocus molecular systematic studies of this charismatic genus.

Reference databases for taxonomic assignment in metagenomics

Metagenomics is providing an unprecedented access to the environmental microbial diversity. The amplicon-based metagenomics approach involves the PCR-targeted sequencing of a genetic locus fitting different features. Namely, it must be ubiquitous in the taxonomic range of interest, variable enough to discriminate between different species but flanked by highly conserved sequences, and of suitable size to be sequenced through next-generation platforms. The internal transcribed spacers 1 and 2 (ITS1 and ITS2) of the ribosomal DNA operon and one or more hyper-variable regions of 16S ribosomal RNA gene are typically used to identify fungal and bacterial species, respectively. In this context, reliable reference databases and taxonomies are crucial to assign amplicon sequence reads to the correct phylogenetic ranks. Several resources provide consistent phylogenetic classification of publicly available 16S ribosomal DNA sequences, whereas the state of ribosomal internal transcribed spacers reference databases is notably less advanced. In this review, we aim to give an overview of existing reference resources for both types of markers, highlighting strengths and possible shortcomings of their use for metagenomics purposes. Moreover, we present a new database, ITSoneDB, of well annotated and phylogenetically classified ITS1 sequences to be used as a reference collection in metagenomic studies of environmental fungal communities. ITSoneDB is available for download and browsing at http://itsonedb.ba.itb.cnr.it/.

Asynchronous origins of ectomycorrhizal clades of Agaricales

The ectomycorrhizal (ECM) symbiosis is the most widespread biotrophic nutritional mode in mushroom-forming fungi. ECM fungi include, though are not limited to, about 5000 described species of Agaricales from numerous, independently evolved lineages. Two central hypotheses suggest different explanations for the origin of ECM fungal diversity: (i) dual origins, initially with the Pinaceae in the Jurassic and later with angiosperms during the Late Cretaceous, and (ii) a simultaneous and convergent radiation of ECM lineages in response to cooling climate during the Palaeogene and advancing temperate ECM plant communities. Neither of these hypotheses is supported here. While we demonstrate support for asynchronous origins of ECM Agaricales, the timing of such events appears to have occurred more recently than suggested by the first hypothesis, first during the Cretaceous and later during the Palaeogene. We are also unable to reject models of rate constancy, which suggests that the diversity of ECM Agaricales is not a consequence of convergent rapid radiations following evolutionary transitions from saprotrophic to ECM habits. ECM lineages of Agaricales differ not only in age, but also in rates of diversification and rate of substitution at nuclear ribosomal RNA loci. These results question the biological uniformity of the ECM guild.

DNA sequence analyses reveal abundant diversity, endemism and evidence for Asian origin of the porcini mushrooms

The wild gourmet mushroom Boletus edulis and its close allies are of significant ecological and economic importance. They are found throughout the Northern Hemisphere, but despite their ubiquity there are still many unresolved issues with regard to the taxonomy, systematics and biogeography of this group of mushrooms. Most phylogenetic studies of Boletus so far have characterized samples from North America and Europe and little information is available on samples from other areas, including the ecologically and geographically diverse regions of China. Here we analyzed DNA sequence variation in three gene markers from samples of these mushrooms from across China and compared our findings with those from other representative regions. Our results revealed fifteen novel phylogenetic species (about one-third of the known species) and a newly identified lineage represented by Boletus sp. HKAS71346 from tropical Asia. The phylogenetic analyses support eastern Asia as the center of diversity for the porcini sensu stricto clade. Within this clade, B. edulis is the only known holarctic species. The majority of the other phylogenetic species are geographically restricted in their distributions. Furthermore, molecular dating and geological evidence suggest that this group of mushrooms originated during the Eocene in eastern Asia, followed by dispersal to and subsequent speciation in other parts of Asia, Europe, and the Americas from the middle Miocene through the early Pliocene. In contrast to the ancient dispersal of porcini in the strict sense in the Northern Hemisphere, the occurrence of B. reticulatus and B. edulis sensu lato in the Southern Hemisphere was probably due to recent human-mediated introductions.

Are similarity- or phylogeny-based methods more appropriate for classifying internal transcribed spacer (ITS) metagenomic amplicons?

• The internal transcribed spacer (ITS) of the nuclear ribosomal DNA region is a widely used species marker for plants and fungi. Recent metagenomic studies using next-generation sequencing, however, generate only partial ITS sequences. Here we compare the performance of partial and full-length ITS sequences with several classification methods. • We compiled a full-length ITS data set and created short fragments to simulate the read lengths commonly recovered from current next-generation sequencing platforms. We compared recovery, erroneous recovery, and coverage for the following methods: best BLAST hit classification, MEGAN classification, and automated phylogenetic assignment using the Statistical Assignment Program (SAP). • We found that summarizing results with more inclusive taxonomic ranks increased recovery and reduced erroneous recovery. The similarity-based methods BLAST and MEGAN performed consistently across most fragment lengths. Using a phylogeny-based method, SAP runs with queries 400 bp or longer worked best. Overall, BLAST had the highest recovery rates and MEGAN had the lowest erroneous recovery rates. • A high-throughput ITS classification method should be selected, taking into consideration read length, an acceptable tradeoff between maximizing the total number of classifications and minimizing the number of erroneous classifications, and the computational speed of the assignment method.

Tidying up international nucleotide sequence databases: ecological, geographical and sequence quality annotation of its sequences of mycorrhizal fungi

Sequence analysis of the ribosomal RNA operon, particularly the internal transcribed spacer (ITS) region, provides a powerful tool for identification of mycorrhizal fungi. The sequence data deposited in the International Nucleotide Sequence Databases (INSD) are, however, unfiltered for quality and are often poorly annotated with metadata. To detect chimeric and low-quality sequences and assign the ectomycorrhizal fungi to phylogenetic lineages, fungal ITS sequences were downloaded from INSD, aligned within family-level groups, and examined through phylogenetic analyses and BLAST searches. By combining the fungal sequence database UNITE and the annotation and search tool PlutoF, we also added metadata from the literature to these accessions. Altogether 35,632 sequences belonged to mycorrhizal fungi or originated from ericoid and orchid mycorrhizal roots. Of these sequences, 677 were considered chimeric and 2,174 of low read quality. Information detailing country of collection, geographical coordinates, interacting taxon and isolation source were supplemented to cover 78.0%, 33.0%, 41.7% and 96.4% of the sequences, respectively. These annotated sequences are publicly available via UNITE (http://unite.ut.ee/) for downstream biogeographic, ecological and taxonomic analyses. In European Nucleotide Archive (ENA; http://www.ebi.ac.uk/ena/), the annotated sequences have a special link-out to UNITE. We intend to expand the data annotation to additional genes and all taxonomic groups and functional guilds of fungi.

Tasting soil fungal diversity with earth tongues: phylogenetic test of SATé alignments for environmental ITS data

An abundance of novel fungal lineages have been indicated by DNA sequencing of the nuclear ribosomal ITS region from environmental samples such as soil and wood. Although phylogenetic analysis of these novel lineages is a key component of unveiling the structure and diversity of complex communities, such analyses are rare for environmental ITS data due to the difficulties of aligning this locus across significantly divergent taxa. One potential approach to this issue is simultaneous alignment and tree estimation. We targeted divergent ITS sequences of the earth tongue fungi (Geoglossomycetes), a basal class in the Ascomycota, to assess the performance of SATé, recent software that combines progressive alignment and tree building. We found that SATé performed well in generating high-quality alignments and in accurately estimating the phylogeny of earth tongue fungi. Drawing from a data set of 300 sequences of earth tongues and progressively more distant fungal lineages, 30 insufficiently identified ITS sequences from the public sequence databases were assigned to the Geoglossomycetes. The association between earth tongues and plants has been hypothesized for a long time, but hard evidence is yet to be collected. The ITS phylogeny showed that four ectomycorrhizal isolates shared a clade with Geoglossum but not with Trichoglossum earth tongues, pointing to the significant potential inherent to ecological data mining of environmental samples. Environmental sampling holds the key to many focal questions in mycology, and simultaneous alignment and tree estimation, as performed by SATé, can be a highly efficient companion in that pursuit.

Dealing with incomplete taxon sampling and diversification of a large clade of mushroom-forming fungi

The absence of an adequate fossil record can hinder understanding the process of diversification that underlies the evolutionary history of a given group. In such cases, investigators have used ultrametric trees derived from molecular data from extant taxa to gain insights into processes of speciation and extinction over time. Inadequate taxon sampling, however, impairs such inferences. In this study, we use simulations to investigate the effect of incomplete taxon sampling on the accumulation of lineages through time for a clade of mushroom-forming fungi, the Hebelomateae. To achieve complete taxon sampling, we use a new Bayesian approach that incorporates substitute lineages to estimate diversification rates. Unlike many studies of animals and plants, we find no evidence of a slowdown in speciation. This indicates the Hebelomateae has not undergone an adaptive radiation. Rather, these fungi have evolved under a relatively constant rate of diversification since their most recent common ancestor, which we date back to the Eocene. The estimated net diversification rate (0.08-0.19 spp./lineage/Ma) is comparable with that of many plants and animals. We suggest that continuous diversification in the Hebelomateae has been facilitated by climatic and vegetation changes throughout the Cenozoic. We also caution against modeling multiple genes as a single partition when performing phylogenetic dating analyses.

Differentiating between effects of invasion and diversity: impacts of aboveground plant communities on belowground fungal communities

Exotic plant species can affect soil microbial communities with the potential for community and ecosystem feedbacks. Yet, separating the effects of exotics from confounded changes in plant community diversity still remains a challenge. We focused on how plant diversity and native or exotic life history affected root fungi because of their significant roles in community and ecosystem processes. Specifically, we examined how fungi colonizing plant roots were affected by plant richness (one, two or four species) replicated across a range of plant community mixtures (natives, exotics, native-exotic mixtures). Fungal biomass inside roots was affected independently by plant richness and mixture, while root fungal community composition was affected only by plant richness. Extraradical networks also increased in size with plant richness. By contrast, plant biomass was a function of plant mixture, with natives consistently smaller than exotics and native-exotic mixtures intermediate. Plant invasions may have an impact on the belowground community primarily through their effects on diversity, at least in the short-term. Disentangling the effects of diversity and invasion on belowground microbial communities can help us to understand both the controllers of belowground resilience and mechanisms of successful colonization and spread of exotic plants.

Molecular characterization of Morchella species from the Western Himalayan region of India

The molecular diversity of thirty-two different Morchella cultures/fruiting bodies, collected from the Western Himalayan region was studied in this investigation. Considerable taxonomic confusion exists regarding many species of Morchella. Although classical taxonomy is helpful in identification for many ascomycetes, morels exhibit considerable morphological diversity and there is disagreement in the identification of morel species. Phylogenetic analyses based on DNA sequences could help in sorting out morel taxonomy which is essential to better define the morel diversity. In this study, sequence analysis revealed that in the Western Himalayan region of India, both yellow (M. crassipes, M. spongiola) and black morels (M. elata, M. angusticeps, and M. gigas) were prominent along with two Verpa species. Phylogenetic analysis by maximum parsimony, maximum likelihood and Bayesian inference revealed two different clades and a clear distinction between yellow and black morels.

PlutoF—a Web Based Workbench for Ecological and Taxonomic Research, with an Online Implementation for Fungal ITS Sequences

DNA sequences accumulating in the International Nucleotide Sequence Databases (INSD) form a rich source of information for taxonomic and ecological meta-analyses. However, these databases include many erroneous entries, and the data itself is poorly annotated with metadata, making it difficult to target and extract entries of interest with any degree of precision. Here we describe the web-based workbench PlutoF, which is designed to bridge the gap between the needs of contemporary research in biology and the existing software resources and databases. Built on a relational database, PlutoF allows remote-access rapid submission, retrieval, and analysis of study, specimen, and sequence data in INSD as well as for private datasets though web-based thin clients. In contrast to INSD, PlutoF supports internationally standardized terminology to allow very specific annotation and linking of interacting specimens and species. The sequence analysis module is optimized for identification and analysis of environmental ITS sequences of fungi, but it can be modified to operate on any genetic marker and group of organisms. The workbench is available at http://plutof.ut.ee.

Community analysis reveals close affinities between endophytic and endolichenic fungi in mosses and lichens

Endolichenic fungi live in close association with algal photobionts inside asymptomatic lichen thalli and resemble fungal endophytes of plants in terms of taxonomy, diversity, transmission mode, and evolutionary history. This similarity has led to uncertainty regarding the distinctiveness of endolichenic fungi compared with endophytes. Here, we evaluate whether these fungi represent distinct ecological guilds or a single guild of flexible symbiotrophs capable of colonizing plants or lichens indiscriminately. Culturable fungi were sampled exhaustively from replicate sets of phylogenetically diverse plants and lichens in three microsites in a montane forest in southeastern Arizona (USA). Intensive sampling combined with a small spatial scale permitted us to decouple spatial heterogeneity from host association and to sample communities from living leaves, dead leaves, and lichen thalli to statistical completion. Characterization using data from the nuclear ribosomal internal transcribed spacer and partial large subunit (ITS-LSU rDNA) provided a first estimation of host and substrate use for 960 isolates representing five classes and approximately 16 orders, 32 families, and 65 genera of Pezizomycotina. We found that fungal communities differ at a broad taxonomic level as a function of the phylogenetic placement of their plant or lichen hosts. Endolichenic fungal assemblages differed as a function of lichen taxonomy, rather than substrate, growth form, or photobiont. In plants, fungal communities were structured more by plant lineage than by the living vs. senescent status of the leaf. We found no evidence that endolichenic fungi are saprotrophic fungi that have been "entrapped" by lichen thalli. Instead, our study reveals the distinctiveness of endolichenic communities relative to those in living and dead plant tissues, with one notable exception: we identify, for the first time, an ecologically flexible group of symbionts that occurs both as endolichenic fungi and as endophytes of mosses.

Rangewide analysis of fungal associations in the fully mycoheterotrophic Corallorhiza striata complex (Orchidaceae) reveals extreme specificity on ectomycorrhizal Tomentella (Thelephoraceae) across North America

Fully mycoheterotrophic plants offer a fascinating system for studying phylogenetic associations and dynamics of symbiotic specificity between hosts and parasites. These plants frequently parasitize mutualistic mycorrhizal symbioses between fungi and trees. Corallorhiza striata is a fully mycoheterotrophic, North American orchid distributed from Mexico to Canada, but the full extent of its fungal associations and specificity is unknown. Plastid DNA (orchids) and ITS (fungi) were sequenced for 107 individuals from 42 populations across North America to identify C. striata mycobionts and test hypotheses on fungal host specificity. Four largely allopatric orchid plastid clades were recovered, and all fungal sequences were most similar to ectomycorrhizal Tomentella (Thelephoraceae), nearly all to T. fuscocinerea. Orchid-fungal gene trees were incongruent but nonindependent; orchid clades associated with divergent sets of fungi, with a clade of Californian orchids subspecialized toward a narrow Tomentella fuscocinerea clade. Both geography and orchid clades were important determinants of fungal association, following a geographic mosaic model of specificity on Tomentella fungi. These findings corroborate patterns described in other fully mycoheterotrophic orchids and monotropes, represent one of the most extensive plant-fungal genetic investigations of fully mycoheterotrophic plants, and have conservation implications for the >400 plant species engaging in this trophic strategy worldwide.

454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity

* Soil fungi play a major role in ecological and biogeochemical processes in forests. Little is known, however, about the structure and richness of different fungal communities and the distribution of functional ecological groups (pathogens, saprobes and symbionts). * Here, we assessed the fungal diversity in six different forest soils using tag-encoded 454 pyrosequencing of the nuclear ribosomal internal transcribed spacer-1 (ITS-1). No less than 166 350 ITS reads were obtained from all samples. In each forest soil sample (4 g), approximately 30 000 reads were recovered, corresponding to around 1000 molecular operational taxonomic units. * Most operational taxonomic units (81%) belonged to the Dikarya subkingdom (Ascomycota and Basidiomycota). Richness, abundance and taxonomic analyses identified the Agaricomycetes as the dominant fungal class. The ITS-1 sequences (73%) analysed corresponded to only 26 taxa. The most abundant operational taxonomic units showed the highest sequence similarity to Ceratobasidium sp., Cryptococcus podzolicus, Lactarius sp. and Scleroderma sp. * This study validates the effectiveness of high-throughput 454 sequencing technology for the survey of soil fungal diversity. The large proportion of unidentified sequences, however, calls for curated sequence databases. The use of pyrosequencing on soil samples will accelerate the study of the spatiotemporal dynamics of fungal communities in forest ecosystems.

The ITS region as a target for characterization of fungal communities using emerging sequencing technologies

The advent of new high-throughput DNA-sequencing technologies promises to redefine the way in which fungi and fungal communities--as well as other groups of organisms--are studied in their natural environment. With read lengths of some few hundred base pairs, massively parallel sequencing (pyrosequencing) stands out among the new technologies as the most apt for large-scale species identification in environmental samples. Although parallel pyrosequencing can generate hundreds of thousands of sequences at an exceptional speed, the limited length of the reads may pose a problem to the species identification process. This study explores whether the discrepancy in read length between parallel pyrosequencing and traditional (Sanger) sequencing will have an impact on the perceived taxonomic affiliation of the underlying species. Based on all 39,200 publicly available fungal environmental DNA sequences representing the nuclear ribosomal internal transcribed spacer (ITS) region, the results show that the two approaches give rise to quite different views of the diversity of the underlying samples. Standardization of which subregion from the ITS region should be sequenced, as well as a recognition that the composition of fungal communities as depicted through different sequencing methods need not be directly comparable, appear crucial to the integration of the new sequencing technologies with current mycological praxis.