ARTICLE

Complete Genome Analysis of Probiotic Levilactobacillus brevis NS2301G3

Kiyeop Kim1https://orcid.org/0000-0002-1072-4032, Byung Wook Son2https://orcid.org/0009-0002-6104-3330, Jong Hoon Kim2https://orcid.org/0009-0006-1336-9722, Sejong Oh1,*https://orcid.org/0000-0002-5870-3038
Author Information & Copyright
1Division of Animal Science, Chonnam National University, Gwangju, Korea
2Research & Development Center, Nong Shim Co. Ltd., Seoul, Korea
*Corresponding author : Sejong Oh, Division of Animal Science, Chonnam National University, Gwangju, Korea, Tel : +82-62-530-0822, Fax : +82-62-530-2129, E-mail : soh@jnu.ac.kr

© Copyright 2024, Korean Society of Dairy Science and Biotechnology. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Dec 10, 2024; Revised: Dec 16, 2024; Accepted: Dec 17, 2024

Published Online: Dec 31, 2024

Abstract

Levilactobacillus brevis NS2301G3, a strain isolated from fermented foods, was analyzed for its complete genome sequenceto identify genes related to probiotic characteristics and functional metabolic pathways. Genome annotation revealed 2,789 protein-coding sequences, highlighting genes involved in carbohydrate metabolism, stress tolerance, and vitamin biosynthesis. Genes responsible for lactic acid production, antimicrobial peptide synthesis, and gamma-aminobutyric acid (GABA) synthesis were identified, emphasizing the potential applications of this strain in fermented food production and probiotic use. CRISPR-associated proteins, exopolysaccharide production genes, and stress response mechanisms further support its adaptability and beneficial characteristics, suggesting its potential as a valuable microbe for human health and food industry applications.

Keywords: Levilactobacillus brevis NS2301G3; lactic acid bacteria; probiotics; whole genome sequencing

Introduction

Fermented foods have long been recognized for their health benefits, with lactic acid bacteria (LAB) playing a crucial role in the fermentation process, contributing to the flavor, texture, and nutritional value of these foods [1]. Among LAB, Lactobacillus brevis is a significant species commonly found in various fermented products, such as kimchi, sauerkraut, and fermented beverages [2]. The ability of L. brevis to adapt to different environments and metabolize a wide range of substrates makes it an important microorganism for both food fermentation and potential probiotic applications [3]. Although L. brevis is not commonly used in fermented dairy products, recent studies have highlighted its potential as a starter culture in fermented milk applications due to its ability to produce gamma-aminobutyric acid (GABA), a bioactive compound with health-promoting properties, including stress reduction and blood pressure regulation [4]. Additionally, there are reports of its application in dairy fermentation processes as a complementary starter culture, contributing to both functional and sensory properties [5,6]. These studies suggest that L. brevis can play a dual role as a fermentative microorganism and a probiotic agent in dairy products. In this study, we analyzed the gene content and functional roles of L. brevis NS2301G3, a strain isolated from traditional fermented food, focusing on its contributions to carbohydrate metabolism, stress response, antimicrobial activity, and probiotic potential. The findings from this study highlight the versatility of this strain and its potential applications not only in traditional fermentation processes but also in the production of functional dairy products and other health-promoting foods.

Materials and Methods

1. Kimchi sample collection and isolation of Lactobacillus

Home-made kimchi samples were homogenized in sterile saline and plated on MRS (pH 5.0) agar [7]. Plates were then incubated anaerobically 37°C for 48–72 hours to allow for the growth of white colonies. Biochemical tests such as catalase test, Gram staining, and 16S rRNA gene analysis (Macrogen, Korea) were performed to initially identify the isolated strain. The strain identified as L. brevis, exhibiting high acid and bile tolerance, was designated as NS2301G3, and its whole genome was analyzed.

2. Hybrid sequencing

Two separate genomic DNA libraries were prepared according to the requirements of the Illumina and Oxford Nanopore systems. A combination of long-read Nanopore GridION and short-read Illumina Nextseq2000 platforms was used to generate the complete genome sequence of L. brevis NS2301G2. For Illumina sequencing, the extracted genomic DNA was fragmented by sonication using a Covaris M220 (Covaris, USA). The sheared DNA were then used to prepare a WGS library with an average insert size of 450 bp using a TruSeq Nano DNA Sample Prep kit (Illumina, USA). The library was sequenced on an Illumina Nextseq2000 platform (Illumina) using the 300 bp paired-end sequencing mode. For Nanopore sequencing, a MinION sequencing library was prepared using the Nanopore Ligation Sequencing Kit (SQK-LSK114; Oxford Nanopore, UK). The library was sequenced with an R10.4.1 GridION flow cell (Flongle) for a 24 h run using MinKNOW with the default settings (MinKNOW core 5.0.0, Guppy 6.0.6).

