Accurate, rapid taxonomic classification of fungal large-subunit rRNA genes

Taxonomic and phylogenetic fingerprinting based on sequence analysis of gene fragments from the large-subunit rRNA (LSU) gene or the internal transcribed spacer (ITS) region is becoming an integral part of fungal classification. The lack of an accurate and robust classification tool trained by a validated sequence database for taxonomic placement of fungal LSU genes is a severe limitation in taxonomic analysis of fungal isolates or large data sets obtained from environmental surveys. Using a hand-curated set of 8,506 fungal LSU gene fragments, we determined the performance characteristics of a naïve Bayesian classifier across multiple taxonomic levels and compared the classifier performance to that of a sequence similarity-based (BLASTN) approach. The naïve Bayesian classifier was computationally more rapid (>460-fold with our system) than the BLASTN approach, and it provided equal or superior classification accuracy. Classifier accuracies were compared using sequence fragments of 100 bp and 400 bp and two different PCR primer anchor points to mimic sequence read lengths commonly obtained using current high-throughput sequencing technologies. Accuracy was higher with 400-bp sequence reads than with 100-bp reads. It was also significantly affected by sequence location across the 1,400-bp test region. The highest accuracy was obtained across either the D1 or D2 variable region. The naïve Bayesian classifier provides an effective and rapid means to classify fungal LSU sequences from large environmental surveys. The training set and tool are publicly available through the Ribosomal Database Project (http://rdp.cme.msu.edu/classifier/classifier.jsp).     

Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys

Taxonomic classification of the thousands–millions of 16S rRNA gene sequences generated in microbiome studies is often achieved using a naïve Bayesian classifier (for example, the Ribosomal Database Project II (RDP) classifier), due to favorable trade-offs among automation, speed and accuracy. The resulting classification depends on the reference sequences and taxonomic hierarchy used to train the model; although the influence of primer sets and classification algorithms have been explored in detail, the influence of training set has not been characterized. We compared classification results obtained using three different publicly available databases as training sets, applied to five different bacterial 16S rRNA gene pyrosequencing data sets generated (from human body, mouse gut, python gut, soil and anaerobic digester samples). We observed numerous advantages to using the largest, most diverse training set available, that we constructed from the Greengenes (GG) bacterial/archaeal 16S rRNA gene sequence database and the latest GG taxonomy. Phylogenetic clusters of previously unclassified experimental sequences were identified with notable improvements (for example, 50% reduction in reads unclassified at the phylum level in mouse gut, soil and anaerobic digester samples), especially for phylotypes belonging to specific phyla (Tenericutes, Chloroflexi, Synergistetes and Candidate phyla TM6, TM7). Trimming the reference sequences to the primer region resulted in systematic improvements in classification depth, and greatest gains at higher confidence thresholds. Phylotypes unclassified at the genus level represented a greater proportion of the total community variation than classified operational taxonomic units in mouse gut and anaerobic digester samples, underscoring the need for greater diversity in existing reference databases.     

metaxa2: improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data

The ribosomal rRNA genes are widely used as genetic markers for taxonomic identification of microbes. Particularly the small subunit (SSU; 16S/18S) rRNA gene is frequently used for species‐ or genus‐level identification, but also the large subunit (LSU; 23S/28S) rRNA gene is employed in taxonomic assignment. The metaxa software tool is a popular utility for extracting partial rRNA sequences from large sequencing data sets and assigning them to an archaeal, bacterial, nuclear eukaryote, mitochondrial or chloroplast origin. This study describes a comprehensive update to metaxa – metaxa 2 – that extends the capabilities of the tool, introducing support for the LSU rRNA gene, a greatly improved classifier allowing classification down to genus or species level, as well as enhanced support for short‐read (100 bp) and paired‐end sequences, among other changes. The performance of metaxa 2 was compared to other commonly used taxonomic classifiers, showing that metaxa 2 often outperforms previous methods in terms of making correct predictions while maintaining a low misclassification rate. metaxa 2 is freely available from http://microbiology.se/software/metaxa2/ .     

Phylogenetic reference data for systematics and phylotaxonomy of arbuscular mycorrhizal fungi from phylum to species level

Although the molecular phylogeny, evolution and biodiversity of arbuscular mycorrhizal fungi (AMF) are becoming clearer, phylotaxonomically reliable sequence data are still limited. To fill this gap, a data set allowing resolution and environmental tracing across all taxonomic levels is provided. Two overlapping nuclear DNA regions, totalling c. 3 kb, were analysed: the small subunit (SSU) rRNA gene (up to 1800 bp) and a fragment spanning c. 250 bp of the SSU rDNA, the internal transcribed spacer (ITS) region (c. 475-520 bp) and c. 800 bp of the large subunit (LSU) rRNA gene. Both DNA regions together could be analysed for 35 described species, the SSU rDNA for c. 76 named and 18 as yet undefined species, and the ITS region or LSU rDNA, or a combination of both, for c. 91 named and 16 as yet undefined species. Present phylogenetic analyses, based on the three rDNA markers, provide reliable and robust resolution from phylum to species level. Altogether, 109 named species and 27 cultures representing as yet undefined species were analysed. This study provides a reference data set for molecular systematics and environmental community analyses of AMF, including analyses based on deep sequencing.     

Tubeufia asiana, the teleomorph of Aquaphila albicans in the tubeufiaceae, pleosporales, based on cultural and molecular data

The teleomorph of Aquaphila albicans was discovered on submerged wood collected in Thailand. Its black, soft-textured, setose ascomata, bitunicate asci and hyaline to pale brown, multiseptate ascospores indicated an affinity to Tubeufiaceae (Dothideomycetes). After morphological or molecular comparisons with related species in Tubeufia, Acanthostigma and Taphrophila, it is described and illustrated as a new species, T. asiana Sivichai & K.M. Tsui, sp. nov. Finding this Tubeufia teleomorph was surprising, given the falcate conidia of its A. albicans anamorph, which superficially resemble the conidia of Fusarium and not the coiled, helicosporous conidia of other species in Tubeufiaceae. We assessed the phylogenetic relationships of A. albicans-T. asiana with ribosomal sequences from SSU and ITS and partial LSU regions by parsimony and Bayesian analysis. An initial set of 40 taxa representing a wide range of ascomycete families and their SSU sequences from GenBank showed A. albicans-T. asiana to be nested within the Tubeufiaceae with 100% bootstrap support. Their placement was inferred with ITS and partial LSU ribosomal sequences. The nearly identical ITS sequences of two isolates of A. albicans and one isolate of Tubeufia asiana united these fungi as a monophyletic group with 100% bootstrap support and further nested them, with 88% bootstrap support, in a clade containing Helicoon gigantisporum and Helicoma chlamydosporum. This is the first molecular phylogenetic study to place a nonhelicosporous species within the Tubeufiaceae and to show that helical conidia were lost at least once within the family.     

Factors that affect large subunit ribosomal DNA amplicon sequencing studies of fungal communities: classification method, primer choice, and error

Nuclear large subunit ribosomal DNA is widely used in fungal phylogenetics and to an increasing extent also amplicon-based environmental sequencing. The relatively short reads produced by next-generation sequencing, however, makes primer choice and sequence error important variables for obtaining accurate taxonomic classifications. In this simulation study we tested the performance of three classification methods: 1) a similarity-based method (BLAST + Metagenomic Analyzer, MEGAN); 2) a composition-based method (Ribosomal Database Project na?ve bayesian classifier, NBC); and, 3) a phylogeny-based method (Statistical Assignment Package, SAP). We also tested the effects of sequence length, primer choice, and sequence error on classification accuracy and perceived community composition. Using a leave-one-out cross validation approach, results for classifications to the genus rank were as follows: BLAST + MEGAN had the lowest error rate and was particularly robust to sequence error; SAP accuracy was highest when long LSU query sequences were classified; and, NBC runs significantly faster than the other tested methods. All methods performed poorly with the shortest 50-100 bp sequences. Increasing simulated sequence error reduced classification accuracy. Community shifts were detected due to sequence error and primer selection even though there was no change in the underlying community composition. Short read datasets from individual primers, as well as pooled datasets, appear to only approximate the true community composition. We hope this work informs investigators of some of the factors that affect the quality and interpretation of their environmental gene surveys.     

