- Open Access
Approaching the taxonomic affiliation of unidentified sequences in public databases – an example from the mycorrhizal fungi
BMC Bioinformaticsvolume 6, Article number: 178 (2005)
During the last few years, DNA sequence analysis has become one of the primary means of taxonomic identification of species, particularly so for species that are minute or otherwise lack distinct, readily obtainable morphological characters. Although the number of sequences available for comparison in public databases such as GenBank increases exponentially, only a minuscule fraction of all organisms have been sequenced, leaving taxon sampling a momentous problem for sequence-based taxonomic identification. When querying GenBank with a set of unidentified sequences, a considerable proportion typically lack fully identified matches, forming an ever-mounting pile of sequences that the researcher will have to monitor manually in the hope that new, clarifying sequences have been submitted by other researchers. To alleviate these concerns, a project to automatically monitor select unidentified sequences in GenBank for taxonomic progress through repeated local BLAST searches was initiated. Mycorrhizal fungi – a field where species identification often is prohibitively complex – and the much used ITS locus were chosen as test bed.
A Perl script package called emerencia is presented. On a regular basis, it downloads select sequences from GenBank, separates the identified sequences from those insufficiently identified, and performs BLAST searches between these two datasets, storing all results in an SQL database. On the accompanying web-service http://emerencia.math.chalmers.se, users can monitor the taxonomic progress of insufficiently identified sequences over time, either through active searches or by signing up for e-mail notification upon disclosure of better matches. Other search categories, such as listing all insufficiently identified sequences (and their present best fully identified matches) publication-wise, are also available.
The ever-increasing use of DNA sequences for identification purposes largely falls back on the assumption that public sequence databases contain a thorough sampling of taxonomically well-annotated sequences. Taxonomy, held by some to be an old-fashioned trade, has accordingly never been more important. emerencia does not automate the taxonomic process, but it does allow researchers to focus their efforts elsewhere than countless manual BLAST runs and arduous sieving of BLAST hit lists. The emerencia system is available on an open source basis for local installation with any organism and gene group as targets.
Mycorrhiza is the term used to denote the root-associated symbiosis between fungus and plant where resources otherwise unattainable or very costly (chiefly carbohydrates and mineral nutrients, respectively) are exchanged. Whereas the identification of the plant partner of the symbiosis often is comparatively straightforward, the identity of the fungal component is typically much more elusive, largely so due to the dearth of information obtainable from the root samples of the plant (Figure 1). The traditional approaches to identification of mycorrhizae include studies based on light microscopy, isozyme assays, mating behaviour experiments, and somatic compatibility tests. However, all of these are associated with drawbacks such as low to moderate precision, high time consumption, or the requirement that the fungus be isolated and grown in culture, which is impossible for many mycorrhizal fungi [1–4]. Alternatively, it is sometimes possible to establish a hyphal connection between the fungal mycelium of the root tips and nearby fungal fruiting-bodies ; traditional fungal taxonomy rests to a large extent on the morphology of fruiting-bodies or other spore producing structures, and there is abundant literature information available for many groups of fungi. The factors triggering the formation of fruiting-bodies of mycorrhizal fungi are, however, poorly understood. This is reflected in the large number of root-associated fungi for which fruiting-bodies have never been found, suggesting that any attempt to characterize the below-ground mycoflora through collection and identification of above-ground fruiting-bodies is likely to give a skewed and incomplete picture [6–8].
The advent of PCR-based methodologies has in many ways revolutionized the field of mycorrhizal research by providing a means by which, at least in theory, every fungus could be identified to species level. It has also served to highlight further the patchiness in our understanding of mycorrhizal fungi; DNA-based studies of mycorrhizal fungi typically contain a considerable proportion of sequences that lack identified matches in public sequence repositories such as GenBank and UNITE , and that hence cannot be identified to species or even genus level [10–13]. The prerequisite of most journals that all sequences used in a publication be made public in GenBank naturally leads authors to submit such unidentified sequences under names like "Fungus: Environmental sample", "Unidentified mycorrhizal basidiomycete", "Uncultured root-tip fungus", and "Cortinarius sp.". Thus, even if an unidentified mycorrhizal sequence turns out to have one or more identical (or nearly so) matches in GenBank, the unidentified nature of the latter still precludes identification to species level.
As more and more sequences from well-identified fruiting-bodies are submitted to and accumulate in GenBank, one can expect some of the new sequences to hint and even resolve the identity of previously unmatched sequences. Any such relation is likely to impart important information to the studies that employed the unidentified sequences in question. Yet there is no generic mechanism on GenBank to alert the sequences' authors – or anyone else who might be interested – of these new matches. In addition, many sequence authors are inefficient in keeping the annotations of their submitted sequences up-to-date, even when new taxonomic progress and knowledge have been attained. As a result, important information goes by undetected, and obsolete annotations remain and are propagated through the databases and literature via subsequent BLAST runs . The present study introduces emerencia, a Perl script package to facilitate the keeping track of the identity of insufficiently identified GenBank sequences over time. A web-service demonstrating emerencia http://emerencia.math.chalmers.se for mycorrhizal fungi and the nuclear ITS region is presented; emerencia is also available on an open source-basis and can be downloaded from the web page for local installation with any gene and organism group as targets.
The web-service provided at http://emerencia.math.chalmers.se is set to monitor the identity of fungal ITS sequences whose taxonomic annotations are poorly resolved. The ITS region offers three sub-loci of very different conservation levels (ITS1- very variable; 5.8S- very conserved; ITS2- semi-conserved ) and is as such the prevalent region of choice when fungi are to be identified [16, 17, 3].
The main script of the emerencia package is written in Perl  and runs in a UNIX environment. On a regular basis, the script is evoked and connects to GenBank through BioPerl routines  to retrieve all fungal ITS sequences deposited since the last update using the query phrase
(("Fungi" [Organism] AND (200 [SLEN] : 3000 [SLEN])) AND (((ITS1 [titl] OR ITS2 [titl]) OR 5.8S [titl]) OR "internal transcribed spacer" [titl] OR "internal transcribed spacers" [titl]))
For each downloaded sequence, the GenBank format annotation is parsed to determine whether the sequence should be regarded as fully identified (i.e., identified to species level) or insufficiently identified (i.e., identified only to phylum Fungi (e.g., "Unidentified fungus"), identified to ordinal level (e.g., "Thelephorales sp."), or identified to generic level (e.g., "Amanita sp.")). The criteria for this decision can be found in full detail at http://emerencia.math.chalmers.se/criteria.html; for example, a sequence is regarded as insufficiently identified if its species annotation (GenBank format reference: SOURCE) contains words like "unidentified", "uncultured", "sp.", "mycorrhizal isolate", or "unknown" (this section of the script was repeatedly fine-tuned to minimize the number of false positives and negatives). All new sequences are appended to one of two tables of a local MySQL database  such that one table holds the identified sequences and the other the insufficiently identified ones. The structure of the database is provided at http://emerencia.math.chalmers.se/criteria.html. Local BLAST search indices  are computed to allow for similarity searches in both tables.
Two sets of local BLAST searches are then run using default settings:
1. all insufficiently identified sequences are BLAST:ed against all identified sequences
2. all insufficiently identified sequences are BLAST:ed against all insufficiently identified sequences.
For the first BLAST run (1), details of the best BLAST match are inspected. If, for any insufficiently identified sequence, the best BLAST match to the table of identified sequences is found to have changed (i.e., a change in the accession number of the best BLAST match in combination with an improved E-value or identical E-value and an improved score) – or if the unidentified sequence lacks any previous significant match – details of the new best BLAST match (accession number, date, and BLAST score and E-value) are appended to the entry of the insufficiently identified sequence. Similarly, for the second BLAST run (2), the best non-self match of each insufficiently identified sequence to the table of insufficiently identified sequences is noted and saved. To retain a BLAST history, the former best BLAST match of each insufficiently identified sequence is also kept. The end product of the main script of emerencia is, thus, two updated, mutually exclusive MySQL tables – one with identified fungal ITS sequences and one with insufficiently identified ITS sequences, and both with cross-linked indices on best and former best BLAST matches in each table.
The web-service enables visitors to interact with the database in a number of ways. Four major search categories are provided (Table 1):
SEARCH FOR INSUFFICIENTLY IDENTIFIED SEQUENCE BY ACCESSION NUMBER
CHECK SPECIFIC PUBLICATION FOR INSUFFICIENTLY IDENTIFIED SEQUENCES AND THEIR IDENTITY
SEARCH FOR INSUFFICIENTLY IDENTIFIED SEQUENCES MATCHING ACCESSION NUMBER OF IDENTIFIED TAXA
SEARCH FOR INSUFFICIENTLY IDENTIFIED SEQUENCES BY KEY WORD
The different search functions produce output pages that are extensively hyperlinked to facilitate further queries against the web-service itself as well as external information resources such as GenBank and the Tree of Life project . Apart from querying the database through a web browser and bookmarking or saving the results, users can subscribe to accession numbers of insufficiently identified sequences. This enables immediate notification by e-mail when the best BLAST matches of those sequences change.
The web-service is hosted on a MacOS X server running the Apache web server . The databases are queried using dedicated CGI scripts written in Perl; parts of the source code of the galaxieEST and mor packages [24, 17] were used for this purpose. The extensively annotated source code of the emerencia core is freely available at the web-service. Local installation and additional technical aspects are described in the documentation.
Results and Discussion
The last decade has seen a dramatic improvement of our understanding of mycorrhizal diversity, largely due to the advent of fast and comparatively cheap PCR-based methods. The discovery of new, previously unsequenced mycorrhizal fungi poses something of a taxonomic problem, particularly when fruiting-bodies and other distinguishing characteristics are absent. Many of these sequences are submitted ad tempus as "environmental samples". Unfortunately, the absence of generic mechanisms – and the apparent lack of motivation of the authors of the sequences – to refine the identity of these sequences as more information is amassed force other researchers to put in a great deal of effort (typically countless manual BLAST runs) in order to make sense of the sequences and the relation of their own sequences to those. With more than a handful of sequences to monitor, the task quickly becomes unreasonably time-consuming. The authors present a prototype software package to minimize the amount of work needd to stay updated on the identity of such insufficiently identified sequences. The web-service provided allows users to subscribe to accession numbers (sequences) with automatic email notification upon identity changes; alternatively, the same – and additional – information can be obtained through the search functions of the web-service. Furthermore, to install emerencia locally and modify it to run with other organism and gene groups should not pose any problem to anyone with a reasonable experience of UNIX-type environments. Such a local installation can stretch from a private, shell window-only tool to a public, user-oriented web-service such as the one presented here. For a local installation, parameters such as how often the script should be started, the BLAST settings, and what information to store locally, can be set as seen fit.
As of May 2005, emerencia has fetched about 29000 identified and 7500 insufficiently identified ITS sequences (Table 2). The identified sequences belong to some 8000 distinct species, which corresponds to approximately 0.5% of the estimated 1.5 million extant species of fungi . While the number of fungal sequences in GenBank is expected to increase drastically over the next few years, it will take a long time before all gaps are filled , leaving taxon sampling a tangible problem for emerencia as well as for other tools used for sequence identification. In addition, the poor state of many taxonomic annotations in GenBank and other databases [27, 28] complicates the above percentage estimates and poses a challenge to users of emerencia. As with other identification tools, it is crucial that the results obtained be viewed and used as guidance for further studies rather than accepted as true and final; emerencia is a tool to refine iteratively the identity of insufficiently identified sequences in public databases and to promote the flow of information pertaining to those sequences. It is not intended – and should never be used – as a shortcut to unequivocally correct species names and annotations.
As with BLAST searches in general, several factors impede the interpretation of the result. The aforementioned problem with taxonomic annotations in GenBank calls, in itself, for subsequent hands-on verification of the results. Furthermore, BLAST explores – and tries to expand – local regions of sequence similarity, and it takes manual inspection of the BLAST results to find out whether the entire, or only a portion of, the query sequence was successfully matched to anything in the database. A match to only a part of the target sequence (such as the very conserved 5.8S sub-locus of the ITS region) is, for identification purposes, tantamount to no match at all . It is also important to keep in mind that BLAST provides a measure of similarity, but similarity does not in turn provide a sound measure of relatedness [29, 30]. Finally, it is notoriously difficult to tell an identified sequence apart from an insufficiently identified one on an automated basis; indeed, the present set-up is likely to yield a small proportion of false positives as well as false negatives (presently less than 0.1%). Such problems would largely have been avoided had there been an accepted standard for annotation of unidentified – and identified – sequences.
emerencia bears some resemblance to tools like Sequence Alerting System , Swiss-Shop , and ReHab , but a number of features set emerencia apart from these. emerencia is primarily a taxonomic utility designed to add an integrative aspect to GenBank data; its automated separation of identified and insufficiently identified sequences paves the way for researchers seeking reliable identification of species rather than merely the best possible match scores. emerencia can be installed locally or accessed over the Internet; in the latter case, the user will need nothing but a web browser. The data structure of emerencia allows many types of interesting questions to be asked; for instance, insufficiently identified sequences and their present best identified matches can be listed publication-wise, or all insufficiently identified sequences that constitute the best BLAST match of some given identified sequence can be listed (essentially amounting to a BLAST run in reverse) (Table 1). emerencia is tailored for variable nucleotide sequences, whereas proteins represent the target for Swiss-Shop and ReHab. Finally, the e-mail subscription utility provides a convenient way for users to stay taxonomically updated on select insufficiently identified sequences with a minimum of effort.
Insufficiently identified sequences generally add little to the studies in which they are included, and it is important to estimate their identity as correctly as possible. The magnitude of this manual task increases with the number of sequences, but this process can fortunately be automated. However, in spite of computational advances, the taxonomic process itself lies beyond automation, alluding to the importance of both good species knowledge and the inherent need to always approach hypothesized identifications in a critical way.
Availability and requirements
Project name: emerencia
Project home page: http://emerencia.math.chalmers.se
Operating system(s): Primarily UNIX type platforms
Programming language: PERL, SQL
Other requirements: BLAST, Apache httpd, BioPerl, Clustal W (optional)
License: GNU GPL version 2
Any restrictions to use by non-academics: None other than those imposed by GNU GPL version 2
Basic Local Alignment Search Tool
- ITS :
Internal Transcribed Spacers
Structured Query Language
Gardes M, White TJ, Fortin JA, Bruns TD, Taylor JW: Identification of indigenous and introduced symbiotic fungi in ectomycorrhizae by amplification of nuclear and mitochondrial ribosomal DNA. Can J Botany 1991, 69: 180–190.
Marmaisse R, Debaud JC, Casselton LA: DNA probes for species and strain identification in the ectomycorrhizal fungus Hebeloma . MycolRes 1992, 96: 161–165.
Horton TR, Bruns TD: The molecular revolution in ectomycorrhizal ecology: peeking into the black-box. Mol Ecol 2001, 10: 1855–1871. 10.1046/j.0962-1083.2001.01333.x
Kõljalg U, Tammi H, Timonen S, Agerer R, Sen R: ITS rDNA sequence-based phylogenetic analysis of Tomentellopsis species from boreal and temperate forests, and the identification of pink-type ectomycorrhizas. Mycol Prog 2002, 1: 81–92.
Agerer R: Characterization of ectomycorrhiza. Method Microbiol 1991, 23: 25–73.
Gardes M, Bruns TD: Community structure of ectomycorrhizal fungi in a Pinus muricata forest: above-and below-ground views. Can J Botany 1996, 74: 1572–1583.
Dahlberg A: Community ecology of ectomycorrhizal fungi: an advancing interdisciplinary field. New Phytol 2001, 150: 555–562. 10.1046/j.1469-8137.2001.00142.x
Valentine LL, Fiedler TL, Hart AN, Petersen CA, Berninghausen HK, Southworth D: Diversity of ectomycorrhizas associated with Quercus garryana in southern Oregon. Can J Botany 2004, 82: 123–135. 10.1139/b03-117
Kõljalg U, Larsson K-H, Abarenkov K, Nilsson RH, Alexander IJ, Eberhardt U, Erland S, Høiland K, Kjøller R, Larsson E, Pennanen T, Sen R, Taylor AFS, Vrålstad T, Tedersoo L, Ursing BM: UNITE – a database providing web based methods for the molecular identification of ectomycorrhizal fungi. New Phytol 2005, 166: 1063–1068. 10.1111/j.1469-8137.2005.01376.x
Rosling A, Landeweert R, Lindahl BD, Larsson K-H, Kuyper TW, Taylor AFS, Finlay RD: Vertical distribution of ectomycorrhizal fungal taxa in a podzol soil profile. New Phytol 2003, 159: 775–783. 10.1046/j.1469-8137.2003.00829.x
Tedersoo L, Kõljalg U, Hallenberg N, Larsson K-H: Fine scale distribution of ectomycorrhizal fungi and roots across substrate layers including coarse woody debris in a mixed forest. New Phytol 2003, 159: 153–165. 10.1046/j.1469-8137.2003.00792.x
Nielsen KB, Kjøller R, Olsson PA, Schweiger PF, Andersen F∅, Rosendahl S: Colonisation and molecular diversity of arbuscular mycorrhizal fungi in the aquatic plants Littorella uniflora and Lobelia dortmanna in southern Sweden. Mycol Res 2004, 108: 616–625. 10.1017/S0953756204000073
Kaldorf M, Renker C, Fladung M, Buscot F: Characterization and spatial distribution of ectomycorrhizas colonizing aspen clones released in an experimental field. Mycorrhiza 2004, 14: 295–306. 10.1007/s00572-003-0266-1
Schüßler A, Schwarzott D, Walker C: Glomeromycota rRNA genes – the diversity of myths? Mycorrhiza 2003, 13: 233–236. 10.1007/s00572-003-0250-9
Álvarez I, Wendel JF: Ribosomal ITS sequences and plant phylogenetic inference. Mol Phylogenet Evol 2003, 29: 417–434. 10.1016/S1055-7903(03)00208-2
Bruns TD, Shefferson RP: Evolutionary studies of ectomycorrhizal fungi: recent advances and future directions. Can J Botany 2004, 82: 1122–1132. 10.1139/b04-021
Hibbett DS, Nilsson RH, Snyder M, Fonseca M, Costanzo J, Shonfeld M: Automated phylogenetic taxonomy: An example in the Homobasidiomycetes (mushroom-forming fungi). Syst Biol, in press.
The Perl Documentation Project[http://www.perldoc.com]
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25: 3389–3402. 10.1093/nar/25.17.3389
The Tree of Life Project[http://tolweb.org/tree/phylogeny.html]
The Apache httpd server project[http://httpd.apache.org/]
Nilsson RH, Rajashekar B, Larsson K-H, Ursing BM: galaxieEST: addressing EST identity through automated phylogenetic analysis. BMC Bioinformatics 2004, 5: 87. 10.1186/1471-2105-5-87
Hawksworth DL: The magnitude of fungal diversity: the 1.5 million species estimate revisited. Mycol Res 2001, 105: 1422–1432.
Berney C, Fahrni J, Pawlowski J: How many novel eukaryotic 'kingdoms'? Pitfalls and limitations of environmental DNA surveys. BMC Biol 2004, 2: 13. 10.1186/1741-7007-2-13
Bridge PD, Roberts PJ, Spooner BM, Panchal G: On the unreliability of published DNA sequences. New Phytol 2003, 160: 43–48. 10.1046/j.1469-8137.2003.00861.x
Vilgalys R: Taxonomic misidentification in public DNA databases. New Phytol 2003, 160: 4–5. 10.1046/j.1469-8137.2003.00894.x
de Queiroz K: Phylogenetic definitions and taxonomic philosophy. Biol Philos 1992, 7: 295–313. 10.1007/BF00129972
Nilsson RH, Larsson K-H, Ursing BM: galaxie – CGI scripts for sequence identification through automated phylogenetic analysis. Bioinformatics 2004, 20: 1447–1452. 10.1093/bioinformatics/bth119
Sequence Altering System[http://www.bork.embl-heidelberg.de/Alerting/]
Whitney J, Esteban DJ, Upton C: Recent hits acquired by BLAST (ReHAB): A tool to identify new hits in sequence similarity searches. BMC Bioinformatics 2005, 6: 23. 10.1186/1471-2105-6-23
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22: 4673–4680.
Tom Bruns, University of California at Berkeley, is acknowledged for initiating the project. The photos of Figure 1 were kindly provided by Leif & Anita Stridvall and Ellen Larsson, respectively. Anders Sjögren is acknowledged for valuable advice on technical matters. The manuscript benefited from helpful comments by Brandon Matheny and Manfred Binder. Financial support to RHN from the Lars Hierta Foundation and from the Helge Ax: son Johnson Foundation are gratefully acknowledged. emerencia was created using only freely available software http://www.gnu.org.
All authors contributed to the structure and functions of emerencia. RHN and EK wrote most parts of the emerencia core, the CGI scripts, and the database handlers. MR was responsible for the mycological part, including literature comparison, integrity testing, and data verification. KHL contributed with advice on fungal taxonomy and systematics. All authors drafted the manuscript and approved the final version.