3. Preassembly

Illumina and Nanopore data were prepared for assembly, respectively with different options. Illumina Sequencing data were processed to remove low quality bases and adapter sequences with the optimized settings using Trimmomatic v0.39 (LEADING:10 TRAILING:10 SLIDINGWINDOW:4:20 MINLEN:200) [8]. Subsequently, additional phiX control were removed from pre-assembled data [9,10]. Trimmed sequences were aligned against phiX genome with bowtie2 v2.3.5.1 with the default options and filtered out by samtools v1.9 [11]. Nanopore sequencing data was basecalled with guppy basecaller v3.1.5. NanoFilt v2.8.0 was used to filter obtained reads with average Phred quality score lower than 7 and length lower than 1,000 [12,13].

1) Genome assembly and annotation

Unicycler v0.4.8 was used to construct genome combined with Filtered NextSeq2000 and GridION data. After, genome was annotated using Prokka v1.14.6 and their coding sequences (CDS) were identified [14].

Results and Discussion

1. Genome structure and size

The basic genome statistics are provided in Table 1. The complete genome of L. brevis NS2301G3 consists of one circular chromosome (2,449,247 bp) with a GC content of 46.01%. According to the genomic results, L. brevis KL251 contains 2,412 CDSs, 67 tRNAs, and 15 rRNAs (Fig. 1).

Table 1. Genomic features of Levilactobacillus brevis NS2301G3
Genomic features L. brevis NS2301G3 (chromosome)
Genome size (bp) 2,449,247
GC content (%) 46.01
rRNA genes 15
tRNA genes 67
CDS 2,412

GC, guanine and cytosine bases; rRNA, ribosomal RNA; tRNA, transfer RNA; CDS, coding sequences.

Download Excel Table
jmsb-42-4-161-g1
Fig. 1. Circular chromosome map of Levilactobacillus brevis NS2301G3. CDS, coding sequences; tRNA, transfer RNA; rRNA, ribosomal RNA; tmRNA, transfer messenger RNA, GC, guanine and cytosine bases.
Download Original Figure
2. Gene content and functional annotation

The complete genome of L. brevis NS2301G3 was analyzed to determine the gene content and its functional roles. Genome annotation was performed using the RAST server and Prokka annotation tools, providing an in-depth overview of the functional repertoire of this strain. The annotation revealed a total of 2,789 protein-CDSs, including genes involved in various metabolic pathways, stress response mechanisms, and cell structure maintenance. The annotated genome also identified several genes encoding carbohydrate-active enzymes (CAZymes), suggesting a strong potential for carbohydrate metabolism, which aligns with the species’ role in fermented foods [15]. Notably, genes responsible for the synthesis of lactic acid, such as lactate dehydrogenase (ldh), were prominently identified, emphasizing its importance in lactic acid fermentation [16]. Additionally, genes related to the catabolism of sugars such as glucose, fructose, and galactose were found, supporting the strain’s adaptability to diverse carbohydrate sources [17]. Key genes involved in carbohydrate metabolism include glucokinase (glk), phosphofructokinase (pfk), aldolase (ald ), beta-galactosidase (lacZ ), lactose permease (lacY ), and galactoside acetyltransferase (lacA) [18]. Functional annotation further revealed genes associated with stress tolerance, including those encoding heat shock proteins (DnaA, DnaN, DnaK, GroS, GroL) and oxidative stress response enzymes, which likely contribute to the strain's ability to survive and adapt to various environmental conditions, such as the acidic and anaerobic environment of fermented foods [19,20]. Additionally, genes involved in the biosynthesis of vitamins, such as riboflavin (ribBA, ribD) and folate (folT, folD, folC, folE, folK, folB), were detected, suggesting potential health-promoting properties [21,22]. In terms of antimicrobial properties, genes encoding CRISPR-associated proteins (casC, cas3), bacteriocins, and other antimicrobial peptides were found, indicating that L. brevis NS2301G3 may play a role in inhibiting pathogenic microorganisms in fermented food systems [23,24]. Genes associated with exopolysaccharide (EPS) production, such as pspA1, pspA2, pspB, were also identified, which are known to enhance the texture and viscosity of fermented products, thus contributing to the overall quality of the final product [25,26]. The genome also contains genes related to fatty acid metabolism regulation (fadR), arabinose metabolism repression (araR), and malolactic enzyme (mleS ), highlighting its versatility in metabolizing various substrates [27,28]. Genes involved in glutathione biosynthesis (gshAB), associated with antioxidant functions, were detected, further supporting the strain’s stress tolerance [29,30]. Additionally, genes related to the GABA pathway were identified, including glutamate/GABA antiporter (gadC ), glutamate decarboxylase (gadB), succinic semialdehyde dehydrogenase (gabD1), glutamate-tRNA ligase (gltX ), and proton/sodium-glutamate symporter (gltT ), suggesting the strain’s potential in GABA production, which is linked to health benefits such as stress reduction [3133]. The presence of mobile genetic elements, including transposases and prophage sequences, was also noted, which may suggest genomic plasticity and the ability to acquire new traits through horizontal gene transfer [34]. This feature could provide an adaptive advantage in diverse ecological niches. Genes related to cell wall and mucosal adhesion, such as Eno, Eno2, LDH1, LDH2, LDH3, were also detected, which are likely important for the strain’s ability to interact with the host environment, contributing to its potential probiotic effects [35,36]. Overall, the gene content and functional annotation of L. brevis NS2301G3 reveal a versatile metabolic capacity, stress tolerance, and potential probiotic properties, supporting its application in fermented food production and its potential as a beneficial microbe for human health.