Toward a better understanding of the infrageneric relationships in Cortinarius (Agaricales, Basidiomycota)

Research on the molecular systematics of Cortinarius, a species-rich mushroom genus with nearly global distribution, is just beginning. The present study explores infrageneric relationships using rDNA ITS and LSU sequence data. One large dataset of 132 rDNA ITS sequences and one combined da-taset with 54 rDNA ITS and LSU sequences were generated. Hebeloma was used as outgroup. Bayesian analyses and maximum-likelihood (ML) analyses were carried out. Bayesian phylogenetic inference performed equally well or better than ML, especially in large datasets. The phylogenetic analysis of the combined dataset with species representing all currently recognized subgenera recovered seven well-supported clades (Bayesian posterior probabilities BPP > 90%). These major clades are: /Myxacium s.l., /subg. Cortinarius, the /phlegmacioid clade (including the subclades /Phlegmacium and /Delibuti), the /calochroid clade (/Calochroi, /Ochroleuci and /Allutus), the /telamonioid clade (/Telamonia, /Orellani, /Anomali), /Dermocybe s.l. and /Myxotelamonia. Our results show that Cortinarius consists of many lineages, but the relationships among these clades could not be elucidated. On one hand, the low divergence in rDNA sequences can be held responsible for this; on the other hand, taxon sampling is problematic in Cortinarius phylogeny. Because of the incredibly high diversity (~2000 Cortinarius species), our sampling included <5% of the known species. By choosing type species of subgenera and sections, our sampling is strongly biased toward Northern Hemisphere taxa. More extensive taxon sampling, especially of species from the Southern Hemisphere, is essential to resolve the phylogeny of this important genus of ectomycorrhizal fungi.     

