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.     

Towards a unified paradigm for sequence‐based identification of fungi

The nuclear ribosomal internal transcribed spacer (ITS) region is the formal fungal barcode and in most cases the marker of choice for the exploration of fungal diversity in environmental samples. Two problems are particularly acute in the pursuit of satisfactory taxonomic assignment of newly generated ITS sequences: (i) the lack of an inclusive, reliable public reference data set and (ii) the lack of means to refer to fungal species, for which no Latin name is available in a standardized stable way. Here, we report on progress in these regards through further development of the UNITE database ( http://unite.ut.ee ) for molecular identification of fungi. All fungal species represented by at least two ITS sequences in the international nucleotide sequence databases are now given a unique, stable name of the accession number type (e.g. Hymenoscyphus pseudoalbidus|GU586904|SH133781.05FU), and their taxonomic and ecological annotations were corrected as far as possible through a distributed, third‐party annotation effort. We introduce the term ‘species hypothesis’ (SH) for the taxa discovered in clustering on different similarity thresholds (97–99%). An automatically or manually designated sequence is chosen to represent each such SH. These reference sequences are released ( http://unite.ut.ee/repository.php ) for use by the scientific community in, for example, local sequence similarity searches and in the QIIME pipeline. The system and the data will be updated automatically as the number of public fungal ITS sequences grows. We invite everybody in the position to improve the annotation or metadata associated with their particular fungal lineages of expertise to do so through the new Web‐based sequence management system in UNITE.     

Description and identification of Alnus acuminata ectomycorrhizae from Argentinean alder stands

The objective of this study was to describe the morphological and anatomical features of five unidentified ectomycorrhizal types of Alnus acuminata and to complement their identification based on ITS-rDNA sequence analysis. The combined approach of morphotyping and sequence analysis based on ITS sequence comparison with sequences contained in GenBank and the UNITE database let us assign three of the five field-collected ectomycorrhiza morphotypes to the tomentella-thelephora lineage that closely matched European and North American species. The sequencing results within Tomentella point toward alder specific clades within T. sublilacina, T. ellisii and T. stuposa sensu lato. The two other EcM morphotypes matched Lactarius omphaliiformis and a Russula sp. Better focused, concomitant fruit body surveys are needed for accurate identification of South American ectomycorrhizal fungi because of the evidence of cryptic speciation in both agaricoid and resupinate mycobionts.     

ITS1 versus ITS2 as DNA metabarcodes for fungi

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

基于形态特征和ITS序列对7个鹅膏菌属菌株的分类鉴定

以采自浙江省丽水地区的7个鹅膏菌属菌株作为研究材料,在基于形态特征进行初步鉴定的基础上,对7种鹅膏菌的rDNAITS区段进行克隆测序和序列特征比较分析。进一步对ITS序列进行核酸序列数据库GenBank同源性检索比对,将从GenBank检索获得的9个最相似物种的ITS序列连同7种鹅膏菌的ITS序列一起作系统发育分析。结果表明:基于ITS序列对f6、f9和f493个菌株的分子鉴定支持了基于形态特征的鉴定结果,对f5的分子鉴定不支持形态鉴定的结果,f8为鹅膏菌属内某种,f66为鹅膏菌属内某种,并与Amanitafulva,A.atrofusca,A.orientifulva3种鹅膏菌的亲缘关系较近,f7与另外6种鹅膏菌的亲缘关系相差甚远。研究结果提示基于分子水平上的ITS序列分析不能单方面作为大型真菌分类鉴定的可靠依据,可以作为基于传统形态学分类鉴定的辅助参考依据。     

PHYMYCO-DB: A Curated Database for Analyses of Fungal Diversity and Evolution

Background

In environmental sequencing studies, fungi can be identified based on nucleic acid sequences, using either highly variable sequences as species barcodes or conserved sequences containing a high-quality phylogenetic signal. For the latter, identification relies on phylogenetic analyses and the adoption of the phylogenetic species concept.Such analysis requires that the reference sequences are well identified and deposited in public-access databases. However, many entries in the public sequence databases are problematic in terms of quality and reliability and these data require screening to ensure correct phylogenetic interpretation.

Methods and Principal Findings

To facilitate phylogenetic inferences and phylogenetic assignment, we introduce a fungal sequence database. The database PHYMYCO-DB comprises fungal sequences from GenBank that have been filtered to satisfy stringent sequence quality criteria. For the first release, two widely used molecular taxonomic markers were chosen: the nuclear SSU rRNA and EF1-α gene sequences. Following the automatic extraction and filtration, a manual curation is performed to remove problematic sequences while preserving relevant sequences useful for phylogenetic studies. As a result of curation, ∼20% of the automatically filtered sequences have been removed from the database. To demonstrate how PHYMYCO-DB can be employed, we test a set of environmental Chytridiomycota sequences obtained from deep sea samples.

Conclusion

PHYMYCO-DB offers the tools necessary to: (i) extract high quality fungal sequences for each of the 5 fungal phyla, at all taxonomic levels, (ii) extract already performed alignments, to act as ‘reference alignments’, (iii) launch alignments of personal sequences along with stored data. A total of 9120 SSU rRNA and 672 EF1-α high-quality fungal sequences are now available.The PHYMYCO-DB is accessible through the URL http://phymycodb.genouest.org/.     

Nucleotide sequence diversity of rDNA internal transcribed spacers extracted from conidia and cleistothecia of several powdery mildew fungi

An improved protocol, including DNA extraction with Chelex, two amplifications with a nested primer set, and DNA purification by electrophoresis, made it possible to analyze nuclear rDNA sequences of powdery mildew fungi using at most several hundred conidia or 20 cleistothecia. Nucleotide sequence diversity of the nuclear rDNA region containing the two internal transcribed spacers (ITS1 and ITS2) and 5.8S rRNA gene derived from conidia and cleistothecia was investigated for four kinds of powdery mildew fungi including two isolates of the same species. The results showed that the nucleotide sequences of the nuclear rDNA region were highly conserved between the teleomorph and the anamorph. Thus, the nucleotide sequence data obtained from either developmental stage can be used for phylogenetic studies of powdery mildew fungi. The nucleotide sequences of the 5.8S rRNA genes of the four species were highly conserved, but those of their ITS regions were variable. This suggests that the nuclear rDNA region is not suitable for phylogenetic studies of distantly related powdery mildew fungi, because too much sequence diversity exists, within the ITS, and too little phylogenetic information is contained within the 5.8S rRNA gene. However, the ITS region will be useful for phylogenetic comparison of closely related species or intraspecies. Contribution No. 132 from the Laboratory of Plant Pathology, Mie University.     

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.     

CORE: a phylogenetically-curated 16S rDNA database of the core oral microbiome

Comparing bacterial 16S rDNA sequences to GenBank and other large public databases via BLAST often provides results of little use for identification and taxonomic assignment of the organisms of interest. The human microbiome, and in particular the oral microbiome, includes many taxa, and accurate identification of sequence data is essential for studies of these communities. For this purpose, a phylogenetically curated 16S rDNA database of the core oral microbiome, CORE, was developed. The goal was to include a comprehensive and minimally redundant representation of the bacteria that regularly reside in the human oral cavity with computationally robust classification at the level of species and genus. Clades of cultivated and uncultivated taxa were formed based on sequence analyses using multiple criteria, including maximum-likelihood-based topology and bootstrap support, genetic distance, and previous naming. A number of classification inconsistencies for previously named species, especially at the level of genus, were resolved. The performance of the CORE database for identifying clinical sequences was compared to that of three publicly available databases, GenBank nr/nt, RDP and HOMD, using a set of sequencing reads that had not been used in creation of the database. CORE offered improved performance compared to other public databases for identification of human oral bacterial 16S sequences by a number of criteria. In addition, the CORE database and phylogenetic tree provide a framework for measures of community divergence, and the focused size of the database offers advantages of efficiency for BLAST searching of large datasets. The CORE database is available as a searchable interface and for download at http://microbiome.osu.edu.