Molecular adaptation of Lactobacillus plantarum WCFS1 to gallic acid revealed by genome-scale transcriptomic signature and physiological analysis
© Reverón et al. 2015
Received: 16 April 2015
Accepted: 23 September 2015
Published: 9 October 2015
Gallic acid (GA) is a model hydroxybenzoic acid that occurs esterified in the lignocellulosic biomass of higher plants. GA displays relevant biological activities including anticancer properties. Owing to its antimicrobial and cellulase-inhibiting activities, GA also imposes constraints to the fermentability of lignocellulosic hydrolysates. In depth-knowledge of the mechanisms used by tolerant microorganisms to adapt to hydroxybenzoic acids would be a step forward to improve the bioavailability of GA or select/engineer production hosts with improved metabolic traits for the bioconversion of pretreated lignocellulosic biomass.
Whole genome transcriptional profiling using DNA microarrays was used to characterize the molecular response of Lactobacillus plantarum WCFS1 to GA. Expression levels of 14 and 40 genes were differentially regulated at 1.5 and 15 mM GA, respectively. The transcriptomic analysis identified a marked induction of genes with confirmed or related roles to gastrointestinal survival, the repression of genes coding for certain ABC-type transporters and modulation of genes involved in the control of intracellular ammonia levels, among other responses. Most notably, a core set of genes dedicated to produce GA from polyphenols (tanB Lp ), decarboxylate GA to pyrogallol (lpdB, lpdC and lpdD) and transport functions (lp_2943) was highly overexpressed at both GA concentrations. Correspondingly, resting cells of strain WCFS1 induced by GA, but not their non-induced controls, produced pyrogallol. Gene expression and organization of genes involved in GA metabolism suggested a chemiosmotic mechanism of energy generation. Resting cells of L. plantarum induced by GA generated a membrane potential and a pH gradient across the membrane immediately upon addition of GA. Altogether, transcriptome profiling correlated with physiological observations indicating that a proton motive force could be generated during GA metabolism as a result of electrogenic GA uptake coupled with proton consumption by the intracellular gallate decarboxylase.
The combination of transcriptome and physiological analyses revealed versatile molecular mechanisms involved in the adaptation of L. plantarum to GA. These data provide a platform to improve the survival of Lactobacillus in the gut. Our data may also guide the selection/engineering of microorganisms that better tolerate phenolic inhibitors present in pretreated lignocellulosic feedstocks.
Edible plants are the main sources of dietary phenolic acids, particularly hydroxycinnamic and hydroxybenzoic acids, which are micronutrients with significant biological properties. Gallic acid (GA) is a model hydroxybenzoic acid widely distributed in edible plants and occurs in several legumes, fruits, vegetables, nuts and beverages of plant origin . GA displays relevant biological activities including anti-inflammatory  or antiviral ones  and has attracted interest because of its reported anticancer properties in animal and in vitro studies (reviewed in ). In the context of symbiosis with human host, the gut microbiota has evolved biochemical pathways for the bioactivation/degradation of dietary polyphenols . Hence, the bioaccesibility and health effects of GA depend on the activity of the subset of gut microbes that resist the anti-bacterial activity reported for tannins  and GA itself . Therefore, a better understanding of the tolerance of gut microbiota to these health-relevant compounds will be a step forward to improve the persistence of beneficial microorganisms in the intestine.
In addition to its significant biological activities, the GA esterified in the lignocellulosic biomass can be released during the pretreatment required for the efficient industrial uses of cellulosic materials by microorganisms . On account of its anti-microbial activity  , released GA and tannins (precursors of GA) have inhibitory effects on these microorganisms. In addition, these compounds are phenolic inhibitors of enzymes used for cellulose hydrolysis to produce ethanol  and consequently impose constraints to the fermentability of lignocellulosic hydrolysates. Therefore, in depth-knowledge of the mechanisms used by microorganisms to overcome the toxicity of hydroxybenzoic acids is of interest to assist in the selection and/or engineering of production hosts with improved metabolic traits for the bioconversion of lignocellulosic biomass.
Lactobacilli are suitable model microorganisms to study resistance mechanisms to phenolic acids as they display higher tolerance to phenolic compounds than other bacterial groups . In fact, some Lactobacillus ssp. have been chosen as models to obtain datasets of specific expression profiles in response to model hydroxycinnamic acids such as ferulic  and p-coumaric acids . Regarding the hydroxybenzoic acids, a proteomic approach revealed some of the molecular adaptive responses to the GA precursor, tannic acid . Among the molecular mechanisms used to persist under tannic acid-stress Lactobacillus plantarum relies on tannase (tannin acyl hydrolase) , an enzyme that transform the gallate esters of tannins into GA and glucose. Recently, the elusive gallate decarboxylase activity (GDC), which decarboxylates GA to yield pyrogallol (PG) as final product of tannin metabolism, has been uncovered in L. plantarum WCFS1 . Despite this crucial advance in the understanding of GA metabolism, knowledge on how gut microorganisms respond to hydroxybenzoic acids is not completely understood.
To provide insight into the microbial mechanisms involved in the tolerance to hydroxybenzoic acids, the current work describes the molecular adaptive responses of the model bacterium L. plantarum WCFS1 to GA as studied by whole-genome transcription profiling. Based on this transcriptional analysis, several mechanisms involved in the response to GA are proposed. The main response identified by the transcriptional datasets, the GA-inducible catabolism of GA to PG, was corroborated by specific metabolic analysis. The transcriptome-based results and the organization of genes involved in GA decarboxylation pointed towards a chemiosmotic mechanism of energy generation associated to GA metabolism, which was experimentally supported by membrane potential and internal pH measurements.
Global transcriptomic responses during adaptation to GA
To investigate the adaptive response of L. plantarum WCFS1 to GA, the transcriptomic profile of L. plantarum WCFS1 was defined in cells exponentially growing in medium devoid of GA after 10 min of exposure to 1.5 or 15 mM of this compound. The time of exposure was chosen considering the short half-life of mRNAs reported for genes involved in stress responses induced by phenolic acids in L. plantarum . The concentrations of GA used (1.5 and 15 mM GA) cover a range which could be representative of the amounts of GA present in the diet, provided that an estimated dietary intake of ≈6 mmol (1 g) GA/day has been reported by some authors .
The impact of GA on the transcriptomic profile of L. plantarum WCFS1 was evaluated by sorting all genes whose transcript level showed changes (log2ratio) of at least ±1.5 (p < 0.05). Overall, 14 genes were significantly differentially expressed at 1.5 mM GA whereas 40 transcripts were affected at 15 mM GA (15 upregulated and 25 downregulated). The differentially expressed genes, functionally distributed according to their Clusters of Orthologous Groups (COGs) categories, are shown (http://www.ncbi.nlm.nih.gov/COG/ and ftp://ftp.ncbi.nih.gov/pub/COG/COG2014/static/lists/listCOGs.html) (Additional file 1: Table S1). The induction of global stress responses was not detected in the transcriptome analysis. In addition, no downregulation of genes coding for ribosomal proteins or translation factors or associated with transcription was evident (Additional file 1: Table S1), indicating that the stringent response was not triggered by GA.
A core set of genes, including lpdB (lp_0271), lpdD (lp_0272), lpdC (lp_2945), lp_2943 (ion transporter), lp_2956 (tanB Lp ) (tannase), lp_2940 (surface protein which has been reportedly demonstrated to play a key role in the persistence and survival of L. plantarum WCFS1 in the GI-tract of mice ) and lp_0274 (transcriptional regulator), were highly overexpressed (Additional file 1: Table S1). In addition, the gene amtB (NH4 + transport protein involved in regulation of nitrogen metabolism) was downregulated. These genes were considered as the pivotal response to GA, as their expression showed the same trend and was roughly conserved at both GA concentrations.
Beside this response, other responses involving the carbohydrate and nitrogen metabolisms were observed at 15 mM GA. At this higher GA concentration several L. plantarum genes coding for ABC-type transporters were significantly downregulated (see below), whereas only two genes putatively involved in stress response pathways were upregulated (Additional file 1: Table S1). In addition, some regulatory networks were activated. These and other variations in the transcriptomic response of L. plantarum to GA are detailed in the following sections.
Correlation between gene expression, GA metabolism and the generation of a proton motive force
GA-mediated induction of genes associated with the transport and metabolism of GA
The GA-responsive genes most strongly upregulated in this study code for LpdC (lp_2945) (272.1-fold induction at 15 mM GA), LpdB (lp_0271) (26.3-fold induction at 15 mM GA) and LpdD (lp_0272) (24.0-fold induction at 15 mM GA). These are the three subunits of the recently uncovered intracellular gallate decarboxylase (GDC) of L. plantarum WCFS1  which decarboxylates GA into PG. The second most upregulated gene upon exposure to GA, lp_2943 (73-fold induction at 15 mM GA), was the sole gene in this study displaying homology with genes coding transporters of the IT (ion transporter) superfamily, which suggest the involvement of the lp_2943 product in GA transport. This transcript profile strongly suggests that the presence of GA is required to induce the transport of GA and its intracellular decarboxylation to PG.
Correlation between GA metabolism and transcriptome profiles
Generation of a proton motive force during GA metabolism
The genetic organization displayed by lp_2943, which clusters downstream lp_2945 (catalytic subunit of the GDC activity ), resembled that of some decarboxylation pathways built from a precursor/product exchanger and a decarboxylase that are activated to supply chemiosmotic energy by the generation of a proton motive force (PMF) across the membrane . This observation, besides the strong induction of lp_2943 and genes coding for the GDC mediated by GA, led us to ask whether a PMF was generated during GA metabolism.
The membrane potential was qualitatively measured with the fluorescent probe DiSC3. Variations in the ∆Ψ showed the same trend as those obtained for the transmembrane pH gradient. Thus, addition of 1.5 mM GA only increased the ∆Ψ in resting cells induced with GA (Fig. 2d) while non-induced resting cells did not generate a ∆Ψ upon the same conditions (Fig. 2b).
The experiments show that L. plantarum generates a proton motive force during GA transport and metabolism.
Effects of GA on central metabolic functions
Microarray data showed the upregulation (2.0-fold) of genes encoding for glucosamine-6-phosphate deaminase/isomerase (nagB), which is involved in N-acetyl glucosamine (GlcNAc) utilization. GA also induced the transcription of lp_2954 (~3.0-fold) which, according to LocateP (a genome-scale subcellular-location predictor ), codes for a hypothetical membrane protein of unknown function. Since the structure of a protein is better conserved than the amino acid sequence during evolution, the putative lp_2954 product was submitted to the Phyre2 modelling server (http://www.sbg.bio.ic.ac.uk/~phyre/) to probe the likely extent of structural similarity to other proteins. This search showed an integral membrane protein EIIC of a PTS system which transports disaccharides with predicted specificity for N, N´-diacetylchitobiose ([GlcNAc]2), as the best hit of the predicted protein. The lp_2954 gene coding for this putative PTS transporter lacks the accompanying genes coding for the structural components EIIA and EIIB, suggesting that it is one of the so called “orphan” PTS-EIIC transport components. The coincident upregulation of nagB and lp_2954 suggests that GlcNAc utilization is promoted by the presence of GA.
Genes coding for a ribose transport protein (lp_3659) and a ribose kinase (lp_3660) were downregulated and a gene coding for a phosphoketolase (lp_2659), involved in boosting the carbon flow from the pentose phosphate pathway to glycolytic intermediates, was also underexpressed. Furthermore, lp_3015 encoding for a lytic transglycosidase that hydrolyzes the glycan strands of the cell wall peptidoglycan was downregulated (3.4-fold).
Microarray data showed the downregulation of glnR (lp_1580), which codes for a major transcription factor that acts as the global repressor of nitrogen metabolism in L. plantarum . Downregulation of glnR was accompanied by the concurrent underexpression of genes encoding functions involved in the assimilation and reallocation of nitrogen within the cell. More particularly, genes coding for functions involved in influx and production of glutamine (Gln) and ammonium, i.e. glnA (glutamine synthetase), lp_0822 (glutamine-fructose-6-phosphate transaminase), genes coding for a Gln ABC transporter (lp_0802 and lp_0803), lp_2830 (aspartate ammonia-lyase which produces fumarate from Asp with the production of ammonium), amtB (NH4 + transport protein), were all downregulated.
Among the genes involved in nitrogen metabolism, only glnH (lp_2312) was upregulated. This gene codes for the substrate binding GlnH subunit of a L. plantarum putative glnPQH ABC transport system. According to Interproscan (http://www.ebi.ac.uk/InterProScan), the putative GlnH protein is predicted to contain a domain with similarity to an ionotropic glutamate receptor, suggesting that this ABC transport system could transport glutamate, which is essential for L. plantarum growth.
GA represses the expression of predicted ABC-type transporters
The adaptive response to GA in L. plantarum included the downregulation (2 to 7.2-fold at 15 mM GA) of six clustered genes coding for components of two ABC-type transport systems (lp_2739-lp_2740 and lp_2743-lp_2744), a putative transcriptional regulator (lp_2742) and a small peptide (lp_2741). The lp_2739 and lp_2743 genes encode the putative ATPases of their respective ABC systems, whereas lp_2740 and lp_2744 display transmembrane domains and encode the putative permeases. The ABC transporter encoded by lp_2739-lp_2740 is homologous to the BceAB transport system from B. subtilis, which functions as a detoxification pump of peptide antibiotics and antibiotics unrelated to antimicrobial peptides .
The predicted lp_2742 product (123 aa) displays the signatures of YtrA subfamily of the GntR superfamily regulators: a reduced C-terminal domain (which tends to be shorter in this subfamily) with only two α-helices . The regulators of the YtrA subfamily often respond to antibiotic stress  and typically modulate the expression of ABC transport systems . Accordingly, the predicted lp_2743 product is homologous to an ABC component of the YtrA operon from B. subtilis, which responds to cell wall antibiotics .
Redox and hydrolase related functions responsive to GA
We noticed the overexpression of two clustered genes, lp_1424 and lp_1425. The gene lp_1425 encodes for one of the six redundant copies of the fumarate reductase subunit A (FrdA) from L. plantarum WCFS1  and lp_1424 encodes for a putative NADPH-dependent FMN oxidoreductase, both encoding proteins predicted to be located in the cytoplasm . The best Phyre hits of the putative lp_1424 product were found with members of the NADH_dh2 family of putative flavin-binding quinone reductases. Fittingly, the predicted lp_1424 product contains the signature sequence LPVTPEYNXXXXXXLKNAIDXXS typical from the members of this family.
Interestingly, the presence of GA induced a marked decrease (4.1-fold) in lp_2229 transcripts encoding for a protein with homology to hydrolases of the metallo β-lactamase (MBL) superfamily. The predicted lp_2229 product displays the HXHXDH motif (HYHHDH; residues 59–64) absolutely conserved among all the MBLs  and presented the best Phyre hits with Zn-dependent MBL-like hydrolases which use tRNA and mRNA substrates. One of these MBL families, the UlaG family of proteins , is functionally specialized to process substrates distinct from nucleic acids.
Some genes involved in the response of L. plantarum to GA clustered together or were grouped in operon structures, suggesting co-regulatory relations. Downregulation of glnR (lp_1580) coding for GlnR, the master regulator of nitrogen metabolism in L. plantarum, was concurrent with the downregulation of glnA (lp_1581), ansB (lp_2830), amtB (lp_0349) and a glnPHQ glutamine ABC transport system (lp_0802 and lp_0803) fitting with the previously predicted GlnR regulon architecture .
Besides GlnR, another GA-responsive regulons of L. plantarum were the regulon Lp_2742 (http://regprecise.lbl.gov/RegPrecise/regulon.jsp) encompassing genes lp_2739 to lp_2744, which is predicted to be controlled by the transcription factor encoded by lp_2742 (GntR family, YtrA subfamily) and the RbsR regulon which encompass rbsK (lp_3660) and rbsD (lp_3659) genes.
Here, we have explored the mechanisms of adaptation of L. plantarum to GA. Gene expression-wide analyses did not reveal stringent or general stress responses. The most striking response to this model hydroxybenzoic acid, the inducible decarboxylation of GA to PG, was revealed by the combination of transcriptome and physiological analyses.
The marked downregulation of genes coding for certain ABC transporters, and its homology-based function prediction related to the response to cell wall drugs, suggests that they could also be involved in the traffic of GA and PG across the membrane, especially as it has been reported that these transporters are usually stress-responsive proteins extruding a wide range of drugs . A gene (lp_2229) coding for a putative MBL-like hydrolase homologous to the UlaG family, which uses tRNA and mRNA substrates, was notably downregulated. The UlaG family is poorly characterized, but it has been suggested that the ancestral RNase fold of this protein had been subjected to extensive modifications in order to acquire new catalytic activities which appear to be very diverse in bacteria . For example, UlaG is a MBL-like hydrolase that has been functionally specialized to act in the cytoplasm as l-ascorbate 6-phosphate lactonase involved in l-ascorbate metabolism . The emergence of these novel metabolic functions, coupled with the fact that LP_2229 is predicted to be located in the cytoplasm (according to LocateP), could provide a rationale for the functional specialization of the lp_2229 product to hydrolyse the aromatic ring of GA, although this remains to be proven. If these were the functions of the MBL-like hydrolase and the downregulated ABC-transporters, the observed adaptation (downregulation) could be directed to preserve GA availability and integrity in order to foster the PMF-generating effects linked to the GA metabolism mentioned above.
The induction of genes coding for redox related functions, such as lp_1424 and lp_1425, would be expected since GA is a quinol and it tends to autoxidize and/or undergo vicarious reduction by cellular one-electron reducers  thus generating highly oxygen reactive semiquinones. According to its intracellular location and predicted functions, the GA-mediated induction of the gene pair lp_1424/lp_1425 could be targeted to activating one-step reduction of intracellular quinols to its corresponding quinone (lp_1424) and vice versa (lp_1425). If active, this redox cycle would preempt the autoxidation and/or vicarious reduction of quinols that, such as p-coumaric and gallic acids, have been transported into the cytoplasm. In support of this hypothesis, the lp_1424/lp_1425 gene pair was also previously shown to be strongly induced by p-coumaric acid stress .
The transcriptional analysis also identified the differential expression of genes involved in nitrogen metabolism. Beyond the response observed at 1.5 mM GA, which only affected the amtB gene, transcriptional changes related to nitrogen metabolism were expanded to genes of the GlnR regulon at a higher GA concentration (15 mM). The concurrent downregulation of glnR (coding for the GlnR repressor) and the genes encompassing the GlnR regulon indicates that GlnR-mediated regulation in L. plantarum occurs in a similar manner to that of Bacillus subtilis . Supporting this hypothesis the L. plantarum GlnR and glutamine synthetase (GlnA) encoding genes are organized, like in B. subtilis, as glnRA operon. In B. subtilis GlnR represses the transcription of the glnRA operon (negative autoregulation) . In this bacterium the presence of the feed-back inhibited glutamine synthase promotes GlnR dimerization (activation) and increases the affinity of GlnR for the binding site thus acting as a chaperone which stabilizes the GlnR-DNA complexes . The downregulation of genes encompassing the GlnR regulon, fits with the recently proposed role of GlnR: to prevent the influx, and simultaneously limit the intracellular production, of ammonia and glutamine . In addition, the downregulation of genes involved in ribose transport (lp_3660 (transporter)) and utilization (lp_3659 (kinase)) is probably related to the tight control of ammonia levels, as ribose catabolism is closely linked to nucleoside metabolism. Thus, purine biosynthesis from ribose 5-P contributes to regenerate 2-ketoglutarate which is interconverted with glutamine and glutamate in central metabolic conversions to assimilate and re-distribute nitrogen within the cell . This tight control of ammonia levels avoids diffusion of ammonia out of the cell and ensures sufficient biosynthetic nitrogen assimilation flux required for growth . The fact that GA act as a signal to accumulate nitrogen could be related with the low nutrient availability encountered in plants, the natural habitat of L. plantarum. In fact, ammonium efflux is limited in L. plantarum under conditions of extremely low growth mimicking that encountered in plant environments .
The insights on the molecular adaptive responses of L. plantarum to GA revealed by the transcriptional datasets identified a marked induction of functions dedicated to metabolize tannic acid and GA, two phenolic compounds with antimicrobial activity and potent inhibitory effects towards cellulases . This knowledge may be used to guide the design of improved detoxification pretreatments of cellulosic hydrolysates by selection/engineering of production hosts well suited to remove these fermentation inhibitors.
Interestingly, genes encoding functions involved in the adaptation of L. plantarum to the GI-tract conditions, such as lp_2940 and lp_2956 (tannase), were highly overexpressed by GA. The crucial role played by lp_2940 in the persistence and survival of L. plantarum WCFS1 in the GI-tract has been already reported , albeit the molecular mechanism involved in this function remains to be elucidated. Regarding tannase, it may confer advantages to microbes to persist in the GI-tract as tannins display antimicrobial activity  and reduce the species richness in the gut . Other responses to GA related to an increased survival of L. plantarum in the GI-tract include the induction of genes involved in GlcNAc utilization, which is part of the human intestinal mucus glycoproteins and can be used as carbon source by Lactobacillus, particularly under bile stress . Altogether, these expression profiles indicate that contact with GA could improve the persistence of L. plantarum in the GI-tract. This hypothesis is in line with the accumulating evidence that the response of L. plantarum to some dietary polyphenol provides crossprotection against gastrointestinal stress [11–13]. Further, the induction by GA of genes required for tannin and GA metabolism provides valuable information to develop strategies aimed at increasing the bioactivity and health effects of these metabolites. For instance, tannase has been recently observed as an underlying function of the microbiota that selectively colonizes colorectal cancer tissues . Considering the anti-carcinogenic roles displayed by GA and PG (end products of tannin metabolism), successful strategies could be designed to kill colorectal cancer cells, more in view of the strong GA-mediated up-regulation of genes required for tannin and GA metabolism observed in this study.
This report expands our insight into the molecular mechanisms underlying the response of Lactobacillus to hydroxybenzoic acids. The transcriptome signature of GA response coincided with physiological analyses to show that the GA catabolic pathway required induction by the substrate (GA), and that the energetic consequence of GA metabolism in L. plantarum was the generation of a proton motive force (PMF). Expression profiling at genome scale provided insights on how GA operates globally suggesting that it acts as a signal to trigger versatile responses including nitrogen storage, tannin metabolism or GlcNAc utilization. The ability of L. plantarum to synthesize energy from GA in the form of PMF and the GA-mediated induction of nitrogen storage, utilization of carbon sources such as GlcNAc (which is used in the GI-tract under bile stress) or tannins (which exerts antimicrobial properties on intestinal bacteria), potentially confer competitive advantages to this microorganism in the GI-tract where nutrients are not in constant supply. This view is supported by the GA-mediated induction of genes playing a crucial role in the survival in the GI-tract. The strong GA-mediated induction of tannase and the enzymes involved in the GA catabolic pathway, highlight the ability of some lactic acid bacteria to better tolerate phenolic inhibitors that are present in pretreated lignocellulosic feedstocks, and further support the suggested potential of these microorganisms to become cell factories for the efficient conversion of cellulosic biomass .
Bacterial strain and culture conditions
Lactobacillus plantarum WCFS1 kindly provided by Dr. Michiel Kleerebezem (NIZO Food Research, The Netherlands) was grown in Man-Rogosa-Sharpe (MRS) broth (Difco Laboratories, Madrid, Spain)  at 30 °C without shaking. This strain is a colony isolate of L. plantarum NCIMB 8826, which was isolated from human saliva. It survives the passage through the human stomach  and persists in the digestive tract of mice and humans better than other Lactobacillus spp. isolated from the human intestine . GA (Sigma) stock solution (200 mM) was prepared in MRS. Appropriate dilutions were used to adjust GA final concentration to 1.5 or 15 mM in MRS.
Twelve paired independent L. plantarum WCFS1 batch cultures (50 mL each) were grown in MRS lacking GA to an OD660 ≈ 0.8–0.9. Then a culture of each pair (twelve cultures) was induced with GA to bring their final concentration to 1.5 mM or 15 mM (twelve biological replicates per each GA concentration). The induced cells and their corresponding controls were centrifuged at 4 °C after 10 min of exposure to this phenolic acid. The pellet was mixed with 2 mL of quenching buffer (60 % methanol, 66.7 mM HEPES, pH 6.5, −40 °C). Following quenching, the cells were centrifuged at 9000×g for 10 min at −10 °C and suspended in an extraction mixture (500 μL 1:4 chloroform-acid phenol, 30 μL of 10 % SDS, 30 μL Na-acetate 3 M pH 5.2, 400 μL Tris–EDTA buffer [10 mM Tris(hydroxymethyl)amino methane, 1 mM EDTA] pH 7.4, 15 mg of polyvinylpoly-pyrrolidone, and 500 mg of glass beads (ϕ, 75–150 μm). The cells were broken under frozen conditions in a FastPrep™ Fp120 (SAVANT) using three treatments of 5000 rpm for 40 s and chilled 1 min between cycles. The suspension was then centrifuged at 4 °C at 10,000×g for 2 min. After two extractions with 500 μL of chloroform the supernatant containing the RNA was immediately frozen in liquid nitrogen, and stored at −80 °C . NanoDrop ND1000 instrument was used for RNA quantification. The A260/A280 and A260/A230 ratios were measured to check RNA purity. Integrity and quality of RNA samples were determined by electrophoresis on agarose gels. Two treatments with DNase I (Ambion) were applied and the absence of genomic DNA was confirmed by PCR .
Microarray: cDNA synthesis, purification and hybridization
Before first-strand cDNA synthesis, RNA integrity was evaluated using the Agilent 2100 Bioanalyzer (Agilent Technologies). Fluorescently labelled cDNA was obtained by using the SuperScript Indirect cDNA Labeling System (Invitrogen). After, the Cy3 and HyPer5 fluorescent dyes (Amersham Biosciences) were coupled to the aminoallyl-modified first-strand cDNA, and purification of probes was carried out with the CyScribe GFX Purification Kit. Labeling efficiency was assessed using a NanoDrop ND1000 spectrophotometer. Preparation of probes and hybridization at 65 °C during 17 h was performed as described on the Two-Color Microarray Based Gene Expression Analysis Manual (Quick Amp Labeling) with Tecan HS Pro Hybridization (V. 5.7/Agilent Technologies). Slide L. plantarum WCFS1 8x15 K microarray GE Agilent G2509F Oligo Microarrays (No. 026636) was custom designed and contains 60-mer probes that were taken at the gene expression omnibus database (GEO Accession No.GPL5874). The oligo-microarray contained an average of three probes per transcript.
Real-time quantitative RT-PCR assays (qRT-PCR)
Real-time qRT-PCR was used to validate the microarray data. Amplification was carried out using a 7500 Fast System (Applied Biosystems). RNA was reverse transcribed using High Capacity cDNA Reverse Transcription Kits (Applied Biosystems). The specific primers used for the qRT-PCR assays are listed (see Additional file 2: Table S2). The SYBR Green method was used and each assay was performed in triplicate using SYBR Green real-time PCR Master Mix (Applied Biosystems). Amplification was initiated at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. Control PCRs were included to confirm the absence of primer dimer formation (no-template control), and to verify that there was no DNA contamination (without RT enzyme negative control). All real-time PCR assays amplified a single product as determined by melting curve analysis and by electrophoresis. A standard curve was plotted with cycle threshold (Ct) values obtained from amplification of known quantities of cDNAs and used to determine the efficiency (E) as E = 10−1/slope. The expression levels of target genes were normalized. The Bestkeeper analysis  was applied, and the geometric mean of the most stably expressed housekeeping genes (16S rRNA, ldhD, gapB, dnaG, and gyrA) was used as a normalization factor. The expression ratios measured by microarrays and by qRT-PCR assay were plotted, and the linear correlation coefficient was calculated (y = 0.973x + 0.107; R2 = 0.97) (see Additional file 3: Table S3).
Images were captured with a GenePix 4000B (Axon) and spots quantified using GenPix software (Axon). Background correction and normalization of expression data were performed using the methods normexp and loess in LIMMA, respectively . The expected False Discovery Rate (FDR) was controlled to be less than 5 %. Genes were considered differentially expressed when nominal p values were <0.05 and had a fold change (FC) equal or higher than ±1.5. FC was calculated as the average of the fold change between significantly regulated probes. Hybridizations and statistical analysis were performed by the Genomics Facility at Centro Nacional de Biotecnología, CSIC, Spain.
Microarray data accession number
The microarray data provided in this study have been deposited in NCBI Gene Expression Omnibus  genomics data repository and are accessible through GEO Series accession numbers GSE56997 and GSE63728.
Pyrogallol production by resting cells
To prepare resting cells, cultures of L. plantarum WCFS1 were grown in MRS medium at 30 °C to an OD660 ≈ 0.3. Cultures were then supplemented with 1.5 mM GA (final concentration) and cultures devoid of GA were used as controls. The GA-induced cultures and their non-induced controls continued to grow at 30 °C to mid-exponential phase (OD660 ≈ 0.6) and then centrifuged at 4 °C for 10 min at 3000 rpm. The cells were washed twice with 50 mM potassium phosphate buffer pH 5.8 and resuspended in the same buffer at 4 °C.
The assay for pyrogallol (PG) production was carried out with 1.5 ml of resting cells resuspended to an OD660 of 1.5 in 50 mM potassium phosphate buffer pH 5.8. Prior to assays the resting cells were incubated at 30 °C for 10 min. The assay was initiated at time zero with the addition of 1.5 mM GA. Aliquots of 250 µl were removed at fixed time intervals (every 15 or 30 min) and immediately centrifuged for 0.5 min at 14,000 rpm. The supernatant was stored on ice until further analysis by high-performance liquid chromatography (HPLC).
Phenolic compounds were extracted from the supernatants by a standard protocol, involving two extraction steps with one-third of the reaction mixture volume of ethyl acetate. The samples were then analyzed by high-pressure liquid chromatography (HPLC) with a diode array detector. A Thermo chromatograph (Thermo Electron Corporation, Waltham, MA) equipped with a P400 SpectraSystem pump, an AS3000 autosampler, and a UV6000LP photodiode array detector was used. A gradient of solvent A (water and acetic acid, 98:2, vol/vol) and solvent B (water, acetonitrile, and acetic acid, 78:20:2, vol/vol/vol) was applied to a reversed-phase Nova-pack C18 cartridge (25 cm by 4.0 mm [inner diameter]; particle size, 4.6 μm) at room temperature as follows: 0–55 min, 0–80 % solvent B, linear, 1.1 ml/min; 55–57 min, 80–90 % solvent B, linear, 1.2 ml/min; 57–70 min, 90–95 % solvent B, isocratic, 1.2 ml/min; 70–80 min, 95–100 % solvent B, linear, 1.2 ml/min; 80–90 min, 100 % linear, 1.2 ml/min; 100–120 min, washing with methanol 1.0 ml/min; and reequilibration of the column under initial gradient conditions. Detection was performed by scanning from 220 to 380 nm. Samples were injected onto the cartridge in duplicate, after being filtered through a 0.45-μm-pore-size polyvinylidene difluoride filter. The identification of phenolic compounds was carried out by comparing the retention times and spectral data of each peak with those of standards from commercial suppliers. GA and PG concentrations were estimated from the peak areas of the recorded chromatograms according to the “Integration Timed Events and Calibration Finnigan ChromQuestTM 4.2 Data System” software (CHROM-97202; Thermo Electron Corporation-2005).
Measurement of internal pH (∆pH) and membrane potential (∆Ψ)
The components of the proton motive force were measured as previously described . GA-induced and non-induced resting cells were prepared as described above (when the culture reached an OD660 ≈ 0.6) but resuspended at high densities, as outlined below. The internal pH was measured by loading resting cells of L. plantarum with the soluble fluorescent pH indicator 2′,7′-bis-(2-carboxyethyl)-5-(and 6)-carboxyfluorescein (BCECF) (Invitrogen Molecular Probes) as previously described . Briefly, 20 µL of resting cells suspended at high density (containing aprox. 50 mg/ml of protein) in 50 mM potassium phosphate (Kpi) buffer pH 5.8 were mixed with 1 µL of a 10 mM BCECF solution and 2.5 µL of 0.5 N HCl and the mixture was incubated 5 min at room temperature to load BCECF. Then the acid shock was stopped by addition of 1 mL of 50 mM Kpi buffer pH 5.8. The BCECF-loaded cells were washed 3 times with 50 mM Kpi buffer pH 5.8, resuspended in 200 µL of the same buffer and kept on ice until use. Fluorescence measurements were performed in 1-cm cuvettes containing 3 mL of 50 mM KPi pH 5.8 buffer equilibrated at 30 °C and 10 µL of BCECF-loaded cells. The mixture in the cuvette was stirred with a magnetic stirring bar. Fluorescence was measured using excitation and emission wavelengths of 502 and 525 nm with slit widths of 4 and 15 nm, respectively. The fluorescence signal was sampled every second in a Fluorolog-3 spectrofluorimeter (Jobin–Yvon-Spex). Opening of the measurement compartment caused loss of data during the first 5–6 s after an addition to the cuvette was made. The internal pH was calculated as described previously .
Membrane potential was measured qualitatively with the fluorescent probe 3,3′-dipropylthiocarbocyanine iodide (DiSC3) (Invitrogen Molecular Probes) . Increases in electrical potential across the membrane correlate with decreases in fluorescence intensity. A volume of 10 µL of resting cells (OD660 ≈ 6 in 50 mM KPi pH 5.8), was added to 3 mL of the same buffer. DiSC3 was added from a stock solution to a final concentration of 2 µM to the cuvette containing the cells in 50 mM KPi pH 5.8 buffer. The system was left to equilibrate for 10 min at 30 °C. Fluorescence measurements were performed using excitation and emission wavelengths of 500 and 705 nm, respectively, and slit widths of 8 nm.
IR, BR, RMu and FLF conceived and designed the experiments. IR, RMa and FLF performed the experiments. IR, RMa, BR, RMu and FLF analyzed the data. FLF drafted the manuscript. IR, BR, RMu helped to draft the manuscript. All authors read and approved the final manuscript.
This work was supported by grants AGL2011-22745 and AGL2014-52911 (MINECO).
The authors have declared no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Serrano J, Puupponen-Pimia R, Dauer A, Aura AM, Saura-Calixto F. Tannins: current knowledge of food sources, intake, bioavailability and biological effects. Mol Nutr Food Res. 2009;53:S310–29.View ArticleGoogle Scholar
- Kim SH, Jun CD, Suk K, Choi BJ, Lim H, Park S, Lee SH, Shin HY, Kim DK, Shin TY. Gallic acid inhibits histamine release and pro-inflammatory cytokine production in mast cells. Toxicol Sci. 2006;91:123–31.View ArticleGoogle Scholar
- Kang MS, Oh JS, Kang IC, Hong SJ, Choi CH. Inhibitory effect of methyl gallate and gallic acid on oral bacteria. J Microbiol. 2008;46:744–50.View ArticleGoogle Scholar
- Verma S, Singh A, Mishra A. Gallic acid: molecular rival of cancer. Environ Toxicol Phar. 2013;35:473–85.View ArticleGoogle Scholar
- Bolca S, Van de Wiele T, Possemiers S. Gut metabotypes govern health effects of dietary polyphenols. Curr Opin Biotechnol. 2013;24:220–5.View ArticleGoogle Scholar
- Scalbert A. Antimicrobial properties of tannins. Phytochemistry. 1991;30:3875–83.View ArticleGoogle Scholar
- Yang S-J, Kataeva I, Hamilton-Brehm SD, Engle NL, Tschaplinski TJ, Doeppke C, Davis M, Westpheling J, Adams MWW. Efficient degradation of lignocellulosic plant biomass, without pretreatment, by the thermophilic anaerobe “Anaerocellum thermophilum” DSM 6725. Appl Environ Microbiol. 2009;75:4762–9.View ArticleGoogle Scholar
- Ximenes E, Kim Y, Mosier N, Dien B, Ladisch M. Deactivation of cellulases by phenols. Enzyme Microb Technol. 2011;48:54–60.View ArticleGoogle Scholar
- Hervert-Hernandez D, Goñi I. Dietary polyphenols and human gut microbiota: a review. Food Rev Int. 2011;27:154–69.View ArticleGoogle Scholar
- Winkler J, Kao KC. Transcriptional analysis of Lactobacillus brevis to N-butanol and ferulic acid stress responses. PLoS One. 2011;6(8):e21438. doi:10.1371/journal.pone.0066473.View ArticleGoogle Scholar
- Reverón I, de las Rivas B, Muñoz R, López de Felipe F. Genome-wide transcriptomic responses of a human isolate of Lactobacillus plantarum exposed to p-coumaric acid stress. Mol Nutr Food Res. 2012;56:1848–59.View ArticleGoogle Scholar
- Curiel J, Rodríguez H, de las Rivas B, Anglade P, Baraige F, Zagorec M, Champomier-Vergès M, Muñoz R, López de Felipe F. Response of a Lactobacillus plantarum human isolate to tannic acid challenge assessed by proteomic analyses. Mol Nutr Food Res. 2011;55:1454–65.View ArticleGoogle Scholar
- Reverón I, Rodríguez H, Campos G, Curiel JA, Ascaso C, Carrascosa AV, Prieto A, de las Rivas B, Muñoz R, López de Felipe F. Tannic acid-dependent modulation of selected Lactobacillus plantarum traits linked to gastrointestinal survival. PLoS One. 2013;8:e66473. doi:10.1371/journal.pone.0066473.View ArticleGoogle Scholar
- Jiménez N, Curiel JA, Reverón I, de las Rivas B, Muñoz R. Uncovering the Lactobacillus plantarum WCFS1 gallate decarboxylase involved in tannin degradation. Appl Environ Microbiol. 2013;79:4253–63.View ArticleGoogle Scholar
- Licandro-Seraut H, Gury J, Tran NP, Barthelmebs L, Cavin JF. Kinetics and intensity of the expression of genes involved in the stress response tightly induced by phenolic acids in Lactobacillus plantarum. J Mol Microb Biotech. 2007;14:41–7.View ArticleGoogle Scholar
- Niho N, Shibutani M, Tamura T, Toyoda K, Uneyama C, Takahashi N, Hirose M. Subchronic toxicity study of gallic acid by oral administration in F344 rats. Food Chem Toxicol. 2001;39:1063–70.View ArticleGoogle Scholar
- Bron PA, Meijer M, Bongers RS, de Vos WM, Kleerebezem M. Dynamics of competitive population abundance of Lactobacillus plantarum ivi gene mutants in faecal samples after passage through the gastrointestinal tract of mice. J Appl Microbiol. 2007;103:1424–34.View ArticleGoogle Scholar
- Sobczak I, Lolkema JS. The 2-hydroxycarboxylate transporter family: physiology, structure, and mechanism. Microbiol Mol Biol Rev. 2005;69:665–95.View ArticleGoogle Scholar
- Zhou M, Boekhorst J, Francke Ch, Siezen R. LocateP: Genome-scale subcellular-location predictor for bacterial proteins. BMC Bioinform. 2008;9:173.View ArticleGoogle Scholar
- Kormelink TG, Koenders E, Hagemeijer Y, Overmars L, Siezen RJ, de Vos WM, Francke C. Comparative genome analysis of central nitrogen metabolism and its control by GlnR in the class Bacilli. BMC Genom. 2012;13:191–206.View ArticleGoogle Scholar
- Rietkötter E, Hoyer D, Mascher T. Bacitracin sensing in Bacillus subtilis. Mol Microbiol. 2008;68:768–85.View ArticleGoogle Scholar
- Rigali S, Derouaux A, Giannotta F, Dusart J. Subdivision of the helix-turn-helix GntR family of bacterial regulators in the FadR, HutC, MocR, and YtrA subfamilies. J Biol Chem. 2002;277:12507–15.View ArticleGoogle Scholar
- Salzberg LI, Luo Y, Hachmann AB, Mascher T, Helmann JD. The Bacillus subtilis GntR family repressor YtrA responds to cell wall antibiotics. J Bacteriol. 2011;14:5793–801.View ArticleGoogle Scholar
- Kleerebezem M, Boekhorst J, van Kranenburg R, Molenaar D, Kuipers OP, Leer R, Tarchini R, Peters SA, Sandbrink HM, Fiers MWEJ, Stiekema W, Lankhorst RMK, Bron PA, Hoffer SM, Groot MSN, Kerkhoven R, de Vries M, Ursing B, de Vos WM, Siezen RJ. Complete genome sequence of Lactobacillus plantarum WCFS1. Proc Natl Acad Sci USA. 2003;2003(100):1990–5.View ArticleGoogle Scholar
- Rasmussen BA, Bush K. Carbapenem-hydrolyzing beta-lactamases. Antimicrob Agents Chemother. 1997;41:223–32.Google Scholar
- Garces F, Fernández FJ, Montellà C, Penya-Soler E, Prohens R, Aguilar J, Baldomà L, Coll M, Badia J, Vega MC. Molecular architecture of the Mn2+-dependent lactonase UlaG reveals an RNase-like metallo-β-lactamase fold and a novel quaternary structure. J Mol Biol. 2010;398:715–29.View ArticleGoogle Scholar
- Davidson AL, Dassa E, Orelle C, Chen J. Structure, function, and evolution of bacterial ATP-binding cassette systems. Microbiol Mol Biol Rev. 2008;72:317–64.View ArticleGoogle Scholar
- Fernández FJ, Garces F, López-Estepa M, Aguilar J, Baldoma L, Coll M, Badia J, Vega MC. The UlaG protein family defines novel structural and functional motifs grafted on an ancient RNase fold. BMC Evol Biol. 2011;11:273.View ArticleGoogle Scholar
- Ackerley DF, Gonzalez CF, Park CH, Blake R, Keyhan M, Matin A. Chromate-reducing properties of soluble flavoproteins from Pseudomonas putida and Escherichia coli. Appl Environ Microbiol. 2004;70:873–82.View ArticleGoogle Scholar
- Schreier HJ, Brown SW, Hirschi KD, Nomellini JF, Sonenshein AL. Regulation of Bacillus subtilis glutamine synthetase gene expression by the product of the glnR gene. J Mol Biol. 1989;210:51–63.View ArticleGoogle Scholar
- Wray LV, Fisher SH. Bacillus subtilis GlnR contains an autoinhibitory C-terminal domain required for the interaction with glutamine synthetase. Mol Microbiol. 2008;68:277–85.View ArticleGoogle Scholar
- Goffin P, van de Bunt B, Giovane M, Leveau J, Höppener-Ogawa S, Teusink B, Hugenholtz J. Understanding the physiology of Lactobacillus plantarum at zero growth. Mol Syst Biol. 2010;6:413.View ArticleGoogle Scholar
- Boogerd FC, Ma H, Bruggeman FJ, van Heeswijk WC, García-Contreras R, Molenaar D, Krab K, Westerhoff HV. AmtB-mediated NH3 transport in prokaryotes must be active and as a consequence regulation of transport by GlnK is mandatory to limit futile cycling of NH4(+)/NH3. FEBS Lett. 2011;585:23–8.View ArticleGoogle Scholar
- Marco ML, De Vries MC, Wels M, Molenaar D, Mangell P, Ahrne S, De Vos WM, Vaughan EE, Kleerebezem M. Convergence in probiotic Lactobacillus gut-adaptive responses in humans and mice. ISME J. 2010;4:1481–4.View ArticleGoogle Scholar
- Chung KT, Lu Z, Chou MW. Mechanism of inhibition of tannic acid and related compounds on the growth of some intestinal bacteria. Food Chem Toxicol. 1998;36:1053–60.View ArticleGoogle Scholar
- Tan HY, Sieo CC, Abdullah N, Liang JB, Huang XD, Ho YW. Effect of condensed tannins on bovine rumen protist diversity based on 18S rRNA gene sequences. J Euk Microbiol. 2013;60:98–100.View ArticleGoogle Scholar
- Koskenniemi K, Laakso K, Koponen J, Kankainen M, Greco D, Auvinen P, Savijoki K, Nyman TA, Surakka A, Salusjärvi T, de Vos WM, Tynkkynen S, Kalkkinen N, Varmanen P. Proteomics and transcriptomics characterization of bile stress response in probiotic Lactobacillus rhamnosus GG. Mol Cell Proteom. 2011;10(2):M110-002741. doi:10.1074/mcp.M110.002741.View ArticleGoogle Scholar
- López de Felipe F, de las Rivas B, Muñoz R. Bioactive compounds produced by gut microbial tannase: implications for colorectal cancer development. Front Microbiol. 2014;5:684. doi:10.3389/fmicb.2014.00684.Google Scholar
- Boguta AM, Bringel F, Martinussen J, Jensen PR. Screening of lactic acid bacteria for their potential as microbial cell factories for bioconversion of lignocellulosic feedstocks. Microb Cell Fact. 2014;13:97.View ArticleGoogle Scholar
- de Man J, Rogosa M, Sharpe M. A medium for the cultivation of lactobacilli. J Appl Microbiol. 1960;23:130–5.Google Scholar
- Vesa T, Pochart P, Marteau P. Pharmacokinetics of Lactobacillus plantarum NCIMB 8826, Lactobacillus fermentum KLD, and Lactococcus lactis MG1363 in the human gastrointestinal tract. Aliment Pharmacol Ther. 2003;14:823–8.View ArticleGoogle Scholar
- Pavan S, Desreumaux P, Mercenier A. Use of mouse models to evaluate the persistence, safety, and immune modulation capacities of lactic acid bacteria. Clin Diagn Lab Inmunol. 2003;10:696–701.Google Scholar
- Saulnier DM, Molenaar D, de Vos WM, Gibson GR, Kolida S. Identification of prebiotic fructooligosaccharide metabolism in Lactobacillus plantarum WCFS1 through microarrays. Appl Environ Microbiol. 2007;73:1753–65.View ArticleGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26:509–15.View ArticleGoogle Scholar
- Smyth GK, Speed T. Normalization of cDNA microarray data. Methods. 2003;31:265–73.View ArticleGoogle Scholar
- Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.View ArticleGoogle Scholar
- Magni C, de Mendoza D, Konings WN, Lolkema JS. Mechanism of citrate metabolism in Lactococcus lactis: resistance against lactate toxicity at low pH. J Bacteriol. 1999;181:1451–7.Google Scholar
- Molenaar D, Abee T, Konings WN. Continuous measurement of the cytoplasmic pH in Lactococcus lactis with a fluorescence pH indicator. Biochim Biophys Acta. 1991;1115:75–83.View ArticleGoogle Scholar
- Šíp M, Heøman P, Plášek J, Hrouda V. Transmembrane potential measurement with carbocyanine dye diS-C3-(5): fast fluorescence decay studies. J Photochem Photobiol B. 1990;4:321–8.View ArticleGoogle Scholar