Skip to main content

Transcriptome-phenotype matching analysis of how nitrogen sources influence Lacticaseibacillus rhamnosus tolerance to heat stress and oxidative stress

Abstract

Background

Spray drying is the most cost-effective production method for lactic acid bacteria starters, but heat and oxidative stresses result in low survival rates. The heat stress and oxidative stress tolerance of Lacticaseibacillus rhamnosus cultured in tryptone-free MRS (NP-MRS) broth was much stronger than that in MRS or tryptone-free MRS broth supplemented with phenylalanine (Phe-MRS). Here, multiple transcriptome-phenotype matching was performed on cells cultured in NP-MRS, MRS and Phe-MRS broths to reveal the mechanism by which nitrogen sources influence L. rhamnosus tolerance to heat stress and oxidative stress.

Results

Compared with cells cultured in NP-MRS broth, 83 overlapping differentially expressed genes (DEGs) were downregulated by either tryptone or phenylalanine. The overlapping DEGs were mainly classified into carbohydrate metabolism and membrane transport pathways, which are often repressed by glucose during carbon catabolite repression (CCR). In the presence of glucose, the heat stress or oxidative stress tolerance of L. rhamnosus hsryfm 1301 was not strengthened by supplementation with secondary carbohydrates. Replacing glucose with mannose, fructose or ribose improved the heat stress and oxidative stress tolerance of L. rhamnosus hsryfm 1301 (5 to 46-fold).

Conclusions

Alleviation of CCR might be a reason for the resistance of L. rhamnosus hsryfm 1301 to heat stress and oxidative stress in a low-nitrogen environment. The survival rate of L. rhamnosus during spray drying will hopefully be improved by relieving CCR. It is a new discovery that nitrogen sources influence CCR in L. rhamnosus.

Background

Probiotics provide a number of health benefits. Viable count is essential for lactic acid bacteria (LAB) to achieve fermentation and beneficial effects. However, various stresses that LAB encounter during the production, transport and storage processes can seriously decrease their viable count. Several attempts, such as freezing, drying and microencapsulation, have been made to improve the viability of probiotics in different food products during their production until the time of consumption [1].

An organism survives stress conditions through multiple physiological responses. Comparative genomics, DNA microarrays, transcriptome and proteome profiling, and genotype–phenotype matching have been employed to identify genetic factors regulated during single stress exposure [2]. Regulators such as HrcA, CtsR, GlnR, and CcpA, and the classical heat shock proteins GroEL and DnaK, play an important role in the stress resistance of LAB [3,4,5].

However, there are disadvantages in the study of single stress tolerance. LAB usually encounter various stress conditions during production and storage processes, but the study of single stress tolerance can only provide information for treating one form of stress. On the other hand, hundreds of functional genes and regulatory elements respond to one stress in LAB [6, 7], so it is difficult to screen novel major genes by single stressome-phenotype matching. Therefore, cross-adaptation which results from stress response overlaps occurring among different stress treatments, has been applied to improve resistance and reveal the underlying mechanisms of stress responses [8,9,10].

Spray drying is cost-effective and highly flexible and has been applied in the production of LAB powder [11, 12]. Heat stress and oxidative stress are the major stresses resulting in a low survival rate of starter cultures during spray drying [13, 14]. In our previous study, mutual cross-adaptation of Lacticaseibacillus rhamnosus hsryfm 1301 to oxidative stress and heat stress was found [15]. Transcriptome-phenotype matching based on RNA sequencing was implemented, and it was revealed that the high proportion of transcriptional homogenization, especially the decrease in abundance of nitrogen source transporter and metabolism enzyme genes was a reason for the cross-adaptation [10]. Tryptone-free MRS broth universally enhanced the heat stress and oxidative stress tolerance (H-OST) of L. rhamnosus cells, with the survival rates under heat stress and oxidative stress being increased by 130 and 40 times, respectively (OD600, 2.0). The spray-drying survival rate of L. rhamnosus hsryfm 1301 cultured in tryptone-free MRS broth rose to 75% (from 30%), and the spray-dried powder also exhibited higher survival in simulated gastrointestinal digestion. After supplementing tryptone-free MRS broth with phenylalanine, isoleucine, glutamate, valine, histidine, or tryptophan, the H-OST of L. rhamnosus hsryfm 1301 exhibited a sharp drop [16]. Amino acids often played a positive role in stress resistance of LAB [9, 17]. Therefore, it is of interesting to determine how the decrease in abundance of nitrogen sources improves the survival of L. rhamnosus hsryfm 1301 under heat stress and oxidative stress.

In the present study, L. rhamnosus hsryfm 1301 was cultured in MRS, tryptone-free MRS and tryptone-free MRS supplemented with phenylalanine to enable transcriptome-phenotype matching and screen gene transcription profiles which were correlated to survival characteristics, revealing the mechanism of the cotolerance of L. rhamnosus to heat stress and oxidative stress through the downregulation of nitrogen metabolism.

Results

Correlation of expression patterns in MRS, NP-MRS and Phe-MRS broths

