Next Article in Journal
Role of Visual Assessment of High-Quality Cakes in Emotional Response of Consumers
Previous Article in Journal
A Targeted Metabolomics Approach to Study Secondary Metabolites and Antioxidant Activity in ‘Kinnow Mandarin’ during Advanced Fruit Maturity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Foods
Volume 11
Issue 10
10.3390/foods11101411
foods-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ribosome Profiling Reveals Genome-Wide Cellular Translational Regulation in Lacticaseibacillus rhamnosus ATCC 53103 under Acid Stress

by
Xuejing Fan
1,
Kenan Zhang
1,
Zongcai Zhang
1,
Zhen Zhang
1,
Xue Lin
1,
Xin Liu
1,
Zhen Feng
1,2,* and
Huaxi Yi
3,*
1
Key Laboratory of Dairy Science, Ministry of Education, College of Food Science, Northeast Agricultural University, 600 Changjiang Road, Harbin 150030, China
2
Spice and Beverage Research Institute, Chinese Academy of Tropical Agricultural Sciences, Wanning 571533, China
3
College of Food Science and Engineering, Ocean University of China, Qingdao 266100, China
*
Authors to whom correspondence should be addressed.
Foods 2022, 11(10), 1411; https://0-doi-org.brum.beds.ac.uk/10.3390/foods11101411
Submission received: 29 March 2022 / Revised: 10 May 2022 / Accepted: 11 May 2022 / Published: 13 May 2022
(This article belongs to the Section Food Microbiology)
Browse Figures
Versions Notes

Abstract

:
During fermentation and food processing, Lacticaseibacillus rhamnosus ATCC 53103 can encounter many adverse conditions, and acid stress is one of them. The purpose of the present study was to investigate the influence of acid stress on the global translational and transcriptional regulation of Lacticaseibacillus rhamnosus ATCC 53103. Two pH values (pH 6.0 vs. pH 5.0) were applied, the effects of which were studied via ribosome profiling and RNA sequencing assay. Under acid stress, many genes showed differential changes at the translational and transcriptional levels. A total of 10 genes showed different expression trends at the two levels. The expression of 337 genes—which mainly participated in the ABC transporters, amino acid metabolism, and ribosome functional group assembly pathways—was shown to be regulated only at the translational level. The translational efficiency of a few genes participating in the pyrimidine and amino acid metabolism pathways were upregulated. Ribosome occupancy data suggested that ribosomes accumulated remarkably in the elongation region of open reading frame regions under acid stress. This study provides new insights into Lacticaseibacillus rhamnosus ATCC 53103 gene expression under acid stress, and demonstrates that the bacterium can respond to acid stress with synergistic translational and transcriptional regulation mechanisms, improving the vitality of cells.

1. Introduction

In recent years, Lacticaseibacillus (formerly Lactobacillus) rhamnosus Gorbach Goldin (L. rhamnosus GG) has been widely used in fermented food and dairy products, such as yogurt, fermented coffee, juice, etc. [1,2,3,4,5,6]. During fermentation and food processing, the continuous accumulation of organic acid causes a decrease in environmental pH, which affects the synthesis of nucleic acids and proteins, the degradation of substrates, and the utilization of nutrients, and also inhibits the growth and metabolism of cells [7,8]. Many acid adaptation mechanisms have been developed in lactic acid bacteria (LAB), such as the proton-translocating ATPase, regulatory factors, and transcription factors, which have been used to maintain cytoplasmic pH homeostasis and improve the resistance ability of LAB under acid stress [9,10,11]. Previous studies on the effect of acid stress on LAB have solely investigated the regulation of gene expression at the transcriptional level [12,13]. At present, translational regulations have been widely demonstrated to be involved in cells’ responses to environmental stresses, and mRNA translation plays an important role in gene expression regulation, as has also has increasingly recognized by researchers [14]. Meanwhile, the way in which LAB regulates gene expression at the translational level in response to acid stress is not fully understood.
Ribosomal profiling (Ribo-Seq) addresses the need for global expression assessment by estimating the ribosome occupancy and single-codon resolution, and quantifying the translation to the genome-wide level. This can precisely delineate translated regions and reveal the full coding potential of the genome [15,16]. Ribosome-protected mRNA fragments (RPFs) and total mRNAs can be quantified simultaneously by sequencing. The distribution of RPF reads illustrates ribosome occupation, while their abundance reveals the ribosome density. They are also used to measure the translational efficiency (TE), which ultimately determines the protein expression level [17,18]. Even since the advent of ribosome profiling, this important and powerful technique has been widely used in many studies of cellular activities in various organisms under different environmental stresses, such as the decrease in the number of RPFs (21 nt) and the translational regulation of Saccharomyces cerevisiae in response to high osmotic pressure and oxidative stresses [19], or the adaptation of yeasts’ responses to oxidative and amino acid starvation stresses [14,20]. The ribosomal footprints gathered in the open reading frames (ORFs) could increase the viability of Escherichia coli under heat stress [21]. Studies have also been conducted on the translational regulation of retinal pigment epithelium cells’ responses to hypoxic stress [22], and the translational regulation of Eubacterium limosum under autotrophic growth conditions [23]. However, the study of translational regulation in LAB has scarcely been reported.
In this study, simultaneous Ribo-Seq and RNA-Seq were used to establish quantitative translational and transcriptional regulation of L. rhamnosus ATCC 53103 under acid stress. The results showed that about 40% of identified differentially expressed genes were shared by the translational and transcriptional levels, and the ribosomal occupancy of the AS group was higher than that of the CG group in the elongation region of ORFs, which affected the TE of L. rhamnosus ATCC 53103 under acid stress. These results provide preliminary insights and new strategies for increasing the viability of L. rhamnosus under acid stress.

