Skip to main content

CoreGenes3.5: a webserver for the determination of core genes from sets of viral and small bacterial genomes

Abstract

Background

CoreGenes3.5 is a webserver that determines sets of core genes from viral and small bacterial genomes as an automated batch process. Previous versions of CoreGenes have been used to classify bacteriophage genomes and mine data from pathogen genomes.

Findings

CoreGenes3.5 accepts as input GenBank accession numbers of genomes and performs iterative BLASTP analyses to output a set of core genes. After completion of the program run, the results can be either displayed in a new window for one pair of reference and query genomes or emailed to the user for multiple pairs of small genomes in tabular format.

Conclusions

With the number of genomes sequenced increasing daily and interest in determining phylogenetic relationships, CoreGenes3.5 provides a user-friendly web interface for wet-bench biologists to process multiple small genomes for core gene determinations. CoreGenes3.5 is available at http://binf.gmu.edu:8080/CoreGenes3.5.

Findings

Background

Genes that are common between a set of genomes are known as core genes. Core sets of genes have been used to understand better bacterial genome evolution [1], orthology in viral genomes [2], viral evolutionary complexity [3], and to mine pathogen genomes [4]. Core genes have also been used to investigate the origins of photosynthesis [5], as well as to classify and untangle the taxonomy of bacteriophages [68]. With such a myriad of uses for core genes and the growing numbers of whole genome sequences, it is important to provide user-friendly and validated software tools for the determination of these genes from sets of genomes. Originally developed in 2002 [9], CoreGenes, a tool for the identification of shared and unique genes among (small) genomes, has been continually updated and refined in response to user demands [10]. These changes include increased robustness of the tool, as well as the ability to upload custom and proprietary data not deposited in GenBank. The major update to this version is the ability and versatility to batch process multiple pairs of small genomes, freeing the user from repetitive and time-consuming manual entry of genome sets. This is of benefit to users who have several large sets of genomes that they wish to analyze, for example a family of bacteriophages.

Other software tools have been developed for the determination of core genes including mGenomeSubtractor [11], CEGMA [12], nWayComp [13], and GenomeBlast [14]. mGenomeSubtractor and GenomeBlast both use BLAST-based algorithms to identify core genes. Of these, mGenomeSubtractor is primarily intended for use with bacterial genomes, whilst CEGMA is intended primarily for eukaryotic genomes; nWayComp and GenomeBlast are no longer accessible online, as is another genome comparison tool called GOAT [15]. In contrast, CoreGenes has been continuously available online since 2002, and shown to be invaluable in characterizing and re-determining the taxonomy and relationships of bacteriophages based on coding sequences [6, 7, 1619]. It is anticipated that this timely update of CoreGenes will enable the analysis of shared proteins among viral and small bacterial genomes in a faster and more efficient manner.

Implementation

CoreGenes3.5 is implemented using Java, Javascript, and HTML. It uses an iterative BLASTP algorithm that processes a reference genome and multiple query genomes. This is based on the GeneOrder algorithm described previously [20, 21]. Briefly, these genomes are retrieved directly from GenBank, or custom-entered by the user, and the gene translations are parsed from the files. Each protein from the first query genome is analyzed against the reference genome proteins using BLASTP, creating a new reference genome which is a subset of the original reference genome and which contains those proteins that meet or exceed the user-defined similarity threshold (BLASTP score). The second query genome is then BLASTP-analyzed against this new reference genome, creating another new reference genome. This iterative process continues until no further query genomes remain. The final output consists of proteins that are common to all the input genomes.

Batch processing of genomes in CoreGenes3.5 is implemented using Javascript and provides the facility to perform comparisons between many pairs of reference and query genomes. Given two lists of genome accession numbers, the script iteratively submits genomes in pairs to the main CoreGenes program such that each reference genome is individually compared with each query genome. When two identical input lists of genome accession numbers are provided, the script results in an ‘all versus all’ comparison, where all genomes are individually compared against each other as both reference and query. Additionally, each genome is compared against itself, providing an internal control where the number of paired proteins should equal the number of protein encoding genes denoted within the genome annotation.

As the BLASTP comparisons are performed ab-initio and not pre-computed, CoreGenes3.5 is limited to genome sizes of 2 Mb or less. While CoreGenes3.5 can take larger genomes as input, the time taken to process them also increases. Therefore, it is recommended that users submit genomes with the aforementioned limit.

Results and discussion

