Marisol Sánchez-García

Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data

Summary 1. The nuclear ribosomal internal transcribed spacer (ITS) region is the primary choice for molecular identification of fungi. Its two highly variable spacers (ITS1 and ITS2) are usually species specific, whereas the intercalary 5.8S gene is highly conserved. For sequence clustering and BLAST searches, it is often advantageous to rely on either one of the variable spacers but not the conserved 5.8S gene. To identify and extract ITS1 and ITS2 from large taxonomic and environmental data sets is, however, often difficult, and many ITS sequences are incorrectly delimited in the public sequence databases. 2. We introduce ITSx, a Perl-based software tool to extract ITS1, 5.8S and ITS2 – as well as full-length ITS sequences – from both Sanger and high-throughput sequencing data sets. ITSx uses hidden Markov models computed from large alignments of a total of 20 groups of eukaryotes, including fungi, metazoans and plants, and the sequence extraction is based on the predicted positions of the ribosomal genes in the sequences. 3. ITSx has a very high proportion of true-positive extractions and a low proportion of false-positive extractions. Additionally, process parallelization permits expedient analyses of very large data sets, such as a one million sequence amplicon pyrosequencing data set. ITSx is rich in features and written to be easily incorporated into automated sequence analysis pipelines. 4. ITSx paves the way for more sensitive BLAST searches and sequence clustering operations for the ITS region in eukaryotes. The software also permits elimination of non-ITS sequences from any data set. This is particularly useful for amplicon-based next-generation sequencing data sets, where insidious non-target sequences are often found among the target sequences. Such non-target sequences are difficult to find by other means and would contribute noise to diversity estimates if left in the data set.

Improving ITS sequence data for identification of plant pathogenic fungi