The reads obtained from sequencing were mapped to the genome of L. rhamnosus hsryfm 1301. The pairwise correlation between each biological replicate was shown based on normalized expression levels. A high reproducibility among biological replicates from the same medium was observed, indicated by a Pearson correlation > 0.98. Moreover, when comparing the results among different media, a Pearson correlation < 0.97 can be obtained (Fig. 1a). Both the Pearson correlation and the PCA (Fig. 1b) of RNA-sequencing indicated the high reproducibility among technical replicates of RNA-seq and the differences in transcription profiles induced by tryptone and Phe. Interestingly, although the H-OST of cells from the Phe-MRS were more similar to those from the MRS, the transcription profile correlation between the Phe-MRS group and the NP-MRS group was higher than that between the Phe-MRS group and the MRS group (Fig. 1).

Fig. 1
figure 1

Correlation of expression patterns in different media. a Pairwise correlation of different biological replicates from the MRS, NP-MRS and Phe-MRS groups. The color intensities (scale in the side bar) and the numbers indicate the degree of pairwise correlation. b PCA of RNA-seq

Transcriptional regulation overlaps in response to tryptone and Phe

In our previous study, it was found that the H-OST of L. rhamnosus hsryfm 1301 cells cultured in MRS broth and Phe-MRS broth was significantly lower than those in NP-MRS broth. Overlapping transcriptional regulation in response to tryptone and Phe would be helpful to understand this phenomenon. To facilitate data presentation, the NP-MRS group was described as a control. Compared with the NP-MRS group, 90 genes were upregulated (log2FC > 1) in the MRS group, and 14 genes were upregulated (log2FC > 1) in the Phe-MRS group (Fig. 2a). There was no upregulation overlap between the MRS vs. NP-MRS comparison and the Phe-MRS vs. NP-MRS comparison. There were 206 genes downregulated (log2FC < -1) in the MRS group. Thirty-seven genes were downregulated (log2FC < -1) by Phe, among which 29 genes were overlapping DEGs (log2FC < -1) responding to both tryptone and Phe (Fig. 2b), and the remaining 8 genes were also slightly downregulated (log2(2/3) < log2FC < -1) by tryptone.

Fig. 2
figure 2

Overlapping DEGs of L. rhamnosus hsryfm 1301 between the MRS group vs. NP-MRS group comparison and the Phe-MRS group vs. NP-MRS group comparison. a |log2FC|> 1, upregulation. b |log2FC|> 1, downregulation. c FC < 2/3, downregulation

To find more overlapping DEGs, the downregulated genes were also analyzed under the standard of |log2FC|> 0.585 (FC < 2/3). There were 671 genes downregulated (log2FC < -0.585) in cells cultured in MRS broth and 112 genes were downregulated (log2FC < -0.585) by Phe, among which 83 genes were overlapping DEGs (log2FC < -0.585) responding to both tryptone and Phe (Fig. 2c). It was suggested that the lower H-OST of L. rhamnosus hsryfm 1301 in the MRS and Phe-MRS groups resulted from the genes downregulated by tryptone and Phe.

KEGG analysis of overlapping DEGs in response to tryptone and Phe

Under the standard of |log2FC|> 1, 22 overlapping DEGs were involved in carbohydrate metabolism (10 DEGs), membrane transport (ABC transporters, 9 DEGs), drug resistance (1 DEG), cell motility (1 DEG), and glycan biosynthesis and metabolism (1 DEG) (Fig. 3a). All 9 overlapping membrane transport DEGs were related to carbohydrate transport. The 10 overlapping carbohydrate metabolism DEGs and 9 overlapping membrane transport DEGs were distributed in the metabolism pathways of 6 common sugars.

Fig. 3
figure 3

KEGG pathway analysis of the overlapping DEGs. a |log2FC|> 1. b FC < 2/3 (|log2FC|> 0.585)

Under the standard of |log2FC|> 0.585, 83 overlapping DEGs were enriched 63 times (the functions of some DEGs were unknown and some DEGs were enriched in more than 1 pathway) and involved in 16 KEGG pathways (Fig. 3b). The enriched overlapping DEGs were mainly involved in carbohydrate metabolism (30 DEGs), membrane transport (ABC transporters, 15 DEGs), glycan biosynthesis and metabolism (3 DEGs), transport and catabolism (2 DEGs), and the endocrine system (2 DEGs) (Fig. 3b). All 15 overlapping membrane transport DEGs were presumed to encode proteins transporting carbohydrates according to NCBI NR annotation. The 30 overlapping carbohydrate metabolism DEGs and 15 overlapping membrane transport DEGs were distributed in the metabolism pathways of more than 10 common carbohydrates.

Overlapping DEGs related to classical stress tolerance

Although the H-OST of cells from the Phe-MRS and MRS groups were far weaker than those from the NP-MRS group, only a few overlapping DEGs were found to be classical stress tolerance genes. One heat shock response (HSR) gene, hsp20 was downregulated significantly by both tryptone and Phe, while groL and sigI/σ70 were downregulated by only Phe (Fig. 4a). Several genes considered to be relevant to the oxidative shock response (OSR), such as pox/spxB and the ascorbic acid operon, were downregulated by both tryptone and Phe. The transcription levels of most classical stress tolerance genes, such as regulators (ctsR, hrcA, and ccpA), chaperonins (dnaK and dnaJ), DNA repair genes (recA, uvrB, and uvrC) and proteases (clpP, clpB, and ftsH), were not influenced by tryptone or Phe. The regulation of classical HSR and OSR genes was unable to explain the cotolerance of L. rhamnosus to heat stress and oxidative stress in a low nitrogen source environment.

