1. Introduction
Lactiplantibacillus plantarum (L. plantarum) is one of the most significant lactic acid bacteria (LAB) and is found in various environments, from different fermented foods[
1] to the human gastrointestinal tract [
2,
3]. Most published studies have focused on
L. plantarum as the starter culture in food fermentation.
L. plantarum strains have many functional properties, such as greater shelf life, flavor properties, nutritional quality, antioxidant activities, and antimicrobial activities [
4]. Low pH, which acts as a major environmental stress factor, significantly affects the growth and reproduction of
L. plantarum. Acids can cause several detrimental effects, such as enzyme inactivation, damage to DNA and proteins, and changes in the growth and metabolism of LAB [
5]. To cope with acid stress, LAB strains have developed various responses to stress in the natural environment, including the maintenance of pH homeostasis, alkali production, metabolic regulation, exopolysaccharide production, and macromolecule repair [
6]. High efficiency and accuracy of gene expression under complex and changeable real-world environments are preconditions for a strain to survive, develop, and even perform specific functions. Transcriptional regulator (TF) genes serve as probes to sense changes in environmental conditions and then adjust the expression level of genes inside the cell quickly and precisely; therefore, they are important for the survival of the strain in response to acid stress.
The tolerance of
L. plantarum to acidic pH is a popular research topic [
2,
7]. Thus, further investigations are needed to analyze the molecular mechanisms by which
L. plantarum is acid-tolerant.
Sugar metabolism is important in the context of abiotic stress [
8]. Sugars can finely regulate the balance between stress resistance and plant growth throughout the life cycle [
9]. Some soluble sugars produced in this pathway can regulate osmotic pressure, provide oxidation resistance, and protect the integrity of the cell membrane by increasing the concentration of cell fluid to enhance stress resistance [
10,
11]. Xu et al. reported that
A. frutescens reduces cold damage by increasing the accumulation of soluble sugar content in the cells [
12]. Another study revealed that drought affects sugar metabolism by increasing the activities of invertase and amylase [
13]. However, the molecular response mechanism underlying sugar metabolism in
L. plantarum under acid stress is unclear.
In our previous study [
14], we identified the OEOE_RS01070 gene, which encodes the AraR transcriptional regulator and harbors 85 nonsynonymous SNPs in
Oenococcus oeni mutants, exhibiting divergent acid resistance phenotypes. Its homolog, LP_RS14895, in
L. plantarum WCFS1 shares 49.41% amino acid identity. In this study, we investigated the regulatory mechanism of acid resistance mediated by the LP_RS14895 gene in
L. plantarum WCFS1.
2. Materials and Methods
2.1. Bacterial Strains, Plasmids, and Culture Conditions
The knockout and expression plasmids were maintained in
Escherichia coli DH5α (Tiangen, China) grown at 200 rpm on Luria-Bertani broth (LB).
Lactiplantibacillus plantarum WCFS1 was propagated statically at 37 °C in deMan Rogosa Sharpe (MRS) broth (HB0384-1, Qingdao Hope Bio-Technology Co., Ltd., China). The experiments were performed in MRSC (modified MRS medium containing 4 g/L yeast extract
, 2 g/L (NH
4)
2SO
4, 5 g/L NaCl, 1 mL/L Tween 80, 0.2 g/L MgSO
4·7H
2O, 0.04 g/L MnSO
4·4H
2O, and 50 mM/L monosaccharide (25 mM/L disaccharide)). The selected sugar sources [
15] included arabinose, ribose, glucose, mannose, fructose, and trehalose (Aladdin Biotechnology Co., Ltd. Shanghai), and the pH was adjusted to 3.8 with HCl. Appropriate antibiotics were added to the medium when needed for plasmid maintenance: 50 µg/mL kanamycin for
E. coli and 50 µg/mL erythromycin for
L. plantarum.
All bacterial strains and plasmids used and constructed in this study are listed in
Table S1. The plasmids were constructed via standard molecular cloning and PCR fusion (overlap extension) techniques. Restriction enzymes (XbaI and ApaI) were purchased from Takara. PCR was performed using a C1000 Touch PCR System (Bio-Rad) following standard procedures.
2.2. Construction of Lp_rs14895 Plasmids and Strains
Construction of LP_RS14895 knockout and expression plasmids
The skeleton of pLCNICK was obtained via double digestion with XbaI and ApaI [
16]. Two 1.0 kb fragments flanking LP_RS14895 (14895-up and 14895-down) were amplified from the genomic DNA of
L. plantarum WCFS1 using the primers 14895-up-1/14895-up-2 and 14895-down-1/14895-down-2, respectively. A 122 bp sgRNA framework that targets 14895 (14895-sgRNA) was obtained via PCR using the primers sgRNA-1/14895-sgRNA-2 with pLCNICK as the template. These three fragments were then assembled with 14895-up-1 and 14895-sgRNA-2 by overlap extension PCR, which yielded a new fragment, 14895-uds. The backbone of pLCNICK and the fragment 14895-uds were assembled to produce a new plasmid, pLCNICK-Δ14895, using the One-Step Cloning Kit. Next, positive clones were verified by PCR amplification using the primers pLCNICK-text-1 and pLCNICK-text-2.
The plasmid pMG36ek11 derived from a study by Yang et al. was used as a starting point and linearized by PCR amplification [
17]. Two fragments were connected to form the overexpression vector pMG36ek11-14895, the backbone of pMG36ek11 (amplified using primers P11-2 and pMG36e-Express-3), and the 14895 fragment (obtained by PCR using primers ek11-14895-1 and ek11-14895-2, with the
L. plantarum WCFS1 strain used as the template). Positive clones were verified by PCR amplification using the primers pMG36e-text-1 and pMG36e-text-2.
Transformation
Heat shock transformation was performed following the instructions for competent E. coli DH5α cells. After shaking at appropriate temperatures (pLCNICK-Δ14895 at 30 °C, pMG36ek11-14895 at 37 °C) and 200 rpm for 30 min, the corresponding antibiotic plates were coated. One colony was selected for colony PCR verification using the primers pLCNICK-test-1, pLCNICK-test-2, pMG36e-test-1, and pMG36e-test-2.
These plasmids were subsequently delivered into
L. plantarum WCFS1 via electroporation. Electrotransformed and electrocompetent
L. plantarum WCFS1 cells were prepared as described in other studies [
16,
17]. Electroporation was performed with a Gene Pulser X-cell (Bio-Rad, USA) and a 2 mm cuvette (BTX, USA) with the following parameters: 2 kV, 4 ms. Then, 1 mL of the recovered SMRS medium (MRS with 0.5 M sucrose and 0.1 M MgCl
2) was added to a cuvette, and the mixture was recovered within 4 h and then plated on MRS supplemented with erythromycin [
18].
Identification of mutants
The mutants were screened as described in another study, with some modifications [
19]. Colony PCR (cPCR) was performed to screen the correct mutants. Successful transformants were selected from MRS plates supplemented with erythromycin. Positive clones were verified via PCR amplification with the primers 14895-in-1 and 14895-in-2 or with the primers 14895-out-1 and 14895-out-2 to obtain the LP_RS14895 knockout mutants. Similarly, the primers pMG36e-test-1 and pMG36e-test-2 were used to verify the LP_RS14895-overexpressing mutants. The PCR products were sequenced to confirm the deletions and expression (Yangling Aoke Biotech Co. Ltd., Shanxi, China).
2.3. Acid Resistance
To evaluate acid resistance, WCFS1, WCFS1-Δ14895-pMG36ek11, and WCFS1-Δ14895-pMG36ek11-14895 were applied. The early exponential phase (OD
600 value of 0.3–0.5) was washed with 0.85% sodium chloride and inoculated at a 1% inoculaiton rate into MRS and MRSC liquid medium at pH 3.8, followed by culture for 40 h at 37 °C. Samples were collected at appropriate time points (
Figure 1A). Cell growth was monitored by measuring the optical density at 600 nm with a spectrophotometer (AOE instruments A380).
2.4. Comparison of Sugar Concentrations Between the Wild-Type and Mutant Strains in Mrsc Medium Via Hplc
We used WCFS1, WCFS1-Δ14895-pMG36ek11, and WCFS1-Δ14895-
pMG36ek11-14895 for the experiment. Overnight cultures were transferred to fresh MRS medium and cultured. When the OD600 reached 0.3, the cells were collected, washed three times with 0.85% (w/v) sterile saline, and then incubated in MRSC medium at pH 3.8 for 36 h. The supernatant was collected by centrifugation and filtered through a 0.45 µm filter membrane for further analysis.
The residual sugar content was analyzed by HPLC (1260 Infinity II, Agilent, United States) as described by Rossouw et al. [
20], with slight modifications. A 300 × 7.8 mm i.d. Aminex HPX-87H column (Bio-Rad, United States) was used with a column temperature of 60 °C and a flow rate of 0.6 mL/min. The mobile phase was 5 mM H
2SO
4, which was prepared by diluting reagent-grade H
2SO
4 with distilled water.
2.5. Real-Time Quantitative Pcr Analysis of the Lp_rs14895 Gene and Capsular Polysaccharide Biosynthesis
The strains were activated and collected as described above, and then, 1% (v/v) inoculum was introduced into MRS and MRSC media at different pH (3.8 and 6.2). Cells were collected for RNA extraction. Total RNA was extracted from the cell pellets using a miRNeasy kit (Qiagen, Hilden, Germany) with DNase digestion following the manufacturer’s protocol. The qualified RNA was stored at –80 °C.
The mRNA was extracted and reverse-transcribed to produce cDNA using a cDNA synthesis kit. The abundance of mRNAs was measured by amplifying the genes via the corresponding cDNAs as PCR templates. To conduct microRNA expression analysis, the miScript PCR starter kit (Tiangen, China) was used following the manufacturer’s instructions and measured on a CFX96 Touch instrument (Bio-Rad Laboratories, Hercules, CA, USA). PCR was performed with the following conditions: 50 cycles of 95 °C for 15 s and 60 °C for 1 min. The Ct (threshold value) values were compared to values obtained from the calibration sample (16S rRNA). The values for these genes were normalized to those of 16S rRNA to estimate the relative copy numbers of the genes. Relative gene expression was calculated using the 2
–ΔΔCt method [
18]. The primer sequences can be found in
Table S2, and the experiments were conducted three times.
2.6. Dna Affinity Purification Sequencing (Dap-Seq)
DAP-seq was conducted following the methods described in another study [
21]. After DAP DNA was extracted from the WCFS1 strain, the enriched DNA was fragmented into short fragments via ultrasonication. Next, the DNA fragments were end-repaired, 3’A was added, and then the fragments were ligated to Illumina sequencing adapters. DNA fragments of appropriate lengths (usually 100–300 bp, including the adapter sequence) were selected for PCR amplification. Finally, a qualified library was obtained for sequencing.
The DNA-binding site for specific transcription factors or histone modifications is not random, as they present a conserved DNA sequence pattern. The MEME suite (
http://meme-suite.org/) was used to detect the motifs. MEME (
http://meme suite.org
/tools/meme) and DREME (
http://memesuite.org/tools/dreme) were used to detect the sequence motif, which was determined to detect long and short consensus sequences.
2.7. Rna-Seq Analysis
The strains WCFS1 and WCFS1-Δ14895 were used. Overnight cultures were transferred into fresh MRS medium and incubated, followed by culture for 18 h at 37 °C, after which the cells were collected for RNA-seq analysis. Total RNA was extracted using TRIzol reagent. RNA-seq was performed at BerryGenomics (Guangzhou, China). The mRNA obtained was isolated using the NEBNext Poly (A) mRNA Magnetic Isolation Module (NEB, E7490), after which they were sequenced on an Illumina NovaSeq 6000 platform, and paired-end reads were generated. Transcriptome analysis was performed using reference genome-based read mapping. Gene expression levels were estimated using FPKM values software [
22]. Genes with
q < 0.05 and |log2-fold change| > 1.5 were identified as significantly differentially expressed genes (DEGs). Differential expression analysis between the WCFS1 and WCFS1-Δ14895 strains was performed using the software package DESeq2 [
23]. GO analysis of these DEGs was conducted by agriGO version 2.0 (
http://systemsbiology.cau.edu.cn/agriGOv2/). Venn diagrams were created using VENNY version 2.1 (
https://bioinfogp.cnb.csic.es/tools/venny/index.html).
3. Results
3.1. Lp_rs14895 Functions as an Arar Transcriptional Regulator to Increase Acid Stress Tolerance in L. Plantarum Wcfs1
In our previous study, we found the OEOE_RS01070 gene, which has 85 nonsynonymous SNPs identified through comparative genomic analysis of three strain groups with different acid-resistant phenotypes [
14]. However, efforts to detect the role of the OEOE_RS1070 gene in
O. oeni have failed because of technical limitations. Studies have indicated that the more similar the protein sequences are, the more likely they are to have similar functions [
24,
25]. Therefore, we considered the role of genes from the perspective of similar genes in other bacteria, especially
L. plantarum WCFS1, which also appears in wine and is a reference strain in other food industries. To test this hypothesis, we analyzed the gene associations in
L. plantarum WCFS1 and identified the gene LP_RS14895, which shares 49.41% amino acid identity with the OEOE_RE1070 gene. LP_RS14895 also shares 41.20% similarity with the arabinose repressor AraR in
Bacillus subtilis (PDB: 3tb6.1. A) [
26].
We then created several mutants, including a gene-deletion strain (WCFS1-Δ14895), a complementation strain (WCFS1-Δ14895-pMG36ek11-14895), and a deletion strain with an empty vector (WCFS1-Δ14895-pMG36ek11). The growth capacity matched our expectations, as indicated by the significantly lower OD
600 value of the knockout mutant WCFS1-Δ14895 than of the wild-type strain WCFS1 at pH 3.8 (
Figure 1). No significant difference in cell growth was observed between the WCFS1 strain and the complementation strain (WCFS1-Δ14895-pMG36ek11-14895). These findings indicated that the WCFS1-Δ14895 knockout mutant is more sensitive to acid stress than the wild-type strain, suggesting that the transcriptional regulator AraR enhances acid tolerance in
L. plantarum WCFS1.
3.2. Arar: A Central Regulatory Hub for Carbohydrate Metabolic Reprogramming Under Acid Stress
Some studies suggest that AraR has limited involvement in L-arabinose utilization [
27,
28]. However, given the existence of pH-dependent compensatory regulatory mechanisms in sugar metabolism pathways [
15], we deliberately expanded the scope of the substrate beyond conventional arabinase-related substrates. We selected arabinose (Ara), glucose (Glu), fructose (Fru), ribose (Rib), mannose (Man), and trehalose (Tre) to investigate the potential involvement of the AraR in the metabolic regulation of these sugars at pH 3.8.
These findings showed that supplementation significantly enhances the growth of strains under highly acidic conditions (pH 3.8) (
Figure 2A). The growth of WCFS1-Δ14895-pMG36ek11-14895 and the wild-type strain WCFS1 in the presence of carbohydrates was significantly greater than that of the mutant strain WCFS1-Δ14895-pMG36ek11. The only prominent difference between WCFS1-Δ14895-pMG36ek11-14895 and WCFS1 was observed in the presence of Ara.
The HPLC data (
Figure 2B) showed the residual sugar levels after 36 h, revealing that the WCFS1-Δ14895-pMG36ek11-14895 strain metabolizes sugars more effectively than the WCFS1-Δ14895-pMG36ek11 knockout strain, except when Glu is added. These findings revealed that LP_RS14895 plays an important role in carbohydrate metabolism at low pH.
The RT-qPCR analysis (
Figure 2C) revealed that the expression of the LP_RS14895 gene was more than four times greater at pH 3.8 than at pH 6.2 and increased 60.79-fold after Rib was added, further confirming its role in carbohydrate metabolism under acidic conditions.
The experimental results revealed that LP_RS14895, which is significantly upregulated (>4-fold) under low-pH conditions, drives efficient carbohydrate metabolism, thereby conferring a survival advantage to the strain in extremely acidic environments.
3.3. Rna-Seq Analysis Revealed That the Arar Regulates Various Carbohydrate Pathways
Transcriptome analysis of the WT and WCFS1-Δ14895 mutant strains based on RNA-seq was performed at pH 3.8. A total of 40 DEGs were obtained, of which 29 genes were upregulated, and 11 genes were downregulated (
Table S3). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of these genes revealed significant enrichment in the pentose and glucuronate interconversions and amino sugar and nucleotide sugar metabolism pathways (
Figure 3A). Additionally, gene set enrichment analysis (GSEA) was performed to further characterize the transcriptome changes caused by the LP_RS14895-deficient mutant. The results revealed different expression patterns between the wild-type and WCFS1-Δ14895 mutant strains, with five gene sets significantly enriched in the mutant strain (FDR<0.25 and nominal p < 0.05) (
Figure 3B). Moreover, pathways related to sugar metabolism, including galactose metabolism, starch and sucrose metabolism, and amino sugar and nucleotide sugar metabolism, were significantly enriched.
3.4. Dap-Seq Provides A Global Overview of Potential Tf Binding Sites Genome-Wide
To define the direct transcriptional targets of AraR in response to acid stress, we subsequently performed DAP-seq to identify the binding sites on genomic DNA [
29]. All reads were counted in the 2000 bp interval upstream of the TSS and 2000 bp downstream of the trans-TES using the DeepTools software. A total of 1279 peaks were identified, which presented the greatest distribution in the adjacent region of the TSS (
Figure 4A) and were distributed mostly in the promoter and exon regions, accounting for 41.9% and 57.0%, respectively (
Figure 4B). These findings confirmed that LP_RS14895 is a typical TF with DNA binding ability and gene regulatory activity. Then, the sequences around the peaks were analyzed via motif elicitation (MEME) to detect significant motif sequences in the peak sequences [
30]. The results revealed that the most significantly enriched motif sequence was “ARCCMATMAHC” (motif 1, E value = 8.7e+001). Three additional significantly enriched LP_RS14895 binding sites were also identified. The sequences were “TGGKGGCW” (motif 2, E value = 5.5e+003), “CRTCGGGTTGCAGG” (motif 3, E value = 1.1e+005), and “GCCGTTACCGC” (motif 4, E value = 2.6e+005) (
Figure 4C). The KEGG cluster analysis of the potential target genes revealed that starch and sucrose metabolism, o-antigen nucleotide sugar biosynthesis, amino sugar, and nucleotide sugar metabolism, polyketide sugar unit biosynthesis and Gal metabolism (involved in sugar metabolism) were enriched among the top biological processes (
Figure 4D).
3.5. Integrated Multiomics Identification of Lp_rs14895 Potential Targets Under Acid Stress
To systematically evaluate the regulatory role of LP_RS14895 in the acid resistance mechanisms of
L. plantarum WCFS1, we integrated DAP-seq, RNA-seq and GSEA to identify putative target genes under acid stress. Crossomics analysis revealed four genes with conserved binding motifs in promoter regions (
Table S4), among which three exhibited significant upregulation in the ΔLP_RS14895 mutant (including LP_RS14880, LP_RS14950, and LP_RS12915, and the log2(FC) values were 5.69, 2.39, and 1.08, respectively). We also conducted DAP-seq and GSEA. Through these analyses, we identified 27 genes interacting with the transcription factor LP_RS14895 (
Table S5). The 31 candidate target genes are likely to represent the primary transcriptional targets through which LP_RS14895 mediates the rapid adaptation of
L. plantarum to lethal acidic environments, warranting further investigation.
4. Discussion
We emphasized the importance of AraR transcriptional regulator, in the acid tolerance of L. plantarum WCFS1 by conducting DAP-seq and RNA-seq profiling. Our findings offered new insights into the acid resistance regulatory network in LAB and identified a potential genetic engineering target to improve acid tolerance in industrial fermentation strains.
4.1. Arar Enhances Microbial Adaptation in Acidic Environments: Transitioning From Arabinose Regulator to Ph Modulator
The regulation of microbial acid tolerance is a complex process governed by various factors, especially transcription factors, which play a crucial role in enabling microorganisms to survive and thrive in acidic environments. This regulation is critical for both pathogenic and industrially relevant microorganisms.
In the context of industrial applications, the engineering of transcription factors has been employed to enhance acid tolerance in various microorganisms. For example, in
Lactococcus lactis, the positive regulation of the DLT operon by the two-component system transcriptional regulator TCSR7 has been shown to enhance acid tolerance. This regulation increases the positive charge on the cell membrane surface, thereby improving the bacterium's resistance to acidic conditions [
31]. Furthermore, the role of transcription factors in acid tolerance is not limited to bacteria. In yeast, such as
Saccharomyces cerevisiae, transcription factors have been identified as key players in conferring tolerance to acetic acid, a common stressor in biofuel production. The use of a zinc-finger-based artificial transcription factor library has led to the identification of novel genes involved in acetic acid tolerance, highlighting the importance of transcriptional regulation in yeast stress responses [
32]. Additionally, the optimization of transcription factor expression levels through strategies like cocktail δ-integration has been shown to improve lactic acid tolerance in yeast, further demonstrating the versatility of transcription factor manipulation in enhancing microbial stress tolerance [
33].
In our study, LP_RS14895 in
L. plantarum WCFS1 and OEOE_RS01070 in
O. oeni both encode an AraR transcriptional regulator that has historically been viewed as a narrow-domain arabinose repressor [
27,
28]. By integrating RNA-seq-guided transcriptomics with growth-kinetic profiling, we now position AraR as a global 'acid-fitness' determinant (
Figure 1). Evolutionary evidence is provided by the
O. oeni mutants, in which three groups of acid-tolerant variants, differentiated by 85 non-synonymous SNPs within the OEOE_RS01070 gene, exhibit a precise genotype-to-phenotype correlation at pH 3.0 [
14]. Moreover, our results support the hypothesis that the OEOE_RS01070 gene contributes to the distinct acid-resistance phenotypes observed among the three sets of strains. Collectively, these results bolster the hypothesis that conserved protein sequence motifs are indicative of conserved regulatory functions across bacterial species [
24,
25].
AraR transcription factors are central to the regulation of microbial acid tolerance, with significant implications for industrial applications. The ability to modulate AraR through genetic engineering offers promising opportunities for developing probiotic therapies and optimizing microbial strains for industrial processes.
4.2. Arar Governs Acid-Driven Sugar Utilization in L. Plantarum
The inhibition of sugar utilization by microbes under high-acid conditions is a critical area of study, particularly in understanding microbial metabolism and its implications for various biological processes [
34,
35]. Our results demonstrated that AraR plays a significant role in influencing the utilization of carbohydrates such as glucose (Glu), fructose (Fru), ribose (Rib), mannose (Man), and trehalose (Tre) at low pH (
Figure 2). Additionally, RNA-seq revealed significant enrichment in carbohydrate metabolism pathways, highlighting the central role of AraR in sugar metabolic reprogramming (
Table S3,
Figure 3A). These findings suggest that AraR represents a key target for elucidating acid-sugar metabolic crosstalk, which could be targeted to alleviate metabolic inhibition induced by low pH.
The phenomenon where acids inhibit sugar utilization is not unique to high-acid environments but is also observed in other contexts, such as in rhizobia, a class of symbiotic diazotrophic bacteria. In rhizobia, the presence of C4 acids like succinate leads to a preference for these acids over sugars, a process known as Succinate Mediated Catabolite Repression (SMCR) [
36]. This mechanism of catabolite repression is crucial as it dictates the hierarchy of carbon source utilization, which is a fundamental aspect of microbial metabolism.
The study of AraR-mediated enhancement of sugar utilization at low pH in L. plantarum WCFS1 offers a valuable model for figuring out how to solve inhibitory processes under high-acid conditions. The parallels in metabolic regulation and genetic control mechanisms underscore the broader applicability of these findings. Further research into the regulatory pathways involved in catabolite repression/enhancement will improve our understanding of microbial adaptation and could lead to innovative strategies for managing microbial processes in various industrial and environmental settings.
4.3. Arar: A Promising Target for Engineering Microbial Carbohydrate Biocatalysts
The exploration of transcription factors as promising targets for engineering microbial carbohydrate biocatalysts is a burgeoning field that holds significant potential for advancing biotechnological applications.
Transcription factors are integral to the hierarchical regulatory networks that govern microbial gene expression. Engineered variants of the XylR protein (R121C/P363S) and increased intracellular concentrations of this regulator have been demonstrated to eliminate glucose- and arabinose-mediated repression, thereby enabling the concurrent fermentation of xylose-arabinose mixtures in
E. coli [
37,
38]. This exemplifies how the reprogramming of a single transcription factor can remove diauxic growth and enhance the conversion of mixed sugars. Liu et al. (2013) highlight the potential of engineering specific, middle-level, and global regulators to improve microbial resilience to environmental stresses, thereby enhancing their utility in industrial processes [
39]. Additionally, transcription factors play a critical role in metabolic regulation, which is essential for the microbial production of fine chemicals, including organic acids. The work by Liu et al. (2020) underscores the importance of TFs in constructing biosensors that can monitor intracellular metabolite concentrations, facilitating high-throughput strain evolution and optimization of metabolic pathways [
40]. Moreover, the application of transcription factor-based biosensors in biotechnology further exemplifies their potential in microbial engineering. As reviewed by Dietrich et al. (2015), TF-based biosensors can be integrated into synthetic regulatory circuits to control gene expression in response to specific stimuli, enabling dynamic regulation of metabolic pathways [
41].
Similarly, our findings identify the AraR protein as a comparable metabolic switch in L. plantarum. Under conditions of low pH stress, this homolog globally activates genes involved in the utilization of glucose, fructose, ribose, mannose, and trehalose, and it also coordinates membrane-proton homeostasis. Therefore, the targeted manipulation of AraR or its orthologs offers a promising avenue for the development of efficient microbial carbohydrate biocatalysts.
5. Conclusions
This study demonstrates that the AraR transcriptional regulator serves as both the central 'commander' of the acid-tolerance network and a global modulator of carbohydrate metabolism in L. plantarum WCFS1 strain. The findings provide a theoretical framework for elucidating microbial acid-stress regulation, optimizing carbohydrate utilization in acidic environments, and guiding the targeted engineering of robust industrial strains.
Supplementary Materials
The following supporting information can be downloaded at the website of this paper posted on
Preprints.org. Supplemental materials are available online only.
Supplementary Table S1. Bacterial strains and plasmids used in this study.
Supplementary Table S2. Primers used in this study.
Supplementary Table S3. Genes that are differentially expressed between the WCFS1-Δ14895 mutant and wild-type WCFS1 at pH 3.8.
Supplementary Table S4. Predicting and screening target genes with binding sites for the transcription factor LP_RS14895 in the promoter region via DAP-seq and RNA-seq.
Supplementary Table S5. Predicting and screening target genes with binding sites for transcription factor LP_RS14895 in the promoter region using DAP-Seq and GSEA.
Author Contributions
Conceptualization, C.Q.; Methodology, Z.L. C.M. F.C.and P.T.; Software, C.M. P.T.; Validation, Z.L. C.M. and F.C.; Formal Analysis, Z.L. C.M. F.C.and P.T; Investigation, Z.L. C.M. F.C.and P.T; Resources, Z.L. C.M. F.C.and P.T; Data Curation, Z.L. C.M. and C.Q.; Writing-Original Draft Preparation, Z.L. and C.M.; Writing-Review&Editing, C.Q.; Visualization, Z.L. and C.M.; Supervision, C.Q.; Project Administration, C.Q.; Funding Acquisition, C.Q.
Funding
This work was financially supported by the National Natural Science Foundation of China Youth Science Fund Project (No. 32302049) and the Key Research and Development Program of the Xinjiang Uygur Autonomous Region (No. 2023B02029).
Ethics approval and consent to participate
No applicable.
Informed Consent Statement
All authors have read and agreed with the manuscript.
Conflicts of Interest
The authors declare that they have no competing interests.
References
- Devi, A.; Anu-Appaiah, K.A.; Lin, T.F. , Timing of inoculation of Oenococcus oeni and Lactobacillus plantarum in mixed malo-lactic culture along with compatible native yeast influences the polyphenolic, volatile and sensory profile of the Shiraz wines. Lwt-Food Sci. Technol. 2022, 158. [Google Scholar] [CrossRef]
- Nath, S.; Sikidar, J.; Roy, M.; Deb, B. , Screening of probiotic properties of isolated from fermented milk product. Food Qual. Saf. 2020, 4, 213–223. [Google Scholar] [CrossRef]
- Liu, H.F.; Li, P.H.; Xin, J.L.; Huang, H.C.; Yang, Y.X.; Deng, H.C.; Zhou, Z.Y.; Zhong, Z.J.; Peng, G.N.; Chen, D.C.; He, C.L. , Probiotic characteristics and whole genome analysis of PM8 from giant panda milk. Probiotics Antimicro. 2025. [Google Scholar]
- Echegaray, N.; Yilmaz, B.; Sharma, H.; Kumar, M.; Pateiro, M.; Ozogul, F.; Lorenzo, J.M. , A novel approach to From probiotic properties to the omics insights. Microbiolo. Res. 2023, 268. [Google Scholar] [CrossRef] [PubMed]
- Lyu, C.; Zhao, W.; Peng, C.; Hu, S.; Fang, H.; Hua, Y.; Yao, S.; Huang, J.; Mei, L. , Exploring the contributions of two glutamate decarboxylase isozymes in Lactobacillus brevis to acid resistance and gamma-aminobutyric acid production. Microb. Cell Fact. 2018, 17, 180. [Google Scholar] [CrossRef]
- Mallick, S.; Das, S. , Acid-tolerant bacteria and prospects in industrial and environmental applications. Appl. Microbiol. Biotechnol. 2023, 107, 3355–3374. [Google Scholar] [CrossRef] [PubMed]
- Tang, W.; Li, C.; He, Z.; Pan, F.; Pan, S.; Wang, Y. , Probiotic properties and cellular antioxidant activity of Lactobacillus plantarum MA2 isolated from Tibetan kefir grains. Probiotics Antimicro. 2018, 10, 523–533. [Google Scholar] [CrossRef]
- Thalmann, M.; Santelia, D. , Starch as a determinant of plant fitness under abiotic stress. New Phytol. 2017, 214, 943–951. [Google Scholar] [CrossRef]
- Jeandet, P.; Formela-Luboinska, M.; Labudda, M.; Morkunas, I. , The role of sugars in plant responses to stress and their regulatory function during fevelopment. Int. J. Mol. Sci. 2022, 23. [Google Scholar] [CrossRef]
- La, V.H.; Lee, B.R.; Islam, M.T.; Park, S.H.; Lee, H.; Bae, D.W.; Kim, T.H. , Antagonistic shifting from abscisic acid- to salicylic acid-mediated sucrose accumulation contributes to drought tolerance in. Environ. Exp. Bot. 2019, 162, 38–47. [Google Scholar] [CrossRef]
- Cao, S.F.; Yang, Z.F.; Zheng, Y.H. , Sugar metabolism in relation to chilling tolerance of loquat fruit. Food Chem. 2013, 136, 139–143. [Google Scholar] [CrossRef]
- Xu, H.Y.; Li, J.J.; Wang, L.J.; Li, X.Y.; Liu, Y.Q.; Wang, X.; Gao, T.T.; Ma, Y.P. , Integrated transcriptomic and metabolomics analysis reveals abscisic acid signal transduction and sugar metabolism pathways as defense responses to cold stress in Argyranthemum frutescens. Environ. Exp. Bot. 2023, 205. [Google Scholar] [CrossRef]
- Han, Y.Y.; Xu, T.Y.; Chen, H.; Tang, M. , Sugar metabolism and 14-3-3 protein genes expression induced by arbuscular mycorrhizal fungi and phosphorus addition to response drought stress in. J. Plant Physiol. 2023, 288. [Google Scholar] [CrossRef]
- Chen, Q.L.; Hao, N.; Zhao, L.L.; Yang, X.K.; Yuan, Y.X.; Zhao, Y.Z.; Wang, F.; Qiu, Z.B.; He, L.; Shi, K.; Liu, S.W. , Comparative functional analysis of malate metabolism genes in and at low pH and their roles in acid stress response. Food Res. Int. 2022, 157. [Google Scholar] [CrossRef] [PubMed]
- Cibrario, A.; Peanne, C.; Lailheugue, M.; Campbell-Sills, H.; Dols-Lafargue, M. , Carbohydrate metabolism in : a genomic insight. Bmc Genomics 2016, 17. [Google Scholar] [CrossRef]
- Huang, H.; Song, X.; Yang, S. , Development of a RecE/T-Assisted CRISPR-Cas9 Toolbox for. Biotechnol. J. 2019, 14. [Google Scholar] [CrossRef]
- Yang, X.; Zhao, L.; Chen, Q.; Wang, N.; Shi, K.; Liu, S. , Functional verification of the citrate transporter gene in a wine lactic acid bacterium, Lactiplantibacillus plantarum. Front Bioeng. Biotechnol. 2022, 10, 894870. [Google Scholar] [CrossRef]
- Meng, Q.; Yuan, Y.X.; Li, Y.Y.; Wu, S.W.; Shi, K.; Liu, S.W. , Optimization of electro transformation parameters and engineered promoters for from Wine. Acs Synth Biol. 2021, 10, 1728–1738. [Google Scholar] [CrossRef] [PubMed]
- Wang, G.Q.; Yu, H.I.; Feng, X.; Tang, H.Y.; Xiong, Z.Q.; Xia, Y.J.; Ai, L.Z.; Song, X. , Specific bile salt hydrolase genes in Lactobacillus plantarum AR113 and relationship with bile salt resistance. Lwt-Food Sci. Technol. 2021, 145. [Google Scholar] [CrossRef]
- Rossouw, D.; Du Toit, M.; Bauer, F.F. , The impact of co-inoculation with Oenococcus oeni on the trancriptome of Saccharomyces cerevisiae and on the flavour-active metabolite profiles during fermentation in synthetic must. Food Microbiol. 2012, 29, 121–31. [Google Scholar] [CrossRef]
- Wang, W.J.; Li, T.; Chen, Q.; Yao, S.X.; Zeng, K.F. , Transcriptional regulatory mechanism of a variant transcription factor CsWRKY23 in citrus fruit resistance to. Food Chem. 2023, 413. [Google Scholar] [CrossRef]
- Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. , Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–5. [Google Scholar] [CrossRef]
- Love, M.I.; Huber, W.; Anders, S. , Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15. [Google Scholar] [CrossRef]
- Li, H.Y.; Liu, T.R.; Wang, B.; Li, H.B. , Genome-wide characterization and expression analysis of GATA transcription factors in response to methyl jasmonate in. Genes 2022, 13. [Google Scholar] [CrossRef]
- Nam, H.J.; Han, S.K.; Bowie, J.U.; Kim, S. , Rampant exchange of the structure and function of extramembrane domains between membrane and water soluble proteins. PLoS Comput. Biol. 2013, 9, e1002997. [Google Scholar] [CrossRef] [PubMed]
- Procházková, K.; Cermáková, K.; Pachl, P.; Sieglová, I.; Fábry, M.; Otwinowski, Z.; Rezácová, P. , Structure of the effector-binding domain of the arabinose repressor AraR from Bacillus subtilis. Acta Crystallogr. D. 2012, 68, 176–185. [Google Scholar] [CrossRef] [PubMed]
- Kuge, T.; Teramoto, H.; Inui, M. , AraR, an L-arabinose-responsive transcriptional regulator in ATCC 31831, exerts different degrees of repression depending on the location of its binding sites within the three target promoter regions. J. Bacteriol. 2015, 197, 3788–3796. [Google Scholar] [CrossRef]
- Mota, L.J.; Tavares, P.; Sa-Nogueira, I. , Mode of action of AraR, the key regulator of L-arabinose metabolism in Bacillus subtilis. Mol. Microbiol. 1999, 33, 476–89. [Google Scholar] [CrossRef] [PubMed]
- Bartlett, A.; O'Malley, R.C.; Huang, S.S.C.; Galli, M.; Nery, J.R.; Gallavotti, A.; Ecker, J.R. , Mapping genome-wide transcription-factor binding sites using DAP-seq. Nat. Protoc. 2017, 12, 1659–1672. [Google Scholar] [CrossRef]
- Machanick, P.; Bailey, T.L. , MEME-ChIP: motif analysis of large DNA datasets. Bioinformatics 2011, 27, 1696–7. [Google Scholar] [CrossRef]
- Wu, H.; Zhang, Y.L.; Li, L.; Li, Y.N.; Yuan, L.; Yue, E.; Qiao, J.J. , Positive regulation of the operon by TCSR7 enhances acid tolerance of F44. J. Dairy Sci. 2022, 105, 7940–7950. [Google Scholar] [CrossRef]
- Ma, C.; Wei, X.W.; Sun, C.H.; Zhang, F.; Xu, J.R.; Zhao, X.Q.; Bai, F.W. , Improvement of acetic acid tolerance of using a zinc-finger-based artificial transcription factor and identification of novel genes involved in acetic acid tolerance. Appl. Microbiol. Biotechnol. 2015, 99, 2441–2449. [Google Scholar] [CrossRef]
- Yamada, R.; Kumata, Y.; Mitsui, R.; Matsumoto, T.; Ogino, H. , Improvement of lactic acid tolerance by cocktail δ-integration strategy and identification of the transcription factor PDR3 responsible for lactic acid tolerance in yeast. World J. Microbiol. Biotechnol. 2021, 37. [Google Scholar] [CrossRef]
- Lu, Y.Q.; Song, S.C.; Tian, H.X.; Yu, H.Y.; Zhao, J.X.; Chen, C. , Functional analysis of the role of CcpA in grown on fructooligosaccharides or glucose: a transcriptomic perspective. Microb. Cell Fact. 2018, 17. [Google Scholar] [CrossRef]
- Xie, S.; Ma, J.Z.; Lu, Z. , Enhances oxidative stress tolerance through rhamnose-dependent mechanisms. Front. Microbiol. 2024, 15. [Google Scholar] [CrossRef]
- Iyer, B.; Rajput, M.S.; Jog, R.; Joshi, E.; Bharwad, K.; Rajkumar, S. , Organic acid mediated repression of sugar utilization in rhizobia. Microbiol. Res. 2016, 192, 211–220. [Google Scholar] [CrossRef]
- Wei, W.P.; Shang, Y.Z.; Zhang, P.; Liu, Y.; You, D.; Yin, B.C.; Ye, B.C. , Engineering prokaryotic transcriptional activator XylR as a xylose-inducible biosensor for transcription activation in Yeast. Acs Synth. Biol. 2020, 9, 1022–1029. [Google Scholar] [CrossRef]
- Martinez, R.; Flores, A.D.; Dufault, M.E.; Wang, X. , The XylR variant (R121C and P363S) releases arabinose-induced catabolite repression on xylose fermentation and enhances coutilization of lignocellulosic sugar mixtures. Biotechnol. Bioeng. 2019, 116, 3476–3481. [Google Scholar] [CrossRef] [PubMed]
- Lin, Z.L.; Zhang, Y.; Wang, J.Q. , Engineering of transcriptional regulators enhances microbial stress tolerance. Biotechnol. Adv. 2013, 31, 986–991. [Google Scholar] [CrossRef] [PubMed]
- Li, J.W.; Zhang, X.Y.; Wu, H.; Bai, Y.P. , Transcription factor engineering for high-throughput strain evolution and organic acid bioproduction: a review. Front. Bioeng. Biotech. 2020, 8. [Google Scholar] [CrossRef] [PubMed]
- Mahr, R.; Frunzke, J. , Transcription factor-based biosensors in biotechnology: current state and future prospects. Appl. Microbiol. Biotechnol. 2016, 100, 79–90. [Google Scholar] [CrossRef] [PubMed]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).