Articles Service
Research article
Metabolic Regulation of Longevity and Immune Response in Caenorhabditis elegans by Ingestion of Lacticaseibacillus rhamnosus IDCC 3201 Using Multi-Omics Analysis
1Department of Agricultural Biotechnology and Research Institute of Agriculture and Life Science, Seoul National University, Seoul 08826, Rrepublic of Korea
2Department of Animal Science and Biotechnology, Chungnam National University, Daejeon 34134, Republic of Korea
3Department of Food Science & Nutrition, Dongseo University, Busan 47011, Republic of Korea
4Department of Animal Science and Technology, Chung-Ang University, Anseong 17546, Republic of Korea
5Department of Microbiology, College of Medicine, Dongguk University, Gyeongju 38066,Republic of Korea
6College of Life Sciences, Sejong University, Seoul 05006, Republic of Korea
7Department of Functional Food and Biotechnology, Jeonju University, Jeonju 55069, Republic of Korea
J. Microbiol. Biotechnol. 2024; 34(5): 1109-1118
Published May 28, 2024 https://doi.org/10.4014/jmb.2402.02025
Copyright © The Korean Society for Microbiology and Biotechnology.
Abstract
Keywords
Graphical Abstract
Introduction
Interest and demand for probiotics are increasing due to their long history of safe use in fermented products and the recognition of their profitable effect on gut and health [1]. According to the Food and Agriculture Organization and the World Health Organization, the current definition of probiotics is “Live microorganisms which confer a health benefit when administered in adequate amounts [2]. Probiotics are acknowledged for their beneficial effects on host health, including their ability to compete with pathogenic bacteria, produce antimicrobial compounds, enhance the integrity of the host’s mucosal barrier, and modulate the immune system [3-6].
The whole transcriptomic analysis is widely used to examine the change in multiple genetic pathways in
Several studies in
Material and Methods
Bacterial and C. elegans Strains and Culture Conditions
In vivo Adhesive Assay Using C. elegans
To evaluate whether 3201 can colonize in
Lifespan and Killing Assay
To evaluate whether 3201 influences
About the killing assay, L4 stage
Worm Size and Locomotive Activity
Worm size and locomotive activity were measured using Wormlab software (MBF Bioscience, USA) as previously reported with slight modifications [19]. Briefly, L4 stage
RNA Isolation and Transcriptome Analysis
Transcriptome analysis was conducted with a slight modification of the previous study [20]. Briefly, L4 stage
Metabolites Extraction and Metabolome Analysis
To measure the change of metabolites, metabolome analysis was performed with a slight modification to a previous study [25]. L4 stage
GC-MS analysis was conducted using a TRACE 1310 Gas Chromatograph (Thermo Fisher Scientific, USA) with an ISQ LT single quadrupole mass spectrometer (Thermo Fisher Scientific), and separation was conducted using a DB-5MS column (60 m × 0.25 mm, 0.25 μm film thickness, Agilent, USA). The oven temperature was held at 50°C for 2 min, constantly increased at 5°C/min, held at 180°C for 8 min, gradually increased at 2.5°C/min to 325°C, and then held for 10 min. The sample was injected at 300°C, and helium was used as the carrier gas, with a flow rate of 1.5 ml/min and a split ratio of 1:60. For GC-MS detection, an electron ionization system with an ionization voltage of 70 eV was used, and the temperature of the ion source was 270°C. The mass scan range was set at 30–450 (m/z), and the acquisition rate was 5 spectra/s. NIST Mass Spectral Search Program (version 2.0; National Institute of Standards and Technology [NIST], USA) was used to identify detected metabolites. Further analyses were performed using MetaboAnalyst 6.0.
Statistics
Data Availability
All the data needed to assess the conclusions of the paper are included in the manuscript and deposited in the NCBI SRA database under the Bioproject number PRJNA1092403. Additional data can be obtained from the authors upon request.
Results
Evaluation of the Adhesive Ability of L. rhamnosus IDCC 3201 in C. elegans
We first evaluated the adhesive ability of probiotic candidate bacteria
-
Fig. 1. Evaluation of the adhesive ability of
Lacticaseibacillus rhamnosus IDCC 3201 inC. elegans . The adhesive ability of OP50, 3201, or LGG inC. elegans strainfer-15;fem-1 after 48 h conditioning period. OP50,E. coli ; 3201,L. rhamnosus IDCC 3201; LGG,L. rhamnosus GG. Statistical analysis was performed using one-way ANOVA, and statistical significance was considered when thep values were below 0.05 (*), 0.01 (**), 0.001 (***), and 0.0001 (****). Statistics compared with 3201:p < 0.0001 andp = 0.3699 for OP50 and LGG, respectively. Data are expressed as means ± SEM.
Evaluation of Longevity and Immune Response of L. rhamnosus IDCC 3201 Using C. elegans
We examined whether 3201 can improve the lifespan of
-
Fig. 2. Evaluation of longevity and immune response of
Lacticaseibacillus rhamnosus IDCC 3201 usingC. elegans . The lifespan ofC. elegans strainfer-15;fem-1 using OP50, 3201, and LGG. The killing assay ofC. elegans strain fer- 15;fem-1 preconditioned with OP50, 3201, or LGG for 48 h then infected with each of four foodborne pathogenic bacteria (two of each gram-negative and positive bacteria). (A) Lifespan assay. (B) Killing assay usingEscherichia coli O157:H7 EDL933. (C) Killing assay usingSalmonella Typhimurium SL1344. (D) Killing assay usingStaphylococcus aureus Newman. (E) Killing assay usingListeria monocytogenes EGD-e. OP50,E. coli ; 3201,L. rhamnosus ; LGG,L. rhamnosus . Statistical analysis was performed using the Kaplan–Meier method, and differences were considered significant when thep value was below 0.05 (*) and 0.01 (**). Survival statistic in lifespan assay compared with 3201:p = 0.0000 andp = 0.0000 for OP50 and LGG, respectively. Survival statistics in killing assay compared with 3201:Escherichia coli O157:H7 EDL933,p = 0.0114 andp = 0.5001 for OP50 and LGG, respectively;Salmonella Typhimurium SL1344,p = 0.0000 andp = 0.2732 for OP50 and LGG, respectively;Staphylococcus aureus Newman,p = 0.0000 andp = 0.0000 for OP50 and LGG, respectively;Listeria monocytogenes EGD-e,p = 0.0000 andp = 0.0496 for OP50 and LGG, respectively. Data are expressed as means ±SEM .
Subsequently, we investigated the effect of 3201 on the host’s immune response against four different kinds of foodborne pathogenic bacteria.
Evaluation of Phenotypic Change of C. elegans after Exposure to L. rhamnosus IDCC 3201
We measured the body size to evaluate whether feeding 3201 changed the phenotype of worms. The body size was measured after feeding with OP50, 3201, or LGG for 48 h. The results showed that the 3201 group demonstrated a significant increase in both length and width compared with both the OP50 group (
-
Fig. 3. Evaluation of phenotypic change of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201. Worm size and locomotive activity ofC. elegans strainfer-15;fem-1 after 48 h exposure period with OP50, 3201, or LGG. (A) Length, (B) Width, (C) Peristaltic speed, and (D) Pumping rate. OP50,E. coli ; 3201,L. rhamnosus ; LGG,L. rhamnosus . Statistical analysis was performed using one-way ANOVA, and statistical significance was considered when thep values were below 0.05 (*), 0.01 (**), 0.001 (***), and 0.0001 (****). Statistics compared with 3201: length,p < 0.0001 andp < 0.0001 for OP50 and LGG, respectively; width,p < 0.0001 andp = 0.0028 for OP50 and LGG, respectively; peristaltic speed,p = 0.5366 andp = 0.6957 for OP50 and LGG, respectively; pumping rate,p < 0.0001 andp = 0.0002 for OP50 and LGG, respectively. Data are expressed as means ±SEM .
Evaluation of Longevity and Immune Response of L. rhamnosus IDCC 3201 Using C. elegans
We examined whether 3201 can improve the lifespan of
-
Fig. 2. Evaluation of longevity and immune response of
Lacticaseibacillus rhamnosus IDCC 3201 usingC. elegans . The lifespan ofC. elegans strainfer-15;fem-1 using OP50, 3201, and LGG. The killing assay ofC. elegans strain fer- 15;fem-1 preconditioned with OP50, 3201, or LGG for 48 h then infected with each of four foodborne pathogenic bacteria (two of each gram-negative and positive bacteria). (A) Lifespan assay. (B) Killing assay usingEscherichia coli O157:H7 EDL933. (C) Killing assay usingSalmonella Typhimurium SL1344. (D) Killing assay usingStaphylococcus aureus Newman. (E) Killing assay usingListeria monocytogenes EGD-e. OP50,E. coli ; 3201,L. rhamnosus ; LGG,L. rhamnosus . Statistical analysis was performed using the Kaplan–Meier method, and differences were considered significant when thep value was below 0.05 (*) and 0.01 (**). Survival statistic in lifespan assay compared with 3201:p = 0.0000 andp = 0.0000 for OP50 and LGG, respectively. Survival statistics in killing assay compared with 3201:Escherichia coli O157:H7 EDL933,p = 0.0114 andp = 0.5001 for OP50 and LGG, respectively;Salmonella Typhimurium SL1344,p = 0.0000 andp = 0.2732 for OP50 and LGG, respectively;Staphylococcus aureus Newman,p = 0.0000 andp = 0.0000 for OP50 and LGG, respectively;Listeria monocytogenes EGD-e,p = 0.0000 andp = 0.0496 for OP50 and LGG, respectively. Data are expressed as means ±SEM .
Subsequently, we investigated the effect of 3201 on the host’s immune response against four different kinds of foodborne pathogenic bacteria.
Evaluation of Phenotypic Change of C. elegans after Exposure to L. rhamnosus IDCC 3201
We measured the body size to evaluate whether feeding 3201 changed the phenotype of worms. The body size was measured after feeding with OP50, 3201, or LGG for 48 h. The results showed that the 3201 group demonstrated a significant increase in both length and width compared with both the OP50 group (
-
Fig. 3. Evaluation of phenotypic change of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201. Worm size and locomotive activity ofC. elegans strainfer-15;fem-1 after 48 h exposure period with OP50, 3201, or LGG. (A) Length, (B) Width, (C) Peristaltic speed, and (D) Pumping rate. OP50,E. coli ; 3201,L. rhamnosus ; LGG,L. rhamnosus . Statistical analysis was performed using one-way ANOVA, and statistical significance was considered when thep values were below 0.05 (*), 0.01 (**), 0.001 (***), and 0.0001 (****). Statistics compared with 3201: length,p < 0.0001 andp < 0.0001 for OP50 and LGG, respectively; width,p < 0.0001 andp = 0.0028 for OP50 and LGG, respectively; peristaltic speed,p = 0.5366 andp = 0.6957 for OP50 and LGG, respectively; pumping rate,p < 0.0001 andp = 0.0002 for OP50 and LGG, respectively. Data are expressed as means ±SEM .
Transcriptomic Analysis of C. elegans after Exposure to L. rhamnosus IDCC 3201
We conducted a transcriptome analysis to evaluate the altered mechanism through feeding 3201 compared with OP50. We selected genes that showed a significant upregulation of more than a 2-fold change in
-
Table 1 . Transcriptomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201.Term Gene count % p valueInnate immune response 75 4.8 0.000 Glutathione metabolic process 18 1.1 0.000 Peptidoglycan catabolic process 8 0.5 0.000 Cell wall macromolecule catabolic process 9 0.6 0.000 Defense response to gram-positive bacterium 17 1.1 0.000 Lipid metabolic process 32 2 0.000 Proteolysis 42 2.7 0.000 Defense response to gram-negative bacterium 19 1.2 0.000 Dephosphorylation 13 0.8 0.000 Fatty acid beta-oxidation using acyl-CoA oxidase 5 0.3 0.000 Top 10 pathways closely associated with genes significantly upregulated more than 2.0 folds in
C. elegans strainfer-15;fem-1 following 48 h of exposure toL. rhamnosus 3201 compared withE. coli OP50.
-
Table 2 . Transcriptomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201.Group and gene Gene number Fold change p valueDescription C-type lectin-related clec-41 CELE_B0365.6 8.443320 0.000 C-type lectin clec-62 CELE_F35C5.5 2.348795 0.000 C-type lectin clec-86 CELE_C54D1.2 9.677637 0.000 C-type lectin clec-186 CELE_ZK896.7 6.833677 0.000 C-type lectin clec-187 CELE_ZK896.6 4.984468 0.000 C-type lectin Lysozyme-related lys-1 CELE_Y22F5A.4 3.182571 0.000 Lysozyme lys-2 CELE_Y22F5A.5 8.078466 0.000 Lysozyme lys-5 CELE_F22A3.6 12.360469 0.000 Lysozyme lys-7 CELE_C02A12.4 5.542044 0.000 Lysozyme lys-8 CELE_C17G10.5 2.706994 0.000 Lysozyme List of genes that were significantly upregulated by more than 2.0 folds in
C. elegans when fed withL. rhamnosus 3201 compared withE. coli OP50 and associated with the innate immune response pathway.
-
Fig. 4. Transcriptomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201. The identification of KEGG pathways related to genes significantly upregulated by more than 3-fold inC. elegans when fed with 3201 compared with OP50. We utilized Cytoscape for the analysis.
Metabolomic Analysis of C. elegans after Exposure to L. rhamnosus IDCC 3201
The metabolomic analysis was performed to assess the impact of 3201 on the metabolite composition of
-
Fig. 5. Metabolomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201. Change of metabolites composition inC. elegans strainfer-15;fem-1 after 48 h exposure period with OP50 or 3201. (A) PLS-DA. (B) Volcano plot. (C) Top 25 enriched heatmap. (D) The quantitative graph using the metabolites which changed more than 4.0 folds.
Discussion
According to the previous study, probiotics are known to exert their beneficial effects on the host by adhering to host epithelial cells, modulating the host’s immune response, and producing specific metabolites [26].
According to a previous study, the adhesive ability to the host’s gastrointestinal is one of the classical selection criteria for potential probiotic bacteria. This adhesive capability can lead to colonization, which can help improve immunomodulatory effects and stimulate the gut barrier and metabolic functions [32]. Our study revealed that 3201 exhibits a significantly higher adhesive ability than OP50 and shows no significant difference compared with LGG. In lifespan assay, 3201 significantly improved the longevity of
According to previous studies, the quality of food and some types of bacteria are known to affect the phenotype of worms such as showing better growth performance [34-36]. Therefore, we sought to assess the impact of feeding 3201 on various phenotypic traits. Measurements of body size, peristaltic speed, and pumping rate were employed as indicators of potential changes in the worms’ phenotype. The results revealed a significant increase in both length and width when
Our previous findings revealed that preconditioning with 3201 enhanced resistances against both gram-negative and gram-positive foodborne pathogens. Therefore, we speculated that 3201 may modulate specific gene expression related to the immune response. Transcriptome analysis using genes that significantly increased more than 2-fold in
The metabolome analysis revealed increased levels of fumaric acid and succinic acid, metabolites of the citric acid cycle, in
In conclusion, our study focused on evaluating the probiotic properties of
Acknowledgements
This study was supported by the Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, Forestry and Fisheries (IPET-321037-5) and by a National Research Foundation of Korea Grant funded by the Korean government (MEST) (NRF-2021R1A2C3011051) and “Cooperative Research Program for Agriculture Science & Technology Development (Project No. RS-2023-00225838)” Rural Development Administration, Republic of Korea.
Conflict of Interest
The authors have no financial conflicts of interest to declare.
References
- Kim B, Kim K, Xu X, Lee H, Pathiraja D, Park DJ,
et al . 2023. Complete genome and two plasmids sequences ofLactiplantibacillus plantarum L55 for probiotic potentials.J. Anim. Sci. Technol. 65 : 1341-1344. - Morelli L, Capurso L. 2012. FAO/WHO guidelines on probiotics: 10 years later.
J. Clin. Gastroenterol. 46 : S1-S2. - Lee D, Goh TW, Kang MG, Choi HJ, Yeo SY, Yang J,
et al . 2022. Perspectives and advances in probiotics and the gut microbiome in companion animals.J. Anim. Sci. Technol. 64 : 197-217. - Yang S, Deng C, Li Y, Li W, Wu Q, Sun Z,
et al . 2022. Complete genome sequence ofLactiplantibacillus plantarum ST, a potential probiotic strain with antibacterial properties.J. Anim. Sci. Technol. 64 : 183-186. - Kim B, Meng Z, Xu X, Baek S, Pathiraja D, Choi IG,
et al . 2023. Complete genome sequence ofLimosilactobacillus fermentum JNU532 as a probiotic candidate for the functional food and feed supplements.J. Anim. Sci. Technol. 65 : 271-274. - Kim JY, Kim JY, Kim H, Moon EC, Heo K, Shim JJ,
et al . 2022. Immunostimulatory effects of dairy probiotic strainsBifidobacterium animalis ssp. HY8002 andLactobacillus plantarum HY7717.J. Anim. Sci. Technol. 64 : 1117-1131. - Kang AN, Mun D, Ryu S, Jae Lee J, Oh S, Kyu Kim M,
et al . 2022. Culturomic-, metagenomic-, and transcriptomic-based characterization of commensal lactic acid bacteria isolated from domestic dogs usingCaenorhabditis elegans as a model for aging.J. Anim. Sci. 100 : skac323. - Jaafar M, Xu P, Mageswaran UM, Balasubramaniam SD, Solayappan M, Woon J,
et al . 2024. Dairy based-LAB improved constipation via increasing fecal bulk and decreasing concentration of fecal threonine while preserving colonic goblet cell count.J. Anim. Sci. Technol. 66 : 178-203. - Kaletta T, Hengartner MO. 2006. Finding function in novel targets:
C. elegans as a model organism.Nat. Rev. Drug discov. 5 : 387-398. - Mahesh R, Ilangovan P, Nongbri D, Suchiang K. 2021. Probiotics interactions and the modulation of major signalling pathways in host model organism
Caenorhabditis elegans .Indian J. Microbiol. 61 : 404-416. - Clark LC, Hodgkin J. 2014. Commensals, probiotics and pathogens in the
Caenorhabditis elegans model.Cell. Microbiol. 16 : 27-38. - Choi H, Mun D, Ryu S, Kwak Mj, Kim BK, Park DJ,
et al . 2023. Molecular characterization and functionality of rumen-derived extracellular vesicles using aCaenorhabditis elegans animal model.J. Anim. Sci. Technol. 65 : 652-663. - Kang A, Kwak MJ, Lee DJ, Lee JJ, Kim MK, Song M,
et al . Dietary supplementation with probiotics promotes weight loss by reshaping the gut microbiome and energy metabolism in obese dogs. Microbiol.Spectr. 0 : e02552-02523. - Anderson SM, Pukkila-Worley R. 2020. Immunometabolism in
Caenorhabditis elegans .PLoS Pathog. 16 : e1008897. - Wan QL, Shi X, Liu J, Ding AJ, Pu YZ, Li Z,
et al . 2017. Metabolomic signature associated with reproduction-regulated aging inCaenorhabditis elegans .Aging (Albany NY) 9 : 447-474. - Kim H, Shin M, Ryu S, Yun B, Oh S, Park D-J,
et al . 2021. Evaluation of probiotic characteristics of newly isolated lactic acid bacteria from dry-aged Hanwoo beef.Food Sci. Anim. Resour. 41 : 468-480. - Park MR, Ryu S, Maburutse BE, Oh NS, Kim SH, Oh S,
et al . 2018. ProbioticLactobacillus fermentum strain JDFM216 stimulates the longevity and immune response ofCaenorhabditis elegans through a nuclear hormone receptor.Sci. Rep. 8 : 7441. - Kim Y, Mylonakis E. 2012.
Caenorhabditis elegans immune conditioning with the probiotic bacteriumLactobacillus acidophilus strain NCFM enhances gram-positive immune responses.Infect. Immun. 80 : 2500-2508. - Shen P, Kershaw JC, Yue Y, Wang O, Kim KH, McClements DJ,
et al . 2018. Effects of conjugated linoleic acid (CLA) on fat accumulation, activity, and proteomics analysis inCaenorhabditis elegans .Food Chem. 249 : 193-201. - Ryu S, Shin M, Yun B, Lee W, Choi H, Kang M,
et al . 2021. Bacterial quality, prevalence of pathogens, and molecular characterization of biofilm-producingStaphylococcus aureus from Korean dairy farm environments.Animals 11 : 1306. - Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data.
Bioinformatics 30 : 2114-2120. - Kim D, Langmead B, Salzberg SL. 2015. HISAT: a fast spliced aligner with low memory requirements.
Nat. Methods 12 : 357-360. - Liao Y, Smyth GK, Shi W. 2014. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features.
Bioinformatics 30 : 923-930. - Robinson MD, McCarthy DJ, Smyth GK. 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.
Bioinformatics 26 : 139-140. - Yoo J, Lee J, Zhang M, Mun D, Kang M, Yun B,
et al . 2022. Enhanced γ-aminobutyric acid and sialic acid in fermented deer antler velvet and immune promoting effects.J. Anim. Sci. Technol. 64 : 166-182. - Williams NT. 2010. Probiotics.
Am. J. Health Syst. Pharm. 67 : 449-458. - Oh YJ, Lee J, Lim SK, Kwon MS, Lee S, Choi SP,
et al . 2023. Complete genome sequence of probioticLactobacillus johnsonii 7409N31 isolated from a healthy Hanwoo calf.J. Anim. Sci. Technol. 65 : 890-893. - Hammes WP, Vogel RF. 1995. The genus lactobacillus, pp. 19-54.
In: The genera of lactic acid bacteria . Ed. Springer,. - Vasquez R, Oh JK, Song JH, Kang DK. 2022. Gut microbiome-produced metabolites in pigs: a review on their biological functions and the influence of probiotics.
J. Anim. Sci. Technol. 64 : 671-695. - Grompone G, Martorell P, Llopis S, González N, Genovés S, Mulet AP,
et al . 2012. Anti-inflammatoryLactobacillus rhamnosus CNCM I-3690 strain protects against oxidative stress and increases lifespan inCaenorhabditis elegans .PLoS One 7 : e52493. - Oh S, Park MR, Son SJ, Kim Y. 2015.
Lactobacillus rhamnosus strain GG influences on the longevity and immune response ofCaenorhabditis elegans to foodborne pathogens with stimulating microRNAs.FASEB J. 29 : 853. - Monteagudo-Mera A, Rastall RA, Gibson GR, Charalampopoulos D, Chatzifragkou A. 2019. Adhesion mechanisms mediated by probiotics and prebiotics and their potential impact on human health.
Appl. Microbiol. Biotechnol. 103 : 6463-6472. - Poupet C, Saraoui T, Veisseire P, Bonnet M, Dausset C, Gachinat M,
et al . 2019.Lactobacillus rhamnosus Lcr35 as an effective treatment for preventingCandida albicans infection in the invertebrate modelCaenorhabditis elegans : first mechanistic insights.PLoS One 14 : e0216184. - Shtonda BB, Avery L. 2006. Dietary choice behavior in
Caenorhabditis elegans .J. Exper. Biol. 209 : 89-102. - Avery L, You YJ.
C. elegans feeding. WormBook: The Online Review ofC. elegans Biology [Internet], 2018. - Kim YJ, Cho SB, Song MH, Lee SI, Hong SM, Yun W,
et al . 2022. Effects of differentBacillus licheniformis andBacillus subtilis ratios on nutrient digestibility, fecal microflora, and gas emissions of growing pigs.J. Anim. Sci. Technol. 64 : 291-301. - Pees B, Yang W, Kloock A, Petersen C, Peters L, Fan L,
et al . 2021. Effector and regulator: diverse functions ofC. elegans C-type lectinlike domain proteins.PLoS Pathog. 17 : e1009454. - Schulenburg H, Hoeppner MP, Weiner III J, Bornberg-Bauer E. 2008. Specificity of the innate immune system and diversity of Ctype lectin domain (CTLD) proteins in the nematode
Caenorhabditis elegans .Immunobiology 213 : 237-250. - O'Rourke D, Baban D, Demidova M, Mott R, Hodgkin J. 2006. Genomic clusters, putative pathogen recognition molecules, and antimicrobial genes are induced by infection of
C. elegans withM. nematophilum .Genome Res. 16 : 1005-1016. - Leippe M. 1999. Antimicrobial and cytolytic polypeptides of amoeboid protozoa-effector molecules of primitive phagocytes.
Dev. Comp. Immunol. 23 : 267-279. - Ciancio A. 2016. Defense and Immune Systems, pp. 205-239.
In: Invertebrate Bacteriology . Ed. Springer,. - Amador‐Noguez D, Yagi K, Venable S, Darlington G. 2004. Gene expression profile of long‐lived Ames dwarf mice and Little mice.
Aging Cell 3 : 423-441. - Perron JT, Tyson RL, Sutherland GR. 2000. Maintenance of tricarboxylic acid cycle kinetics in Brown-Norway Fischer 344 rats may translate to longevity.
Neurosc. Lett. 281 : 91-94. - Edwards CB, Copes N, Brito AG, Canfield J, Bradshaw PC. 2013. Malate and fumarate extend lifespan in
Caenorhabditis elegans .PLoS One 8 : e58345. - Zarse K, Schulz TJ, Birringer M, Ristow M. 2007. Impaired respiration is positively correlated with decreased life span in
Caenorhabditis elegans models of Friedreich Ataxia.FASEB J. 21 : 1271-1275. - Ryu S, Kim K, Cho D-Y, Kim Y, Oh S. 2022. Complete genome sequences of
Lactococcus lactis JNU 534, a potential food and feed preservative.J. Anim. Sci. Technol. 64 : 599-602. - Fernández MF, Boris S, Barbes C. 2003. Probiotic properties of human
lactobacilli strains to be used in the gastrointestinal tract.J. Appl. Microbiol. 94 : 449-455. - Dinev T, Beev G, Tzanova M, Denev S, Dermendzhieva D, Stoyanova A. 2018. Antimicrobial activity of
Lactobacillus plantarum against pathogenic and food spoilage microorganisms: a review.Bulg. J. Vet. Med. 21 . DOI:10.15547/bjvm.1084. - Shu Q, Gill HS. 2002. Immune protection mediated by the probiotic
Lactobacillus rhamnosus HN001 (DR20™) againstEscherichia coli O157: H7 infection in mice.FEMS Immunol. Med. Microbiol. 34 : 59-64.
Related articles in JMB
Article
Research article
J. Microbiol. Biotechnol. 2024; 34(5): 1109-1118
Published online May 28, 2024 https://doi.org/10.4014/jmb.2402.02025
Copyright © The Korean Society for Microbiology and Biotechnology.
Metabolic Regulation of Longevity and Immune Response in Caenorhabditis elegans by Ingestion of Lacticaseibacillus rhamnosus IDCC 3201 Using Multi-Omics Analysis
Daniel Junpyo Lee1†, Ju Young Eor1†, Min-Jin Kwak1, Junbeom Lee1, An Na Kang1, Daye Mun1, Hyejin Choi1, Minho Song2, Jong Nam Kim3, Jun-Mo Kim4, Jungwoo Yang5, Hyung Wook Kim6, Sangnam Oh7*, and Younghoon Kim1*
1Department of Agricultural Biotechnology and Research Institute of Agriculture and Life Science, Seoul National University, Seoul 08826, Rrepublic of Korea
2Department of Animal Science and Biotechnology, Chungnam National University, Daejeon 34134, Republic of Korea
3Department of Food Science & Nutrition, Dongseo University, Busan 47011, Republic of Korea
4Department of Animal Science and Technology, Chung-Ang University, Anseong 17546, Republic of Korea
5Department of Microbiology, College of Medicine, Dongguk University, Gyeongju 38066,Republic of Korea
6College of Life Sciences, Sejong University, Seoul 05006, Republic of Korea
7Department of Functional Food and Biotechnology, Jeonju University, Jeonju 55069, Republic of Korea
Correspondence to:Sangnam Oh, osangnam@jj.ac.kr
Younghoon Kim, ykeys2584@snu.ac.kr
†These authors contributed equally to this study.
Abstract
Probiotics, specifically Lacticaseibacillus rhamnosus, have garnered attention for their potential health benefits. This study focuses on evaluating the probiotic properties of candidate probiotics L. rhamnosus IDCC 3201 (3201) using the Caenorhabditis elegans surrogate animal model, a well-established in vivo system for studying host–bacteria interactions. The adhesive ability to the host’s gastrointestinal tract is a crucial criterion for selecting potential probiotic bacteria. Our findings demonstrated that 3201 exhibits significantly higher adhesive capabilities compared with Escherichia coli OP50 (OP50), a standard laboratory food source for C. elegans and is comparable with the widely recognized probiotic L. rhamnosus GG (LGG). In lifespan assay, 3201 significantly increased the longevity of C. elegans compared with OP50. In addition, preconditioning with 3201 enhanced C. elegans immune response against four different foodborne pathogenic bacteria. To uncover the molecular basis of these effects, transcriptome analysis elucidated that 3201 modulates specific gene expression related to the innate immune response in C. elegans. C-type lectin-related genes and lysozyme-related genes, crucial components of the immune system, showed significant upregulation after feeding 3201 compared with OP50. These results suggested that preconditioning with 3201 may enhance the immune response against pathogens. Metabolome analysis revealed increased levels of fumaric acid and succinic acid, metabolites of the citric acid cycle, in C. elegans fed with 3201 compared with OP50. Furthermore, there was an increase in the levels of lactic acid, a well-known antimicrobial compound. This rise in lactic acid levels may have contributed to the robust defense mechanisms against pathogens. In conclusion, this study demonstrated the probiotic properties of the candidate probiotic L. rhamnosus IDCC 3201 by using multi-omics analysis.
Keywords: Lacticaseibacillus rhamnosus, C. elegans, immune response, multi-omics analysis
Introduction
Interest and demand for probiotics are increasing due to their long history of safe use in fermented products and the recognition of their profitable effect on gut and health [1]. According to the Food and Agriculture Organization and the World Health Organization, the current definition of probiotics is “Live microorganisms which confer a health benefit when administered in adequate amounts [2]. Probiotics are acknowledged for their beneficial effects on host health, including their ability to compete with pathogenic bacteria, produce antimicrobial compounds, enhance the integrity of the host’s mucosal barrier, and modulate the immune system [3-6].
The whole transcriptomic analysis is widely used to examine the change in multiple genetic pathways in
Several studies in
Material and Methods
Bacterial and C. elegans Strains and Culture Conditions
In vivo Adhesive Assay Using C. elegans
To evaluate whether 3201 can colonize in
Lifespan and Killing Assay
To evaluate whether 3201 influences
About the killing assay, L4 stage
Worm Size and Locomotive Activity
Worm size and locomotive activity were measured using Wormlab software (MBF Bioscience, USA) as previously reported with slight modifications [19]. Briefly, L4 stage
RNA Isolation and Transcriptome Analysis
Transcriptome analysis was conducted with a slight modification of the previous study [20]. Briefly, L4 stage
Metabolites Extraction and Metabolome Analysis
To measure the change of metabolites, metabolome analysis was performed with a slight modification to a previous study [25]. L4 stage
GC-MS analysis was conducted using a TRACE 1310 Gas Chromatograph (Thermo Fisher Scientific, USA) with an ISQ LT single quadrupole mass spectrometer (Thermo Fisher Scientific), and separation was conducted using a DB-5MS column (60 m × 0.25 mm, 0.25 μm film thickness, Agilent, USA). The oven temperature was held at 50°C for 2 min, constantly increased at 5°C/min, held at 180°C for 8 min, gradually increased at 2.5°C/min to 325°C, and then held for 10 min. The sample was injected at 300°C, and helium was used as the carrier gas, with a flow rate of 1.5 ml/min and a split ratio of 1:60. For GC-MS detection, an electron ionization system with an ionization voltage of 70 eV was used, and the temperature of the ion source was 270°C. The mass scan range was set at 30–450 (m/z), and the acquisition rate was 5 spectra/s. NIST Mass Spectral Search Program (version 2.0; National Institute of Standards and Technology [NIST], USA) was used to identify detected metabolites. Further analyses were performed using MetaboAnalyst 6.0.
Statistics
Data Availability
All the data needed to assess the conclusions of the paper are included in the manuscript and deposited in the NCBI SRA database under the Bioproject number PRJNA1092403. Additional data can be obtained from the authors upon request.
Results
Evaluation of the Adhesive Ability of L. rhamnosus IDCC 3201 in C. elegans
We first evaluated the adhesive ability of probiotic candidate bacteria
-
Figure 1. Evaluation of the adhesive ability of
Lacticaseibacillus rhamnosus IDCC 3201 inC. elegans . The adhesive ability of OP50, 3201, or LGG inC. elegans strainfer-15;fem-1 after 48 h conditioning period. OP50,E. coli ; 3201,L. rhamnosus IDCC 3201; LGG,L. rhamnosus GG. Statistical analysis was performed using one-way ANOVA, and statistical significance was considered when thep values were below 0.05 (*), 0.01 (**), 0.001 (***), and 0.0001 (****). Statistics compared with 3201:p < 0.0001 andp = 0.3699 for OP50 and LGG, respectively. Data are expressed as means ± SEM.
Evaluation of Longevity and Immune Response of L. rhamnosus IDCC 3201 Using C. elegans
We examined whether 3201 can improve the lifespan of
-
Figure 2. Evaluation of longevity and immune response of
Lacticaseibacillus rhamnosus IDCC 3201 usingC. elegans . The lifespan ofC. elegans strainfer-15;fem-1 using OP50, 3201, and LGG. The killing assay ofC. elegans strain fer- 15;fem-1 preconditioned with OP50, 3201, or LGG for 48 h then infected with each of four foodborne pathogenic bacteria (two of each gram-negative and positive bacteria). (A) Lifespan assay. (B) Killing assay usingEscherichia coli O157:H7 EDL933. (C) Killing assay usingSalmonella Typhimurium SL1344. (D) Killing assay usingStaphylococcus aureus Newman. (E) Killing assay usingListeria monocytogenes EGD-e. OP50,E. coli ; 3201,L. rhamnosus ; LGG,L. rhamnosus . Statistical analysis was performed using the Kaplan–Meier method, and differences were considered significant when thep value was below 0.05 (*) and 0.01 (**). Survival statistic in lifespan assay compared with 3201:p = 0.0000 andp = 0.0000 for OP50 and LGG, respectively. Survival statistics in killing assay compared with 3201:Escherichia coli O157:H7 EDL933,p = 0.0114 andp = 0.5001 for OP50 and LGG, respectively;Salmonella Typhimurium SL1344,p = 0.0000 andp = 0.2732 for OP50 and LGG, respectively;Staphylococcus aureus Newman,p = 0.0000 andp = 0.0000 for OP50 and LGG, respectively;Listeria monocytogenes EGD-e,p = 0.0000 andp = 0.0496 for OP50 and LGG, respectively. Data are expressed as means ±SEM .
Subsequently, we investigated the effect of 3201 on the host’s immune response against four different kinds of foodborne pathogenic bacteria.
Evaluation of Phenotypic Change of C. elegans after Exposure to L. rhamnosus IDCC 3201
We measured the body size to evaluate whether feeding 3201 changed the phenotype of worms. The body size was measured after feeding with OP50, 3201, or LGG for 48 h. The results showed that the 3201 group demonstrated a significant increase in both length and width compared with both the OP50 group (
-
Figure 3. Evaluation of phenotypic change of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201. Worm size and locomotive activity ofC. elegans strainfer-15;fem-1 after 48 h exposure period with OP50, 3201, or LGG. (A) Length, (B) Width, (C) Peristaltic speed, and (D) Pumping rate. OP50,E. coli ; 3201,L. rhamnosus ; LGG,L. rhamnosus . Statistical analysis was performed using one-way ANOVA, and statistical significance was considered when thep values were below 0.05 (*), 0.01 (**), 0.001 (***), and 0.0001 (****). Statistics compared with 3201: length,p < 0.0001 andp < 0.0001 for OP50 and LGG, respectively; width,p < 0.0001 andp = 0.0028 for OP50 and LGG, respectively; peristaltic speed,p = 0.5366 andp = 0.6957 for OP50 and LGG, respectively; pumping rate,p < 0.0001 andp = 0.0002 for OP50 and LGG, respectively. Data are expressed as means ±SEM .
Evaluation of Longevity and Immune Response of L. rhamnosus IDCC 3201 Using C. elegans
We examined whether 3201 can improve the lifespan of
-
Figure 2. Evaluation of longevity and immune response of
Lacticaseibacillus rhamnosus IDCC 3201 usingC. elegans . The lifespan ofC. elegans strainfer-15;fem-1 using OP50, 3201, and LGG. The killing assay ofC. elegans strain fer- 15;fem-1 preconditioned with OP50, 3201, or LGG for 48 h then infected with each of four foodborne pathogenic bacteria (two of each gram-negative and positive bacteria). (A) Lifespan assay. (B) Killing assay usingEscherichia coli O157:H7 EDL933. (C) Killing assay usingSalmonella Typhimurium SL1344. (D) Killing assay usingStaphylococcus aureus Newman. (E) Killing assay usingListeria monocytogenes EGD-e. OP50,E. coli ; 3201,L. rhamnosus ; LGG,L. rhamnosus . Statistical analysis was performed using the Kaplan–Meier method, and differences were considered significant when thep value was below 0.05 (*) and 0.01 (**). Survival statistic in lifespan assay compared with 3201:p = 0.0000 andp = 0.0000 for OP50 and LGG, respectively. Survival statistics in killing assay compared with 3201:Escherichia coli O157:H7 EDL933,p = 0.0114 andp = 0.5001 for OP50 and LGG, respectively;Salmonella Typhimurium SL1344,p = 0.0000 andp = 0.2732 for OP50 and LGG, respectively;Staphylococcus aureus Newman,p = 0.0000 andp = 0.0000 for OP50 and LGG, respectively;Listeria monocytogenes EGD-e,p = 0.0000 andp = 0.0496 for OP50 and LGG, respectively. Data are expressed as means ±SEM .
Subsequently, we investigated the effect of 3201 on the host’s immune response against four different kinds of foodborne pathogenic bacteria.
Evaluation of Phenotypic Change of C. elegans after Exposure to L. rhamnosus IDCC 3201
We measured the body size to evaluate whether feeding 3201 changed the phenotype of worms. The body size was measured after feeding with OP50, 3201, or LGG for 48 h. The results showed that the 3201 group demonstrated a significant increase in both length and width compared with both the OP50 group (
-
Figure 3. Evaluation of phenotypic change of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201. Worm size and locomotive activity ofC. elegans strainfer-15;fem-1 after 48 h exposure period with OP50, 3201, or LGG. (A) Length, (B) Width, (C) Peristaltic speed, and (D) Pumping rate. OP50,E. coli ; 3201,L. rhamnosus ; LGG,L. rhamnosus . Statistical analysis was performed using one-way ANOVA, and statistical significance was considered when thep values were below 0.05 (*), 0.01 (**), 0.001 (***), and 0.0001 (****). Statistics compared with 3201: length,p < 0.0001 andp < 0.0001 for OP50 and LGG, respectively; width,p < 0.0001 andp = 0.0028 for OP50 and LGG, respectively; peristaltic speed,p = 0.5366 andp = 0.6957 for OP50 and LGG, respectively; pumping rate,p < 0.0001 andp = 0.0002 for OP50 and LGG, respectively. Data are expressed as means ±SEM .
Transcriptomic Analysis of C. elegans after Exposure to L. rhamnosus IDCC 3201
We conducted a transcriptome analysis to evaluate the altered mechanism through feeding 3201 compared with OP50. We selected genes that showed a significant upregulation of more than a 2-fold change in
-
Table 1 . Transcriptomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201..Term Gene count % p valueInnate immune response 75 4.8 0.000 Glutathione metabolic process 18 1.1 0.000 Peptidoglycan catabolic process 8 0.5 0.000 Cell wall macromolecule catabolic process 9 0.6 0.000 Defense response to gram-positive bacterium 17 1.1 0.000 Lipid metabolic process 32 2 0.000 Proteolysis 42 2.7 0.000 Defense response to gram-negative bacterium 19 1.2 0.000 Dephosphorylation 13 0.8 0.000 Fatty acid beta-oxidation using acyl-CoA oxidase 5 0.3 0.000 Top 10 pathways closely associated with genes significantly upregulated more than 2.0 folds in
C. elegans strainfer-15;fem-1 following 48 h of exposure toL. rhamnosus 3201 compared withE. coli OP50..
-
Table 2 . Transcriptomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201..Group and gene Gene number Fold change p valueDescription C-type lectin-related clec-41 CELE_B0365.6 8.443320 0.000 C-type lectin clec-62 CELE_F35C5.5 2.348795 0.000 C-type lectin clec-86 CELE_C54D1.2 9.677637 0.000 C-type lectin clec-186 CELE_ZK896.7 6.833677 0.000 C-type lectin clec-187 CELE_ZK896.6 4.984468 0.000 C-type lectin Lysozyme-related lys-1 CELE_Y22F5A.4 3.182571 0.000 Lysozyme lys-2 CELE_Y22F5A.5 8.078466 0.000 Lysozyme lys-5 CELE_F22A3.6 12.360469 0.000 Lysozyme lys-7 CELE_C02A12.4 5.542044 0.000 Lysozyme lys-8 CELE_C17G10.5 2.706994 0.000 Lysozyme List of genes that were significantly upregulated by more than 2.0 folds in
C. elegans when fed withL. rhamnosus 3201 compared withE. coli OP50 and associated with the innate immune response pathway..
-
Figure 4. Transcriptomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201. The identification of KEGG pathways related to genes significantly upregulated by more than 3-fold inC. elegans when fed with 3201 compared with OP50. We utilized Cytoscape for the analysis.
Metabolomic Analysis of C. elegans after Exposure to L. rhamnosus IDCC 3201
The metabolomic analysis was performed to assess the impact of 3201 on the metabolite composition of
-
Figure 5. Metabolomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201. Change of metabolites composition inC. elegans strainfer-15;fem-1 after 48 h exposure period with OP50 or 3201. (A) PLS-DA. (B) Volcano plot. (C) Top 25 enriched heatmap. (D) The quantitative graph using the metabolites which changed more than 4.0 folds.
Discussion
According to the previous study, probiotics are known to exert their beneficial effects on the host by adhering to host epithelial cells, modulating the host’s immune response, and producing specific metabolites [26].
According to a previous study, the adhesive ability to the host’s gastrointestinal is one of the classical selection criteria for potential probiotic bacteria. This adhesive capability can lead to colonization, which can help improve immunomodulatory effects and stimulate the gut barrier and metabolic functions [32]. Our study revealed that 3201 exhibits a significantly higher adhesive ability than OP50 and shows no significant difference compared with LGG. In lifespan assay, 3201 significantly improved the longevity of
According to previous studies, the quality of food and some types of bacteria are known to affect the phenotype of worms such as showing better growth performance [34-36]. Therefore, we sought to assess the impact of feeding 3201 on various phenotypic traits. Measurements of body size, peristaltic speed, and pumping rate were employed as indicators of potential changes in the worms’ phenotype. The results revealed a significant increase in both length and width when
Our previous findings revealed that preconditioning with 3201 enhanced resistances against both gram-negative and gram-positive foodborne pathogens. Therefore, we speculated that 3201 may modulate specific gene expression related to the immune response. Transcriptome analysis using genes that significantly increased more than 2-fold in
The metabolome analysis revealed increased levels of fumaric acid and succinic acid, metabolites of the citric acid cycle, in
In conclusion, our study focused on evaluating the probiotic properties of
Acknowledgements
This study was supported by the Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, Forestry and Fisheries (IPET-321037-5) and by a National Research Foundation of Korea Grant funded by the Korean government (MEST) (NRF-2021R1A2C3011051) and “Cooperative Research Program for Agriculture Science & Technology Development (Project No. RS-2023-00225838)” Rural Development Administration, Republic of Korea.
Conflict of Interest
The authors have no financial conflicts of interest to declare.
Fig 1.
Fig 2.
Fig 3.
Fig 4.
Fig 5.
-
Table 1 . Transcriptomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201..Term Gene count % p valueInnate immune response 75 4.8 0.000 Glutathione metabolic process 18 1.1 0.000 Peptidoglycan catabolic process 8 0.5 0.000 Cell wall macromolecule catabolic process 9 0.6 0.000 Defense response to gram-positive bacterium 17 1.1 0.000 Lipid metabolic process 32 2 0.000 Proteolysis 42 2.7 0.000 Defense response to gram-negative bacterium 19 1.2 0.000 Dephosphorylation 13 0.8 0.000 Fatty acid beta-oxidation using acyl-CoA oxidase 5 0.3 0.000 Top 10 pathways closely associated with genes significantly upregulated more than 2.0 folds in
C. elegans strainfer-15;fem-1 following 48 h of exposure toL. rhamnosus 3201 compared withE. coli OP50..
-
Table 2 . Transcriptomic analysis of
C. elegans after exposure toLacticaseibacillus rhamnosus IDCC 3201..Group and gene Gene number Fold change p valueDescription C-type lectin-related clec-41 CELE_B0365.6 8.443320 0.000 C-type lectin clec-62 CELE_F35C5.5 2.348795 0.000 C-type lectin clec-86 CELE_C54D1.2 9.677637 0.000 C-type lectin clec-186 CELE_ZK896.7 6.833677 0.000 C-type lectin clec-187 CELE_ZK896.6 4.984468 0.000 C-type lectin Lysozyme-related lys-1 CELE_Y22F5A.4 3.182571 0.000 Lysozyme lys-2 CELE_Y22F5A.5 8.078466 0.000 Lysozyme lys-5 CELE_F22A3.6 12.360469 0.000 Lysozyme lys-7 CELE_C02A12.4 5.542044 0.000 Lysozyme lys-8 CELE_C17G10.5 2.706994 0.000 Lysozyme List of genes that were significantly upregulated by more than 2.0 folds in
C. elegans when fed withL. rhamnosus 3201 compared withE. coli OP50 and associated with the innate immune response pathway..
References
- Kim B, Kim K, Xu X, Lee H, Pathiraja D, Park DJ,
et al . 2023. Complete genome and two plasmids sequences ofLactiplantibacillus plantarum L55 for probiotic potentials.J. Anim. Sci. Technol. 65 : 1341-1344. - Morelli L, Capurso L. 2012. FAO/WHO guidelines on probiotics: 10 years later.
J. Clin. Gastroenterol. 46 : S1-S2. - Lee D, Goh TW, Kang MG, Choi HJ, Yeo SY, Yang J,
et al . 2022. Perspectives and advances in probiotics and the gut microbiome in companion animals.J. Anim. Sci. Technol. 64 : 197-217. - Yang S, Deng C, Li Y, Li W, Wu Q, Sun Z,
et al . 2022. Complete genome sequence ofLactiplantibacillus plantarum ST, a potential probiotic strain with antibacterial properties.J. Anim. Sci. Technol. 64 : 183-186. - Kim B, Meng Z, Xu X, Baek S, Pathiraja D, Choi IG,
et al . 2023. Complete genome sequence ofLimosilactobacillus fermentum JNU532 as a probiotic candidate for the functional food and feed supplements.J. Anim. Sci. Technol. 65 : 271-274. - Kim JY, Kim JY, Kim H, Moon EC, Heo K, Shim JJ,
et al . 2022. Immunostimulatory effects of dairy probiotic strainsBifidobacterium animalis ssp. HY8002 andLactobacillus plantarum HY7717.J. Anim. Sci. Technol. 64 : 1117-1131. - Kang AN, Mun D, Ryu S, Jae Lee J, Oh S, Kyu Kim M,
et al . 2022. Culturomic-, metagenomic-, and transcriptomic-based characterization of commensal lactic acid bacteria isolated from domestic dogs usingCaenorhabditis elegans as a model for aging.J. Anim. Sci. 100 : skac323. - Jaafar M, Xu P, Mageswaran UM, Balasubramaniam SD, Solayappan M, Woon J,
et al . 2024. Dairy based-LAB improved constipation via increasing fecal bulk and decreasing concentration of fecal threonine while preserving colonic goblet cell count.J. Anim. Sci. Technol. 66 : 178-203. - Kaletta T, Hengartner MO. 2006. Finding function in novel targets:
C. elegans as a model organism.Nat. Rev. Drug discov. 5 : 387-398. - Mahesh R, Ilangovan P, Nongbri D, Suchiang K. 2021. Probiotics interactions and the modulation of major signalling pathways in host model organism
Caenorhabditis elegans .Indian J. Microbiol. 61 : 404-416. - Clark LC, Hodgkin J. 2014. Commensals, probiotics and pathogens in the
Caenorhabditis elegans model.Cell. Microbiol. 16 : 27-38. - Choi H, Mun D, Ryu S, Kwak Mj, Kim BK, Park DJ,
et al . 2023. Molecular characterization and functionality of rumen-derived extracellular vesicles using aCaenorhabditis elegans animal model.J. Anim. Sci. Technol. 65 : 652-663. - Kang A, Kwak MJ, Lee DJ, Lee JJ, Kim MK, Song M,
et al . Dietary supplementation with probiotics promotes weight loss by reshaping the gut microbiome and energy metabolism in obese dogs. Microbiol.Spectr. 0 : e02552-02523. - Anderson SM, Pukkila-Worley R. 2020. Immunometabolism in
Caenorhabditis elegans .PLoS Pathog. 16 : e1008897. - Wan QL, Shi X, Liu J, Ding AJ, Pu YZ, Li Z,
et al . 2017. Metabolomic signature associated with reproduction-regulated aging inCaenorhabditis elegans .Aging (Albany NY) 9 : 447-474. - Kim H, Shin M, Ryu S, Yun B, Oh S, Park D-J,
et al . 2021. Evaluation of probiotic characteristics of newly isolated lactic acid bacteria from dry-aged Hanwoo beef.Food Sci. Anim. Resour. 41 : 468-480. - Park MR, Ryu S, Maburutse BE, Oh NS, Kim SH, Oh S,
et al . 2018. ProbioticLactobacillus fermentum strain JDFM216 stimulates the longevity and immune response ofCaenorhabditis elegans through a nuclear hormone receptor.Sci. Rep. 8 : 7441. - Kim Y, Mylonakis E. 2012.
Caenorhabditis elegans immune conditioning with the probiotic bacteriumLactobacillus acidophilus strain NCFM enhances gram-positive immune responses.Infect. Immun. 80 : 2500-2508. - Shen P, Kershaw JC, Yue Y, Wang O, Kim KH, McClements DJ,
et al . 2018. Effects of conjugated linoleic acid (CLA) on fat accumulation, activity, and proteomics analysis inCaenorhabditis elegans .Food Chem. 249 : 193-201. - Ryu S, Shin M, Yun B, Lee W, Choi H, Kang M,
et al . 2021. Bacterial quality, prevalence of pathogens, and molecular characterization of biofilm-producingStaphylococcus aureus from Korean dairy farm environments.Animals 11 : 1306. - Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data.
Bioinformatics 30 : 2114-2120. - Kim D, Langmead B, Salzberg SL. 2015. HISAT: a fast spliced aligner with low memory requirements.
Nat. Methods 12 : 357-360. - Liao Y, Smyth GK, Shi W. 2014. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features.
Bioinformatics 30 : 923-930. - Robinson MD, McCarthy DJ, Smyth GK. 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.
Bioinformatics 26 : 139-140. - Yoo J, Lee J, Zhang M, Mun D, Kang M, Yun B,
et al . 2022. Enhanced γ-aminobutyric acid and sialic acid in fermented deer antler velvet and immune promoting effects.J. Anim. Sci. Technol. 64 : 166-182. - Williams NT. 2010. Probiotics.
Am. J. Health Syst. Pharm. 67 : 449-458. - Oh YJ, Lee J, Lim SK, Kwon MS, Lee S, Choi SP,
et al . 2023. Complete genome sequence of probioticLactobacillus johnsonii 7409N31 isolated from a healthy Hanwoo calf.J. Anim. Sci. Technol. 65 : 890-893. - Hammes WP, Vogel RF. 1995. The genus lactobacillus, pp. 19-54.
In: The genera of lactic acid bacteria . Ed. Springer,. - Vasquez R, Oh JK, Song JH, Kang DK. 2022. Gut microbiome-produced metabolites in pigs: a review on their biological functions and the influence of probiotics.
J. Anim. Sci. Technol. 64 : 671-695. - Grompone G, Martorell P, Llopis S, González N, Genovés S, Mulet AP,
et al . 2012. Anti-inflammatoryLactobacillus rhamnosus CNCM I-3690 strain protects against oxidative stress and increases lifespan inCaenorhabditis elegans .PLoS One 7 : e52493. - Oh S, Park MR, Son SJ, Kim Y. 2015.
Lactobacillus rhamnosus strain GG influences on the longevity and immune response ofCaenorhabditis elegans to foodborne pathogens with stimulating microRNAs.FASEB J. 29 : 853. - Monteagudo-Mera A, Rastall RA, Gibson GR, Charalampopoulos D, Chatzifragkou A. 2019. Adhesion mechanisms mediated by probiotics and prebiotics and their potential impact on human health.
Appl. Microbiol. Biotechnol. 103 : 6463-6472. - Poupet C, Saraoui T, Veisseire P, Bonnet M, Dausset C, Gachinat M,
et al . 2019.Lactobacillus rhamnosus Lcr35 as an effective treatment for preventingCandida albicans infection in the invertebrate modelCaenorhabditis elegans : first mechanistic insights.PLoS One 14 : e0216184. - Shtonda BB, Avery L. 2006. Dietary choice behavior in
Caenorhabditis elegans .J. Exper. Biol. 209 : 89-102. - Avery L, You YJ.
C. elegans feeding. WormBook: The Online Review ofC. elegans Biology [Internet], 2018. - Kim YJ, Cho SB, Song MH, Lee SI, Hong SM, Yun W,
et al . 2022. Effects of differentBacillus licheniformis andBacillus subtilis ratios on nutrient digestibility, fecal microflora, and gas emissions of growing pigs.J. Anim. Sci. Technol. 64 : 291-301. - Pees B, Yang W, Kloock A, Petersen C, Peters L, Fan L,
et al . 2021. Effector and regulator: diverse functions ofC. elegans C-type lectinlike domain proteins.PLoS Pathog. 17 : e1009454. - Schulenburg H, Hoeppner MP, Weiner III J, Bornberg-Bauer E. 2008. Specificity of the innate immune system and diversity of Ctype lectin domain (CTLD) proteins in the nematode
Caenorhabditis elegans .Immunobiology 213 : 237-250. - O'Rourke D, Baban D, Demidova M, Mott R, Hodgkin J. 2006. Genomic clusters, putative pathogen recognition molecules, and antimicrobial genes are induced by infection of
C. elegans withM. nematophilum .Genome Res. 16 : 1005-1016. - Leippe M. 1999. Antimicrobial and cytolytic polypeptides of amoeboid protozoa-effector molecules of primitive phagocytes.
Dev. Comp. Immunol. 23 : 267-279. - Ciancio A. 2016. Defense and Immune Systems, pp. 205-239.
In: Invertebrate Bacteriology . Ed. Springer,. - Amador‐Noguez D, Yagi K, Venable S, Darlington G. 2004. Gene expression profile of long‐lived Ames dwarf mice and Little mice.
Aging Cell 3 : 423-441. - Perron JT, Tyson RL, Sutherland GR. 2000. Maintenance of tricarboxylic acid cycle kinetics in Brown-Norway Fischer 344 rats may translate to longevity.
Neurosc. Lett. 281 : 91-94. - Edwards CB, Copes N, Brito AG, Canfield J, Bradshaw PC. 2013. Malate and fumarate extend lifespan in
Caenorhabditis elegans .PLoS One 8 : e58345. - Zarse K, Schulz TJ, Birringer M, Ristow M. 2007. Impaired respiration is positively correlated with decreased life span in
Caenorhabditis elegans models of Friedreich Ataxia.FASEB J. 21 : 1271-1275. - Ryu S, Kim K, Cho D-Y, Kim Y, Oh S. 2022. Complete genome sequences of
Lactococcus lactis JNU 534, a potential food and feed preservative.J. Anim. Sci. Technol. 64 : 599-602. - Fernández MF, Boris S, Barbes C. 2003. Probiotic properties of human
lactobacilli strains to be used in the gastrointestinal tract.J. Appl. Microbiol. 94 : 449-455. - Dinev T, Beev G, Tzanova M, Denev S, Dermendzhieva D, Stoyanova A. 2018. Antimicrobial activity of
Lactobacillus plantarum against pathogenic and food spoilage microorganisms: a review.Bulg. J. Vet. Med. 21 . DOI:10.15547/bjvm.1084. - Shu Q, Gill HS. 2002. Immune protection mediated by the probiotic
Lactobacillus rhamnosus HN001 (DR20™) againstEscherichia coli O157: H7 infection in mice.FEMS Immunol. Med. Microbiol. 34 : 59-64.