- 1Department of Bioscience and Technology for Food, Agriculture and Environment, University of Teramo, Teramo, Italy
- 2Istituto Zooprofilattico Sperimentale dell'Abruzzo e del Molise “G. Caporale”, Teramo, Italy
Preventing L. monocytogenes infection is crucial for food safety, considering its widespread presence in the environment and its association with contaminated RTE foods. The pathogen's ability to persist under adverse conditions, for example, in food processing facilities, is linked to virulence and resistance mechanisms, including biofilm formation. In this study, the protein expression patterns of two L. monocytogenes 1/2a strains, grown under environmental stressors (mild acidic pH, thermal abuse, and high concentration of NaCl), were investigated. Protein identification and prediction were performed by nLC-ESI-MS/MS and nine different bioinformatic software programs, respectively. Gene enrichment analysis was carried out by STRING v11.05. A total of 1,215 proteins were identified, of which 335 were non-cytosolic proteins and 265 were immunogenic proteins. Proteomic analysis revealed differences in protein expression between L. monocytogenes strains in stressful conditions. The two strains exhibited unique protein expression profiles linked to stress response, virulence, and pathogenesis. Studying the proteomic profiles of such microorganisms provides information about adaptation and potential treatments, highlighting their genetic diversity and demonstrating the utility of bioinformatics and proteomics for a broader analysis of pathogens.
1. Introduction
Listeria monocytogenes is a pathogenic bacterium that is generally recognized as an important hazard in the food industry, as it causes outbreaks of listeriosis mainly in pregnant women, newborns, elderly, and individuals with compromised immune systems (Shevchuk, 2021). This microorganism is ubiquitous, and therefore the most effective control strategy for listeriosis involves the prevention of food contamination (Carpentier and Cerf, 2011; Dhama et al., 2015). The crucial niches for its transmission to animals, plants, and the food chain are mainly soil and water (Linke et al., 2014), and the outbreaks are often related to the consumption of contaminated ready-to-eat (RTE) foods, such as RTE meat products, soft cheeses, and fresh produce (Zhang et al., 2021). The pathogen is able to survive and persist under unfavorable conditions (i.e., desiccation, heat, high salt content, and refrigerated temperatures) in food processing facilities for extended periods of time due to ineffective cleaning and sanitation practices in difficult access sites, insufficient controls of movements of employees or equipment in areas with different risks of contamination, and a lack of design of preventive strategies (Rodríguez-López et al., 2018; Belias et al., 2022).
Foods are complex matrices in which bacteria can be exposed to natural stressful conditions with consequent physiological and structural changes in the cells, which allow them not only to withstand but also to grow and to express virulence-associated genes (Alvarez-Ordóñez et al., 2015; Fang et al., 2016). It is assumed that biofilm production by L. monocytogenes is one of the ways that permits increased resistance and persistence in the food chain (Lee et al., 2019), as cells resident in biofilm remain adhered to and protected by the extracellular polymeric matrix (Hobley et al., 2015). Furthermore, Li et al. (2020) demonstrated how some protein expressions, such as Lmo2193, encoded by the oppD gene, can affect positively both biofilm formation and virulence. However, the prevalence of L. monocytogenes in food facilities may not automatically be comparable to its occurrence in food products (Buchanan et al., 2017).
Therefore, rapid, accurate, sensitive, and high-throughput methods to detect and identify such microorganisms are important requirements for both the scientific community and the food industry (André et al., 2017). Traditionally, foodborne pathogen identification was carried out by culture-based approaches and the study of morphological, biochemical, and physiological microorganism characteristics. However, such techniques require some days or weeks before obtaining results, and dead cells or toxic substances cannot be revealed. Then, molecular techniques (i.e., DNA-based methods or antibody-based assays) have been developed to provide faster and more reliable data.
Currently, proteomics, genomics, metabolomics, and transcriptomics, denoted as “foodomics” when used together with bioinformatics tools, can give significant information on the biomolecules produced by food or organisms (Abril et al., 2022). These techniques have allowed scientific researchers to gain insights into the virulence mechanisms of L. monocytogenes, the genes and proteins involved in its pathogenesis, as well as the host–pathogen interactions that underlie listeriosis (Palumbo et al., 2022).
Proteomics enables the identification, characterization, and quantification of a comprehensive range of proteins upregulated or downregulated by a live organism when exposed to specific growth conditions (Chen et al., 2016), revealing protein function and expression (Böhme et al., 2013). Furthermore, proteomics methods are useful to identify structural and non-structural proteins involved in virulence as well as to investigate metabolic and physiological characteristics. Utilizing the combination of mass spectrometry analysis with different gel-based (i.e., SDS-PAGE, one- or two-dimensional gel electrophoresis) or shotgun proteomics approaches, researchers can detect bacterial infections more effectively, providing information on the proteins and the virulence factors involved in pathogenicity. In fact, understanding protein functions is crucial for the prevention and control of food pathogens, as well as for improving diagnosis (Zubair et al., 2022). Proteomic analyses of L. monocytogenes allow us to identify the proteins involved in the various aspects of its pathogenesis, including adhesion, invasion, intracellular survival, and virulence (Karthikeyan et al., 2019). Some studies have shown that L. monocytogenes-expressed proteins are implicated in the formation of cell surface structures to adhere to and invade host cells, but also in intracellular survival to change the host cell environment and promote bacterial growth (Birk et al., 2021).
In the present study, we investigated the protein expression patterns of L. monocytogenes involved in the adaptation to various stressful conditions using proteomic methods. Notably, previous studies have also explored similar approaches to gain insights into proteomic studies on the growth of such microorganisms and their adaptation to mild acid, thermal abuse, and high salt concentrations (D'Onofrio et al., 2023). For example, Mata et al. (2015) conducted a comprehensive proteomic analysis to examine the influence of strain genetics and temperature. Furthermore, Bowman et al. (2012) used similar techniques to study acid tolerance responses, while Abram et al. (2008) identified components contributing to acid and salt tolerance. However, proteomic studies on the growth and adaptation of L. monocytogenes to mild acid, low temperature, and high salt concentration at the same time are lacking. Moreover, profiling the whole protein expression to investigate systemic modifications in the microbial world is still a pioneering research field.
Therefore, the aim of the present study is the investigation of the protein expression patterns involved in the adaptation of two L. monocytogenes strains to various stressful conditions. A combination set of cardinal growth parameters (12°C as thermal abuse, maximum and minimum values of pH and NaCl), considered realistic for pork products, was carried out to observe the survival and multiplication of L. monocytogenes serotype 1/2a strains. Furthermore, this study aimed to identify the proteins involved in the virulent pathways during the exponential growth phase, comparing in vitro data with optimal growth conditions.
2. Materials and methods
2.1. Cultivation of L. monocytogenes serotype 1/2a strains
Two L. monocytogenes serotype 1/2a strains (268M “ATCC® BAA-679™” as reference strain and 3178-CB-2018 as wild type, isolated from pork sausage) were grown at different conditions (C1-opt: 37°C, pH 7.0, NaCl 0.5%; C2: 37°C, pH 5.5, NaCl 7%; C3: 12°C, pH 7.0, NaCl 0.5%; C4: 12°C, pH 5.5, NaCl 7%). Biological and technical triplicates were performed for each experimental condition using Brain Heart Infusion (BHI) broth (Oxoid Thermo Fisher Scientific, Rodano, Italy) (NaCl 0.5%, pH 7.0) and modified BHI broth (NaCl 7%, pH 5.5), according to D'Onofrio et al. (2022). Bacterial cells were grown in a 15 mL tube (OD600 = 0.9), diluted 1:100 in BHI and modified BHI broth, and then cultivated under continuous agitation. Bacterial cells were collected at the late exponential growth phase (OD600 = 0.6), after 6 h, 3 days, 24 h, and 14 days for C1-opt, C2, C3, and C4, respectively and washed with sterile ice-cold 0.01 M phosphate-buffered saline (PBS) pH 7.0. The pellets were centrifuged without adding any cryoprotective agent (Eppendorf, Hamburg, Germany) at 5,600 × g at 4°C for 10 min and stored at −80°C until use.
2.2. Protein extraction
The proteins were extracted by CelLytic B Cell Lysis Reagent (Sigma-Aldrich, Milan, Italy) and CelLytic IB Inclusion Body Solubilization Reagent (Sigma-Aldrich, Milan, Italy). The lysates were precipitated by precipitation buffer (TCA 6.1 N), incubated for 45 min at 4°C, and centrifuged (10,000 g × 15 min at 4°C), washed with 100 μl of cold acetone, and dried under vacuum (Sanchez, 2001). The samples were resuspended in a solubilization solution (10 mM Tris–HCl; pH 7.5, 200 mM NaCl; 1 mM PMSF). The PierceTM BCA Protein Assay Kit (Thermo Fisher Scientific, Rodano, Italy) was used to evaluate protein concentrations. The protein extracts were resolved by SDS-PAGE loading 12 μg of proteins per well using NuPAGETM 4–12% Bis-Tris pre-cast gels (Life Technologies Thermo Fisher Scientific, Monza, Italy) at 200 V.
2.3. Mass spectrometry analysis
L. monocytogenes purified lysates were in-gel digested, excising the lanes with a scalpel (Shevchenko et al., 1996).
The gel lanes (Supplementary material) were subjected to a reduction-alkylation process by 10 mM dithiothreitol (DTT) and 55 mM iodoacetamide (IAA), respectively, and trypsin digested at 37°C overnight. Peptides were desalted by StageTip C18 (Rappsilber et al., 2003), dried, and resuspended in 5% formic acid. A peptide mixture (5 μL) was injected on a quadrupole Orbitrap Q-Exactive HF mass spectrometer coupled with a UHPLC Easy-nLC 1200 (Thermo Fisher Scientific, Waltham, MA, USA) with a 25 cm fused-silica emitter of 75 μm inner diameter (New Objective), packed in-house with ReproSil-Pur C18-AQ 1.9 μm beads (Dr. Maisch GmbH, Ammerbuch, Germany). Peptide separation was conducted with a single run time of 33 min using a linear gradient (23 min) from 95% solvent A (2% ACN and 0.1% formic acid) to 50% solvent B (80% acetonitrile and 0.1% formic acid) and from 50 to 100% solvent B (2 min) at a constant flow rate of 0.25 μl/min.
A data-dependent acquisition (DDA) top 15 method was performed. Survey full scan MS spectra (300–1,750 Th) were acquired in the Orbitrap with 60,000 resolution, AGC target 1e6, and IT 120 ms. For HCD spectra, resolution was set to 15,000, AGC target 1e5, IT 120 ms; normalized collision energy 28%, isolation width of 3.0 m/z, and a dynamic exclusion of 5 s.
2.4. MS data analysis
Raw MS files were processed using Proteome Discoverer (version 1.4.1.14, Thermo Fisher Scientific, Waltham, MA, USA). MS/MS peak lists were searched by Mascot engine (version 2.6.0, Matrix Science, Boston, MA, USA) against the database “uniprot_listeria_monocytogenes” [setting parameters: enzyme trypsin; maximum missed cleavage 2; fixed modification carbamidomethylation (C); variable modification oxidation (M); and protein N-terminal acetylation; peptide mass tolerance 10 ppm; MS/MS tolerance 20 mmu]. Scaffold (version Scaffold_4.3.3, Proteome Software Inc., Portland, OR) was used to validate peptide and protein identification by the Scaffold Local FDR algorithm (peptide and protein threshold >95.0 and 99.0%, respectively). Furthermore, protein identifications were accepted if they contained at least two identified peptides. Protein probabilities were assigned by the Protein Prophet algorithm (Nesvizhskii et al., 2003). The proteins were grouped to satisfy the principles of parsimony, and those sharing significant peptide evidence were grouped into clusters. The proteins identified in at least two out of three biological replicates were included in the protein prediction analysis. The unique and commonly identified proteins among experimental conditions were grouped by Venn Diagram using Venny 2.0 software.
2.5. Protein prediction analysis
The proteins were analyzed by five software programs to evaluate subcellular localization (SCL): PSORTb version 3.0.2 and CELLO version 2.5 were used as bacterial protein SCL predictors (Krogh et al., 2001; Nesvizhskii et al., 2003; Rahman et al., 2008). LipoP 1.0 server was used to predict lipoproteins and to differentiate lipoprotein signal peptides (SPs) (Shevchenko et al., 1996); TMHMM Server version 2.0 was used to evaluate the prediction of transmembrane helices (Rappsilber et al., 2003; Perez-Riverol et al., 2022), and SignalP 4.1 server predicts the presence and location of SP cleavage sites in amino acid sequences (Nesvizhskii et al., 2003). The immunogenic proteins were examined to evaluate the immunogenic candidates by means of VirulentPred and VaxiJen 2.0 servers (Garp and Gupta, 2008; Petersen et al., 2009; Jespersen et al., 2017). VirulentPred is predictive software for the identification of bacterial virulent proteins. This method is founded on a two-layer cascade support vector machine (SVM) approach, employing a preset threshold of 0.0. VaxiJen operates as an autonomous server focused on foretelling the physicochemical traits of antigens and proteins without relying on sequence alignment. In terms of data evaluation through VaxiJen, the selection of proteins is influenced by their “Minimum Adhesin Probability” (MAP) score. In VaxiJen, a MAP varying between 0.4 and 0.5 indicates that a protein is a potential adhesin or possesses adhesin-like characteristics (D'Onofrio et al., 2022). The immunogenic proteins were analyzed to predict the protein immunogenic regions by NetSurfP version 1.1 (Srivastava et al., 2021) and BepiPred version 1.0 server, the former to predict the secondary structure and the relative surface exposure of the individual amino acid residues, and the latter to predict the B-cell epitopes (Lin et al., 2009).
2.6. STRING v.11.05 analysis
The in silico enrichment analysis was performed using STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) software version 11.05. The proteins with epitopes exposed to solvent (EETS), identified among each condition, were listed in the Supplementary material. The proteins were analyzed, setting a minimum required interaction score of 0.700, and visualized by the network. The proteins were clustered by the K-means clustering algorithm to highlight the different enriched pathways. Different colors were used to evaluate the functional characteristics of clustered proteins. Gene Ontology Biological Processes (GOBPs), Gene Ontology Molecular Functions (GOMFs), and text mining functions were used to describe the network functional enrichment analysis. Moreover, GOBPs were obtained from ShinyGO 0.77 for 3178-CB-2018 C2 and C3 conditions. Data visualization was performed using R version 4.3.0, utilizing the ggplot2 package.
3. Results
In the present study, the protein expression profiles of L. monocytogenes in response to different environmental conditions, considered realistic for pork products, and specifically mild acidic pH, thermal abuse, and high osmolarity level (C1-opt: 37°C, pH 7.0, NaCl 0.5%; C2: 37°C, pH 5.5, NaCl 7%; C3: 12°C, pH 7.0, NaCl 0.5%; C4: 12°C, pH 5.5, NaCl 7%), were examined.
L. monocytogenes serotype 1/2a 268M (accession number: ASM19603v1; https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_000196035.1/) was chosen because its sequence was already analyzed in a previous study (Palma et al., 2022), and it is the reference strain used in many publications on the genomic organization of this species. L. monocytogenes serotype 1/2a 3178-CB-2018 (BioSample accession SAMN37347693, BioProject: PRJNA1015181) was studied for its genetic characteristics, which are similar to a strain that caused an outbreak in Italy (Orsini et al., 2018).
The primary focus was to discern the proteins implicated in the virulence pathways of L. monocytogenes during its exponential growth phase and compare them with those expressed under optimal growth conditions (C1-opt). Through the integration of proteomic and immunoinformatic results, it was possible to identify mechanisms enabling L. monocytogenes to endure stressful environments that could potentially enhance its pathogenicity.
The total number of proteins identified by nLC-ESI-MS/MS was 1,215, considering all experimental conditions. Unique and overlapping proteins for each cultivated strain (268M and 3178-CB-2018) were grouped by Venn diagram (Figure 1). The reference (268M) and the wild-type (3178-CB-2018) strains presented 126 vs. 142 proteins in C1-opt, 158 vs. 104 in C2, 234 vs. 102 in C3, and 125 vs. 120 in C4, respectively. Out of these proteins, 87 and 72 were common in 268M and 3178-CB-2018, respectively. The SCL was performed by five software programs (CELLO version 2.5, PSORTb version 3.0.2, LipoP 1.0 Server, TMHMM Server version 2.0, and SignalP 4.1 Server). The SCL analysis led to the identification of 335 non-cytosolic proteins (NCPs). The immunogenic proteins predicted among NCP by VirulentPred and VaxiJen 2.0 servers were 265. The candidates with EETS predicted via NetSurfP v.1.1 and BepiPred v.1.0 were 259 (Figure 2). All results are available as Supplementary material.
Figure 1. Venn diagram of identified proteins for two L. monocytogenes serotype 1/2a strains: (A) 268M and (B) 3178-CB-2018. C1-opt: 37°C, pH 7.0, NaCl 0.5%; C2: 37°C, pH 5.5, NaCl 7%; C3: 12°C, pH 7.0, NaCl 0.5%; and C4: 12°C, pH 5.5, NaCl 7%.
Figure 2. Immunoinformatic workflow for immunogenic protein identification of 268M and 3178-CB-2018 L. monocytogenes strains after exposure to C1-opt: 37°C, pH 7.0, NaCl 0.5%; C2: 37°C, pH 5.5, NaCl 7%; C3: 12°C, pH 7.0, NaCl 0.5%; and C4: 12°C, pH 5.5, NaCl 7%.
GOBPs and GOMFs enriched for each experimental condition are shown in Figures 3, 4. GOBPs and GOMFs represent the biological processes in which genes are involved and the molecular functions of genes (i.e., enzyme activity, receptor binding, or transporter activity), respectively. All the proteins identified among C1-opt, C2, C3, and C4 for both strains were used to perform the enrichment analysis. The most enriched GOBPs in C1-opt, C2, C3, and C4 were protein metabolic process (26 genes), transport (46 genes), localization (81 genes), and macromolecule metabolic process (42 genes) for the 268M strain, whereas protein metabolic process (27 genes), membrane (31 genes), membrane (33 genes), and macromolecule biosynthetic process (23 genes) were the most enriched GOBPs for 3178-CB-2018. On the other side, the most enriched GOMFs for the 268M strain in C1-opt and C2 were 44 and 50 genes in small molecule binding, respectively, while in C3 and C4, they were transmembrane transporter activity (54 genes) and organic and heterocyclic compound binding (58 genes), respectively. GOMFs results for the wild-type strain showed that hydrolase activity (46 genes), organic and heterocyclic compound binding (51 genes), organic and heterocyclic compound binding (54 genes), and small molecule binding (40 genes) were enriched for C1-opt, C2, C3, and C4.
Figure 3. Data visualization obtained by R version 4.3.0 (ggplot2 package) of in silico enrichment analysis performed by ShinyGO 0.77, of GOBPs expressed by the 268M and 3178-CB-2018 L. monocytogenes strains after exposure to C1-opt: 37°C, pH 7.0, NaCl 0.5%; C2: 37°C, pH 5.5, NaCl 7%; C3: 12°C, pH 7.0, NaCl 0.5%; and C4: 12°C, pH 5.5, NaCl 7%. The blue bars indicate the number of genes (x-axes) that are codified for the proteins involved in each biological process.
Figure 4. Data visualization obtained by R version 4.3.0 (ggplot2 package) of in silico enrichment analysis performed by ShinyGO 0.77 of GOMFs expressed by the 268M and 3178-CB-2018 L. monocytogenes strains after exposure to C1-opt: 37°C, pH 7.0, NaCl 0.5%; C2: 37°C, pH 5.5, NaCl 7%; C3: 12°C, pH 7.0, NaCl 0.5%; and C4: 12°C, pH 5.5, NaCl 7%. The blue bars indicate the number of genes (x-axes) that are codified for the proteins involved in each biological process.
The in silico enrichment analysis led to visualize eight different networks based on the experimental growth conditions for 268M and 3178-CB-2018. The proteins are represented by nodes and identified by their code genes. The edges represent the protein–protein interactions and the different interactions among them. The nodes of the same colors are related to the same molecular pathway.
In 268M, 126 nodes and 338 edges for C1-opt (Figure 5A) were displayed. A total of 94 proteins were identified and annotated in 6 GOBPs, while 289 proteins were identified and included in 9 GOMFs. A total of 158 nodes and 152 edges were identified in response to high osmolarity level and low pH (C2), as shown in Figure 5B. A total of 127 and 406 proteins were identified and annotated across 6 GOBPs and 12 GOMFs, respectively. The analysis of the effects of thermal abuse (C3) identified 234 nodes and 314 edges (Figure 5C). In addition, a total of 456 proteins were annotated across 14 GOBPs, while 406 proteins were explained across 10 GOMFs. Against thermal abuse, mild acid conditions, and high osmolarity levels, 125 nodes and 124 edges were identified (Figure 5D). Furthermore, the analysis identified a total of 204 proteins annotated across 11 GOBPs, along with 623 proteins noted across 23 GOMFs.
Figure 5. (A) Network STRING v11.05 diagram showing the expressed proteins by the 268M L. monocytogenes 1/2a strain after exposure to C1-opt: 37°C, pH 7.0, NaCl 0.5%. (B) Network STRING v11.05 diagram showing the expressed proteins by the 268M L. monocytogenes 1/2a strain after exposure to C2: 37°C, pH 5.5, NaCl 7%. (C) Network STRING v11.05 diagram showing the expressed proteins by 268 M L. monocytogenes 1/2a strain after exposure to C3: 12°C, pH 7.0, NaCl 0.5%. (D) Network STRING v11.05 diagram showing the expressed proteins by 268 M L. monocytogenes 1/2a strain after exposure to C4: 12°C, pH 5.5, NaCl 7%. Nodes = protein; edges = protein–protein interactions; known interactions = pink and light blue colors; predicted interactions = green, red, and blue colors; other interactions = yellow, black, and gray colors; empty nodes = proteins of unknown 3D structure; filled nodes = some 3D structure is known or predicted.
Considering the 3178-CB-2018 strain, 142 nodes and 167 edges were identified for C1-opt (Figure 6A). The gene enrichment analysis highlighted 10 GOBPs and 13 GOMFs, including 150 and 284 proteins, respectively. To face the high osmolarity level combined with the low pH (C2), 6 GOBFs and 17 GOMFs were enriched among a network of 104 nodes and 72 edges (Figure 6B). The analysis also identified 109 and 503 proteins annotated across GOBPs and GOMFs, respectively.
Figure 6. (A) Network STRING v11.05 diagram showing the expressed proteins by the 3178-CB-2018 L. monocytogenes 1/2a strain after exposure to C1-opt: 37°C, pH 7.0, NaCl 0.5%. (B) Network STRING v11.05 diagram showing the expressed proteins by the 3178-CB-2018 L. monocytogenes 1/2a strain after exposure to C2: 37°C, pH 5.5, NaCl 7%. (C) Network STRING v11.05 diagram showing the expressed proteins by the 3178-CB-2018 L. monocytogenes 1/2a strain after exposure to C3: 12°C, pH 7.0, NaCl 0.5%. (D) Network STRING v11.05 diagram showing the expressed proteins by the 3178-CB-2018 L. monocytogenes 1/2a strain after exposure to C4: 12°C, pH 5.5, NaCl 7%. Nodes = protein; edges = protein–protein interactions; known interactions = pink and light blue colors; predicted interactions = green, red, and blue colors; other interactions = yellow, black, and gray colors; empty nodes = proteins of unknown 3D structure; filled nodes = some 3D structure is known or predicted.
Thermal abuse (C3) analysis revealed a network containing 102 nodes and 65 edges (Figure 6C). Furthermore, the analysis identified 112 and 463 proteins that were annotated in 6 GOBPs and 15 GOMFs, respectively. Considering C4, 120 nodes and 133 edges were identified (Figure 6D). Furthermore, 5 GOBFs and 11 GOMFs were significantly enriched by 81 and 347 proteins, respectively.
4. Discussion
Several enriched BPs were identified across the different growth conditions in both L. monocytogenes 1/2a strains. In C1-opt, 268M and 3178-CB-2018 exhibited enrichment in the BPs related to protein metabolic processes and cellular component organization. Such shared enrichment suggests that the strains focused on maintaining cellular homeostasis rather than survival under stress conditions, and this feature was common in both strains. In C2, 268M prominently enriched the BPs associated with transport, particularly organic substance transport, cell wall organization, and biogenesis. In contrast, 3178-CB-2018 showed significant enrichment in the pathways related to membrane, transmembrane, and transmembrane helix. This may indicate that both L. monocytogenes strains (268M and 3178-CB-2018) modified cell wall and membrane structures, probably as an adaptive strategy to maintain cellular homeostasis and resist the acidic pH of the C2 condition (Cheng et al., 2015). To face thermal abuse, both strains activated genes related to the synthesis and localization of flagella. Furthermore, to withstand lower temperatures, the strains modified the cellular structures, not only at the level of the cell wall but also inside the membrane, both below 25°C. Under the most stressful conditions, the investigated strains prioritized cell wall organization and biogenesis of BPs as a strategy to improve resistance to cold temperatures and defend cellular components. Additionally, peptidoglycan-based cell wall biogenesis, associated with penicillin-binding proteins (Pbp), plays a significant role in these adaptations (Bucur et al., 2018). Notably, both strains extended their modifications beyond the cell wall, also affecting membrane structures (Guariglia-Oropeza et al., 2018; Quereda et al., 2021). Although the L. monocytogenes strains used in this study belonged to the same serogroup (1/2a), it needs to be pointed out that they are different clonal complexes (CCs): in fact, 268M and 3178-CB-2018 are CC9 and CC7, respectively. This phylogenetic distance may have led to distinct protein expression profiles. Among all the experimental conditions, the same clusters were identified for both investigated strains.
The oligopeptide permease (Opp) complex and the ribosomal protein cluster (RPC) are essential for the survival and growth of bacterial cells. The Opp complex and the RPC are involved in the uptake of nutrients and in the fundamental process of protein synthesis, respectively (Calamita et al., 2018). The former is constituted by OppA (lmo2196), OppC (lmo2194), OppD (lmo2193), and OppF (lmo2192): they are peptide transporter subunits implicated in virulence gene expression as part of cell-to-cell signaling regulation. In fact, OppD seems to be involved in biofilm formation and virulence pathways. Moreover, the Opp complex includes Lmo2360 and Lmo2361, which are transmembrane permeases as in the 268M strain (Disson et al., 2021). The RPC is related to translation and is conserved across bacterial species (rplU, rpmB, rpsT, rpmI, rpsU, rpsZ, rpmE2, rplB, rpsI, rpsL, and rpmA) (Hain et al., 2008). These ribosomal proteins play a critical role in ensuring accurate and efficient translation of mRNA into proteins. Mutations in RPC proteins have been shown to affect ribosome biogenesis and function, leading to a variety of cellular defects and diseases (Calamita et al., 2018). The locus of enterocyte effacement proliferation A (LepA) and signal recognition particle (SRP) protein Ffh (ffh) are involved in signaling pathways to induce molecular switches during stressor exposure. LepA is required to face stress conditions such as nutrient deprivation, heat shock, and oxidative stress. Ffh binds to the signal peptide of the secreted protein, ensuring proper insertion and folding. Ffh has also been implicated in the regulation of cell division (Pleitner et al., 2014). The energy coupling factors A2 (EcfA2) and ATP synthase F (AtpF) are part of the cluster in the presence of a low pH and high osmolarity (C2); these proteins are connected to RPC. EcfA2 clusters with RPC also in C4; it is a subunit transporter of ThiT, which is an integral membrane complex required for thiamine recognition and transport. L. monocytogenes needs thiamine as a cofactor because its genome lacks a partial list of genes involved in thiamine biosynthesis (Madeo et al., 2012). Moreover, ThiT is necessary to enhance acid tolerance and intracellular growth during pathogenesis (Schauer et al., 2009). Overall, the presence of EcfA2 and AtpF in this cluster suggests that they are involved in the energy metabolism and protein synthesis required for the survival and adaptation of L. monocytogenes under different environmental conditions. AtpF is a subunit of plasma membrane proton-transporting ATP synthase complex coupling factor (F) 0, which is responsible for generating ATP by using the proton motive force across the bacterial membrane and contributes to facing acidity stress in L. monocytogenes by maintaining membrane potential and pH homeostasis (Bucur et al., 2018). In C3, AtpF was connected to AtpB, which seems to be overexpressed in response to various stress conditions, such as low pH combined with anaerobiosis, and in this study was expressed to face thermal abuse (Chakravarty et al., 2021).
In C2, the connection between LepA and YidC2 suggests a potential role for LepA in protein folding and insertion during nutrient transport (Burg-Golani et al., 2013). The fact that YidC2 was identified in response to thermal abuse and was clustered in the translocase complex further highlights its importance in the adaptation of L. monocytogenes to different environmental conditions. ParA and FtsK, which were present in all investigated conditions, are implicated in the optimal replication of L. monocytogenes in macrophages, as these proteins are important for cell division and the control of protein complex localization (Fischer et al., 2022). In the same cluster, genes involved in DNA metabolism and transport, such as lmo1606, DNA repair ruvB, and transcriptional regulator birA, are included. The presence of the first two genes listed above in the same cluster as lepA and yidC2 in C2 suggests that these proteins could also play a role in DNA maintenance and repair during nutrient transport and adaptation to different environmental conditions. Then, the transcriptional regulator birA in this cluster may also be involved in the regulation of gene expression related to nutrient transport and stress response (Burg-Golani et al., 2013; Fischer et al., 2022).
In all experimental conditions, one cluster related to sugar uptake was present. Lmo0097, Lmo0098, Lmo0781, and Lmo0782 are kinase proteins involved in the carbohydrate phosphotransferase system (PTS) and are related to fructose and mannose-specific PTS. Moreover, they are positively regulated by the stress response factor σB, which also contributes to host entry by L. monocytogenes (Mattila et al., 2020). The fructose permease IIC component (FruA) protein is connected to this cluster and is positively regulated by the agr gene. Such genes have a main role in the quorum sensing system, which is associated with cell survival, adhesion, biofilm formation, invasion, and the infectivity of host cells (Pinheiro et al., 2018).
Under the optimal growth condition (C1-opt), surface proteins that contribute to L. monocytogenes host cell adhesion clustered together. Lmo1068 is involved in attachment and biofilm formation (Tiong and Muriana, 2016), Lmo0723 in chemotaxis, and Lmo2713 is upregulated by the σB factor during host cell invasion (Neuhaus et al., 2013; Piercey et al., 2016). These proteins are also visualized in cluster 3 (C3), while Iap protein was identified in all studied conditions along with autolysin protein (Lmo1076), muramidase protein MurA (lmo2691) (Tiong and Muriana, 2016), and Lmo0129, which is involved in the host gastrointestinal tract colonization (Desvaux et al., 2010).
Chaperones and proteases contribute to virulence and stress tolerance. Caseinolytic protein-P (ClpP) is a subunit of the ATP-dependent Clp protease complex, which is involved in the adaptive response of L. monocytogenes during the infection process and upregulates listeriolysin O (LLO) production at 37°C. Moreover, a similar study on a different Gram-positive pathogen, Staphylococcus aureus, showed that it is required for motility and growth at high temperatures (Savijoki et al., 2020). Lmo0292 protein, also known as “high temperature requirement A” (HtrA), acts as a chaperone and serine protease (Abfalter et al., 2019). ClpP and Lmo0292 were present in all conditions. In such clusters, proteinases of the Fts family (FtsY, FtsH, FtsE, and FtsK) and other proteinases (Lmo1851 and Lmo1077) are also included, and they are involved in cell wall metabolism and in wall teichoic acid (WTA) biosynthesis and modification, respectively. Those proteins seem to be upregulated in the presence of stress factors, such as carnocyclin A (CCLA), which is an antimicrobial peptide used against L. monocytogenes in RTE meat products (Camejo et al., 2009; Liu et al., 2014).
The penicillin-binding protein (Pbp) family is involved in cell wall metabolism, antimicrobial resistance, and infections. Lmo0441, PbpA (lmo1855), Pbp5 (lmo2754), and Pbp3 (lmo1438) were present in all experimental conditions, except for PbpB (lmo2039) that was expressed only in C3 and C4 of 268M strain, indicating a different cell wall structure when exposed to various stressors. As regards the 3178-CB-2018 strain, the above-cited Pbp were all expressed at the optimal condition (C1-opt), while Lmo2754, Lmo2229, and Lmo1855 were identified in C2, C3, and C4, respectively (Shin et al., 2010; Lima et al., 2011; Fischer et al., 2020). Guinae et al. (2006) demonstrated that strains lacking Pbp3 and Pbp5 exhibit attenuated virulence. The upregulation of some of these proteins in response to stressors, such as antimicrobial peptides, suggests that L. monocytogenes has evolved mechanisms to resist host defense systems, and this may contribute to its ability to cause infection and illnesses. The importance of Pbp3 and Pbp5 in virulence further highlights the role of cell wall metabolism in the pathogenicity of the microorganism.
The CBS domain is found in both prokaryotes and eukaryotes, with an important role in protein–protein interactions for processes, such as cell adhesion, cell activation, and molecular recognition (Oliver et al., 2010). The CBS domain (Lmo2232) was present in C1-opt both for 268M clustered to Lmo1422 and for 3178-CB-2018 to Lmo1422 and OpuCC. Lmo1422 and OpuCC are subunits of an ABC transporter complex: the former is involved in glycine betaine transport, protection by oxidative damage, and the SOS response of L. monocytogenes (Zhang et al., 2022), and the latter is concerned with the carnitine accumulation pathway (Van Der Veen et al., 2010). Moreover, both proteins are part of the σB-dependent osmotic stress resistance systems (Liu et al., 2019). These differences suggest that 3178-CB-2018 in C1-opt has a distinct surface protein profile that may contribute to its ability to adhere to both surfaces and host cells.
Proteins associated with the oxidative stress response were expressed in different conditions, suggesting L. monocytogenes adaptation to oxidative environments. In C1-opt for both 268M and 3178-CB-2018, superoxide dismutase (SOD) was connected to Pbp3 in the network. This protein is classified as an ROS scavenger and, together with both thioredoxin reductase (TrxB) and glutathione reductase (Lmo1433); it aims to address oxidative stress in bacterial cells (Savijoki et al., 2020). Hain et al. (2008) determined that Lmo1433 is σB-dependent in L. monocytogenes. PdhB is a stress response protein involved in the repair of DNA during oxidative stress (D'Onofrio et al., 2022). The same cluster was found in C4 for both strains. In C2 for 268M, more proteins related to oxidative stress were identified, as the MntA protein is linked to SOD by MntH; these proteins are expressed in this experimental condition to cope with low pH and high osmolarity (Zeng et al., 2019). In fact, manganese (Mn) transporters (MntH and MntB) were found to be highly expressed after exposure to mild acid stress (pH 5) (Wu et al., 2023).
ZurA, Lmo0153, and Lmo0154, involved as complexes in zinc uptake systems (ZurAM and ZinABC), were both expressed in C2 and C3 of 268M.
Only in C3 of 268M, there was a cluster formed by four thermosensors (Lmo2417, Lmo2418, Lmo2419, and MetN1), upregulated by the prfA gene. These proteins appear to be upregulated during the exponential growth phase and when the pathogen invades the host intestinal lumen. Among the ABC transporter systems regulated to cope with low temperatures, the Lmo0607–Lmo0608 complex, identified in all experimental conditions of 268M, acts as a multidrug transport system (Krawczyk-Balska et al., 2021). Arpj is an ABC transporter subunit for arginine, expressed only in C3 during intracellular replication and acid shock, while Lmo0847–Lmo0848 complex subunits act as glutamine ABC transporters and were found in all conditions (Neuhaus et al., 2013). Lmo0193–Lmo0195 complex acts as an ABC transporter and is upregulated during osmotic stress (C4) (Bergholz et al., 2012).
Lmo1388 (tcsA)–Lmo1389 are two-component system (TCS) proteins, positively regulated by σB regulon (Oliver et al., 2010) and present in all conditions for both strains.
Lmo0530 is part of a five-gene operon (lmo0527–lmo0531), which seems to be σBand σA dependent. It is involved in the enhanced synthesis of exopolysaccharides in response to osmolarity changes, facilitating bacterial survival under this stress condition (Ribeiro et al., 2014). In fact, Lmo0530 was expressed not only at 12°C but also in C2 and C3 of the 268M strain. Chen et al. (2014) and Köseoglu et al. (2015) demonstrated the involvement of cyclic diguanylate (c-di-GMP) signaling in L. monocytogenes and its impact on various bacterial behaviors, including the production of an exopolysaccharide, tolerance to stressors, and virulence.
To cope with low temperatures, L. monocytogenes synthesizes type II branched-chain fatty acids. Lmo0970 (FabI) is the regulatory protein reductase of type II in the fatty acid synthesis of L. monocytogenes and is essential together with Lmo2201 (FabF-B-ketoacyl-acyl-carrier-protein-synthase complex) and Lmo0814 (FabK) (Yao et al., 2016). In fact, this cluster was observed in C3 to face thermal abuse and also in C2 and C4, to overcome high salt content and low pH.
TagB, lmo1074, and lmo1091 are specific genes of the L. monocytogenes serovar 1/2a strain: the first two encode for cell wall-associated proteins (tail-anchored translocation permeases), while the second is for a glycosyl transferase that is involved in cell wall biosynthesis. Being strain-specific, these proteins could have a definite impact on host–pathogen interaction (Chatterjee et al., 2006) and were present in C3 and C4 for 268M and only in C4 for 3178-CB-2018.
Genes involved in flagellar biosynthesis, such as flhA, were only observed in C3, in response to thermal abuse. In particular, induced by lmo0681, flaA codes for flagellin, which is the principal component of the bacterial flagellum and is usually expressed below 25°C and under specific stress conditions, such as the presence of detergents (Way et al., 2004; Hain et al., 2012). FliM is a complex subunit of the flagellar motor switch, essential for the transmembrane transport of protons and rotation of the flagella, thus enabling the bacterial cell to become motile with the help of FliF proteins, which are involved in flagellar assembly (Bigot et al., 2005). Flagellar structural genes are also part of the cluster, such as flgK (hook-associated protein), flgE (flagellar hook protein), and flgG (basal body rod proteins).
Lmo0553 and Lmo1466 proteins are part of the cyclic di-adenosine monophosphate (c-di-AMP) complex, which is a largely conserved second messenger essential for L. monocytogenes growth, stress responses, antibiotic resistance, cellular morphology, and virulence (Witte et al., 2013). C-di-AMP acts on the intracellular signaling pathway and is an important component of innate immune detection during infection (Sureka et al., 2014). Witte et al. (2013) highlighted the importance of c-di-AMP as a critical signaling molecule for L. monocytogenes, essential for bacterial replication, cell wall stability, and pathogenicity.
In C1-opt, 3178-CB-2018 expressed cell surface proteins that contribute to both surface and host cell adhesion. Specifically, in contrast to C1-opt for 268M, 3178-CB-2018 expressed InlA and FlaA, which are linked to the proteins Iap and Lmo0723, respectively.
On the other hand, there were several proteins present in 268M under the optimal condition that were absent in 3178-CB-2018. These include EzrA involved in cell division, Lmo1593 in iron–sulfur cofactor synthesis, RuvB and BirA in DNA metabolism and transcription, respectively, and Lmo0653, which is a hypothetical protein that is likely elaborated in lipoteichoic acid metabolism and could contribute to nisin resistance in L. monocytogenes, according to Pang et al. (2022).
To face high osmolarity and low pH (C2), the 3178-CB-2018 L. monocytogenes strain appears to produce fewer proteins than those expressed by the strain 268M grown under the same conditions. Despite this, InlA was only found in the wild-type strain and not in the reference strain, and therefore the presence of InlA in 3178-CB-2018 suggests that it may have a greater ability to infect host cells than 268M under high osmolarity and low pH conditions.
The same response was observed in the presence of thermal abuse (C3). However, several proteins clustered only in 3178-CB-2018, including LepA involved in signaling pathways, the fts proteinases family, Lmo1077 involved in wall teichoic acid (WTA) biosynthesis and modification, ClpP and ParA related to cell division, PdhB, Lmo1606 implicated in DNA metabolism, and Lmo0653. Moreover, thermal abuse induced the expression of several proteins that appear to be involved in L. monocytogenes pathogenesis, such as Lmo1526 (Knudsen et al., 2012) and Lmo2215; the latter is expressed by the lmo2215 gene, which is one of the genes most associated with survival and stress response of L. monocytogenes, especially in meat and meat products (Tanui et al., 2022).
Lmo1529 is a YajC translocase preprotein subunit associated with protein secretion and codified by genes important for L. monocytogenes virulence (Muchaamba et al., 2022).
The lmo2637 gene encodes for an EET-linked lipoprotein, PplA, while the lmo2638 gene encodes for an EET-linked NADH dehydrogenase involved in ethanolamine (EA) metabolism (Zeng et al., 2021). EA is a source of carbon and nitrogen derived from the eukaryotic cell membrane; it can be used by a specific microbial redox metabolic subsystem called the bacterial microcompartment (BMC), inducing the activation of EA utilization (eut) genes, which are also linked to L. monocytogenes pathogenesis and growth (Zeng et al., 2021).
5. Conclusion
Investigating the proteomic profiles of L. monocytogenes under different stress conditions is of great importance for several reasons. In fact, understanding the molecular mechanisms that allow the pathogen to adapt and survive under stressful conditions is crucial for the development of effective strategies to control and prevent its spread. This study provides valuable insights into the expression patterns of specific genes and proteins under different stress conditions, which can be used to identify potential targets to discriminate between hypovirulent and hypervirulent strains of L. monocytogenes.
Furthermore, the study highlights the differences between the proteome profiles of two strains of L. monocytogenes grown under exposure to different stressors. This information can be used to better understand the genetic and phenotypic diversity of such a pathogen and to identify factors that contribute to its virulence and pathogenicity.
Finally, our results demonstrate the power of bioinformatic and proteomic approaches to analyze large datasets and identify important patterns and associations that would be difficult to detect using traditional laboratory methods. This methodology can be applied to other pathogens and stress conditions to gain a better understanding of their biology and pathogenesis.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
FD'O: Methodology, Writing—original draft, Data curation, Formal analysis, Writing—review and editing. MS: Conceptualization, Supervision, Writing—original draft, Writing—review and editing. IK: Data curation, Formal analysis, and Methodology. MTi: Investigation and Visualization. LI: Investigation, Resources, Supervision, Writing—review and editing. FP: Funding acquisition, Investigation, and Visualization. MTo: Supervision, Writing—review and editing. AP: Supervision, Writing—review and editing. ND'A: Investigation, Visualization and Funding. ML: Conceptualization, Supervision, Writing—original draft, Writing—review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The Italian Ministry of Health funded research project IZS AM 03/17 RC (study of the phenotypic and genotypic characteristics of Listeria monocytogenes subpopulations resistant to environmental stress and to disinfectants—MSRCTE0719) covered part of the analyses that were carried out for this study.
Acknowledgments
The authors are thankful to the European Union—Next Generation EU. Project Code: ECS00000041. Project Title: Innovation, digitalization, and sustainability for the diffused economy in Central Italy—VITALITY.
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.2023.1271787/full#supplementary-material
References
Abfalter, C. M., Bernegger, S., Jarzab, M., Posselt, G., Ponnuraj, K., and Wessler, S. (2019). The proteolytic activity of Listeria monocytogenes HtrA. BMC Microbiol. 19, 255. doi: 10.1186/s12866-019-1633-1
Abram, F., Starr, E., Karatzas, K. A. G., Matlawska-Wasowska, K., Boyd, A., Wiedmann, M., et al. (2008). Identification of components of the Sigma B regulon in Listeria monocytogenes that contribute to acid and salt tolerance. App. Environ. Microbiol. 74, 6848–6858. doi: 10.1128/AEM.00442-08
Abril, A. G., Villa, T. G., Calo-Mata, P., Barros-Velázquez, J., and Carrera, M. (2022). “Chapter 11 – Application of proteomics to the identification of foodborne pathogens,” in Food Proteomics: Technological Advances, Current Applications and Future Perspectives, eds M. L. Pedrouso, D. F. Ruiz, and J. M. Lorenzo (Academic Press), 337–362.
Alvarez-Ordóñez A. Broussolle V. Colin P. Nguyen-The C. Prieto M. (2015). The adaptative response of bacterial food-borne pathogens in the environment, host and food: Implications for food safety. Int. J. Food Microbiol. 213, 99–109. doi: 10.1016/j.ijfoodmicro.2015.06.004
André, S., Vallaeys, T., and Planchon, S. (2017). Spore-forming bacteria responsible for food spoilage. Res. Microbiol. 168, 379–387. doi: 10.1016/j.resmic.2016.10.003
Belias, A., Sullivan, G., Wiedmann, M., and Ivanek, R. (2022). Factors that contribute to persistent Listeria in food processing facilities and relevant interventions: a rapid review. Food Control 133, 108579. doi: 10.1016/j.foodcont.2021.108579
Bergholz, T. M., Bowen, B., Wiedmann, M., and Boor, K. J. (2012). Listeria monocytogenes shows temperature-dependent and -independent responses to salt stress, including responses that induce cross-protection against other stresses. Appl. Environ. Microbiol. 78, 2602–2612. doi: 10.1128/AEM.07658-11
Bigot, A., Pagniez, H., Botton, E., Fréhel, C., Dubail, I., Jacquet, C., et al. (2005). Role of FliF and FliI of Listeria monocytogenes in flagellar assembly and pathogenicity. ASM J. Infect. Immun. 73, 5530–5539. doi: 10.1128/IAI.73.9.5530-5539.2005
Birk, M. S., Ahmed-Begrich, R., Tran, S., Elsholz, A. K. W., Frese, C. K., and Charpentier, E. (2021). Time-resolved proteome analisis of Listeria monocytogenes during infection reveals the role of the AAA+ chaperone ClpC for host cell adaptation. mSystems 6, e0021521. doi: 10.1128/mSystems.00215-21
Böhme, K., Fernández-No, I. C., Pazos, M., Gallardo, J. M., and Barros-Velázquez, J. (2013). Identification and classification of seafood-borne pathogenic and spoilage bacteria: 16S rRNA sequencing versus MALDI-TOF MS fingerprinting. Electrophoresis 34, 877–887. doi: 10.1002/elps.201200532
Bowman, J. P., Hages, E., Nilsson, R. E., Kocharunchitt, C., and Ross, T. (2012). Investigation of the Listeria monocytogenes Scott A acid tolerance response and associated physiological and phenotypic features via whole proteome analysis. J. Proteome Res. 11, 2409–2426. doi: 10.1021/pr201137c
Buchanan, R. L., Gorris, L. G. M., Hayman, M. M., Jackson, T. C., and Whiting, R. C. (2017). A review of Listeria monocytogenes: an update on outbreaks, virulence, dose-response, ecology, and risk assessments. Food Control 75, 1–13. doi: 10.1016/j.foodcont.2016.12.016
Bucur, F. I., Grigore-Gurgu, L., Crauwels, P., Riedel, C. U., and Nicolau, A. I. (2018). Resistance of Listeria monocytogenes to stress conditions encountered in food and food processing environments. Front. Microbiol. 9, 2700. doi: 10.3389/fmicb.2018.02700
Burg-Golani, T., Pozniak, Y., Rabinovich, L., Sigal, N., Nir Paz, R., and Herskovits, A. A. (2013). Membrane chaperone SecDF plays a role in the secretion of Listeria monocytogenes major virulence factors. J. Bacteriol. 195, 5262–5272. doi: 10.1128/JB.00697-13
Calamita, P., Gatti, G., Miluzio, A., Scagliola, A., and Biffo, S. (2018). Translating the game: Ribosomes as active players. Front. Genet. 9, 533. doi: 10.3389/fgene.2018.00533
Camejo, A., Buchrieser, C., Couvé, E., Carvalho, F., Reis, O., Ferreira, P., et al. (2009). In vivo transcriptional profiling of Listeria monocytogenes and mutagenesis identify new virulence factors involved in infection. PLoS Pathog. 5, e1000449. doi: 10.1371/journal.ppat.1000449
Carpentier, B., and Cerf, O. (2011). Review-Persistence of Listeria monocytogenes in food industry equipment and premises. Int. J. Food Microbiol. 145, 1–8. doi: 10.1016/j.ijfoodmicro.2011.01.005
Chakravarty, D., Sahukhal, G., Arick, M., Davis, M. L., and Donaldson, J. R. (2021). Transcriptomic analysis of Listeria monocytogenes in response to bile under aerobic and anaerobic conditions. Front. Microbiol. 12, 754748. doi: 10.3389/fmicb.2021.754748
Chatterjee, S. S., Hossain, H., Otten, S., Kuenne, C., Kuchmina, K., Machata, S., et al. (2006). Intracellular gene expression profile of Listeria monocytogenes. Infect. Immun. 74, 1323–1338. doi: 10.1128/IAI.74.2.1323-1338.2006
Chen, L.-H., Köseoglu, V. K., Güvener, Z. T., Myers-Morales, T., Reed, J. M., D'Orazio, S. E. F., et al. (2014). Cyclic di-GMP_dependent signaling pathways in the pathogenic Firmicute Listeria monocytogenes. PLoS Pathog. 10, e1004301. doi: 10.1371/journal.ppat.1004301
Chen, T. Y., Kuo, S. H., Chen, S. T., and Hwang, D. F. (2016). Differential proteomics to explore the inhibitor effects of acidic slightly acidic electrolysed water and sodium hypochlorite solution on Vibrio parahaemolyticus. Food Chem. 194, 529–537. doi: 10.1016/j.foodchem.2015.08.019
Cheng, C., Yang, Y., Dong, Z., Wang, X., Fang, C., Yang, M., et al. (2015). Listeria monocytogenes varies among strains to maintain intracellular pH homeostasis under stresses by different acids as analyzed by a high-throughput microplate-based fluorometry. Front. Microbiol. 6, 15. doi: 10.3389/fmicb.2015.00015
Desvaux, M., Dumas, E., Chafsey, I., Chambon, C., and Hébraud, M. (2010). Comprehensive appraisal of the extracellular proteins from a monoderm bacterium: Theoretical and empirical exoproteomes of Listeria monocytogenes EGD-e by secretomics. J. Proteome Res. 9, 5076–5092. doi: 10.1021/pr1003642
Dhama, K., Karthik, K., Tiwari, R., Shabbir, M. Z., Barbuddhe, S., Malik, S. V. S., et al. (2015). Listeriosis in animals, its public health significance (food-borne zoonosis) and advances in diagnosis and control: a comprehensive review. Vet. Quart. 35, 211–235. doi: 10.1080/01652176.2015.1063023
Disson, O., Moura, A., and Lecuit, M. (2021). Making sense of the biodiversity and virulence of Listeria monocytogenes. Trends Microbiol. 29, 811–822. doi: 10.1016/j.tim.2021.01.008
D'Onofrio, F., Schirone, M., Paparella, A., Krasteva, I., Tittarelli, M., Pomilio, F., et al. (2023). Stress adaptation responses of a Listeria monocytogenes 1/2a strain via proteome profiling. Foods 12, 2166. doi: 10.3390/foods12112166
D'Onofrio, F., Visciano, P., Krasteva, I., Torresi, M., Tittarelli, M., Pomilio, F., et al. (2022). Immunoproteome profiling of Listeria monocytogenes under mild acid and salt stress conditions. Proteomics 22:2200082. doi: 10.1002/pmic.202200082
Fang, F. C., Frawley, E. R., Tapscott, T., and Vázquez-Torres, A. (2016). Bacterial stress responses during host infection. Cell Host Microbe 20, 133–143. doi: 10.1016/j.chom.2016.07.009
Fischer, M. A., Engelgeh, T., Rothe, P., Fuchs, S., Thürmer, A., and Halbedel, S. (2022). Listeria monocytogenes genes supporting growth under standard laboratory cultivation conditions and during macrophage infection. Genome Res. 32, 1711–1726. doi: 10.1101/gr.276747.122
Fischer, M. A., Wamp, S., Fruth, A., Allerberger, F., Flieger, A., and Halbedel, S. (2020). Population structure-guided profiling of antibiotic resistance patterns in clinical Listeria monocytogenes isolates from Germany identifies pbpB3 alleles associated with low levels of cephalosporin resistance. Emerg. Microbes Infec. 9, 1804–1813. doi: 10.1080/22221751.2020.1799722
Garp, A., and Gupta, D. (2008). VirulentPred: a SVM based prediction method for virulent proteins in bacterial pathogens. BMC Bioinformat. 9, 62. doi: 10.1186/1471-2105-9-62
Guariglia-Oropeza, V., Orsi, R. H., Guldimann, C., Wiedmann, M., and Boor, K. J. (2018). The Listeria monocytogenes bile stimulon under acidic conditions is characterized by strain-specific patterns and the upregulation of motility, cell wall modification functions, and the PrfA regulon. Front. Microbiol. 9, 120. doi: 10.3389/fmicb.2018.00120
Guinae, C. M., Cotter, P. D., Ross, R. P., and Hill, C. (2006). Contribution of penicillin-binding protein homologs to antibiotic resistance, cell morphology, and virulence of Listeria monocytogenes EGDe. Antimicrob. Agents Chemother. 50, 2824–2828. doi: 10.1128/AAC.00167-06
Hain, T., Ghai, R., Billion, A., Kuenne, C. T., Steinweg, C., Izar, B., et al. (2012). Comparative genomics and transcriptomics of lineages I, II, and III strains of Listeria monocytogenes. BMC Genom. 13, 144. doi: 10.1186/1471-2164-13-144
Hain, T., Hossain, H., Chatterjee, S. S., Machata, S., Volk, U., Wagner, S., et al. (2008). Temporal transcriptomic analysis of the Listeria monocytogenes EGD-e σB regulon. BMC Microbiol. 8, 20. doi: 10.1186/1471-2180-8-20
Hobley, L., Harkins, C., MacPhee, C. E., and Stanley-Wall, N. R. (2015). Giving structure to the biofilm matrix: an overview of individual strategies and emerging common themes. FEMS Microbiol. Rev. 39, 649–669. doi: 10.1093/femsre/fuv015
Jespersen, M. C., Peters, B., Nielsen, M., and Marcantili, P. (2017). BepiPred-2.0: improving sequence-based B-cell epitope prediction using conformational epitopes. Nucleic Acids Res. 45, W24–W29. doi: 10.1093/nar/gkx346
Karthikeyan, R., Gayathri, P., Gunasekaran, P., Jagannadham, M. V., and Rajendhran, J. (2019). Comprehensive proteomic analysis and pathogenic role of membrane vesicles of Listeria monocytogenes serotype 4b reveals proteins associated with virulence and their possible interaction with host. Int. J. Med. Microbiol. 309, 199–212. doi: 10.1016/j.ijmm.2019.03.008
Knudsen, G. M., Holch, A., and Gram, L. (2012). Subinhibitory concentrations of antibiotics affect stress and virulence gene expression in Listeria monocytogenes and cause enhanced stress sensitivity but do not affect Caco-2 cell invasion. J. Appl. Microbiol. 113, 1273–1286. doi: 10.1111/j.1365-2672.2012.05435.x
Köseoglu, V. K., Heiss, C., Azadi, P., Topchiy, E., Güvener, Z. T., Lehmann, T. E., et al. (2015). Listeria monocytogenes exopolysaccharide: origin, structure, biosynthetic machinery and c-di-GMP-dependent regulation. Mol. Microbiol. 96, 728–743. doi: 10.1111/mmi.12966
Krawczyk-Balska, A., Ładziak, M., Burmistrz, M., Scibek, K., and Kallipolitis, B. H. (2021). RNA-mMediated cControl in Listeria monocytogenes: insights into regulatory mechanisms and roles in metabolism and virulence. Front. Microbiol. 12, 622829. doi: 10.3389/fmicb.2021.622829
Krogh, A., Larsson, B., von Heijne, G., and Sonnhammer, E. L. L. (2001). Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J. Mol. Biol. 305, 567–580. doi: 10.1006/jmbi.2000.4315
Lee, B.-H., Cole, S., Badel-Berchoux, S., Guillier, L., Felix, B., Krezdorn, N., et al. (2019). Biofilm formation of Listeria monocytogenes strains under food processing environments and pan-genome-wide association study. Front. Microbiol. 10, 2698. doi: 10.3389/fmicb.2019.02698
Li, H., Qiao, Y., Du, D., Wang, J., and Ma, X. (2020). Deletion of the oligopeptide transporter Lmo2193 decreases the virulence of Listeria monocytogenes. J. Vet. Sci. 21, e88. doi: 10.4142/jvs.2020.21.e88
Lima, A., Durán, R., Schujman, G. E., Marchissio, M. J., Portela, M. M., Obal, G., et al. (2011). Serine/threonine protein kinase PrkA of the human pathogen Listeria monocytogenes: Biochemical characterization and identification of interacting partners through proteomic approaches. J. Proteom. 74, 1720–1734. doi: 10.1016/j.jprot.2011.03.005
Lin, H. N., Chen, C. T., Sung, T. Y., Ho, S. Y., and Hsu, W. L. (2009). Protein subcellular localization prediction of eukaryotes using a knowledge-based approach. BMC Bioinformat. 10, S8. doi: 10.1186/1471-2105-10-S15-S8
Linke, K., Rückerl, I., Brugger, K., Karpiskova, R., Walland, J., Muri-Klinger, S., et al. (2014). Reservoirs of Listeria species in three environmental ecosystems. Appl. Environ. Microbiol. 80, 5583–5592. doi: 10.1128/AEM.01018-14
Liu, X., Basu, U., Miller, P., and Mcmullen, L. M. (2014). Stress response and adaptation of Listeria monocytogenes 08-5923 exposed to a sublethal dose of carnocyclin A. App. Environ. Microbiol. 80, 3835–3841. doi: 10.1128/AEM.00350-14
Liu, Y., Orsi, R. H., Gaballa, A., Wiedmann, M., Boor, K. J., and Guariglia-Oropeza, V. (2019). Systematic review of the Listeria monocytogenes σ B regulon supports a role in stress response, virulence and metabolism. Fut. Microbiol. 14, 801–828. doi: 10.2217/fmb-2019-0072
Madeo, M., O'Riordan, N., Fuchs, T. M., Utratna, M., Karatzas, K. A. G., and O'Byrne, C. P. (2012). Thiamine plays a critical role in the acid tolerance of Listeria monocytogenes. FEMS Microbiol. Lett. 326, 137–143. doi: 10.1111/j.1574-6968.2011.02442.x
Mata, M. M., da Silva, W. P., Wilson, R., Lowe, E., and Bowman, J. (2015). Attached and planktonic Listeria monocytogenes global proteomic responses and associated influence of strain genetics and temperature. J. Proteome Res. 14, 1161–1173. doi: 10.1021/pr501114e
Mattila, M., Somervuo, P., Korkeala, H., Stephan, R., and Tasara, T. (2020). Transcriptomic and phenotypic analyses of the Sigma B-dependent characteristics and the synergism between Sigma B and Sigma L in Listeria monocytogenes EGD-e. Microorganisms 8, 1644. doi: 10.3390/microorganisms8111644
Muchaamba, F., Eshwar, A. K., Stevens, M. J. A., Stephan, R., and Tasara, T. (2022). Different shades of Listeria monocytogenes: Strain, serotype, and lineage-based variability in virulence and stress tolerance profiles. Front. Microbiol. 12, 792162. doi: 10.3389/fmicb.2021.792162
Nesvizhskii, A. I., Keller, A., Kolker, E., and Aebersold, R. (2003). A statistical model for identifying proteins by tandem mass spectrometry. Anal. Chem. 75, 4646–4658. doi: 10.1021/ac0341261
Neuhaus, K., Satorhelyi, P., Schauer, K., Scherer, S., and Fuchs, T. M. (2013). Acid shock of Listeria monocytogenes at low environmental temperatures induces prfA, epithelial cell invasion, and lethality towards Caenorhabditis elegans. BMC Genom. 14, 285. doi: 10.1186/1471-2164-14-285
Oliver, H. F., Orsi, R. H., Wiedmann, M., and Boor, K. J. (2010). Listeria monocytogenes σ B has a mall core regulon and a conserved role in virulence but makes differential contributions to stress tolerance across a diverse collection of strains. Appl. Environ. Microbiol. 76, 4216–4232. doi: 10.1128/AEM.00031-10
Orsini, M., Cornacchia, A., Patavino, C., Torresi, M., Centorame, P., Acciari, V. A., et al. (2018). Whole-genome sequences of two Listeria monocytogenes serovar 1/2a strains responsible for a severe listeriosis outbreak in Central Italy. Genome Announc. 6, e00236–e00218. doi: 10.1128/genomeA.00236-18
Palma, F., Mangone, I., Janowicz, A., Moura, A., Chiaverini, A., Torresi, M., et al. (2022). In vitro and in silico parameters for precise cgMLST typing of Listeria monocytogenes. BMC Genom. 23, 335. doi: 10.1186/s12864-022-08437-4
Palumbo, M., Sosa, E., Castello, F., Schottlender, G., Serral, F., Turjanski, A., et al. (2022). Integrating diverse layers of omic data to identify novel drug targets in Listeria monocytogenes. Front. Drug Discov. 2, 969415. doi: 10.3389/fddsv.2022.969415
Pang, X., Wu, Y., Liu, X., Wu, Y., Shu, Q., Niu, J., et al. (2022). The lipoteichoic acid-related proteins YqgS and LafA contribute to the resistance of Listeria monocytogenes to nisin. Microbiol. Spectr. 10, e02095–e02021. doi: 10.1128/spectrum.02095-21
Perez-Riverol, Y., Bai, J., Bandla, C., García-Seisdedos, D., Hewapathirana, S., Kamatchinathan, S., et al. (2022). The PRIDE database resources in 2022: a hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 50, D543–D552. doi: 10.1093/nar/gkab1038
Petersen, B., Petersen, T. N., Andersen, P., Nielsen, M., and Lundegaard, C. A. (2009). A generic method for assignment of reliability scores applied to solvent accessibility predictions. BMC Struct. Biol. 9, 1–10. doi: 10.1186/1472-6807-9-51
Piercey, M. J., Hingston, P. A., and Truelstrup Hansen, L. (2016). Genes involved in Listeria monocytogenes biofilm formation at a simulated food processing plant temperature of 15°C. Int. J. Food Microbiol. 223, 63–74. doi: 10.1016/j.ijfoodmicro.2016.02.009
Pinheiro, J., Lisboa, J., Pombinho, R., Carvalho, F., Carreaux, A., Brito, C., et al. (2018). MouR controls the expression of the Listeria monocytogenes Agr system and mediates virulence. Nucleic Acids Res. 46, 9338–9352. doi: 10.1093/nar/gky624
Pleitner, A. M., Trinetta, V., Morgan, M. T., Linton, R. L., and Oliver, H. F. (2014). Transcriptional and phenotypic responses of Listeria monocytogenes to chlorine dioxide. Appl. Environ. Microbiol. 80, 2951–2963. doi: 10.1128/AEM.00004-14
Quereda, J. J., Morón-García, A., Palacios-Gorba, C., Dessaus, C., García-del-Portillo, F., Pucciarelli, M. G., et al. (2021). Pathogenicity and virulence of Listeria monocytogenes: a trip from environmental to medical microbiology. Virulence 12, 2509–2545. doi: 10.1080/21505594.2021.1975526
Rahman, O., Cummings, S. P., Harrington, D. J., and Sutcliffe, I. C. (2008). Methods for the bioinformatic identification of bacterial lipoproteins encoded in the genomes of Gram-positive bacteria. World J. Microb. Biot. 24, 2377–2382. doi: 10.1007/s11274-008-9795-2
Rappsilber, J., Ishihama, Y., and Mann, M. (2003). Stop and go extraction tips for matrix-assisted laser desorption/ionization, nanoelectrospray, and LC/MS sample pretreatment in proteomics. Anal. Bioanal. Chem. 75, 663–670. doi: 10.1021/ac026117i
Ribeiro, V. B., Mujahid, S., Orsi, R. H., Bergholz, T. M., Wiedmann, M., Boor, K. J., et al. (2014). Contributions of σB and PrfA to Listeria monocytogenes salt stress under food relevant conditions. Int. J. Food Microbiol. 177, 98–108. doi: 10.1016/j.ijfoodmicro.2014.02.018
Rodríguez-López, P., Rodríguez-Herrera, J. J., Vázquez-Sánchez, D., and López Cabo, M. (2018). Current knowledge on Listeria monocytogenes biofilms in food-related environments: Incidence, resistance to biocides, ecology and biocontrol. Foods 7, 85. doi: 10.3390/foods7060085
Sanchez, L. (2001). TCA Protein Precipitation Protocol. Available online at: http://www.its.caltech.edu/bjorker/TCA_ppt_protocol.pdf
Savijoki, K., Miettinen, I., Nyman, T. A., Kortesoja, M., Hanski, L., Varmanen, P., et al. (2020). Growth mode and physiological state of cells prior to biofilm formation affect immune evasion and persistence of Staphylococcus aureus. Microorganisms 8, 106. doi: 10.3390/microorganisms8010106
Schauer, K., Stolz, J., Scherer, S., and Fuchs, T. M. (2009). Both thiamine uptake and biosynthesis of thiamine precursors are required for intracellular replication of Listeria monocytogenes. J. Bacteriol. 191, 2218–2227. doi: 10.1128/JB.01636-08
Shevchenko, A., Wilm, M., Vorm, O., and Mann, M. (1996). Mass spectrometric sequencing of proteins silver-stained polyacrylamide gels. Anal. Chem. 68, 850–858. doi: 10.1021/ac950914h
Shevchuk, O. (2021). “Chapter 11 - Targeted proteomics for rapid and sensitive detection of foodborne pathogens,” in Comprehensive Foodomics, A. Cifuentes ed. (Elsevier Inc.), 123–136. doi: 10.1016/B978-0-08-100596-5.22841-2
Shin, J.-H., Kim, J., Kim, S.-M., Kim, S., Lee, J.-C., Ahn, J.-M., et al. (2010). σB-dependent protein induction in Listeria monocytogenes during vancomycin stress: σB and vancomycin stress in L. monocytogenes. FEMS Microbiol. Lett. 308, 94–100. doi: 10.1111/j.1574-6968.2010.01998.x
Srivastava, S., Sharma, S. K., Srivastava, V., and Kumar, A. (2021). Proteomic exploration of Listeria monocytogenes for the purpose of vaccine designing using a reverse vaccinology approach. Int. J. Pept. Res. Ther. 27, 779–799. doi: 10.1007/s10989-020-10128-1
Sureka, K., Choi, P. H., Precit, M., Delince, M., Pensinger, D. A., Huynh, T. N., et al. (2014). The cyclic dinucleotide c-di-amp is an allosteric regulator of metabolic enzyme function. Cell 158, 1389–1401. doi: 10.1016/j.cell.2014.07.046
Tanui, C. K., Benefo, E. O., Karanth, S., and Pradhan, A. K. (2022). A machine learning model for food source attribution of Listeria monocytogenes. Pathogens 11, 691. doi: 10.3390/pathogens11060691
Tiong, H., and Muriana, P. (2016). RT-qPCR analysis of 15 Genes encoding putative surface proteins involved in adherence of Listeria monocytogenes. Pathogens 5, 60. doi: 10.3390/pathogens5040060
Van Der Veen, S., Van Schalkwijk, S., Molenaar, D., De Vos, W. M., Abee, T., and Wells-Bennik, M. H. J. (2010). The SOS response of Listeria monocytogenes is involved in stress resistance and mutagenesis. Microbiology 156, 374–384. doi: 10.1099/mic.0.035196-0
Way, S. S., Thompson, L. J., Lopes, J. E., Hajjar, A. M., Kollmann, T. R., Freitag, N. E., et al. (2004). Characterization of flagellin expression and its role in Listeria monocytogenes infection and immunity. Cell. Microbiol. 6, 235–242. doi: 10.1046/j.1462-5822.2004.00360.x
Witte, C. E., Whiteley, A. T., Burke, T. P., Sauer, J.-D., Portnoy, D. A., and Woodward, J. J. (2013). Cyclic di-AMP is critical for Listeria monocytogenes growth, cell wall homeostasis, and establishment of infection. MBio 4, e00282–e00213. doi: 10.1128/mBio.00282-13
Wu, J., McAuliffe, O., and O'Byrne, C. P. (2023). Manganese uptake mediated by the NRAMP-type transporter MntH is required for acid tolerance in Listeria monocytogenes. Int. J. Food Microbiol. 399, 110238. doi: 10.1016/j.ijfoodmicro.2023.110238
Yao, J., Ericson, M. E., Frank, M. W., and Rock, C. O. (2016). Enoyl-Acyl Carrier protein reductase I (FabI) is essential for the intracellular growth of Listeria monocytogenes. Infect. Immun. 84, 3597–3607. doi: 10.1128/IAI.00647-16
Zeng, Z., Boeren, S., Bhandula, V., Light, S. H., Smid, E. J., Notebaart, R. A., et al. (2021). Bacterial microcompartments coupled with extracellular electron transfer drive the anaerobic utilization of ethanolamine in Listeria monocytogenes. mSystems 6, e01349–e01320. doi: 10.1128/mSystems.01349-20
Zeng, Z., Smid, E. J., Boeren, S., Notebaart, R. A., and Abee, T. (2019). Bacterial microcompartment-dependent 1,2-propanediol utilization stimulates anaerobic growth of Listeria monocytogenes EGDe. Front. Microbiol. 10, 2660. doi: 10.3389/fmicb.2019.02660
Zhang, H., Wang, J., Chang, Z., Liu, X., Chen, W., Yu, Y., et al. (2021). Listeria monocytogenes contamination characteristics in two ready-to-eat meat plants from 2019 to 2020 in Shanghai. Front. Microbiol. 12, 729114. doi: 10.3389/fmicb.2021.729114
Zhang, Y., Gonzalez-Gutierrez, G., Legg, K. A., Walsh, B. J. C., Pis Diez, C. M., Edmonds, K. A., et al. (2022). Discovery and structure of a widespread bacterial ABC transporter specific for ergothioneine. Nat. Comm. 13, 7586. doi: 10.1038/s41467-022-35277-3
Keywords: Listeria monocytogenes, proteome, bioinformatic, STRING database, stressors, pathogen, virulence
Citation: D'Onofrio F, Schirone M, Krasteva I, Tittarelli M, Iannetti L, Pomilio F, Torresi M, Paparella A, D'Alterio N and Luciani M (2023) A comprehensive investigation of protein expression profiles in L. monocytogenes exposed to thermal abuse, mild acid, and salt stress conditions. Front. Microbiol. 14:1271787. doi: 10.3389/fmicb.2023.1271787
Received: 02 August 2023; Accepted: 19 September 2023;
Published: 09 October 2023.
Edited by:
Didier Cabanes, Universidade do Porto, PortugalReviewed by:
Ricardo Monteiro, Universidade do Porto, PortugalJialun Wu, University of Galway, Ireland
Copyright © 2023 D'Onofrio, Schirone, Krasteva, Tittarelli, Iannetti, Pomilio, Torresi, Paparella, D'Alterio and Luciani. 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: Maria Schirone, bXNjaGlyb25lJiN4MDAwNDA7dW5pdGUuaXQ=