The input to CoreGenes3.5 consists of a reference genome and a query genome. These are entered into the text fields on the web interface as GenBank accession numbers (Figure 1). As noted, the BLASTP threshold score can be changed from its default of “75”. An email address must be entered to which results are sent. For genomes that are not available in GenBank, a link is provided to a custom genome entry page where the user can upload proprietary and unpublished data (these data and results are scrubbed from the server nightly to provide a level of confidentiality). In the batch entry mode (Figure 2), the input consists of two “comma-delineated” lists of GenBank accession numbers in addition to the threshold field. Additionally, CoreGenes3.5 works best with the Mozilla Firefox web browser.

Figure 1
figure 1

Web interface for CoreGenes3.5. A reference genome accession number is entered into the first text field and query genome accession numbers are entered into the subsequent fields. Options for modifying the BLASTP threshold score and an email address are provided.

Figure 2
figure 2

Web interface for batch processing in CoreGenes3.5. Reference and query accession numbers are entered in ‘comma delimited’ format. An email address is required so that output tables are emailed to the user following completion of the analysis.

The output of CoreGenes3.5 is a list of core genes in tabular format (Figure 3). Each gene is linked to its corresponding entry in GenBank. This allows users to identify particular genes of interest for further investigation, provide valuable insights for annotation of function and inform the design of wet-bench studies. In batch processing mode, results tables are emailed to the user in HTML format which can be viewed using any web browser.

Figure 3
figure 3

Partial output of the analysis between the Enterobacteria phages lambda and P22, generated using the batch mode input. Links (blue) are provided, pointing to whole genome data and specific proteins archived in GenBank. Statistics of how many homologs are found by CoreGenes are summarized at the top of the table.

The advent and continued development of next generation technologies has substantially increased the throughput and fidelity of genome sequence data. With reducing costs, the number of viral and bacterial genomes deposited in the International Nucleotide Sequence Databases/GenBank has grown rapidly (and continues to do so). It is therefore crucial to continue the development and improvement of novel and existing software tools that can efficiently mine this expanding wealth of sequence data and facilitate comparisons of multiple closely or distantly related genomes.

CoreGenes3.5 is the latest and most versatile update to a user-friendly tool for locating and identifying core genes from viral and small bacterial genomes. Like previous versions of CoreGenes, this newest version will be continually updated in response to demands from the user community. The ability of CoreGenes to deal with larger bacterial genomes is actively being addressed.

Conclusions

The batch processing feature of CoreGenes3.5 enables researchers to analyze multiple small genomes expeditiously using a web interface. This allows users to data mine the increasing numbers of genomes in sequence databases and to determine quickly the phylogenetic relationships amongst them.

Availability and requirements

Project name: CoreGenes3.5

Project home page: http://binf.gmu.edu:8080/CoreGenes3.5

Operating system(s): Platform independent

Programming language: Java

Any restrictions to use by non-academics: License required for commercial usage

