- Research article
- Open Access
Comparative genome analysis of Pediococcus damnosus LMG 28219, a strain well-adapted to the beer environment
© Snauwaert et al.; licensee BioMed Central. 2015
- Received: 19 May 2014
- Accepted: 6 March 2015
- Published: 3 April 2015
Pediococcus damnosus LMG 28219 is a lactic acid bacterium dominating the maturation phase of Flemish acid beer productions. It proved to be capable of growing in beer, thereby resisting this environment, which is unfavorable for microbial growth. The molecular mechanisms underlying its metabolic capabilities and niche adaptations were unknown up to now. In the present study, whole-genome sequencing and comparative genome analysis were used to investigate this strain’s mechanisms to reside in the beer niche, with special focus on not only stress and hop resistances but also folate biosynthesis and exopolysaccharide (EPS) production.
The draft genome sequence of P. damnosus LMG 28219 harbored 183 contigs, including an intact prophage region and several coding sequences involved in plasmid replication. The annotation of 2178 coding sequences revealed the presence of many transporters and transcriptional regulators and several genes involved in oxidative stress response, hop resistance, de novo folate biosynthesis, and EPS production. Comparative genome analysis of P. damnosus LMG 28219 with Pediococcus claussenii ATCC BAA-344T (beer origin) and Pediococcus pentosaceus ATCC 25745 (plant origin) revealed that various hop resistance genes and genes involved in de novo folate biosynthesis were unique to the strains isolated from beer. This contrasted with the genes related to osmotic stress responses, which were shared between the strains compared. Furthermore, transcriptional regulators were enriched in the genomes of bacteria capable of growth in beer, suggesting that those cause rapid up- or down-regulation of gene expression.
Genome sequence analysis of P. damnosus LMG 28219 provided insights into the underlying mechanisms of its adaptation to the beer niche. The results presented will enable analysis of the transcriptome and proteome of P. damnosus LMG 28219, which will result in additional knowledge on its metabolic activities.
- Pediococcus damnosus
- Illumina sequencing
- Beer environment
- Draft genome
- Comparative genomics
Beer is a fermented beverage that is high in ethanol, carbon dioxide and flavorful yeast metabolites, contains hop-derived flavor and antimicrobial compounds, and is low in pH, oxygen, and residual nutrients . This environment has selected for unique groups of bacteria specialized in growth in beer, including several species of lactic acid bacteria (LAB). Overall, one of the key metabolic actions of LAB is the reduction of pyruvate into lactate to regenerate NAD+ by means of lactate dehydrogenase activity. Depending on the beer type produced, lactate production may or may not be desired [2-5]. As most LAB species tolerate high ethanol concentrations and a low pH, hop resistance is the major factor limiting their growth in beer . Hop-derived antimicrobial compounds, known as iso-alpha acids, cause permeability changes in the bacterial cell wall , leakage of the cytoplasmic membrane and subsequent inhibition of respiration and protein, DNA, and RNA syntheses , as well as changes in leucine uptake and proton ionophore activity . Additionally, iso-alpha acids alter the redox properties of the bacterial cell, causing oxidative stress . Accordingly, hop resistance is a multifactorial property that requires different mechanisms to counteract the action of iso-alpha acids. A well-known factor mediating hop resistance is the ATP-binding cassette multidrug transporter, HorA, which extrudes iso-alpha acids  and is plasmid-encoded .
LAB species frequently encountered in the beer environment belong to the genera Lactobacillus and Pediococcus, which are both Gram-positive, facultative anaerobic, chemo-organotrophic bacteria of the family Lactobacillaceae of the Firmicutes phylum. The species with the highest capacity to grow in beer are Lactobacillus brevis and Pediococcus damnosus , although the ability of bacteria to grow in beer is a strain- rather than species-specific characteristic . Recent studies have investigated the mechanisms of L. brevis for overcoming stresses in beer by means of reverse transcription quantitative polymerase chain reaction (RT-qPCR) and proteomic analyses [13,14]. Less is known about the adaptations of P. damnosus to the beer environment. So far, Pediococcus claussenii ATCC BAA-344T originating from spoiled beer, Pediococcus pentosaceus ATCC 25745 of plant origin, and P. pentosaceus IE-3 originating from a dairy effluent sample, are the only members of the genus Pediococcus for which the genome has been sequenced completely [15-17]. However, draft genome sequences are available for several strains of Pediococcus acidilactici, namely strain MA 18/5 M originating from animal food, strain DSM 20284 isolated from barley, strain 7_4 from a human fecal sample, and strain NGRI 0510Q, formerly known as Pediococcus lolii , isolated from ryegrass silage (http://genomesonline.org). In general, Pediococcus species possess rather small genomes (approximately 2 Mb), encoding a broad repertoire of transporters for efficient carbon and nitrogen acquisition and reflecting a limited range of biosynthetic capabilities. This suggests both extensive gene loss as well as acquisitions via horizontal gene transfer during the evolution of pediococci within their habitats . Several strains possess plasmids containing genes regulating the fermentation of carbohydrates and encoding different types of resistances (e. g., resistance to stress, hop, or antibiotics) [17,19].
The strain P. damnosus LMG 28219 was isolated in 2013 from a Flemish acid beer, i.e. a sour (containing both lactic acid and acetic acid) and ethanolic [5-6% (vol/vol)] beer at the end of its maturation phase. A microbial diversity analysis of this mature beer proved that this strain was present in high numbers and is therefore able to replicate in this environment (unpublished results). Flemish acid beers are produced by a mixed-culture fermentation and represent culturally important products, for which microbial activities play critical roles in beer production and quality formation [3,5]. In the present study, the draft genome sequence of P. damnosus LMG 28219 is presented and analysed to obtain insights into its genome-based metabolic features. A better understanding of the molecular mechanisms underlying its metabolic capabilities enabled detailed insights into the mechanisms of adaptation of this strain to the beer environment. Furthermore, comparison of P. damnosus LMG 28219 with other sequenced genomes of members of the genus Pediococcus addressed the potentially unique properties of this strain and other strains adapted to the beer environment.
General architecture and annotation of the Pediococcus damnosus LMG 28219 draft genome
Overview of the Pediococcus damnosus LMG 28219 draft genome and the best BLAST hits
Draft genome of P. damnosus LMG 28219
Mean GC content (%)
Mean size consensus (bp)
Best BLAST hit
82, 98, 145, 116, 77, 119, 47, 53, 5, 89, 56, 76, 166, 84, 81, 136, 102, 141, 72, 30, 52, 129, 74, 54, 19, 110, 147, 49, 26, 85, 137, 61, 146, 177, 7, 68, 131, 104, 87, 32, 17, 80, 66, 165, 2, 20, 101, 96, 4, 59, 45, 164, 118, 64, 176, 69, 90, 24, 127, 10, 46, 139, 114, 60, 124, 23, 83, 51, 151, 169, 106, 11, 121, 78, 55, 43, 130, 48, 65, 168, 27, 62, 60, 79, 6, 152, 120, 73, 35, 99, 29, 71*
chromosomal DNA (NC_016605)
P. claussenii ATCC BAA-344 T
21, 39, 41, 183
12, 21, 94, 95, 112, 135
1,21, 39, 40, 95, 97, 107, 132, 135, 170
21, 107, 135, 178
8,21, 39, 41, 86, 94, 182, 183
Bacillus megaterium QM B1551 (pBM400, NC_004604) and WSH-002 (WSH-002_p1, CP003018)
123, 134, 140
L. brevis 925A (pLB925A04, NC_012551), KB290 (pKB290-4, AP012171), and ATCC 367 (plasmid 1, NC_008498)
L. plantarum WCFS1 (pWCFS103, NC_006377)
L. buchneri CD034 (pCD034-3, CP003044) and NRRL B-30929 (pLBUC01, NC_015420)
18, 22, 93, 108, 109, 113, 115, 143, 158, 162, 163
Several Lactobacillus spp.
A BLAST search of the remaining non-chromosomally encoded contigs against a plasmid-specific database (i.e., PATRIC) revealed many similarities towards plasmids of different bacterial species (see below). Many of these contigs (as listed in Table 1) were similar to the plasmids of P. claussenii ATCC BAA-344T, with the exception of plasmids pPECL-1 and pPECL-2, onto which no contigs mapped. These two plasmids are small and cryptic (1.85Kb and 2.45Kb in size, respectively), whereas the other six plasmids of P. claussenii ATCC BAA-344T range from 16 to 36 kb and contribute to roughly 7% of the strain’s coding capacity . Four, six, ten, four, eight, and three contigs mapped onto plasmids pPECL-3, pPECL-4, pPECL-5, pPECL-6, pPECL-7, and pPECL-8 of P. claussenii ATCC BAA-344T, respectively. A total of 19 contigs did not show high similarities towards the plasmids of P. claussenii ATCC BAA-344T. Contigs 16 and 57 were similar to plasmids pBM400 and WSH-002_p1 of Bacillus megaterium strains QM B1551 and WSH-002, respectively, whereas others were similar to plasmids of L. brevis strains 925A, KB290, and ATCC 367 (contigs 123, 134, and 140, respectively), Lactobacillus plantarum WCFS1 (contig 37), Lactobacillus buchneri strains CD034 and NRRL B-30929 (contigs 100 and 122, respectively), or other Lactobacillus species (all remaining contigs). In addition, several CDSs involved in plasmid replication were found, which encoded replication initiation proteins (AH70_10315 on contig 8, AH70_03600 on contig 21, AH70_00575 on contig 107, AH70_02015 on contig 138, and AH70_03140 on contig 182).
Gene finding and annotation of the P. damnosus LMG 28219 draft genome with the Integrated Microbial Genomics Expert review (IMG-RE) software resulted in a total gene count of 2266, among which 96.12% were CDSs. A total of 79.08% of the CDSs could be assigned to a protein, among which 23.39% and 12.80% were predicted to be enzymes and transporters, respectively. This agrees with the property of LAB to encode a broad repertoire of transporters for efficient nutrient uptake and reflects their limited biosynthetic capabilities . Additionally, 2.60% of the genes encoded signal peptides, whereas 26.04% encoded transmembrane proteins. Furthermore, a total of 56 tRNA and three 5S, 16S, and 23S rRNA genes were predicted. As expected, most CDSs were involved in carbohydrate, protein, DNA, and RNA metabolism and in cell wall and capsule construction. These findings agreed with those of Makarova and colleagues , who analyzed 12 genomes belonging to the order Lactobacillales and defined a core of Lactobacillales-specific clusters of orthologous groups (LaCOGs). The functional distribution of the conserved core of 567 LaCOGs showed that the majority encodes components of the information processing systems (translation, transcription, and replication), which are likely to perform essential functions. With the exception of 9 LaCOGs (i.e., COGs 0230, 0419, 0454, 0457, 0470, 0762, 2815, 2855, and 4608), all core LaCOGs were found in the draft genome of P. damnosus LMG 28219. Genome closure could possibly lead to the identification of these missing core functions, although an update of the core LaCOGs is mandatory, because only 12 Lactobacillales species out of more than 500 species of this order were included; in addition, only one of these 12 species belonged to the genus Pediococcus.
Next to these core functions, the draft genome sequence of P. damnosus LMG 28219 comprised 50 predicted CDSs that were associated with virulence, disease, and defense, including three multidrug resistance efflux pumps. Also, 42 CDSs related to stress response were detected, seven of which were involved in oxidative stress response. Four of the non-core LaCOGs present in the P. damnosus LMG 28219 draft genome comprised 12 or more CDSs in a COG, namely arabinose efflux permeases (COG2814), transcriptional regulators (COG1309), predicted transcriptional regulators (COG0789), and transposases and their inactivated derivatives (COG2826). Arabinose efflux permeases belong to the major facilitator superfamily, a large and diverse group of secondary transporters that includes uniporters, symporters, and antiporters . A total of 12 and 10 CDSs of COG2814 were annotated as arabinose efflux permeases and drug resistance transporters from the EmrB/QacA subfamily, respectively. In addition, one CDS was predicted to be a drug resistance transporter of the Bcr/CflA subfamily and one was predicted to be a transporter belonging to the carbohydrate transporter family. The COG1309 members belonged to the TetR family of transcriptional repressors, which are involved in transcriptional control of multidrug efflux pumps, pathways for the biosynthesis of antibiotics, responses to osmotic stress and toxic chemicals, control of catabolic pathways, differentiation processes, and pathogenicity. Analysis of the CDSs of COG0789 did not reveal insights into their function. Finally, the CDSs in COG2826 were all transposases belonging to the insertion sequence IS30 family. The presence of IS elements may lead to genomic instability due to its role in gene loss and gene gain .
Comparison of genome characteristics of P. pentosaceus ATCC 25745, P. claussenii ATCC BAA-344 T , and P. damnosus LMG 28219 (based on the IMG-ER server)
P. pentosaceus ATCC 25745
P. claussenii ATCC BAA-344 T
P. damnosus LMG 28219
Accession number *
Genes involved in hop resistance
Interestingly, several studies have revealed that the horA gene is consistently surrounded by the same set of genes in bacteria growing in beer, suggesting that these entire regions are acquired via horizontal gene transfer [21,28] and could possibly play a role in their adaptation to the beer niche. Most of these horA-surrounding genes are involved in phospholipid or cell wall biosynthesis. Orthologs of PECL_1950, PECL_1952, and PECL_1954 flanking the horA gene were found in the genome of P. damnosus LMG 28219 (AH70_09615 on contig 70 and AH70_02080 and AH70_02070 on contig 14, respectively) but not in the P. pentosaceus ATCC 25745 genome sequence. AH70_09615 and AH70_02080 were annotated as an acyl-phosphate glycerol 3-phospate and a glycosyl transferase, respectively. Both PECL_1950 and PECL_1952 are members of the lysophospholipid acyltransferase superfamily, which contains acyltransferases of de novo and remodeling pathways of glycerophospholipid biosynthesis. Finally, AH70_02070 contained a family 8 glycosyl transferase domain that catalyzes the transfer of sugar moieties from activated donor molecules to specific acceptor molecules, forming glycosidic bonds.
In addition to the horA gene, the hitA gene encodes a divalent cation transporter and confers hop resistance by importing manganese and thereby counteracting proton gradient dissipation . Orthologs of the hitA gene of L. brevis L5784 (AB035808) were found in the genomes of P. damnosus LMG 28219 (AH70_05540), P. claussenii ATCC BAA-344T (AEV94625.1), and P. pentosaceus ATCC 25745 (YP_805164). The reason for the presence of the hitA gene in P. pentosaceus ATCC 25745 of plant origin remains unclear. Furthermore, HorC has been suggested to be a proton motive force-dependent multidrug transporter, whose gene expression is under control of the HorB transcriptional regulator . Orthologs of the horC gene of Lactobacillus backii LMG 23555 (BAF56899.1) were found on contig 18 (AH70_03090) in the P. damnosus LMG 28219 and P. claussenii ATCC BAA-344T genome sequences (AEV95756.1), but not in that of P. pentosaceus ATCC 25745. Finally, horB, bsrA, and bsrB homologs were not found in the draft genome sequence of P. damnosus LMG 28219 .
The beer isolates P. damnosus LMG 28219 and P. claussenii ATCC BAA-344T contained a set of genes involved in folate biosynthesis that were absent in the genome of P. pentosaceus ATCC 25745 (plant origin). These included a 5,10-methylenetetrahydrofolate reductase (EC 220.127.116.11, AH70_04245 on contig 29), a dihydropteroate synthase (EC 18.104.22.168, AH70_05560 on contig 4), a 2-amino-4-hydroxy-6-hydroxymethyldihydropteridine pyrophosphokinase (EC 22.214.171.124, AH70_05580 on contig 4), a GTP cyclohydrolase I type 1 (EC 126.96.36.199, AH70_05575 on contig 4), and a dihydroneopterin aldolase (EC 188.8.131.52, AH70_05585 on contig 4).
Laboratory experiments indicated that P. damnosus LMG 28219 produces EPS. The ability of a strain to produce EPS is not directly correlated to its ability to reside in beer but probably has importance in biofilm formation , thereby enabling persistence in the brewery environment. The EPS produced by P. claussenii ATCC BAA-344T is a high-molecular-mass β-glucan produced by the action of a transmembrane glycosyl transferase (gtf) gene. This gene is located on plasmid pPECL-7, which is not essential for growth in beer . The fibrillar polymer consists of a trisaccharide repeating unit with a β-1,3-linked glucose backbone and branches made up of single β-1,2-linked D-glucopyranosyl residues. Walling and colleagues  reported a glucosyl transferase gene (dps) in P. damnosus IOEB8801, originating from wine, that produces a linear backbone of 3-β-D-glucose-1 moieties. Surprisingly, no homologies towards the gtf and dps genes were found in the draft genome sequence of P. damnosus LMG 28219. Yet, contig 56 harbored a CDS (AH70_07835), containing a glycosyl transferase (group 2) domain, which may be involved in EPS production. OrthoMCL analysis indicated that this CDS was unique to P. damnosus LMG 28219. Other CDS candidates involved in EPS production by P. damnosus LMG 28219 were AH70_01405, AH70_01410, and AH70_01415, which were all predicted to be present on contig 122. These proteins did not show amino acid sequence homologies of more than 70% towards the proteomes of P. claussenii ATCC BAA-344T and P. pentosaceus ATCC 25745. AH70_01405 is an EPS biosynthesis protein, consisting of an AAA domain containing a P-loop, and showed similarities towards the capsular EPS family protein (EHO53752) found in the genome of Lactobacillus kisonensis F0435, originating from the human oral cavity . AH70_01410 is a lipopolysaccharide biosynthesis domain that is involved in the biosynthesis of EPS. This protein showed more than 70% protein sequence homology towards the chain length determinant protein (EEI18212) of L. buchneri ATCC 11577, isolated from the human oral cavity . Finally, AH70_01415 harbored a cell envelope-related transcriptional attenuator domain, which describes a domain of unknown function that is found in the predicted extracellular domain of a number of putative membrane-bound proteins . One of those is CpsA, a putative regulatory protein involved in EPS biosynthesis . AH70_01415 showed more than 70% protein sequence homology towards a cell envelope-like function transcriptional attenuator common domain protein (EEI18211) of L. buchneri ATCC 11577 and a biofilm regulatory protein of L. kisonensis F0435 (see above). Besides the presence of proteins potentially involved in EPS production in P. damnosus LMG 28219 (as discussed above), the mechanism of EPS production in P. damnosus LMG 28219 remains unclear.
Genes involved in oxidative stress response
Pittet and coworkers  found a set of highly transcribed genes in P. claussenii ATCC BAA-344T during growth in beer as a response to the oxidative stress imposed by hops. These genes are manganese transport proteins, methionine sulfoxide reductases MsrA and MsrB as well as other metal transport and homeostasis proteins. Orthologs of these genes were found in the genomes of P. damnosus LMG 28219 and P. pentosaceus ATCC 25745, indicating that the presence of these genes is not unique to strains capable of growing in beer. These results indicate that growth in beer is a multifactorial reaction of the cell towards a challenging environment, not only involving the presence of specific genes but also the up-and down-regulation of specific sets of genes.
COGs enriched in the genomes of bacteria originating from beer
Overview of the functions enriched in the genomes of lactic acid bacteria (LAB) of beer origin compared to those of LAB not originating from beer (based on IMG-ER)
No beer origin
L. delbrueckii bulgaricus ATCC 11842
P. acidilactici 7_4
P. acidilactici DSM 20284
P. acidilactici MA 18/5 M
P. acidilactici NGRI 0510Q
P. pentosaceus ATCC 25745
P. pentosaceus IE-3
P. claussenii ATCC BAA-344 T
P. damnosus LMG 28219
L. malefermentans KCTC 3548
L. rhamnosus ATCC 8530
Arabinose efflux permease
Predicted transcriptional regulator
Transposase and inactivated derivatives, IS30 family
Predicted hydrolases or acyltransferases (alpha/beta hydrolase superfamily)
Response regulators consisting of a CheY-like receiver domain and a winged-helix DNA-binding domain
The draft genome of P. damnosus LMG 28219 and its comparative analysis with the genomes of other pediococci provided insights into the adaptation of this strain to the beer environment. These adaptations included the presence of the horA gene and its surrounding genes, genes involved in de novo folate biosynthesis, genes involved in the production of EPS, and the enrichment of functions related to transcriptional regulation. Furthermore, the results presented in this study will enable future transcriptome analysis of P. damnosus LMG 28219, which can provide additional insights into its metabolic activities.
Bacterial strain and growth conditions
A mature Belgian red-brown acidic ale brew sample was serially diluted in 0.86% (wt/vol) saline and 50 μL of each dilution was plated on de Man-Rogosa-Sharpe (MRS, Oxoid, Basingstoke, Hampshire, United Kingdom) agar medium, supplemented with 5 ppm amphotericin B (Sigma-Aldrich) and 200 ppm cycloheximide (Sigma-Aldrich) to favor the growth of LAB and inhibit yeast growth. Agar plates were incubated at 28°C in an aerobic atmosphere for 4–6 days. LAB isolates were identified through sequence analysis of the 16S rRNA and phenylalanyl-tRNA synthase alpha subunit (pheS) genes . Strain P. damnosus LMG 28219 was present in large numbers in the mature brew sample (data not shown) and was used for genome sequencing. To obtain cell pellets, P. damnosus LMG 28219 was propagated in 750 ml de Man-Rogosa-Sharpe broth (MRS-B; Oxoid) at 28°C for 3 days, followed by microcentrifugation (11000 rpm, 15 min, 4°C).
DNA extraction and Illumina sequencing
Total DNA was extracted using a procedure applied by Gevers and colleagues , with several modifications. (i) The thawed pellet was washed in 1 ml TES buffer (6.7% sucrose, 50 mM Tris–HCl, pH 8.0, 1 mM EDTA) and resuspended in 275 μl STET buffer (8.0% sucrose, 5.0% Triton X-100, 50 mM Tris–HCl, pH 8.0, 50 mM EDTA). (ii) Ninety μl lysozyme-mutanolysin-proteinase K solution [TES buffer containing 1667 U/ml mutanolysine (Sigma-Aldrich, St. Louis, MO, USA), 33 mg/ml lysozyme (Serva, Heidelberg, Germany), and 2.78 mg/ml proteinase K (Merck, Darmstadt, Germany)] was added and the suspension was incubated at 37°C for 1 h. (iii) Prior to extraction with a phenol/chloroform/isoamylalcohol (49.5:49.5:1.0) solution, proteins were precipitated by adding 250 μl ammonium acetate. (iv) Five μl RNase (2 mg/l, Sigma-Aldrich) was added to the DNA solution, which was incubated at 37°C for 60 min and stored at −20°C. The integrity, concentration, and purity of the DNA isolated was evaluated using 1.0% (wt/vol) agarose gels, stained in ethidium bromide, and by spectrophotometric measurements at 234, 260, and 280 nm. A fluorescent stain-based kit (Qubit® dsDNA Broad range assay kit; Life Technologies, Carlsbad, CA, USA) was used for an accurate determination of the DNA concentration. Because of the low DNA concentration obtained, the extraction procedure was performed in triplicate and the resulting DNAs were pooled, after which an additional purification [by adding 12.5 μl sodium acetate (3.0 M, pH 4.8) and 200 μl 95% ice-cold ethanol to 100 μl DNA product] and quantification step was performed.
Illumina sequencing and genome assembly
Library preparation and genome sequencing was performed by BaseClear BV (Leiden, The Netherlands). A paired-end DNA library with a mean gap length size between 230 and 360 bp was sequenced, with average reads of 100 bp, on an Illumina HiSeq2500 apparatus (Illumina Inc., San Diego, CA, USA). Quality trimmed reads were obtained by quality trimming in CLC Genomics Workbench v6.5.1 (CLC Inc., Aarhus, Denmark) based on a quality score limit of 0.05 and a maximum number of ambiguous nucleotides of 2. Quality trimming was based on the Phred quality scores (Q) from the fastq read files as follows. The Q value was used to calculate an error probability (perror = 10Q/-10). Next, for every base a new value was calculated (limit – perror). This value was negative for low-quality bases. For every base, the CLC workbench calculated the running sum of this value. If the sum dropped below zero, it was set to zero. The part of the sequence not trimmed was the region ending at the highest value of the running sum and started at the last zero value before this highest score. Everything in front and behind this region was trimmed. With the ambiguous trim set to 2, the algorithm found the maximum length region containing 2 or fewer ambiguities and then trimmed away the ends not included in this region. Trimmed reads shorter than 15 bp were discarded. After trimming, an initial de novo assembly was performed in CLC Genomics Workbench v6.5.1, using the quality-trimmed paired reads. All contigs shorter than 500 bp were discarded. Contigs were ordered automatically with MAUVE  by comparing the draft genome sequence of P. damnosus LMG 28219 with the complete chromosomal DNA of P. claussenii ATCC BAA-344T. Contigs that did not show similarities towards the P. claussenii ATCC BAA-344T chromosomal DNA were blasted against a plasmid and phage database, using the PATRIC (Pathosystems Resource Integration Center)  and PHAST (PHAge Search Tool)  websites, respectively. In addition, the average GC content of the contigs was calculated using an in house developed Python script.
Functional annotation and metabolic reconstruction were performed with (i) the Rapid Annotation Subsystem Technology (RAST) server , using GLIMMER  for gene calling and allowing frameshift corrections, backfilling of gaps, and automatic fixing of errors; (ii) the Integrated Microbial Genomics Expert review (IMG-ER) annotation pipeline, using GenePRIMP for gene prediction ; and (iii) the National Center for Biotechnology Information’s Prokaryotic Genome Automatic Annotation Pipeline (PGAAP), which uses GeneMark and GeneMark.HMM for gene calling . The automated gene prediction and annotation of PGAAP was followed by manual curation of the CDSs of interest, using the results of the IMG-ER and RAST server, BLASTp , UniProt (http://www.uniprot.org/), and InterProScan .
Comparative genomics of the draft genome sequence of P. damnosus LMG 28219 was performed using the available complete genome sequences of P. claussenii ATCC BAA-344T and P. pentosaceus ATCC 25745. The average nucleotide identities (ANIs) between the draft genome (P. damnosus LMG 28219) and both reference genomes were calculated using the in silico DNA-DNA hybridization method implemented in the Jspecies software , using BLAST, as proposed by Goris and colleagues . The OrthoMCL tool was used with default settings for ortholog finding in the genomes of P. damnosus LMG 28219, P. claussenii ATCC BAA-344T, and P. pentosaceus ATCC 25745 (http://orthomcl.org/). The Kyoto encyclopedia of genes and genomes (KEGG) database was used for the reconstruction of the folate biosynthesis metabolic map (http://www.genome.jp/kegg/).
Availability of supporting data
The raw sequence data received from BaseClear BV were deposited at the Short Read Archive (SRA) of GenBank (accession number SRP035530). This Whole Genome Shotgun project was deposited at DDBJ/EMBL/GenBank under the accession number JANK00000000 after automatic annotation by PGAAP and manual curation of the genes of interest.
The authors acknowledge the financial support of the Research Foundation-Flanders (FWO-Vlaanderen), the Research Council of Ghent University (BOF project) and the Vrije Universiteit Brussel (SRP, IRP, and IOF projects), and of the Hercules Foundation. In addition, the authors would like to acknowledge the brewery involved and Charlotte Peeters and Bart Verheyde for the bioinformatics support.
- Bokulich NA, Bamforth CW. The microbiology of malting and brewing. Microbiol Mol Biol Rev. 2013;77:157–72.View ArticlePubMed CentralPubMedGoogle Scholar
- Bokulich NA, Bamforth CW, Mills DA. Brewhouse-resident microbiota are responsible for multi-stage fermentation of American coolship ale. PLoS One. 2012;7:e35507.View ArticlePubMed CentralPubMedGoogle Scholar
- Martens H, Iserentant D, Verachtert H. Microbiological aspects of a mixed yeast-bacterial fermentation in the production of a special Belgian acidic ale. J Inst Brew. 1997;103:85–91.View ArticleGoogle Scholar
- Van Oevelen D, Spaepen M, Timmermans P, Verachtert H. Microbiological aspects of spontaneous wort fermentation in the production of lambic and gueuze. J Inst Brew. 1977;83:356–60.View ArticleGoogle Scholar
- Spitaels F, Wieme AD, Janssens M, Aerts M, Daniel H-M, Van Landschoot A, et al. The microbial diversity of traditional spontaneously fermented lambic beer. PLoS One. 2014;9:e95384.View ArticlePubMed CentralPubMedGoogle Scholar
- Shimwell JL. On the relation between the staining properties of bacteria and their reaction towards hop antiseptic. J Inst Brew. 1937;43:111–8.View ArticleGoogle Scholar
- Teuber M, Schmalreck AF. Membrane leakage in Bacillus subtilis 168 induced by the hop constituents lupulone, humulone, isohumulone and humulinic acid. Arch Mikrobiol. 1973;94:159–71.View ArticlePubMedGoogle Scholar
- Simpson WJ. Ionophoric action of trans-isohumulone on Lactobacillus brevis. J Gen Microbiol. 1993;139:1041–5.View ArticleGoogle Scholar
- Behr J, Vogel RF. Mechanisms of hop inhibition include the transmembrane redox reaction. Appl Environ Microbiol. 2010;76:142–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Sakamoto K, Margolles A, van Veen HW, Konings WN. Hop resistance in the beer spoilage bacterium Lactobacillus brevis is mediated by the ATP-binding cassette multidrug transporter HorA. J Bacteriol. 2001;183:5371–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Sakamoto K, Konings WN. Beer spoilage bacteria and hop resistance. Int J Food Microbiol. 2003;89:105–24.View ArticlePubMedGoogle Scholar
- Pittet V, Morrow K, Ziola B. Ethanol tolerance of lactic acid bacteria, including relevance of the exopolysaccharide gene gtf. J Am Soc Brew Chem. 2011;69:57–61.Google Scholar
- Bergsveinson J, Pittet V, Ziola B. RT-qPCR analysis of putative beer-spoilage gene expression during growth of Lactobacillus brevis BSO 464 and Pediococcus claussenii ATCC BAA-344T in beer. Appl Microbiol Biotechnol. 2012;96:461–70.View ArticlePubMedGoogle Scholar
- Behr J, Israel L, Gänzle MG, Vogel RF. Proteomic approach for characterization of hop-inducible proteins in Lactobacillus brevis. Appl Environ Microbiol. 2007;73:3300–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Midha S, Ranjan M, Sharma V, Kumari A, Singh PK, Korpole S, et al. Genome sequence of Pediococcus pentosaceus strain IE-3. J Bacteriol. 2012;194:4468.View ArticlePubMed CentralPubMedGoogle Scholar
- Pittet V, Abegunde T, Marfleet T, Haakensen M, Morrow K, Jayaprakash T, et al. Complete genome sequence of the beer spoilage organism Pediococcus claussenii ATCC BAA-344T. J Bacteriol. 2012;194:1271–2.View ArticlePubMed CentralPubMedGoogle Scholar
- Makarova K, Slesarev A, Wolf Y, Sorokin A, Mirkin B, Koonin E, et al. Comparative genomics of the lactic acid bacteria. Proc Natl Acad Sci U S A. 2006;103:15611–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Wieme A, Cleenwerck I, Van Landschoot A, Vandamme P. Pediococcus lolii DSM 19927T and JCM 15055T are strains of Pediococcus acidilactici. Int J Syst Evol Microbiol. 2012;62:3105–8.View ArticlePubMedGoogle Scholar
- Haakensen M, Ziola B. Identification of novel horA-harbouring bacteria capable of spoiling beer. Can J Microbiol. 2008;54:321–5.View ArticlePubMedGoogle Scholar
- Yoon B, Jang S, Chang H-I. Sequence analysis of the Lactobacillus temperate phage Sha1. Arch Virol. 2011;156:1681–4.View ArticlePubMedGoogle Scholar
- Pittet V, Phister TG, Ziola B. Transcriptome sequence and plasmid copy number analysis of the brewery isolate Pediococcus claussenii ATCC BAA-344T during growth in beer. PLoS One. 2013;8:e73627.View ArticlePubMed CentralPubMedGoogle Scholar
- Reddy VS, Shlykov MA, Castillo R, Sun EI, Saier MH. The major facilitator superfamily (MFS) revisited. FEBS J. 2012;279:2022–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Darmon E, Leach DRF. Bacterial genome instability. Microbiol Mol Biol Rev. 2014;78:1–39.View ArticlePubMed CentralPubMedGoogle Scholar
- Klappenbach JA, Dunbar JM, Schmidt TM. rRNA operon copy number reflects ecological strategies of bacteria. Appl Environ Microbiol. 2000;66:1328–33.View ArticlePubMed CentralPubMedGoogle Scholar
- Di Mattia E, Grego S, Cacciari I. Eco-physiological characterization of soil bacterial populations in different states of growth. Microb Ecol. 2002;43:34–43.PubMedGoogle Scholar
- De Bruyne K, Franz CMAP, Vancanneyt M, Schillinger U, Mozzi F, de Valdez GF, et al. Pediococcus argentinicus sp. nov. from Argentinean fermented wheat flour and identification of Pediococcus species by pheS, rpoA and atpA sequence analysis. Int J Syst Evol Microbiol. 2008;58:2909–16.View ArticlePubMedGoogle Scholar
- Suzuki K, Ozaki K, Yamashita H. Comparative analysis of conserved genetic markers and adjacent DNA regions identified in beer-spoilage lactic acid bacteria. Lett Appl Microbiol. 2004;39:240–5.View ArticlePubMedGoogle Scholar
- Iijima K, Suzuki K, Asano S, Kuriyama H, Kitagawa Y. Isolation and identification of potential beer-spoilage Pediococcus inopinatus and beer-spoilage Lactobacillus backi strains carrying the horA and horC gene clusters. J Inst Brew. 2007;113:96–101.View ArticleGoogle Scholar
- Hayashi N, Ito M, Horiike S, Taguchi H. Molecular cloning of a putative divalent-cation transporter gene as a new genetic marker for the identification of Lactobacillus brevis strains capable of growing in beer. Appl Microbiol Biotechnol. 2001;55:596–603.View ArticlePubMedGoogle Scholar
- Iijima K, Suzuki K, Ozaki K, Yamashita H. HorC confers beer-spoilage ability on hop-sensitive Lactobacillus brevis ABBC45cc. J Appl Microbiol. 2006;100:1282–8.View ArticlePubMedGoogle Scholar
- Suzuki K, Iijima K, Sakamoto K, Sami M, Yamashita H. A review of hop resistance in beer spoilage lactic acid bacteria. J Inst Brew. 2006;112:173–91.View ArticleGoogle Scholar
- de Crecy-Lagard V, El Yacoubi B, de la Garza RD, Noiriel A, Hanson AD. Comparative genomics of bacterial and plant folate synthesis and salvage: predictions and validations. BMC Genomics. 2007;8:245.View ArticlePubMed CentralPubMedGoogle Scholar
- Walling E, Gindreau E, Lonvaud-Funel A. A putative glucan synthase gene dps detected in exopolysaccharide-producing Pediococcus damnosus and Oenococcus oeni strains isolated from wine and cider. Int J Food Microbiol. 2005;98:53–62.View ArticlePubMedGoogle Scholar
- Chen T, Yu W-H, Izard J, Baranova OV, Lakshmanan A, Dewhirst FE. The human oral microbiome database: a web accessible resource for investigating oral microbe taxonomic and genomic information. Database 2010:baq013.Google Scholar
- Tilden EB, Svec M. Further studies of a differential culture technique for estimations of acidogenic bacteria in saliva. II. Species of lactobacilli isolated from saliva and their distribution in a group of children. J Dent Res. 1952;31:831–8.View ArticlePubMedGoogle Scholar
- Zdobnov EM, Apweiler R. InterProScan – an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001;17:847–8.View ArticlePubMedGoogle Scholar
- Gevers D, Huys G, Swings J. Applicability of rep-PCR fingerprinting for identification of Lactobacillus species. FEMS Microbiol Lett. 2001;205:31–6.View ArticlePubMedGoogle Scholar
- Darling ACE, Mau B, Blattner FR, Perna NT. MAUVE: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14:1394–403.View ArticlePubMed CentralPubMedGoogle Scholar
- Wattam AR, Abraham D, Dalay O, Disz TL, Driscoll T, Gabbard JL, et al. PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic Acids Res. 2014;42:D581–91.View ArticlePubMed CentralPubMedGoogle Scholar
- You Z, Yongjie L, Karlene HL, Jonathan JD, David SW. PHAST: A fast phage search tool. Nucleic Acids Res. 2011;39:W347–52.View ArticleGoogle Scholar
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75.View ArticlePubMed CentralPubMedGoogle Scholar
- Delcher AL, Harmon D, Kasif S, White O, Salzberg SL. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999;27:4636–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Mavromatis K, Chu K, Ivanova N, Hooper SD, Markowitz VM, Kyrpides NC. Gene context analysis in the integrated microbial genomes (IMG) data management system. PLoS One. 2009;4:e7979.View ArticlePubMed CentralPubMedGoogle Scholar
- Lukashin AV, Borodovsky M. GeneMark.HMM: new solutions for gene finding. Nucleic Acids Res. 1998;26:1107–15.View ArticlePubMed CentralPubMedGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.View ArticlePubMed CentralPubMedGoogle Scholar
- Rosselló-Mora R. DNA-DNA reassociation methods applied to microbial taxonomy and their critical evaluation. In: Stackebrandt E, editor. In Molecular Identification, Systematics, and Population Structure of Prokaryotes. Berlin, Germany: Springer; 2006. p. 23–50.View ArticleGoogle Scholar
- Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. DNA–DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol. 2007;57:81–91.View ArticlePubMedGoogle Scholar
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.