2. Materials and Methods

2.1. Strain, Culture Conditions, and Acid Stress

L. rhamnosus ATCC 53103 was obtained from the American Type Culture Collection (Manassas, VA, USA). The strain was revitalized in de Man, Rogosa, and Sharpe (MRS, HuanKai Microbial Sci. and Tech. Co., Ltd., Guangdong, China) medium at 37 °C. For acid stress, the cells were collected via centrifugation when the OD600 reached 1.8 (6000× g, 10 min, 4 °C). Then, the strain with the inoculation of 2% was resuspended in MRS medium (pH 5.0). The OD600 and the cell number of the cultures were 0.212 and 6.09 log10CFU/mL at t0, respectively. The growth of the strain was measured with an automatic growth curve analyzer (Bioscreen C Pro; OY Growth Curves, Helsinki, Finland) at 600 nm (OD600).

2.2. Total RNA Extraction and Library Construction

Cultures were harvested by centrifugation when the OD600 reached 1.0 (6000× g for 15 min at 4 °C). Total RNA was extracted, and the library was constructed as described in a previous study [24].

2.3. Ribo-Seq

The RPFs were prepared and the Ribo-Seq library was constructed as described in a previous study [24]. We added chloramphenicol (200 µg/mL) to the cultures and shook them for 2 min when the OD600 reached 1.0. The cells were washed using 20 mL of resuspension buffer and collected via centrifugation. Mixing the cells with cell lysis buffer and incubating them on ice for 10 min, the cells were triturated and the lysate was collected. Then, 5 µL of DNase I (Roche) and 7.5 µL of RNase I (Abgene) were added to 300 µL of lysate to cultivate for 45 min, and 10 µL of RNase inhibitor (Ambion) was added to the lysate to stop the digestion. RPFs were isolated using the Concentrator-25 Kit (R1017, Zymo Research, Orange County, CA, USA) and RNAClean. The Ribo-Seq library was constructed using the NEBNext® Multiple Small RNA Library Prep Set for Illumina® (catalog no. E7300S, E7300L). The 140–160 bp PCR products were enriched to generate a cDNA library, and sequenced using Illumina HiSeqTM 2500.

2.4. Analysis of Differentially Expressed Genes