Fig. 4
figure 4

Overlapping DEGs related to classical stress tolerance. a Classical stress tolerance genes. b Genes in the ascorbic acid operon. a−.cValues within each column with different superscripts are significantly different (P < 0.05)

Overlapping DEGs related to sugar utilization and pyruvate flow

Most of the overlapping DEGs were related to the utilization of carbohydrates. Specifically, 4 genes in the mlt operon (mannose metabolism) were upregulated by more than threefold in a low nitrogen environment (Fig. 5a). Operons related to the utilization of fructose (Fig. 5b), ribose (Fig. 5c), trehalose (Fig. 5d), galactose (Fig. 5e), cellobiose (Fig. 5f) and isomaltose (Fig. 5g) were downregulated by both tryptone and Phe. These results suggested that L. rhamnosus hsryfm 1301 consumes a diversity of carbon sources under low nitrogen conditions. On the other hand, genes related to mixed acid fermentation, such as the pdh operon (acetate fermentation, Fig. 6a) and pfl operon (formate fermentation, Fig. 6a), were also downregulated by both tryptone and Phe. Moreover, most of these genes were similarly expressed in the MRS and NP-MRS groups. Both multiple carbon sources consumption and mixed acid fermentation were indicators of carbon catabolite repression (CCR).

Fig. 5
figure 5

Overlapping DEGs related to secondary carbohydrates. a Genes in the mannose operon. b Genes in the fructose operon. c Genes in the ribose operon. d Genes in the trehalose operon. e. Genes in the galactose operon. f Genes in the cellobiose operon. g Genes in the isomaltose operon. a−.cValues within each column with different superscripts are significantly different (P < 0.05)

Fig. 6
figure 6

Overlapping DEGs related to pyruvate flow. a Genes related to acetate fermentation. b Genes related to formate fermentation. a−.cValues within each column with different superscripts are significantly different (P < 0.05)

Heat stress and oxidative stress tolerance of L. rhamnosus hsryfm 1301 under conditions with different sugars

Several genes related to multiple carbon source consumption and mixed acid fermentation were upregulated in a low nitrogen environment. These genes were often regulated by glucose during CCR. To identify whether this response involved intake of secondary carbon sources or was just performance of the carbon catabolite derepression, the H-OST of L. rhamnosus hsryfm 1301 with different carbon sources were investigated. In Man-MRS, Fru-MRS, and Gal-MRS broths, L. rhamnosus hsryfm 1301 could grow as well as in Glc-MRS broth, in which the OD600 value reached 1.8 in 6 h (Fig. 7a). However, L. rhamnosus hsryfm 1301 grew much slower in Rib-MRS broth and its OD600 reached 0.45 at 6 h, so the survival rate of L. rhamnosus hsryfm 1301 in Glc-MRS broth at OD600 0.45 was also measured (Fig. 7b).

Fig. 7
figure 7

Growth characteristics of L. rhamnosus hsryfm 1301 under different carbon source conditions. a Growth curves. b Measurement points A, B and C were determined according to the Materials and methods section; Measurement points B and C, OD600, 0.45

In the case of heat shock, the survival rate of L. rhamnosus hsryfm 1301 (OD600 1.80) was 2.32% in Glc-MRS broth. After supplementation with mannose, fructose, galactose, or ribose (in Glc + Man-MRS, Glc + Fru-MRS, Glc + Gal-MRS, or Glc + Rib-MRS broths), the heat stress tolerance of L. rhamnosus hsryfm 1301 was not enhanced. However, when the glucose of Glc-MRS was replaced by mannose, galactose, or ribose (changing the broth to Man-MRS broth, Gal-MRS broth, or Rib-MRS broth, respectively), the thermotolerance of L. rhamnosus hsryfm 1301 was significantly enhanced (Fig. 8a). Specifically, the survival rate of L. rhamnosus hsryfm 1301 in Rib-MRS broth (OD600 0.45) reached 105% (Fig. 8b). In the case of oxidative shock, treatment with H2O2 killed 97.21% of L. rhamnosus hsryfm 1301 cells in Glc-MRS broth (OD600 1.80). The survival rate of L. rhamnosus hsryfm 1301 in Man-MRS, Fru-MRS, and Gal-MRS broths rose to 15.15%, 20.78%, and 13.63%, respectively (Fig. 8c). The survival rates of L. rhamnosus hsryfm 1301 at OD600 0.45 in Glc-MRS and Rib-MRS broths were 13.40% and 78.74%, respectively (Fig. 8d). However, aerotolerance was not improved in Glc + Man-MRS broth, Glc + Fru-MRS broth, Glc + Gal-MRS broth, or Glc + Rib-MRS broth, which were the broths supplemented with different sugars (Fig. 8c). These results suggested that sugar replacement, not sugar supplementation improved the H-OST of L. rhamnosus hsryfm 1301.

