PHYLOViZ: phylogenetic inference and data visualization for sequence based typing methods
© Francisco et al.; licensee BioMed Central Ltd. 2012
Received: 23 December 2011
Accepted: 8 May 2012
Published: 8 May 2012
With the decrease of DNA sequencing costs, sequence-based typing methods are rapidly becoming the gold standard for epidemiological surveillance. These methods provide reproducible and comparable results needed for a global scale bacterial population analysis, while retaining their usefulness for local epidemiological surveys. Online databases that collect the generated allelic profiles and associated epidemiological data are available but this wealth of data remains underused and are frequently poorly annotated since no user-friendly tool exists to analyze and explore it.
PHYLOViZ is platform independent Java software that allows the integrated analysis of sequence-based typing methods, including SNP data generated from whole genome sequence approaches, and associated epidemiological data. goeBURST and its Minimum Spanning Tree expansion are used for visualizing the possible evolutionary relationships between isolates. The results can be displayed as an annotated graph overlaying the query results of any other epidemiological data available.
PHYLOViZ is a user-friendly software that allows the combined analysis of multiple data sources for microbial epidemiological and population studies. It is freely available at http://www.phyloviz.net.
Microbial typing methods are fundamental tools for epidemiological studies and for the study of bacterial populations. These techniques allow the characterization of bacteria at the strain level providing researchers with important information for the surveillance of infectious diseases, outbreak investigation, pathogenesis studies, natural history of an infection and bacterial population genetics. The recent advances in sequencing technologies and the resulting decrease in costs, are promoting sequence based typing methods over traditional mole- cular methodologies, due to their reproducibility and portability of results, which allows for straightforward inter-laboratory comparison of microbial strains.
Multilocus Sequence Typing (MLST)  and Multilocus Variable Number of Tandem Repeats Analysis (MLVA)  are examples of such methods that are available for a multitude of bacterial species and are being used globally in epidemiological microbial typing and bacterial population studies . Over seventy different microorganisms have public MLST databases available at five different online locations, namely MLST.net , PubMLST , Institut Pasteur MLST Databases , EcMLST  and MLST Databases at the ERI . With the advent of Next Generation Sequencing (NGS), Single Nucleotide Polymorphism (SNP) analysis of entire bacterial genomes is being increasingly used [9–12] and may eventually in the future replace MLST as the typing method of choice in global epidemiological studies.
These advances created the need for algorithms and tools to make sense of this wealth of data in the context of epidemiology, population genetics and evolution. Existing free software tools for data analysis and visualization of microbial typing data are focused on MLST analysis, including the algorithms START , eBURST  and goeBURST . However, these tools lack the capacity to integrate epidemiological data, which is crucial for the correct inference and analysis of relationships between strains. Other well known software tools, provide the ability for network visualization and data integration, usually of generic use, such as Gephi  and Tulip  or more focused on biological networks such as regulatory, protein interaction or metabolic networks, as is the case of Cytoscape . Neither these nor any other current tools for network visualization  are specifically directed towards population or evolutionary analysis and all depend on already defined trees or networks, lacking the ability to infer them directly from data origi- nating from sequence-based typing methods. Moreover, given their focus on general use, their configuration and customization for epidemiological data analysis requires expertise on graph theory and network mining and frequently calls for programming skills, so that their use for this purpose is restricted only to a few advanced users.
PHYLOViZ allows users to manipulate multiple data sets of sequence-based typing data for analysis, which are presented within a window with a tree based navigation menu. Each data set is described by a file containing typing data, i.e., allelic profiles or nucleotide sequences generated by a sequence based typing method such as MLST, MLVA or SNP data. Moreover, a file containing ancillary data, sometimes also refered to as metadata, such as epidemiological, demographic information or any other relevant data may also be added and linked to the allelic profile data, representing information available for different isolates. Both typing data profiles and isolate data can be visualized in a table or a tree view. These views allow the user to query and filter the data using regular expressions or by using a simple point and click interface as described in detail below. Figure 1 displays the basic interface functionalities of the PHYLOViZ interface.
Once the data files are loaded, the user can infer probable patterns of evolutionary descent between allelic profiles using the goeBURST algorithm or its expansion to generalized Minimum Spanning Tree (MST) like structures using different distance metrics. The unrooted trees resulting from the algorithms are displayed using a force-directed layout, where they can be interactively visualized, freely manipulated by the user to explore the existing links and interrogated, either through clicking on the nodes or links of interest to display additional information or by searching the graph for specific nodes.
Users can then combine or query specific characteristics of either typing or ancillary data and overlay the results onto the graphical display. Queries can be easily performed in a tree menu view, by selecting character values to be displayed onto the tree using selection boxes. A regular expression composition box is also provided to allow more complex queries. The resulting annotated trees can then be exported in a variety of common graphic formats, including both vector and bitmap formats. A complete tutorial for PHYLOViZ and a description of its features is available at http://www.phyloviz.net/wiki/tutorial.
The software architecture is based on the NetBeans Platform (available at http://platform.netbeans.org), which includes a plugin system. Users can develop and use new plugins according to their requirements and make them available for the community through their own website or by depositing them at the PHYLOViZ repository.
It is our objective to keep PHYLOViZ as modular and extensible as possible, allowing a seamless integration of new data analysis methods and visualization features. PHYLOViZ is written in Java and is compatible with any operating system running Java 1.6 or higher. The software uses several well-known open source libraries for data export and visualization, notably the Prefuse toolkit for interactive information visualization (available at http://prefuse.org) and the FreeHEP Java libraries for image exporting (available at http://java.freehep.org). The code for PHYLOViZ core and most plugins is available under the GNU General Public License (GPL) version 3 or any later version, with a Classpath Exception. The Classpath Exception means that the copyright holders of a library give you permission to link this library with independent modules to produce an executable, regardless of the license terms of these independent modules. More complete information on the license is provided at http://www.phyloviz.net/wiki/download/. Both the source code and binary versions for major operating systems, as well as the documentation are available at http://www.phyloviz.net.
The current version of PHYLOViZ supports the inference of patterns of evolutionary descent using the goeBURST algorithm  and its extension to a full MST-like approach.
The first improvement over the original implementation of the goeBURST algorithm is that PHYLOViZ allows any possible distance for linkage comparison. The user can choose among available distances, or add new distances by implementing the interface AbstractDistance (see Data Model section) and providing a new plugin as described in Supplemental Material. In order to allow such extensibility, we have identified basic distance requisites. Thus, any implemented distance should provide a way to compare two links (i.e., a link comparator), and it should also compute a level for each link. Note that only a link comparator is required for most MST like approaches and, in particular, for goeBURST. The level is only relevant for link filtering, as is traditionally done for filtering links at different number of locus variants in the case of MLST data analysis. When link levels are not applicable for a given distance measure, we can define the level as being equal for all links.
Secondly, our implementation extends the goeBURST rules to any number of variable loci in order to create a single tree. This approach uses the goeBURST rules to identify, out of all the possible MST-like trees, the one that is concordant with the evolutionary model supporting the rules, guaranteeing a unique and consistent solution. The MST analysis allows the user to dynamically visualize the groups formed at any linkage level, including those formed at single-locus variant (SLV) level (the standard result of goeBURST). The information of the group in which each isolate is classified, for each linkage level, can be exported to the ancillary data file and combined with existing data for later analysis.
The capacity to create groups at any linkage level is particularly useful to analyze MLVA and SNP typing data, since schemata can have much larger number of loci than traditional MLST schemes, frequently relying only on seven loci. The higher number of loci could be expected to yield a larger number of differences between individuals, motivating the exploration of higher linkage levels to generate useful trees translating evolutionary hypotheses.
Note that by allowing any possible distance for link comparison and thereby fully generalizing link comparison, any problem of tree construction becomes a graphical matroid instance [15, 21], that can be solved by well known algorithms, such as the Borǔvka algorithm , the Kruskal algorithm  or the Prim algorithm . For goeBURST, we enumerated explicitly all possible links, because we computed tie statistics for them, and we built the tree using the Kruskal algorithm. However, for the full MST computation, we are now using the Prim algorithm. The main reason for this is to avoid enumerating explicitly all links, the number of which, in the generalized version, grows quadratically with the number of nodes. The Prim algorithm only requires us to generate links implicitly, involving at most to keep in memory as many links as the number of nodes and, thus, making the approach more scalable for large data sets. Moreover, since we may have dense graphs, we also implemented a more efficient priority queue for Prim’s algorithm, based on rank relaxed heaps as proposed by Driscoll et al. . These improvements are part of the PHYLOViZ Algorithms plugin. The visualization of the full MST was also optimized to support large data sets, with thousands of nodes. This was made possible by incrementally adding nodes to the force directed display following a Breadth First Search approach starting from the node with greater number of links. This lead to increased optimization speed and lowered the memory requirements for the MST display.
Creating new plugins
The PHYLOViZ plugin system allows new algorithms and visualization tools to be easily added. The implementation of these features must also use the NetBeans Platform respecting the documentation regarding the creation of new modules available at http://netbeans.org/kb/trails/platform.html.
In order to have access to the visualization and querying capabilities of PHYLOViZ, the developer must add the intended PHYLOViZ plugins as library dependencies of the new plugin. Extensive documentation, describing the data model and algorithms API, is available on the website and in the Data Model section. A plugin implementation example is provided in Supplemental Material.
The graphical user interface for the incorporation of newly developed plugins is found in the Plugins entry of the Tools menu. The user has the option of adding a plugin stored locally, or through the use of a remote Update Center. The latter is periodically checked for new versions of the installed plugins and the user notified in case of the availability of updates. The update center can also be used to install new publicly available plugins from within PHYLOViZ.
Results and Discussion
PHYLOViZ aims to be an efficient and user-friendly software, providing an interface where users can analyze and visualize their data, and export the results in several commonly used graphic formats.
To achieve these goals, it allows users to create data sets by loading the data from tab separated files for any kind of allelic profile data (be it MLST, MLVA or SNP data, i.e. Typing Data), that can be easily exported from Excel worksheets, or FASTA formatted files for sequence data, a common format easily obtained from any sequence analysis editor software. Nevertheless, as previously pointed out , the use of sequence information to infer relatedness between microbial isolates can be confounded by recombination that may introduce multiple nucleotide changes in a single event. On the other hand for microbial species where recombination is infrequent or for chromosomic regions where recombination is rare, this may be an approach that can be pursued. The users can also use the available plugins to download data directly from public databases that provide their data through a Web service interface such as PubMLST.org  or as simple flat files. This flexibility, coupled with the ability to use the goeBURST algorithm and its expansion to a MST, allows the user to load multiple data sets simultaneously for the analysis, visual exploration and comparison of multiple typing results for the same isolate collection. Any available additional ancillary clinical or typing data can be similarly easily exported from Excel worksheets or from other database software, or directly retrieved from PubMLST.org from within PHYLOViZ (see below).
Also, by applying the same analysis algorithms to different typing techniques, namely MLST, MLVA or SNP typing, PHYLOViZ provides researchers with a more direct comparison of the different typing methodologies results, since using the same underlying evolutionary models and assumptions avoids possible artifacts or differences resulting from using multiple data analysis techniques for each typing method and then comparing the results. For the analysis of MLVA data, an alternative distance based on goeBURST is also made available in the software. In this implementation, the algorithm still creates trees only linking SLVs, however instead of simply using the goeBURST rules to decide which links should be drawn, it starts by considering the euclidean distance between two profiles. This distance is calculated as the sum of the differences across all loci in the MLVA profile. For instance in a two locus scheme the difference between profiles 1-3 and 1-1 would be calculated as 2 in this metric whereas the profiles would be SLVs in the classic MLST distance. This assumes that the genetic distance between two MLVA profiles is proportional to the difference in the number of repeats at each locus. This is in accordance with models of accumulation of tandem repeats in each loci [26–28], and for some data sets could create a more realistic and meaningful representation than using only the goeBURST rules.
The ability to analyze data sets comprising thousands of isolates, with thousands of target loci, also allows PHYLOViZ to efficiently handle the large volumes of data generated by NGS when applied to molecular epidemiology. This, coupled with its ability for data integration with a simple query system and dynamic visualization capabilities provide a user friendly platform aimed at non-bioinformaticians such as clinical microbiologists, epidemiologists and infection control specialists.
Phyloviz usage example
To demonstrate a typical usage of PHYLOViZ for exploring sequence-based typing methods databases and related epidemiological data, we show the analysis of the MLST database from Streptococcus pneumoniae retrieved from https://spneumoniae.mlst.net. This data set (retrieved in 28 October 2011) contains 7146 unique sequence types and information on 15329 isolates. The available data on penicilin susceptibility was curated and recoded for R (resistant), I (intermediate susceptibility) and S (susceptible) based on the Clinical and Laboratory Standards Institute breakpoints for oral penicillin  and saved on a tab delimited file. In subsequent analysis isolates classified as R or I are collectively referred to as non-susceptible. A new data set was created on PHYLOViZ were 7146 unique STs and respective allelic profiles were used as Typing Data and the penicillin resistance information was used as Isolate Data. The goeBURST algorithm  was computed at SLV level for the data set dividing it in 506 groups of two or more STs and 1514 singletons (STs unrelated to any other in the collection at SLV level). The largest clonal complex (CC), currently having 986 unique STs and representing 2260 isolates is presented in Figure 3.
A screencast to demonstrate how fast and easy it is to produce this comparison is available at http://www.phyloviz.net/wiki/videos/.
Remote MLST public repositories
The PHYLOViZ LoadMLST DBs plugin allows PHYLOViZ to load allelic profiles and sequence data from public MLST repositories. It is currently capable of retrieving data from five public repositories: PubMLST , MLST.net , MLST Databases at ERI , Institut Pasteur MLST Databases  and EcMLST . These repositories make their data available through flat files. Additionally, this plugin allows the integration of a local file, containing the isolate ancillary data, to complement the retrieved allelic profiles.
As described in Additional file 1, each developed plugin must specify which are the necessary dependencies. In this case, the PHYLOViZ Core, PHYLOViZ MLST and PHYLOViZ Aligned Sequences need to be specified, since the plugin permits the retrieval of allelic profiles and sequence data.
The PHYLOViZ LoadMLST DBs plugin creates an entry under the File menu, to launch a wizard-like interface guiding the user through the loading process. The interface is composed of four distinct steps. The first, presents the user with the list of all the possible data sets available in the five public repositories, together with a mandatory text field where the user specifies the name of the selected data set. The second, presents the full list of genes constituting the allelic profile of the selected data set, as well as the number of STs composing it. The complete set of allelic profiles can be downloaded by a single click. In the third step, the user can integrate with the retrieved allelic profiles, a local file containing isolate ancillary data. Finally, in the last step of the wizard, the user is presented with two choices: the immediate creation of a data set without any sequence data, or to add the complete sequence data for each loci. In the latter, the user can import the sequence data independently for each locus, either by downloading it directly from the public repository or by specifying a local file in FASTA format. By pressing Finish, the new data set is created with the selected typing data.
Since both allelic profiles and sequence data are specific instances of TypingData, both the goeBURST algorithm and its MST-like approach can be directly applied. This means that PHYLOViZ permits the inference of patterns of evolutionary descent, not only from the MLST data but also directly from the sequence data of the selected data set (see Algorithms subsection).
A similar plugin (PHYLOViZ PubMLST) takes advantage of the PubMLST Web services to also automatically retrieve isolate ancillary data from the database. Both the PHYLOViZ LoadMLST DBs and the PHYLOViZ PubMLST plugins allow PHYLOViZ to directly access, analyse and visualize the latest epidemiological data available in remote MLST public repositories.
PHYLOViZ provides the research community with a freely available and easy to use software tool for the analysis of sequence-based typing methods. It is in continuous development, focusing on user-friendliness, platform independence, modularity, reusability and efficiency and offering a plugin architecture gives users the possibility to create modules tailored to specific needs. Its flexibility in terms of epidemiological typing data analysis, relies on its ability to use any data format currently used for the established methods, such as MLST and MLVA, as well as those of upcoming methodologies based on NGS.
We are currently working on new analysis methods and improved visualization capabilities in order to support and efficiently analyze the continuously growing data sets. We are also extending the software to directly interact with online databases, with the ultimate goal of providing a unified and simple interface for users to interact and explore currently available data.
Availability and requirements
∙Project name: PHYLOViZ ∙Project home page:http://www.phyloviz.net/∙Operating system(s): Platform independent ∙Programming language: Java ∙Other requirements: Java 1.6 or higher ∙License: GNU GPL v3 w/ Classpath Exception ∙Any restrictions to use by non-academics: Only those imposed already by the license.
This work was partly supported by Fundação para a Ciencia e Tecnologia, Portugal, projects PPTDC/SAU-ESA/71499/2006 and PTDC/EIA-CCO/118533/2010, and INESC-ID multiannual funding through the PIDDAC Program funds, by Fundação Calouste Gulbenkian and by an unrestricted grant from GlaxoSmithKline.
- Maiden MCJ, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA, et al.: Multilocus sequence typing: A portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Nat Acad Sci 1998, 95(6):3140. 10.1073/pnas.95.6.3140PubMed CentralView ArticlePubMedGoogle Scholar
- Lindstedt BA: Multiplelocus variable number tandem repeats analysis for genetic fingerprinting of pathogenic bacteria. Electrophoresis 2005, 26(13):2567. 10.1002/elps.200500096View ArticlePubMedGoogle Scholar
- van Belkum A, Struelens M, de Visser A, Verbrugh H, Tibayrenc M: Role of genomic typing in taxonomy, evolutionary genetics, and microbial epidemiology. Clin Microbiol Rev 2001, 14(3):547. 10.1128/CMR.14.3.547-560.2001PubMed CentralView ArticlePubMedGoogle Scholar
- MLST: Multi Locus Sequence Typing. Imperial College of London;http://www.mlst.net
- PubMLST. University of Oxford (UK);http://pubmlst.org/
- Institut Pasteur MLST Databases. Pasteur Institute. . [http://www.pasteur.fr/mlst/] Pasteur Institute.
- EcMLST. Michigan State University. . [http://www.shigatox.net/cgi-bin/mlst7/index] Michigan State University.
- MLST Databases at ERI. University College Cork. . [http://mlst.ucc.ie/] University College Cork.
- Roumagnac P, Weill F, Dolecek C, Baker S, Brisse S, Chinh N, Le T, Acosta C, Farrar J, Dougan G: Evolutionary history of Salmonella typhi. Science 2006, 314(5803):1301. 10.1126/science.1134933PubMed CentralView ArticlePubMedGoogle Scholar
- Harris SR, Feil EJ, Holden MTG, Quail MA, Nickerson EK, Chantratita N, Gardete S, Tavares A, Day N, Lindsay JA, Edgeworth JD, de Lencastre H, Parkhill J, Peacock SJ, Bentley SD: Evolution of MRSA During Hospital Transmission and Intercontinental Spread. Science 2010, 327(5964):469. 10.1126/science.1182395PubMed CentralView ArticlePubMedGoogle Scholar
- Holt KE, Parkhill J, Mazzoni CJ, Roumagnac P, Weill FX, Goodhead I, Rance R, Baker S, Maskell DJ, Wain J, Dolecek C, Achtman M, Dougan G: High-throughput sequencing provides insights into genome variation and evolution in Salmonella Typhi. Nat Genet 2008, 40(8):987. 10.1038/ng.195PubMed CentralView ArticlePubMedGoogle Scholar
- Croucher N, Harris S, Fraser C, Quail M: Rapid pneumococcal evolution in response to clinical interventions. Science 2011, 331: 430. 10.1126/science.1198545PubMed CentralView ArticlePubMedGoogle Scholar
- Jolley KA, Feil EJ, Chan MS, Maiden MCJ: Sequence type analysis and recombinational tests (START). Bioinformatics 2001, 17(12):1230. 10.1093/bioinformatics/17.12.1230View ArticlePubMedGoogle Scholar
- Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG: eBURST: Inferring Patterns of Evolutionary Descent among Clusters of Related Bacterial Genotypes from Multilocus Sequence Typing Data. J Bacteriol 2004, 186(5):1518. 10.1128/JB.186.5.1518-1530.2004PubMed CentralView ArticlePubMedGoogle Scholar
- Francisco AP, Bugalho MF, Ramirez M, Carriço JA: Global optimal eBURST analysis of multilocus typing data using a graphic matroid approach. BMC Bioinf 2009., 10(152):Google Scholar
- Bastian M, Heymann S, Jacomy M: Gephi: An open source software for exploring and manipulating networks. International AAAI Conference Weblogs Social Media 2009. [http://www.aaai.org/ocs/index.php/ICWSM/09/paper/view/154]Google Scholar
- Auber D: Tulip-a huge graph visualization framework. In Graph Drawing Software. Springer Springer, ; 2003:105–126.Google Scholar
- Shannon P, Markiel A, Ozier O, Baliga N, Wang J, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 13(11):2498. 10.1101/gr.1239303PubMed CentralView ArticlePubMedGoogle Scholar
- Suderman M, Hallett M: Tools for visually exploring biological networks. Bioinformatics 2007, 23(20):2651. 10.1093/bioinformatics/btm401View ArticlePubMedGoogle Scholar
- Shneiderman B: The Eyes Have It: A Task by Data Type Taxonomy for Information Visualizations. In IEEE Symposium on Visual Languages. , ; 1996:336–343.Google Scholar
- Papadimitriou CH, Steiglitz K: Combinatorial Optimization. Dover, ; 1998.Google Scholar
- Borǔvka O: On a minimal problem. Prace Moraské Pridovedecké Spolecnosti 1926., 3:Google Scholar
- Kruskal JB: On the shortest spanning subtree of a graph and the traveling salesman problem. Proc Am Math Soc 1956, 7: 48. 10.1090/S0002-9939-1956-0078686-7View ArticleGoogle Scholar
- Prim RC: Shortest connection networks and some generalizations. Bell Syst Tech J 1957, 36: 1389.View ArticleGoogle Scholar
- Driscoll JR, Gabow HN, Shrairman R, Tarjan RE: Relaxed Heaps: An Alternative to Fibonacci Heaps with Applications to Parallel Computation. Commun ACM 1988, 31(11):1343. 10.1145/50087.50096View ArticleGoogle Scholar
- Xu X, Peng M, Fang Z: The direction of microsatellite mutations is dependent upon allele length. Nat Genet 2000, 24: 396. 10.1038/74238View ArticlePubMedGoogle Scholar
- Vogler A, Keys C, Allender C, Bailey I, Girard J, Pearson T, Smith K, Wagner D, Keim P: Mutations, mutation rates, and evolution at the hypervariable VNTR loci of Yersinia pestis. Mutation Research/Fundamental Mol Mech Mutagenesis 2007, 616(1–2):145. 10.1016/j.mrfmmm.2006.11.007View ArticleGoogle Scholar
- Cooley MB, Carychao D, Nguyen K, Whitehand L, Mandrell R: Effects of Environmental Stress on Stability of Tandem Repeats in Escherichia coli O157:H7. Appl Environ Microbiol 2010, 76(10):3398. 10.1128/AEM.02373-09PubMed CentralView ArticlePubMedGoogle Scholar
- Clinical, Institute LS, Clinical: Performance Standards for Antimicrobial Susceptibility Testing - Sixteenth Informational Supplement, Volume CLSI document M100-S18. Wayne, PA: Clinical and Laboratory Standards Institute 2008 Wayne, PA: Clinical and Laboratory Standards Institute 2008Google Scholar
- Bergmann C, Chi F, Rachid S, Hakenbeck R: Mechanisms for penicillin resistance in Streptococcus pneumoniae: penicillin binding proteins, gene transfer, and cell wall metabolism. In The pneumococcus. Edited by: Tuomanen EI, Mitchell TJ, Morrison DA, Spratt BG. ASM Press, Washington, DC; 2004:339–349.Google Scholar
- Enright MC, Spratt BG: Extensive variation in the ddl gene of penicillin-resistant Streptococcus pneumoniae results from a hitchhiking effect driven by the penicillin-binding protein 2b gene. Mol Biol Evol 1999, 16(12):1687. 10.1093/oxfordjournals.molbev.a026082View 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.