Galerina Earle: A polyphyletic genus in the consortium of dark-spored agarics

The basidiomycete genus Galerina Earle accommodates more than 300 small brown-spored agarics worldwide, predominantly described from the Northern hemisphere. The delimitation of species and infrageneric units hitherto has been based on morphological and, to some extent, ecological characters. In this study we have analyzed nuclear ribosomal LSU and ITS sequences to reveal infrageneric phylogeny and the phylogenetic placement of Galerina among the dark-spored agarics. Sequences from 36 northern hemisphere Galerina species and 19 other dark-spored taxa were analyzed, some of them obtained from EMBL/GenBank. Our results, received from Bayesian and distance methods, strongly suggest that Galerina is a polyphyletic genus. The LSU analysis shows that Galerina is composed of three or four separate monophyletic main groups. In addition, a few species cluster together with other dark-spored agarics. The same groups are recognized in the ITS tree and they correspond roughly to previously recognized subgenera or sections in Galerina. With high support our LSU analysis suggests that Gymnopilus is a monophyletic genus and that Gymnopilus and one of the Galerina lineages ("mycenopsis") are sister groups. The analyses further indicate that the Galerina lineages, as well as the genus Gymnopilus, could be referred to a strongly emendated family Strophariaceae, which corresponds largely to the family as circumscribed by Kühner (1980). Our results affirm that morphological characters often are highly homoplastic in the agarics. At the present stage formal taxonomic consequences or nomenclatural changes are not proposed.     

Multiple Independent Losses of Photosynthesis and Differing EvolutionaryRates in the Genus Cryptomonas (Cryptophyceae): Combined Phylogenetic Analyses of DNA Sequences of the Nuclear and the Nucleomorph Ribosomal Operons

In this study, evidence for at least three independent losses of photosynthesis in the freshwater cryptophyte genus Cryptomonas is presented. The phylogeny of the genus was inferred by molecular phylogenetic analyses of the nuclear internal transcribed spacer 2 (nuclear ITS2), partial nuclear large subunit ribosomal DNA (LSU rDNA), and nucleomorph small subunit ribosomal DNA (SSU rDNA, NM). Both concatenated and single data sets were used. In all data sets, the colorless Cryptomonas strains formed three different lineages, always supported by high bootstrap values (maximum parsimony, neighbor joining and maximum likelihood) and posterior probabilities (Bayesian analyses). The three leukoplast-bearing lineages displayed differing degrees of accelerated evolutionary rates in nuclear and nucleomorph rDNA. Also an increase in A+T-content in highly variable regions of the nucleomorph SSU rDNA was observed in one of the leukoplast-bearing lineages.This article contains three online-only supplementary tables.Reviewing Editor: Dr. Yves Van de Peer     