References

  1. Liang W, Zhao Y, Chen C, Cui X, Yu J, Xiao J, Kan B: Pan-Genomic analysis provides insights into the genomic variation and evolution of Salmonella Paratyphi A. PLoS One. 2012, 7: e45346-

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  2. Garavaglia MJ, Miele SAB, Iserte JA, Belaich MN, Ghiringhelli PD: The ac53, ac78, ac101 and ac103 are newly discovered core genes in the family Baculoviridae. J Virol. 2012, 86: 12069-12079.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  3. Yutin N, Koonin EV: Hidden evolutionary complexity of Nucleo-Cytoplasmic Large DNA viruses of eukaryotes. Virol J. 2012, 9: 161-

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Mahadevan P, King JF, Seto D: Data mining pathogen genomes using GeneOrder and CoreGenes and CGUG: gene order, synteny and in silico proteomes. Int J Comput Biol Drug Des. 2009, 2: 100-114.

    Article  PubMed  CAS  Google Scholar 

  5. Mulkidjanian AY, Koonin EV, Makarova KS, Mekhedov SL, Sorokin A, Wolf YI, Dufresne A, Partensky F, Burd H, Kaznadzey D, Haselkorn R, Galperin MY: The cyanobacterial genome core and the origin of photosynthesis. Proc Natl Acad Sci USA. 2006, 103: 13126-13131.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  6. Lavigne R, Darius P, Summer EJ, Seto D, Mahadevan P, Nilsson AS, Ackermann HW, Kropinski AM: Classification of Myoviridae bacteriophages using protein sequence similarity. BMC Microbiol. 2009, 9: 224-

    Article  PubMed  PubMed Central  Google Scholar 

  7. Lavigne R, Seto D, Mahadevan P, Ackermann H-W, Kropinski AM: Unifying classical and molecular taxonomic classification: analysis of the Podoviridae using BLASTP-based tools. Res Microbiol. 2008, 159: 406-414.

    Article  PubMed  CAS  Google Scholar 

  8. Mahadevan P, Seto D: Taxonomic parsing of bacteriophages using core genes and in silico proteome-based CGUG and applications to small bacterial genomes. Adv Exp Med Biol. 2010, 680: 379-385.

    Article  PubMed  Google Scholar 

  9. Zafar N, Mazumder R, Seto D: CoreGenes: a computational tool for identifying and cataloging “core” genes in a set of small genomes. BMC Bioinforma. 2002, 3: 12-

    Article  Google Scholar 

  10. Mahadevan P, King JF, Seto D: CGUG: in silico proteome and genome parsing tool for the determination of “core” and unique genes in the analysis of genomes up to ca. 1.9 Mb. BMC Res Notes. 2009, 2: 168-

    Article  PubMed  PubMed Central  Google Scholar 

  11. Shao Y, He X, Harrison EM, Tai C, Ou H-Y, Rajakumar K, Deng Z: mGenomeSubtractor: a web-based tool for parallel in silico subtractive hybridization analysis of multiple bacterial genomes. Nucleic Acids Res. 2010, 38: W194-200.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  12. Parra G, Bradnam K, Korf I: CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007, 23: 1061-1067.

    Article  PubMed  CAS  Google Scholar 

  13. Yao J, Lin H, Doddapaneni H, Civerolo EL: nWayComp: a genome-wide sequence comparison tool for multiple strains/species of phylogenetically related microorganisms. In Silico Biol (Gedrukt). 2007, 7: 195-200.

    CAS  Google Scholar 

  14. Lu G, Jiang L, Helikar RMK, Rowley TW, Zhang L, Chen X, Moriyama EN: GenomeBlast: a web tool for small genome comparison. BMC Bioinforma. 2006, 7 (Suppl 4): S18-

    Article  Google Scholar 

  15. Kaluszka A, Gibas C: Interactive gene-order comparison for multiple small genomes. Bioinformatics. 2004, 20: 3662-3664.

    Article  PubMed  CAS  Google Scholar 

  16. Chibeu A, Lingohr EJ, Masson L, Manges A, Harel J, Ackermann H-W, Kropinski AM, Boerlin P: Bacteriophages with the ability to degrade uropathogenic Escherichia coli biofilms. Viruses. 2012, 4: 471-487.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Kropinski AM, Van den Bossche A, Lavigne R, Noben J-P, Babinger P, Schmitt R: Genome and proteome analysis of 7-7-1, a flagellotropic phage infecting Agrobacterium sp H13-3. Virol J. 2012, 9: 102-

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Lehman SM, Kropinski AM, Castle AJ, Svircev AM: Complete genome of the broad-host-range Erwinia amylovora phage phiEa21-4 and its relationship to Salmonella phage felix O1. Appl Environ Microbiol. 2009, 75: 2139-2147.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Villegas A, She Y-M, Kropinski AM, Lingohr EJ, Mazzocco A, Ojha S, Waddell TE, Ackermann H-W, Moyles DM, Ahmed R, Johnson RP: The genome and proteome of a virulent Escherichia coli O157:H7 bacteriophage closely resembling Salmonella phage Felix O1. Virol J. 2009, 6: 41-

    Article  PubMed  PubMed Central  Google Scholar 

  20. Celamkoti S, Kundeti S, Purkayastha A, Mazumder R, Buck C, Seto D: GeneOrder3.0: software for comparing the order of genes in pairs of small bacterial genomes. BMC Bioinforma. 2004, 5: 52-

    Article  Google Scholar 

  21. Mahadevan P, Seto D: Rapid pair-wise synteny analysis of large bacterial genomes using web-based GeneOrder4.0. BMC Res Notes. 2010, 3: 41-

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Chris Ryan for maintaining the server on which CoreGenes3.5 is hosted and Jason Seto for critical comments and software validation. We also thank Andrew Kropinski for suggestions and comments over the years to improve these software tools. Publication of this article was funded in part by the George Mason University Libraries Open Access Publishing Fund.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Padmanabhan Mahadevan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DT and PM implemented the batch processing in CoreGenes3.5. DT, DR, DS, and PM wrote the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Turner, D., Reynolds, D., Seto, D. et al. CoreGenes3.5: a webserver for the determination of core genes from sets of viral and small bacterial genomes. BMC Res Notes 6, 140 (2013). https://doi.org/10.1186/1756-0500-6-140

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1756-0500-6-140

Keywords