PanACEA: a bioinformatics tool for the exploration and visualization of bacterial pan-chromosomes
BMC Bioinformatics volume 19, Article number: 246 (2018)
Bacterial pan-genomes, comprised of conserved and variable genes across multiple sequenced bacterial genomes, allow for identification of genomic regions that are phylogenetically discriminating or functionally important. Pan-genomes consist of large amounts of data, which can restrict researchers ability to locate and analyze these regions. Multiple software packages are available to visualize pan-genomes, but currently their ability to address these concerns are limited by using only pre-computed data sets, prioritizing core over variable gene clusters, or by not accounting for pan-chromosome positioning in the viewer.
We introduce PanACEA (Pan-genome Atlas with Chromosome Explorer and Analyzer), which utilizes locally-computed interactive web-pages to view ordered pan-genome data. It consists of multi-tiered, hierarchical display pages that extend from pan-chromosomes to both core and variable regions to single genes. Regions and genes are functionally annotated to allow for rapid searching and visual identification of regions of interest with the option that user-supplied genomic phylogenies and metadata can be incorporated. PanACEA’s memory and time requirements are within the capacities of standard laptops. The capability of PanACEA as a research tool is demonstrated by highlighting a variable region important in differentiating strains of Enterobacter hormaechei.
PanACEA can rapidly translate the results of pan-chromosome programs into an intuitive and interactive visual representation. It will empower researchers to visually explore and identify regions of the pan-chromosome that are most biologically interesting, and to obtain publication quality images of these regions.
Next-generation sequencing technologies and a realization that single reference genomes are insufficient to grasp species-level diversity have resulted in a phenomenal rise in the number of publicly available bacterial genome sequences. A comparison of just six strains of Streptococcus agalactiae demonstrated that many more isolates are needed to capture strain diversity and helped define the concept of the bacterial pan-genome: the set of genes (core and variable) that are encoded within a bacterial species . Tools have been developed to perform multiple genome comparisons by computing orthologous gene clusters and the resulting sets of core and variable genes [2,3,4,5,6,7,8,9,10]. Chan et al. extended the pan-genome concept to the “pan-chromosome”, where the order and orientation of core genes produce a consensus circular scaffold; thus, providing the framework for placing variable genes into discrete “flexible genomic regions (fGRs)” . It is these fGRs that help define phenotypic subspecies differences  and provide the means for survival under iron limiting conditions, host immune pressure, and antibiotics .
To facilitate the interpretation of results for biological discovery, visualization tools have been developed, but still suffer from a number of caveats. A subset of pan-genome visualization tools are web-based (which is good for human intuitive data representation, but poses costly overhead), but only work with pre-computed and/or static data and do not allow user-supplied sequence data [13,14,15,16,17]. Pan-Tetris  and PanViz  are both interactive, but do not easily display variable (a.k.a., flexible) genomic islands (fGIs) . Some visualization tools focus on alignments of core regions , require complicated database dependencies or produce complicated network diagrams . None of the existing pan-genome visualization tools are geared toward a standalone (i.e., client side), intuitive, pan-chromosome-based interactive browser that will enable researchers to navigate to those parts of the pan-genome that are most relevant to understanding strain-specific differences that may impact pathogenesis, antimicrobial resistance, and general fitness in a given environment.
PanACEA uses PERL scripts and a tab-delimited human-readable flat file that contains the following necessary information for the script to generate platform-independent visualizations: the gene order of the pan-chromosome “assemblies”, including the flexible and core regions (such as output of gene_order.pl ); detailed information about each gene; and the location of the sequences of the genes. Though this file can be recreated ad hoc and the user manual does provide descriptions, the PanACEA software package includes a script designed to translate the output of pan-genome software packages to the PanACEA flat file (Fig. 1). Currently, PanACEA must be downloaded or cloned from the GitHub site and run locally. As such, the flat file input provides flexibility for the user independent of which pan-genome generation software they wish to use, both current and future programs. Currently, PanACEA optimally works with PanOCT  and gene_order.pl  output (both are availible at https://sourceforge.net/projects/panoct/). An example dataset consisting of the PanOCT and gene_order.pl derived pangenome of 19 Acinetobacter baumannii genomes along with GO term and ARO term based gene annotations is also available at the PanACEA GitHub repository.
Beyond generic input requirements, PanACEA is highly configurable, allowing for customization of input features specific to the needs and available data of the researcher. Additional information, such as that describing the functionality of the genes or the relationship between genomes, can be incorporated (Fig. 1). Any functional annotation (i.e., Gene Ontology (GO) [23, 24] or Antibiotic Resistance Ontology (ARO)  terms) can be added modularly through a configuration file that will associate colors with functional annotation as well as ontology information. Included with the package are scripts that will add annotation to the gene clusters in a format that PanACEA can read. For sets of genomes with a known evolutionary relationship, a Newick-formatted phylogenetic tree file can also be added, along with metadata information about the genomes such as isolation date, host, serotype, pathogen/non-pathogen, etc.
The PanACEA interface enables the interactive exploration of pan-genomic data through multiple spatial views, from broad pan-chromosome/scaffold context through multi-gene regions to single gene details (Additional file 1: Figure S1). Pan-scaffold representations can be cyclic or linear and highlight flexible and core regions, with core genes individually colored by protein function. For cyclic representations, the nucleotide position coordinate system of the consensus pan-chromosome is used. The pan-scaffolds are shown at identical heights, independent of the number of genomes found in each region. For ease of differentiating short flexible and core regions, the flexible regions are all shown at staggered instances of three-quarters height, again regardless of how many genomes are contained in that region. Regions of interest, such as those involved in antibiotic resistance, virulence, bacteriophage, plasmid, or any other user-configured high-level feature can be preferentially displayed. Likewise, the pan-scaffold (main) page contains a table listing regions, genes, and specific functional terms and can be selected to also highlight the location of the genes. The main page includes a text search function to facilitate identifying specific genes and regions in the table and a zoom function on the top of the main page. The user can scale from the pan-scaffold to a more detailed view of single regions, whether a set of core genes or a fGR, either by clicking on the region on the pan-scaffold map or in the table. On separate pages, PanACEA provides a linear representation of gene context, associated functional annotation, and prevalence of the region in each genome. Given the possible complexity of a fGR, the display can be trimmed to focus on a reduced set of fGIs of interest. Additionally, when included, the genomic phylogeny, accessible from the fGR and core region pages, as well as the gene pages, enables phylogenomic analysis of any region of interest overlaid with user-provided metadata. This functionality can be extended to individual gene summary pages, which display gene annotation and provide access to sequence data and single gene analysis tools such as multiple sequence alignments. All PanACEA displays can be exported as publication-quality SVGs or preview graphics files in other formats (e.g., PNG) and the gene and region lists in tabular data as text files.
A more detailed description of both the PanACEA software package and the web pages with the visualization, complete with examples and help pages, is available in the PanACEA manual on the GitHub site.
The biological utility and output of PanACEA is illustrated using the Enterobacter hormaechei pan-genome data generated from PanOCT from 219 genomes where PanACEA helped to visualize fGIs responsible for the known metabolic differences historically used to classify E. hormaechei subspecies . The time to generate all necessary files from the PanOCT output to the final web pages was 466 s. In addition to the pan-genome, annotation files for each of the gene clusters calculated using GO terms and anti-microbial resistance genes from the CARD database using RGI were used [24, 25]. All the E. hormaechei PanACEA files are available on the GitHub site. The fGR depicted contains two GIs (one flexible and one core between core gene clusters 3936 and 3949) and encodes metabolic pathways historically used to define phenotypic differences between E. hormaechei subspecies (Fig. 2). E. hormaechei subsp. hormaechei is distinguishable from E. hormaechei subsp. oharae and E. hormaechei subsp. steigerwaltii by growth on dulcitol (a.k.a. galactitol) as the sole carbon source via the gat operon . In contrast, E. hormaechei subsp. oharae and subsp. steigerwaltii both encode a different fGI (the aga operon) for the metabolism of N-acetylgalactosamine  (Fig. 2). We readily identified and located the genes and regions of interest by inputting “N-acetylgalactosamine” in the text search and selecting the highlighted regions and genes of interest in the main pan-chromosome view as shown in Fig. 2, thus allowing for analysis of the positional context. The output demonstrates the capability of PanACEA to highlight differences between strains in a visually informative manner and present the users with publication-ready images.
The memory and time usage required by the PanACEA scripts to run does not exceed the capabilities of most laptops, as shown in Additional file 1: Table S1. We compared runs of pan-chromosomes generated from between 20 and 219 genomes. The compute times ranged from 80 to 456 s, while the memory usage varied from 208 Mb to 3.16 Gb. We further found that increasing the number of fGR paths also lead to an increase in these requirements - surprisingly somewhat independent of number of genomes. For instance, the 193 E. coli genome pan-chromosome has almost twice as many fGR paths compared to a 219 E. hormaechei genome pan-chromosome and showed relative increases in time and memory usage. However, this increase is limited to a few minutes in terms of the CPU and a few gigabytes in terms of memory usage.
The modularity of PanACEA also allows for more functionality to be added. Further possible functions that can be included in future versions of PanACEA may include: multiple region views where genomes can be compared across neighboring fG and Core regions; additional gene annotation on the core region images, such as three letter gene names; graphs and text demonstrating the prevalence of different gene order and gene prevalence in clusters of genomes with the available metadata; and finally, to write additional scripts to transform the output from other pan-genome tools such as Roary  so that it can be used as input for PanACEA.
PanACEA is an interactive visualization tool that leverages bacterial genomic data for the analysis of pan-genomes in the context of a consensus pan-chromosome. Its browser interface displays customizable annotation features such as the anti-microbial resistance and gene ontologies, which expedite the point-and-click exploration of pan-chromosomes when compared to text files and previous visualizations that lacked contextual browsing of variable regions. Its hierarchical design enables the navigation of both detailed and high level views of the data. The search and zoom functions permit users to identify genes and regions of interest and view these regions in the context of the full pan-chromosome, zoomed in close, or in the detail views in another window, as shown in our use case. PanACEA is database independent and browser agnostic, easy to install, and works off generalized flat files promoting interoperability across pan-genome software.
Availability and requirements
Project name: PanACEA.
Project home page: https://github.com/JCVenterInstitute/PanACEA/
Operating system(s): Platform independent.
Other requirements: PERL v5.22.1, BioPerl v1.007001.
License: GNU GPL.
Any restrictions to use by non-academics: none.
Antibiotic Resistance Ontology
flexible genomic island
flexible genome region
Resistance Gene Identifier
Tettelin H, Masignani V, Cieslewicz MJ, Donati C, Medini D, Ward NL, et al. Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial pan-genome. Proc Natl Acad Sci U S A. 2005;102:13950–5.
Remm M, Storm CEV, Sonnhammer ELL. Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001;314:1041–52.
Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of Ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.
Fouts DE, Brinkac L, Beck E, Inman J, Sutton G. PanOCT: automated clustering of orthologs using conserved gene neighborhood for pan-genomic analysis of bacterial strains and closely related species. Nucleic Acids Res. 2012;40:e172.
Ozer EA, Allen JP, Hauser AR. Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools spine and AGEnt. BMC Genomics. 2014;15:737.
Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, et al. Roary: Rapid large-scale prokaryote pan genome analysis. bioRxiv. 2015:019315.
Contreras-Moreira B, Vinuesa P. GET_HOMOLOGUES, a versatile software package for scalable and robust microbial Pangenome analysis. Appl Environ Microbiol. 2013;79:7696–701.
Blom J, Kreis J, Spänig S, Juhre T, Bertelli C, Ernst C, et al. EDGAR 2.0: an enhanced software platform for comparative gene content analyses. Nucleic Acids Res. 2016;44:W22–8.
Chaudhari NM, Gupta VK, Dutta C. BPGA- an ultra-fast pan-genome analysis pipeline. Sci Rep. 2016;6:srep24373.
Zhao Y, Wu J, Yang J, Sun S, Xiao J, Yu J. PGAP: pan-genomes analysis pipeline. Bioinformatics. 2012;28:416–8.
Chan AP, Sutton G, DePew J, Krishnakumar R, Choi Y, Huang X-Z, et al. A novel method of consensus pan-chromosome assembly and large-scale comparative analysis reveal the highly flexible pan-genome of Acinetobacter baumannii. Genome Biol. 2015;16:143.
Chavda KD, Chen L, Fouts DE, Sutton G, Brinkac L, Jenkins SG, et al. Comprehensive genome analysis of Carbapenemase-producing Enterobacter spp.: new insights into phylogeny, population structure, and resistance mechanisms. MBio. 2016;7:e02093–16.
Blom J, Albaum SP, Doppmeier D, Pühler A, Vorhölter F-J, Zakrzewski M, et al. EDGAR: a software framework for the comparative analysis of prokaryotic genomes. BMC Bioinformatics. 2009;10:154.
Laing C, Buchanan C, Taboada EN, Zhang Y, Kropinski A, Villegas A, et al. Pan-genome sequence analysis using Panseq: an online tool for the rapid analysis of core and accessory genomic regions. BMC Bioinformatics. 2010;11:461.
Brittnacher MJ, Fong C, Hayden HS, Jacobs MA, Radey M, Rohmer L. PGAT: a multistrain analysis resource for microbial genomes. Bioinformatics. 2011;27:2429–30.
Ding W, Baumdicker F, Neher RA pX. pan-genome analysis and exploration. bioRxiv. 2017:072082.
Pantoja Y, Pinheiro K, Veras A, Araújo F, de SAL, Guimarães LC, et al. PanWeb: a web interface for pan-genomic analysis. PLoS One. 2017;12:e0178154.
Hennig A, Bernhardt J, Nieselt K. Pan-Tetris: an interactive visualisation for Pan-genomes. BMC Bioinformatics. 2015;16:S3.
Pedersen TL, Nookaew I, Wayne Ussery D, Månsson M. PanViz: interactive visualization of the structure of functionally annotated pangenomes. Bioinformatics. 2017;33:1081–2.
Treangen TJ, Ondov BD, Koren S, Phillippy AM. The harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol. 2014;15:524.
Sheikhizadeh S, Schranz ME, Akdel M, de Ridder D, Smit S. PanTools: representation, storage and exploration of pan-genomic data. Bioinformatics. 2016;32:i487–93.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium Nat Genet. 2000;25:25–9.
Gene Ontology Consortium. Gene ontology consortium: going forward. Nucleic Acids Res. 2015;43:D1049–56.
Jia B, Raphenya AR, Alcock B, Waglechner N, Guo P, Tsang KK, et al. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017;45:D566–73.
Hoffmann H, Stindl S, Ludwig W, Stumpf A, Mehlen A, Monget D, et al. Enterobacter hormaechei subsp. oharae subsp. nov., E. hormaechei subsp. hormaechei comb. nov., and E. hormaechei subsp. steigerwaltii subsp. nov., three new subspecies of clinical importance. J Clin Microbiol. 2005;43:3297–303.
Reizer J, Ramseier TM, Reizer A, Charbit A, Saier MH. Novel phosphotransferase genes revealed by bacterial genome sequencing: a gene cluster encoding a putative N-acetylgalactosamine metabolic pathway in Escherichia coli. Microbiology. 1996;142:231–50.
The authors would like to thank Chris Greco, Pratap,Venepally, and Harinder Singh for his assistance in software testing and critical review of the help manual, and Matthew LaPointe for his suggestions on software and image generation.
This project has been funded in whole or part with federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services under Award Number U19AI110819.
Availability of data and materials
All software and example data sets are available on GitHub at https://github.com/JCVenterInstitute/PanACEA/
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Clarke, T.H., Brinkac, L.M., Inman, J.M. et al. PanACEA: a bioinformatics tool for the exploration and visualization of bacterial pan-chromosomes. BMC Bioinformatics 19, 246 (2018). https://doi.org/10.1186/s12859-018-2250-y