Complete sequence and structure of ribosomal RNA gene of Heterosporis anguillarum

The ribosomal RNA (rRNA) gene region of the microsporidium Heterosporis anguillarum has been examined. Complete DNA sequence data (4060 bp, GenBank Accession No. AF402839) of the rRNA gene of H. anguillarum are presented for the small subunit gene (SSU rRNA: 1359 bp), the internal transcribed spacer (ITS: 37 bp), and the large subunit gene (LSU rRNA: 2664 bp). The secondary structures of the H. anguillarum SSU and LSU rRNA genes are constructed and described. This is the first complete sequence of an rRNA gene published for a fish-infecting microsporidian species. In the phylogenetic analysis, the sequences, including partial SSU rRNA, ITS, and partial LSU rRNA sequences of the fish-infecting microsporidia, were aligned and analysed. The taxonomic position of H. anguillarum as suggested by Lom et al. (2000; Dis Aquat Org 43:225-231) is confirmed in this paper.     

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.     

Evolutionary relationships of the cup-fungus genus Peziza and Pezizaceae inferred from multiple nuclear genes: RPB2, beta-tubulin, and LSU rDNA

To provide a robust phylogeny of Pezizaceae, partial sequences from two nuclear protein-coding genes, RPB2 (encoding the second largest subunit of RNA polymerase II) and beta-tubulin, were obtained from 69 and 72 specimens, respectively, to analyze with nuclear ribosomal large subunit RNA gene sequences (LSU). The three-gene data set includes 32 species of Peziza, and 27 species from nine additional epigeous and six hypogeous (truffle) pezizaceous genera. Analyses of the combined LSU, RPB2, and beta-tubulin data set using parsimony, maximum likelihood, and Bayesian approaches identify 14 fine-scale lineages within Pezizaceae. Species of Peziza occur in eight of the lineages, spread among other genera of the family, confirming the non-monophyly of the genus. Although parsimony analyses of the three-gene data set produced a nearly completely resolved strict consensus tree, with increased confidence, relationships between the lineages are still resolved with mostly weak bootstrap support. Bayesian analyses of the three-gene data, however, show support for several more inclusive clades, mostly congruent with Bayesian analyses of RPB2. No strongly supported incongruence was found among phylogenies derived from the separate LSU, RPB2, and beta-tubulin data sets. The RPB2 region appeared to be the most informative single gene region based on resolution and clade support, and accounts for the greatest number of potentially parsimony informative characters within the combined data set, followed by the LSU and the beta-tubulin region. The results indicate that third codon positions in beta-tubulin are saturated, especially for sites that provide information about the deeper relationships. Nevertheless, almost all phylogenetic signal in beta-tubulin is due to third positions changes, with almost no signal in first and second codons, and contribute phylogenetic information at the "fine-scale" level within the Pezizaceae. The Pezizaceae is supported as monophyletic in analyses of the three-gene data set, but its sister-group relationships is not resolved with support. The results advocate the use of RPB2 as a marker for ascomycete phylogenetics at the inter-generic level, whereas the beta-tubulin gene appears less useful.     

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.     

Phylogenetic relationships among 28 spirotrichous ciliates documented by rDNA

