Phylometrics: a pipeline for inferring phylogenetic trees from a sequence relationship network perspective
© Smits and Ouverney; licensee BioMed Central Ltd. 2010
Published: 7 October 2010
Comparative sequence analysis of the 16S rRNA gene is frequently used to characterize the microbial diversity of environmental samples. However, sequence similarities do not always imply functional or evolutionary relatedness due to many factors, including unequal rates of change and convergence. Thus, relying on top BLASTN hits for phylogenetic studies may misrepresent the diversity of these constituents. Furthermore, attempts to circumvent this issue by including a large number of BLASTN hits per sequence in one tree to explore their relatedness presents other problems. For instance, the multiple sequence alignment will be poor and computationally costly if not relying on manual alignment, and it may be difficult to derive meaningful relationships from the resulting tree. Analyzing sequence relationship networks within collective BLASTN results, however, reveal sequences that are closely related despite low rank.
We have developed a web application, Phylometrics, that relies on networks of collective BLASTN results (rather than single BLASTN hits) to facilitate the process of building phylogenetic trees in an automated, high-throughput fashion while offering novel tools to find sequences that are of significant phylogenetic interest with minimal human involvement. The application, which can be installed locally in a laboratory or hosted remotely, utilizes a simple wizard-style format to guide the user through the pipeline without necessitating a background in programming. Furthermore, Phylometrics implements an independent job queuing system that enables users to continue to use the system while jobs are run with little or no degradation in performance.
Phylometrics provides a novel data mining method to screen supplied DNA sequences and to identify sequences that are of significant phylogenetic interest using powerful analytical tools. Sequences that are identified as being similar to a number of supplied sequences may provide key insights into their functional or evolutionary relatedness. Users require the same basic computer skills as for navigating most internet applications.
The breadth of knowledge of microbial diversity continues to rapidly expand as 16S rRNA genes are sequenced from environmental samples and comparisons to existing data are drawn. Culture-independent methods have enabled the application of a molecular-phylogenetic approach to discover that the environmental microbes constitute the majority of evolutionary diversity . Phylogenetic studies have further recognized that over 99% of all prokaryotes are uncultivated [2, 3]. Therefore, in attempts to continue to explore the evolution and composition of microbial communities it is now standard practice to sequence the 16S rRNA gene .
These DNA sequences are being added to public databases rapidly since they have become the most cost effective, if not the only method available to identify and quantify the uncultivated microbes. Deriving meaningful knowledge from the inundation of sequences, then, necessitates methods to judge which sequences are of particular interest to the study. A common method in screening DNA sequences derived in a study is to assign each sequence to a taxonomic group by comparing it to the closest relative in publically available databases, such as Greengene’s Simrank  and NCBI’s BLASTN . Although this method is rapid, it has been recognized that neither the top ranking hits nor the most similar sequences are always the most closely related phylotype [7, 8].
Phylotypes may also be determined to be closely related by examining the network of sequence relationships based on their similarities . A network is established using the following principle: if sequences A, B, and C each return sequence a as a BLASTN hit (with varying ranks), sequence a is considered to be “associated” with the three sequences. In the context of screening the sequences derived from environmental samples, analyzing the associations between the BLASTN hits of a collection will unveil sequences of great phylogenetic interest despite the possibility of having scored low in the BLASTN ranking system . Recognizing these associations before aligning large collections of sequences reduces both the computational expense and the obfuscation that many sequences will inherently contribute to its analysis.
Here we introduce a web application, Phylometrics, which automates many key tasks for DNA sequence analysis in a high-throughput method. Instead of single BLASTN inquiries, Phylometrics automatically analyzes multiple BLASTN results, cross-comparing each BLASTN result against others within seconds. A graphical tool for mining these sequence relationships is also provided and includes them in a multiple sequence alignment process according to user-defined parameters. This approach can discover sequences of great interest for further phylogenetic examination which would otherwise be missed without significant human involvement. Finally, the application generates phylogenetic trees according to numerous methods of tree inference which may be selected and configured using a wizard-style form. As the results are being generated, the user may continue to interact and queue additional analyses without interrupting the ongoing processes.
Phylometrics offers a number of key advantages: simplicity in the user interface, powerful analytics to examine sequence relationships, the ability to install the application locally and the ability to queue batch jobs. The code is also released as open source allowing it to be extended and integrated into workflows as desired. Finally, due to the architectural decision to make the application driven by job queues (a list of tasks to be run asynchronously), cloud computing services may be integrated and directed by the queue to perform the more CPU-intensive processes such as multiple sequence alignments.
Implement a simple, widely usable interface that does not involve user programming, by implementing a wizard-styled approach.
Provide rich, interactive visual representations of the data generated by the pipeline.
Allow for the import of sequences in multiple formats and the export of data at every stage of the pipeline so as to give the user the flexibility to use alternative software solutions for those steps. Furthermore, the system must maintain a record of the parameters the user set for each step in the pipeline.
Allow batch jobs to be queued and/or running while being able to simultaneously create new jobs or able to access completed reports.
User interface and pipeline wizard
A common web browser serves as the user interface. It is advantageous due to the low learning curve and accessibility from a number of devices. It also incorporates simple yet useful features such as printing and searching from within the browser experience.
The second step of the Wizard is “BLAST Setup.” The software queries NCBI’s live database by default; this may be modified, however, to search a local BLASTN database instead. After entering BLASTN parameter options, the user has the ability to “Generate a Screening Report,” or to continue with the pipeline. Both processes are queued and allow the user to continue to use the system, while the option to generate the screening report will halt the pipeline at this second step before continuing the progression through the pipeline. A Screening Report will be generated irrespective of the choice the user makes at this point, allowing one to revisit this report to create any number of alignments and trees based on their filters.
The following step is “Alignment” with the options to select parameters based on ClustalW , a multiple sequence alignment program. Once the Screening Report is generated, users may choose to create alignments using various filters which set thresholds, gate and also specifically toggle individual sequences. In fact, one can use the same dataset to create many alignments using different parameters and compare the results in the generated trees.
The final step of the Wizard is “Build Tree”, where one can choose between various methods of tree inferences based on the PHYLIP  and PhyML software packages. The NEWICK file format outputs of these packages are then rendered into images using TreeVector . Similar to the many-to-one relationship of alignments to collections, many trees can be generated from one alignment.
Interactive visual analytics
Sequences are initially imported by the user into a “Collection” from the Wizard in FASTA, Genbank or Raw format. On Wizard completion, the user is be able to download the BLASTN report in XML format, the multiple sequence alignments in PHYLIP format, and the tree inferences in NEWICK, SVG and PNG file formats. The user may also click on the generated results to review the parameters that were used to generate the reports within the Wizard process.
Asynchronous job queues
The pipeline processing and job queues are managed by Perl scripts which utilize the BioPerl library . These scripts are decoupled from the system in such a way so as to allow them to be run on a completely separate hardware platform, receiving job instructions from the database layer. This setup allows the user to continue to access and analyze generated reports while jobs are processing, in addition to being able to queue an unlimited number of jobs. As the tasks progress through the pipeline an unobtrusive panel displays the progress of each submitted job. Furthermore, this setup can be exploited to process jobs which are computationally expensive on other servers simultaneously or even utilize cloud computing services such as Amazon EC2.
Results and discussion
Comparison of Phylometrics to existing commonly-used tools for generating phylogenetic trees.
Where is it installed?
Local or Remote
Offers Manual Alignment
Visualization of Sequence Networks
Phylometrics makes it possible to explore sequence relationships within a collection in order to quickly build phylogenetically meaningful trees. Using a visual and interactive interface, sequences that may provide key insights into functional or evolutionary relatedness are identified much faster than the current standard manual approach. Furthermore, the identification of homologous sequences previous to multiple sequence alignments and tree building significantly reduces computational costs and time by minimizing the inclusion of unrelated sequences.
Once a dataset has been refined and pertinent sequences have been identified, Phylometrics enables one to explore the differences between alignments and each method of tree inference to determine which set of methods and parameters may suit the data best. The parameters used in each of the processes are also permanently stored, allowing one to replicate the results manually.
Results can be exported from any stage in the pipeline for personal collection databases or to be further analyzed in other popular sequence analysis packages, such as ARB . The formats of the output include: Phylip for alignments, Newick for tree inferences and JPG and SVG formats for tree visualizations. The multiple sequence alignments are especially useful for export, as manually aligning sequences may be necessary to improve the quality of the resulting tree. The visual representations are also printable from within the Flash movie by accessing the menu with a right-click.
Phylometrics offers an automated, high-throughput pipeline for phylogenetic inferences in a simple wizard-style web interface, requiring basic computer skills and no user programming experience. It incorporates novel analytical tools to mine sequence collections for biologically significant associations generating trees that provide insights into their respective functional or evolutionary implications.
Availability and requirements
Project name: Phylometrics
Project home page: http://www.phylometrics.com
Operating system: Platform independent
We thank the four anonymous reviewers for their time in providing constructive comments on the manuscript. This work has been funded in part by NIH Grant #SC3GM082291 and the Howard Hughes Medical Institute Science Education Grant # 52006312.
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 6, 2010: Proceedings of the Seventh Annual MCBIOS Conference. Bioinformatics: Systems, Biology, Informatics and Computation. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S6.
- Pace NR: A Molecular View of Microbial Diversity and the Biosphere. Science 1997, 276: 734–740. 10.1126/science.276.5313.734View ArticlePubMedGoogle Scholar
- Hugenholtz P: Exploring prokaryotic diversity in the genomic era. Genome Biol 2002, 3: REVIEWS0003. 10.1186/gb-2002-3-2-reviews0003PubMed CentralView ArticlePubMedGoogle Scholar
- Rappé MS, Giovannoni SJ: The uncultured microbial majority. Annu. Rev. Microbiol 2003, 57: 369–394. 10.1146/annurev.micro.57.030502.090759View ArticlePubMedGoogle Scholar
- Imhoff JF, Caumette P: Recommended standards for the description of new species of anoxygenic phototrophic bacteria. Int J Syst Evol Microbiol 2004, 54: 1415–1421. 10.1099/ijs.0.03002-0View ArticlePubMedGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol 2006, 72: 5069–5072. 10.1128/AEM.03006-05PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul S, Madden T, Schaffer A, Zhang J, Zhang Z, Miller W, Lipman D: 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.3389PubMed CentralView ArticlePubMedGoogle Scholar
- Koski LB, Golding GB: The Closest BLAST Hit Is Often Not the Nearest Neighbor. Journal of Molecular Evolution 2001, 52: 540–542.View ArticlePubMedGoogle Scholar
- Eisen JA: Phylogenomics: Improving Functional Predictions for Uncharacterized Genes by Evolutionary Analysis. Genome Research 1998, 8: 163–167.View ArticlePubMedGoogle Scholar
- Suderman M, Hallett M: Tools for visually exploring biological networks. Bioinformatics 2007, 23: 2651–2659. 10.1093/bioinformatics/btm401View ArticlePubMedGoogle Scholar
- Fokkens L, Botelho S, Boekhorst J, Snel B: Enrichment of homologs in insignificant BLAST hits by co-complex network alignment. BMC Bioinformatics 2010, 11: 86. 10.1186/1471-2105-11-86PubMed CentralView ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics 1989, 5: 166. 164 164Google Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol 2003, 52: 696–704. 10.1080/10635150390235520View ArticlePubMedGoogle Scholar
- Larkin M, Blackshields G, Brown N, Chenna R, McGettigan P, McWilliam H, Valentin F, Wallace I, Wilm A, Lopez R, Thompson J, Gibson T, Higgins D: Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23: 2947–2948. 10.1093/bioinformatics/btm404View ArticlePubMedGoogle Scholar
- Pethica R, Barker G, Kovacs T, Gough J: TreeVector: Scalable, Interactive, Phylogenetic Trees for the Web. PLoS ONE 2010, 5: e8934. 10.1371/journal.pone.0008934PubMed CentralView ArticlePubMedGoogle Scholar
- Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C, Fuellen G, Gilbert JGR, Korf I, Lapp H, Lehväslaiho H, Matsalla C, Mungall CJ, Osborne BI, Pocock MR, Schattner P, Senger M, Stein LD, Stupka E, Wilkinson MD, Birney E: The Bioperl toolkit: Perl modules for the life sciences. Genome Res 2002, 12: 1611–1618. 10.1101/gr.361602PubMed CentralView ArticlePubMedGoogle Scholar
- Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar BuchnerA, Lai T, Steppi S, Jobb G, Forster W, Brettske I, Gerber S, Ginhart AW, Gross O, Grumann S, Hermann S, Jost R, Konig A, Liss T, Lussmann R, May M, Nonhoff B, Reichel B, Strehlow R, Stamatakis A, Stuckmann N, Vilbig A, Lenke M, Ludwig T, Bode A, Schleifer K: ARB: a software environment for sequence data. Nucleic Acids Res 2004, 32: 1363–1371. 10.1093/nar/gkh293PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol. Biol. Evol 2007, 24: 1596–1599. 10.1093/molbev/msm092View ArticlePubMedGoogle Scholar
- Wu D, Hartman A, Ward N, Eisen JA: An Automated Phylogenetic Tree-Based Small Subunit rRNA Taxonomy and Alignment Pipeline (STAP). PLoS ONE 2008, 3: e2566. 10.1371/journal.pone.0002566PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.