- 1Department of Molecular Biology and Genetics, Democritus University of Thrace, Alexandroupolis, Greece
- 2Department of Medicine, Faculty of Health Sciences, Democritus University of Thrace, Alexandroupolis, Greece
Probiotic lactobacilli exert antimicrobial and antibiofilm activity through a plethora of mechanisms, including the production of inhibitory compounds and competition for nutrients and adhesion sites on the host mucosa. Their ability to interfere with pathogen-host interactions may mitigate the harmful effects of infection. Lactiplantibacillus pentosus L33 (L33) and Lactiplantibacillus plantarum L125 (L125) are two potential probiotic lactic acid bacteria (LAB) strains, previously shown to reduce pathogen viability and biofilm formation in vitro. The present study aims to investigate their ability to limit cell death induced by Staphylococcus aureus and Escherichia coli in the human colon adenocarcinoma cell line HT-29. To this end, we examined the protective effects of the two strains using in silico, in vitro and omic approaches, with Lacticaseibacillus rhamnosus GG (LGG) serving as a reference strain, due to its well-documented antimicrobial properties. Based on the findings of our study, direct contact of HT-29 cells with L125 for 4 h prior to the addition of S. aureus or E. coli prevented pathogen-induced cell death at rates comparable to LGG. In contrast, L33 failed to exert a protective effect. Moreover, L125 significantly reduced adherence of S. aureus to HT-29 cells, and the internalization capacity of both pathogens (>1.5 Log CFU/mL). Dual RNA-seq and protein microarrays were used to determine expression changes in L125 and host cells during co-incubation. L125 expressed high levels of adhesins and moonlighting proteins, homologous to those encoded by the pathogens. Pathways involved in pathogen adhesion and internalization, endocytosis, cell–cell and cell-extracellular matrix (ECM) adhesion, were downregulated in HT-29 cells. Finally, L125 reduced the secretion of various pro-inflammatory mediators. Our findings highlight the strain-specific protective effects of LAB against pathogen-induced cell death achieved through competitive exclusion and priming of host cell responses. Future studies will focus on elucidating the specific surface components of L125 involved in these events, paving the way for targeted interventions at the host-pathogen interface.
1 Introduction
Probiotics are live microorganisms that, when administered in sufficient quantities, provide health benefits to the host (Hill et al., 2014). A well-documented feature of probiotic strains is their ability to exert antimicrobial and antibiofilm activity against foodborne and clinically relevant pathogens (Silva et al., 2020). These effects are often attributed to lactic acid production and matrix acidification, and secretion of bacteriocins or other small molecules (Guo et al., 2020). In addition, cell-surface components may exert direct antimicrobial activity (Muscariello et al., 2020) or participate in competitive exclusion events (Tuo et al., 2018). Key components involved in these processes include pilins, adhesins and moonlighting proteins, which mediate adhesion to the ECM and host cell surface, such as fibronectin, collagen, and mucins (Sengupta et al., 2013; Jeffery, 2019). The close proximity of lactobacilli to host cells facilitates trans-kingdom signaling events, influencing critical cellular processes, such as cell death and survival, immune responses, and adhesion to the ECM and neighboring cells (Ye et al., 2023).
Staphylococcus aureus and Escherichia coli are two common bacteria that can cause localized or invasive infections in hospital and community settings (Ramos et al., 2020; Raineri et al., 2022). Pathogen adhesion and internalization are critical steps during the early stages of infection. S. aureus and E. coli code for a large repertoire of cell-surface proteins involved in these events, including microbial surface component recognizing adhesive matrix molecules (MSCRAMMs) (Josse et al., 2017), moonlighting adhesins (Furuya and Ikeda, 2011; Harvey et al., 2019), pilins and fimbriae (Croxen and Finlay, 2010). These proteins enable pathogens to attach to mucins and the ECM components, triggering the production of enterotoxins (Younes et al., 2016; Josse et al., 2017). Interactions between pathogens and host proteins such as fibronectin and actin, facilitate their internalization into non-phagocytotic host cells (Sinha and Fraunholz, 2010). This allows pathogens to evade immune surveillance, thus providing a fertile ground for the establishment of chronic infections (Croxen and Finlay, 2010; Foster et al., 2014). Limiting these interactions during the early phases of infection could effectively limit host colonization, and increase pathogen clearance. Therefore, commensals and probiotic microorganisms with the ability to interfere at the pathogen-host interface offer a promising complementary strategy for the prevention or management of life-threatening infections.
In the present study, we examined the ability of two potential probiotic strains, Lactiplantibacillus pentosus L33 and Lactiplantibacillus plantarum L125, isolated from fermented products (Pavli et al., 2016), to protect HT-29 human colon adenocarcinoma cells against S. aureus and E. coli infection. As we have previously shown, both strains exerted antimicrobial and antibiofilm effects against these pathogens in vitro (Kiousi et al., 2023). In particular, viable bacteria and cell-free culture supernatants (CFCS) reduced pathogen viability in suspension and inhibited biofilm formation (Kiousi et al., 2023). Furthermore, their WGS have been published and are available in public repositories (Stergiou et al., 2021; Tegopoulos et al., 2021). Here, we investigated the ability of L33 and L125 to prevent pathogen cytotoxicity with various experimental setups. Lacticaseibacillus rhamnosus GG was used as a reference strain due to its well-documented ability to limit adhesion, biofilm formation, virulence and pathogen-induced cell death in vitro and in vivo, through secreted metabolites and cell-surface molecules (De Keersmaecker et al., 2006; Spacova et al., 2020). LGG is, also, known to modulate pathogen clearance (Capurso, 2019). Trans-well systems were used to determine whether direct contact between the lactobacilli and host cells is required to prevent pathogen-induced cell death. Subsequently, the ability of the strains to adhere to HT-29 cells and inhibit pathogen adhesion and invasion was evaluated in vitro. WGS analysis was performed to identify proteins potentially involved in competitive exclusion events. Finally, dual RNA-seq and protein microarrays were employed to reveal genes and pathways involved in lactobacilli-host interactions and the associated protective phenotype.
2 Materials and methods
2.1 Bacterial strains and culture conditions
LAB strains L33 and L125 (Pavli et al., 2016) and LGG (DMSZ, Braunschweig, Germany) were cultivated in de Man, Rogosa and Sharp (MRS, Applichem, Darmstadt, Germany) broth at 37°C, under anaerobic conditions. S. aureus and E. coli were obtained from the microbial collection of the Laboratory of Clinical Microbiology, University Hospital of Alexandroupolis and were routinely cultured in Tryptic Soy Broth (TSB, Condalab, Madrid, Spain) at 37°C under aerobic conditions.
2.2 Human cell lines
HT-29 human colorectal adenocarcinoma cell line (ATCC, Manassas, VA, United States) was maintained in Roswell Park Memorial Institute GlutaMAX™ (RPMI)-1640 medium supplemented with 10% fetal bovine serum (FBS), 100 μg/mL streptomycin and 100 U/mL penicillin (all from Thermo Fisher Scientific, Waltham, MA, United States). For HT-29-bacteria co-incubations, a modified cell culture medium consisted of 10% FBS and 20 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES; Thermo Fisher Scientific) was used. Cells were incubated in a humidified atmosphere at 37°C, 5% CO2 under sterile conditions.
2.3 Cell and bacteria viability assays
The cytotoxic capacity of pathogens was initially evaluated using the Sulforhodamine B (SRB) (Invitrogen, Waltham, MA, United States) assay. HT-29 cells were seeded at a density of 7,500 cells per well in a 96-well plate (SPL Life Sciences, Pochon, South Korea) and incubated O/N in standard cell culture conditions. The following day, cells were washed with phosphate-buffered saline (PBS, Thermo Fischer Scientific) and incubated with viable pathogens or conditioned media (CM). For viable cell treatments, pathogens were added at a concentration of 107 CFU/mL for 1, 2 or 4 h. CM was prepared as follows: pathogens were incubated in RPMI-1640 medium supplemented with 10% FBS and 20 mM HEPES for 24 h. Cell pellets were discarded after centrifugation at 8,000 g for 5 min, and the supernatants were sterile filtered using a 0.2 μm filter (Merck, Rahway, NJ, United States). Cells were then treated with pathogen CM for 2 or 4 h. HT-29 cell viability was determined using a previously published protocol (Chondrou et al., 2018). Cell survival (%) is expressed as: [(Sample OD570 – Blank OD570)/(Control mean OD570 – Blank OD570)] × 100. Cells treated with standard culture medium served as the untreated control.
The cytotoxic capacity of pathogen CM was also evaluated using the HoloMonitor Live Imaging System (PHI, Boston, MA, United States). In this case, 2 × 105 cells per well were seeded into 24-well plates. The following day, cells were treated with undiluted, sterile-filtered pathogen CM. Cells treated with cell culture medium were included as an untreated control. Cell proliferation was assessed over a 24 h period with the Kinetic Cell Proliferation Assay software (PHI).
The potential protective effects of viable lactobacilli were evaluated against pathogen-induced toxicity, under two different conditions: (i) co-treatment of cells with lactobacilli (107 CFU/mL) and pathogens (108 CFU/mL), (ii) pre-treatment (2 or 4 h) of cells with lactobacilli (107 CFU/mL) prior to the addition of pathogens (108 CFU/mL). Pathogen challenge lasted for 4 or 2 h for S. aureus and E. coli, respectively. Cell survival was calculated as described above. Untreated cells served as a control, while cells exposed only to pathogens were used as a positive control.
2.4 Flow cytometry
Propidium iodide (PI) staining, and flow cytometry were used to examine membrane permeability and cell death. HT-29 cells were seeded at a density of 2 × 105 cells per well in 6-well plates. The following day, cells were pre-treated with lactobacilli (107 CFU/mL) for 4 h after which S. aureus (108 CFU/mL for 4 h) or E. coli (108 CFU/mL for 2 h) were added. After incubation, cells were collected by trypsinization and were centrifuged at 600 g for 5 min. Then, cell pellets were washed twice with PBS, and stained with 50 μg/mL PI (Biotium, Fremont, CA, United States) for 3 min. Untreated cells served as a negative control, and pathogen-only treated cells as a positive control. Flow cytometry was performed using the Attune NxT Flow Cytometer (Thermo Fischer Scientific) and results were analyzed with FlowJo V10 software (BD Biosciences, San Jose, CA, United States).
2.5 Trans-well system
Cell culture inserts were used to determine whether direct contact between L125 and HT-29 cells is required for protection against pathogen-induced cytotoxicity. Specifically, HT-29 cells were seeded at a density of 5 × 104 cells per well into 24-well plates. The following day, L125 cells (107 CFU/mL) were added to either side of a 0.4 μm trans-well insert (SPL Life Sciences) for 4 h. Then, pathogens were inoculated in the lower compartment, in direct contact with HT-29 cells. Cell survival was determined using the SRB assay, after 4 h challenge with S. aureus or 2 h challenge with E. coli.
2.6 Adherence and competitive exclusion assays
HT-29 were seeded in 24-well plates at a density of 4 × 105 cells per well and incubated until reaching 100% confluency, following a previously published protocol (Plessas et al., 2020). Lactobacilli at a concentration of 107 CFU/mL were then added for 2 or 4 h. After incubation, the monolayers were washed twice with sterile PBS, and cells were detached using 1× Trypsin (Thermo Fischer Scientific). The resulting cell suspensions were serially diluted in Ringer’s solution and plated onto MRS agar for enumeration. Plates were incubated at 37°C until visible colonies formed. Attached bacteria are expressed as Log CFU/mL. Lactobacilli-HT-29 interactions were further visualized using an inverted microscope (ZEISS, Jena, Germany).
For the competitive exclusion assay, HT-29 cells were seeded into 24-well plates (4 × 105 cells/well). The following day, cells were pretreated with lactobacilli for 4 h prior to pathogen challenge. Then, pathogens (108 CFU/mL) were added and incubated for 2 h. At the end of the incubation period, the monolayers were washed twice with sterile PBS and cells were disassociated with 1× Trypsin. The suspension was serially diluted in Ringer’s solution and plated onto TSA for colony enumeration of pathogens. The plates were incubated at 37°C, under aerobic conditions, until visible colonies formed. Attached bacteria were expressed as Log CFU/mL. Controls included pathogen-only treated cells and untreated cells to account for potential contamination.
2.7 Gentamicin protection assay
HT-29 cells were seeded in 24-well plates at a density of 5 × 104 cells per well. The following day, cells were pre-treated with lactobacilli at a concentration of 107 CFU/mL for 4 h. Pathogens were, then, added at a concentration of 108 CFU/mL for 1 h. Monolayers were washed twice with PBS and incubated in RPMI medium supplemented with 10% FBS and 100 μg/mL gentamicin (all from Thermo Fischer Scientific) for 1 h. Cells were then lysed using 1% (v/v) Triton-X (Thermo Fischer Scientific), the resulting suspension was serially diluted in Ringer’s solution and plated onto TSA plates for colony enumeration of internalized pathogens. Plates were incubated at 37°C, under aerobic conditions, until visible colonies formed. Internalized bacteria were expressed as Log CFU/mL.
2.8 Dual RNA-seq
To examine global changes in the expression of genes in L125 and HT-29 cells during co-incubation, dual RNA-seq experiments were performed. Specifically, HT-29 cells were seeded at a density of 8 × 105 cells in 100 mm plates and incubated O/N. The following day, L125 cells were added at a concentration of 107 CFU/mL for 4 h. HT-29 cells treated with cell culture medium were utilized as non-treated control samples. After treatments, plates were washed with PBS, and cells and attached bacteria were collected by trypsinization. Samples were centrifuged at 8,000 g for 5 min and incubated in a lysis buffer containing 1 M Τris-Cl pH 8, 0.5 M EDTA pH 8, 0.1% Triton-X and 100 mg/mL lysozyme, at 37°C for 30 min to maximize the efficacy of RNA isolation from the Gram-positive bacterium L125. Trizol reagent (Sigma-Aldrich, Saint Louis, MO, United States) was finally added, and RNA isolation was performed following manufacturer’s instructions. RNA quality was determined via agarose gel electrophoresis and spectrophotometrically (Thermo Fisher Scientific, NanoDrop 1000 Spectrophotometer). RNA concentration was measured with Qubit according to the Qubit™ RNA Broad Range protocol (Thermo Fisher Scientific). Dual RNA-seq (metatranscriptomics) was performed with Illumina NovaSeq X Plus (Strategy PE150). Bioinformatic analysis was performed on the Galaxy server (Abueg et al., 2024). The quality of the obtained reads was determined with FASTQC (v0.74; Andrews, 2010), and low-quality reads were discarded with Trimmomatic (v0.39; Bolger et al., 2014). HISAT2 (v2.2.1; Kim et al., 2019) was used for read assembly against the L125 or human (GRCh38/hg38) genomes. Gene counts and lengths were determined with featureCounts (v2.0.3; Liao et al., 2014) and differential expression with Limma (v3.25.1; Ritchie et al., 2015) and DESeq2 (v1.1.0; Love et al., 2014). Goseq (v1.50.0; Young et al., 2010) and EGSEA (v1.20.0; Alhamdoosh et al., 2017) were used to classify the genes into KEGG pathways and functional categories and for gene ontology analysis, respectively. Transcripts per million (TPM) were calculated for L125 transcripts with StringTie (v2.1.5; Pertea et al., 2015). Based on TPM values expression was characterized as high (>1,000 TPM), medium (10 < TPM < 1,000), low (0.5 < TPM < 10) or below cut-off (TPM <0.5).1
2.9 In silico analysis
The WGS of L125 (Accession number: JAKJPP000000000.1) was searched for S-layer proteins (SLAPS) and proteins containing cell-surface associated domains (LysM, SH3, WxL) and motifs (LPxTG and LPxTG-like) using EggNOG (Huerta-Cepas et al., 2019) and InterPro (Blum et al., 2021). Proteins containing signals for translocation to the bacterial cell surface were predicted using SignalP 6.0 (Teufel et al., 2022). Their localization and topology were also validated with PSORTb (v3.0; Yu et al., 2010) and DeepTMHMM (Hallgren et al., 2022), respectively. Putative surface proteins were selected for further analysis based on the presence of motifs and domains implicated in interactions with host cells and ECM, including mucin-binding domains (MuBPs), fibronectin-binding domains (FnBPs), collagen-binding domains (CnBPs) or proteins containing Ig-like folds and leucine-rich repeats (LRRs). Furthermore, known moonlighting proteins (Jeffery, 2019) were also included. Sequence homology between the annotated cell surface proteins and proteins encoded by S. aureus and E. coli was examined using Blastp (Camacho et al., 2009).
2.10 Docking study
The structure of the putative cell-surface associated proteins containing adhesin-related domains and motifs was predicted with ColabFold (Mirdita et al., 2022). Protein structures of human receptors and proteins participating in the interactions were downloaded from PDB or AlphaFold and used for docking experiments. These include: human fibrinogen chain A (FnA, 1BBR), fibrinogen fragment D (FnB, 1FZA), fibrinogen chain G (FnG, 1DUG), mucin 1 (MUC1, 1SM3), mucin 3b (MUC3b, AF-Q9H195), mucin 4 (MUC4, AF-Q99102), mucin 7 (MUC7, AF-Q8TAX7), mucin 13 (MUC13, AF-Q9H3R2), mucin 16 (MUC16, 7SA9), Toll-like receptor 2 (TLR2, AF-O606063), Toll-like receptor 4 (TLR4, AF-O00206), Toll-like receptor 5 (TLR5, AF-O60602), Toll-like receptor 6 (TLR6, AF-Q6Y2C9), carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1, AF-P13688), carcinoembryonic antigen-related cell adhesion molecule 7 (CEACAM7, Q14002), cadherin 1 (CDH1, AF-P12830) and collagen alpha-1(I) chain (COL1A1, AF-P02452). Docking experiments were performed on the HDOCK server (Yan et al., 2020).
2.11 Human cytokine array
High throughput determination of cytokines produced by HT-29 cells after co-incubation with L125 was performed using the Proteome Profiler Human Cytokine Array Kit (R&D Systems, Minneapolis, MN, United States). HT-29 cells were treated with the LAB strain at a concentration of 107 CFU/mL for 4 h. Then, culture supernatants were collected, centrifuged at 11,000 g for 5 min and stored at −80°C until further analysis. Supernatants from three independent experiments were pooled and analyzed using the protein microarray, following manufacturer’s instructions. Briefly, the membranes spotted with antibodies against 36 different cytokines, chemokines and acute phase proteins were blocked for 1 h while cell culture supernatants were mixed and incubated with the detection antibody cocktail for 1 h at RT. Then, the membranes were incubated with the sample/antibody mixtures overnight at 4°C on a rocking platform. The following day, membranes were washed three times, incubated with Streptavidin-HRP for 30 min at RT and washed again three times prior to chemiluminescence detection. Membranes were imaged using the ChemiDoc Imaging System (Bio-Rad, Hercules, CA, United States). Mean pixel density was calculated using a protein array analyzer plugin for ImageJ (National Institute of Health, United States).
2.12 Statistical analysis
For the statistical analysis of the experimental data, Student’s t-test was performed using GraphPad PRISM 9 (GraphPad Software Inc., CA, United States). All experiments were performed in triplicate unless otherwise stated. Results are represented as mean ± standard deviation. A p-value <0.05 was considered statistically significant.
3 Results
3.1 Staphylococcus aureus and Escherichia coli induce cell death of HT-29 cells in a time-dependent manner
First, we sought to determine the virulence effects of clinically isolated S. aureus and E. coli strains by their ability to induce cytotoxic effects on the colon adenocarcinoma cell line HT-29 at 1, 2 and 4 h (Figure 1). As shown in Figures 1A,C, a time-dependent reduction in cell viability was recorded for both pathogens. After 4 h-treatment of HT-29 cells with S. aureus or 2 h-treatment with E. coli, cell viability was reduced by >80% (p < 0.005). These timepoints were selected for subsequent experiments. These findings were further validated by PI staining and flow cytometry (Figures 1B,D). The effect of pathogen-derived CM on cell viability after 4 h (S. aureus, SACM) or 2 h (E. coli, EC CM), was initially investigated with the SRB cytotoxicity assay. A smaller reduction (p < 0.005) in cell survival was recorded (Figure 1E), compared to treatments with viable pathogens. These effects were further monitored for 24 h with the HoloMonitor live cell imaging system, using the Kinetic Cell Proliferation Assay software. It was shown that CM derived from either pathogen reduced cell confluency in a time dependent manner (Figure 1F).
Figure 1. Cytotoxic effects of S. aureus and E. coli against HT-29 cells. Co-incubation of HT-29 cells with 108 CFU/mL (A,B) S. aureus (SA) or (C,D) E. coli (EC) for 1, 2 or 4 h, limited cell viability as measured by (A,C) the SRB assay, (B,D) PI staining and flow cytometry. (E) S. aureus-derived CM (SA CM), and E. coli-derived CM (EC CM) limited cell viability after 4 h or 2 h of co-incubation, respectively. (F) The effect of pathogen-derived CM was monitored for 24 h using the HoloMonitor live cell imaging system. SA CM and EC CM decreased confluency in HT-29 cells in a time-dependent manner. The data presented are the mean ± standard deviation of three independent experiments. *p < 0.05 and **p < 0.005 compared to control (CTRL) untreated cells.
3.2 L125 pretreatment protects cells from pathogen-induced cytotoxicity
The ability of the potential probiotic LAB strains to limit pathogen-induced cell death was then assessed using various experimental designs (Figures 2A,3A). Initially, cells were co-incubated with lactobacilli and pathogens for 4 h (S. aureus-treated cells) or 2 h (E. coli-treated cells), and viability was determined with the SRB assay. No difference in cell survival (%) was observed between pathogen-only treated cells and those treated with any of the three lactobacilli (Figures 2B,C). Then, we investigated whether lactobacilli pretreatment could prevent pathogen-induced cell death (Figure 3A). As shown in Figure 3B, L125 and LGG significantly limited cell death, while L33 had no significant effect. These results were also confirmed through flow cytometry analysis (Figure 3C). It should be noted that lactobacilli treatment did not affect cell survival (Supplementary Figure S1).
Figure 2. Investigation of the protective effects of viable lactobacilli against pathogen-induced cytotoxicity in co-incubation experiments. (A) Schematic representation of the experimental design. Cells were simultaneously treated with 107 CFU/mL of lactobacilli and 108 CFU/mL of (B) S. aureus (SA) or (C) E. coli (EC) and cell viability was determined after 4 or 2 h, respectively, with the SRB assay. The data presented are the mean ± standard deviation of three independent experiments.
Figure 3. Investigation of the protective effects of viable lactobacilli against pathogen-induced cytotoxicity. (A) Schematic representation of the experimental design. (B) HT-29 cells were treated with lactobacilli at a concentration of 107 CFU/mL for 2 or 4 h before the addition of the pathogens (108 CFU/mL) (pretreatment). At the end of the incubation period, cell viability was determined using the SRB assay. (C) Determination of the protective effect of 4 h-pretreatment of the HT-29 cells against S. aureus (SA) or E. coli (EC) cytotoxicity via PI staining and flow cytometry. The data presented are the mean ± standard deviation of three independent experiments. *p < 0.05 and **p < 0.005 compared to control infected cells (SA or EC).
3.3 Direct contact of L125 with HT-29 cells confers protection against pathogen-induced cell death
A trans-well cell culture system was employed to determine whether direct contact of L125 with the epithelial cells is necessary for protection against pathogen-induced cytotoxicity (Figure 4A). Indeed, direct contact for 4 h prior to the addition of the pathogens resulted in a 2-fold reduction in cell death (p < 0.05). However, pretreatment with L125 in the upper compartment of the trans-well system, where no direct contact occurred, did not confer any protective effect (Figure 4B).
Figure 4. Use of a trans-well cell culture system to determine whether metabolites or cell surface molecules are responsible for the protective effects exerted by L125 against pathogen-induced cytotoxicity. (A) HT-29 cells at a density of 5 × 104 cells/well were seeded into the 24-well plates. The next day, L125 was added in the upper (separated) or lower (in contact with the cells) compartment for 4 h. Then, pathogens were added in the lower compartment for 4 (S. aureus, SA) or 2 h (E. coli, EC). (B) Cell viability was determined at the end of the incubation period using the SRB assay. The data presented are the mean ± standard deviation of three independent experiments. *p < 0.05 and **p < 0.005 compared to control infected cells (SA or EC).
3.4 L125 limits pathogen adhesion and invasion in HT-29 cells
To further explore the protective phenotype, we measured the capacity of L125 to adhere to HT-29 cells and limit pathogen attachment and invasion, using established microbiological assays. As shown in Figure 5A, L125 and LGG exhibited strong binding to HT-29 cells, while L33 demonstrated markedly lower adherence capacity (p < 0.05). Notably, L125 formed aggregates and adopted a net-like morphology over the HT-29 cells (Figure 5B). Pretreatment with L125 or LGG for 4 h reduced the attachment of S. aureus (p < 0.05), having no effect against E. coli (Figure 5C). L33 did not limit the adhesion ability of either pathogen, but did limit internalized S. aureus counts. L125 significantly inhibited the invasion of both pathogens in HT-29 cells, reducing counts by more than 1.5 Log CFU/mL (Figure 5D). A slight reduction in pathogen viability was also recorded at the same timepoints (Supplementary Figure S2).
Figure 5. Investigation of lactobacilli-HT-29-pathogen interactions. (A) Cells were treated with cell culture medium (CTRL) or lactobacilli at a concentration of 107 CFU/mL for 2 or 4 h. Attached bacteria are expressed as Log CFU/mL. (a, b) Significance was not achieved for variables with the same letter (p > 0.05). (B) Visualization of the lactobacilli-HT-29 interactions using an inverted optical microscope (ZEISS). (C) Competitive exclusion, and (D) inhibition of pathogen invasion was determined in cells pretreated with the lactobacilli for 4 h prior to the addition of pathogens. Then, cells were treated with the pathogens for 2 or 1 h for the competitive exclusion or gentamicin protection assay, respectively. Pathogen counts are expressed as Log CFU/mL. The data presented is the mean ± standard deviation of three independent experiments. *p < 0.05 and **p < 0.005 compared to cells incubated with the pathogens (SA or EC).
3.5 Lactobacilli-HT-29 cells interactions support the expression of adhesins and moonlighting proteins in L125
Dual RNA seq was utilized to analyze the gene expression changes in L125 and HT-29 cells during a 4 h co-incubation period. Specifically, for L125, a total of 108 genes exhibited high expression levels (>1,000 TPM), 2.568 genes displayed medium expression levels (10 < TPM < 1,000), and 497 genes showed low expression levels (0.5 < TPM < 10). Additionally, 68 genes were not expressed in neither of the two independent experiments (TPM: 0). The 108 highly expressed genes clustered into 14 KEGG functional categories, 23 KEGG pathways (Figure 6) and 16 clusters of orthologous groups (COGs) (Table 1). The most represented KEGG functional category was “genetic information processing” with “translation” being most prominent KEGG pathway and COG category.
Figure 6. Assignment of L125 transcripts with high expression levels (>1,000 TPM) in (A) KEGG functional categories, and (B) KEGG pathways using BlastKOALA.
To identify transcriptional changes in cell-surface exposed proteins, the WGS of L125 was, firstly, re-annotated in silico. Genome mining was performed to identify proteins potentially involved in the competitive exclusion phenotype, and the strain’s ability to limit pathogen internalization. By filtering the EggNOG output, 77 proteins containing motifs and domains indicative of cell-wall exposure were identified (Supplementary Table S1). Specifically, these proteins included: 17 proteins carrying the WxL domain, 19 proteins with a Gram-positive anchor (i.e., LPxTG or LPxTG-like motifs), 2 proteins with an S-layer homology (SLH) domain and 7 proteins with a SlaP (S-layer assembly protein) domain, 13 proteins containing the LysM domain and 9 proteins with the SH3 domain. Among these proteins: 7 possess Ig-like motifs (Big_2 and Big_3 groups), 5 carry LRR motifs, 3 carry collagen-binding domains (Cna_B), 13 contain mucin-binding domains (MucBP), 2 possess a fibronectin-binding domain (FbpA), 2 carry a FIVAR domain, and 1 has a domain belonging to the concanavalin A-like lectin/glucanases superfamily. Additionally, a protein belonging to the invasin/intimin bacterial adhesion mediator protein family was annotated. Out of these, 29 proteins were predicted to have extracellular localization or be exposed on the cell wall (Table 2), 15 were predicted to be localized to the cytoplasmic membrane, 6 to the cytoplasm, and 31 proteins to have unknown or multiple localization sites. Extracellular and cell wall proteins were selected for further analysis. Most predicted surface proteins were found to utilize a Sec/SPI signal recognized by signal peptidase I (Table 2), and one protein (LP125_RS11350) carries a YSIRK processing signal for transport outside the cell. Finally, 16 moonlighting proteins with adhesin function were identified in the L125 genome (Table 3).
Table 2. Putative cell-surface proteins encoded by L125 and their expression levels (TPM) during co-incubation with HT-29 cells.
Table 3. Moonlighting proteins with adhesin-like functions encoded by L125 and their expression levels (TPMs) during co-incubation with HT-29 cells.
To investigate possible exclusion events, the homology of L125 surface exposed proteins to those encoded by S. aureus or E. coli was analyzed in silico (Supplementary Tables S2, S3). Several proteins encoded by L125 contain domains with high alignment scores to those found in the two pathogens. Specifically, LP125_RS10075, LP125_RS11095 contain LPXTG anchor domains that show homology to those present in S. aureus (HDP6314765.1, HDM8613333.1), LP125_RS11095, LP125_RS11335, LP125_RS01150 and LP125_RS06780 carry MucBP domains that present high similarity to domains encoded by S. aureus (HDM8613333.1, EZX22480.1, MDU3924263.1). Furthermore, LRR repeats are contained in proteins encoded by L. plantarum L125 (LP125_RS02875, LP125_RS02875) and a BspA family leucine-rich repeat surface protein encoded by S. aureus (MDF4035866.1). Accordingly, LRR repeats found in L125 protein LP125_RS11350 share similarity to BspA family leucine-rich repeat surface protein (MBL0960876.1) and a bacterial Ig-like domain-containing protein (WP_142456104.1), respectively. Accordingly, LP125_RS11385, LP125_RS01150, LP125_RS07170 contain MucBP domains with similarity to MucBP domain-containing proteins encoded by E. coli (MBC8929979.1, MBC8921713.1). Furthermore, two proteins (LP125_RS01490 and LP125_RS09530) exhibited similarity to collagen-binding proteins (PPI97406.1, HCC5232179.1). The majority of moonlighting proteins annotated in the genome of L125, except for maltose phosphorylase (MapA, LP125_RS10155 and LP125_RS06445), share high sequence similarity with proteins encoded by S. aureus. On the contrary, L125 proteins generally exhibit lower query cover and/or similarity with E. coli proteins (Table 4).
Table 4. Sequence identity of moonlighting proteins encoded by L125 to proteins encoded by S. aureus and E. coli.
Concerning the expression at the transcriptional level, most of the putative cell surface proteins exhibited medium expression levels (100 < TPM < 1,000), with the exception of two extracellular proteins LP125_RS14880 and LP125_RS09840, which had high expression levels (>1,000 TPMs). No transcripts were detected for LP125_RS11935 and very low expression levels were recorded for LP125_RS08445 (TPM < 0.5) (Table 3). As expected, several of the moonlighting adhesins exhibited particularly high expression, reflecting their crucial role in bacterial homeostasis (Table 3). The highest expression levels were recorded for EF-Tu (LP125_RS08655), GAPDH (LP125_RS02930), and Hsp60 (LP125_RS03215). Notably, no transcripts were measured for ornithine carbamoyltransferase (LP125_RS10650).
The ability of the annotated, putative cell-exposed proteins to interact with host proteins was determined with in silico docking experiments. To this end, docking partners were determined based on the presence of domains and motifs involved in host–microbe interactions (Table 2). As shown in Figures 7, 8 proteins from L125 demonstrated the lowest docking scores (<−300) when interacting with host receptors and ECM components indicating strong interaction potential. More specifically, the LRR or Ig-like domains contained in LP125_RS10075, LP125_RS11335, LP125_RS12380 and LP125_RS04395 interact with CEACAM1 with a high confidence score (>95%). Accordingly, LP125_RS10075 and LP125_RS11350 interact with FnB (confidence score >95%) and LP125_RS12380 with FnA (confidence score 95.3). Furthermore, the MucBP and LRR domains present in LP125_RS11350 interact with MUC1, 4 and MUC16, as well as with LRR domains present in TLR2, 4, 5, 6 with a confidence score >97%. Additionally, the Ig-like folds of LP125_RS07170 demonstrated interaction with MUC4, TLR4 and 5 with a confidence score >94%. Interestingly, the highly expressed LP125_RS14880 and LP125_RS09840 proteins, which do not contain specific adhesin-related domains or motifs, may still contact extracellular domains of TLR4 and −5. Moreover, LP125_RS09840 may also form interactions with MUC1 and LP125_RS14880 with MUC3b and 7, as indicated by the low docking scores and high confidence scores (Figures 7A,B). Protein–protein interactions of partners with the lowest free energy are presented in Figure 7C.
Figure 7. Docking analysis of putative cell surface proteins encoded by L125 with host cell receptors and ECM proteins. (A) Docking scores, and (B) confidence scores derived from in silico docking experiments. (C) De novo docking of L125 cell surface proteins that present the highest scores with host receptors and ECM proteins. The 3D structure of the annotated proteins was predicted using Collabfold, and host protein structures were derived from PDB or AlphaFold. Docking was performed with HDOCK and the interactions of the proteins presenting the highest docking and confidence were visualized using PyMoL.
Figure 8. Differential expression analysis in cells treated with L125 for 4 h. (A) MD and volcano plots indicating differentially expressed genes in cells treated with L125 (L125-Tr) compared to control, untreated cells (CTRL): 600 genes were upregulated (red dot), 1,137 were significantly downregulated (blue dot), and no change was recorded for 26,658 genes (grey dot). (B) Heatmap of the top differentially expressed (topDEs) genes. Z-scores were calculated in rows. Top over-represented categories in (C) biological processes, and (D) molecular function determined with Goseq.
3.6 Microbe-host interactions alter gene expression in HT-29 cells
RNA seq was utilized to explore the capacity of L125 to prime antimicrobial responses in the host cell. As shown in Figure 8, treatment of cells with L125 for 4 h led to the upregulation of 600 genes and the downregulation of 1,137 genes in HT-29 cells, while 26,658 genes remained unaffected. The top over-represented categories in “biological function” and in “molecular function” are visualized in Figures 8C,D. Interestingly, L125 did not significantly modulate immune-related pathways, including TLR- or NOD- signaling cascades, nor did it affect the production of cytokines and chemokines (Figure 9A). At the protein level, L125 was shown to limit the secretion of immunological markers in pooled cell culture supernatants, collected from three independent experiments (Figures 9B–D). Pathway analysis revealed that L125 primarily exerted inhibitory effects (Supplementary Table S4), negatively impacting bacterial invasion (Figure 10), adherens junction (Supplementary Figure S3), and endocytosis (Supplementary Figure S4). These phenotypes are mediated by central processes related to the regulation of actin cytoskeleton, internalization and vacuole formation. In this context, genes involved on the mitotic spindle assembly pathway and G2/M progression were also downregulated, further emphasizing L125’s regulatory role in cellular processes.
Figure 9. Investigation of the effect of L125 on the expression of immune-related genes in HT-29 cells. (A) The expression levels of cytokines and chemokines were determined at the transcriptional level with RNA-seq. Heatmaps of normalized counts of two independent experiments calculated with Limma-voom were constructed on GENE-E. (B,C) The effect of lactobacilli treatment (4 h) (L125-Tr) on the expression of 36 immunological markers in HT-29 cells was estimated using the Proteome Profiler Human Cytokine Array Kit. The supernatants of three independent experiments were pooled and utilized in the assay. Untreated cells were used as control (CTRL). (D) Mean pixel density of spots was calculated using the protein array a protein array analyzer plugin for ImageJ. Pairs of spots in the upper-left, upper-right and down-left corners are positive controls.
Figure 10. Pathway enrichment analysis of cells treated with L. plantarum L125 vs. untreated cells. Pathway “Bacterial invasion of epithelial cells” (hsa05100) was downregulated. Pathway enrichment analysis was performed with EGSEA and the visualization of KEGG pathways with Pathview.
4 Discussion
Probiotic LAB employ various direct and indirect antimicrobial mechanisms to interfere with the host-pathogen interface, thereby limiting pathogen viability, adhesion and biofilm formation. These mechanisms include the production of antimicrobial metabolites and bacteriocins, competitive exclusion and competition for nutrients and resources (Arqués et al., 2015). Concomitantly, certain strains can also modulate cellular signaling to promote pathogen clearance (Llewellyn and Foey, 2017). S. aureus and E. coli are pathogens of clinical significance, that can invade the epithelial barrier and cause serious systemic and blood-borne infections. The production of membrane-damaging toxins (i.e., hemolysins, leukotoxins, leucocidin) and toxins that interfere with receptor function (i.e., enterotoxins and enterotoxin-like toxins) is upregulated by adhesion and biofilm formation in the host niche, supporting pathogen internalization and dissemination (Younes et al., 2016; Josse et al., 2017). Therefore, limiting these early interactions is a key to preventing the deleterious effects of infection.
We have previously shown that L. pentosus L33 and L. plantarum L125 possess potent antimicrobial and antibiofilm potential against S. aureus and E. coli (Kiousi et al., 2023). In this study, we aimed to elucidate the potential of these strains to inhibit infection and prevent loss of cell viability. Based on our findings, L125 was capable of limiting pathogen-induced cytotoxicity in HT-29 cells, when in direct contact for 4 h prior to S. aureus or E. coli addition. These inhibitory effects induced by L125 were comparable to those observed with LGG, a strain known for its ability to inhibit S. aureus-induced cytotoxicity (Spacova et al., 2020; El-Chami et al., 2022). Notably, this is the first study to report that LGG can also prevent E. coli-induced cell death in HT-29 cells. In contrast, pretreatment with L33 did not prevent pathogen-induced cell death, but significantly decreased internalized S. aureus counts. We previously demonstrated that L33 promotes adhesion of S. aureus in co-incubation with HT-29 cells (Kiousi et al., 2023), suggesting that this phenotype could be attributed to steric hindrance or other specific interactions. Additionally, the cytoprotective effects of L125 and LGG were also recorded in the melanoma cell line A375 (data not shown). Similar direct contact requirements have been reported for other potential probiotic strains, such as L. jensenii RC-28 or L. reuteri RC-14, which limited adhesion, pathogen-induced inflammation and production of S. aureus toxin TSST-1 in the vaginal epithelial cell line VK2-E6E7 (Younes et al., 2016). Likewise, L. fermentum 8,711 increased the viability of Caco-2 infected with methicillin-resistant S. aureus (Jayashree et al., 2018). In this context, we hypothesized that the incubation period allows for the establishment of a protective net that limits pathogen adhesion and invasion and/or primes host responses towards pathogen clearance. To test this hypothesis, we employed in silico and in vitro approaches.
First, the capacity of L125 to exclude pathogen adhesion and invasion was examined, after a 4 h pretreatment period. A modest, yet significant, reduction in S. aureus attachment to HT-29 cells was recorded (~0.2 Log CFU/mL). Νotably, L125 demonstrated a higher capacity to limit invasion of both S. aureus and E. coli in HT-29 cells by more than 1.5 Log CFU/mL. Subsequently, we sought to investigate potential lactobacilli-host and lactobacilli-pathogen interactions, in silico, by re-annotating the genome of L125. The predicted surface proteins of L125 include those anchored on the cell wall through covalent interactions (LPxTG and LPxTG-like motifs) or non-covalent interactions (LysM, SH3 and WXL domains), contain domains involved in adhesion to host components like mucins (MucBPs) and fibronectin (FbpA), as well as interactions with TLRs via LRR domains and CEACAMs through Ig-like domains. Specifically, for S. aureus, which showed reduced cell attachment and internalization in HT-29 cells pretreated with L125, these motifs and domains are conserved in MSCRAMMs and utilized during host colonization and epithelial invasion. Previous studies have showed that LGG exerts cytoprotective effects via pilli encoded by the spaCBA cluster, which block S. aureus attachment to keratinocytes (Spacova et al., 2020). Similarly, moonlighting proteins with adhesin function, enolase, GAPDH, EF-Tu, and TpiA have been found to form transient interactions with host proteins, like fibrinogen, fibronectin, and cytoskeletal proteins (Jeffery, 2019). Among the annotated proteins in this study, only moonlighting adhesins exhibited significant similarity to proteins encoded by the two pathogens, alluding to their central role in mutual exclusion phenomena. In this context, enolase and TpiA encoded by LGG were found to inhibit S. aureus binding and cytotoxic effects on keratinocytes, in a dose-dependent manner (El-Chami et al., 2022). In addition to genome analysis, we examined the expression levels of adhesins and moonlighting proteins with adhesin function in L125 during co-incubation with HT-29 cells. As anticipated, moonlighting proteins participating in crucial housekeeping functions were highly expressed. Most adhesins were expressed in medium levels, with the exception of two extracellular proteins (LP125_RS14880, LP125_RS09840), which were expressed at high levels. These findings suggest that during the pre-incubation period these adhesins and moonlighting proteins may help anchor L125 to the HT-29 cell surface, promoting competitive exclusion via specific interactions and steric hindrance. Interestingly, we previously showed that L125 presents high co-aggregation capacity with S. aureus, which could further contribute to its ability to limit pathogen adhesion and invasion (Kiousi et al., 2023). These interactions may play a key role in the observed protective effects of L125 in cellular models.
Bacterial cell-surface proteins with adhesion capacity can participate in trans-kingdom signaling events. In this context, lactobacilli are shown to induce strain-specific immunomodulatory effects that are mainly mediated by cell-surface MAMPs or excreted metabolites (Wells, 2011). Several annotated proteins in the genome of L125 were shown to possess LRR domains that can interact with TLRs. TLR-2 activation by L. plantarum strains has been linked to strengthening the gut barrier by upregulating and enhancing trafficking of zona occludens-1 (ZO-1) and occludin (Karczewski et al., 2010), while also stimulating pro-inflammatory responses (Fidanza et al., 2021). On the other hand, L. plantarum glycolipids were shown to trigger TLR-2/TLR-6 signaling, leading to tolerogenic IL-10 responses and the induction of regulatory T (Treg) populations (Karczewski et al., 2010). MUCs, also, play a role in regulating cellular survival and immune responses. For instance, the formation of the MUC4-ERBB2-ERBB3-NRG1 complex can suppress apoptosis (Carraway et al., 2002), while MUC13 has been identified as an oncogenic glycoprotein that promotes cancer cell growth (Sheng et al., 2016; Pang et al., 2022). Regarding the immunomodulatory capacity of these proteins, MUC1 and MUC16 were found to mask TLR signaling, resulting in the reduction of effector T cell functions, therefore enhancing cancer cell immune evasion and metastasis (Bhatia et al., 2019). Pro-inflammatory cytokines participate in pathogen clearance via the induction of local immune response and modulate the production of antimicrobial peptides and tight junction proteins from the epithelium, contributing to gut homeostasis. Indeed, the IL-1 family cytokines (i.e., IL-1α, IL-1β and IL-18) have been shown to upregulate transcription of defensins, iron-sequester lipocalin, as well as psoriasin and calprotectin (Kolls et al., 2008). Meanwhile, IL-8 and MCP-1 can exhibit direct antimicrobial effects (Yount et al., 2006). Moreover, short-term exposure to IL-10, IL-4, IFN-γ, IL-39γ, IL-6 and IL-22, and low levels of TNF-α can stimulate tissue regeneration and proliferation (Andrews et al., 2018). In this study, dual RNA-seq and protein microarrays were used to study immunomodulatory events that may contribute to pathogen clearance and resistance to pathogen-induced cell death. At the transcriptional level, L125 did not significantly alter cytokine or chemokine levels, but decreased the secretion of specific immunological markers in the cell culture medium. This aligns with our previous findings that immune response induction by potential probiotic lactobacilli in non-professional antigen-presenting cells requires prolonged incubation (Chondrou et al., 2020), suggesting that the short incubation period in this study may not have been sufficient for significant changes in immunological marker levels. Finally, L125 did not affect the production of antimicrobial peptides, as cathelicidin and defensin levels remained unchanged between treated and untreated cells. In the future, metabolomic and proteomic approaches will be used to further characterize the bidirectional L125-host interactions and validate the mechanisms involved in the recorded protective phenotype. In this context, Q Exactive-Based Quantitative Proteomics shown that Lactobacillus mucosae strain LM1 could affect tight junction assembly, and cellular and metabolic processes in IPEC-J2 intestinal epithelial cells (Pajarillo et al., 2017). At the same time, similar approaches have been applied to probe the Salmonella spp. (Birk et al., 2024), Pseudomonas auruginosa (Liu et al., 2024), and Acinetobacter baumanii (Li et al., 2021) - host crosstalk, providing significant insight into the early and late stages of infection. The application of proteomics in the context of the lactobacilli-host-pathogen interface can provide novel insight into the capacity of lactobacilli to limit the detrimental effects of infection (Stastna, 2024).
Pathway enrichment analysis revealed a significant downregulation of key signaling pathways involved in endocytosis, adherens junction formation, and the attachment and invasion of epithelial cells by pathogens. To the best of our knowledge, this is the first study to use dual RNA-seq to investigate the priming effects of lactobacilli on cellular function, in co-culture experiments. The affected modules are predominantly exploited by pathogens that colonize and internalize into host epithelia. In this context, S. aureus utilizes Fn-binding proteins to engage integrin α5β1, re-organize the actin cytoskeleton and generate focal contacts by regulating protein tyrosine kinase signaling and focal adhesion kinase activity (Hauck and Ohlsen, 2006). Similar integrin-mediated interactions that result in internalization are also utilized by Streptococcus spp. (Wang et al., 2007), Yersinia spp. (Matsumoto and Young, 2009), and certain E. coli pathotypes (Pokharel et al., 2023). Additionally, E. coli employs proteins with Ig domains to interact with host cell surface molecules like CEACAMs (van Sorge et al., 2021). Indeed, CEACAM-1 and -6 binding is used by pathogens to invade target cells (Pakbin et al., 2021). Several lactobacilli strains have been shown to reduce pathogen internalization, including L. salivarius W24, L. rhamnosus W71 (Campana et al., 2017), L. gasseri LG-7528, L. crispatus LCR-A21, L. paracasei LPP-A16 and L. rhamnosus LR-B5 (Poimenidou et al., 2023), and L. helveticus R0052 (Wine et al., 2009), in studies utilizing the gentamicin protection assay. Curiously, L125 upregulated the production Mad2L2, a protein involved in cell cycle arrest, an effect previously reported for Shigella, which exploits the epithelial renewal mechanisms to promote perseverance in the intestinal niche (Iwai et al., 2007).
Overall, our findings highlight the strain-specific ability of lactobacilli to enhance epithelial resistance against pathogenic insults via interference with host-pathogen interactions and regulation of pathogen internalization pathways (Figure 11). Future studies will focus on the identification of the specific surface components of L. plantarum L125 that orchestrate these events. Understanding these novel interactions may provide deeper insights into microbe-host interactions under homeostatic conditions and offer new strategies for targeting the host-pathogen interface.
Figure 11. Schematic representation of a putative mechanism by which L. plantarum L125 prevents pathogen-induced cell death. L125 cell surface is decorated by C-terminal anchored proteins containing LPxTG, WxL, SH3 or LysM motifs, S-layer proteins carrying the SLH domain, and moonlighting proteins with adhesin function. Adhesins form covalent bonds with cell wall or membrane components, while moonlighting proteins form reversible interactions based on their charge and hydrophobicity. L125 pretreatments resulted in the formation of a protective layer on epithelial cells, preventing the invasion of S. aureus and E. coli. At the same time, L125 downregulated pathways involved in pathogen adhesion and internalization, endocytosis, cell–cell adherence, and actin cytoskeleton formation.
5 Conclusion
Probiotics possess several mechanisms to limit pathogen viability and infectivity. In this study, the cytoprotective effects of two potential probiotic LAB strains against S. aureus and E. coli-induced cytotoxicity were examined using in vitro and in silico approaches. Among the strains tested, L125 exhibited a strong ability to limit cell death via direct contact with HT-29 cells for 4 h prior to pathogen exposure. Concomitantly, L125 significantly reduced pathogen internalization (>1.5 log reduction) and limited S. aureus attachment. Genome annotation of L125 revealed cell surface-associated and extracellular proteins likely involved in competitive exclusion events. Finally, microbe-host interactions triggered the expression of proteins with adhesin function in L125, and downregulated pathways related to endocytosis, adherence and pathogen internalization in epithelial cells. Future studies will focus on the identification and characterization of bacterial components that mediate the recorded cytoprotective effects.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, JAIGOE000000000.1; https://www.ncbi.nlm.nih.gov/, PRJNA1162724.
Ethics statement
Ethical approval was not required for the studies on humans in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used.
Author contributions
DK: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. MP: Formal analysis, Methodology, Supervision, Validation, Writing – original draft, Writing – review & editing. AP: Formal analysis, Supervision, Validation, Writing – original draft, Writing – review & editing. AG: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The research project was supported by the Hellenic Foundation for Research and Innovation (H.F.R.I.) under the “1st Call for H.F.R.I. Research Projects to support Faculty Members and Researchers and the Procurement of High-Cost Research Equipment” (Project Number: HFRIFM17C3-2007).
Acknowledgments
The authors acknowledge the support of the Biomedical Data Science and Bioinformatics Facility of the Department of Molecular Biology and Genetics, Democritus University of Thrace and the Institute of Technology of Agricultural Products, Hellenic Agricultural Organization DIMITRA (Likovrisi, Attiki, Greece), for kindly providing the strains L. pentosus L33 and L. plantarum L125. Figures 2A, 3A, 4A, 11 were created with Biorender.com.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1501119/full#supplementary-material
Footnotes
1. ^EMBL-EBI, Expression Atlas. Available on: https://www.ebi.ac.uk/gxa/FAQ.html.
References
Abueg, L. A. L., Afgan, E., Allart, O., Awan, A. H., Bacon, W. A., Baker, D., et al. (2024). The Galaxy platform for accessible, reproducible, and collaborative data analyses: 2024 update. Nucleic Acids Res. 52, W83–W94. doi: 10.1093/NAR/GKAE410
Alhamdoosh, M., Law, C. W., Tian, L., Sheridan, J. M., Ng, M., and Ritchie, M. E. (2017). Easy and efficient ensemble gene set testing with EGSEA. F1000Res 6:2010. doi: 10.12688/F1000RESEARCH.12544.1
Andrews, S. (2010), Babraham Bioinformatics—FastQC A Quality Control tool for High Throughput Sequence Data. Available at: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/. (Accessed August 10, 2024)
Andrews, C., McLean, M. H., and Durum, S. K. (2018). Cytokine tuning of intestinal epithelial function. Front. Immunol. 9:368738. doi: 10.3389/fimmu.2018.01270
Arqués, J. L., Rodríguez, E., Langa, S., Landete, J. M., and Medina, M. (2015). Antimicrobial activity of lactic acid bacteria in dairy products and gut: effect on pathogens. Biomed. Res. Int. 2015:584183. doi: 10.1155/2015/584183
Bhatia, R., Gautam, S. K., Cannon, A., Thompson, C., Hall, B. R., Aithal, A., et al. (2019). Cancer-associated mucins: role in immune modulation and metastasis. Cancer Metastasis Rev. 38, 223–236. doi: 10.1007/s10555-018-09775-0
Birk, M. S., Walch, P., Baykara, T., Sefried, S., Amelang, J., Buerova, E., et al. (2024). Salmonella infection impacts host proteome thermal stability. Eur. J. Cell Biol. 103:151448. doi: 10.1016/j.ejcb.2024.151448
Blum, M., Chang, H. Y., Chuguransky, S., Grego, T., Kandasaamy, S., Mitchell, A., et al. (2021). The InterPro protein families and domains database: 20 years on. Nucleic Acids Res. 49, D344–D354. doi: 10.1093/NAR/GKAA977
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/BIOINFORMATICS/BTU170
Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC Bioinformatics 10:421. doi: 10.1186/1471-2105-10-421
Campana, R., van Hemert, S., and Baffone, W. (2017). Strain-specific probiotic properties of lactic acid bacteria and their interference with human intestinal pathogens invasion. Gut Pathog. 9:12. doi: 10.1186/S13099-017-0162-4
Capurso, L. (2019). Thirty years of Lactobacillus rhamnosus GG: a review. J. Clin. Gastroenterol. 53, S1–S41. doi: 10.1097/MCG.0000000000001170
Carraway, K. L., Perez, A., Idris, N., Jepson, S., Arango, M., Komatsu, M., et al. (2002). Muc4/sialomucin complex, the intramembrane ErbB2 ligand, in cancer and epithelia: to protect and to survive. Prog. Nucleic Acid Res. Mol. Biol. 71, 149–185. doi: 10.1016/S0079-6603(02)71043-X
Chondrou, P., Karapetsas, A., Kiousi, D. E., Tsela, D., Tiptiri-Kourpeti, A., Anestopoulos, I., et al. (2018). Lactobacillus paracasei K5 displays adhesion, anti-proliferative activity and apoptotic effects in human colon cancer cells. Benef. Microbes 9, 975–983. doi: 10.3920/BM2017.0183
Chondrou, P., Karapetsas, A., Kiousi, D. E., Vasileiadis, S., Ypsilantis, P., Botaitis, S., et al. (2020). Assessment of the immunomodulatory properties of the probiotic strain Lactobacillus paracasei K5 in vitro and in vivo. Microorganisms 8:709. doi: 10.3390/MICROORGANISMS8050709
Croxen, M. A., and Finlay, B. B. (2010). Molecular mechanisms of Escherichia coli pathogenicity. Nat. Rev. Microbiol. 8, 26–38. doi: 10.1038/NRMICRO2265
De Keersmaecker, S. C. J., Verhoeven, T. L. A., Desair, J., Marchal, K., Vanderleyden, J., and Nagy, I. (2006). Strong antimicrobial activity of Lactobacillus rhamnosus GG against Salmonella typhimurium is due to accumulation of lactic acid. FEMS Microbiol. Lett. 259, 89–96. doi: 10.1111/J.1574-6968.2006.00250.X
El-Chami, C., Choudhury, R., Mohammedsaeed, W., McBain, A. J., Kainulainen, V., Lebeer, S., et al. (2022). Multiple proteins of Lacticaseibacillus rhamnosus GG are involved in the protection of keratinocytes from the toxic effects of Staphylococcus aureus. Front. Microbiol. 13:875542. doi: 10.3389/fmicb.2022.875542
Fidanza, M., Panigrahi, P., and Kollmann, T. R. (2021). Lactiplantibacillus plantarum-nomad and ideal probiotic. Front. Microbiol. 12:712236. doi: 10.3389/fmicb.2021.712236
Foster, T. J., Geoghegan, J. A., Ganesh, V. K., and Höök, M. (2014). Adhesion, invasion and evasion: the many functions of the surface proteins of Staphylococcus aureus. Nat. Rev. Microbiol. 12, 49–62. doi: 10.1038/nrmicro3161
Furuya, H., and Ikeda, R. (2011). Interaction of triosephosphate isomerase from Staphylococcus aureus with plasminogen. Microbiol. Immunol. 55, 855–862. doi: 10.1111/J.1348-0421.2011.00392.X
Guo, X., Chen, J., Sun, H., Luo, L., Gu, Y., Yi, Y., et al. (2020). Mining, heterologous expression, purification and characterization of 14 novel bacteriocins from Lactobacillus rhamnosus LS-8. Int. J. Biol. Macromol. 164, 2162–2176. doi: 10.1016/J.IJBIOMAC.2020.08.067
Hallgren, J., Tsirigos, K. D., Pedersen, M. D., Almagro Armenteros, J. J., Marcatili, P., Nielsen, H., et al. (2022). DeepTMHMM predicts alpha and beta transmembrane proteins using deep neural networks. bioRxiv. Available at: https://doi.org/10.1101/2022.04.08.487609. [Epub ahead of preprint]
Harvey, K. L., Jarocki, V. M., Charles, I. G., and Djordjevic, S. P. (2019). The diverse functional roles of elongation factor Tu (EF-Tu) in microbial pathogenesis. Front. Microbiol. 10:2351. doi: 10.3389/fmicb.2019.02351
Hauck, C. R., and Ohlsen, K. (2006). Sticky connections: extracellular matrix protein recognition and integrin-mediated cellular invasion by Staphylococcus aureus. Curr. Opin. Microbiol. 9, 5–11. doi: 10.1016/J.MIB.2005.12.002
Hill, C., Guarner, F., Reid, G., Gibson, G. R., Merenstein, D. J., Pot, B., et al. (2014). The International Scientific Association for Probiotics and Prebiotics consensus statement on the scope and appropriate use of the term probiotic. Nat. Rev. Gastroenterol. Hepatol. 11, 506–514. doi: 10.1038/nrgastro.2014.66
Huerta-Cepas, J., Szklarczyk, D., Heller, D., Hernández-Plaza, A., Forslund, S. K., Cook, H., et al. (2019). eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314. doi: 10.1093/NAR/GKY1085
Iwai, H., Kim, M., Yoshikawa, Y., Ashida, H., Ogawa, M., Fujita, Y., et al. (2007). A bacterial effector targets Mad2L2, an APC inhibitor, to modulate host cell cycling. Cell 130, 611–623. doi: 10.1016/J.CELL.2007.06.043
Jayashree, S., Karthikeyan, R., Nithyalakshmi, S., Ranjani, J., Gunasekaran, P., and Rajendhran, J. (2018). Anti-adhesion property of the potential probiotic strain Lactobacillus fermentum 8711 against methicillin-resistant Staphylococcus aureus (MRSA). Front. Microbiol. 9:411. doi: 10.3389/fmicb.2018.00411
Jeffery, C. J. (2019). Mini-review intracellular/surface moonlighting proteins that aid in the attachment of gut microbiota to the host. AIMS Microbiol. 5, 77–86. doi: 10.3934/microbiol.2019.1.77
Josse, J., Laurent, F., and Diot, A. (2017). Staphylococcal adhesion and host cell invasion: fibronectin-binding and other mechanisms. Front. Microbiol. 8:2433. doi: 10.3389/fmicb.2017.02433
Karczewski, J., Troost, F. J., Konings, I., Dekker, J., Kleerebezem, M., Brummer, R. J. M., et al. (2010). Regulation of human epithelial tight junction proteins by Lactobacillus plantarum in vivo and protective effects on the epithelial barrier. Am. J. Physiol. Gastrointest. Liver Physiol. 298, G851–G859. doi: 10.1152/AJPGI.00327.2009
Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4
Kiousi, D. E., Efstathiou, C., Tzampazlis, V., Plessas, S., Panopoulou, M., Koffa, M., et al. (2023). Genetic and phenotypic assessment of the antimicrobial activity of three potential probiotic lactobacilli against human enteropathogenic bacteria. Front. Cell. Infect. Microbiol. 13:1127256. doi: 10.3389/fcimb.2023.1127256
Kolls, J. K., McCray, P. B., and Chan, Y. R. (2008). Cytokine-mediated regulation of antimicrobial proteins. Nat. Rev. Immunol. 8, 829–835. doi: 10.1038/NRI2433
Li, X., Liu, X., Horvatovich, P., Hu, Y., and Zhang, J. (2021). Proteomics landscape of host-pathogen interaction in Acinetobacter baumannii infected mouse lung. Front. Genet. 12:563516. doi: 10.3389/fgene.2021.563516
Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/BIOINFORMATICS/BTT656
Liu, Y. S., Zhang, C., Khoo, B. L., Hao, P., and Chua, S. L. (2024). Dual-species proteomics and targeted intervention of animal-pathogen interactions. J. Adv. Res. S2090-1232, 00383–00387. doi: 10.1016/j.jare.2024.08.038
Llewellyn, A., and Foey, A. (2017). Probiotic modulation of innate cell pathogen sensing and signaling events. Nutrients 9:1156. doi: 10.3390/NU9101156
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Matsumoto, H., and Young, G. M. (2009). Translocated effectors of Yersinia. Curr. Opin. Microbiol. 12, 94–100. doi: 10.1016/J.MIB.2008.12.005
Mirdita, M., Schütze, K., Moriwaki, Y., Heo, L., Ovchinnikov, S., and Steinegger, M. (2022). ColabFold: making protein folding accessible to all. Nat. Methods 19, 679–682. doi: 10.1038/s41592-022-01488-1
Muscariello, L., De Siena, B., and Marasco, R. (2020). Lactobacillus cell surface proteins involved in interaction with mucus and extracellular matrix components. Curr. Microbiol. 77, 3831–3841. doi: 10.1007/S00284-020-02243-5
Pajarillo, E. A. B., Kim, S. H., Valeriano, V. D., Lee, J. Y., and Kang, D. K. (2017). Proteomic view of the crosstalk between Lactobacillus mucosae and intestinal epithelial cells in co-culture revealed by Q Exactive-based quantitative proteomics. Front. Microbiol. 8:2459. doi: 10.3389/fmicb.2017.02459
Pakbin, B., Brück, W. M., and Rossen, J. W. A. (2021). Virulence factors of enteric pathogenic Escherichia coli: a review. Int. J. Mol. Sci. 22:9922. doi: 10.3390/IJMS22189922
Pang, Y., Zhang, Y., Zhang, H. Y., Wang, W. H., Jin, G., Liu, J. W., et al. (2022). MUC13 promotes lung cancer development and progression by activating ERK signaling. Oncol. Lett. 23:37. doi: 10.3892/OL.2021.13155
Pavli, F. G., Argyri, A., Papadopoulou, O. S., Nychas, G.-J. E., Chorianopoulos, N. G., and Tassou, C. C. (2016). Probiotic potential of lactic acid bacteria from traditional fermented dairy and meat products: assessment by in vitro tests and molecular characterization. J. Probiotics Health 4:3. doi: 10.4172/2329-8901.1000157
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi: 10.1038/nbt.3122
Plessas, S., Kiousi, D. E., Rathosi, M., Alexopoulos, A., Kourkoutas, Y., Mantzourani, I., et al. (2020). Isolation of a Lactobacillus paracasei strain with probiotic attributes from kefir grains. Biomedicines 8, 1–15. doi: 10.3390/BIOMEDICINES8120594
Poimenidou, S. V., Skarveli, A., Saxami, G., Mitsou, E. K., Kotsou, M., and Kyriacou, A. (2023). Inhibition of Listeria monocytogenes growth, adherence and invasion in Caco-2 cells by potential probiotic lactic acid Bacteria isolated from fecal samples of healthy neonates. Microorganisms 11:363. doi: 10.3390/MICROORGANISMS11020363
Pokharel, P., Dhakal, S., and Dozois, C. M. (2023). The diversity of Escherichia coli pathotypes and vaccination strategies against this versatile bacterial pathogen. Microorganisms 11:344. doi: 10.3390/MICROORGANISMS11020344
Raineri, E. J. M., Maaß, S., Wang, M., Brushett, S., Palma Medina, L. M., Sampol Escandell, N., et al. (2022). Staphylococcus aureus populations from the gut and the blood are not distinguished by virulence traits—a critical role of host barrier integrity. Microbiome 10:239. doi: 10.1186/s40168-022-01419-4
Ramos, S., Silva, V., Dapkevicius, M., Caniça, M., Tejedor-Junco, M., Igrejas, G., et al. (2020). Escherichia coli as commensal and pathogenic bacteria among food-producing animals: health implications of extended spectrum β-lactamase (ESBL) production. Animals 10:2239. doi: 10.3390/ANI10122239
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43:e47. doi: 10.1093/NAR/GKV007
Sengupta, R., Altermann, E., Anderson, R. C., McNabb, W. C., Moughan, P. J., and Roy, N. C. (2013). The role of cell surface architecture of lactobacilli in host-microbe interactions in the gastrointestinal tract. Mediat. Inflamm. 2013:237921. doi: 10.1155/2013/237921
Sheng, Y. H., He, Y., Hasnain, S. Z., Wang, R., Tong, H., Clarke, D. T., et al. (2016). MUC13 protects colorectal cancer cells from death by activating the NF-κB pathway and is a potential therapeutic target. Oncogene 36, 700–713. doi: 10.1038/onc.2016.241
Silva, D. R., de Cássia Orlandi Sardi, J., de Souza Pitangui, N., Roque, S. M., da Silva, A. C. B., and Rosalen, P. L. (2020). Probiotics as an alternative antimicrobial therapy: current reality and future directions. J. Funct. Foods 73:104080. doi: 10.1016/J.JFF.2020.104080
Sinha, B., and Fraunholz, M. (2010). Staphylococcus aureus host cell invasion and post-invasion events. Int. J. Med. Microbiol. 300, 170–175. doi: 10.1016/J.IJMM.2009.08.019
Spacova, I., O’Neill, C., and Lebeer, S. (2020). Lacticaseibacillus rhamnosus GG inhibits infection of human keratinocytes by Staphylococcus aureus through mechanisms involving cell surface molecules and pH reduction. Benef. Microbes 11, 703–715. doi: 10.3920/BM2020.0075
Stastna, M. (2024). The role of proteomics in identification of key proteins of bacterial cells with focus on probiotic bacteria. Int. J. Mol. Sci. 25:8564. doi: 10.3390/ijms25168564
Stergiou, O. S., Tegopoulos, K., Kiousi, D. E., Tsifintaris, M., Papageorgiou, A. C., Tassou, C. C., et al. (2021). Whole-genome sequencing, phylogenetic and genomic analysis of Lactiplantibacillus pentosus L33, a potential probiotic strain isolated from fermented sausages. Front. Microbiol. 12:746659. doi: 10.3389/FMICB.2021.746659
Tegopoulos, K., Stergiou, O. S., Kiousi, D. E., Tsifintaris, M., Koletsou, E., Papageorgiou, A. C., et al. (2021). Genomic and phylogenetic analysis of Lactiplantibacillus plantarum L125, and evaluation of its anti-proliferative and cytotoxic activity in cancer cells. Biomedicines 9:1718. doi: 10.3390/BIOMEDICINES9111718
Teufel, F., Almagro Armenteros, J. J., Johansen, A. R., Gíslason, M. H., Pihl, S. I., Tsirigos, K. D., et al. (2022). SignalP 6.0 predicts all five types of signal peptides using protein language models. Nat. Biotechnol. 40, 1023–1025. doi: 10.1038/s41587-021-01156-3
Tuo, Y., Song, X., Song, Y., Liu, W., Tang, Y., Gao, Y., et al. (2018). Screening probiotics from Lactobacillus strains according to their abilities to inhibit pathogen adhesion and induction of pro-inflammatory cytokine IL-8. J. Dairy Sci. 101, 4822–4829. doi: 10.3168/JDS.2017-13654
van Sorge, N. M., Bonsor, D. A., Deng, L., Lindahl, E., Schmitt, V., Lyndin, M., et al. (2021). Bacterial protein domains with a novel Ig-like fold target human CEACAM receptors. EMBO J. 40:e106103. doi: 10.15252/EMBJ.2020106103
Wang, B., Li, S., Dedhar, S., and Cleary, P. P. (2007). Paxillin phosphorylation: bifurcation point downstream of integrin-linked kinase (ILK) in streptococcal invasion. Cell. Microbiol. 9, 1519–1528. doi: 10.1111/J.1462-5822.2007.00889.X
Wells, J. M. (2011). Immunomodulatory mechanisms of lactobacilli. Microb. Cell Fact. 10:S17. doi: 10.1186/1475-2859-10-S1-S17
Wine, E., Gareau, M. G., Johnson-Henry, K., and Sherman, P. M. (2009). Strain-specific probiotic (Lactobacillus helveticus) inhibition of Campylobacter jejuni invasion of human intestinal epithelial cells. FEMS Microbiol. Lett. 300, 146–152. doi: 10.1111/J.1574-6968.2009.01781.X
Yan, Y., Tao, H., He, J., and Huang, S. Y. (2020). The HDOCK server for integrated protein–protein docking. Nat. Protocols 15, 1829–1852. doi: 10.1038/s41596-020-0312-x
Ye, Q., Lao, L., Zhang, A., Qin, Y., Zong, M., Pan, D., et al. (2023). Multifunctional properties of the transmembrane LPxTG-motif protein derived from Limosilactobacillus reuteri SH-23. J. Dairy Sci. 106, 8207–8220. doi: 10.3168/JDS.2023-23440
Younes, J. A., Reid, G., van der Mei, H. C., and Busscher, H. J. (2016). Lactobacilli require physical contact to reduce staphylococcal TSST-1 secretion and vaginal epithelial inflammatory response. Pathog. Dis. 74:ftw029. doi: 10.1093/FEMSPD/FTW029
Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11:R14. doi: 10.1186/gb-2010-11-2-r14
Yount, N. Y., Waring, A. J., Gank, K. D., Welch, W. H., Kupferwasser, D., and Yeaman, M. R. (2006). Structural correlates of antimicrobial efficacy in IL-8 and related human kinocidins. Biochim. Biophys. Acta 1768, 598–608. doi: 10.1016/j.bbamem.2006.11.011
Yu, N. Y., Wagner, J. R., Laird, M. R., Melli, G., Rey, S., Lo, R., et al. (2010). PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics 26, 1608–1615. doi: 10.1093/BIOINFORMATICS/BTQ249
Keywords: Lactiplantibacillus plantarum , probiotics, pathogen-induced cytotoxicity, pathogen invasion, competitive exclusion, dual RNA-seq, transcriptomics, proteomics
Citation: Kiousi DE, Panopoulou M, Pappa A and Galanis A (2024) Lactobacilli-host interactions inhibit Staphylococcus aureus and Escherichia coli-induced cell death and invasion in a cellular model of infection. Front. Microbiol. 15:1501119. doi: 10.3389/fmicb.2024.1501119
Edited by:
Natasa Golic, University of Belgrade, SerbiaReviewed by:
Becky Hess, Pacific Northwest National Laboratory (DOE), United StatesAnita Tewari, Nanaji Deshmukh Veterinary Science University, India
Copyright © 2024 Kiousi, Panopoulou, Pappa and Galanis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alex Galanis, YWdhbGFuaXNAbWJnLmR1dGguZ3I=