Fig. 8
figure 8

Heat stress and oxidative stress tolerance of L. rhamnosus hsryfm 1301 under different carbon source conditions. a Heat stress tolerance at OD600 1.80. b Heat stress tolerance at OD600 0.45. c Oxidative stress tolerance at OD600 1.80. d Oxidative stress tolerance at OD600 0.45. a−.cValues within each column with different superscripts are significantly different (P < 0.05)

Discussion

L. rhamnosus is one of the most commonly used probiotics [18,19,20]. To keep L. rhamnosus cells alive to perform health promotion functions, many studies have been carried out to improve its stress tolerance. Bacteria can withstand harsh conditions and sudden environmental changes through stress-sensing systems [10]. Omics data-phenotype matching approaches are practical methods for identifying novel genes that respond to stresses in L. rhamnosus [10, 21,22,23]. However, the stress tolerance of LAB is achieved by a complex regulatory network containing hundreds of functional genes and regulatory elements [6, 7], so it is difficult to screen novel genes by single stressome-phenotype matching. Compared with single stress tolerance studies, the study of the response of LAB to complex stresses was more powerful to solve the loss of activity caused by multiple stress conditions.

In our previous study, it was found that L. rhamnosus hsryfm 1301 was able to cotolerate heat stress and oxidative stress by decreasing the abundance of nitrogen sources [10]. The H-OST of L. rhamnosus hsryfm 1301 cultured in NP-MRS broth was significantly stronger than that from MRS broth or NP-MRS broth supplemented with phenylalanine, isoleucine, glutamate, valine, histidine, or tryptophan. The spray drying survival rate of L. rhamnosus hsryfm 1301 cultured in NP-MRS broth rose to 75% (from 30% in MRS broth) [16]. Therefore, these results helped reveal the underlying mechanisms of stress responses. In the present study, the mechanism of the cotolerance of L. rhamnosus to heat and oxidative stresses through the downregulation of nitrogen metabolism was investigated.

In the MRS and Phe-MRS groups, L. rhamnosus hsryfm 1301 performed similarly and had much weaker H-OST than that in the NP-MRS group. After culturing L. rhamnosus hsryfm 1301 in MRS broth, Phe-MRS broth, or NP-MRS broth, the stressomes in the MRS and Phe-MRS groups were compared with the stressomes in the NP-MRS group, and then the regulation overlaps facilitated the screening of novel gene functions in response to stress.

Compared with the NP-MRS group, the number of regulated genes (|log2FC|> 1) in the MRS group was much greater than that in the Phe-MRS group, which was consistent with the composition complexity of the media. There was no upregulation overlap between the MRS vs. NP-MRS comparison and the Phe-MRS vs. NP-MRS comparison. Among the 206 genes downregulated in the MRS group and 37 genes downregulated by Phe, 29 genes were overlapping DEGs. It was suggested that the weaker H-OST of L. rhamnosus hsryfm 1301 in the presence of tryptone and Phe were due to gene down-regulation, and less attention should be given to the upregulated DEGs. Through multiple transcriptome-phenotype matching, the scope in gene selection was decreased. If single transcriptome-phenotype matching is implemented, a large number of genes will be found. For example, 746 genes of Limosilactobacillus fermentum were regulated under bile stress [6], and 506 differently expressed proteins were identified in Lactiplantibacillus plantarum under cold stress [7]. Because the FC of many DEGs was close to 2.0, the standard that |log2FC|> 1 might result in omissions of underlying DEGs. So, the downregulated genes were also analyzed under the standard FC < 2/3. As a result, 83 genes were overlapping DEGs responding to both tryptone and Phe.

Although amino acids can influence the stress resistance of LAB [9, 17], almost no overlapping DEGs (FC < 2/3) belonged to the amino acid metabolism pathway, and only few overlapping DEGs were related to classical stress tolerance. Some common stress response genes, such as regulators (ctsR, hrcA, and ccpA), chaperonins (dnaK, dnaJ, groES, and groEL) and proteases (clp and ftsH), respond to heat stress in LAB [2]. A theoretical analysis showed that HSR repressors (HrcA and CtsR), alternative σ factors, protein repair systems (chaperonins and proteases) and proteins involved in DNA repair and export of toxic substances were related to HSR in Lacticaseibacillus casei and L. rhamnosus [4]. On the other hand, genes related to detoxification of ROS (sod, peroxidase gene, noxE, and so on) and DNA repair (recA, uvrA, and uvrB) are often found in the OSR of LAB [2]. However, only one HSR gene, hsp20, was downregulated significantly by both tryptone and Phe, and several genes considered to be relevant to OSR, such as pox/spxB [24] and the ascorbic acid operon, were downregulated by both tryptone and Phe. The classical stress tolerance genes groL and sigI/σ70 were only downregulated by Phe. The transcription of most classical stress tolerance genes was not influenced. Therefore, the regulation of classical HSR and OSR genes was insufficient to explain the cotolerance of L. rhamnosus to heat stress and oxidative stress in a low nitrogen source environment.

