- Data note
- Open Access
Draft genome sequence of Enterococcus faecium SP15, a potential probiotic strain isolated from spring water
BMC Research Notesvolume 12, Article number: 99 (2019)
Enterococci are Gram-positive lactic acid bacteria and common inhabitants of the gastrointestinal tract of mammals, including humans. They are also widely distributed in diverse environments such as soil, water, vegetables and food. Enterococcus faecium is able to produce antimicrobial compounds (enterocins) and thus can act as a probiotic. E. faecium SP15 is a newly identified enterocin-producing strain from spring water that has been subjected to genome sequence analysis to provide understanding of its antimicrobial and probiotic properties.
The draft genome of E. faecium SP15 comprises of 2,783,033 bp with a G+C content of 38.08%. Five genetic loci predicted to specify enterocin production were identified, but no virulence factors could be detected and only two potential antibiotic resistance genes were noted.
Enterococci are Gram-positive lactic acid bacteria with a wide environmental distribution, encompassing many species from a wide variety of ecological niches . Enterococcus faecium is a major nosocomial pathogen often causing neonatal meningitis or endocarditis . However, certain strains of E. faecium have beneficial effects on human health due to their probiotic activity . For example, E. faecium T-110 is a syndicate member in several probiotic products including BIO-THREER which is widely prescribed for human, animal and aquacultural use  and E. faecium strain L-3 is the principle organism in the probiotic Laminolact . E. faecium is well known for its ability to produce bacteriocins, but there are relatively few reports in the literature on the genome sequence of E. faecium from non-clinical sources [5, 6] and there are no current NCBI database genome entries for E. faecium isolated from natural water sources.
Enterococcus faecium SP15 was isolated from spring water in Rawalakot (Azad Kashmir, Pakistan), a site of relative isolation with little waste water contamination. The strain exhibited strong antimicrobial activity against a panel of seven indicator strains, including Listeria monocytogenesis, indicative of enterocin production. E. faecium SP15 genomic DNA was extracted using a GeneJET genomic DNA purification kit (Thermofisher Scientific) as recommended by the vendor and assessed using a NanoDrop ND-1000 spectrophotometer and 0.7% agarose gel electrophoresis. Genome sequencing was performed by MicrobesNG (University of Birmingham, UK) using Illumina MiSeq and HiSeq 2500 platforms (Illumina, UK) with 2 × 250 bp paired-end reads. The reads were trimmed using Trimmomatic version 0.30  and the quality was assessed using in-house scripts combined with BWA-MEM software 0.7.16 . De novo assembly was performed with SPAdes software version 3.9.0  and assembly metrics were calculated using QUAST version 2.0 . Gene prediction and annotation were carried out using the Pathosystems Resources Integration Center (PATRIC) web server , RAST version 2.0  and the NCBI, PGAP version 4.6 . The tRNA genes were predicted by tRNA scan-SE 2.0 . Antimicrobial mechanisms (e.g. enterocin production) were explored with BAGEL 3  and anti-SMASH V4 . Virulence factors were identified using the virulence factor database, VFDB . Antibiotic resistance gene were identified using the Comprehensive Antibiotic Resistance Database, CARD , and acquired resistance genes were predicted by Resistance Finder 3.0 . Contigs were ordered by alignment against the most closely related sequence in GenBank  (E. faecium T110, CP006030; 99% identity) using progressive Mauve version 2.4.0 . Intact and incomplete prophage regions were identified through the integrated search and annotation tool, PHAST . Clustered regularly interspaced short palindromic repeat (CRISPR) arrays were identified using CRISPR finder .
The draft genome assembly consisted of 121 contigs with a total size of 2,783,033 bp (Table 1). The genome sequence data was at 30× coverage with an N50 of 102,590 bp and mean GC content 38.08%. A total of 2900 protein-encoding genes were predicted of which 2063 were assigned putative functions while 837 remain hypothetical. A total of 63 tRNA structural genes were identified. BAGLE 3 predicted five bacteriocin biosynthetic gene clusters (enterocin-HF, enterocin-P, enterocin SE-K4, enterocin L50A/L50B and enterolysin). No virulence factors (AS, Ace, Acm, Scm, EfaA, EcbA, Esp, Cyl, GelE and SprE) were detected. Two antibiotic resistance genes aac (6′)-li and msrC (98 and 97% identity, respectively) were identified conferring resistance to aminoglycosides, and macrolides and streptogramin B antibiotics, but no acquired resistance determinants were found. Four prophage loci were predicted of which three were intact (Strept_9871, Lactob_phig1e and Staphy_SPbeta (40, 36.9 and 32.7 kb, respectively) and one incomplete (Salmon_SJ46, 17.4 kb). One CRISPR array was identified of 190 bp, containing three spacers with a highly conserved 24 bp DR region, and although two cas gene clusters (cas3_typeI, cas4_typeI-II) were found, these were not associated with the CRISPR array.
Two major replicons are apparent: a chromosome of ~ 2,545,000 bp and a plasmid of ~ 149,300 bp (related to plasmid pNB2354 from E. faecium NRRL-B-2354, CP004064).
Current data is based on the draft level genome such that the exact length of the genome, and the number of rRNA genes and repetitive elements, cannot be absolutely determined. Furthermore, the genome includes extrachromosomal elements that cannot be predicted precisely.
prokaryotic genome annotation pipeline
- BAGEL 3:
BActeriocin GEnome mining tooL version 3
antibiotic and secondary Metabolite Analysis SHell
Franz CM, Holzapfel WH, Stiles ME. Enterococci at the crossroads of food safety. Int J Food Microbiol. 1999;47:1–24.
Linden SM, Miller CB. Vancomycin-resistant enterococci: the clinical effect of a common nosocomial pathogen. Diagn Microbiol Infect Dis. 1999;33:113–20.
Purushothaman M, Parani M. First complete genome sequence of a probiotic Enterococcus faecium strain T-110 and its comparative genome analysis with pathogenic and non-pathogenic Enterococcus faecium genome. J Genet Genomics. 2015;42:43–6. https://doi.org/10.1016/j.jgg.2014.07.002.
Karaseva A, Tsapieva A, Pachebat J, Suvorov A. Draft genome sequence of probiotic Enterococcus faecium strain L-3. Genome Announ. 2016;4:e01622.
Qin X, Galloway-Peña JR, Sillanpaa J, Roh JH, Nallapareddy SR, Chowdhury S, Bourgogne A, Choudhury T, Muzny DM, Buhay CJ, Ding Y, Dugan-Rocha S, Liu W, Kovar C, Sodergren E, Highlander S, Petrosino JF, Worley KC, Gibbs RA, Weinstock GM, Murray BE. Complete genome sequence of Enterococcus faecium strain TX16 and comparative genomic analysis of Enterococcus faecium genomes. BMC Microbiol. 2012;12:135.
Van Hal SJ, Ip CLC, Ansari MA, Wilson DJ, Espedido BA, Jensen SO, Bowden R. Evolutionary dynamics of Enterococcus faecium reveals complex genomic relationships between isolates with independent emergence of vancomycin resistance. Microbial Genomics. 2016. https://doi.org/10.1099/mgen.0.000048.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics. 2014;30(15):2114–200.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler Transform. Bioinformatics. 2009;25:1754–60.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29:1072–5.
Wattam AR, Davis JJ, Assaf R, Boisvert S, Brettin T, Bun C, Conrad N, Dietrich EM, Disz T, Gabbard JL, Gerdes S, Henry CS, Kenyon RW, Machi D, Mao C, Nordberg EK, Olsen GJ, Murphy-Olson DE, Olson R, Overbeek R, Parrello B, Pusch GD, Shukla M, Vonstein V, Warren A, Xia F, Yoo H, Stevens RL. Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center. Nucleic Acids Res. 2017;45:D535–D542542.
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O. The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75.
Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, Lomsadze A, Pruitt KD, Borodovsky M, Ostell J. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016;44(14):6614–24. https://doi.org/10.1093/nar/gkw569.
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.
van Heel AJ, de Jong A, Montalban-Lopez M, Kok J, Kuipers OP. BAGEL3: automated identification of genes encoding bacteriocins and (non-bactericidal post translationally modified peptides. Nucleic Acids Res. 2013;41:W448–W45353. https://doi.org/10.1093/nar/gkt391.
Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, et al. antiSMASH: rapid identification, annotation and analysis of secondary metabolite bio-synthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011;39:W339–W34646. https://doi.org/10.1093/nar/gkr466.
Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y, Jin Q. VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 2005;33:D325–D328328.
McArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57(7):3348–57. https://doi.org/10.1128/AAC.00419-13.
Zankari E, Hasman H, Consentino S, Vestergaard M, Rasmussen S, Lund O, Aarestrup FM, Larsen MV. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012. https://doi.org/10.1093/jac/dks261.
Agarwala R, Barrett T, Beck J, Benson DA, Bollin C, Bolton E, et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2018;46(D1):D8–D13. https://doi.org/10.1093/nar/gkx1095.
Darling AE, Mau B, Perna NT. Progressive mauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS ONE. 2010;5(6):e11147. https://doi.org/10.1371/journal.pone.0011147.
Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: a fast phage search tool. Nucleic Acids Res. 2011;39:W347–W352352.
Grissa I, Vergnaud G, Pource LC. CRISPR Finder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35:W52–W5757. https://doi.org/10.1093/nar/gkm360.
FA obtained the samples and completed the sequence analysis. MNK and SA advised on the study design and SCA advised and assisted the experimental work. All authors read and approved the final manuscript.
FA gratefully acknowledges an IRSIP award from the Higher Education Commission of Pakistan
The authors declare that they have no competing interests.
Availability of data materials
This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession number RDQA00000000. The version described in this paper is version RDQA01000000. The data described in this Data Note can be freely and openly accessed at DDBJ/ENA/GenBank. Please see Table 1 for details and links to the data.
Consent for publication
Ethics approval and consent to participate
Higher Education Commission of Pakistan (International research support initiative program, I-8HEC/HRD/2017/7639 (IRSIP 36 BMS 02)). The funding body had no role in the experimental design, collection, analysis and interpretation of data or in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.