The gene expression level was normalized using the fragments per kilobase of transcript per million (FPKM) mapped reads method to eliminate the influence of different gene lengths and the amount of sequencing data on the calculation of gene expression. The edgeR package (http://www.rproject.org/ 28 December 2020) was used to ascertain differentially expressed genes (DEGs). Gene Ontology (GO) annotation and the Kyoto Encyclopedia for Genes and Genomes (KEGG) pathways of DEGs were analyzed. DEGs were divided into five groups as described in a previous study [24]: unchanged, homodirectional, opposite, translational, and transcriptional.

2.5. Analysis of DEGs’ Translational Efficiency, and Correlation with DEGs’ Expression

The analysis of the TE of all DEGs was carried out using RiboDiff [24,25]. According to the TE of the DEGs, and their expression at the transcriptional level, all DEGs were divided into five groups: homodirectional, opposite, TE, transcriptional, and unchanged.

2.6. Data Analysis

Raw data were filtered as described in a previous study [26]. The trimmed reads were mapped to the reference transcriptomic data of L. rhamnosus (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/nuccore/NC_013198.1 28 December 2020). Retained reads were aligned, and DEGs were identified and calculated as described in a previous study [24]; we chose genes with a false discovery rate (FDR) < 0.05 and ∣log2(fold change)∣ > 1, which were defined as DEGs.

3. Results

3.1. Effects of Acid Stress on DEGs’ Expression and Function

The regulation of translation and transcription under acid stress was monitored by sequencing of RPFs and total mRNAs in the control group (CG, pH 6.0) and the acid stress group (AS, pH 5.0). Gene expression changes induced by acid stress at the translational and transcriptional levels are shown in Figure 1. There were 132 and 372 upregulated DEGs at the transcriptional and translational levels, respectively, and there were 104 and 117 downregulated DEGs at these two levels, respectively (Figure 1a–c). Around 40% of the regulated genes were shared between the translational and transcriptional levels, including 12.2% of upregulated genes and 27.0% of downregulated genes (Figure 1d,e). GO annotation and KEGG pathway analysis of DEGs at these two levels were performed (Figure 2). The results indicated that DEGs largely overlapped in GO at the two levels. Most DEGs were prominently upregulated in GO at the translational level, and were mainly involved in metabolic and single-organism processes, catalytic activity, and the membrane components (Figure 2a). Many KEGG pathways appeared to largely overlap, such as ABC transporters, quorum sensing, two-component systems, biosynthesis of amino acids, and metabolism of purine pathways. However, large numbers of KEGG pathways did not overlap at these two levels, such as aminoacyl-tRNA biosynthesis and RNA degradation pathways (Figure 2b).

3.2. Effects of Acid Stress on the Dynamic Profiles of Translation and Transcription

To access the overall changes in DEGs’ expression trends at the two levels simultaneously, DEGs were classified into distinct groups (Figure 3a).
There were 102 co-regulated DEGs at these two levels, and their expression had similar trends (yellow dots); 55 of these DEGs were co-upregulated, and 47 were co-downregulated. These DEGs mainly participated in the purine metabolism, ABC transporters, two-component system, pyrimidine metabolism, and biosynthesis of amino acids pathways. The DEGs purQ, purS, purM, purN, purH, and purD were downregulated by 193.03-, 151.33-, 57.73-, 47.84-, 16.14-, and 6.32-fold at the transcriptional level, respectively; they also were downregulated from 50- to 200-fold at the translational level. They mainly participated in the purine metabolism pathway. The DEGs dltB, dltD, bceA, and agrA were also downregulated by about 2- to 3-fold at the two levels, respectively; they mainly participated in the two-component system pathway.
There were 10 co-regulated DEGs at these two levels, and their expression had opposite trends (blue dots). ArgG was upregulated by 2.11-fold at the translational level and downregulated by 2.14-fold at the transcriptional level; it participated in the amino acid biosynthesis pathway. CitC was upregulated by 3.18-fold at the translational level and downregulated by 2.33-fold at the transcriptional level; it participated in the two-component system pathway. PyrR1 was upregulated by 10.95-fold at the translational level and downregulated by 2.89-fold at the transcriptional level; it participated in the pyrimidine metabolism pathway.
There were 124 regulated DEGs at the transcriptional level only (red dots). These DEGs participated in the fatty acid metabolism, microbial metabolism in diverse environments, biosynthesis of secondary metabolites, and pyruvate metabolism pathways. The DEGs accA, accC1, fabZ, accB, fabF, bkr2, and pksA were each downregulated by about 2-fold. They mainly participated in the fatty acid metabolism pathway.
There were 337 regulated DEGs at the translation level only (green dots). The DEGs dppE, livH, braF, tcyJ, ftsX, potD, nlpA, dppE, gbuA, and oppA were upregulated, and they mainly participated in the ABC transporter pathway. The DEGs trpF, trpC, and trpD were upregulated by 14.51-, 12.16-, and 107.45- fold, respectively, and they mainly participated in the tyrosine, phenylalanine, and tryptophan biosynthesis pathways. The DEGs carB, carA, pyrB, gltB, argH, and mtnU were upregulated by 3.32-, 4.48, 11.59-, 4.11-, 4.53-, and 3.02-fold, respectively, and they participated in the alanine, aspartate, and glutamate metabolism pathways. The DEGs rpsP, rpmI, rplQ, rplR, rpsH, and rpmH were upregulated, and they mainly participated in the assembly of ribosome functional groups. Many DEGs participated in the phosphotransferase system (PTS), fructose, mannose, pyrimidine, starch, sucrose, amino sugar, nucleotide sugar metabolism, and microbial metabolism in diverse environments pathways; these DEGs also were upregulated. Meanwhile, the DEGs coaD, ppdK, and ldh participated in metabolic pathways, and they were each downregulated by about 2-fold. mprF, encoding lysyl-tRNA synthetase, was also downregulated by only 2.32-fold.

3.3. Effects of Acid Stress on DEGs’ Translational Efficiency

With FDR < 0.05 and∣log2FC∣ > 1, the TE of 27 DEGs was upregulated, while the TE of 5 DEGs was downregulated (Figure 3b and Table 1). The KEGG pathway enrichment analysis determined that these DEGs primarily participated in the amino acid, pyrimidine, and purine metabolism pathways, as well as the biosynthesis of secondary metabolites pathway. The DEGs trpF and trpD participated in the tyrosine, phenylalanine, and tryptophan biosynthesis pathways, while the DEGs carA, pyrB, and argH participated in the alanine, aspartate, and glutamate metabolism pathways. The DEGs pyrC and pyrR1 participated in the pyrimidine metabolism pathway. The DEG acpP participated in the biosynthesis of secondary metabolites pathway. The TE of these DEGs was upregulated. Meanwhile, the DEGs purD and purH participated in the biosynthesis of purine metabolism pathway, and their TE was downregulated (Table 2).

3.4. Effects of Acid Stress on the Dynamic Profiles of Translational Efficiency and Transcription

Based on the expression of genes at the transcriptional level and the regulation of TE, the DEGs were divided into five groups (Figure 3c). There were 21 DEGs that were not regulated at the transcriptional level; however, their TE was changed significantly. The DEGs carA, pyrB, and argH participated in the alanine, aspartate, and glutamate metabolism pathways, and their TE was upregulated. The DEGs trpF and trpD participated in the tyrosine, phenylalanine, and tryptophan biosynthesis pathways, and their TE was also upregulated. The TE of pyrC and accp was upregulated under acid stress; these DEGs mainly participated in the metabolism of pyrimidine and the biosynthesis of secondary metabolites, respectively. The DEGs hsp18 and HSP17.6C encode heat-shock proteins, and their TE was upregulated. The TE of fmnP was also upregulated; this gene encodes predicted membrane proteins.
There were 225 DEGs that were only regulated at the transcriptional level; however, their TE was not changed significantly. The DEGs accA, accC1, accD, fabZ, accB, fabF, bkr2, and pksA participated in the biosynthesis and metabolism of fatty acids pathway, and their expression was downregulated significantly. The DEGs purA, guaC, purK, purB, purN, purM, purF, purL, purQ, purS, purC, purK, and purE mainly participated in the purine metabolism pathway, and their expression was downregulated significantly. The DEGs dltB, dltC, dltD, citC, bceA, and agrA mainly participated in the two-component system pathway; their expression was also downregulated. The DEGs tauA, potA, potB, pstS1, ssuB1, macB2, nodI, ugpB, ugpA, and ugpC mainly participated in the ABC transporter pathway, and their expression was upregulated. In addition, 22 DEGs were identified in the biosynthesis of secondary metabolites pathway, and they were all regulated only at the transcriptional level.
Only 4 DEGs’ expression had regulation trends similar to their TE. These included the DEGs purD and purH, involved in the purine metabolism and biosynthesis of secondary metabolites pathways, which were downregulated at the transcriptional level, and their TE was also downregulated. A total of 7 DEGs’ expression had opposite regulation trends to their TE. The DEGs glnP and nodI participated in the ABC transporter pathway, and their expression was downregulated at the transcriptional level; however their TE was upregulated. The expression of pyrR1 was also downregulated at the transcriptional level, while its TE was upregulated; this gene was involved in the pyrimidine metabolism pathway.

3.5. Effects of Acid Stress on Ribosome Accumulation in the ORF Regions

In order to quantify the translational regulation, the ribosomes’ average occupancy was computed around the start and stop codons. The distribution of averaged reads in the 5’ and 3’UTRs is shown in Figure 3d,e. Two groups of footprints mapped to the initiation region showed a clear and significant increase, with the 5’ end roughly at the nucleotide −15 position. The averaged read occupancy of the CG was higher than that of the AS group. The results showed that acid stress did not reduce the averaged read occupancy of the ribosomes, which increased from the region of nucleotide 20 to 120. The average ribosome occupancy of the AS group at nucleotide 120 after the start codon was about double the reads of the CG. In addition, the average ribosome occupancy of the AS group at the nucleotides from −100 to −20 in the stop codon regions was about 50% higher than that of the CG.

4. Discussion

Not only were the changes of gene expression explored on a genome-wide scale, the interactions between the translational and transcriptional regulation were also revealed in L. rhamnosus ATCC 53103 under acid stress using Ribo-Seq and RNA-Seq. There were 488 regulated DEGs at the translational level, while the number of regulated DEGs at the transcriptional level was only 236. These results indicate that translational regulation can be a more rapid and direct means of acid stress response compared to transcriptional regulation, because no new mRNA needs to be produced [27]. Hence, translational regulation plays an integral role in acid stress responses for L. rhamnosus ATCC 53103. There were 337 regulated DEGs at the translational level only, which mainly participated in the ABC transporter, amino acids biosynthesis and metabolism, and ribosome assembly pathways. The translational regulation of ABC transporter genes could regulate the synthesis of ABC transporters, which mediate the active transport of molecules across the cell membranes, and this process is helpful for L. rhamnosus ATCC 53103 to resist acid shock rapidly and maintain intercellular homeostasis [28,29]. Amino-acid-dependent acid-tolerance systems have been approved as a sophisticated defense mechanism against damage caused by acid in LAB [30]. The upregulated expression of DEGs participating in the biosynthesis and metabolism of amino acids at the translational level may help to increase the production of ATP and sequester protons, which are important for pH homeostasis [31]. Based on the high-level synthesis of proteins and the requirements of growing cells, the formation of ribosomes in the bacteria is an important and efficiently regulated process. The assembly of an intact ribosome requires many complex and regulated functions, such as the coordinated synthesis of ribosomal proteins [32]. In this study, some DEGs involved in the assembly of ribosome functional groups were upregulated only at the translational level under acid stress. This might help the bacteria to adjust the number of ribosomes in proportion to the metabolic state and the growth rate of cells in order to meet the demand for protein synthesis under acid stress [32]. Hence, translational regulation appeared to play an important role in L. rhamnosus ATCC 53103 under acid stress.
About 40% of DEGs were shared between the translational and transcriptional levels. There were 112 co-regulated DEGs at the two levels, and they mainly participated in the purine metabolism, ABC transporter, two-component system, pyrimidine metabolism, and amino acid biosynthesis pathways. The DEGs purH, purQ, purM, purN, purD, and purS were downregulated at both levels, and they were mainly involved in purine metabolism. As the signal mediator, the downregulation of pur genes could help bacteria to save energy and accelerate the growth of cells under acid stress, indicating that purine metabolism involves the bacterial stringent response, and could conceivably enhance resistance to acid stress in L. rhamnosus ATCC 53103 [33]. Normally, the two-component system can help bacteria to respond to environmental stresses, because it contains a target gene expression regulator and a membrane-associated histidine kinase [21]. The overlap between transcriptional and translational regulation could increase the flexibility of gene expression, which could improve bacterial adaptability to acid stress [34].
TE is an important index that can be used to identify previously unknown translational events [35]. Only 21 DEGs were not regulated at the transcriptional level, and their TE was changed significantly under acid stress. These DEGs participated in the amino acid metabolism pathway, and their TE increased significantly. Amino acids participate in the physiological processes of L. rhamnosus ATCC 53103, such as intracellular pH regulation, stress resistance, protein synthesis, and metabolic energy generation [31]. Although the genes’ expression was not changed significantly at the transcriptional level, the upregulation of TE still improved the biosynthesis of amino acids, and further increased the tolerance under acid stress. Many DEGs participated in the purine metabolism, fatty acid biosynthesis and metabolism, and two-component system pathways; their expression was downregulated at the transcriptional level. These results are different from previous results that showed that these genes’ expression was upregulated under different environmental stresses [30], and the TE of these DEGs was not changed significantly, indicating that translational regulation might control the production of effector proteins in response to certain stresses [21]. Therefore, we hypothesized that translational regulation might play different roles in cells under different stresses, and that the TE remaining relatively constant would be helpful for L. rhamnosus ATCC 53103 to respond to acid shock, while the downregulated expression of related DEGs might help to reduce cells’ energy consumption [36].
ORFs are important mediators of transcript-specific translational control—especially from the upstream ORFs to their downstream coding sequences [22]. In this study, the averaged read occupancy of the CG was higher than that of the AS group around the start codons. This might have been due to the reduction in ribosome abundance at the initiation region, which could partially rescue protein production under acid stress [37]. Meanwhile, the average ribosome occupancy of the AS group was higher than that of the CG in the elongation region. Typically, upstream ORFs would reduce the TE of the main downstream ORFs under normal conditions [38]. Another study also reported that when the translational inhibition by an upstream ORF was reduced, the ribosome occupancy and TE of numerous mRNAs were increased in Arabidopsis thaliana under hypoxic stress [39]. Stress-induced repression of ribosomes and the reduction in the initiation rates of genes have been thought to be adaptive, because they not only increase the pool of free ribosomes and tRNAs, but also minimize resource waste. These responses could have direct benefits in increasing the elongation rate of genes, rescuing protein production, and promoting cell growth under stress [37].

5. Conclusions

The present work revealed the landscape of translational and transcriptional regulation in L. rhamnosus ATCC 53103 under acid stress by using Ribo-Seq and RNA-Seq. At the translational and transcriptional levels, about 40% of DEGs were co-regulated, indicating that the overlap between the translational and transcriptional regulations improved the adaptability of L. rhamnosus ATCC 53103 under acid stress. Many DEGs were regulated only at the transcriptional level; however, their TE was not changed significantly. For some DEGs, L. rhamnosus ATCC 53103 appeared to possess specialized mechanisms to offset or enhance the effects of transcriptional regulation under acid stress. The ribosomal footprints of the AS group were higher than those of the CG in the elongation region of the ORFs, which also maintained the TE in L. rhamnosus ATCC 53103 under acid stress. These results indicate that ORFs play an important role in the regulation of gene expression in response to acid stress; other mechanisms of the effects of ORFs still require further study.

Author Contributions

Methodology, X.F., K.Z. and Z.Z. (Zongcai Zhang); software, Z.Z. (Zhen Zhang) and X.L. (Xue Lin); data curation, X.F. and X.L. (Xin Liu); writing-original draft preparation, X.F.; writing-review and editing, Z.F. and H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by grants from the National Natural Science Foundation of China, grant number 31771989.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zheng, J.; Wittouck, S.; Salvetti, E.; Franz, C.M.A.P.; Harris, H.M.B.; Mattarelli, P.; O’Toole, P.W.; Pot, B.; Vandamme, P.; Walter, J.; et al. A taxonomic note on the genus Lactobacillus: Description of 23 novel genera, emended description of the genus Lactobacillus Beijerinck 1901, and union of Lactobacillaceae and Leuconostocaceae. Int. J. Syst. Evol. Micr. 2020, 70, 2782–2858. [Google Scholar] [CrossRef] [PubMed]
  2. Hekmat, S.; Soltani, H.; Reid, G. Growth and survival of Lactobacillus reuteri RC-14 and Lactobacillus rhamnosus GR-1 in yogurt for use as a functional food. Innov. Food Sci. Emerg. Technol. 2009, 10, 293–296. [Google Scholar] [CrossRef]
  3. Chan, M.Z.A.; Toh, M.; Liu, S.Q. Growth, survival, and metabolic activities of probiotics Lactobacillus rhamnosus GG and Saccharomyces cerevisiae var. boulardii CNCM-I745 in fermented coffee brews. Int. J. Food Microbiol. 2021, 350, 109229. [Google Scholar] [CrossRef] [PubMed]
  4. Lei, W.; Luo, J.; Wu, K.; Chen, Q.; Hao, L.; Zhou, X.; Wang, X.; Liu, C.; Zhou, H. Dendrobium candidum extract on the bioactive and fermentation properties of Lactobacillus rhamnosus GG in fermented milk. Food Biosci. 2021, 41, 100987. [Google Scholar] [CrossRef]
  5. Wang, Y.; Huo, K.; Gao, L.; Cai, D.; Wang, B.; Zhao, G.; Liu, J.; Hao, J. Improvement of l-lactic acid fermentation by Lactobacillus rhamnosus immobilized in polyvinyl alcohol/Fe3O4 composite using sweet sorghum juice. Ind. Crop. Prod. 2022, 182, 114922. [Google Scholar] [CrossRef]
  6. Sengun, I.Y.; Kirmizigul, A.; Atlama, K.; Yilmaz, B. The viability of Lactobacillus rhamnosus in orange juice fortified with nettle (Urtica dioica L.) and bioactive properties of the juice during storage. LWT-Food Sci. Technol. 2020, 118, 108707. [Google Scholar] [CrossRef]
  7. Zhu, Z.M.; Ji, X.M.; Wu, Z.M.; Zhang, J.; Du, G.C. Improved acid-stress tolerance of Lactococcus lactis NZ9000 and Escherichia coli bl21 by overexpression of the anti-acid component rect. J. Ind. Microbiol. Biotechnol. 2018, 45, 1091–1101. [Google Scholar] [CrossRef]
  8. Yanez, R.; Marques, S.; Girio, F.M.; Roseiro, J.C. The effect of acid stress on lactate production and growth kinetics in Lactobacillus rhamnosus cultures. Process. Biochem. 2008, 43, 356–361. [Google Scholar] [CrossRef]
  9. Wu, C.D.; Huang, J.; Zhou, R.Q. Progress in engineering acid stress resistance of lactic acid bacteria. Appl. Microbiol. Biotechnol. 2014, 98, 1055–1063. [Google Scholar] [CrossRef]
  10. Koponen, J.; Laakso, K.; Koskenniemi, K.; Kankainen, M.; Savijoki, K.; Nyman, T.A.; de Vos, W.M.; Tynkkynen, S.; Kalkkinen, N.; Varmanen, P. Effect of acid stress on protein expression and phosphorylation in Lactobacillus rhamnosus GG. J. Proteom. 2012, 75, 1357–1374. [Google Scholar] [CrossRef]
  11. Zhu, Z.M.; Yang, J.H.; Yang, P.S.; Wu, Z.M.; Zhang, J.; Du, G.C. Enhanced acid-stress tolerance in lactococcus lactis nz9000 by overexpression of abc transporters. Microb. Cell Factories 2019, 18, 136. [Google Scholar] [CrossRef] [PubMed]
  12. Broadbent, J.R.; Larsen, R.L.; Deibel, V.; Steele, J.L. Physiological and transcriptional response of Lactobacillus casei ATCC 334 to acid stress. J. Bacteriol. 2010, 192, 2445–2458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wei, Y.X.; Gao, J.; Liu, D.B.; Li, Y.; Liu, W.L. Adaptational changes in physiological and transcriptional responses of Bifidobacterium longum involved in acid stress resistance after successive batch cultures. Microb. Cell. Fact. 2019, 18, 156. [Google Scholar] [CrossRef] [PubMed]
  14. Li, G.W.; Oh, E.; Weissman, J.S. The anti-shine-dalgarno sequence drives translational pausing and codon choice in bacteria. Nature 2012, 484, 538–541. [Google Scholar] [CrossRef] [PubMed]
  15. Ingolia, N.T.; Ghaemmaghami, S.; Newman, J.R.S.; Weissman, J.S. Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. Science 2009, 324, 218–223. [Google Scholar] [CrossRef] [Green Version]
  16. Ingolia, N.T. Ribosome footprint profiling of translation throughout the genome. Cell 2016, 165, 22–33. [Google Scholar] [CrossRef] [Green Version]
  17. Ingolia, N.T.; Brar, G.A.; Rouskin, S.; McGeachy, A.M.; Weissman, J.S. The ribosome profiling strategy for monitoring translation in vivo by deep sequencing of ribosome-protected mrna fragments. Nat. Protoc. 2012, 7, 1534–1550. [Google Scholar] [CrossRef]
  18. Waudby, C.A.; Dobson, C.M.; Christodoulou, J. Nature and regulation of protein folding on the ribosome. Trends Biochem. Sci. 2019, 44, 914–926. [Google Scholar] [CrossRef] [Green Version]
  19. Wu, C.C.C.; Zinshteyn, B.; Wehner, K.A.; Green, R. High-resolution ribosome profiling defines discrete ribosome elongation states and translational regulation during cellular stress. Mol. Cell 2019, 73, 1–12. [Google Scholar] [CrossRef] [Green Version]
  20. Gerashchenko, M.V.; Lobanov, A.V.; Gladyshev, V.N. Genome-wide ribosome profiling reveals complex translational regulation in response to oxidative stress. Proc. Natl. Acad. Sci. USA 2012, 109, 17394–17399. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, Y.Q.; Xiao, Z.T.; Zou, Q.; Fang, J.H.; Wang, Q.F.; Yang, X.R.; Gao, N. Ribosome profiling reveals genome-wide cellular translational regulation upon heat stress in Escherichia coli. Genom. Proteom. Bioinf. 2017, 15, 324–330. [Google Scholar] [CrossRef] [PubMed]
  22. Jiang, Z.W.; Yang, J.Q.; Dai, A.M.; Wang, Y.M.; Li, W.; Xie, Z. Ribosome profiling reveals translational regulation of mammalian cells in response to hypoxic stress. BMC Genom. 2017, 18, 638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Song, Y.; Shin, J.; Jin, S.; Lee, J.K.; Kim, D.R.; Kim, S.C.; Cho, S.; Cho, B.K. Genome-scale analysis of syngas fermenting acetogenic bacteria reveals the translational regulation for its autotrophic growth. BMC Genom. 2018, 19, 837. [Google Scholar] [CrossRef] [PubMed]
  24. Fan, X.; Bao, T.; Yi, H.; Zhang, Z.; Zhang, K.; Liu, X.; Lin, X.; Zhang, Z.; Feng, Z. Ribosome profiling and RNA sequencing reveal genome-wide cellular translation and transcription regulation under osmotic stress in Lactobacillus rhamnosus ATCC 53103. Front. Microbiol. 2021, 12, 781454. [Google Scholar] [CrossRef] [PubMed]
  25. Zhong, Y.; Karaletsos, T.; Drewe, P.; Sreedharan, V.T.; Kuo, D.; Singh, K.; Wendel, H.G.; Ratsch, G. Ribodiff: Detecting changes of mRNA translation efficiency from ribosome footprints. Bioinformatics 2017, 33, 139–141. [Google Scholar] [CrossRef]
  26. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  27. Sonenberg, N.; Hinnebusch, A.G. Regulation of translation initiation in eukaryotes: Mechanisms and biological targets. Cell 2009, 136, 731–745. [Google Scholar] [CrossRef] [Green Version]
  28. Shi, M.Y.; Wang, S.S.; Zhang, Y.; Wang, S.; Zhao, J.; Feng, H.; Sun, P.P.; Fang, C.B.; Xie, X.B. Genome-wide characterization and expression analysis of ATP-binding cassette (abc) transporters in strawberry reveal the role of fvabcc11 in cadmium tolerance. Sci. Hortic. 2020, 271, 109464. [Google Scholar] [CrossRef]
  29. Qin, J.Y.; Wang, X.W.; Wang, L.D.; Zhu, B.B.; Zhang, X.H.; Yao, Q.S.; Xu, P. Comparative transcriptome analysis reveals different molecular mechanisms of bacillus coagulans 2–6 response to sodium lactate and calcium lactate during lactic acid production. PLoS ONE 2015, 10, e0124316. [Google Scholar] [CrossRef] [Green Version]
  30. Guan, N.Z.; Liu, L. Microbial response to acid stress: Mechanisms and applications. Appl. Microbiol. Biotechnol. 2020, 104, 51–65. [Google Scholar] [CrossRef] [Green Version]
  31. Wu, H.; Liu, J.G.; Miao, S.; Zhao, Y.; Zhu, H.J.; Qiao, M.Q.; Saris, P.E.J.; Qiao, J.J. Contribution of YthA, a PspC family transcriptional regulator of Lactococcus lactis F44 acid tolerance and nisin yield: A transcriptomic approach. Appl. Environ. Microb. 2018, 84, e02483-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Piersimoni, L.; Giangrossi, M.; Marchi, P.; Brandi, A.; Gualerzi, C.O.; Pon, C.L. De novo synthesis and assembly of rRNA into ribosomal subunits during cold acclimation in Escherichia coli. J. Mol. Biol. 2016, 428, 1558–1573. [Google Scholar] [CrossRef] [PubMed]
  33. Goncheva, M.I.; Flannagan, R.S.; Sterling, B.E.; Laakso, H.A.; Friedrich, N.C.; Kaiser, J.C.; Watson, D.W.; Wilson, C.H.; Sheldon, J.R.; McGavin, M.J.; et al. Stress-induced inactivation of the Staphylococcus aureus purine biosynthesis repressor leads to hypervirulence. Nat. Commun. 2019, 10, 775. [Google Scholar] [CrossRef] [PubMed]
  34. Lei, L.; Shi, J.P.; Chen, J.; Zhang, M.; Sun, S.L.; Xie, S.J.; Li, X.J.; Zeng, B.A.; Peng, L.Z.; Hauck, A.; et al. Ribosome profiling reveals dynamic translational landscape in maize seedlings under drought stress. Plant J. 2015, 84, 1206–1218. [Google Scholar] [CrossRef]
  35. Ingolia, N.T. Ribosome profiling: New views of translation, from single codons to genome scale applications of next-generation sequencing-innovation. Nat. Rev. Genet. 2014, 15, 205–213. [Google Scholar] [CrossRef] [PubMed]
  36. Zhou, A.L.; Cao, Y.F.; Zhou, D.G.; Hu, S.F.; Tan, W.J.; Xiao, X.L.; Yu, Y.G.; Li, X.F. Global transcriptomic analysis of Cronobacter Sakazakii CICC 21544 by RNA-seq under inorganic acid and organic acid stresses. Food Res. Int. 2020, 130, 108963. [Google Scholar] [CrossRef] [PubMed]
  37. Shah, P.; Ding, Y.; Niemczyk, M.; Kudla, G.; Plotkin, J.B. Rate-limiting steps in yeast protein translation. Cell 2013, 153, 1589–1601. [Google Scholar] [CrossRef] [Green Version]
  38. Hsu, M.K.; Chen, F.C. Selective constraint on the upstream open reading frames that overlap with coding sequences in animals. PLoS ONE 2012, 7, e48413. [Google Scholar] [CrossRef] [Green Version]
  39. Juntawong, P.; Girke, T.; Bazin, J.; Bailey-Serres, J. Translational dynamics revealed by genome-wide profiling of ribosome footprints in Arabidopsis. Proc. Natl. Acad. Sci. USA 2014, 111, E203–E212. [Google Scholar] [CrossRef] [Green Version]
Foods 11 01411 g001 550
Figure 1. (a) The volcano map of DEGs in L. rhamnosus ATCC 53103 at the transcriptional level under acid stress (∣log2fold change∣ > 1 and FDR < 0.05). (b) The volcano map of DEGs in L. rhamnosus ATCC 53103 at the translational level under acid stress (∣log2fold change∣ > 1 and FDR < 0.05). (c) The number of DEGs in L. rhamnosus ATCC 53103 under acid stress at the transcriptional or translational levels. (d) The relationships between upregulated DEGs at these two levels. (e) The relationships between downregulated DEGs at these two levels.
Figure 1. (a) The volcano map of DEGs in L. rhamnosus ATCC 53103 at the transcriptional level under acid stress (∣log2fold change∣ > 1 and FDR < 0.05). (b) The volcano map of DEGs in L. rhamnosus ATCC 53103 at the translational level under acid stress (∣log2fold change∣ > 1 and FDR < 0.05). (c) The number of DEGs in L. rhamnosus ATCC 53103 under acid stress at the transcriptional or translational levels. (d) The relationships between upregulated DEGs at these two levels. (e) The relationships between downregulated DEGs at these two levels.
Foods 11 01411 g001
Foods 11 01411 g002 550
Figure 2. (a) GO analysis of DEGs in the AS and CG groups at the translational and transcriptional levels. (b) Significantly enriched KEGG pathways of DEGs in the AS and CG groups at the translational and transcriptional levels.
Figure 2. (a) GO analysis of DEGs in the AS and CG groups at the translational and transcriptional levels. (b) Significantly enriched KEGG pathways of DEGs in the AS and CG groups at the translational and transcriptional levels.
Foods 11 01411 g002
Foods 11 01411 g003 550
Figure 3. (a) The gene expression changes at the translational and transcriptional levels simultaneously in L. rhamnosus ATCC 53103 under acid stress. (b) The volcano plot of the TE responses of DEGs under acid stress. (c) Simultaneous monitoring of the distribution of log2fold change TE and corresponding log2fold change mRNA abundance in L. rhamnosus ATCC 53103 under acid stress. The four dashed lines indicate log2fold change values of −1 or +1. (d) The density of ribosomal footprints around the start codon under acid stress. (e) The density of ribosomal footprints around the stop codon under acid stress.
Figure 3. (a) The gene expression changes at the translational and transcriptional levels simultaneously in L. rhamnosus ATCC 53103 under acid stress. (b) The volcano plot of the TE responses of DEGs under acid stress. (c) Simultaneous monitoring of the distribution of log2fold change TE and corresponding log2fold change mRNA abundance in L. rhamnosus ATCC 53103 under acid stress. The four dashed lines indicate log2fold change values of −1 or +1. (d) The density of ribosomal footprints around the start codon under acid stress. (e) The density of ribosomal footprints around the stop codon under acid stress.
Foods 11 01411 g003
Table
Table 1. Genes with differential translational efficiency regulation under acid stress.
Table 1. Genes with differential translational efficiency regulation under acid stress.
CategoryGene_IDGene_NameTE_CGTE_ASlog2FCFDR
UpregulatedLGG_RS00110--6.7543.022.670.0396827
LGG_RS00155--0.674.932.880.0153195
LGG_RS00490trpF0.2813.155.550.0303543
LGG_RS00500trpD0.6141.486.092.95 × 10−6
LGG_RS00715TTE16500.010.072.810.0409188
LGG_RS02170--0.011.77.411.666 × 10−9
LGG_RS02790TTE16500.010.072.810.0409188
LGG_RS03170HSP17.6C0.010.355.132.827 × 10−8
LGG_RS04550TTE16500.010.072.810.0409188
LGG_RS04690--0.10.792.980.0112087
LGG_RS04695nodI0.090.672.90.0327759
LGG_RS06620fmnP0.241.472.610.0184817
LGG_RS06985carA4.4325.522.530.0251229
LGG_RS06990pyrC1.7612.292.80.0081897
LGG_RS06995pyrB0.4810.294.421.185 × 10−7
LGG_RS07000pyrP0.3811.494.922.265 × 10−9
LGG_RS07005pyrR10.154.684.962.265 × 10−9
LGG_RS08210TTE16500.010.072.810.0409188
LGG_RS08915--0.864.822.490.0409188
LGG_RS09795clpC1.1519.984.123.804 × 10−6
LGG_RS10200acpP5.5633.812.60.0251229
LGG_RS10665--1.9410.392.420.0396827
LGG_RS13400Cndp11.9215.493.010.0219037
LGG_RS13405oppA1.7273.855.422.659E-06
LGG_RS13420hsp180.789.473.60.0001423
LGG_RS13465glnP0.020.122.580.0327759
LGG_RS13475argH0.241.62.740.0278452
DownregulatedLGG_RS05995--9032.64−8.420.0251229
LGG_RS08685purD15.220.44−5.115.176 × 10−10
LGG_RS08690purH22.291.58−3.825.257 × 10−6
LGG_RS09365oppA3.010.56−2.430.0219037
LGG_RS13815--40.174.72−3.090.002218
Note: FC, fold change; TE, translational efficiency.
Table
Table 2. DEG enrichment results of TE with KEGG pathway analysis.
Table 2. DEG enrichment results of TE with KEGG pathway analysis.
CategoryPathway IDPathway Namep-ValueGenesGenes
Upregulatedko00250Alanine, aspartate, and glutamate metabolism0.002LGG_RS06985, LGG_RS06995, LGG_RS13475carA, pyrB, argH
ko00240Pyrimidine metabolism0.002LGG_RS06985, LGG_RS06990, LGG_RS06995, LGG_RS07005carA, pyrC, pyrB, pyrR1
ko00400Phenylalanine, tyrosine, and tryptophan biosynthesis0.008LGG_RS00490, LGG_RS00500trpF, trpD
ko00999Biosynthesis of secondary metabolites—unclassified0.029LGG_RS10200acpP
ko01110Biosynthesis of secondary metabolites0.057LGG_RS00490, LGG_RS00500, LGG_RS13475trpF, trpD, argH
ko01100Metabolic pathways0.068LGG_RS00490, LGG_RS00500, LGG_RS06985, LGG_RS06990, LGG_RS06995, LGG_RS07005, LGG_RS13475trpF, trpD, carA, pyrC, pyrB, pyrR1, argH
ko00220Arginine biosynthesis0.071LGG_RS13475argH
ko01230Biosynthesis of amino acids0.108LGG_RS00490, LGG_RS00500, LGG_RS13475trpF, trpD, argH
Downregulatedko02024Quorum sensing0.151LGG_RS09365oppA
ko00230Purine metabolism0.187LGG_RS08685 LGG_RS08690purD, purH
ko02010ABC transporters0.432LGG_RS09365oppA
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fan, X.; Zhang, K.; Zhang, Z.; Zhang, Z.; Lin, X.; Liu, X.; Feng, Z.; Yi, H. Ribosome Profiling Reveals Genome-Wide Cellular Translational Regulation in Lacticaseibacillus rhamnosus ATCC 53103 under Acid Stress. Foods 2022, 11, 1411. https://0-doi-org.brum.beds.ac.uk/10.3390/foods11101411

AMA Style

Fan X, Zhang K, Zhang Z, Zhang Z, Lin X, Liu X, Feng Z, Yi H. Ribosome Profiling Reveals Genome-Wide Cellular Translational Regulation in Lacticaseibacillus rhamnosus ATCC 53103 under Acid Stress. Foods. 2022; 11(10):1411. https://0-doi-org.brum.beds.ac.uk/10.3390/foods11101411

Chicago/Turabian Style

Fan, Xuejing, Kenan Zhang, Zongcai Zhang, Zhen Zhang, Xue Lin, Xin Liu, Zhen Feng, and Huaxi Yi. 2022. "Ribosome Profiling Reveals Genome-Wide Cellular Translational Regulation in Lacticaseibacillus rhamnosus ATCC 53103 under Acid Stress" Foods 11, no. 10: 1411. https://0-doi-org.brum.beds.ac.uk/10.3390/foods11101411

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop