title: The large-scale blast score ratio (LS-BSR) pipeline: a method to rapidly compare genetic content between bacterial genomes.
abstract: Background. As whole genome sequence data from bacterial isolates becomes cheaper to generate, computational methods are needed to correlate sequence data with biological observations. Here we present the large-scale BLAST score ratio (LS-BSR) pipeline, which rapidly compares the genetic content of hundreds to thousands of bacterial genomes, and returns a matrix that describes the relatedness of all coding sequences (CDSs) in all genomes surveyed. This matrix can be easily parsed in order to identify genetic relationships between bacterial genomes. Although pipelines have been published that group peptides by sequence similarity, no other software performs the rapid, large-scale, full-genome comparative analyses carried out by LS-BSR. Results. To demonstrate the utility of the method, the LS-BSR pipeline was tested on 96 Escherichia coli and Shigella genomes; the pipeline ran in 163 min using 16 processors, which is a greater than 7-fold speedup compared to using a single processor. The BSR values for each CDS, which indicate a relative level of relatedness, were then mapped to each genome on an independent core genome single nucleotide polymorphism (SNP) based phylogeny. Comparisons were then used to identify clade specific CDS markers and validate the LS-BSR pipeline based on molecular markers that delineate between classical E. coli pathogenic variant (pathovar) designations. Scalability tests demonstrated that the LS-BSR pipeline can process 1,000 E. coli genomes in 27-57 h, depending upon the alignment method, using 16 processors. Conclusions. LS-BSR is an open-source, parallel implementation of the BSR algorithm, enabling rapid comparison of the genetic content of large numbers of genomes. The results of the pipeline can be used to identify specific markers between user-defined phylogenetic groups, and to identify the loss and/or acquisition of genetic information between bacterial isolates. Taxa-specific genetic markers can then be translated into clinical diagnostics, or can be used to identify broadly conserved putative therapeutic candidates.
authors.affiliation: Division of Pathogen Genomics, Translational Genomics Research Institute , Flagstaff, AZ , USA ; Department of Biological Sciences, Northern Arizona University , Flagstaff, AZ , USA ; Center for Microbial Genetics and Genomics, Northern Arizona University , Flagstaff, AZ , USA.
authors.forename: Jason W
authors.initials: JW
authors.lastname: Sahl
authors.affiliation: Department of Biological Sciences, Northern Arizona University , Flagstaff, AZ , USA.
authors.forename: J Gregory
authors.initials: JG
authors.lastname: Caporaso
authors.affiliation: Department of Microbiology and Immunology, Institute for Genome Sciences, University of Maryland School of Medicine , Baltimore, MD , USA.
authors.forename: David A
authors.initials: DA
authors.lastname: Rasko
authors.affiliation: Division of Pathogen Genomics, Translational Genomics Research Institute , Flagstaff, AZ , USA ; Department of Biological Sciences, Northern Arizona University , Flagstaff, AZ , USA ; Center for Microbial Genetics and Genomics, Northern Arizona University , Flagstaff, AZ , USA.
authors.forename: Paul
authors.initials: P
authors.lastname: Keim
keywords: Bioinformatics; Comparative genomics; Genomics; Microbiology; Pathogens
references.citation: Bioinformatics. 2012 Feb 1;28(3):416-8
references.pmid: 22130594
references.citation: PLoS One. 2013;8(1):e53786
references.pmid: 23342000
references.citation: Proc Natl Acad Sci U S A. 2013 Jul 30;110(31):12810-5
references.pmid: 23858472
references.citation: Appl Environ Microbiol. 2013 Dec;79(24):7696-701
references.pmid: 24096415
references.citation: J Bacteriol. 2008 Oct;190(20):6881-93
references.pmid: 18676672
references.citation: Proc Natl Acad Sci U S A. 2011 May 17;108(20):8329-34
references.pmid: 21536867
references.citation: Bioinformatics. 2010 Oct 1;26(19):2460-1
references.pmid: 20709691
references.citation: Genome Res. 2003 Sep;13(9):2178-89
references.pmid: 12952885
references.citation: Bioinformatics. 2011 Feb 1;27(3):334-42
references.pmid: 21148543
references.citation: Nucleic Acids Res. 2002 Apr 1;30(7):1575-84
references.pmid: 11917018
references.citation: BMC Bioinformatics. 2010 Mar 08;11:119
references.pmid: 20211023
references.citation: Infect Immun. 2011 Feb;79(2):950-60
references.pmid: 21078854
references.citation: Appl Environ Microbiol. 2009 Dec;75(23):7537-41
references.pmid: 19801464
references.citation: Nucleic Acids Res. 2012 Jan;40(Database issue):D48-53
references.pmid: 22144687
references.citation: BMC Genomics. 2014 Jan 03;15:8
references.pmid: 24387194
references.citation: BMC Bioinformatics. 2005 Jan 05;6:2
references.pmid: 15634352
references.citation: Nat Genet. 2011 May;43(5):491-8
references.pmid: 21478889
references.citation: Mol Syst Biol. 2008;4:198
references.pmid: 18523433
references.citation: BMC Genomics. 2008 Feb 08;9:75
references.pmid: 18261238
references.citation: Curr Opin Microbiol. 2008 Oct;11(5):472-7
references.pmid: 19086349
references.citation: Methods Enzymol. 2006;411:134-93
references.pmid: 16939790
references.citation: Genome Res. 2002 Apr;12(4):656-64
references.pmid: 11932250
references.citation: Nucleic Acids Res. 1997 Sep 1;25(17):3389-402
references.pmid: 9254694
references.citation: PLoS One. 2010 Mar 10;5(3):e9490
references.pmid: 20224823
references.citation: PLoS One. 2013;8(1):e54287
references.pmid: 23365658
references.citation: Bioinformatics. 2014 May 1;30(9):1297-9
references.pmid: 24420766