Genome sequence and description of the mosquitocidal and heavy metal tolerant strain Lysinibacillus sphaericus CBAM5
© Peña-Montenegro et al.; licensee BioMed Central. 2015
Received: 12 June 2014
Accepted: 21 November 2014
Published: 20 January 2015
Lysinibacillus sphaericus CBAM5, was isolated from subsurface soil of oil well explorations in the Easter Planes of Colombia. This strain has potential in bioremediation of heavy-metal polluted environments and biological control of Culex quinquefasciatus. According to the phylogenetic analysis of 16S rRNA gene sequences, the strain CBAM5 was assigned to the Lysinibacillus sphaericus taxonomic group 1 that comprises mosquito pathogenic strains. After a combination assembly-integration, alignment and gap-filling steps, we propose a 4,610,292 bp chromosomal scaffold. The whole genome (consisting of 5,146,656 bp long, 60 contigs and 5,209 predicted-coding sequences) revealed strong functional and syntenial similarities to the L. sphaericus C3-41 genome. Mosquitocidal (Mtx), binary (Bin) toxins, cereolysin O, and heavy metal resistance clusters from nik, ars, czc, mnt, ter, cop, cad, and znu operons were identified.
KeywordsLysinibacillus sphaericus CBAM5 DNA homology Binary toxins Mosquitocidal toxins S-layer proteins Heavy metal tolerance
Lysinibacillus sphaericus is one of the bacteria used as a bio-insecticide as part of vector control programs against tropical diseases, such as malaria, filariasis, yellow fever, dengue fever and West Nile virus . L. sphaericus isolates may be classified according to their larvicidal activity into high and low toxicity strains. High- and low-toxicity strains synthesize mosquitocidal toxins (Mtx) in vegetative growth cells . Highly toxic strains produce a binary toxin coded by binA and binB genes in sporulating stages . In addition, L. sphaericus larvicidal toxicity may be explained due to expression of Cry48/Cry49 toxin  and the S-layer protein . Vegetative and sporulated cells of L. sphaericus CBAM5 are pathogenic towards Culex quinquefasciatus larvae . LC50 (50% lethal concentration) toward C. quinquefasciatus larvae of strain CBAM5 is 1400 cells/mL from sporulated cultures, being this isolate assigned as a high-toxicity strain .
The biotechnological application of L. sphaericus is not limited to biological control. L. sphaericus biomass has been applied in the bioremediation of heavy metals, such as cobalt, copper, chromium and lead  with specific metal binding in the cell surface . Native Colombian isolates L. sphaericus OT4b.31 and IV(4)10 showed heavy metal biosorption in living and dead biomass, both strains expressing the S-layer proteins . L. sphaericus strain CBAM5, along with other 24 native isolates, shown effective growth in arsenate, hexavalent chromium and/or lead [6, 10].
Considering that Lysinibacillus sphaericus CBAM5 is a relevant native strain, not only by its highly toxic larvicidal activity but also by its heavy metal tolerance, we have chosen this strain to analyze its genomic sequence. In this report, we present a summary classification, and set of general features for Lysinibacillus sphaericus strain CBAM5 including previously unreported aspects of its phenotype, together with the description of its genome sequence and annotation.
Lysinibacillus sphaericus is an aerobic, mesophilic, spore-forming and Gram-positive bacterium, commonly isolated from soil and water . Formerly known as Bacillus sphaericus, the species was later reassigned to the genus Lysinibacillus because of its distinctive peptidoglycan membrane composition, and physiological features [12, 13]. Lysinibacillus sphaericus strains have been classified into five DNA homology groups, where mosquito larvicidal strains were placed into DNA subgroup IIA  while the subgroup IIB was reclassified as Lysinibacillus fusiformis . Later, according to 16S rDNA and lipid profile comparisons, Lysinibacillus sphaericus strains were classified into seven similarity subgroups, of which only four retained the previous description by Krych et al. . Groups VI and VII were later reclassified as new species . Because of the phenotypic and genetic diversity summarized above, most of the groups remain designated as Lysinibacillus sphaericus sensu lato.
Classification and general features of Lysinibacillus sphaericus CBAM5 according to the MIGS recommendations 
Evidence code a
Species Lysinibacillus sphaericus
Positive in vegetative cells, variable in sporulating stages
15 – 40°C
Growth in Luria-Bertani broth (5% NaCl)
Pathogenic toward Culex quinquefasciatus larvae
Eastern Planes oil basins, Colombia
Sample collection time
350 m above sea level
Genome sequencing information
Genome project history
Genome sequencing project information
Improved high-quality draft
One paired end tags 90:90 bp with 500 bp insert
Illumina Hi-Seq 2000
CISA version 1.3, SOAPdenovo version 2.04, Velvet version 1.2.10, ABySS version 1.3.7, CLC Assembly Cell version 4.0.10
Gene calling method
Glimmer3, tRNAscan-SE, RNAmmer
Genbank Date of Release
February 1, 2014
Biotechnology, metabolic pathway
Growth conditions and DNA isolation
Lysinibacillus sphaericus strain CBAM5 was grown in nutrient broth for 16 hours at 30°C and 150 rev/min. High molecular weight DNA was isolated using the EasyDNA® Kit (Carlsbad, CA, USA. Invitrogen) as indicated by the manufacturer. DNA purity and concentration were determined in a NanoDrop spectrophotometer (Wilmington, DE, USA. Thermo Scientific).
Genome sequencing and assembly
After DNA extraction, samples were sent to the Beijing Genome Institute (BGI) Americas Laboratory (Tai Po, Hong Kong). Purified DNA sample was first sheared into smaller fragments with a desired size by a Covaris E210 ultrasonicator. Then the overhangs resulting from fragmentation were converted into blunt ends by using T4 DNA polymerase, Klenow Fragment and T4 polynucleotide kinase. After adding an “A” base to the 3’ end of the blunt phosphorylated DNA fragments, adapters were ligated to the ends of the DNA fragments. The desired fragments were purified though gel-electrophoresis, then selectively enriched and amplified by PCR. The index tag was introduced into the adapter at the PCR stage as appropriate, and a library quality test was performed. Lastly, qualified, short, paired-ends of 90:90 bp length with 500 bp insert libraries were used to cluster preparation and to conduct whole-shotgun sequencing in Illumina Hi-Seq 2000 sequencers.
Using the FASTX-Toolkit version 0.6.1  and clean_reads version 0.2.3 from the ngs_backbone pipeline  reads were trimmed and quality filtered. Four preliminary assemblies were obtained by using: SOAPdenovo version 2.04 , Velvet version 1.2.10 , ABySS version 1.3.7 , and CLC Assembly Cell version 4.0.10 . Those assemblies were integrated in the CISA pipeline resulting in a consensus assembly . SOAPdenovo and CLC Assembly Cell packages included automatic scaffolding and k-mer/overlapping optimization steps. To obtain structural insight of a chromosomal scaffold, we used CONTIGuator.2 , using the Lysinibacillus sphaericus strain C3-41 chromosome (accession number: CP000817.1) as reference. Some gaps were successfully filled by using GapFiller . Gap-filling steps were applied over each one of the preliminary assemblies and over the final consensus assembly. Quality assessment of the assembly was performed with iCORN . The error rate of the final assembly is less than 1 in 1,000,000 bp. We compared the chromosomal assembly of L. sphaericus CBAM5 with the chromosome sequences of L. sphaericus C3-41 and L. sphaericus OT4b.31 by maximal unique matching of translated sequences with PROmer , and a read mapping single nucleotide polymorphism (SNP) effect analysis with SnpEff package .
The Glimmer 3 gene finder was used to identify and extract sequences for potential coding regions. To achieve the functional annotation steps, the RAST server  and Blast2GO pipelines  were used. Blast2GO performed the blasting, GO-mapping and annotation steps; which included a description according to the ProDom, FingerPRINTScan, PIR-PSD, Pfam, TIRGfam, PROSITE, ProDom, SMART, SuperFamily, Pattern, Gene3D, PANTHER, SignalIP and TM-HMM databases. The results were summarized with InterPro . Additionally, a GO-EnzymeCode mapping step was used to retrieve KEGG pathway-maps. tRNA genes were identified by using tRNAscan-SE  and rRNA genes by using RNAmmer . The possible orthologs of the genome were identified based on the COG database and classified accordingly . Prophage region prediction was also conducted by using the PHAST tool .
Summary of genome
Nucleotide content and gene count levels of the genome
% of total a
Genome size (bp)
DNA GC content (bp)
DNA coding region (bp)
Number of replicons
Genes in paralog clusters
Genes assigned to COGs
1 or more conserved domains
2 or more conserved domains
3 or more conserved domains
Genes with function prediction
Genes assigned Pfam domains
Genes with signal peptides
Genes with transmembrane helices
Number of genes associated with the 25 general COG functional categories
% age a
RNA processing and modification
Replication, recombination and repair
Chromatin structure and dynamics
Cell cycle control, mitosis and meiosis
Signal transduction mechanisms
Cell wall/membrane biogénesis
Intracellular trafficking and secretion
Posttranslational modification, protein turnover, chaperones
Energy production and conversión
Carbohydrate transport and metabolism
Amino acid transport and metabolism
Nucleotide transport and metabolism
Coenzyme transport and metabolism
Lipid transport and metabolism
Inorganic ion transport and metabolism
Secondary metabolites biosynthesis, transport and catabolism
General function prediction only
Not in COGs
Insights into the genome
The origin of replication of the chromosome of L. sphaericus CBAM5 was estimated by similarities to several features of the corresponding regions in L. sphaericus C3-41, Bacillus sp. B-14905 and other close related bacteria, including colocalization of the genes: dnaA, dnaN, dnaX, recR and recF; and GC nucleotide skew [(G–C)/(G + C)] analysis. In the contig 19 (EWH31640:EWH31645) we found a typical cluster consisting of dnaA, dnaN, recF, gyrA and gyrB boxes. The predicted genes dnaB, dnaI, dnaG, dnaE, holA, holB, priAB, polA and recA were also found spread in the chromosomal and extrachosmomal sequences. The replication termination site of the chromosomal scaffold is believed to be localized near 2.92 Mbp in the contig 14. According to GC skew analysis, the coding bias for the two strands of the chromosome is for the majority of CDSs to be on the outer strand from 0 to ~2.92 Mbp, and on the inner strand from ~2.92 Mbp to the end of the chromosomal scaffold (contig 19, Figure 2). This was also confirmed by the presence of parC (EWH32537) and parE (EWH32538), which encode the subunits of the chromosome-partitioning enzyme topoisomerase IV. Similar to previous reports [50, 51], we did not find the homolog of rtp (replication terminator protein-encoding gene) in the chromosomal assembly of CBAM5.
Lysinibacillus sphaericus CBAM5 displays 28 CDSs annotated as transposases, including three allocated in the extrachromosomal sequences. The most frequent families are IS1182, IS3 and IS4. In addition, four incomplete prophage regions were identified as follows: Thermus phage φOH2 (contig 12), Burkholderia phage ST79 (contig 14), and two regions comprising the Clostridium phage φSM101 (contigs 14 and 28). Prophage regions φOH2 and ST79 included putative encoding sequences for tail, lysis and baseplate proteins. None of the reported phages has been described in the Colombian strain L. sphaericus OT4b.31 .
The genome of L. sphaericus CBAM5 shows a wide repertoire of potential encoding sequences in terms of mosquitocidal toxins. In the contig 11, we found Mtx1 (EWH35097) and Mtx2 (EWH35034) CDSs located in an identical cluster as Hu et al.  described in the genome of L. sphaericus C3-41. This cluster includes two insertion sequences, one of them consisting of a disrupted transposase between the mtx1 and mtx2, as well. One Mtx3 CDSs (EWH32377) was found in contig 14. Upstream of this sequence, we could identify some IS3 family mobile elements and putative DeoR family transcriptional regulators. In addition, we found one hypothetical toxin from the family Mtx2 (PFam PF03318) in contig 11 (EWH35106) and a putative cereolysin O CDS (EWH31995) being described to be active against the German cockroach Blattela germanica  in contig 15.
Surface (S) layer proteins and toxic metal resistances
L. sphaericus CBAM5 exhibits 21 CDSs described as surface (S) layer proteins or S-layer homologs in its genome. The fragment covered from EWH35069 to EWH35072 includes four CDSs encoding for a variable protein, a putative S-layer associated protein, a P60 invasion-associated protein and a N-acetylmuramoyl-L-alanine amidase. Probably the genes located in this fragment may participate in the larvicidal activity of the strain CBAM5, given that the same genes have been described as differentially expressed in virulent infections of Lysteria monocytogenes . A total of 14 CDSs show three SLHs motifs near to the N terminal region, similarly to the slpC gene previously described in native strains . In addition, we found two S-layer surface array proteins in the chromosomal scaffold and another in extracromosomal sequences.
A total of 64 CDSs corresponded to encoding sequences involved in responses against toxic metal(oid)s. Among those coding sequences, we found the following operons: arsRBCDA, arsRBC, cadCA, mntABCD, nikABCDE, terD-terD-terD, zurR-znuBC, and czrA-czcD-csoR-copZA. We could identify various genes probably involved in metal(loid) resistances spread across the genome (Additional file 4: Table S2). The chrA gene seems to be the only representative of the chr operon in the genome of L. sphaericus CBAM5. Previous reports have shown that microorganisms bearing chrA homologues display highly variable resistance levels against Cr(VI) . However, two superoxide dismutase putative proteins (EWH33050, EWH30224) and several CDS ascribed as flavin reductases (EC 22.214.171.124), nitroreductases (EC 126.96.36.199) and quinone reductases (EC 188.8.131.52) could cooperate in the Cr(VI) resistance, in agreement with previous reports [55, 56].
Given the heavy metal resistance of L. sphaericus CBAM5 in polluted environments, and supported by the identification of genes in Additional file 4: Table S2, we could expect the assistance of efflux pump systems and heavy-metal resistance proteins specific to As, Sb, Ni, Zn, Cu, Cd, Te, Cr, Mn and Co. By the evaluation of coalescent models, Villegas-Torres et al.  proposed that L. sphaericus CBAM5 may have acquired the arsC gene through recent events of horizontal gene transfer as a possible adaptation to polluted environments. However, we found highly similar homologues of heavy metal resistance proteins of the CBAM5 strain in microorganisms isolated from non-polluted environments (i.e. czrA-czcD-csoR-copZA, cadCA, and arsRBC in L sphaericus OT4B.31 ). Further analysis on plasmids, prophage content, or conjugation factors may clarify the origin of resistance (as well as larvicidal) genes. Finally, based in the KEGG analysis, some predicted proteins might participate in peripheral pathways for the degradation of geraniol, chlorocyclohexane, chlorobenzene, benzoate, bisphenol, fluorobenzoate, toluene, chloroalkane, chloroalkene, naphthalene, aminobenzoate, styrene, atrazine, limonene and pinene.
Lysinibacillus sphaericus CBAM5 was isolated from drilling mineral base oil samples at the subsurface soil level. By comparing the chromosomal sequences between L. sphaericus strains CBAM5 and C3-41, we identified distinctive similarities of the DNA homology group IIA. The evidence of the binary toxins allocated in a conserved cluster delimited by mobile elements, resembles a probable phage invasion in the DNA subgroup IIA of the Lysinibacillus sphaericus species. Along with the biological control potential given by the Mtx, Bin and cerolysin toxins, L. sphaericus CBAM5 displays encoding sequences for S-layer proteins and heavy-metal efflux pumps, which may confer resistance to As, Sb, Ni, Zn, Cu, Cd, Te, Cr, Mn and Co in polluted environments.
This work was performed under the auspices of the Grant 1204-452-21129 from the Instituto Colombiano para el fomento de la Investigación Francisco José de Caldas – Colciencias, the Research Fund from the Faculty of Sciences at Universidad de los Andes, and the Centro de Investigaciones Microbiológicas (CIMIC).
- Berry C: The bacterium, Lysinibacillus sphaericus , as an insect pathogen. J Invertebr Pathol 2012, 109:1–10. 10.1016/j.jip.2011.11.008View ArticlePubMedGoogle Scholar
- Thanabalu T, Porter AG: Efficient expression of a 100-kilodalton mosquitocidal toxin in protease-deficient recombinant Bacillus sphaericus . Appl Environ Microbiol 1995, 61:4031–6.PubMed CentralPubMedGoogle Scholar
- Baumann P, Clark MA, Baumann L, Broadwell AH: Bacillus sphaericus as a mosquito pathogen: properties of the organism and its toxins. Microbiol Rev 1991, 55:425–36.PubMed CentralPubMedGoogle Scholar
- Jones GW, Nielsen-LeRoux C, Yang Y, Yuan Z, Dumas VF, Monnerat RG, Berry C: A new Cry toxin with a unique two-component dependency from Bacillus sphaericus . FASEB J 2007, 21:4112–20. 10.1096/fj.07-8913comView ArticlePubMedGoogle Scholar
- Lozano LC, Ayala JA, Dussán J: Lysinibacillus sphaericus S-layer protein toxicity against Culex quinquefasciatus . Biotechnol Lett 2011, 33:2037–41. 10.1007/s10529-011-0666-9View ArticlePubMedGoogle Scholar
- Lozano LC, Dussán J: Metal tolerance and larvicidal activity of Lysinibacillus sphaericus . World J Microbiol Biotechnol 2013, 29:1383–9. 10.1007/s11274-013-1301-9View ArticlePubMedGoogle Scholar
- Tuzen M, Uluozlu OD, Usta C, Soylak M: Biosorption of copper(II), lead(II), iron(III) and cobalt(II) on Bacillus sphaericus -loaded Diaion SP-850 resin. Anal Chim Acta 2007, 581:241–6. 10.1016/j.aca.2006.08.040View ArticlePubMedGoogle Scholar
- Pollmann K, Raff J, Merroun M, Fahmy K, Selenska-Pobell S: Metal binding by bacteria from uranium mining waste piles and its technological applications. Biotechnol Adv 2006, 24:58–68. 10.1016/j.biotechadv.2005.06.002View ArticlePubMedGoogle Scholar
- Velásquez L, Dussan J: Biosorption and bioaccumulation of heavy metals on dead and living biomass of Bacillus sphaericus . J Hazard Mater 2009, 167:713–6. 10.1016/j.jhazmat.2009.01.044View ArticlePubMedGoogle Scholar
- Villegas-Torres MF, Bedoya-Reina OC, Salazar C, Vives-Florez MJ, Dussan J: Horizontal arsC gene transfer among microorganisms isolated from arsenic polluted soil. Int Biodeterior Biodegradation 2011, 65:147–52. 10.1016/j.ibiod.2010.10.007View ArticleGoogle Scholar
- Claus D, Berkeley RCW: Genus Bacillus Cohn 1872, 174AL. In Bergey’s Man Syst Bacteriol. Volume 2. Edited by: Sneath PHA, Mair NS, Sharpe ME, Holt JG. Baltimore: The Williams and Wilkins Co; 1986:1105–39.Google Scholar
- White PJ, Lotay HK: Minimal nutritional requirements of Bacillus sphaericus NCTC9602 and 26 other strains of this species: the majority grow and sporulate with acetate as sole major source of carbon. J Gen Microbiol 1980, 118:13–9.Google Scholar
- Ahmed I, Yokota A, Yamazoe A, Fujiwara T: Proposal of Lysinibacillus boronitolerans gen. nov. sp. nov., and transfer of Bacillus fusiformis to Lysinibacillus fusiformis comb. nov. and Bacillus sphaericus to Lysinibacillus sphaericus comb. nov. Int J Syst Evol Microbiol 2007, 57:1117–25. 10.1099/ijs.0.63867-0View ArticlePubMedGoogle Scholar
- Krych VK, Johnson JL, Yousten AA: Deoxyribonucleic acid homologies among strains of Bacillus sphaericus . Int J Syst Bacteriol 1980, 30:476–84. 10.1099/00207713-30-2-476View ArticleGoogle Scholar
- Nakamura LK: Phylogeny of Bacillus sphaericus -like organisms. Int J Syst Evol Microbiol 2000, 50:1715–22.PubMedGoogle Scholar
- Nakamura LK, Shida O, Takagi H, Komagata K: Bacillus pycnus sp. nov. and Bacillus neidei sp. nov., round-spored bacteria from soil. Int J Syst Evol Microbiol 2002,52(Pt 2):501–5.PubMedGoogle Scholar
- Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 1993, 10:512–26.PubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S: MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 2013, 30:2725–9. 10.1093/molbev/mst197View ArticlePubMed CentralPubMedGoogle Scholar
- Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, Tatusova T, Thomson N, Allen MJ, Angiuoli SV, Ashburner M, Axelrod N, Baldauf S, Ballard S, Boore J, Cochrane G, Cole J, Dawyndt P, De Vos P, DePamphilis C, Edwards R, Faruque N, Feldman R, Gilbert J, Gilna P, Glöckner FO, Goldstein P, Guralnick R, Haft D, Hancock D, et al.: The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol 2008, 26:541–7. 10.1038/nbt1360View ArticlePubMed CentralPubMedGoogle Scholar
- Woese CR, Kandler O, Wheelis ML: Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci 1990, 87:4576–9. 10.1073/pnas.87.12.4576View ArticlePubMed CentralPubMedGoogle Scholar
- Gibbons NE, Murray RGE: Proposals Concerning the Higher Taxa of Bacteria. Int J Syst Bacteriol 1978, 28:1–6. 10.1099/00207713-28-1-1View ArticleGoogle Scholar
- Garrity GM, Holt JG: The Road Map to the Manual. In Bergey’s Man Syst Bacteriol. Volume 1. Second edition. Edited by: Garrity GM, Boone DR, Castenholz RW. New York: Springer; 2001:119–69.View ArticleGoogle Scholar
- Murray RGE: The Higher Taxa, or, a Place for Everything…? In Bergey’s Man Syst Bacteriol. Volume 1. First edition. Edited by: Holt JG. Baltimore: The Williams and Wilkins Co; 1984:31–4.Google Scholar
- Ludwig W, Schleifer KH, Whitman WB, Krieg NR, Ludwig W, Rainey EA, Schleifer KH, Withman WB: Class I. Bacilli Class Nov . In Bergey’s Manual of Systematic Bacteriology. In: Vol 3 (The Firmicutes). Second. Volume 3. Edited by: De Vos P, Garrity GM, Jones D. Dordrecht, Heidelberg, London, New York: Springer; 2009:19–20.Google Scholar
- List of new names and new combinations previously effectively, but not validly, published Int J Syst Evol Microbiol 2010, 60:469–72.
- Skerman VBD, McGowan V, Sneath PHA: Approved Lists of Bacterial Names. Int J Syst Bacteriol 1980, 30:225–420. 10.1099/00207713-30-1-225View ArticleGoogle Scholar
- Prévot AR: Hauduroy P, Ehringer G, Guillot G, Magrou J, Prévot AR, Rosset, Urbain A (eds) Dictionnaire Des Bactéries Pathogènes. 2nd edition. Paris: Masson; 1953:1–692.Google Scholar
- Fischer A: Untersuchungen über bakterien. Jahrbücher für Wissenschaftliche Bot 1895, 27:1–163.Google Scholar
- Jung MY, Kim J-S, Paek WK, Styrak I, Park I-S, Sin Y, Paek J, Park KA, Kim H, Kim HL, Chang Y-H: Description of Lysinibacillus sinduriensis sp. nov., and transfer of Bacillus massiliensis and Bacillus odysseyi to the genus Lysinibacillus as Lysinibacillus massiliensis comb. nov. and Lysinibacillus odysseyi co. Int J Syst Evol Microbiol 2012, 62:2347–55. 10.1099/ijs.0.033837-0View ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology. Nat Genet 2000, 25:25–9. 10.1038/75556View ArticlePubMed CentralPubMedGoogle Scholar
- FASTX-Toolkit. http://hannonlab.cshl.edu/fastx_toolkit/
- Blanca JM, Pascual L, Ziarsolo P, Nuez F, Cañizares J: ngs_backbone: a pipeline for read cleaning, mapping and SNP calling using Next Generation Sequence. BMC Genomics 2011., 12: Google Scholar
- Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, Li Y, Li S, Shan G, Kristiansen K, Li S, Yang H, Wang J, Wang J: De novo assembly of human genomes with massively parallel short read sequencing. Genome Res 2010, 20:265–72. 10.1101/gr.097261.109View ArticlePubMed CentralPubMedGoogle Scholar
- Zerbino DR, Birney E: Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 2008, 18:821–9. 10.1101/gr.074492.107View ArticlePubMed CentralPubMedGoogle Scholar
- Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJM, Birol I: ABySS: A parallel assembler for short read sequence data. Genome Res 2009, 19:1117–23. 10.1101/gr.089532.108View ArticlePubMed CentralPubMedGoogle Scholar
- CLC Genomics Workbench. http://www.clcbio.com/products/clc-genomics-workbench/
- Lin S-H, Liao Y-C: CISA: contig integrator for sequence assembly of bacterial genomes. PLoS One 2013, 8:e60843. 10.1371/journal.pone.0060843View ArticlePubMed CentralPubMedGoogle Scholar
- Galardini M, Biondi EG, Bazzicalupo M, Mengoni A: CONTIGuator: a bacterial genomes finishing tool for structural insights on draft genomes. Source Code Biol Med 2011., 6: Google Scholar
- Nadalin F, Vezzi F, Policriti A: GapFiller: a de novo assembly approach to fill the gap within paired reads. BMC Bioinformatics 2012,13(Suppl 14):S8. 10.1186/1471-2105-13-S14-S8View ArticlePubMed CentralPubMedGoogle Scholar
- Otto TD, Sanders M, Berriman M, Newbold C: Iterative Correction of Reference Nucleotides (iCORN) using second generation sequencing technology. Bioinformatics 2010, 26:1704–7. 10.1093/bioinformatics/btq269View ArticlePubMed CentralPubMedGoogle Scholar
- Delcher AL, Phillippy A, Carlton J, Salzberg SL: Fast algorithms for large-scale genome alignment and comparison. Nucleic Acids Res 2002, 30:2478–83. 10.1093/nar/30.11.2478View ArticlePubMed CentralPubMedGoogle Scholar
- Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM: A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 2012, 6:80–92. 10.4161/fly.19695View ArticlePubMed CentralPubMedGoogle Scholar
- Aziz R, Bartels D, Best A, DeJongh M, Disz T, Edwards R, Formsma K, Gerdes S, Glass E, Kubal M, Meyer F, Olsen G, Olson R, Osterman A, Overbeek R, McNeil L, Paarmann D, Paczian T, Parrello B, Pusch G, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST Server: Rapid Annotations using Subsystems Technology. BMC Genomics 2008., 9: Google Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res 2008, 36:3420–35. 10.1093/nar/gkn176View ArticlePubMed CentralPubMedGoogle Scholar
- Apweiler R, Attwood TK, Bairoch A, Bateman A, Birney E, Biswas M, Bucher P, Cerutti L, Corpet F, Croning MDR, Durbin R, Falquet L, Fleischmann W, Gouzy J, Hermjakob H, Hulo N, Jonassen I, Kahn D, Kanapin A, Karavidopoulou Y, Lopez R, Marx B, Mulder NJ, Oinn TM, Pagni M, Servant F, Sigrist CJA, Zdobnov EM: The InterPro database, an integrated documentation resource for protein families, domains and functional sites. Nucleic Acids Res 2001, 29:37–40. 10.1093/nar/29.1.37View ArticlePubMed CentralPubMedGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 1997, 25:955–64. 10.1093/nar/25.5.0955View ArticlePubMed CentralPubMedGoogle Scholar
- Lagesen K, Hallin P, Rødland EA, Staerfeldt H-H, Rognes T, Ussery DW: RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res 2007, 35:3100–8. 10.1093/nar/gkm160View ArticlePubMed CentralPubMedGoogle Scholar
- Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, Kiryutin B, Galperin MY, Fedorova ND, Koonin EV: The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res 2001, 29:22–8. 10.1093/nar/29.1.22View ArticlePubMed CentralPubMedGoogle Scholar
- Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS: PHAST: A Fast Phage Search Tool. Nucleic Acids Res 2011.Google Scholar
- Peña-Montenegro TD, Dussán J: Genome sequence and description of the heavy metal tolerant bacterium Lysinibacillus sphaericus strain OT4b.31. Stand Genomic Sci 2013, 9:42–56. 10.4056/sigs.4227894View ArticlePubMed CentralPubMedGoogle Scholar
- Hu X, Fan W, Han B, Liu H, Zheng D, Li Q, Dong W, Yan J, Gao M, Berry C, Yuan Z: Complete genome sequence of the mosquitocidal bacterium Bacillus sphaericus C3–41 and comparison with those of closely related Bacillus species. J Bacteriol 2008, 190:2892–902. 10.1128/JB.01652-07View ArticlePubMed CentralPubMedGoogle Scholar
- Nishiwaki H, Nakashima K, Ishida C, Kawamura T, Matsuda K: Cloning, functional characterization, and mode of action of a novel insecticidal pore-forming toxin, sphaericolysin, produced by Bacillus sphaericus . Appl Environ Microbiol 2007, 73:3404–11. 10.1128/AEM.00021-07View ArticlePubMed CentralPubMedGoogle Scholar
- Okinaka RT, Cloud K, Hampton O, Hoffmaster AR, Hill KK, Keim P, Koehler TM, Lamke G, Kumano S, Mahillon J, Manter D, Martinez Y, Ricke D, Svensson R, Jackson PJ: Sequence and organization of pXO1, the large Bacillus anthracis plasmid harboring the anthrax toxin genes. J Bacteriol 1999, 181:6509–15.PubMed CentralPubMedGoogle Scholar
- Camejo A, Buchrieser C, Couvé E, Carvalho F, Reis O, Ferreira P, Sousa S, Cossart P, Cabanes D: In vivo transcriptional profiling of Listeria monocytogenes and mutagenesis identify new virulence factors involved in infection. PLoS Pathog 2009, 5:e1000449. 10.1371/journal.ppat.1000449View ArticlePubMed CentralPubMedGoogle Scholar
- Viti C, Marchi E, Decorosi F, Giovannetti L: Molecular mechanisms of Cr(VI) resistance in bacteria and fungi. FEMS Microbiol Rev 2014, 38:633–59. 10.1111/1574-6976.12051View ArticlePubMedGoogle Scholar
- He M, Li X, Liu H, Miller SJ, Wang G, Rensing C: Characterization and genomic analysis of a highly chromate resistant and reducing bacterial strain Lysinibacillus fusiformis ZC1. J Hazard Mater 2011, 185:682–8. 10.1016/j.jhazmat.2010.09.072View ArticlePubMedGoogle Scholar
- Liolios K, Chen I-MA, Mavromatis K, Tavernarakis N, Hugenholtz P, Markowitz VM, Kyrpides NC: The Genomes On Line Database (GOLD) in 2009: status of genomic and metagenomic projects and their associated metadata. Nucleic Acids Res 2010,38(suppl 1):D346–54.View ArticlePubMed CentralPubMedGoogle 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.