ORIGINAL RESEARCH article
New Genes Involved in Mild Stress Response Identified by Transposon Mutagenesis in Lactobacillus paracasei
- Université de Bourgogne Franche-Comté, AgroSup Dijon, PAM UMR A 02.102, Dijon, France
Lactic acid bacteria (LAB) are associated with various plant, animal, and human niches and are also present in many fermented foods and beverages. Thus, they are subjected to several stress conditions and have developed advanced response mechanisms to resist, adapt, and grow. This work aimed to identify the genes involved in some stress adaptation mechanisms in LAB. For this purpose, global reverse genetics was applied by screening a library of 1287 Lactobacillus paracasei transposon mutants for mild monofactorial stresses. This library was submitted independently to heat (52°C, 30 min), ethanol (170 g.L−1, 30 min), salt (NaCl 0.8 M, 24 h), acid (pH 4.5, 24 h), and oxidative (2 mM H2O2, 24 h) perturbations which trigger mild monofactorial stresses compatible with bacterial adaptation. Stress sensitivity of mutants was determined either by evaluating viability using propidium iodide (PI) staining, or by following growth inhibition through turbidity measurement. The screening for heat and ethanol stresses lead respectively to the identification of 63 and 27 genes/putative promoters whose disruption lead to an increased sensitivity. Among them, 14 genes or putative promoters were common for both stresses. For salt, acid and oxidative stresses, respectively 8, 6, and 9 genes or putative promoters were identified as essential for adaptation to these unfavorable conditions, with only three genes common to at least two stresses. Then, RT-qPCR was performed on selected stress response genes identified by mutant screenings in order to evaluate if their expression was modified in response to stresses in the parental strain. Eleven genes (membrane, transposase, chaperone, nucleotide and carbohydrate metabolism, and hypothetical protein genes) were upregulated during stress adaptation for at least two stresses. Seven genes, encoding membrane functions, were upregulated in response to a specific stress and thus could represent potential transcriptomic biomarkers. The results highlights that most of the genes identified by global reverse genetics are specifically required in response to one stress and that they are not differentially transcribed during stress in the parental strain. Most of these genes have not been characterized as stress response genes and provide new insights into the adaptation of lactic acid bacteria to their environment.
Lactobacillus casei/paracasei is one of the most emblematic groups of lactic acid bacteria (LAB), probably because it was one of the first studied and marketed for its benefits on health in the form of dairy products (Saxelin et al., 2005). Besides, strains of this group are nonstarter LAB in ripened cheeses, where they contribute to flavor development. Comparative genomics demonstrated that this versatile group is highly adaptable to various niches, thus it can be present and active in plethora of environments such as dairy and plant products, oral, intestinal, and reproductive tracts of humans and animals (Cai et al., 2010; Douillard et al., 2013; Duar et al., 2017). Moreover, among Lactobacillus genus, L. casei, paracasei, and rhamnosus species are widely used and studied for their probiotic function. For example, it was demonstrated that some strains can inhibit pathogenic bacteria (Ingrassia et al., 2005) and prevent diarrhea (de Roos and Katan, 2000). Considering the importance of these strains for industrial application as well as for health effects, many publications focused on stress response in L. casei/paracasei group, in particular when they are exposed to acid (Broadbent et al., 2010), cold (Beaufils et al., 2007), and bile stresses (Wu et al., 2010; Alcantara and Zuniga, 2012; Hamon et al., 2012). These data are discussed by Hosseini Nezhad et al. (2015).
Although LAB are fastidious because they possess numerous auxotrophies, and that only a few strains are genetically manipulable, many studies have been dedicated to the comprehension of the high adaptability of LAB to environmental stresses, guided by fundamental issues and potential industrial applications. Physiological approaches have highlighted the phenomenon of stress adaptation, i.e., the transient improvement of stress resistance when LAB are submitted to a first stress of mild intensity (Andrade-Linares et al., 2016). This phenomenon, also called priming stress, is a key factor for the development of active starters or probiotics and for their preservation. Then, as transcriptomics and proteomics have progressed, genetic determinants of stress response have been identified (De Angelis and Gobbetti, 2004). Some genes were qualified as “general” stress response genes whereas others were stress-specific. Some of the most studied proteins are Hsp (Heat shock protein; Sugimoto et al., 2008), UspA (Universal stress protein; Kvint et al., 2003; Gury et al., 2009), and those of the ATR (Acid Tolerant Response; Cotter and Hill, 2003). In addition, some genes can be functionally determinant for stress adaptation without, however, exhibiting transcriptomic or proteomic changes. Indeed, comparative genome analyses demonstrated that genes for niche adaptation are mainly carried by plasmids and adaptation islands on the chromosome (Cai et al., 2010; Smokvina et al., 2013). Also, targeted inactivation of genes (allowing reverse genetics) brought new information on LAB stress response. All knowledges generated by these studies were reviewed in a comprehensive and detailed manner by van de Guchte et al. (2002) and Papadimitriou et al. (2015, 2016).
More global approaches such as random mutagenesis or gene expression library have made possible global reverse genetics but they have been only slightly used so far, probably because of the difficulty to develop them. The main reason is that random mutagenesis approaches need the building of a transposon mutant library and then the screening of each mutant impaired for the function of interest in order to identify the corresponding gene(s) responsible for the phenotype. This requires transposon vectors tailored for the bacterial species and is often a limiting factor for LAB which can be difficult to transform, and for which the selective antibiotics panel is limited. To overcome this limitation, the Pjunc-TpaseIS1223 system has been developed specifically for Lactobacillus genus (Licandro-Seraut et al., 2012). We succeeded in this approach by applying transposon mutagenesis in the frame of the adaptation of Lactobacillus pentosus to olive brine, which is a multifactorial stress, and identified five brine sensitive mutants (Perpetuini et al., 2013, 2016). We have demonstrated the strength of this approach by uncovering many new genetic determinants during the early stage of L. paracasei establishment at the intestinal level, when bacteria have to cope with numerous changes of their environment (Licandro-Seraut et al., 2014). For this study, thanks to the sequencing of each transposon target of the mutant library, we have set up a non-redundant library of 1,110 transposon mutants in which the transposon target was identified in a coding region for each mutant (Scornec et al., 2014).
In the present work, we revisit the issue of LAB stress response by a global reverse genetics approach. We took advantage of the annotated random transposon library of L. paracasei cited above. In order to have a fine resolution, we have chosen to apply simple (monofactorial) stress of mild intensities, compatible with an adaptation of the LAB population. We defined this type of perturbation as “mild stress.” For this purpose, we have implemented suitable screening strategies to study separately the responses to five of the major components that can change when L. paracasei is under multifactorial environmental stresses: heat, salt, oxidative, ethanol, and acidic stresses. The cross-checking of the screening results and the transcriptional analysis of selected genes allowed us to identify new genetic determinants of stress response and to qualify them for a “general stress” or “specific stress” response.
Materials and Methods
Strains and Growth Conditions
Wild-type (WT) L. paracasei ATCC 334 (CIP 107868, Institut Pasteur Collection; Judicial Commission of the International Committee on Systematics of Bacteria, 2008) and its corresponding mutants obtained by random transposon mutagenesis (Licandro-Seraut et al., 2014) were grown statically at 37°C in MRS (Difco). The antibiotic used was 5 μg.mL−1 erythromycin (Erm) during mutant preculture. The mutants corresponded to the 1,100 genic mutants already described (Licandro-Seraut et al., 2012) and to 177 new mutants selected among the 8,053 previously annotated mutants (Scornec et al., 2014) because their transposon was inserted at <100 bp upstream of a gene for which no genic mutant had been found. For each mutant, the putative inactivated function was assigned thanks to the genome annotation of this strain (Cai et al., 2010).
Determination of Stress Intensities Based on WT Viability
WT L. paracasei was grown until stationary phase because it was the most convenient for the mutant screening. The culture (1 ml) was centrifuged (4,000 g, 4 min, 25°C) and the pellet was suspended with one volume of buffer modified for the stress. Phosphate buffer (pH 6.5) was supplemented with 1 M to 4 M NaCl for salt stress, with 450 μM to 350 mM H2O2 for oxidative stress and with 150–200 g.L−1 ethanol for ethanol stress. For acid stress, the buffer was a 0.2 M KCl/0.2 M HCl solution with pH ranging from 2.0 to 2.4. Exposure to heat stress (48–60°C) was performed in phosphate buffer in a C1000™ Thermal Cycler (Biorad). WT viability after a 30-min exposure to stress was determined by plate counts (three biological replicates).
Viability Screening by Propidium Iodide
Mutant cultures in stationary phase (20 μl) were put in a white 96-well PCR plate (Biorad), centrifuged, pellets were suspended in 100 μl of stress buffer. Phosphate buffer was supplemented with 170 g.L−1 ethanol for ethanol stress. Exposure to heat stress (52°C) was performed in prewarmed phosphate buffer and in a C1000™ Thermal Cycler (Biorad). After a 30-min incubation, bacterial pellets were suspended in 100 μl of 20 μM propidium iodide (PI, Sigma-Aldrich). All centrifugations were done at 4,000 g for 10 min and all suspensions were performed by vortexing (1 min at 1,550 rpm in a Mix mate, Eppendorf). Fluorescence (λem/λexc: 490/635 nm) was measured in a plate reader (Beckman Coulter) to calculate the viability percent with respect to a calibration range of known proportions of “dead” cells of L. paracasei WT (heat inactivated, 30 min at 80°C). Two biological replicates were made. When necessary, a third biological replicate was made to confirm the phenotype. Mutants were qualified as sensitive if their viability was less or equal to the WT viability minus twice the standard deviation for the two biological replicates. Corresponding genes were aligned using BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) against all bacteria to determine their specificity.
Growth Inhibition Screening by Turbidimetry
Mutant cultures in stationary phase were inoculated at 1/20 (v/v) (OD600 nm ≈ 0.2) in 200 μl of modified MRS as specified hereafter. MRS was supplemented with 0.8M NaCl for salt stress, 2 mM H2O2 for oxidative stress and acidified with HCl until pH 4.5 for acid stress. Microplates were incubated 24 h at 37°C and vortexed (1,000 rpm 1 min in a Mix Mate, Eppendorf) to read the final OD600 nm in a plate reader (Beckman Coulter). Mutants were qualified as sensitive when OD600 nm (24 h) was 1.5-fold lower than that of the WT final OD600 nm (two biological replicates). Corresponding genes were aligned using BLAST against all bacteria to determine their specificity. An alternative method was also used for selected sensitive mutants. After a 30-min stress in modified MRS, bacterial pellets were suspended in MRS to follow the growth at 37°C during 7 h. Mutants were considered sensitive when their growth percentage at 7 h in stress condition compared to the control condition, was inferior to the WT growth percentage minus the value of two standard derivations (two biological repeats).
Growth Inhibition Determined by Diffusion Test Assay for Oxidative Stress
Mutant cultures in stationary phase of growth were inoculated at 1/20 in 20 mL of MRS agar before solidification in a Petri dish. Disks impregnated with 10 μl of 3.5 M H2O2 were placed on MRS agar. After a 24-hincubation at 37°C, inhibition diameters were measured and mutants were considered as sensitive when their inhibition diameters were significantly higher than that of the WT (Student test, p < 0.05, six biological replicates).
For transcriptomic analysis, cells in mid-exponential growth phase (OD600 between 0.6 and 0.8) were centrifuged and cell pellets were suspended with one volume of the dedicated media. Thereby, in function of the applied stress, MRS was supplemented with 1 mM H2O2, 1 M NaCl, 150 g.L−1 ethanol, or acidified with HCl to pH 3.0. Phosphate buffer was maintained in a water bath at 50°C for heat stress (four biological replicates). Total RNA isolation, cDNA synthesis, and qPCR was made as previously described (Licandro-Seraut et al., 2008) using TRI Reagent (Sigma Aldrich), DNase I (Roche), iScript™ Reverse Transcription Supermix (Bio-Rad), and SsoAdvanced™ Universal SYBR® Green Supermix (Bio-Rad). Primers were designed by using Primer3Plus (Untergasser et al., 2007) (Table S1). Quantitative PCR (qPCR) were performed using a CFX96 Touch™ Real-Time PCR Detection System (Bio-Rad) in triplicate, in a 20 μl-reaction mixture. Cq (threshold value) calculation was determined by a regression model of the CFX Manager™ Software. Gene expression was calculated using 2−ΔΔCT method (Livak and Schmittgen, 2001). In order to select appropriate reference genes, 10 housekeeping genes (fusA, ileS, lepA, leuS, mutL, pcrA, pyrG, recA, recG, and rpoB) were tested with all experimental conditions and analyzed using the CFX Manager™ Software. The genes recG and rpoB were selected as the references because they displayed the lowest M-values (0.57) and coefficients of variation (0.20), meaning that they have the most stable expression in the tested stress conditions.
Development of the Mutant Screening Strategy
In addition to the 1,110 mutants that have already been screened in a previous study (Licandro-Seraut et al., 2014), we have selected 177 new mutants in this work. These new mutants allowed us to increase the potential of information from the screening. They were selected because their transposon was inserted in the putative promoter of a gene for which no genic mutant was present in the initial genic library. Thus, a library containing a total number of 1,287 L. paracasei mutants was screened for the ability to resist to mild thermic, ethanol, acid, osmotic, and oxidative stresses in comparison with their parental strain (wild type, WT) as described in section Materials and Methods. The global screening strategy is shown in Figure 1.
Figure 1. Organizational chart of the screening strategy to identify genetic factors involved in L. paracasei response to five monofactorial mild stresses: heat, ethanol, salt, acid, and oxidative.
Initially, we aimed to measure mutant viability using propidium iodide (PI) in 96-well plates for all types of stress since PI presents the advantage of a direct detection and a high sensitivity. For each stress, a range of conditions (intensity, time) was applied to WT resting cells (in phosphate buffer) to determine mild stress conditions, i.e., the highest intensity condition for which L. paracasei WT was at least 80% viable. A correlation between viabilities obtained by PI method and plate count agar method (CFU counts) was made for each stress with the WT. It demonstrated that PI staining was a relevant method for heat and ethanol stresses, but not for salt, acid, and oxidative stresses. Viabilities were 91 ± 7% (CFU counts) vs. 70 ± 1% (PI counts) for a 2 M NaCl stress, 94 ± 5% (CFU counts) and 76 ± 1% (PI counts) for a pH 2.4 stress, 41 ± 9% (UFC counts) and 84 ± 0% (PI counts) for a 3.5 mM H2O2 stress.
The following parameters were retained for the screening: 52°C during 30 min for heat stress and 170 g.L−1 ethanol during 30 min for ethanol stress. The WT viability was 82% (±1%) for heat stress and 84% (±2%) for ethanol stress (six biological repeats).
As an alternative, the three other stresses were applied at the beginning of the WT growth and their sensitivity was determined according to the growth inhibition triggered by the stress. Mild stress conditions were defined as the highest intensity condition for which the WT OD600 nm (24 h) was less or equal to the 80% of OD600 nm (24 h) in absence of stress. The retained parameters were pH 4.5, 0.8 M NaCl, and 2 mM H2O2 for acid, salt and oxidative stresses respectively. When WT strain was exposed to these stresses, OD600 nm (24 h) was 1.33 (±0.03) for acid stress, 1.00 (±0.05) for salt stress, 1.31 (±0.03) for oxidative stress whereas it reached 1.73 (±0.05) in the control condition (MRS). For oxidative stress, as standard derivations were high for some mutants, the 23 selected mutants were subjected to diffusion test with 3.5 M H2O2 in order to confirm and quantify the sensitivity of each mutant for oxidative stress. Mutants were considered sensitive if their inhibition diameters were significantly higher than that of the WT (p < 0.05, six biological repeats).
Identification of Genes Involved in Responses to Heat and Ethanol Mild Stresses by Viability Screening
The mutant library was subjected independently to heat stress and ethanol stress as described before. Viability analysis revealed 104 sensitive mutants including 71 genic and 33 intergenic transposon mutants (Figure 2). Mean viabilities ranged from 49 to 80% for sensitive mutants (Table S2). A large majority of mutants, 61 genic and 29 intergenic, were sensitive to only one stress condition (heat or ethanol). These disrupted genes encode for various predicted functions such as transcriptional regulation, amino acid transport and metabolism, carbohydrate transport and metabolism, membrane biogenesis, and hypothetical proteins. The most sensitive mutants were LSEI_0656 (P) (putative DNA-entry nuclease), LSEI_0824 (P) (putative hypothetical protein), and LSEI_2289 (putative hydrolase) for heat stress, LSEI_0221 (putative D-alanyl-D-alanine carboxypeptidase) and LSEI _2733 (putative L-xylulose-5-phosphate 3-epimerase) for ethanol stress, with viability lower than 70%. The 14 mutants sensitive to the two stresses are disrupted in cell wall/membrane genes such as putative ABC transporters, PTS system, membrane enzymatic function, but also in putative transposase genes, in a putative nucleotide and carbohydrate metabolism gene, and in putative genes of unknown function (Figure 2). LSEI_0040 (putative hypothetical protein) was the most sensitive mutant in this category. These genes could be attributed to general stress response.
Figure 2. Venn diagram of L. paracasei mutants sensitive to heat (52°C, 30 min) and ethanol (170 g.L−1, 30 min) mild stresses. Corresponding viabilities using propidium iodide labeling are indicated only for mutants sensitive to the 2 stresses (at least 2 biological repeats). (P), intergenic insertion mutant (putative promoter). The predictive function of each gene is reported in Table 5.
Identification of Genes Involved in Response to Salt, Acid, and Oxidative Mild Stresses by Growth Inhibition Screening
The mutant library was subjected independently to salt stress (0.8 M NaCl), acid stress (pH 4.5), and oxidative stress (2 mM H2O2) during a 24 h-growth as described before. The screening resulted in detection of eight mutants for salt stress, six mutants for acid stress, and 23 mutants for oxidative stress (Table S3). For the sensitive mutants, DO600 (24 h) were 1.5–7 times lower than that of the WT. In a second round, oxidative mutants were screened on MRS plates with an H2O2 gradient. Among the 23 mutants initially selected, nine mutants displayed a significant sensitivity (p < 0.05) and were selected as sensitive (Figure 3). The mean final inhibition diameter ranged from 31 to 33 mm compared to 29 mm for the WT. Most of mutants were sensitive to only one of the three stresses. Their disrupted genes encode for various predicted functions such as Zn protease, D-alanine activating enzyme, ABC transporter for acid stress, ABC transporter, DNA binding response regulator, phosphoglycerate mutase for salt stress and ABC transporter, cysteine desulfurase, nucleotidase for oxidative stress. The most sensitive mutants were LSEI_1468 (putative ribonucleotide reductase), LSEI_1565 (DnaK), LSEI_2540 (putative Zn protease) with an OD600 (24 h) lower than 50% of that of the WT. Only three genes seemed to be involved in response to several stresses: LSEI_1468 and LSEI_2540 for the three stresses and LSEI_1565 for salt and acid stresses. So genetic responses to salt, acid, and oxidative mild stresses seem rather specific of the type of stress applied.
Figure 3. Venn diagram of L. paracasei mutants sensitive to salt (0.8 M NaCl, 24 h), acid (pH 4.5, 24 h), and/or oxidative (3.5 M H2O2, 24 h) mild stresses, their corresponding OD600 24 h and their inhibition diameters (ID). (P), putative promoter. The predictive function of each gene is reported in Table 5.
Transcriptomic Analysis of Genes Involved in Acid, Salt, or Oxidative Stresses
In order to expand the analysis, transcriptomic analysis of L. paracasei WT was performed on the 18 genes identified after the growth inhibition screening. WT was subjected to 15-min mild stresses (1 M NaCl for salt stress, pH 3.0 for acid stress, and 1 mM H2O2 for oxidative stress) and a culture without stress was used as the reference (Table 1). In parallel, for the 18 corresponding mutants, growth inhibition was measured after 30-min salt, acid and oxidative stresses of same intensity than for transcriptomics (Table 2). The growth of nine mutants was more inhibited than the WT after a stress application, whereas the growth of the 18 mutants was comparable to the WT in absence of stress. Among the 18 genes, six were upregulated for the stress for which the corresponding mutant was sensitive: LSEI_1468 (putative ribonucleotide reductase) and LSEI_2262 (putative hypothetical protein) for salt and oxidative stresses, LSEI_0794 (putative D-alanine-activating enzyme) and LSEI_1565 (DnaK) for salt stress, LSEI_2540 (putative Zn protease) for oxidative stress, LSEI_0938 (putative phosphate ABC transporter) for acid stress. The most up-regulated genes -with values over 3.0- were LSEI_0938 for acid stress (3.35 ± 0.69), LSEI_0794 and LSEI_2262 for salt stress (3.43 ± 0.93 and 3.99 ± 1.12) and LSEI_1565 for oxidative stress (3.21 ± 1.01).
Table 1. Differentially expressed genes in L. paracasei after acid, salt or oxidative stress among the 18 genes for which a corresponding mutant has been identified as sensitive for at least one of these stress conditions.
Table 2. Growth percentages after acid, salt, or oxidative stress, for the 18 mutants identified as sensitive for at least one of these stress conditions.
Transcriptomic Analysis of Genes Involved in Both Heat and Ethanol Stresses
To study general stress response, the identification of the sensitive mutants for both heat and ethanol stresses (14 mutants) was complemented by transcriptomic analysis on the corresponding genes for all stresses. Transcriptomic analysis was performed on the WT subjected to 15-min stresses (150 g.L−1 for ethanol stress, 50°C for heat stress, 1 M NaCl for salt stress, pH 3.0 for acid stress, and 1 mM H2O2 for oxidative stress; Table 3). Among the 14 genes, eight were upregulated for the stress for which the corresponding mutant was sensitive (Table 4): LSEI_0281 (putative cell wall-associated hydrolase) and LSEI_1241 (putative transposase) for heat stress, LSEI_1884 (putative peptide ABC transporter) and LSEI_2320 (putative mucin binding protein) for salt stress, LSEI_1111 and LSEI_2304 (putative hypothetical proteins) for heat and salt stress, LSEI_0389 (putative pseudo gene) for heat and acid stress and LSEI_0631 for all stress conditions. For this last gene, encoding a putative PTS system, expressions were over 4.0 for heat (5.47 ± 3.07), salt (4.45 ± 1.11), acid (4.68 ± 0.77), and oxidative (10.20 ± 3.01) stresses. LSEI_1884 gene, encoding a putative peptide ABC transporter, was strongly upregulated only for salt stress (11.11 ± 3.43). Four genes, LSEI_0040 (putative hypothetical protein), LSEI_0389 (putative pseudo gene), LSEI_0613 and LSEI_1241 (putative transposases) were more upregulated for heat stress (value >3.0). LSEI_2399 (putative metal-dependent membrane protease) was upregulated especially for heat (5.43 ± 1.96) and acid (4.51 ± 1.07) stresses.
Table 3. Differentially expressed genes in L. paracasei after heat, ethanol, acid, salt, or oxidative stress, for the 14 genes for which a corresponding mutants has been identified as sensitive to both heat and ethanol stress conditions.
Table 4. Growth percentages after heat, ethanol, acid, salt, or oxidative mild stress, for the 14 mutants identified as sensitive to both heat and ethanol stress conditions.
Development of the Mutant Screening
In this study, a transposon mutant library of 1,287 mutants was exposed to five monofactorial mild stresses in order to identify genetic determinants involved in bacterial adaptation. The first part of the work was to develop a suitable screening method which should be easy and fast but also and mostly accurate. We observed that viability determinations using propidium iodide as the marker of mortality, although very easy and fast to carry out, were aberrant for studying salt, acid, and oxidative stresses. Thus we developed growth inhibition screening for these stresses. Previous studies have also reported that CFU counts and fluorescence staining methods were sometimes not correlated in particular when adherent bacteria were harvested using sonication, which could lead to a higher number of colonies on plates (Hannig et al., 2007; Tawakoli et al., 2013). For oxidative stress, as observed by Zotta et al. (2012), bacteria were not stained with PI. Maybe this kind of stress does not immediately affect membrane integrity.
Comparison of Mutagenesis and Transcriptomic Approaches for Stress Response
To complete the results obtained after mutant screenings, transcriptomic analyzes were performed on selected genes. Among the 32 genes required in stress response and selected for transcriptomics, 10 genes were not differentially expressed in the stressed WT. This result strengthens the scope of global reverse genetics compared to global transcriptomics since many genes may be essential for a function but constitutively expressed. Other studies have demonstrated that some genes can be required for stress resistance without exhibiting transcriptomic changes. It is particularly the case for genes encoding monocistronic transcriptional regulators (Tran et al., 2008). Recently, Price et al. (2013) have shown that there was only little correlation between gene upregulation and mutant sensitivity in their bacterial models.
Genes Involved in General Stress Response
Mutant library screening led to the identification of 20 genes and 5 putative promoters involved in multiple stress response (at least two stresses; Table 5). For 15 of them, we did not find any literature reporting their role in stress response (Tables S2, S3) thus they can be considered as new multiple stress response determinants. Most are specific of the L. casei/paracasei/rhamnosus group. On the contrary, the putative ABC transporter LSEI_1884 and the putative response regulator of a two component system LSEI_1041 caught our attention because they are well conserved among Gram-positive bacteria. On the other hand, four genes, LSEI_0631 (putative beta-glucoside-specific PTS system), LSEI_1111 (putative hypothetical protein), LSEI_0389 (putative pseudo gene), and LSEI_2320 (putative mucin binding protein) can be considered as important general stress determinants because involved in response to 4 or 5 stress conditions (Table 5). The putative PTS system beta-glucoside-specific transporter (LSEI_0631) was the only one upregulated for all stress conditions, which perfectly matches the screening results. The PTS mannose system is one of the main sugar transporters for LAB. Several authors have reported a link between PTS systems and stress response for Lactobacillus genus. L. plantarum mutant for rpoN displayed an impaired expression of the mannose PTS operon, which increased their sensitivity to peroxide (Stevens et al., 2010). The authors have assumed that the suppression of this transporter led to the diminution of glucose capture and energy that may explain the susceptibility to oxidative stress. A L. casei resistant mutant for acid stress obtained by adaptation had a higher PTS activity than the WT after 1 h at pH 5.0 (Wu et al., 2012). The authors suggest that this results may lead to higher level of ATP via glycolysis which allows maintaining intracellular pH. LSEI_1111 had no function assigned. Transcriptomic analysis of LSEI_0389, encoding a putative pseudo-gene, demonstrated that it was transcribed and upregulated for heat and acid stresses. LSEI_2320 possesses three mucin binding domains (MucBP) which are involved in cell adhesion (Munoz-Provencio et al., 2012; Yamasaki-Yashiki et al., 2017). It has never been reported as determinant of stress response before.
Moreover, we have observed that 11 of the 25 identified genes were also involved in multiple stress response at the transcriptomic level (genes in red in Table 6). This is the case of the chaperone dnaK (LSEI_1565) which is known to be involved in various LAB stress responses (van de Guchte et al., 2002). Also, the expression of genes involved in membrane modification such as LSEI_0281 (cell wall-associated hydrolase) and LSEI_2399 (predicted metal-dependent membrane protease) was increased and highlighted that cell envelope feature was essential for multiple stress response. The phosphoglucose isomerase gene (LSEI_1126, putative glucose-6-phosphate isomerase) was upregulated for heat, ethanol, and oxidative stresses. This gene was upregulated during exposure of L. plantarum to manganese (Tong et al., 2017).
Table 6. Summary of genes involved in general and specific stress response in function of the genetic approach used (mutagenesis or transcriptomics).
Our results have shown that the salt and oxidative stresses increased the expression of a putative ribonucleotide reductase gene (LSEI_1468, RNR). On the contrary, under bile stress, a decrease of RNR expression was observed in other Lactobacillus (Burns et al., 2010; Koskenniemi et al., 2011). The hypothesis would be that this decrease was probably a consequence of a reduced growth rate and a decrease of DNA replication. Mutagenesis and transcriptomics highlight the involvement of two putative transposase genes (LSEI_1241 and LSEI_0613) in several stresses. Similarly, 12 other transposase genes were upregulated during acid adaptation at pH 4.5 in L. casei (Broadbent et al., 2010). Authors have assumed that stimulation of transposition could provide an evolutionary advantage to the host.
Genes Involved in Specific Stress Responses
Among the 1,287 screened mutants, 118 were sensitive to at least one stress condition. Comparison of stress response indicates that genes are generally specific to one stress condition (93 genes, Table 5). All of these specific genes are relevant as indicators of a particular stress (i.e., stress markers). A large part of them have never been reported as required for stress response. Hence, to our knowledge, we have identified 47 L. casei/paracasei/rhamnosus group genes, 22 Lactobacillus specific genes and 23 Gram-positive bacteria genes as new genes for specific stress responses (Tables S2, S3). Remarkably, we have reported the involvement of LSEI_1945 and LSEI_1951 which encode putative phage proteins, a functional family that has never been reported to act during stress response.
Identifications of Stress Markers
Cellular biomarkers are useful tools to predict the impact of environmental changes on bacterial robustness and survival. Generally, applications of stress markers are rather based on upregulations at the transcriptomic and/or proteomic level. For instance, in L. casei, a work has investigated the potential markers of bile tolerance by comparison of the proteomic profiles of six strains with different bile tolerance levels (Hamon et al., 2012). We have analyzed transcription of some stress-specific genes on L. paracasei WT when subjected to monofactorial mild stresses. We have selected genes as stress markers because they were upregulated specifically for one stress condition and their disruption led to an increased sensitivity to the same stress. Three genes (LSEI_0794, LSEI_1884 and LSEI_2320) were identified for NaCl stress, two genes (LSEI_0040 and LSEI_0613) for heat stress, one gene (LSEI_0938) for acid stress, and one gene (LSEI_2540) for oxidative stress (genes in gray in the Table 6). Interestingly, all biomarkers identified except LSEI_0613 (a transposase, IS30 family) are related to cell wall functions or localization (LSEI_0040 is a gene of unknown function with several transmembrane domains). LSEI_0794 (putative D-alanine-activating enzyme), _1884 (putative peptide ABC transporter), _0613 (putative transposase), _0938 (putative phosphate ABC transporter), and _2540 (putative ATP dependant Zn protease) present a higher potential as biomarker since they are not restricted to the L. paracasei group genomes.
Some functions have already been associated to stresses. Firstly, the dlt operon (the LSEI_0793–0797ortholog), responsible for D-alanine incorporation into teichoic acids, is upregulated during bile and heat stresses in L. plantarum (Xie et al., 2004; Bron et al., 2006). Secondly, in L. plantarum, ftsH mutation (LSEI_2540 ortholog, ATP dependent Zn protease) reduced the growth rate in physiological conditions, and the growth defect became more important under stress (Bove et al., 2012). Thirdly, several studies have reported the involvement of ABC transporters in LAB in response to bile or osmotic stresses (Bouvier et al., 2000; Hamon et al., 2011; Alcantara and Zuniga, 2012; Wang et al., 2015). Peptide transport systems do not only play a role in cell nutrition, but are also involved in various signaling processes (Detmers et al., 2001). In this study, LSEI_0938 and LSEI_1884 encode putative ABC transporters, the latter represents a strong marker of salt stress since the measured upregulation was the most important of this work (11.11 ± 3.43).
This work explored the screening of a transposon mutant library of L. paracasei to identify genes required to face five monofactorial mild stresses. It was complemented by targeted transcriptomics on selected genes. Firstly, 118 genes whose mutation led to a sensitive mutant were identified. These genes were generally specific to one stress condition and most of them had never been reported to be involved in stress response or adaptation. They deserve to be studied in a more thorough way to explore LAB genetic stress response. At least half of them are conserved in Lacobacillus genus. Secondly, five of the seven genes identified as potential stress biomarkers are related to the cell wall and could be used as selective stress reporters of membrane injury. Finally, some of these genes whose expression was not regulated by stress would have escaped a transcriptomics or proteomics screening. The identified stress determinants and biomarkers contribute to the comprehension of stress responses and adaptation, including stresses encountered in microbial processes and in food matrices. They could be targets for a better control of growth performance and functional properties of starters in food. Hence, the selection of high performance strains could be achieved by screening strains for the presence of the target gene. In addition, the expression of the target gene could be monitored and would indicate the intensity of stress during food process. In the context of improving starters by stress adaptation, since most of the stress response genes identified here are stress-specific, we assume that the nature of the mild stress (or priming stress) has to be carefully selected in function of the perturbations during the process.
AP, HS, HL, and J-FC: study conception and design. AP: experimentation and acquisition of data. AP and HL: analysis and interpretation of data. AP: drafting of manuscript. AP, HS, HL, and J-FC: critical revision. All authors read and approved the final manuscript.
This work was supported by the Regional Council of Bourgogne—Franche Comté and the Fonds Européen de DEveloppement Régional (FEDER). AP thesis was supported by the Ministry of Higher Education, Research, and Innovation (MESRI).
Conflict of Interest Statement
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.
We thank the other members of the laboratory for their assistance, in particular, Christine Rojas for her skillful technical assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00535/full#supplementary-material
Beaufils, S., Sauvageot, N., Mazé, A., Laplace, J.-M., Auffray, Y., Deutscher, J., et al. (2007). The cold shock response of Lactobacillus casei: relation between HPr phosphorylation and resistance to freeze/thaw cycles. J. Mol. Microbiol. Biotechnol. 13, 65–75. doi: 10.1159/000103598
Bouvier, J., Bordes, P., Romeo, Y., Fourçans, A., Bouvier, I., Gutierrez, C., et al. (2000). Characterization of OpuA, a glycine-betaine uptake system of Lactococcus lactis. J. Mol. Microbiol. Biotechnol. 2, 199–205. Available online at: https://www.caister.com/jmmb/tocs/v2n2toc.html
Bove, P., Capozzi, V., Garofalo, C., Rieu, A., Spano, G., and Fiocco, D. (2012). Inactivation of the ftsH gene of Lactobacillus plantarum WCFS1: effects on growth, stress tolerance, cell surface properties and biofilm formation. Microbiol. Res. 167, 187–193. doi: 10.1016/j.micres.2011.07.001
Broadbent, J. R., Larsen, R. L., Deibel, V., and Steele, J. L. (2010). Physiological and transcriptional response of Lactobacillus casei ATCC 334 to acid Stress. J. Bacteriol. 192, 2445–2458. doi: 10.1128/JB.01618-09
Bron, P. A., Molenaar, D., de Vos, W. M., and Kleerebezem, M. (2006). DNA micro-array-based identification of bile-responsive genes in Lactobacillus plantarum. J. Appl. Microbiol. 100, 728–738. doi: 10.1111/j.1365-2672.2006.02891.x
Burns, P., Sanchez, B., Vinderola, G., Ruas-Madiedo, P., Ruiz, L., Margolles, A., et al. (2010). Inside the adaptation process of Lactobacillus delbrueckii subsp. lactis to bile. Int. J. Food Microbiol. 142, 132–141. doi: 10.1016/j.ijfoodmicro.2010.06.013
Cai, H., Thompson, R., Budinich, M. F., Broadbent, J. R., and Steele, J. L. (2010). Genome sequence and comparative genome analysis of Lactobacillus casei: insights into their niche-associated evolution. Genome Biol. Evol. 1, 239–257. doi: 10.1093/gbe/evp019
de Roos, N. M., and Katan, M. B. (2000). Effects of probiotic bacteria on diarrhea, lipid metabolism and carcinogenesis: a review of papers published between 1988 and 1998. Am. J. Clin. Nutr. 71, 405–411. doi: 10.1093/ajcn/71.2.405
Douillard, F. P., Ribbera, A., Jarvinen, H. M., Kant, R., Pietila, T. E., Randazzo, C., et al. (2013). Comparative genomic and functional analysis of Lactobacillus casei and Lactobacillus rhamnosus strains marketed as probiotics. Appl. Environ. Microbiol. 79, 1923–1933. doi: 10.1128/AEM.03467-12
Duar, R. M., Lin, X. B., Zheng, J., Martino, M. E., Grenier, T., Pérez-Muñoz, M. E., et al. (2017). Lifestyles in transition: evolution and natural history of the genus Lactobacillus. FEMS Microbiol. Rev. 41, S27–S48. doi: 10.1093/femsre/fux030
Gury, J., Seraut, H., Tran, N. P., Barthelmebs, L., Weidmann, S., Gervais, P., et al. (2009). Inactivation of PadR, the repressor of the phenolic acid stress response, by molecular interaction with Usp1, a universal stress protein from Lactobacillus plantarum, in Escherichia coli. Appl. Environ. Microbiol. 75, 5273–5283. doi: 10.1128/AEM.00774-09
Hamon, E., Horvatovich, P., Bisch, M., Bringel, F., Marchioni, E., Aoudé-Werner, D., et al. (2012). Investigation of biomarkers of bile tolerance in Lactobacillus casei using comparative proteomics. J. Proteome Res. 11, 109–118. doi: 10.1021/pr200828t
Hamon, E., Horvatovich, P., Izquierdo, E., Bringel, F., Marchioni, E., Aoudé-Werner, D., et al. (2011). Comparative proteomic analysis of Lactobacillus plantarum for the identification of key proteins in bile tolerance. BMC Microbiol. 11:63. doi: 10.1186/1471-2180-11-63
Hannig, C., Hannig, M., Rehmer, O., Braun, G., Hellwig, E., and Al-Ahmad, A. (2007). Fluorescence microscopic visualization and quantification of initial bacterial colonization on enamel in situ. Arch. Oral Biol. 52, 1048–1056. doi: 10.1016/j.archoralbio.2007.05.006
Ingrassia, I., Leplingard, A., and Darfeuille-Michaud, A. (2005). Lactobacillus casei DN-114 001 inhibits the ability of adherent-invasive Escherichia coli isolated from Crohn's disease patients to adhere to and to invade intestinal epithelial cells. Appl. Environ. Microbiol. 71, 2880–2887. doi: 10.1128/AEM.71.6.2880-2887.2005
Judicial Commission of the International Committee on Systematics of Bacteria (2008). The type strain of Lactobacillus casei is ATCC 393, ATCC 334 cannot serve as the type because it represents a different taxon, the name Lactobacillus paracasei and its subspecies names are not rejected and the revival of the name “Lactobacillus zeae” contravenes Rules 51b (1) and (2) of the International Code of Nomenclature of Bacteria. Opinion 82. Int. J. Syst. Evol. Microbiol. 58, 1764–1765. doi: 10.1099/ijs.0.2008/005330-0
Koskenniemi, K., Laakso, K., Koponen, J., Kankainen, M., Greco, D., Auvinen, P., et al. (2011). Proteomics and transcriptomics characterization of bile stress response in probiotic Lactobacillus rhamnosus GG. Mol. Cell. Proteomics 10:M110.002741. doi: 10.1074/mcp.M110.002741
Licandro-Seraut, H., Brinster, S., van de Guchte, M., Scornec, H., Maguin, E., Sansonetti, P., et al. (2012). Development of an efficient in vivo system (Pjunc-TpaseIS1223) for random transposon mutagenesis of Lactobacillus casei. Appl. Environ. Microbiol. 78, 5417–5423. doi: 10.1128/AEM.00531-12
Licandro-Seraut, H., Gury, J., Tran, N. P., Barthelmebs, L., and Cavin, J. F. (2008). Kinetics and intensity of the expression of genes involved in the stress response tightly induced by phenolic acids in Lactobacillus plantarum. J. Mol. Microbiol. Biotechnol. 14, 41–47. doi: 10.1159/000106081
Licandro-Seraut, H., Scornec, H., Pédron, T., Cavin, J.-F., and Sansonetti, P. J. (2014). Functional genomics of Lactobacillus casei establishment in the gut. Proc. Natl. Acad. Sci. U.S.A. 111, E3101–E3109. doi: 10.1073/pnas.1411883111
Munoz-Provencio, D., Rodriguez-Diaz, J., Collado, M. C., Langella, P., Bermudez-Humaran, L. G., and Monedero, V. (2012). Functional analysis of the Lactobacillus casei BL23 sortases. Appl. Environ. Microbiol. 78, 8684–8693. doi: 10.1128/AEM.02287-12
Papadimitriou, K., Alegría, Á., Bron, P. A., de Angelis, M., Gobbetti, M., Kleerebezem, M., et al. (2016). Stress physiology of lactic acid bacteria. Microbiol. Mol. Biol. Rev. 80, 837–890. doi: 10.1128/MMBR.00076-15
Papadimitriou, K., Zoumpopoulou, G., Foligné, B., Alexandraki, V., Kazou, M., Pot, B., et al. (2015). Discovering probiotic microorganisms: in vitro, in vivo, genetic and omics approaches. Front. Microbiol. 6:58. doi: 10.3389/fmicb.2015.00058
Perpetuini, G., Pham-Hoang, B. N., Scornec, H., Tofalo, R., Schirone, M., Suzzi, G., et al. (2016). In Lactobacillus pentosus, the olive brine adaptation genes are required for biofilm formation. Int. J. Food Microbiol. 216, 104–109. doi: 10.1016/j.ijfoodmicro.2015.10.002
Perpetuini, G., Scornec, H., Tofalo, R., Serror, P., Schirone, M., Suzzi, G., et al. (2013). Identification of critical genes for growth in olive brine by transposon mutagenesis of Lactobacillus pentosus C11. Appl. Environ. Microbiol. 79, 4568–4575. doi: 10.1128/AEM.01159-13
Price, M. N., Deutschbauer, A. M., Skerker, J. M., Wetmore, K. M., Ruths, T., Mar, J. S., et al. (2013). Indirect and suboptimal control of gene expression is widespread in bacteria. Mol. Syst. Biol. 9:660. doi: 10.1038/msb.2013.16
Saxelin, M., Tynkkynen, S., Mattila-Sandholm, T., and de Vos, W. M. (2005). Probiotic and other functional microbes: from markets to mechanisms. Curr. Opin. Biotechnol. 16, 204–211. doi: 10.1016/j.copbio.2005.02.003
Scornec, H., Tichit, M., Bouchier, C., Pédron, T., Cavin, J.-F., Sansonetti, P. J., et al. (2014). Rapid 96-well plates DNA extraction and sequencing procedures to identify genome-wide transposon insertion sites in a difficult to lyse bacterium: Lactobacillus casei. J. Microbiol. Methods 106, 78–82. doi: 10.1016/j.mimet.2014.08.001
Smokvina, T., Wels, M., Polka, J., Chervaux, C., Brisse, S., Boekhorst, J., et al. (2013). Lactobacillus paracasei comparative genomics: towards species pan-genome definition and exploitation of diversity. PLoS ONE 8:e68731. doi: 10.1371/journal.pone.0068731
Stevens, M. J. A., Molenaar, D., de Jong, A., de Vos, W. M., and Kleerebezem, M. (2010). Involvement of the mannose phosphotransferase system of Lactobacillus plantarum WCFS1 in peroxide stress tolerance. Appl. Environ. Microbiol. 76, 3748–3752. doi: 10.1128/AEM.00073-10
Sugimoto, S., and Abdullah-Al-Mahin Sonomoto, K. (2008). Molecular chaperones in lactic acid bacteria: physiological consequences and biochemical properties. J. Biosci. Bioeng. 106, 324–336. doi: 10.1263/jbb.106.324
Tawakoli, P. N., Al-Ahmad, A., Hoth-Hannig, W., Hannig, M., and Hannig, C. (2013). Comparison of different live/dead stainings for detection and quantification of adherent microorganisms in the initial oral biofilm. Clin. Oral Invest. 17, 841–850. doi: 10.1007/s00784-012-0792-3
Tong, Y., Zhai, Q., Lu, W., Tian, F., Zhao, J., Zhang, H., et al. (2017). New insights in integrated response mechanism of Lactobacillus plantarum under excessive manganese stress. Food Res. Int. 102, 323–332. doi: 10.1016/j.foodres.2017.10.014
Tran, N. P., Gury, J., Dartois, V., Nguyen, T. K. C., Seraut, H., Barthelmebs, L., et al. (2008). Phenolic acid-mediated regulation of the padC gene, encoding the phenolic acid decarboxylase of Bacillus subtilis. J. Bacteriol. 190, 3213–3224. doi: 10.1128/JB.01936-07
Untergasser, A., Nijveen, H., Rao, X., Bisseling, T., Geurts, R., and Leunissen, J. A. M. (2007). Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 35, W71–W74. doi: 10.1093/nar/gkm306
van de Guchte, M., Serror, P., Chervaux, C., Smokvina, T., Ehrlich, S. D., and Maguin, E. (2002). Stress responses in lactic acid bacteria. Antonie Van Leeuwenhoek 82, 187–216. doi: 10.1023/A:1020631532202
Wang, G., Li, D., Ma, X., An, H., Zhai, Z., Ren, F., et al. (2015). Functional role of oppA encoding an oligopeptide-binding protein from Lactobacillus salivarius Ren in bile tolerance. J. Ind. Microbiol. Biotechnol. 42, 1167–1174. doi: 10.1007/s10295-015-1634-5
Wu, C., Zhang, J., Chen, W., Wang, M., Du, G., and Chen, J. (2012). A combined physiological and proteomic approach to reveal lactic-acid-induced alterations in Lactobacillus casei Zhang and its mutant with enhanced lactic acid tolerance. Appl. Microbiol. Biotechnol. 93, 707–722. doi: 10.1007/s00253-011-3757-6
Wu, R., Sun, Z., Wu, J., Meng, H., and Zhang, H. (2010). Effect of bile salts stress on protein synthesis of Lactobacillus casei Zhang revealed by 2-dimensional gel electrophoresis. J. Dairy Sci. 93, 3858–3868. doi: 10.3168/jds.2009-2967
Xie, Y., Chou, L., Cutler, A., and Weimer, B. (2004). DNA Macroarray profiling of Lactococcus lactis subsp. lactis IL1403 gene expression during environmental stresses. Appl. Environ. Microbiol. 70, 6738–6747. doi: 10.1128/AEM.70.11.6738-6747.2004
Yamasaki-Yashiki, S., Sawada, H., Kino-oka, M., and Katakura, Y. (2017). Analysis of gene expression profiles of Lactobacillus paracasei induced by direct contact with Saccharomyces cerevisiae through recognition of yeast mannan. Biosci. Microbiota Food Health 36, 17–25. doi: 10.12938/bmfh.BMFH-2016-015
Zotta, T., Guidone, A., Tremonte, P., Parente, E., and Ricciardi, A. (2012). A comparison of fluorescent stains for the assessment of viability and metabolic activity of lactic acid bacteria. World J. Microbiol. Biotechnol. 28, 919–927. doi: 10.1007/s11274-011-0889-x
Keywords: lactic acid bacteria, transposon mutants, mild stresses, stress response genes, bacterial adaptation
Citation: Palud A, Scornec H, Cavin J-F and Licandro H (2018) New Genes Involved in Mild Stress Response Identified by Transposon Mutagenesis in Lactobacillus paracasei. Front. Microbiol. 9:535. doi: 10.3389/fmicb.2018.00535
Received: 12 January 2018; Accepted: 08 March 2018;
Published: 23 March 2018.
Edited by:Vittorio Capozzi, University of Foggia, Italy
Reviewed by:Teresa Zotta, Consiglio Nazionale delle Ricerche (CNR), Italy
Rosanna Tofalo, Università di Teramo, Italy
Copyright © 2018 Palud, Scornec, Cavin and Licandro. 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 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.