Interestingly, 71.4% of the enriched overlapping DEGs were related to transport and metabolism of carbohydrate, such as mannose, fructose, ribose, trehalose, galactose, cellobiose and isomaltose. Moreover, most of these genes were similarly expressed in the MRS and NP-MRS groups. Does a high concentration of secondary carbohydrates help L. rhamnosus hsryfm 1301 to survive heat stress and oxidative stress? Supplementation of Glc-MRS broth with mannose, galactose, fructose, or ribose did not strengthen the heat stress or oxidative stress tolerance. The use of a secondary carbon source is induced by the carbon source and repressed in carbon catabolite repression (CCR) that the use of secondary carbon sources is reduced in the presence of glucose [5]. CCR is important for competition in natural environments, as selection of the preferred carbon source is a major determining factor in growth rate and therefore competitive success with other microorganisms [25]. In LAB, catabolite control protein A (CcpA) is the global regulatory element in the control of CCR. CcpA also has a key role in the regulation of aerobic metabolism [26, 27]. In the present study, it was also found that genes related to acetate fermentation and formate fermentation were downregulated by both tryptone and Phe, which suggested that nitrogen sources influence CCR in L. rhamnosus.

CcpA has been shown to be involved in stress response mechanisms by exerting positive control on the dnaK and groESL operons in L. plantarum [28]. For exponential phase cells, inactivation of ccpA impaired both heat stress and cold/starvation stress, but increased oxidative stress tolerance in L. plantarum [5]. Replacing glucose with mannose, fructose or ribose not only strengthened the heat stress tolerance, but also improved the oxidative stress tolerance of L. rhamnosus hsryfm 1301, which indicated that relieving CCR could help L. rhamnosus hsryfm 1301 resist heat stress and oxidative stress. Unlike L. plantarum, the CCR in L. rhamnosus hsryfm 1301 was conducive to the tolerance of both heat stress and oxidative stress, which could explain the characteristic cotolerance of L. rhamnosus to heat stress and oxidative stress. To our knowledge, catabolite responsive elements (cre) and fructose-1,6-biphosphate are necessary for achieving CCR [25, 29], and nitrogen sources have not been reported to influence the function of CcpA. The finding that a low nitrogen source environment relieved CCR and then strengthened stress tolerance was interesting and worthy of further study.

Conclusion

The resistomes of L. rhamnosus hsryfm 1301 in the MRS, NP-MRS and Phe-MRS groups were identified. Most of the enriched overlapping DEGs were involved in secondary carbon source transport and metabolism. Replacing glucose by a secondary carbon source rather than adding a secondary carbon source helped L. rhamnosus hsryfm 1301 survive heat stress and oxidative stress, suggesting that a low nitrogen source environment strengthened the H-OST of L. rhamnosus hsryfm 1301 by relieving CCR. It is a new point of penetration to investigate the relationship between stress tolerance and carbon metabolism. The survival rate of L. rhamnosus during spray drying will hopefully be improved by relieving CCR.

Materials and methods

Bacterial strain and growth conditions

L. rhamnosus hsryfm 1301 was isolated from the gut of a centenarian from Bama, Guangxi Province, China, in 2013 and preserved in the China General Microbiological Culture Collection Center (CGMCC No. 8545) [10]. The strain was cultured in MRS broth (2% (v/v) inoculation) at 37 °C under static incubation.

MRS, NP-MRS (tryptone-free MRS) and Phe-MRS (addition of Phe to NP-MRS with a final concentration of 0.8 g L−1) broths were prepared according to our previous study [16]. Sugar supplement media were prepared with the addition of mannose, fructose, galactose, or ribose to the glycoprival MRS medium (Table 1 10 g L−1) [30].

Table 1 Concentrations of carbon sources in the modified MRS media

RNA isolation and sequencing

After overnight incubation, L. rhamnosus hsryfm 1301 was diluted 50-fold in 50 mL of fresh MRS broth, NP-MRS broth and Phe-MRS broth, and incubated at 37 °C until the optical density at 600 nm (OD600) reached 1.8–2.0. Each broth group included 3 replicates. The cells of L. rhamnosus hsryfm 1301 were washed with PBS buffer and harvested by centrifugation at 6000×g for 10 min at 4 °C. The resulting pellets were immediately frozen in liquid nitrogen.

For RNA extraction, total RNA was isolated using TRIzol Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA). Quality and integrity were determined using a Nano Drop spectrophotometer (Thermo Scientific, Wilmington, DE, USA) and a Bioanalyzer 2100 system (Agilent, Palo Alto, CA, USA). Library construction and RNA-Seq were performed by Shanghai Majorbio Bio-pharm Technology Co., Ltd. (Shanghai, China) as reported in a previous study [31].

The data generated from the Illumina platform were used for bioinformatics analysis. Analyses such as correlation of expression patterns were performed using the free online Majorbio Cloud Platform (www.majorbio.com) from Shanghai Majorbio Bio-pharm Technology Co., Ltd.

Mapping reads to the reference genome

