- Short genome report
- Open Access
Near complete genome sequence of the animal feed probiotic, Bacillus amyloliquefaciens H57
Standards in Genomic Sciencesvolume 11, Article number: 60 (2016)
Bacillus amyloliquefaciens H57 is a bacterium isolated from lucerne for its ability to prevent feed spoilage. Further interest developed when ruminants fed with H57-inoculated hay showed increased weight gain and nitrogen retention relative to controls, suggesting a probiotic effect. The near complete genome of H57 is ~3.96 Mb comprising 16 contigs. Within the genome there are 3,836 protein coding genes, an estimated sixteen rRNA genes and 69 tRNA genes. H57 has the potential to synthesise four different lipopeptides and four polyketide compounds, which are known antimicrobials. This antimicrobial capacity may facilitate the observed probiotic effect.
Bacillus amyloliquefaciens species have been taxonomically classified as part of the Bacillus subtilis group. Members of this group share substantial morphological similarities and near identical (98.1 %–99.8 %) 16S rRNA gene sequences . Other members of the Bacillus subtilis group include B. subtilis , B. atrophaeus , B. licheniformis , B. sonorensis , B. tequilensis , B. vallismortis , and the B. mojavensis subgroup. The production of bioactive metabolites, the ability to form spores and a lack of pathogenicity make members of the Bacillus subtilis group ideal candidates for use as probiotics. Strains of B. amyloliquefaciens synthesise non-ribosomal bioactive lipopeptides such as surfactin, fengycin, bacillomycin D and members of the iturin family [2–4]. These lipopeptides have demonstrated activity as antimicrobials and inhibit a wide range of bacterial and fungal pathogens [3, 5].
The strain B. amyloliquefaciens H57 (H57 hereafter) was first isolated in the search for a biological control agent to prevent fungal spoilage of hay . Due to its spore forming ability and production of antimicrobial compounds, H57 was revealed as the best candidate of a panel of isolates for commercialisation as a spoilage control agent under the product name HayRite™. Importantly, sheep and cattle fed on HayRite™ treated feed showed an increase in digestibility and nitrogen retention leading to increased live weight gain . This new development into the potential of H57 to act as a probiotic has led to further investigation of this strain.
Here, we present a summary description of the classification and features of H57, along with a sequencing description and annotation summary. The availability of a genome sequence for H57 will facilitate research into the probiotic effects observed in animals treated with this bacterium.
Classification and features
A near-complete 16S rRNA gene was identified in the H57 genome, which by BLAST  is most closely related (99 % identical) to other B. amyloliquefaciens strains including FZB42 ( B. amyloliquefaciens subsp. plantarum ; acc. NR075005.1), HPCAQB14 (acc. KF861603.1) and SB 3200 (acc. GU191911.1). Comparison of the average read coverage of the genome and 16S rRNA gene, suggests that H57 has 13 copies of the rRNA operon. A concatenated alignment of 99 single copy marker genes obtained from publicly available Bacillus genomes using HMMER  confirmed the classification of strain H57 as a member of the species B. amyloliquefaciens (Fig. 1).
H57 is a Gram-positive rod shaped bacterium averaging 2.5 μm in length and 1 μm in width (Fig. 2d). It is an aerobic spore forming bacterium that is motile with peritrichous flagella. H57 spores are centrally located and average 1.25 μm in length (Fig. 2b). Optimum growth occurs at a temperature of 29 °C and pH 7.0 (Table 1). The colony morphology of strain H57 is circular convex with undulate margins. When grown on a nutrient agar plate, colonies are an off-white colour as shown in Fig. 2c.
Genome sequencing information
Genome project history
Strain H57 was selected for sequencing due to its ability to act as a probiotic in agricultural animals. The draft genome was deposited in GenBank under the accession number LMUC00000000. Genome sequencing and assembly was performed at the Australian Centre for Ecogenomics, The University of Queensland. Gene annotation was performed using the AnnotateM script . A summary of the project is shown in Table 2 using MIGS version 2.0  criteria.
Growth conditions and genomic DNA preparation
Genomic DNA of H57 was isolated from a freeze-dried product of H57 spores combined with sodium bentonite (1:1). DNA was extracted from the H57 spores using the ‘Repeated Bead-beating and Column Extraction’ method described by Yu and Forster (2005) . In brief, 0.1 g of sporulated product was added to 1 mL of lysis buffer (2.9 % NaCl, 0.6 % Tris, 0.05 M EDTA pH 8.0 and 4 % SDS) in a cryotube containing 0.5 g zirconia beads (BioSpec Products Inc., Bartlesville, USA). The sample was then homogenised in a mini bead beater 16 (BioSpec Products Inc., Bartlesville, USA) for 2 cycles of 3 min. Between cycles the samples were incubated for 15 min at 70 °C, centrifuged (13,200 rpm for 5 min at 4 °C) and supernatant transferred to a fresh tube. Following bead beating further extraction was performed on the supernatant using the QIAGEN QIAmp DNA Mini Kit as per kit instructions (QIAGEN, Doncaster, VIC).
Genome sequencing and assembly
The genome of H57 was sequenced on an Illumina MiSeq sequencing platform (Illumina, Inc. San Diego, CA). DNA libraries were prepared using the Nextera® XT DNA Library Preparation Kit (Illumina, San Diego, CA) according to the manufacturer’s instructions. An input of 1 ng was used to prepare DNA libraries, which was then cleaned using Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA). The purified PCR product was then size selected for amplicons with a size between 300 bp and 800 bp. Illumina paired-end sequencing was performed, producing a total of 1,351,526 reads. Primer and adaptor sequences were removed using Trimmomatic v0.32  resulting in an average read length of 256 bp. Reads were assembled using SPAdes 3.0.0. . The H57 genome was obtained in 16 contigs ranging in size from 701,147 bp to 10,158 bp with a combined length of 3,958,833 bp. Genome completeness and contamination was estimated using CheckM version 1.0.0, indicating that the genome was near complete (99.51 %) with no detectable contamination (0 %) .
Gene annotation was achieved using a combination of protein databases via AnnotateM Version 6.0 . Open reading frames were initially generated using PROKKA . The resulting protein sequence was then searched against the IMG, Uniref, COG, PFAM and TIGRfam databases [16–20] to identify homologous genes. The software ProteinOrtho  was used to identify orthologous genes to other known B. amyloliquefaciens strains for further comparison. Genes unique to H57 were compared against the KEGG gene database  to identify metabolic functions.
The draft genome assembly of H57 consists of sixteen contigs totalling 3,958,833 bp and a G + C content of 46.42 %, which is likely a slight underestimate of its genome size due to unresolved collapsed repeats, primarily rRNA operons (Table 3). With a coding region of 3,549,557 bp, this assembly represents a total of 3,945 ORFs. Of those genes, 3,836 encode proteins and the remainder encode sixteen rRNAs (7 × 5S, 7 × 16S and 2 × 23S), 69 tRNAs and 24 other RNA genes (Table 3). Of the annotated genes, the majority were assigned a putative function (80.66 %) with 69.81 % assigned into Clusters of Orthologous Groups, presented in Table 4. Of the 3,945 ORFs in the H57 genome, 3,751 were inferred to be orthologous to other B. amyloliquefaciens strains, including strains CC178, DSM7, XH7, TF28, Y2, IT-45, LFB112 and B. amyloliquefaciens subsp plantarum strains UCMB5113, FZB42, NAU-B3, YAU B9601-Y2, and TrigoCor1448. Of the 194 genes unique to H57, several appear to be involved in the degradation of aromatic compounds, more specifically the breakdown of 4-hydroxyphenylacetic acid.
Insights from the genome sequence
Comparative analysis of the H57 genome indicates that its central metabolism is consistent with other strains of B. amyloliquefaciens . The presence of a complete TCA cycle and electron transport chain indicates the potential for aerobic respiration. H57 has a narGHJI operon and the transcriptional regulator fnr, suggesting that it is also capable of growing anaerobically using nitrate as an electron acceptor . This capability would be required for H57 to grow in anoxic environments.
The genome of H57 also encodes a number of enzymes involved in carbohydrate metabolism. A search against the carbohydrate-active enzyme database  reveals that H57 is dominant in glycoside hydrolase families 1, 43 and 13 (Table 5). The GH 1 and GH 43 families comprise enzymes that degrade the various sugar monomers of hemicellulose. This suggests that H57 may contribute to breaking down the less fibrous components of the plant cell wall. The abundance of GH 13 enzymes, which are a family of α-amylases, suggests that H57 also contributes to the breakdown of starch. The presence of these carbohydrate-activated enzymes alludes to the notion that H57 may assist in the digestion of animal feeds by breaking down certain polysaccharides of the plant cell wall.
Consistent with observed anti-fungal activity, the H57 genome encodes a broad range of antimicrobial compounds. These include genes for non-ribosomal synthesis of antimicrobial lipopeptides such as surfactin (srfABCD), iturin (ituABCD), bacillomycin D (bmyABC) and fengycin (fenABCDE). Surfactin is capable of inhibiting a wide range of microorganisms due to its ability to insert itself into the cell wall creating ion pores . Bacillomycin D, iturin and fengycin all have demonstrated antifungal properties primarily based on their ability to disrupt the fungal cell wall [26–28]. The genes for the expression of antibiotic polyketides are also present on the H57 genome. These include the operons mlnABCDEFGHI, dfnABCDEFGHIJ and baeEDLMNJRS, which encode macrolactin, difficidin and bacillaene respectively. These compounds inhibit a wide range of microorganisms acting chiefly on preventing protein synthesis [29–31].
The ~3.96 Mbp genome of B. amyloliquefaciens H57 reveals the basis of its antimicrobial nature and potential to survive and reproduce in anoxic animal gastrointestinal tracts. In common with other B. amyloliquefaciens strains, H57 encodes a wide range of antimicrobial compounds that explain its effectiveness as a biocontrol agent for fungi and other feed spoilage organisms. The production of these compounds may also contribute to the observed probiotic effect by inhibiting potentially pathogenic organisms creating a healthier microbial ecosystem.
Wang L-T, Lee F-L, Tai C-J, Kasai H. Comparison of gyrB gene sequences, 16S rRNA gene sequences and DNA–DNA hybridization in the Bacillus subtilis group. Int J Syst Evol Microbiol. 2007;57(8):1846–50.
Chen XH, Koumoutsi A, Scholz R, Schneider K, Vater J, Süssmuth R, Piel J, Borriss R. Genome analysis of Bacillus amyloliquefaciens FZB42 reveals its potential for biocontrol of plant pathogens. J Biotechnol. 2009;140(1–2):27–37.
Koumoutsi A, Chen X-H, Henne A, Liesegang H, Hitzeroth G, Franke P, Vater J, Borriss R. Structural and functional characterization of gene clusters directing nonribosomal synthesis of bioactive cyclic lipopeptides in Bacillus amyloliquefaciens strain FZB42. J Bacteriol. 2004;186(4):1084–96.
Yu GY, Sinclair JB, Hartman GL, Bertagnolli BL. Production of iturin A by Bacillus amyloliquefaciens suppressing Rhizoctonia solani. Soil Biol Biochem. 2002;34(7):955–63.
Huang X, Suo J, Cui Y. Optimization of antimicrobial activity of surfactin and polylysine against Salmonella enteritidis in milk evaluated by a response surface methodology. Foodborne Pathog Dis. 2011;8(3):439–43.
Dart P, Brown S. Testing hay treated with mould-inhibiting, biocontrol inoculum: microbial inoculant for hay: a report for the Rural Industries Research and Development Corporation. no. 05/103: Rural Industries Research and Development Corporation; 2005.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39(Web Server issue):W29–37.
Haroon MF. AnnotateM. https://github.com/fauziharoon/annotateM. Accessed 23rd April 2014.
Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, Tatusova T, Thomson N, Allen MJ, Angiuoli SV, et al. The minimum information about a genome sequence (MIGS) specification. Nat Biotech. 2008;26(5):541–7.
Yu Z, Forster RJ. Nucleic acid extraction, oligonucleotide probes and PCR methods. In: Makkar HPS, McSweeney CS, editors. Methods in Gut Microbial Ecology for Ruminants. Netherlands: Springer; 2005. p. 81–104.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J Comp Biol. 2012;19(5):455–77.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Peer J PrePrints. 2014;2:e554v1.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42(D1):D222–30.
Haft DH, Selengut JD, White O. The TIGRFAMs database of protein families. Nucleic Acids Res. 2003;31(1):371–3.
Markowitz VM, Chen IMA, Palaniappan K, Chu K, Szeto E, Grechkin Y, Ratner A, Jacob B, Huang J, Williams P, et al. IMG: the integrated microbial genomes database and comparative analysis system. Nucleic Acids Res. 2012;40(Database issue):D115–22.
Suzek BE, Huang H, McGarvey P, Mazumder R, Wu CH. UniRef: comprehensive and non-redundant UniProt reference clusters. Bioinformatics. 2007;23(10):1282–8.
Tatusov R, Fedorova N, Jackson J, Jacobs A, Kiryutin B, Koonin E, Krylov D, Mazumder R, Mekhedov S, Nikolskaya A, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4(1):41.
Lechner M, Findeiss S, Steiner L, Marz M, Stadler P, Prohaska S. Proteinortho: detection of (co-)orthologs in large-scale analysis. BMC Bioinformatics. 2011;12(1):124.
Gonzalez DS. KEGG. Brief Bioinform. 2002;3(3):316.
Nakano MM, Zuber P. Anaerobic growth of a "strict aerobe" (Bacillus subtilis). Annu Rev Microbiol. 1998;52:165–90.
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42(Database issue):D490–5.
Peypoux F, Bonmatin JM, Wallach J. Recent trends in the biochemistry of surfactin. Appl Microbiol Biotechnol. 1999;51(5):553–63.
Deleu M, Paquot M, Nylander T. Effect of fengycin, a lipopeptide produced by Bacillus subtilis, on model biomembranes. Biophys J. 2008;94(7):2667–79.
Maget-Dana R, Peypoux F. Iturins, a special class of pore-forming lipopeptides: biological and physicochemical properties. Toxicology. 1994;87(1–3):151–74.
Nasir MN, Besson F. Conformational analyses of bacillomycin D, a natural antimicrobial lipopeptide, alone or in interaction with lipid monolayers at the air-water interface. J Colloid Interface Sci. 2012;387(1):187–93.
Patel PS, Huang S, Fisher S, Pirnik D, Aklonis C, Dean L, Meyers E, Fernandes P, Mayerl F. Bacillaene, a novel inhibitor of procaryotic protein synthesis produced by Bacillus subtilis: production, taxonomy, isolation, physico-chemical characterization and biological activity. J Antibiot. 1995;48(9):997–1003.
Romero-Tabarez M, Jansen R, Sylla M, Lünsdorf H, Häußler S, Santosa DA, Timmis KN, Molinari G. 7-O-malonyl macrolactin A, a new macrolactin antibiotic from Bacillus subtilis active against methicillin-resistant Staphylococcus aureus, vancomycin-resistant Enterococci, and a small-colony variant of Burkholderia cepacia. Antimicrob Agents Chemother. 2006;50(5):1701–9.
Zweerink MM, Edison A. Difficidin and oxydifficidin: novel broad spectrum antibacterial antibiotics produced by Bacillus subtilis. III. Mode of action of difficidin. J Antibiot. 1987;40(12):1692–7.
Price MN, Dehal PS, Arkin AP. FastTree: Computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26(7):1641–50.
Wolfgang Ludwig OS, Ralf W, Lothar R, Harald Meier Y, Arno B, Tina L, Susanne S, Gangolf J, Wolfram F, Igor B, Stefan G, Ginhart AW, Oliver G, Silke G, Stefan H, Ralf J, Andreas K, Thomas L, Ralph L, Michael M, Björn N, Boris R, Robert S, Alexandros S, Norbert S, Alexander V, Michael L, Thomas L, Arndt B, Karl-Heinz S. ARB: a software environment for sequence data. Nucleic Acids Res. 2004;32(4):1363–71.
Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci U S A. 1990;87(12):4576–9.
Garrity G, Holt J. The road map to the manual. In: Garrity G, Boone D, Castenholz R, editors. Bergey’s manual of systematic bacteriology, vol. 1. 2nd ed. New York: Springer; 2001. p. 119–69.
Gibbons NE, Murray RGE. Proposals concerning the higher taxa of bacteria. Int J Syst Evol Microbiol. 1978;28(1):1–6.
Murray R. The higher taxa, or, a place for everything…? In: Holt J, editor. Bergey's Manual of Systematic Bacteriology, vol. 1. 1st ed. Baltimore: The Williams and Wilkins Co; 1984. p. 31–4.
Oren A, Garrity GM. List of new names and new combinations previously effectively, but not validly, published. Int J Syst Evol Microbiol. 2016;66(3):1603–6.
Ludwig W, Schleifer K, Whitman W. ClassI. Bacilli class nov. In: De Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey F, Schleifer K, Whitman W, editors. Bergey's Manual of Systematic Bacteriology, vol. 3. 2nd ed. New York: Springer; 2009. p. 19–20.
Skerman VBD, McGowan V, Sneath PHA. Approved lists of bacterial names. Int J Syst Bacteriol. 1980;30(1):225–420.
Prévot AR. Bacillales. In: Hauderoy P, Ehringer G, Guillot G, Magrou J, Prévot AR, Rosset D, Urbain A, editors. Dictionnaire des Bactéries Pathogènes. 2nd ed. Paris: Masson et Cie; Paris 1953. p. 1-692.
Fischer A. Untersuchungen über bakterien. Jahrbücher für Wissenschaftliche Botanik. 1985;27:1–163.
Gibson T, Gordon R, Genus I. Bacillus Cohn 1872, 174; Nom. gen. cons. Nomencl. Comm. Intern. Soc. Microbiol. 1937, 28; Opin. A. Jud. Comm. 1955, 39. In: Buchanan R, Gibbons N, editors. Bergey’s Manual of Determinative Bacteriology. 8th ed. Baltimore: The Williams and Wilkins Co; 1974. p. 529–50.
Cohn F. Untersuchungen über Bakterien. Beitr Biol Pflanz. 1872;1:127–224.
Fukomoto J. Studies on the production of bacterial amylase. I. Isolation of bacteria secreting potent amylase and their distribution. Nippon Nogeikagaku Kaishi. 1943;19:487–503.
Priest FG, Goodfellow M, Shute LA, Berkeley RCW. Bacillus amyloliquefaciens sp. nov., nom. rev. Int J Syst Bacteriol. 1987;37(1):69–71.
Wang L-T, Lee F-L, Tai C-J, Kuo H-P. Bacillus velezensis is a later heterotypic synonym of Bacillus amyloliquefaciens. Int J Syst Evol Microbiol. 2008;58(3):671–5.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.
The authors acknowledge the facilities, and the scientific and technical assistance, of the Australian Microscopy & Microanalysis Research Facility at the Centre for Microscopy and Microanalysis, The University of Queensland. We thank Roger Shivas of the Department of Agriculture and Fisheries, QLD for assistance in microscopy. The authors also thank Serene Low and Nicola Angel of the Australian Centre for Ecogenomics for library preparation and sequencing. This work was supported by funding from Ridley AgriProducts and the Australian Research Council (LP120200837).
BJS cultivated the bacterium, contributed to bioinformatic analysis, submitted the genome and drafted the manuscript. AS completed genome assembly, gene annotation and revised the manuscript. NL performed DNA library preparation, participated in the design of the study and revised the manuscript. DO and AVK participated in the design and supervision of the study, assisted with the interpretation of results and helped to draft the manuscript. PD isolated the bacterium, conceived of the study and participated in its design. PH participated in the studies conception, design and coordination. PH provided support with interpretation of the data and helped draft the manuscript. All authors revised and approved of the final manuscript.
The authors declare that they have no competing interests.