- Research article
- Open Access
Comparative genome analysis of the Lactobacillus brevis species
BMC Genomics volume 20, Article number: 416 (2019)
Lactobacillus brevis is a member of the lactic acid bacteria (LAB), and strains of L. brevis have been isolated from silage, as well as from fermented cabbage and other fermented foods. However, this bacterium is also commonly associated with bacterial spoilage of beer.
In the current study, complete genome sequences of six isolated L. brevis strains were determined. Five of these L. brevis strains were isolated from beer (three isolates) or the brewing environment (two isolates), and were characterized as beer-spoilers or non-beer spoilers, respectively, while the sixth isolate had previously been isolated from silage. The genomic features of 19 L. brevis strains, encompassing the six L. brevis strains described in this study and thirteen L. brevis strains for which complete genome sequences were available in public databases, were analyzed with particular attention to evolutionary aspects and adaptation to beer.
Comparative genomic analysis highlighted evolution of the taxon allowing niche colonization, notably adaptation to the beer environment, with approximately 50 chromosomal genes acquired by L. brevis beer-spoiler strains representing approximately 2% of their total chromosomal genetic content. These genes primarily encode proteins that are putatively involved in oxidation-reduction reactions, transcription regulation or membrane transport, functions that may be crucial to survive the harsh conditions associated with beer. The study emphasized the role of plasmids in beer spoilage with a number of unique genes identified among L. brevis beer-spoiler strains.
Lactobacillus brevis is a member of the lactic acid bacteria (LAB), which are catalase-negative, non-sporulating, non-motile, rod or coccus-shaped Gram-positive bacteria. L. brevis grows optimally at 30 °C and within a pH range of 4 to 6 [1,2,3]. It is an obligatory hetero-fermentative bacterium producing lactic acid, carbon dioxide and ethanol and/or acetic acid [1,2,3]. Using phylogenomic and comparative genomic analysis Duar et al. studied the relatedness within the Lactobacillus genus in light of their natural habitat in order to understand their evolutionary history . They assigned lactobacilli species into three main lifestyle categories: free living (environmental and plant isolates), host adapted or as “nomadic” . Sequenced genomes of the Lactobacillus genus range in size from 1.27 (L. iners) to 4.91 (L. parakefiri) Mbp .
L. brevis has been isolated from silage, as well as from fermented cabbage and other fermented foods [5, 6], and is assigned to the free-living lifestyle group of lactobacilli . L. brevis strains, among other lactobacilli, are of particular interest as they have been granted Qualified Presumption of Safety (QPS) status and consequently have been widely used in the production of fermented foods [1, 7]. In addition to their application in food fermentations they are purported to have potential as health-promoting or probiotic bacteria [1, 7]. In contrast to these positive attributes, L. brevis strains have also been reported as the causative agent of food or beverage spoilage, in particular of beer [8, 9]. LAB species are reported to cause approximately 70% of microbial beer-spoilage incidents, and among this group L. brevis isolates are particularly problematic [10,11,12]. They are associated with the production of malodorous compounds, acidity and/or turbidity with negative impacts on the organoleptic properties of the final product. Bacterial spoilage of beer may result in product withdrawal or recall with concomitant economic losses for the brewing industry [10,11,12]. Beer spoilage by Gram-positive bacteria has been studied previously and the main mechanism of hop resistance known so far involves an active extrusion of the toxic compound using transporters identified as: (a) HorA which functions as an ABC-type multidrug transporter to expel hop compounds, in particular iso-α-acids, from the bacterial cytoplasm, (b) HorC a proton motive force-dependent hop excretion transporter, and (c) the H+-ATPase which increases the pumping of protons released from the hop compounds [13,14,15]. The transmembrane protein HitA is also thought to play a role in the transport of divalent cations, where iso-α-acids exchange protons for cellular divalent cations such as Mn2+ .
To date a number of comparative genome studies of the Lactobacillus genus have been described [1, 17,18,19], some of which have provided insights into the taxonomy of the Lactobacillus genus [3, 7], or its fermentation capabilities . Carbohydrate metabolism has been assessed in several Lactobacillus species LAB such as L. casei or L. plantarum . However, a broad comparative genome analysis of the L. brevis species has as yet not been undertaken. Recent advances in next generation sequencing technologies has facilitated a rapid surge in the number of bacterial genomes now available for comparative analysis within a genus or a species.
In the current study, Single-Molecule-Real-Time (SMRT) sequencing technology [20, 21] was employed to generate the complete genome sequence of an additional six L. brevis strains isolated from silage and the brewery environment. Using the dataset of 19 complete chromosomal sequences, a comparative genome analysis of the L. brevis taxon was undertaken through an assessment of the phylogeny, pan- and core-genome, and niche adaptation with particular emphasis on adaptation to the brewing environment. The importance of plasmids was also investigated in relation to beer spoilage ability.
Results and discussion
Isolation of L. brevis strains
Six L. brevis strains were isolated and included as part of the study, with the aim of expanding the collection of L. brevis genome sequences currently available, as well as studying their ability to grow and colonize harsh environments such as beer. Three L. brevis strains were isolated from beer and characterized as beer-spoilers based on their ability to grow in beer. Two additional L. brevis strains were isolated from the brewing environment, yet lack the ability to grow in beer and are thus not classified as beer-spoilers (Table 1 and Fig. 1). In addition, the sixth L. brevis strain sequenced as part of this study originates from silage  and was included as a non-brewing environmental isolate (Table 1). The different L. brevis isolates showed different colony morphology ranging from a dry irregular colony type for L. brevis UCCLBBS449 and UCCLB95 to a slimy and ropy colony type for UCCLBBS124 (Table 1). Plasmid profiling of the different isolates revealed a distinct plasmid content for each isolate. Growth curves in MRS broth demonstrated the unique growth profiles of the individual isolates confirming that the isolates were distinct from each other. Furthermore, only L. brevis UCCLB95, UCCLBBS124 and UCCLBBS449 were characterized as beer-spoilers having the ability to survive and grow in beer, while L. brevis SA-C12, UCCLB521 and UCCLB556 were defined as non-beer-spoilers (Table 1 and Fig. 1).
General genome features
The complete chromosomal sequences of nineteen L. brevis strains were selected for analysis, thirteen of which available at that time were obtained from the NCBI database, while the remaining six were sequenced as part of this study using SMRT sequencing technology (Table 2). These 19 selected L. brevis strains had been isolated from different ecological niches: silage, fermented food, animal’s gut and the brewery environment (Table 2). The general features of the 19 L. brevis genomes are indicated in Table 3 and include an average chromosome length of 2.49 Mbp (ranging from 2.27 to 2.79 Mbp) and a G + C content of 46%. An average of 2338 predicted CDSs (Coding Sequences) per chromosome were identified to which approximately 78.3% could be assigned a function based on in silico predictions using BLAST (Basic Alignment Search Tool), while the remaining 21.7% were annotated as hypothetical proteins (Table 3). A type II CRISPR-Cas (Clustered Regularly Interspaced Short Palindromic Repeats) locus was found in the chromosome of L. brevis BGP6, L. brevis NPS-QW145 and L. brevis SRCM101106 where variability was observed in the spacer region, distinct spacers were observed in each of these three L. brevis strains suggesting an active system acquiring unique and various spacers for protection against invading DNA over time. Conversely, in the chromosome of L. brevis TMW1.2112 and L. brevis TMW1.2113 ten identical spacers were detected suggesting that these two strains are clonal or that this CRISPR-Cas system is inactive, and that these common spacers originate from a common ancestor that acquired genetic material from viruses/plasmids that it encountered in the past . The L. brevis strain ZLB004 chromosome revealed the presence of four CRISPR locus, one was associated to a type I-E CRISPR-Cas system, a second one was associated to a type II CRISPR-Cas system both potentially active systems. The two other CRISPR locus were not associated to any CRISPR-Cas proteins suggesting inactive system.
PacBio SMRT sequencing was used to determine the diversity and frequency of methylated motifs recognized by Restriction Modification (R/M) systems within the six newly sequenced and annotated L. brevis strains as part of this study. R/M systems constitute one of the biological barriers exerted by a strain against foreign DNA . This analysis revealed the presence of various m6A motifs and allowed the identification of three motifs assignable to Type I R/M system and six motifs assignable to Type II R/M system (Table 4). The presence of specific methylated motifs was linked to the presence of specific R/M systems in the corresponding L. brevis strains (Table 4). Somewhat surprisingly, L. brevis UCCLB95 does not appear to encode any R/M systems.
The predicted mobilome of L. brevis
All complete chromosome sequences were analyzed for the presence of mobile elements such as IS elements and genes specified as encoding transposases. This analysis indicates that L. brevis strains UCCLBBS449 and UCCLB95 contain the highest number of insertion sequence (IS) elements/transposases, 114 and 132, respectively (Table 3). The genome sequences were also investigated for prophages, revealing various predicted intact or partial prophage regions (Table 3), displaying in most cases similarity to the published L. brevis temperate bacteriophage LBR48 . The plasmid content of the L. brevis strains is detailed below.
The phylogenetic relationship between the genomes of the nineteen L. brevis strains were investigated by a comparative analysis of their 16S rRNA sequences (Fig. 2a). The resulting phylogenetic tree distinguishes five clades (clades A through to E). Clade A represents two L. brevis beer-spoiling strains UCCLB95 and UCCLBBS449 both isolated from spoiled beer, displaying slow growth in nutritive media or beer (.
Table 1 and Fig. 1). Clade B encompasses three L. brevis strains: L. brevis SRCM101106 isolated from food L. brevis BDGP6 isolated from the gut of a drosophila and L. brevis NCTC13768 from an unknown isolation source. Clade C is represented by six L. brevis strains, of which one was isolated from food (L. brevis SRCM101174), three from silage (L. brevis SA-C12, ATCC 367 and 100D8) and two strains (L. brevis UCCLB521 and UCCLB556), both isolated from the brewing environment, yet unable to survive and grow in beer (Table 1). These latter two strains may have been introduced into the brewery through raw materials such as cereal grains thus explaining the observed phylogenetic relation to the silage L. brevis isolate SA-C12. Clade D includes five L. brevis strains, all retrieved as beer-spoiler strains from the brewing environment, and all exhibiting a slimy, ropy phenotype (Table 1) . Clade E gathers three L. brevis strains, two isolated from fermented food (L. brevis KB290 and NPS-QW-145) as well as L. brevis ZLB004 isolated from pig’s feces.
In order to obtain a more refined view of the phylogeny of the 19 analyzed strains, a so-called phylogenetic supertree was constructed based on 631 conserved orthologous proteins that had been identified as single-copy genes conserved across all investigated chromosomal sequences (19 L. brevis strains and Enterococcus faecalis V583 as an outgroup) [26, 27]. This supertree does not display distinct clades separating the L. brevis strains in different groups as was observed with the 16S rRNA phylogenetic tree, suggesting a close relatedness within the species (Fig. 2b). However, upon close inspection of this phylogenetic tree, it appears that L. brevis strains isolated from food and silage cluster on one branch of the tree, while L. brevis brewery isolates cluster on another. The L. brevis strains isolated from gut microbiota BDGP6 and ZLB004 as well as L. brevis strains NCTC13768 and 100D8 gather in a smaller clade.
A pan-genome analysis was performed in order to determine the total number of distinct genes present on the combined chromosomal sequences of the analyzed L. brevis strains. The pan-genome curve displays an asymptotic trend, growing with an average rate of 136 genes per genome in the first nine iterations, then the number of new genes decreased leading to a total pan-genome content of 3968 genes (Fig. 3). The mathematical function displayed on the graph reveals an exponential value lower than 0.5 indicating that the pan-genome is in a closed state. The core genome was determined to encompass 1428 genes (Fig. 3). As an overall result, both analyses indicate a closed pan-genome for L. brevis species, while also indicating that a sufficient number of strains had been included to adequately describe the genetic repertoire of the L. brevis species.
Comparative analysis of orthologous genes
The comparative analysis used in this study was based on chromosomal sequences. The core genome of 1428 genes is divided in 1170 orthologous gene families (single copy) and 258 paralogous gene families (multi-copy). Unique gene families to each chromosome were also recorded and 246 unique gene families were identified across the nineteen L. brevis strains (Fig. 4a). Functional assignment efforts revealed that 75.2% of the unique gene families encoded proteins of unknown function (hypothetical proteins), while 4.5% encoded (pro)phage-related proteins. The remaining unique gene families encode proteins that could benefit the fitness of the strain such as CRISPR-Cas system (e.g. Type I-E CRISPR Cas system in L. brevis ZLB004), restriction-modification systems (e.g. Type I R/M system in L. brevis strain UCCLBBS124), or cell wall polysaccharide synthesis (e.g. genes predicted to encode glycosyltransferases and a polysaccharide polymerase only found in the L. brevis strain NPS-QW-145).
In order to further investigate the functionality and diversity encoded by the core and dispensable genomes, a Cluster of Orthologous Group (COG) analysis was employed. The genome content of the 19 selected L. brevis strains was classified into different groups depending on their function. More than 75% were predicted to be involved in housekeeping functions, vital for the strain to grow such as those participating in transcription or translation. Approximately 16% of the genes were assigned to COG groups with only a general function predicted or of unknown function (Fig. 4b).
Evolution and adaptation to beer environment
When plotting the number of CDSs as a function of genome size for the different L. brevis strains, the group exhibiting the largest genome size as well as the highest number of CDS are L. brevis strains isolated from beer and characterized as beer-spoilers as well as the L. brevis strain BDGP6 displaying the biggest CDS number. L. brevis strains known to be beer-spoilers possess an average of 2385 CDS, while those isolated from food, silage, animal’s gut and non-beer spoiling brewery isolates display an average of 2311 CDSs (Fig. 5). This observation suggests a link to adaptation to a new environment, i.e. the beer or brewery environment, which may have necessitated the acquisition of novel genes and corresponding functions in order to survive in this harsh environment. To understand if the beer-spoiling strains had acquired a specific set of genes or associated functions, genes that may putatively be associated to beer adaptation were first predicted to be those that would be present in the genomes of at least four beer spoiling strains (Table 5). From this analysis, 58 genes of interest were highlighted as well as 26 hypothetical proteins. Out of these 58 genes, approximately 21% encode proteins related to oxido-reduction reactions (Flavodoxin, oxidoreductases and short-chain dehydrogenases), 22% are linked to transcription (transcriptional regulators, RNA polymerase sigma-24 subunit ECF subfamily), 21% encode membrane and cell surface proteins and 14% are related to membrane transport (MFS transporter, permease, ABC transporters) (Table 5).
When exposed to beer, bacteria are subjected to various stresses, among them a low pH (3.8–4.7) and hop compounds . When iso-α-acids enter the cell cytoplasm, they dissociate into hop anions and protons decreasing the intracellular pH . Therefore, bacteria would have to adapt in order to regulate their internal pH in order to survive [28,29,30]. Furthermore, in beer the presence of ethanol (0.5–10% w/w) causes oxidative stress in bacteria, this results in the production of Reactive Oxygen Species (ROS) such as hydrogen peroxide and free radicals leading to cell damage [31, 32]. Despite the stress and harsh environment imposed by the beer environment, some bacteria have evolved and acclimatized to this new medium. It may thus be possible that some of the L. brevis strains acquired additional functions which allow them to grow and survive in beer and which has led to an increased genome size. The fact that 21% of these chromosomal genes encode proteins related to redox reactions is of interest and suggests a link between L. brevis beer-spoiler strains and oxidative stress response. Six of the 12 genes that encode functions relating to oxido-reduction reactions present in at least four beer-spoiler L. brevis strains are predicted to encode NADH oxidoreductases and short-chain dehydrogenases/reductases (SDRs). These proteins are part of the large family of NAD(P)(H)-dependent oxidoreductases and are believed to behave as scaffold proteins for an NAD(P)(H) redox sensor system . In previous studies, the role of SDRs during oxidative stress was highlighted in species such as Bacillus subtilis where they are required for survival in severe ethanol stress , or in Burkholderia pseudomallei during salt stress .
Furthermore, 22% of the chromosomal genes that seem to be specifically associated with beer-spoiling L. brevis strains are linked to transcriptional regulation, suggesting that these regulators act on specific genes to control their expression and confer an advantage when present in beer. It would be interesting to study which genes are affected by these transcriptional regulators to assess the mechanisms employed to survive in this harsh environment. Of the L. brevis beer-spoiler specific chromosomal genes 21% encode membrane and cell surface proteins suggesting an adaptation to survive in the harsh beer environment. 14% encode proteins associated with membrane transport such as permeases and ABC transporters suggesting exchange between the strain and its environment and possibly a role in extrusion where the L. brevis isolate would expel protons or iso-α-acids in order to survive and thrive in beer, as has been described previously [10, 13, 15].
Interestingly some of the chromosomal genes identified among L. brevis beer-spoiler strains in this analysis had also been highlighted in a previous study as beer-spoilage diagnostic marker genes (DMG) . These genes are predicted to code for an oligogalacturonide transporter, a short chain dehydrogenase and a RNA polymerase sigma factor ECF subfamily, which reinforces the hypothesis for their involvement in beer spoilage adaptation.
The role of plasmids in adaptation to beer environment
Different proteins involved in beer spoilage have been identified on plasmids indicating the importance of plasmids for bacterial strains in beer spoilage. This might suggest a role for plasmid mobilization and transfer between bacterial strains throughout evolution to adapt to a new environment such as beer.
The nineteen analyzed L. brevis strains were predicted to harbour up to nine plasmids with strains L. brevis KB290 and L. brevis UCCLBBS449 exhibiting the largest plasmid complements of the assessed strains. The plasmid size ranges from 2.8 Kb to 107.0 Kb (Table 3). The number of plasmids and their size do not appear to be linked to the isolation source of the L. brevis strains (e.g. four plasmids for L. brevis SRCM101106 versus nine plasmids for L. brevis KB290, both isolated from fermented food) or to the beer spoilage ability of the isolate (two plasmids for L. brevis UCCLB95 versus nine plasmids for L. brevis UCCLBBS449 both characterised as beer-spoilers). Investigating analogies between plasmids among L. brevis beer-spoiler strains revealed that the plasmid content of L. brevis TMW1.2108 and L. brevis TMW1.2111 were very similar. Indeed, the six plasmids of L. brevis TMW1.2111 show at least 90% identity to seven of the eight plasmids contained by strain L. brevis TMW1.2108, with the exception of plasmid TMW1.2108–5. Similarly, L. brevis strains TMW1.2112 and TMW1.2113 present a close plasmid composition as the four plasmids of L. brevis TMW1.2113 are at least 90% identical to four out of five plasmids of L. brevis TMW1.2112 with the exception of plasmid TMW1.2112–1.
Out of the 38 plasmids shared between L. brevis beer-spoiler strains, only three plasmids seem to be unique, sharing less than 10% similarity with any other plasmid. These three plasmids were found in L. brevis UCCLBBS449 (UCCLBBS449_pF, UCCLBBS449_pH and UCCLBBS449_pI) and contain mostly genes coding for hypothetical proteins, replication proteins as well as genes coding for proteins involved in conjugation such as mobilization proteins and a relaxase.
Refined analysis of specific genes shared only between at least three L. brevis beer-spoiler strains, identified only twenty-five genes (Table 6). In this list of unique genes shared only between L. brevis beer-spoiler strains, the gene coding for the membrane protein HorC is noteworthy, as it is known to be involved in beer spoilage  and is present in all L. brevis beer-spoiler strains with the exception of L. brevis TMW1.2113.
Interestingly, the gene encoding the ABC transporter HorA  and present in L. brevis beer-spoiler strains TMW1.2108, TMW1.2111, TMW1.2113, UCCLBBS124 and UCCLBBS449 does not figure in this list as a similar protein can be found in plasmid sequences of the L. brevis strains KB290, SRCM101106 isolated from fermented food and L. brevis UCCLB556 isolated from the brewery and characterized as a non-beer spoiler. Moreover, the transmembrane protein HitA  has been identified only in two of the L. brevis beer-spoiler strains UCCLBBS449 and TMW1.2112. These observations reinforce the statement that involvement of these genes in beer survival and spoilage is not always verified as they are not consistently present in beer-spoiler organisms nor are always corresponding to beer spoilage ability if present in a strain . The list of genes present only in L. brevis beer-spoiler strains shows that strains L. brevis TMW1.2108 and TMW1.2111 possess more than 90% of these genes whereas L. brevis UCCLB95 only possesses one gene coding for the membrane transporter HorC (Table 6). The remainder of the L. brevis beer-spoiler strains carry approximately 50% of these particular genes (Table 6).
Out of these 25 unique genes shared among L. brevis beer-spoiler strains approximately 25% code for hypothetical proteins of unknown function. Meanwhile, ~ 30% of these genes appear to encode cell wall-associated proteins either as membrane transporters (MFS transporter, HorC) or as cell wall biosynthesis (lipopolysaccharide glycosyltransferases, acyltransferases). As mentioned above, a beer-spoiling strain would need to extrude toxic compounds using transporters and adapt its cell wall composition to survive the harsh beer environment. A smaller portion of these unique genes are linked to transcription regulation, replication or mobilome.
Interestingly, some of the plasmid-associated genes identified among L. brevis beer-spoiler strains in this analysis have also been highlighted previously as unique attributes of beer-spoiling strain plasmids . The gene coding for the CrcB like-protein involved in ion transport was found on plasmid BSO 464–2 of the L. brevis beer-spoiler strain BSO 464 as well as a gene coding for enolase involved in glucose metabolism. A gene coding for cytosine deaminase is present in five out of the seven L. brevis beer-spoiler strains used in this analysis (Table 6) and was identified as a unique attribute on the plasmid pPECL-8 of the beer-spoiler Pediococcus claussenii ATCC BAA-344 . Moreover, a plasmid-associated gene coding for a glycosyltransferase family 2 was highlighted in the analysis (Table 6), this protein was associated with excess β-glucan formation leading to a slimy ropy phenotype in the L. brevis beer-spoiler TMW1.2112 . This slimy phenotype was observed in the L. brevis strain UCCLBBS124 first described in this study (Table 1), and the gene coding for the glycosyltransferase family 2 was identified on one of its plasmid UCCLBBS124_pB.
This overall examination of plasmid-associated genes shows the importance of extrachromosomal DNA in beer spoilage adaptation and opens new possibilities for understanding the beer spoilage process with an updated list of potential proteins of interest only present in L. brevis beer-spoiler strains.
The isolation and genome sequencing of six L. brevis strains combined with thirteen additional, publicly available L. brevis genomes allowed a comparative genome analysis of the L. brevis species. The deduced pan-genome of these L. brevis isolates appears to be in a closed state, indicating that the representatives used in this study are sufficient to describe the genetic diversity of the taxon. Throughout evolution, it appears that L. brevis strains specified and differentiated one from another by acquiring plasmids and prophages, despite for the presence of CRISPR-Cas and R/M systems which may have limited such foreign DNA invasion events. These latter systems are of relevance for future functional investigations that may necessitate the development of DNA transfer and/or mutagenesis tools. L. brevis strains represent a significant threat for the brewing industry being the most common cause of beer spoilage; however, this spoiling ability is strain specific. The comparative genome analysis performed here highlights that L. brevis strains with the ability to grow in beer possess a higher number of CDSs in their overall chromosomal sequences. This observation suggests a link to evolution and adaptation to beer in which the strain would have acquired novel genes and functions in order to adapt and survive in the harsh environment that beer represents. The role(s) of the “acquired” or beer-specific CDSs revealed that almost a quarter of these are linked to oxido-reduction reactions, possibly playing a role in the response to oxidative stress. Another 22% are linked to transcription regulation, 21% encode cell surface proteins while 14% are encoding membrane transport related proteins and possibly associated to harmful compound extrusion encountered by the L. brevis strains when surviving and growing in beer. Additional genetic diversification of these L. brevis strains is expected to have occurred through plasmid acquisition that also likely contributes to beer adaptation. The plasmid content analysis of the different L. brevis beer-spoiler strains highlighted the presence of unique proteins shared among these strains. These proteins are mostly hypothetical proteins while approximately 30% are linked to membrane transport, and cell-wall synthesis. These observations demonstrate the complexity of microorganisms’ beer spoilage ability and suggests that adaptation of the L. brevis strain to beer is a complex process, not due to the action of only one specific gene, but more likely the intervention of a complex, multi-factorial response.
Isolation of L. brevis strains
Five distinct L. brevis strains (UCCLB521, UCCLB556, UCCLB95, UCCLBBS124 and UCCLBBS449) were isolated from the brewing environment, while SA-C12 had previously been isolated from silage . The strains were characterized by evaluating their plasmid content as well as growth curve profiles in MRS broth or in beer at 30 °C. Plasmids were isolated after overnight growth of the L. brevis strains in MRS broth at 30 °C, cells were pelleted by centrifugation for 10 min at 5000 rpm followed by lysozyme treatment (30 mg/mL lysozyme in TE + 25% sucrose) at 37 °C for 30 min. Plasmids were extracted using the GeneJET Plasmid Miniprep Kit (Thermo Scientific™). Plasmid profiles of the different L. brevis isolates were analyzed using a 1% agarose gel. Growth curve profiles in MRS broth or in beer were performed at 30 °C by hourly OD600nm measurements for a period of 55 h. Moreover, colony morphology was recorded following growth on MRS agar plate at 30 °C.
Sequencing and annotation
L. brevis strains were streaked on MRS agar plates and grown at 30 °C for 24 h. For each L. brevis strain, a single colony was inoculated into MRS broth and grown overnight at 30 °C. Cells were pelleted by centrifugation at 5000 rpm for 10 min. The supernatant was removed and the pelleted cells were frozen at − 20 °C prior sending for sequencing. Sequencing was performed using the PacBio SMRT next generation sequencing technology (performed by GATC Biotech, Germany). De novo genome assemblies were performed using the Pacific Biosciences SMRT Portal analysis platform. Open Reading Frame (ORF) or coding sequence (CDS) prediction was performed using Prodigal prediction software  and confirmed using BLASTX alignments . Automatic annotations were refined using Artemis v16.0.0 where ORF predictions were manually checked, start codons adjusted and pseudogenes identified. Transfer RNA (tRNA) genes were predicted using tRNA-scan-SE v2.0 , while ribosomal RNA (rRNA) genes were identified using RNAmmer v1.2 . The sixteen L. brevis genomes obtained from NCBI were re-annotated as described above in order to treat identically all sequenced genomes used in this study.
Following de novo genome assembly, the RS_Modification_and_Motif_Analysis.1 protocol of the SMRT Analysis portal was employed for base modification and methylated motif detection. This analysis was performed on L. brevis strains sequenced, assembled and annotated as part of this study. Methylation motifs with a score equal or higher than 40 (corresponding to a P-value of < 0.0001) were considered specific and used for further analysis. ORFs of genomes were investigated for the presence of restriction/modification systems using the BLASTP alignment function of the REBASE database  (cut-off E-value of 0.0001; with at least 30% similarity over at least 80% of the sequence length). A comparative genome analysis was employed to associate the presence of R/M system-encoding genes with the presence of methylation motif(s).
All protein sequence comparisons were performed using all-against-all, bi-directional BLAST alignments . An alignment cut-off value of E-value 0.0001, and a similarity cut-off level of at least 30% amino acid identity across 80% of the sequence length was used. Results were analyzed with the Markov Clustering Algorithm (MCL)  and proteins encoded were categorized in predicted functional groups based on COG (Clusters of Orthologous Groups) assignments .
The supertree was prepared using the BLAST-based comparative approach described above in order to identify chromosomal orthologous proteins. The set of chromosomal orthologous proteins was concatenated for each strain and an ungapped alignment was performed using MUSCLE v3.8.31 . The phylogenetic tree was computed using the maximum-likelihood method in PhyML v3.0 and bootstrapped employing 1000 replicates . The final tree was visualized using MEGA7. A tree based on 16 S rRNA genes was constructed using clustalw and visualized via ITOL (Interactive Tree Of Life) . The chromosome sequence of Enterococcus faecalis V583 (Accession: AE016830) was included as an outgroup.
The pan-core genome analysis of the above-mentioned 19 L. brevis chromosomal sequences, was performed using PGAP v1.0 . ORF content for each chromosome is classified in functional gene clusters using the Gene Family method. From this analysis a pan/core genome profile was generated.
Genome accession numbers
L. brevis 100D8: CP015338, L. brevis ATCC 367: CP000416, L. brevis BDGP6: CP024635, L. brevis KB290: AP012167, L. brevis NCTC13768: LS483405, L. brevis NPS-QW-145: CP015398, L. brevis SA-C12: CP031185, L. brevis SA-C12_pA: CP031186, L. brevis SA-C12_pB: CP031187, L. brevis SRCM101106: CP021674, L. brevis SRCM101174: CP021479, L. brevis TMW 1.2108: CP019734, L. brevis TMW 1.2111: CP019743, L. brevis TMW 1.2112: CP016797, L. brevis TMW 1.2113: CP019750, L. brevis UCCLB521: CP031208, L. brevis UCCLB521_pA: CP031209, L. brevis UCCLB521_pB: CP031210, L. brevis UCCLB521_pC: CP031211, L. brevis UCCLB521_pD: CP031212, L. brevis UCCLB521_pE: CP031213, L. brevis UCCLB556: CP031174, L. brevis UCCLB556_pA: CP031175, L. brevis UCCLB556_pB: CP031176, L. brevis UCCLB556_pC: CP031177, L. brevis UCCLB556_pD: CP031178, L. brevis UCCLB556_pE: CP031179, L. brevis UCCLB556_pF: CP031180, L. brevis UCCLB556_pG: CP031181, L. brevis UCCLB95: CP031182, L. brevis UCCLB95_pA: CP031183, L. brevis UCCLB95_pB: CP031184, L. brevis UCCLBBS124: CP031169, L. brevis UCCLBBS124_pA: CP031170, L. brevis UCCLBBS124_pB: CP031171, L. brevis UCCLBBS124_pC: CP031172, L. brevis UCCLBBS124_pD: CP031173, L. brevis UCCLBBS449: CP031198, L. brevis UCCLBBS449_pA: CP031199, L. brevis UCCLBBS449_pB: CP031200, L. brevis UCCLBBS449_pC: CP031201, L. brevis UCCLBBS449_pD: CP031202, L. brevis UCCLBBS449_pE: CP031203, L. brevis UCCLBBS449_pF: CP031204, L. brevis UCCLBBS449_pG: CP031205, L. brevis UCCLBBS449_pH: CP031206, L. brevis UCCLBBS449_pI: CP031207, L. brevis ZLB004: CP021456 and Enterococcus faecalis V583: AE016830.
Basic Alignment Search Tool
Cluster of Orthologous Group
Clustered Regularly Interspaced Short Palindromic Repeats
Diagnostic Marker Genes
Interactive Tree Of Life
Lactic Acid Bacteria
Markov Clustering Algorithm
Open Reading Frame
Qualified Presumption of Safety
Reactive Oxygen Species
Vos, P., G. Garrity, D. Jones, N.R. Krieg, W. Ludwig, F.A. Rainey, K.-H. Schleifer, and W. Whitman, Bergey's manual of systematic bacteriology: volume 3: the Firmicutes. 2011. 3: p. 464–511.
Kandler O. Carbohydrate metabolism in lactic acid bacteria. Antonie Van Leeuwenhoek. 1983;49(3):209–24.
Felis GE, Dellaglio F. Taxonomy of lactobacilli and Bifidobacteria. Current issues in intestinal microbiology. 2007;8(2):44–61.
Duar RM, Lin XB, Zheng J, Martino ME, Grenier T, Pérez-Muñoz ME, Leulier F, Gänzle M, Walter J. Lifestyles in transition: evolution and natural history of the genus Lactobacillus. FEMS Microbiol Rev. 2017;41(Supp_1):S27–48.
Makarova K, Slesarev A, Wolf Y, Sorokin A, Mirkin B, Koonin E, Pavlov A, Pavlova N, Karamychev V, Polouchine N, Shakhova V, Grigoriev I, Lou Y, Rohksar D, Lucas S, Huang K, Goodstein DM, Hawkins T, Plengvidhya V, Welker D, Hughes J, Goh Y, Benson A, Baldwin K, Lee J, Díaz-Muñiz I, Dosti B, Smeianov V, Wechter W, Barabote R, Lorca G, Altermann E, Barrangou R, Ganesan B, Xie Y, Rawsthorne H, Tamir D, Parker C, Breidt F, Broadbent J, Hutkins R, O'Sullivan D, Steele J, Unlu G, Saier M, Klaenhammer T, Richardson P, Kozyavkin S, Weimer B, Mills D. Comparative genomics of the lactic acid bacteria. Proc Natl Acad Sci U S A. 2006;103(42):15611–6.
Fukao M, Oshima K, Morita H, Toh H, Suda W, Kim SW, Suzuki S, Yakabe T, Hattori M, Yajima N. Genomic analysis by deep sequencing of the probiotic Lactobacillus brevis KB290 harboring nine plasmids reveals genomic stability. PLoS One. 2013;8(3).
Salvetti E, Torriani S, Felis GE. The genus Lactobacillus: a taxonomic update. Probiotics and antimicrobial proteins. 2012;4(4):217–26.
Bergsveinson, J., V. Pittet, E. Ewen, N. Baecker, and B. Ziola, Genome Sequence of Rapid Beer-Spoiling Isolate Lactobacillus brevis BSO 464. 2015. 3(6).
Fraunhofer ME, Geissler AJ, Wefers D, Bunzel M, Jakob F, Vogel RF. Characterization of β-glucan formation by Lactobacillus brevis TMW 1.2112 isolated from slimy spoiled beer. Int J Biol Macromol. 2017.
Sakamoto K, Konings WN. Beer spoilage bacteria and hop resistance. Int J Food Microbiol. 2003;89(2–3):105–24.
Vaughan A, O'Sullivan T, Van Sinderen D. Enhancing the microbiological stability of malt and beer — a review. J Inst Brew. 2005;111(4):355–71.
Bokulich NA, Bamforth CW. The microbiology of malting and brewing. Microbiology and molecular biology reviews : MMBR. 2013;77(2):157–72.
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(18):5371–5.
Iijima K, Suzuki K, Ozaki K, Yamashita H. horC confers beer-spoilage ability on hop-sensitive Lactobacillus brevis ABBC45cc. J Appl Microbiol. 2006;100(6):1282–8.
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(2):173–91.
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(5):596–603.
Kant R, Blom J, Palva A, Siezen RJ, de Vos WM. Comparative genomics of Lactobacillus. Microb Biotechnol. 2011;4(3):323–32.
Claesson MJ, van Sinderen D, O'Toole PW. Lactobacillus phylogenomics--towards a reclassification of the genus. Int J Syst Evol Microbiol. 2008;58(Pt 12:2945–54.
Goh YJ, Klaenhammer TR. Genomic features of Lactobacillus species. Front Biosci (Landmark Ed). 2009;(14):1362–86.
Gupta PK. Single-molecule DNA sequencing technologies for future genomics research. Trends Biotechnol. 2008;26(11):602–11.
McCarthy A. Third generation DNA sequencing: Pacific Biosciences' single molecule real time technology. Chem Biol. 2010;17(7):675–6.
Deasy T, Mahony J, Neve H, Heller KJ, van Sinderen D. Isolation of a virulent Lactobacillus brevis phage and its application in the control of beer spoilage. J Food Prot. 2011;74(12):2157–61.
Horvath P, Barrangou R. CRISPR/Cas, the immune system of Bacteria and archaea. Science. 2010;327(5962):167–70.
Modrich P. Structures and mechanisms of DNA restriction and modification enzymes. Q Rev Biophys. 2009;12(3):315–69.
Jang SH, Yoon BH, Chang HI. Complete nucleotide sequence of the temperate bacteriophage LBR48, a new member of the family Myoviridae. Arch Virol. 2011;156(2):319–22.
Bottacini, F., M. O’Connell-Motherway, J. Kuczynski, K.J. O’Connell, F. Serafini, S. Duranti, C. Milani, F. Turroni, G.A. Lugli, A. Zomer, D. Zhurina, C. Riedel, M. Ventura, and D. Sinderen, Comparative genomics of the Bifidobacterium breve taxon. BMC Genomics, 2014. 15.
Kelleher P, Bottacini F, Mahony J, Kilcawley KN, van Sinderen D. Comparative and functional genomics of the Lactococcus lactis taxon; insights into evolution and niche adaptation. BMC Genomics. 2017;18(1):267.
Booth IR. Regulation of cytoplasmic pH in bacteria. Microbiol Rev. 1985;49(4):359–78.
Cook GM, Russell JB. The effect of extracellular pH and lactic acid on pH homeostasis inLactococcus lactis andStreptococcus bovis. Curr Microbiol. 1994;28(3):165–8.
Siegumfeldt H, Rechinger KB, Jakobsen M. Dynamic changes of intracellular pH in individual lactic acid bacterium cells in response to a rapid drop in extracellular pH. Appl Environ Microbiol. 2000;66(6):2330–5.
Betteridge DJ. What is oxidative stress? Metabolism. 2000;49(2 Suppl 1):3–8.
van de Guchte M, Serror P, Chervaux C, Smokvina T, Ehrlich SD, Maguin E. Stress responses in lactic acid bacteria. Antonie Van Leeuwenhoek. 2002;82(1):187–216.
Kavanagh KL, Jörnvall H, Persson B, Oppermann U. Medium- and short-chain dehydrogenase/reductase gene and protein families: the SDR superfamily: functional and structural diversity within a family of metabolic and regulatory enzymes. Cell Mol Life Sci. 2008;65(24):3895–906.
Höper D, Völker U, Hecker M. Comprehensive characterization of the contribution of individual SigB-dependent general stress genes to stress resistance of Bacillus subtilis. J Bacteriol. 2005;187(8):2810–26.
Pumirat P, Boonyuen U, Vanaporn M, Pinweha P, Tandhavanant S, Korbsrisate S, Chantratita N. The role of short-chain dehydrogenase/oxidoreductase, induced by salt stress, on host interaction of B. pseudomallei. BMC Microbiol. 2014;14(1).
Bergsveinson J, Ziola B. Comparative genomic and plasmid analysis of beer-spoiling and non-beer-spoiling Lactobacillus brevis isolates. Can J Microbiol. 2017;63(12):970–83.
Fraunhofer ME, Geissler AJ, Wefers D, Bunzel M, Jakob F, Vogel RF. Characterization of β-glucan formation by Lactobacillus brevis TMW 1.2112 isolated from slimy spoiled beer. Int J Biol Macromol. 2018;107:874–81.
Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC bioinformatics. 2010;11(1):119.
Altschul S, Gish W, Miller W, Myers E, Lipman D. Basic local alignment search tool. J Mol Biol. 1990;215.
Schattner P, Brooks AN, Lowe TM. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 2005;33(suppl_2):W686–9.
Lagesen K, Hallin P, Rødland EA, Stærfeldt H-H, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35(9):3100–8.
Roberts RJ, Vincze T, Posfai J, Macelis D. REBASE—a database for DNA restriction and modification: enzymes, genes and genomes. Nucleic Acids Res. 2015;(43, Database issue):D298–9.
Enright AJ, Van Dongen S, Ouzounis CA. An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 2002;30(7):1575–84.
Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28(1):33–6.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Guindon S, Gascuel O, Rannala B. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–5.
Zhao Y, Wu J, Yang J, Sun S, Xiao J, Yu J. PGAP: pan-genomes analysis pipeline. Bioinformatics. 2012;28.
Fukao M, Oshima K, Morita H, Toh H, Suda W, Kim S-W, Suzuki S, Yakabe T, Hattori M, Yajima N. Genomic analysis by deep sequencing of the probiotic Lactobacillus brevis KB290 harboring nine plasmids reveals genomic stability. PLoS One. 2013;8(3):e60521.
Wu Q, Law YS, Shah NP. Dairy Streptococcus thermophilus improves cell viability of Lactobacillus brevis NPS-QW-145 and its γ-aminobutyric acid biosynthesis ability in milk. Sci Rep. 2015;5.
Marine Feyereisen is the recipient of an Irish Research Council Enterprise Partnership Scheme postgraduate scholarship (Ref. No. EPSPG/2015/7). Douwe van Sinderen is supported by a Principal Investigator award (Ref. No. 450 13/IA/1953) through Science Foundation Ireland (SFI). Jennifer Mahony is in receipt of Starting Investigator Research Grant (SIRG) (Ref. No. 15/SIRG/3430) funded by Science Foundation Ireland (SFI). The funding bodies were involved in study design.
Availability of data and materials
All genome sequences analyzed during the current study are available from the National Centre for Biotechnology Information (NCBI) genome database. Additionally, all genomes sequenced during the current study have been deposited in the NCBI genome database. Accession numbers are as follow: L. brevis SA-C12: CP031185, L. brevis SA-C12_pA: CP031186, L. brevis SA-C12_pB: CP031187,, L. brevis UCCLB521: CP031208, L. brevis UCCLB521_pA: CP031209, L. brevis UCCLB521_pB: CP031210, L. brevis UCCLB521_pC: CP031211, L. brevis UCCLB521_pD: CP031212, L. brevis UCCLB521_pE: CP031213, L. brevis UCCLB556: CP031174, L. brevis UCCLB556_pA: CP031175, L. brevis UCCLB556_pB: CP031176, L. brevis UCCLB556_pC: CP031177, L. brevis UCCLB556_pD: CP031178, L. brevis UCCLB556_pE: CP031179, L. brevis UCCLB556_pF: CP031180, L. brevis UCCLB556_pG: CP031181, L. brevis UCCLB95: CP031182, L. brevis UCCLB95_pA: CP031183, L. brevis UCCLB95_pB: CP031184, L. brevis UCCLBBS124: CP031169, L. brevis UCCLBBS124_pA: CP031170, L. brevis UCCLBBS124_pB: CP031171, L. brevis UCCLBBS124_pC: CP031172, L. brevis UCCLBBS124_pD: CP031173, L. brevis UCCLBBS449: CP031198, L. brevis UCCLBBS449_pA: CP031199, L. brevis UCCLBBS449_pB: CP031200, L. brevis UCCLBBS449_pC: CP031201, L. brevis UCCLBBS449_pD: CP031202, L. brevis UCCLBBS449_pE: CP031203, L. brevis UCCLBBS449_pF: CP031204, L. brevis UCCLBBS449_pG: CP031205, L. brevis UCCLBBS449_pH: CP031206, L. brevis UCCLBBS449_pI: CP031207.
Ethics approval and consent to participate
Consent for publication
The authors declare that TOS and JG are employees of Heineken and supplied Lb. brevis strains UCCLB521, UCCLB556, UCCLB95, UCCLBBS124 and UCCLBBS449.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Feyereisen, M., Mahony, J., Kelleher, P. et al. Comparative genome analysis of the Lactobacillus brevis species. BMC Genomics 20, 416 (2019) doi:10.1186/s12864-019-5783-1
- Lactobacillus brevis
- SMRT sequencing
- Beer adaptation
- Beer spoilage