Conflict of Interest

The authors declare no potential conflict of interest.

Acknowledgements

This work was supported by Nong Shim Co. Ltd. And this work was also supported by Korea Institute of Planning and Evaluation for Technology in Food, Agriculture and Forestry (IPET) through the Agricultural Microbiome R&D Program for Advancing innovative technology Program, funded by Ministry of Agriculture, Food and Rural Affairs (MAFRA) (RS-2024-0040347740982119420101).

The complete chromosome, plasmid sequences of L. brevis NS2301G3 have been deposited in GenBank under the accession numbers CP095737, CP165721-CP165724, respectively. The BioProject and BioSample accession numbers are PRJNA1141725 and SAMN42962384.

References

1.

Silva CCG, Ribeiro SC. Microorganisms and their importance in the food industry: safety, quality, and health properties. Foods. 2024;13:1452.

2.

Park JM, Moon JW, Zhang BZ, An BK. Antioxidant activity and other characteristics of lactic acid bacteria isolated from Korean traditional sweet potato stalk kimchi. Foods. 2024;20:3261.

3.

Fan X, Yu L, Shi Z, Li C, Zeng X, Wu Z, et al. Characterization of a novel flavored yogurt enriched in γ-aminobutyric acid fermented by Levilactobacillus brevis CGMCC1.5954. J Dairy Sci. 2023;106:852-867.

4.

Wu Q, Shah NP. High γ-aminobutyric acid production from lactic acid bacteria: emphasis on Lactobacillus brevis as a functional dairy starter. Crit Rev Food Sci Nutr. 2017;57:3661-3672.

5.

Kim K, Shin DJ, Lee J, Oh S. Complete genome sequence of Levilactobacillus brevis KL251 isolate from kimchi. J Dairy Sci Biotechnol. 2024;42:18-22.

6.

Seo MJ, Nam YD, Park SL, Lee SY, Yi SH, Lim SI. γ-Aminobutyric acid production in skim milk co-fermented with Lactobacillus brevis 877G and Lactobacillus sakei 795. Food Sci Biotechnol. 2013;22:751-755.

7.

De Man JC, Rogosa M, Sharpe ME. A medium for the cultivation of lactobacilli. J Appl Bacteriol. 1960;23:130-135.

8.

Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114-2120.

9.

Andrews S. FastQC: a quality control tool for high throughput sequence data [Internet]. 2010 [cited 2024 Oct 6]. Available from: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/

10.

Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. 2013. arXiv.

11.

Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, et al. Twelve years of SAMtools and BCFtools. Gigascience. 2021;10:giab008.

12.

De Coster W, D’Hert S, Schultz DT, Cruts M, Van Broeckhoven C. NanoPack: visualizing and processing long-read sequencing data. Bioinformatics. 2018;34: 2666-2669.

13.

GitHub. Oxford nanopore technologies [Internet]. 2024 [cited 2024 Oct 6]. Available from: https://github.com/nanoporetech

14.

Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLOS Comput. 2017;13:e1005595.

15.

Liang T, Jiang T, Liang Z, Zhang N, Dong B, Wu Q, et al. Carbohydrate-active enzyme profiles of Lactiplantibacillus plantarum strain 84-3 contribute to flavor formation in fermented dairy and vegetable products. Food Chem X. 2023;14: 101036.

16.

Tarrah A, Pakroo S, Lemos Junior WJF, Guerra AF, Corich V, Giacomini A. Complete genome sequence and carbohydrates-active EnZymes (CAZymes) analysis of Lactobacillus paracasei DTA72, a potential probiotic strain with strong capability to use inulin. Curr Microbiol. 2020;77:2867-2875.

17.