SummaryPlant pathogenic fungi are a large and diverse assemblage of eukaryotes with substantial impacts on natural ecosystems and human endeavours. These taxa often have complex and poorly understood life cycles, lack observable, discriminatory morphological characters, and may not be amenable to in vitro culturing. As a result, species identification is frequently difficult. Molecular (DNA sequence) data have emerged as crucial information for the taxonomic identification of plant pathogenic fungi, with the nuclear ribosomal internal transcribed spacer (ITS) region being the most popular marker. However, international nucleotide sequence databases are accumulating numerous sequences of compromised or low-resolution taxonomic annotations and substandard technical quality, making their use in the molecular identification of plant pathogenic fungi problematic. Here we report on a concerted effort to identify high-quality reference sequences for various plant pathogenic fungi and to re-annotate incorrectly or insufficiently annotated public ITS sequences from these fungal lineages. A third objective was to enrich the sequences with geographical and ecological metadata. The results – a total of 31,954 changes – are incorporated in and made available through the UNITE database for molecular identification of fungi (http://unite.ut.ee), including standalone FASTA files of sequence data for local BLAST searches, use in the next-generation sequencing analysis platforms QIIME and mothur, and related applications. The present initiative is just a beginning to cover the wide spectrum of plant pathogenic fungi, and we invite all researchers with pertinent expertise to join the annotation effort.

Into and out of the tropics : global diversification patterns in a hyperdiverse clade of ectomycorrhizal fungi

Ectomycorrhizal (ECM) fungi, symbiotic mutualists of many dominant tree and shrub species, exhibit a biogeographic pattern counter to the established latitudinal diversity gradient of most macroflora and fauna. However, an evolutionary basis for this pattern has not been explicitly tested in a diverse lineage. In this study, we reconstructed a mega‐phylogeny of a cosmopolitan and hyperdiverse genus of ECM fungi, Russula, sampling from annotated collections and utilizing publically available sequences deposited in GenBank. Metadata from molecular operational taxonomic unit cluster sets were examined to infer the distribution and plant association of the genus. This allowed us to test for differences in patterns of diversification between tropical and extratropical taxa, as well as how their associations with different plant lineages may be a driver of diversification. Results show that Russula is most species‐rich at temperate latitudes and ancestral state reconstruction shows that the genus initially diversified in temperate areas. Migration into and out of the tropics characterizes the early evolution of the genus, and these transitions have been frequent since this time. We propose the ‘generalized diversification rate’ hypothesis to explain the reversed latitudinal diversity gradient pattern in Russula as we detect a higher net diversification rate in extratropical lineages. Patterns of diversification with plant associates support host switching and host expansion as driving diversification, with a higher diversification rate in lineages associated with Pinaceae and frequent transitions to association with angiosperms.

How to know the fungi: combining field inventories and DNA‐barcoding to document fungal diversity

The fungi kingdom is among the most diverse eukaryotic lineages on Earth with estimates of several million extant species (O’Brien et al., 2005; Blackwell, 2011; Taylor et al., 2014). Fungi play critical roles in carbon andnutrient cycling of terrestrial and aquatic ecosystems, and they are important pathogens and mutualists (Read & Perez-Moreno, 2003; Taylor et al., 2012; Grossart et al., 2016). More than 80% of plant species form symbioses with fungi and these symbioses have been crucial to the colonization of terrestrial ecosystems (Field et al., 2015a; Selosse et al., 2015). Despite their impacts on primary ecosystem functions, assessments of fungal biodiversity estimate that only c. 10% of fungal species have been described (Bass & Richards, 2011; Hibbett et al., 2011). Traditionally, specimen-based taxonomic studies have been the only way to discover new species. Because most fungi have microscopic life-stages and convergent morphological features (Rivas-Plata & Lumbsch, 2011; Wynns, 2015), many fungal groups remain severely undersampled. DNA-barcoding and highthroughput sequencing methods have provided a new framework for studying fungal biodiversity (Fierer et al., 2012; Schoch et al., 2012; Myrold et al., 2014), and diversity estimates based on environmental sequences have increased exponentially. Although these ‘sequence-based classification and identification’ methods are a powerful means to rapidly detect hidden diversity, careful interpretation of these data is needed to make accurate inferences (K~oljalg et al., 2013; Lindahl et al., 2013; Nguyen et al., 2015; Hibbett et al., 2016). In particular, many environmental sequences cannot be associated with a known fungal species or lineage. This remains a major challenge to decipher fungal community composition and understand ecological roles of fungi in leaf litter, soil, or inside plants (Yahr et al., 2016). In some cases, these fungi are truly undescribed and their ecological roles are unknown but in other cases they represent described taxa for which no sequence is available (Nagy et al., 2011; Nilsson et al., 2016). DNA barcoding of herbarium specimens and culture collections is extremely valuable to link unidentified sequences to known taxa (e.g. Brock et al., 2009; Nagy et al., 2011; Osmundson et al., 2013; Garnica et al., 2016).DNA sequences have been generated from fungal type specimens > 200 years old (Larsson & Jacobsson, 2004), but in many cases obtaining sequences from historical material is challenging (Dentinger et al., 2010). Today’s threats to biodiversity from habitat loss and climate change are occurring at an unprecedented scale, and it is possible that many species may become extinct before they have been discovered (Costello et al., 2013; Monastersky, 2014). In the need to describe and protect as many species as possible we addressed the following questions: what are the best methods to rapidly document fungal biodiversity? Are traditional, specimen-based approaches still useful?

Deconstructing the Tricholomataceae (Agaricales) and introduction of the new genera Albomagister, Corneriella, Pogonoloma and Pseudotricholoma

The family Tricholomataceae, contained within the Tricholomatoid clade, has traditionally been one of the largest families of the Agaricales. However, in this sense it is highly polyphyletic and requires emendation. Here, we present a phylogeny of the Tricholomatoid clade based on nucleotide sequence data from two nuclear ribosomal RNA genes (large subunit and small subunit) and the second-largest subunit of RNA polymerase II (rpb2). Our aim is to delimit the Tricholomataceae and identify monophyletic groups within the Tricholomatoid clade. We also infer a separate phylogeny, based on the three genes above, in addition to sequences of the nuclear ribosomal internal transcribed spacers (ITS), in order to evaluate generic-level boundaries within the Tricholomataceae s.str. Based on this analysis we recover seven monophyletic genera within the Tricholomataceae s.str. that correspond to Leucopaxillus, Tricholoma, Pseudotricholoma stat. nov., Porpoloma s.str., Dennisiomyces, Corneriella gen. nov., and Albomagister gen. nov. Of the 98 genera that have been traditionally assigned to the Tricholomataceae sensu Singer, only four can be placed within it ( Tricholoma, Porpoloma, Dennisiomyces, Leucopaxillus). The genus Porpoloma is highly polyphyletic and divided into four genera: Porpoloma s.str., Corneriella gen. nov., Pseudotricholoma stat. nov., and Pogonoloma stat. nov. In all, four new genera are proposed. Taxonomic descriptions, and a key to genera of the Tricholomat - aceae as emended here are also presented.

Is the switch to an ectomycorrhizal state an evolutionary key innovation in mushroom-forming fungi? a case study in the tricholomatineae (agaricales)

Although fungi are one of the most diverse groups of organisms, little is known about the processes that shape their high taxonomic diversity. This study focuses on evolution of ectomycorrhizal (ECM) mushroom‐forming fungi, symbiotic associates of many trees and shrubs, in the suborder Tricholomatineae of the Agaricales. We used the BiSSE model and BAMM to test the hypothesis that the ECM habit represents an evolutionary key innovation that allowed the colonization of new niches followed by an increase in diversification rate. Ancestral state reconstruction (ASR) supports the ancestor of the Tricholomatineae as non‐ECM. We detected two diversification rate increases in the genus Tricholoma and the Rhodopolioid clade of the genus Entoloma. However, no increases in diversification were detected in the four other ECM clades of Tricholomatineae. We suggest that diversification of Tricholoma was not only due to the evolution of the ECM lifestyle, but also to the expansion and dominance of its main hosts and ability to associate with a variety of hosts. Diversification in the Rhodopolioid clade could be due to the unique combination of spore morphology and ECM habit. The spore morphology may represent an exaptation that aided spore dispersal and colonization. This is the first study to investigate rate shifts across a phylogeny that contains both non‐ECM and ECM lineages.