The quality information of raw data in FASTQ format was calculated, and then the raw data were filtered using Cutadapt (v1.15) software [32]. Then, high-quality clean data were obtained by removing reads containing adapters, reads containing poly-N and low-quality reads. The remaining clean reads were mapped to the annotated genome of L. rhamnosus hsryfm 1301 using Bowtie 2 software [10, 33].

Differential expression analysis

Gene and isoform abundances from paired-end RNA-Seq data were quantified using RSEM (http://deweylab.github.io/RSEM/). To compare the gene expression levels of different genes and different samples, FPKM (fragments per kilobase of exon per million fragments mapped) was used to normalize the expression [34]. DESeq2 software was used to quantify all transcripts with differential expression [35]. Transcripts with |log2FoldChange|> 1 (log2FC) and P value < 0.05 were considered differentially expressed genes (DEGs).

KEGG analysis

Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.kegg.jp/) pathway analysis was performed on target genes of differentially expressed mRNAs. KOBAS 2.0 (http://kobas.cbi.pku.edu.cn) was used to identify statistically significantly enriched pathway using Fisher’s exact test. After multiple testing correction, pathways with a P value ≤ 0.05 were considered significantly enriched [36].

Growth of L. rhamnosus hsryfm 1301 in broths with different carbon sources

After overnight incubation, L. rhamnosus hsryfm 1301 was diluted 50-fold in 5 ml of each fresh modified MRS broth in which glucose was replaced by different carbon sources and incubated at 37 °C. The OD600 was measured at 2-h intervals.

Heat stress and oxidative stress tolerance detection

After overnight incubation, L. rhamnosus hsryfm 1301 was diluted 50-fold in 5 mL of fresh modified MRS broths and incubated at 37 °C. For the modified MRS broth in which L. rhamnosus hsryfm 1301 grew as rapidly as in Glc-MRS broth, tolerance detection was implemented when the OD600 reached 1.8–2.0 (at 6 h). For the modified MRS broth in which L. rhamnosus hsryfm 1301 grew more slowly than in the Glc-MRS broth, tolerance detection was implemented at 6 h (Point B, Fig. 9), and the tolerance in the Glc-MRS broth at the same OD600 was detected as a control (Point C, Fig. 9).

Fig. 9
figure 9

The method applied to define the OD600 value at which the stress tolerance was measured. Measurement point A: L. rhamnosus hsryfm 1301 cultured in Glc-MRS broth; time, 6 h; OD600, 1.8–2.0. Measurement point B: L. rhamnosus hsryfm 1301 cultured in modified MRS broth; time, 6 h; OD600, measured. Measurement point C: L. rhamnosus hsryfm 1301 cultured in Glc-MRS broth; time, at which the OD600 was equal to that of measurement point B

The sublethal heat stress and oxidative stress was determined according to the method described previously [15]. (1) For oxidative stress. The samples were supplemented with 4.0 mM H2O2, and then the samples were incubated at 37 °C for 1 h. (2) For heat stress. The samples were incubated at 55 °C in a water bath for 1 h. The colony-forming unit (CFU) counts of the treated samples were calculated according to a previously described drop plate technique [30].

Statistical analysis

Significant differences between two values were evaluated by unpaired Student’s t test (P < 0.05).

Nucleotide sequence accession numbers

The whole genome sequence of L. rhamnosus hsryfm 1301 was deposited into the National Center for Biotechnology Information (NCBI) under the accession numbers CP044228 (chromosome) and CP044229 (plasmid pGu1301). The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive in National Genomics Data Center [37, 38], China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences (accession number: CRA007167, experimental number: CRX455121-CRX455129), which are publicly accessible at https://bigd.big.ac.cn/gsa/browse/CRA007167.

Data Availability Statement

The data that support the findings of this study are openly available in the Genome Sequence Archive in National Genomics Data Center at https://bigd.big.ac.cn/gsa/browse/CRA007167

Abbreviations

LAB:

Lactic acid bacteria

H-OST:

Heat stress and oxidative stress tolerance

OD:

Optical density

DEGs:

Differentially expressed genes

KEGG:

Kyoto Encyclopedia of Genes and Genomes

FPKM:

Fragments per kilobase of exon per million fragments mapped

CFU:

Colony-forming unit

HSR:

Heat shock response

OSR:

Oxidative shock response

CCR:

Carbon catabolite repression

References

  1. Tripathi MK, Giri SK. Probiotic functional foods: survival of probiotics during processing and storage. J Funct Foods. 2014;9:225–41.

    Article  CAS  Google Scholar 

  2. Papadimitriou K, Alegría Á, Bron PA, de Angelis M, Gobbetti M, Kleerebezem M, Lemos JA, Linares DM, Ross P, Stanton C, et al. Stress physiology of lactic acid bacteria. Microbiol Mol Biol R. 2016;80(3):837–90.

    Article  CAS  Google Scholar 

  3. Gong L, Ren C, Xu Y. GlnR negatively regulates glutamate-dependent acid resistance in Lactobacillus brevis. Appl Environ Microbiol. 2020;86(7):e02615-2619.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  5. Zotta T, Ricciardi A, Guidone A, Sacco M, Muscariello L, Mazzeo MF, Cacace G, Parente E. Inactivation of ccpA and aeration affect growth metabolite production and stress tolerance in Lactobacillus plantarum WCFS1. Int J Food Microbiol. 2012;155(1–2):51–9.

    Article  CAS  Google Scholar 

  6. Ali SA, Singh P, Tomar SK, Mohanty AK, Behare P. Proteomics fingerprints of systemic mechanisms of adaptation to bile in Lactobacillus fermentum. J Proteomics. 2020;213: 103600.

    Article  CAS  Google Scholar 

  7. Liu S, Ma Y, Zheng Y, Zhao W, Zhao X, Luo T, Zhang J, Yang Z. Cold-stress response of probiotic Lactobacillus plantarum K25 by iTRAQ proteomic analysis. J Microbiol Biotechn. 2020;30(2):187–95.

    Article  CAS  Google Scholar 

  8. Gao X, Kong J, Zhu H, Mao B, Cui S, Zhao J. Lactobacillus Bifidobacterium and Lactococcus response to environmental stress: mechanisms and application of cross-protection to improve resistance against freeze-drying. J Appl Microbiol. 2022;132(2):802–21.

    Article  CAS  Google Scholar 

  9. Yang H, He M, Wu C. Cross protection of lactic acid bacteria during environmental stresses: Stress responses and underlying mechanisms. LWT. 2021;144: 111203.

    Article  CAS  Google Scholar 

  10. Zhang C, Gui Y, Chen X, Chen D, Guan C, Yin B, Pan Z, Gu R. Transcriptional homogenization of Lactobacillus rhamnosus hsryfm 1301 under heat stress and oxidative stress. Appl Microbiol Biot. 2020;104(6):2611–21.

    Article  CAS  Google Scholar 

  11. Agudelo-Chaparro J, Ciro-Velásquez HJ, Sepúlveda-Valencia JU, Pérez-Monterroza EJ. Microencapsulation of Lactobacillus rhamnosus ATCC 7469 by spray drying using maltodextrin whey protein concentrate and trehalose. Food Sci Technol Int. 2021;1:10820132211020620.

    Google Scholar 

  12. Moayyedi M, Eskandari MH, Rad AHE, Ziaee E, Khodaparast MHH, Golmakani MT. Effect of drying methods (electrospraying freeze drying and spray drying) on survival and viability of microencapsulated Lactobacillus rhamnosus ATCC 7469. J Funct Foods. 2018;40:391–9.

    Article  CAS  Google Scholar 

  13. Dijkstra AR, Setyawati MC, Bayjanov JR, Alkema W, van Hijum SA, Bron PA, Hugenholtz J. Diversity in robustness of Lactococcus lactis strains during heat stress oxidative stress and spray drying stress. Appl Environ Microbiol. 2014;80(2):603–11.

    Article  CAS  Google Scholar 

  14. Simpson PJ, Stanton C, Fitzgerald GF, Ross RP. Intrinsic tolerance of Bifidobacterium species to heat and oxygen and survival following spray drying and storage. J Appl Microbiol. 2005;99(3):493–501.

    Article  CAS  Google Scholar 

  15. Zhang C, Lu J, Yang D, Chen X, Huang Y, Gu R. Stress influenced the aerotolerance of Lactobacillus rhamnosus hsryfm 1301. Biotechnol Lett. 2018;40(4):729–35.

    Article  Google Scholar 

  16. Zhang C, Han Y, Gui Y, Wa Y, Chen D, Huang Y, Yin B, Gu R. Influence of nitrogen sources on the tolerance of Lacticaseibacillus rhamnosus to heat stress and oxidative stress. J Ind Microbiol Biot. 2022;1:20.

    Google Scholar 

  17. Fernández M, Zúñiga M. Amino acid catabolic pathways of lactic acid bacteria. Crit Rev Microbiol. 2006;32(3):155–83.

    Article  Google Scholar 

  18. Capurso L. Thirty years of Lactobacillus rhamnosus GG: a review. J Clin Gastroenterol. 2019;53(Suppl 1):S1-41.

    Article  CAS  Google Scholar 

  19. Wickens K, Barthow C, Mitchell EA, Kang J, van Zyl N, Purdie G, Stanley T, Fitzharris P, Murphy R, Crane J. Effects of Lactobacillus rhamnosus HN001 in early life on the cumulative prevalence of allergic disease to 11 years. Pediatric Allergy Immunol. 2018;29(8):808–14.

    Article  Google Scholar 

  20. Wa Y, Yin B, He Y, Xi W, Huang Y, Wang C, Guo F, Gu R. Effects of single probiotic- and combined probiotic-fermented milk on lipid metabolism in hyperlipidemic rats. Front Microbiol. 2019;10:1312.

    Article  Google Scholar 

  21. Koponen J, Laakso K, Koskenniemi K, Kankainen M, Savijoki K, Nyman TA, de Vos WM, Tynkkynen S, Kalkkinen N, Varmanen P. Effect of acid stress on protein expression and phosphorylation in Lactobacillus rhamnosus GG. J Proteomics. 2012;75(4):1357–74.

    Article  CAS  Google Scholar 

  22. Koskenniemi K, Laakso K, Koponen J, Kankainen M, Greco D, Auvinen P, Savijoki K, Nyman TA, Surakka A, Salusjärvi T, et al. Proteomics and transcriptomics characterization of bile stress response in probiotic Lactobacillus rhamnosus GG. Mol Cell Proteomics. 2011;10(2):M110002741.

    Article  Google Scholar 

  23. Prasad J, McJarrow P, Gopal P. Heat and osmotic stress responses of probiotic Lactobacillus rhamnosus HN001 (DR20) in relation to viability after drying. Appl Environ Microbiol. 2003;69(2):917–25.

    Article  CAS  Google Scholar 

  24. Carvalho SM, Farshchi Andisi V, Gradstedt H, Neef J, Kuipers OP, Neves AR, Bijlsma JJ. Pyruvate oxidase influences the sugar utilization pattern and capsule production in Streptococcus pneumoniae. PLoS ONE. 2013;8(7): e68277.

    Article  CAS  Google Scholar 

  25. Görke B, Stülke J. Carbon catabolite repression in bacteria: many ways to make the most out of nutrients. Nat Rev Microbiol. 2008;6(8):613–24.

    Article  Google Scholar 

  26. de Felipe F, Gaudu P. Multiple control of the acetate pathway in Lactococcus lactis under aeration by catabolite repression and metabolites. Appl Microbiol Biot. 2009;82(6):1115–22.

    Article  Google Scholar 

  27. Goffin P, Muscariello L, Lorquet F, Stukkens A, Prozzi D, Sacco M, Kleerebezem M, Hols P. Involvement of pyruvate oxidase activity and acetate production in the survival of Lactobacillus plantarum during the stationary phase of aerobic growth. Appl Environ Microbiol. 2006;72(12):7933–40.

    Article  CAS  Google Scholar 

  28. Castaldo C, Siciliano RA, Muscariello L, Marasco R, Sacco M. CcpA affects expression of the groESL and dnaK operons in Lactobacillus plantarum. Microb Cell Fact. 2006;5:35.

    Article  Google Scholar 

  29. Zhang C, Guo T, Xin Y, Gao X, Kong J. Catabolite responsive element deficiency of xyl operon resulting in carbon catabolite derepression in Lactobacillus fermentum 1001. J Appl Microbiol. 2016;120(1):126–37.

    Article  CAS  Google Scholar 

  30. Zhang C, Yang L, Ding Z, Yin B, Chen D, Guan C, Gu R. New selective media for isolation and enumeration of Lactobacillus rhamnosus and Streptococcus thermophilus. J Food Meas Charact. 2019;13(2):1431–9.

    Article  Google Scholar 

  31. Liu H, Li M, He X, Cai S, He X, Lu X. Transcriptome sequencing and characterization of ungerminated and germinated spores of Nosema bombycis. Acta Biochim Biophys Sin (Shanghai). 2016;48(3):246–56.

    Article  CAS  Google Scholar 

  32. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.

    Article  Google Scholar 

  33. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.

    Article  CAS  Google Scholar 

  34. Frazee AC, Pertea G, Jaffe AE, Langmead B, Salzberg SL, Leek JT. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat Biotechnol. 2015;33(3):243–6.

    Article  CAS  Google Scholar 

  35. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  Google Scholar 

  36. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(Database):D480-484.

    Article  CAS  Google Scholar 

  37. CNCB-NGDC Members and Partners. Database resources of the national genomics data center china national center for bioinformation in 2022. Nucleic Acids Res. 2022;50(D1):D27-38.

    Article  Google Scholar 

  38. Chen T, Chen X, Zhang S, Zhu J, Tang B, Wang A, Dong L, Zhang Z, Yu C, Sun Y, et al. The genome sequence archive family: toward explosive data growth and diverse data types. Genom Proteom Bioinf. 2021;19(4):578–83.

    Article  Google Scholar 

Download references

Acknowledgements

Sequencing service were provided by Shanghai Majorbio Bio-pharm Technology Co., Ltd., Shanghai, China.

Funding

This work was supported by the National Natural Science Foundation of China (No. 32272333, No. 31972094), the Scientific and Technological Innovation Platform Co-built by Yangzhou City-Yangzhou University (YZ2020265) and the “Lv Yang Jin Feng” Project 2019 in Yangzhou.

Author information

Authors and Affiliations

Authors

Contributions

CCZ and RXG conceived and designed research. HHC and YMH conducted experiments, and analyzed the data. CCZ and HHC wrote the manuscript. YCW, DWC, GRG, and YJH reviewed the manuscript, performed investigation, and analyzed data. All authors read and approved the manuscript.

Corresponding author

Correspondence to Ruixia Gu.

Ethics declarations

Competing interests

All author declares that they have no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, C., Cheng, H., Han, Y. et al. Transcriptome-phenotype matching analysis of how nitrogen sources influence Lacticaseibacillus rhamnosus tolerance to heat stress and oxidative stress. Microb Cell Fact 21, 257 (2022). https://doi.org/10.1186/s12934-022-01985-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12934-022-01985-0

Keywords