Sharma A, Sharma N, Gupta D, Lee HJ, Park YS. Comparative genome analysis of four Leuconostoc strains with a focus on carbohydrate-active enzymes and oligosaccharide utilization pathways. Comput Struct Biotechnol J. 2022;20:4771-4785.

18.

Apostolakos I, Paramithiotis S, Mataragas M. Comparative genomic analysis reveals the functional traits and safety status of lactic acid bacteria retrieved from artisanal cheeses and raw sheep milk. Foods. 2023;12:599.

19.

Zhang C, Gui Y, Chen X, Chen D, Guan C, Yin B, et al. Transcriptional homogenization of Lactobacillus rhamnosus hsryfm 1301 under heat stress and oxidative stress. Appl Microbiol Biotechnol. 2020;104:2611-2621.

20.

Rossi F, Zotta T, Iacumin L, Reale A. Theoretical insight into the heat shock response (HSR) regulation in Lactobacillus casei and L. rhamnosus. J Theor Biol. 2016;390:21-37.

21.

Dijkstra AR, Setyawati MC, Bayjanov JR, Alkema W, van Hijum SA, Bron PA, et al. Diversity in robustness of Lactococcus lactis strains during heat stress, oxidative stress, and spray drying stress. Appl Environ Microbiol. 2014;80:603-611.

22.

Zhang C, Lu J, Yang D, Chen X, Huang Y, Gu R. Stress influenced the aerotolerance of Lactobacillus rhamnosus hsryfm 1301. Biotechnol Lett. 2018;40:729-735.

23.

Mu Y, Zhang C, Li T, Jin FJ, Sung YJ, Oh HM, et al. Development and applications of CRISPR/Cas9-based genome editing in Lactobacillus. Int J Mol Sci. 2022;23: 12852.

24.

Kumar P, Singh S, Sankhyan S, Ray S. Metabolic engineering of lactic acid bacteria for antimicrobial peptides production. Antimicrob Pept Lactic Acid Bact. 2024;3: 67-95.

25.

Plavec TV, Berlec A. Engineering of lactic acid bacteria for delivery of therapeutic proteins and peptides. Appl Microbiol Biotechnol. 2019;103:2053-2066.

26.

Van Pijkeren JP, Barrangou R. Genome editing of food-grade lactobacilli to develop therapeutic probiotics. Microbiol Spectr. 2017;5:BAD-0013-2016.

27.

Franco IS, Mota LJ, Soares CM, de Sá-Nogueira I. Functional domains of the Bacillus subtilis transcription factor AraR and identification of amino acids important for nucleoprotein complex assembly and effector binding. J Bacteriol. 2006;188: 3024-3036.

28.

Kawaguchi H, Sasaki M, Vertès AA, Inui M, Yukawa H. Identification and functional analysis of the gene cluster for L-arabinose utilization in Corynebacterium glutamicum. Appl Environ Microbiol. 2009;75:3419-3429.

29.

Liu X, Zhao J, Zang J, Peng C, Lv L, Li Z. Integrated analysis of physiology, antioxidant activity and transcriptomic of Lactobacillus plantarum 120 in response to acid stress. LWT-Food Sci Technol. 2023:117109.

30.

Pophaly SD, Singh R, Pophaly SD, Kaushik JK, Tomar SK. Current status and emerging role of glutathione in food grade lactic acid bacteria. Microb Cell Fact. 2012;11:114.

31.

Altermann E, Hutkins RW. Gamma-aminobutyric acid and probiotics: multiple health benefits including stress reduction and gut health. J Funct Foods. 2019;57: 94-103.

32.

Mu Y, Zhang C, Li T, Jin FJ, Sung YJ, Oh HM. Production of gamma-aminobutyric acid (GABA) by lactic acid bacteria: enzymatic pathways and health implications. Int J Mol Sci. 2022;23:995.

33.

Yang S, Li X, Zhou W. Comprehensive characterization of GABA production pathways in lactic acid bacteria. Front Microbiol. 2024;15:1408624.

34.

Pei Z, Sadiq FA, Han X, Zhao J, Zhang H, Ross RP, et al. Comprehensive scanning of prophages in Lactobacillus: distribution, diversity, antibiotic resistance genes, and linkages with CRISPR-Cas systems. mSystems. 2021;6:e0121120.

35.

Glenting J, Beck HC, Vrang A, Riemann H, Ravn P, Hansen AM, et al. Anchorless surface associated glycolytic enzymes from Lactobacillus plantarum 299v bind to epithelial cells and extracellular matrix proteins. Microbiol Res. 2013;168:245-253.

36.

Liu W, Wang Z, Wang S, Liu M, Zhang J, Li X, et al. Identification of moonlighting adhesins of highly-adhesive Lactobacillus plantarum PO23 isolated from the intestine of Paralichthys olivaceus. Aquaculture. 2024;558:741044.