The contiguous sequence of the SSU rDNA, ITS 1, 5.8S, ITS 2, and approximately 1370 bp at the 5(') end of the LSU rDNA was determined in 25 stichotrichs, one oligotrich, and two hypotrichs. Maximum parsimony, neighbor-joining, and quartet-puzzling analyses were used to construct individual phylogenetic trees for SSU rDNA, for LSU rDNA, and ITS 1+5.8S+ITS 2, as well as for all these components combined. All trees were similar, with the greatest resolution obtained with the combined components. Phylogenetic relationships were largely consistent with classical taxonomy, with notable disagreements. DNA sequences indicate that Oxytricha granulifera and Oxytricha longa are rather distantly related. The oligotrich, Halteria grandinella, is placed well within the order Stichotrichida. Uroleptus pisces and Uroleptus gallina probably belong to different genera. Holosticha polystylata (family Holostichidae) and Urostyla grandis (family Urostylidae) are rather closely related. These rDNA sequence analyses imply the need for some modifications of classical taxonomic schemes.     

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.     

The novel organization and complete sequence of the ribosomal RNA gene of Nosema bombycis

We present here for the first time the complete DNA sequence data (4301bp) of the ribosomal RNA (rRNA) gene of the microsporidian type species, Nosema bombycis. Sequences for the large subunit gene (LSUrRNA: 2497bp, GenBank Accession No. ), the internal transcribed spacer (ITS: 179bp, GenBank Accession No. ), the small subunit gene (SSUrRNA: 1232bp), intergenic spacer (IGS: 279bp), and 5S region (114bp) are also given, and the secondary structure of the large subunit is discussed. The organization of the N. bombycis rRNA gene is LSUrRNA-ITS-SSUrRNA-IGS-5S. This novel arrangement, in which the LSU is 5' of the SSU, is the reverse of the organizational sequence (i.e., SSU-ITS-LSU) found in all previously reported microsporidian rRNAs, including Nosema apis. This unique character in the type species may have taxonomic implications for the members of the genus Nosema.     

Molecular systematics of Helicoma, Helicomyces and Helicosporium and their teleomorphs inferred from rDNA sequences

Three genera of asexual, helical-spored fungi, Helicoma, Helicomyces and Helicosporium traditionally have been differentiated by the morphology of their conidia and conidiophores. In this paper we assessed their phylogenetic relationships from ribosomal sequences from ITS, 5.8S and partial LSU regions using maximum parsimony, maximum likelihood and Bayesian analysis. Forty-five isolates from the three genera were closely related and were within the teleomorphic genus Tubeufia sensu Barr (Tubeufiaceae, Ascomycota). Most of the species could be placed in one of the seven clades that each received 78% or greater bootstrap support. However none of the anamorphic genera were monophyletic and all but one of the clades contained species from more than one genus. The 15 isolates of Helicoma were scattered through the phylogeny and appeared in five of the clades. None of the four sections within the genus were monophyletic, although species from Helicoma sect. helicoma were concentrated in Clade A. The Helicosporium species also appeared in five clades. The four Helicomyces species were distributed among three clades. Most of the clades supported by sequence data lacked unifying morphological characters. Traditional characters such as the thickness of the conidial filament and whether conidiophores were conspicuous or reduced proved to be poor predictors of phylogenetic relationships. However some combinations of characters including conidium colour and the presence of lateral, tooth-like conidiogenous cells did appear to be predictive of genetic relationships.     

Molecular phylogeny of penaeid shrimps inferred from two mitochondrial markers

Penaeid shrimps are an important resource in crustacean fisheries, representing more than the half of the gross production of shrimp worldwide. In the present study, we used a sample of wide-ranging diversity (41 shrimp species) and two mitochondrial markers (758 bp) to clarify the evolutionary relationships among Penaeidae genera. Three different methodologies of tree reconstruction were employed in the study: maximum likelihood, neighbor joining and Bayesian analysis. Our results suggest that the old Penaeus genus is monophyletic and that the inclusion of the Solenocera genus within the Penaeidae family remains uncertain. With respect to Metapenaeopsis monophyly, species of this genus appeared clustered, but with a nonsignificant bootstrap value. These results elucidate some features of the unclear evolution of Penaeidae and may contribute to the taxonomic characterization of this family.     

Analysis of sequences of ITS1 internal transcribed spacer and 5.8S ribosome gene of Malus species

The nucleotide sequence of the ITS1-5.8S ribosomal DNA spacer fragment was determined for 41 samples of the Malus species. The total length of compared sequences ranged from 389 to 392 bp. The nucleotide sequence of the 5.8S gene within the genus was highly conserved. The level of polymorphism of ITS1 region comprised 14%. Both species- and group-specific substitutions were identified. The analysis of M. orientalis and M. turkmenorum sequences revealed their full identity, which indicates the need to perform more research with a larger number of samples of both species from other collections to clarify the taxonomic status of the M. turkmenorum species. The previous findings on the synonymy of species M. baccata, M. mandshurica, M. pallasiana, and M. sachalinensis were also confirmed.