Chen Chen
Yanqing Lu
Linlin Wang
Haiyan Yu
Huaixiang Tian*- Department of Food Science and Technology, Shanghai Institute of Technology, Shanghai, China
Fructooligosaccharides (FOSs) metabolism in Lactobacillus plantarum is controlled by two gene clusters, and the global regulator catabolite control protein A (CcpA) may be involved in the regulation. To understand the mechanism, this study focused on the regulation relationships of CcpA toward target genes and the binding effects on the catabolite responsive element (cre). First, reverse transcription-PCR analysis of the transcriptional organization of the FOS-related gene clusters showed that they were organized in three independent polycistronic units. Diauxic growth, hierarchical utilization of carbohydrates and repression of FOS-related genes were observed in cultures containing FOS and glucose, suggesting carbon catabolite repression (CCR) control in FOS utilization. Knockout of ccpA gene eliminated these phenomena, indicating the principal role of this gene in CCR of FOS metabolism. Furthermore, six potential cre sites for CcpA binding were predicted in the regions of putative promoters of the two clusters. Direct binding was confirmed by electrophoretic mobility shift assays in vitro and chromatin immunoprecipitation in vivo. The results of the above studies suggest that CcpA is a vital regulator of FOS metabolism in L. plantarum and that CcpA-dependent CCR regulates FOS metabolism through the direct binding of CcpA toward the cre sites in the promoter regions of FOS-related clusters.
Introduction
Lactobacilli have complex nutritional requirements for fermentable carbohydrates and derive metabolic energy from homofermentative or heterofermentative carbohydrate fermentation (Gänzle and Follador, 2012). The broad ecological distribution and diverse habitats of lactobacilli reflect their metabolic flexibility and suggest that they are able to utilize a wide range of carbohydrates (Kant et al., 2011; Goh and Klaenhammer, 2015). Nevertheless, lactobacilli seldom use different carbon sources simultaneously. Rather, they organize their carbohydrate utilization in a hierarchical manner to achieve maximal growth (Zeng et al., 2017). The presence of a preferred carbon source prevents the utilization of the secondary substrate until the more favorable carbon source is exhausted. This hierarchical phenomenon, which was first described in relation to glucose-lactose diauxie in Escherichia coli, has been termed carbon catabolite repression (CCR) (Jankovic and Brückner, 2007; Görke and Stülke, 2008). CCR is important for competitive ability of lactobacilli in natural environments, such as the gastrointestinal tract (GIT). As the selection of the preferred carbon source is a major determining factor of growth rate for the microorganisms in GIT, and is also the competitive success compared with other microorganisms (Görke and Stülke, 2008; Becerra et al., 2015; Cardarelli et al., 2015).
The mechanism of CCR in Gram-positive bacteria involves negative regulation mediated by the catabolite control protein A (CcpA), a protein of the LacI-GalR family (Hueck and Hillen, 1995; Swint-Kruse and Matthews, 2009). It is a pleiotropic regulator involved in many important cellular processes. In the presence of a preferred carbon source, usually glucose, CcpA can effectively bind to its cognate operator sites, termed catabolite responsive elements (cre) with the assistance of seryl-phosphorylated histidine phosphocarrier protein (HPr-Ser∼P) (Muscariello et al., 2001; Deutscher et al., 2006; Wu et al., 2015). The consensus sequence of a cre site typically has a palindromic nucleotide motif, although some variations in nucleotide composition and length have been reported in recent literatures (Kim and Burne, 2017; Yang et al., 2017). The cre site is often located in the promoter region or within open reading frames of the regulated genes and operons (Marciniak et al., 2012). The binding of CcpA to cre sequences represses transcription of non-preferred metabolism genes (Fujita, 2009; Cai et al., 2012). CCR is relieved when decreased availability of the preferred carbohydrate results in decreased flux through glycolysis. CCR has been reported in several lactobacilli species, in which it is involved in the regulation of carbon catabolism, aerobic and anaerobic metabolism (Ana et al., 2011; Zotta et al., 2012), stress tolerance (Li et al., 2016), and metabolite production. During these regulatory processes, CcpA plays a central role that inactivation the ccpA gene leads to partial or complete relief of CCR (Moye et al., 2014; Zeng et al., 2017).
As one of the important prebiotics, fructooligosaccharides (FOSs) are non-digestible food ingredients that can stimulate the growth and activity of beneficial microbes residing in the GIT (Kaplan and Hutkins, 2000; Goh and Klaenhammer, 2015). Despite it have been widely applied in biotechnological applications and food industry, limited information is available about the regulation of FOS metabolism for specific members of the intestinal microbiota. Goh discovered a diauxic growth pattern when Lactobacillus paracasei 1195 was grown on media containing 1% FOS and 0.1% glucose, and predicted that the regulation is controlled by CCR via binding of CcpA to the cre site (Goh et al., 2007). Some cre-like sequences have also been identified in the operons for FOS utilization in L. acidophilus NCFM (Barrangou et al., 2003) and Bifidobacterium breve UCC2003 (Ryan et al., 2005). In our previous study, two gene clusters (sacPTS1 and sacPTS26) were found to participate in the metabolism of FOS in Lactobacillus plantarum ST-III. A number of potential cre sites for CcpA binding have been predicted in the regions of putative promoters of the two clusters (Chen et al., 2015).
To determine whether FOS metabolism in L. plantarum is regulated by CcpA-dependent CCR, the present study investigated the growth of L. plantarum and its mutant strains in CDM medium containing FOS and limited glucose. Growth profiles, sugar consumption, and gene expression were monitored. The binding of CcpA to the cre sites was verified using electrophoretic mobility shift assays (EMSA) in vitro and chromatin immunoprecipitation assay (ChIP) in vivo. Our findings reveal the hierarchical nature of carbohydrate utilization in L. plantarum, and highlight the crucial role of CcpA in controlling FOS utilization.
Materials and Methods
Bacterial Strains and Growth Conditions
The bacterial strains, plasmids and primers used in this study and their relevant features are listed in Table 1 and Supplementary Table 1. E. coli strain (DH5α and BL21), used as an intermediate cloning and expression host, was grown in Luria-Bertani (LB) medium at 37°C with aeration at 200 rpm/min. L. plantarum ST-III and its ccpA mutant strain were routinely cultured in de Man-Rogosa-Sharpe (MRS) broth (Merck, Darmstadt, Germany) at 37°C without agitation. To measure the growth of the wild-type and ccpA mutant strains, cells were grown in chemically defined medium (CDM) (Robert et al., 2000; Teusink et al., 2005) supplemented with FOS and glucose solutions (sterilized by passing through a 0.22 μm filter) as the basic carbon sources. The FOS used in this study was a commercial compound supplied by Meiji Seika Kaisha (Tokyo, Japan), comprising 9.8% (w/w) fructosyl-nystose, 37.3% (w/w) 1-kestose, 2.3% (w/w) sucrose, 49.1% (w/w) nystose, and 1.3% (w/w) glucose and fructose. Appropriate antibiotics were added to the media at the following concentrations: 100 μg/mL kanamycin, 50 μg/mL ampicillin, 30 μg/mL chloramphenicol and 250 μg/mL erythromycin for E. coli, 10 μg/mL chloramphenicol and 10 or 30 μg/mL (for replica plating) erythromycin for L. plantarum.
TABLE 1. Strains and plasmids used in this study.
Growth Analysis, Sampling, and Detection of Sugar Consumption Under Mixed Carbon Sources
Overnight cultures of L. plantarum ST-III or its ccpA mutant strain were transferred with 2% (v/v) inoculum into 500 mL of CDM containing 0.1% glucose and 0.4% FOS. The cultures were incubated without shaking for 16–18 h at 37°C. During the cells’ growth up to the stationary phase, the samples were withdrawn every 2 h to measure OD600 for growth analysis. To quantify sugars consumption and gene expression, culture samples were taken per hour (every 10 min in the diauxic growth period) during fermentation and centrifuged (12,000 × g, 4°C, 15 min). The supernatants were collected for residual sugar detection, and the cell pellets were used for RNA isolation.
For quantification of sugar consumption in the cell supernatants, high performance anion exchange chromatography (HPAEC) was carried out on an ICS5000 chromatograph (Dionex, Corp., Sunnyvale, CA, United States) with a CarboPac PA20 anion exchange column (Dionex, Corp., ID 3 mm × 150 mm) and a pulsed amperometric detector (PAD) as described previously (L’homme et al., 2001; Sánchez et al., 2004). After centrifugation (15,000 × g, 4°C, 15 min), the supernatant of each sample was filtered through a 0.45 μm nylon filter and then analyzed using a gradient elution procedure with H2O-250 mM NaOH as the mobile phase. The column was eluted at a flow rate of 0.5 mL/min and the injection volume of sample was 20 μL. Quantitative analyses were carried out using 5 mg/mL of solutions of the FOS standard set (Wako Pure Chemical Industries, Osaka, Japan), glucose, fructose, and sucrose (Sigma-Aldrich, Co., St. Louis, MO, United States).
RNA Extraction and RT-qPCR Analysis
Total RNA from L. plantarum ST-III cells collected from various conditions was extracted using TRIzol reagent (Invitrogen, Shanghai, China) according to the manufacturer’s instructions. RNA was subjected to RNase-free DNase I digestion and then purified using a PrimeScript RT reagent kit (Takara Bio, Dalian, China). The quantity and quality of total RNA were evaluated using a NanoDrop 2000C (Thermo, Waltham, MA, United States) and agarose gel electrophoresis. RNA preparations were stored at -80°C until use. Total RNA was reverse-transcribed to cDNA with a PrimeScript RT reagent kit (Takara Bio, Dalian, China) following the manufacturer’s instructions.
For RT-qPCR analysis, the generated cDNA was mixed with 0.2 μM gene specific primers (Supplementary Table 1) in a total volume of 25 μL. The PCR cycling conditions were as follows: 95°C for 10 min, followed by 40 cycles of amplification at 95°C for 15 s and 60°C for 30 s. All of the samples were measured in triplicate. Gene expression was normalized by the 2-ΔΔCT method and the 16S rRNA gene was used as the normalized standard.
Target Gene Structure Analysis and Prediction of Cre Sites
Gene structure analysis of sacPTS1 and sacPTS26 clusters (Chen et al., 2015) was performed on L. plantarum ST-III grown to early logarithmic phase in CDM medium supplemented with FOS. Total RNA was extracted and reverse-transcribed as described above and then cDNA products were generated as templates for PCR.
Primers (Supplementary Table 1) were designed on the basis of intergenic regions spanning the potential cotranscribed genes of the sacPTS1 and sacPTS26 clusters, respectively. As positive controls, each region was amplified with the same primers using chromosomal DNA of L. plantarum ST-III as the template. As negative controls, total RNA without the reverse transcription stage was used as the template.
Regulatory Sequence Analysis Tools (RSAT)1 was used to analyze the consensus motif of the cre sites. They were found by scanning all upstream regions in the genome of L. plantarum ST-III based on the profile of binding sites of CcpA in L. plantarum WCSF1 using the RegPrecise database2. A positional frequency matrix (PFM) was constructed for the collection of binding sites and those in the upstream regions of sacPTS1 and sacPTS26 clusters were searched. Scores of candidate sites were defined as the sum of positional nucleotide weights as previously described (Sánchez et al., 2004) and values greater than five were considered as the cre sites of binding of CcpA.
Construction of ccpA Mutant
The Cre-lox-based mutagenesis system was used for gene deletion of ccpA (Lambert et al., 2007). DNA fragments corresponding to the chromosomal regions upstream (1047 bp fragment; primer pairs UpF-UpR) and downstream (1001 bp fragment; primer pairs DownF-DownR) of the ccpA gene were amplified by PCR using a proofreading DNA polymerase (CWBIO, Shanghai, China) with chromosomal DNA used as the template. The amplicons were respectively cloned in the XhoI and Eco53kI restriction sites of suicide vector pNZ5319, and the recombinant mutagenesis vector, pNZ5319-up-down, was electroporated into L. plantarum ST-III cells. Candidate double-crossover clones, were selected and confirmed by PCR. Then, the lox66-P32-cat-lox71 cassette was resolved to a single double-mutant lox72 site through the transient plasmid pNZ5348. A double-crossover ccpA mutant without any resistance gene (ΔccpA strain) was acquired and confirmed by DNA sequencing as described previously (Capozzi et al., 2011; Chen et al., 2015).
Expression and Purification of His6-Tagged CcpA Protein
Expression and purification of His6-tagged CcpA was performed using the pET-28a(+) vector as previously described (Chen et al., 2014) with some modifications. Briefly, a 1,012 bp sequence of the ccpA gene was PCR amplified from L. plantarum genomic DNA using the primer pair ccpA-F and ccpA-R, which includes the HindIII and NheI sites (Supplementary Table 1). The recombinant plasmid pET-28a-ccpA was constructed as described previously (Chen et al., 2014) and the strain harboring this plasmid was named E. coil BL21-ccpA.
The His6-tagged CcpA protein expression was induced by addition of 1 mM isopropyl-β-D-thioIsopropyl-β-D-thio-galactoside (IPTG) when OD600 reached 0.4 to 0.6 at 37°C. The culture was allowed to grow at 25°C for 8 h. Cells were collected by centrifugation at 6000 × g for 10 min, and resuspended in binding buffer (0.2 M sodium phosphate, 5 M NaCl, 1 M imidazole, pH 7.4), then 1 mM PMSF and 0.1 mg/mL lysozyme solution were added, and the cells were then lysed by sonication. The cell debris was removed by centrifugation at 10,000 × g for 10 min at 4°C before RNase and DNase were added. The soluble protein was purified by nickel ion affinity chromatography using a Chelating Sepharose Fast Flow column (GE Healthcare, Waukesha, WI, United States) according to the manufacturer’s instructions. The purified protein was desalted and concentrated by Amicon Ultra-0.5 centrifugal filter devices (Millipore, Billerica, MA, United States).
Electrophoretic Mobility Shift Assay
Double-stranded 200-bp DNA fragments (named PsacK, Ppts1-sacA, Pagl4, and PsacR2) containing the putative cre sites (Supplementary Table 2), which were located in the four promoter regions of the sacPTS1 and sacPTS26 clusters, were amplified by PCR using specific primers, respectively (Supplementary Table 1). Then fluorescent FAM was added to these DNA fragments by PCR reaction with Dpx DNA polymerase (TOLO Biotech, Shanghai, China) using the universal primers M13F-47 (FAM) and M13R-48 (Supplementary Table 1). The FAM-labeled probes were purified using a Wizard® SV Gel and PCR Clean-Up System (Promega, United States) and were quantified with NanoDrop 2000C. A competition assay was performed by addition of 100-fold molar excess of unlabelled probe to the EMSA reaction mixture.
Electrophoretic mobility shift assay was performed in a 20 μL reaction volume that contains 40 ng probe and different amounts of purified His6-tagged CcpA proteins, in a reaction buffer of 50 mM Tris-HCl (pH 8.0), 100 mM KCl, 2.5 mM MgCl2, 0.2 mM DTT, 2000 ng polydIdC and 10% glycerol. After incubation for 30 min at 30°C, the reaction system was loaded onto 2% agarose gels buffered with 0.5× TBE. The gels were scanned with ImageQuant LAS 4000 mini (GE Healthcare) (Almengor et al., 2007). To verify the specific binding of CcpA to the cre sites, each putative cre site (generated from RSAT analysis according to the consensus motif) was mutated and named creMUT. The main principle of the mutation was as follows: the defined base in the consensus motif was mutated to the other three bases, and the “W” that represents A or T was mutated to “S,” which represents G or C. Specific primers were then synthesized for PCR amplification, and the amplicons were purified and self-ligated with an Ezmax One-Step Cloning Kit (TOLO Biotech, Shanghai, China) following the manufacturer’s instructions. The creMUT sequences were verified by DNA sequencing, and the probe preparation and EMSA analysis were the same as described above.
Chromatin Immunoprecipitation Assay
Inducible expression vectors carrying Flag-tagged CcpA in L. plantarum ST-III were constructed using the pSIP409 vector as previously described (Chen et al., 2014) with some modifications. Briefly, the ccpA gene sequence was amplified by PCR using the Flag-ccpA-409F (containing a Flag-tag on the 5′ terminus) and Flag-ccpA-409R primers (Supplementary Table 1). The purified PCR-products were ligated into the same restriction enzyme digested pSIP409 vector to construct the recombinant plasmid. Then, the constructed plasmid pSIP409-Flag-ccpA was transformed into L. plantarum ST-III. For the recombinant strain 409-Flag-ccpA, the cells were induced at an OD600 of 0.3 by adding the inducing peptide pheromone IP-673 (synthesized by Invitrogen, Shanghai, China) to a final concentration of 50 ng/mL. The cells were fixed in 1% formaldehyde and quenched by the addition of glycine. Then the cells were washed twice with 5 mM Tris-HCl (pH 8.0) containing 25 mg/mL lysozyme and sheared to an average fragment size of 300 to 500 bp by sonication. After centrifugation, the 5 μL supernatant was retained for use as the input sample, and the rest of the supernatant was incubated with Flag antibodies. Normal Rabbit IgG was set up as a negative control. Immunoprecipitated and input samples were de-crosslinked by incubation in 1× IP elution buffer and 5 M NaCl and 20 mg/mL proteinase K for 65°C at 1.5 h. The purified ChIP products and genomic input DNA were analyzed by qPCR using specific primers (Supplementary Table 1). The enrichment of DNA fragments was analyzed with the input DNA samples serving as controls. All of the samples were measured in triplicate.
Statistical Analysis
The data shown are representative of at least three independent experiments. Student’s t-test was used to determine statistical differences. Differences between samples with a p-value ≤ 0.05 were considered to be statistically significant.
Results
Analysis of Target Gene Structure and cre Consensus
Previous sequence analysis suggested that the FOS-related gene clusters of L. plantarum ST-III might be organized in polycistronic units (Chen et al., 2015). To verify this prediction, cDNA products generated from total RNA extracted from cells induced by FOS were used as PCR templates in RT-PCRs to test for the amplification of overlapping regions on each cluster (Vastano et al., 2016). The results of specific amplifications demonstrated that the sacPTS1 cluster is transcribed in two polycistronic units: sacK and pts1 are transcribed together and sacA, sacR, and agl2 are divergently oriented and cotranscribed as an operon. As expected, only a single PCR product was obtained from each amplification reaction for the sacPTS26 cluster, confirming it constitutes a single operon (Figure 1).
FIGURE 1. RT-PCR analysis of the transcriptional organization of sacPTS1 and sacPTS26 clusters in L. plantarum ST-III. M, DNA size marker; C+, positive control using chromosomal DNA of the strain as template; A, fragment amplified with the specific primer set using cDNA generated from FOS-induced L. plantarum ST-III cells as a PCR template; C–, negative control using total RNA without the reverse transcription step. Primers used are listed in Supplementary Table 1.
It has been shown that CcpA can directly bind to cre sites to regulate its target genes. However, the cre sequence differs among Gram-positive bacteria (Schumacher et al., 2011; Marciniak et al., 2012). This prompted us to search for a conserved common cre consensus motif within the genome of L. plantarum ST-III based on the RegPrecise database (Figure 2). Using the generated PFM to search in the sacPTS1 and sacPTS26 clusters, six potential cre sites were found. In addition to the five cre sites found in the previous study, a new cre site was discovered in the Ppts1-sacA region. As a result, this region has three putative cre sites (Figure 3). The scores of the six candidate sites were all greater than 5, suggesting the possibility of CcpA binding (Medina-Rivera et al., 2011).
FIGURE 2. Consensus sequence motif of cre sites in L. plantarum ST-III generated using RSAT software. A positional frequency matrix (PFM) was generated according to the occurrence frequency of each base at each location of the consensus sequence. The sequence-logo represents occurrence frequency and the height of each individual symbol reflects its prevalence at a given position.
FIGURE 3. Prediction of the potential cre sites in the sacPTS1 and sacPTS26 clusters of L. plantarum ST-III. Putative cre sites are underlined with red; red backgrounds show scores for each cre sites, defined as the sum of positional nucleotide weight. The presumed start codon of each gene is shown in capital letters, and the putative-10 and-35 promoter regions and possible ribosome-binding sites (RBSs) are marked.
Growth, Sugar Consumption, and Gene Expression of L. plantarum ST-III in Medium With Mixed Carbon Sources
The growth of L. plantarum ST-III in CDM containing both glucose and FOS (0.1 and 0.4%, respectively) was monitored (Figure 4A). A typical diauxic growth pattern with two distinct growth phases separated by a lag phase of about 2 h was observed. During the growth, residual sugars were detected in the supernatant at the same time. In the first 6 h, glucose and traces of fructose in the FOS mixture were consumed and the main components of FOS were kept constant. The diauxic lag is likely to have been caused by the depletion of glucose, as cessation of growth was observed at a similar time. After the diauxic lag phase, the cells resumed growing and entered a second growth phase using FOS as the carbon source. 1-Kestose was rapidly consumed and was undetectable after 10 h of fermentation. Nystose and fructosyl-nystose gradually decreased to an undetectable level at 14 and 16 h, respectively. To examine the kinetics of transcription of the FOS-related genes during the diauxic growth, RT-qPCR was performed to analyze the changes in the expression levels of these genes. RNA samples were collected at different sampling times and CT values obtained at 8 h of incubation were chosen as the control values to calculate the fold changes in gene expression over time. Each gene of the two clusters was analyzed using specific primer pairs (Supplementary Table 1). Similar CT values were obtained for group of sacK and pts1, group of sacA, sacR1, and agl2 and genes in the sacPTS26 clusters respectively (data not shown), which is consistent with the result of the RT-PCR experiments. Then sacK, sacA, and pts26 genes were chosen as the representatives for gene expression in the subsequent studies. As shown in Figure 4C, the expression levels of sacK were very low during the first growth phase when glucose was utilized as the preferred carbon source. On entering into the lag phase, gradually higher expression was observed. The maximum transcription levels were observed at 9 and 10 h when 1-kestose and nystose were actively hydrolyzed. The expression profiles of sacA and pts26 were similar. Expression was repressed during the first growth phase, then, after glucose depletion, the repression began to release and the transcription levels remained almost constant during the second growth phase. After 13 h, the transcription levels of all three genes were sharply decreased and this period was associated with the depletion of FOS in the medium. These observations are consistent with the results of previous transcriptome experiments in which glucose repressed the transcription of FOS-induced genes (Denich et al., 2003; Goh et al., 2006).
FIGURE 4. The growth curves, sugar utilization and gene expression of L. plantarum ST-III in CDM containing 0.1% glucose and 0.4% FOS. (A) The growth curves and sugar consumption of the wild-type strain during fermentation under mixed carbon sources. (B) The growth curves and sugar consumption of the ΔccpA strain during fermentation under mixed carbon sources. (C) The kinetics of transcription of the FOS-related genes during the diauxic growth of L. plantarum ST-III by RT-qPCR. CT values obtained at 8 h of incubation were chosen as the control.
ccpA Deletion Leads to Relief From CCR
To determine whether the inhibition of FOS utilization by glucose is regulated by catabolic repression via CcpA, a ΔccpA strain was constructed using the Cre-lox-based mutagenesis system (Lambert et al., 2007). The growth of the ΔccpA strain on glucose and FOS was compared with its parent strain. The ΔccpA strain exhibited a significantly reduced growth rate compared with the wild-type strain in glucose (a doubling time of 98.4 ± 1.0 min versus 76.2 ± 6 min) and FOS (a doubling time of 102.6 ± 1.5 min versus 81.9 ± 7.2 min). These results are in accordance with findings that CcpA inactivation significantly affected the growth rate of other microorganisms (Cai et al., 2012; Zeng et al., 2017). The deletion of the ccpA gene resulted in a complete loss of diauxic growth when the ΔccpA strain was grown in a mixture of glucose and FOS (Figure 4B). The components in the FOS mixture were used simultaneously with glucose during the whole growth process, although glucose was consumed at a faster speed. All of the sugars were depleted at the end of fermentation.
Gene expression of the FOS-related clusters in the wild-type and ΔccpA strains was also measured by RT-qPCR in the early logarithmic phase of the growth. The expression levels of the sacK, sacA, and pts26 genes in the wild-type and ΔccpA strains were measured with glucose, FOS and the mixture of glucose and FOS as the carbon sources respectively. The CT values of the wild-type strain grown on glucose were chosen as the control for the three genes to calculate the fold changes between the different conditions. As expected, all three genes were significantly down-regulated for the wild-type strains grown on glucose and in the mixture of glucose and FOS, compared with the wild-type grown on FOS (Table 2). In contrast, after ccpA inactivation, repression of these genes by glucose was drastically relieved in the ΔccpA strain compared with the wild-type. For example, the expression levels of the sacK, sacA, and sacPTS26 operons in the ΔccpA strain were respectively 2.9-, 3-, and 3.9-fold higher than the corresponding values for the wild-type strain when the cells were grown in a mixture of glucose and FOS. These results confirm the dominant role of CcpA in the CCR of FOS metabolism in L. plantarum ST-III. Notably, in the ΔccpA strain, sacA and pts26 showed 1.5- to 1.8-fold increases in the presence of FOS (FOS alone or mixed with glucose) versus glucose alone. This result implies that in addition to the CcpA-dependent CCR, these genes may be induced or derepressed by FOS.
TABLE 2. Relative transcript abundances of FOS-related genes in the wild-type and ΔccpA strains grown in different sugarsa.
CcpA Binds to cre Motif
To determine whether CcpA can bind to the six putative cre sites in vitro, His6-tagged CcpA was expressed heterologously in E. coli, and the purified protein was used to perform an EMSA (Almengor et al., 2007; Tiffert et al., 2008). As shown in Figure 5A (lanes 1 to 4), with increasing amounts of His6-tagged CcpA (0 to 1.2 μg), the intensities of the bands representing the shifted CcpA-DNA complex strengthened. In contrast, when labeled and unlabelled probes were used for a specific competitive assay (lane 5), no shift was detected for the labeled probe, demonstrating the binding specificity of CcpA to these DNA fragments.
FIGURE 5. Characterization and verification of the CcpA-DNA binding of the four promoter regions by EMSA. (A) EMSA of His6-tagged CcpA with the DNA fragments of four promoter regions carrying intact cre sites. (B) EMSA of His6-tagged CcpA with the DNA fragments of four promoter regions carrying mutated cre sites. The positions of CcpA-DNA complex (C) or free DNA (F) are indicated to the left of the figure.
Next, to verify that CcpA combined directly with the cre sites within the regions, EMSA experiments were performed to examine whether mutation of cre (Supplementary Table 2) affected the interaction between His6-tagged CcpA and its DNA fragments (Tiffert et al., 2008; Wu et al., 2015). Then DNA fragments of the four promoter regions, which contained the creMUT sites, were generated by PCR and used in EMSA (Willenborg et al., 2014). The results showed that the binding of His6-tagged CcpA protein to the new PsacK, Pagl4, and PsacR2 regions was completely abolished (Figure 5B), indicating that the CcpA did not bind to the three promoter regions that included creMUT. As the Ppts1-sacA region exists at three putative cre sites, the region was mutated three times and EMSA was conducted. The binding affinities of His6-tagged CcpA protein to the Ppts1-sacA region were weakened after one (cresacA1) or two sites (cresacA1 and cresacA1-sacpts1) were mutated (data not shown). When the Ppts1-sacA region was mutated at third time (cresacA2), the binding affinity of His6-tagged CcpA protein to the Ppts1-sacA region vanished. These results indicate that all three cre sites can bind to the CcpA protein.
Chromatin immunoprecipitation was then performed to test for the binding of CcpA to the promoter regions in vivo. The CcpA protein was Flag-tagged at its N terminus, and then the successful expression of 409-Flag-ccpA in L. plantarum was confirmed via western blot analysis (Figure 6A). The cross-linked DNA fragments were analyzed by RT-qPCR. As shown in Figure 6B, it is clear that the four regions were remarkably enriched by CcpA protein compared to the negative control sample, suggesting that CcpA interacts specifically with those promoter regions in vivo. Together, these findings suggest that CcpA protein can specifically bind to the four promoter regions, and six cre sites exist in these regions with high-affinity binding for CcpA.
FIGURE 6. Chromatin immunoprecipitation (ChIP) asssy for analysis of binding of CcpA with the DNA fragment of the four promoter regions. (A) Detection of Flag-tagged CcpA protein by western blotting assay with Flag antibodies. (B) Enrichment of Flag-tagged CcpA at the promoter regions was determined by ChIP-qPCR. Data presented are mean values based on at least three replicates. Error bars indicate standard deviations. Values that implied statistically significant differences (p ≤ 0.05) from negative control (Normal Rabbit IgG) are indicated with asterisks.
Discussion
Lactobacillus plantarum is a normal inhabitant of the human GIT and previous reports have shown that it can utilize FOS as efficiently as glucose, with two gene clusters participating in FOS metabolism (Chen et al., 2015). Due to the complex interspecies competition in the GIT, the regulation of carbohydrate metabolism is crucial to the survival and colonization of L. plantarum. CCR is one of the most important regulatory mechanisms used by both Gram-positive and Gram-negative bacteria to achieve maximum growth in mixed carbohydrate environments (Marasco et al., 2002; Almengor et al., 2007; Zeng et al., 2017). In the present work, the FOS metabolism of L. plantarum ST-III regulated by CcpA-dependent CCR was studied in vivo and in vitro.
When L. plantarum ST-III was grown in media that contained glucose in addition to FOS, diauxic growth was observed (Cai et al., 2012). Diauxic growth is common among bacteria and higher organisms because most of them preferentially use the carbon sources that are most easily accessible and allow the fastest growth. The presence of preferred carbon sources prevents the expression of catabolic systems that enable the use of secondary substrates (Görke and Stülke, 2008). When the preferred carbon source is exhausted, bacteria first need to synthesize the related enzymes (during the lag phase) before they can resume growth (Kearns and Russell, 1996). Goh also reported this phenomenon when cells were grown on FOS together with limited glucose (Yong, 2005). However, this is disadvantageous for completely fermenting all available sugars and efficiently transforming biomass to high value-added products for microbes in complex carbohydrate environments (Chu, 2015).
The diauxic growth and repression of FOS-related genes were observed in cultures containing glucose and FOS, which suggests that FOS utilization in L. plantarum is controlled by CCR. CcpA is believed to play a key role in CCR regulation in Gram-positive bacteria. As expected, in-frame deletion of ccpA resulted in the complete loss of diauxic growth, simultaneous use of both sugars and relief of the repressed gene expression. These phenomena are in accordance with previous reports that mutation of the ccpA gene leads to complete or partial release from CCR (Marasco et al., 2002; Fujita, 2009). These results show that CcpA is crucial in the CCR of FOS metabolism in L. plantarum.
Regulation of the transcription of CcpA-regulated genes involves the binding of CcpA to cre sites of the target genes. The consensus sequence of cre has been determined as a typical 14- to 16-nucleotides sequence, such as “TGWAANCGNTNWCA” (where W is A or T, N is any base) and “WTGNNARCGNWWWCAW” in Bacillus subtilis (Yang et al., 2017) and L. plantarum (Wels et al., 2011), respectively. However, it has recently been found that CcpA has flexible binding site architecture that is highly variable in both length and base composition (Yang et al., 2017). In this study, six putative cre sites were found in the promoter regions of FOS-related clusters based on the consensus motif generated from RSAT analysis. These sites deviate by one or two nucleotides from published cre consensus sequences (Hueck et al., 1994; Miwa et al., 2000; Jankovic and Brückner, 2007), but the binding was verified by EMSA in vivo and ChIP in vitro. These results provide a new insight into the structure of CcpA recognition sites in Gram-positive bacteria. In addition, the repressor role played by CcpA observed in this study is consistent with the location of cre sites within the promoter region. According to Zomer et al. (2007), the repression will occur at cre sites located in or downstream of the putative-35 and -10 sequences).
In the absence of a functional CcpA, other regulatory processes could be studied when CcpA-mediated repression is abolished. In our study, FOS in the medium resulted in activation of the FOS-related genes in the ΔccpA strain, demonstrating that other regulators participate in FOS metabolism of L. plantarum. These results confirm our previous prediction that L. plantarum have a double effect of global and local regulation of FOS metabolism. The local regulators SacR1 and SacR2 may also be involved in metabolic regulation in L. plantarum (Chen et al., 2015). In these two clusters, the local regulators SacR1 and SacR2, which are also members of LacI-GalR family, are cotranscribed with other FOS-related genes. In the presence of FOS, SacR1 and SacR2 may be induced by their substrates and maintain their own expression at a certain level to help bacteria to adjust sugar utilization to their metabolic capacities (Robert et al., 2008; Teixeira et al., 2013). Although the effects of CcpA have been confirmed, the interactions between specific local operators and the putative binding sites in these clusters are not directly proven and the exact binding site is not yet clear. To answer these questions, related experiments are currently being carried out.
Conclusion
The diauxic growth, hierarchical use of sugars and repression of FOS-related genes in L. plantarum grown on FOS in the presence of limiting glucose demonstrated that FOS utilization is subject to CCR. Inactivation of the ccpA gene eliminated these phenomena, proving the dominant role of CcpA in CCR for FOS metabolism. The binding of CcpA to the cre sites in the promoter regions was verified both in vivo and in vitro, which suggests that CcpA regulates CCR through direct regulation of the transcription of FOS-related clusters. As CCR is part of the regulatory network, further analysis of CCR may lead to deeper insights into a complex regulatory network. Our evolving understanding of the mechanistic interactions of probiotics and prebiotics will provide the molecular basis for the design of effective probiotic–prebiotic combinations to maximize host benefits.
Author Contributions
CC wrote the manuscript and the statistical analysis. YL analyzed the growth and expression of related genes in wild-type and mutated strains. LW executed the construction of the expression vector and ΔccpA strain. HY executed the target genes structure analysis and confirmed CcpA binding to the putative cre sites. HT designed the research.
Funding
This work was supported by the National Natural Science Foundation of China (Award No. 31501451).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01114/full#supplementary-material
Footnotes
References
Almengor, A. C., Kinkel, T. L., Day, S. J., and Mciver, K. S. (2007). The catabolite control protein CcpA binds to Pmga and influences expression of the virulence regulator Mga in the Group A Streptococcus. J. Bacteriol. 189, 8405–8416. doi: 10.1128/JB.01038-07
Ana, A., Isabelle, M. V., and Bruno, D. (2011). CcpA-mediated repression of Clostridium difficile toxin gene expression. Mol. Microbiol. 79, 882–899. doi: 10.1111/j.1365-2958.2010.07495.x
Barrangou, R., Altermann, E., Hutkins, R., Cano, R., and Klaenhammer, T. R. (2003). Functional and comparative genomic analyses of an operon involved in fructooligosaccharide utilization by Lactobacillus acidophilus. Proc. Natl. Acad. Sci. U.S.A. 100, 8957–8962. doi: 10.1073/pnas.1332765100
Becerra, J. E., Yebra, M. J., and Monedero, V. (2015). An l-fucose operon in the probiotic Lactobacillus rhamnosus GG is involved in adaptation to gastrointestinal conditions. Appl. Environ. Microbiol. 81, 3880–3888. doi: 10.1128/AEM.00260-15
Cai, J., Tong, H., Qi, F., and Dong, X. (2012). CcpA-dependent carbohydrate catabolite repression regulates galactose metabolism in Streptococcus oligofermentans. J. Bacteriol. 194, 3824–3832. doi: 10.1128/JB.00156-12
Capozzi, V., Weidmann, S., Fiocco, D., Rieu, A., Hols, P., Guzzo, J., et al. (2011). Inactivation of a small heat shock protein affects cell morphology and membrane fluidity in Lactobacillus plantarum WCFS1. Res. Microbiol. 162, 419–425. doi: 10.1016/j.resmic.2011.02.010
Cardarelli, H. R., Martinez, R. C., Albrecht, S., Schols, H., Franco, B. D., Saad, S. M., et al. (2015). In vitro fermentation of prebiotic carbohydrates by intestinal microbiota in the presence of Lactobacillus amylovorus DSM 16998. Benef. Microbes 7, 119–133. doi: 10.3920/BM2014.0151
Chen, C., Zhao, G., Chen, W., and Guo, B. (2015). Metabolism of fructooligosaccharides in Lactobacillus plantarum ST-III via differential gene transcription and alteration of cell membrane fluidity. Appl. Environ. Microbiol. 81, 7697–7707. doi: 10.1128/AEM.02426-15
Chen, C., Zhou, F., Ren, J., Ai, L., Dong, Y., Wu, Z., et al. (2014). Cloning, expression and functional validation of a β-fructofuranosidase from Lactobacillus plantarum. Process Biochem. 49, 758–767. doi: 10.1016/j.procbio.2014.02.013
Chu, D. F. (2015). In silico evolution of diauxic growth. BMC Evol. Biol. 15:211. doi: 10.1186/s12862-015-0492-0
Denich, T. J., Beaudette, L. A., Lee, H., and Trevors, J. T. (2003). Effect of selected environmental and physico-chemical factors on bacterial cytoplasmic membranes. J. Microbiol. Methods 52, 149–182. doi: 10.1016/S0167-7012(02)00155-0
Deutscher, J., Francke, C., and Postma, P. W. (2006). How phosphotransferase system-related protein phosphorylation regulates carbohydrate metabolism in bacteria. Microbiol. Mol. Biol. Rev. 70, 939–1031. doi: 10.1128/MMBR.00024-06
Fujita, Y. (2009). Carbon catabolite control of the metabolic network in Bacillus subtilis. Biosci. Biotechnol. Biochem. 73, 245–259. doi: 10.1271/bbb.80479
Gänzle, M. G., and Follador, R. (2012). Metabolism of oligosaccharides and starch in lactobacilli: a review. Front. Microbiol. 3:340. doi: 10.3389/fmicb.2012.00340
Goh, Y. J., and Klaenhammer, T. R. (2015). Genetic mechanisms of prebiotic oligosaccharide metabolism in probiotic microbes. Annu. Rev. Food Sci. Technol. 6, 137–156. doi: 10.1146/annurev-food-022814-015706
Goh, Y. J., Lee, J. H., and Hutkins, R. W. (2007). Functional analysis of the fructooligosaccharide utilization operon in Lactobacillus paracasei 1195. Appl. Environ. Microbiol. 73, 5716–5724. doi: 10.1128/AEM.00805-07
Goh, Y. J., Zhang, C., Benson, A. K., Schlegel, V., Lee, J. H., and Hutkins, R. W. (2006). Identification of a putative operon involved in fructooligosaccharide utilization by Lactobacillus paracasei. Appl. Environ. Microbiol. 72, 7518–7530. doi: 10.1128/AEM.00877-06
Görke, B., and Stülke, J. (2008). Carbon catabolite repression in bacteria: many ways to make the most out of nutrients. Nat. Rev. Microbiol. 6, 613–624. doi: 10.1038/nrmicro193
Hueck, C. J., and Hillen, W. (1995). Catabolite repression in Bacillus subtilis: a global regulatory mechanism for the Gram-positive bacteria? Mol. Microbiol. 15, 395–401. doi: 10.1111/j.1365-2958.1995.tb02252.x
Hueck, C. J., Hillen, W., and Saier, M. H. Jr. (1994). Analysis of a cis-active sequence mediating catabolite repression in Gram-positive bacteria. Res. Microbiol. 145, 503–518. doi: 10.1016/0923-2508(94)90028-0
Jankovic, I., and Brückner, R. (2007). Carbon catabolite repression of sucrose utilization in: catabolite control protein CcpA ensures glucose preference and autoregulatory limitation of sucrose utilization. Oncology 12, 114–120. doi: 10.1159/000096466
Kant, R., Blom, J., Palva, A., Siezen, R. J., and de Vos, W. M. (2011). Comparative genomics of Lactobacillus. Microb. Biotechnol. 4, 323–332. doi: 10.1111/j.1751-7915.2010.00215.x
Kaplan, H., and Hutkins, R. W. (2000). Fermentation of fructooligosaccharides by Lactic acid bacteria and Bifidobacteria. Appl. Environ. Microbiol. 66, 2682–2684. doi: 10.1128/AEM.66.6.2682-2684.2000
Kearns, D. B., and Russell, J. B. (1996). Catabolite regulation in a diauxic strain and a nondiauxic strain of Streptococcus bovis. Curr. Microbiol. 33, 216–219. doi: 10.1007/s002849900102
Kim, J. N., and Burne, R. A. (2017). CcpA and CodY coordinate acetate metabolism in Streptococcus mutans. Appl. Environ. Microbiol. 83:e03274-16. doi: 10.1128/AEM.03274-16
Lambert, J. M., Bongers, R. S., and Kleerebezem, M. (2007). Cre-lox-based system for multiple gene deletions and selectable-marker removal in Lactobacillus plantarum. Appl. Environ. Microbiol. 73, 1126–1135. doi: 10.1128/AEM.01473-06
L’homme, C., Peschet, J. L., Puigserver, A., and Biagini, A. (2001). Evaluation of fructans in various fresh and stewed fruits by high-performance anion-exchange chromatography with pulsed amperometric detection. J. Chromatogr. A 920, 291–297. doi: 10.1016/S0021-9673(00)01262-0
Li, C., Sun, J. W., Zhang, G. F., and Liu, L. B. (2016). Effect of the absence of the CcpA gene on growth, metabolic production, and stress tolerance in Lactobacillus delbrueckii ssp. bulgaricus. J. Dairy Sci. 99, 104–111. doi: 10.3168/jds.2015-10321
Marasco, R., Muscariello, L., Rigano, M., and Sacco, M. (2002). Mutational analysis of the bglH catabolite-responsive element (cre) in Lactobacillus plantarum. FEMS Microbiol. Lett. 208, 143–146. doi: 10.1111/j.1574-6968.2002.tb11074.x
Marciniak, B. C., Pabijaniak, M., Jong, A. D., Dűhring, R., Seidel, G., Hillen, W., et al. (2012). High- and low-affinity cre boxes for CcpA binding in Bacillus subtilis revealed by genome-wide analysis. BMC Genomics 13:401. doi: 10.1186/1471-2164-13-401
Medina-Rivera, A., Abreu-Goodger, C., Thomas-Chollier, M., Salgado, H., Collado-Vides, J., and van Helden, J. (2011). Theoretical and empirical quality assessment of transcription factor-binding motifs. Nucleic Acids Res. 39, 808–824. doi: 10.1093/nar/gkq710
Miwa, Y., Nakata, A., Ogiwara, A., Yamamoto, M., and Fujita, Y. (2000). Evaluation and characterization of catabolite-responsive elements (cre) of Bacillus subtilis. Nucleic Acids Res. 28, 1206–1210. doi: 10.1093/nar/28.5.1206
Moye, Z. D., Lin, Z., and Burne, R. A. (2014). Fueling the caries process: carbohydrate metabolism and gene regulation by Streptococcus mutans. J. Oral Microbiol. 6, 635–661. doi: 10.3402/jom.v6.24878
Muscariello, L., Marasco, R., De Felice, M., and Sacco, M. (2001). The functional ccpA gene is required for carbon catabolite repression in Lactobacillus plantarum. Appl. Environ. Microbiol. 67, 2903–2907. doi: 10.1128/AEM.67.7.2903-2907.2001
Robert, H., Marrec, C. L., Blanco, C., and Jebbar, M. (2000). Glycine betaine, carnitine, and choline enhance salinity tolerance and prevent the accumulation of sodium to a level inhibiting growth of Tetragenococcus halophila. Appl. Environ. Microbiol. 66, 509–517. doi: 10.1128/AEM.66.2.509-517.2000
Robert, K., Christof, F., Michiel, W., and Siezen, R. J. (2008). A generic approach to identify transcription factor-specific operator motifs; inferences for LacI-family mediated regulation in Lactobacillus plantarum WCFS1. BMC Genomics 9:145. doi: 10.1186/1471-2164-9-145
Ryan, S. M., Fitzgerald, G. F., and Van, S. D. (2005). Transcriptional regulation and characterization of a novel beta-fructofuranosidase-encoding gene from Bifidobacterium breve UCC2003. Appl. Environ. Microbiol. 71, 3475–3482. doi: 10.1128/AEM.71.7.3475-3482.2005
Sánchez, B., Noriega, L., Ruas-Madiedo, P., de los Reyes-Gavilán, C. G., and Margolles, A. (2004). Acquired resistance to bile increases fructose-6-phosphate phosphoketolase activity in Bifidobacterium. FEMS Microbiol. Lett. 235, 35–41. doi: 10.1016/j.femsle.2004.04.009
Schumacher, M. A., Sprehe, M., Bartholomae, M., Hillen, W., and Brennan, R. G. (2011). Structures of carbon catabolite protein A–(HPr-Ser46-P) bound to diverse catabolite response element sites reveal the basis for high-affinity binding to degenerate DNA operators. Nucleic Acids Res. 39, 2931–2942. doi: 10.1093/nar/gkq1177
Sørvig, E., Mathiesen, G., Naterstad, K., Eijsink, V. G., and Axelsson, L. (2005). High-level, inducible gene expression in Lactobacillus sakei and Lactobacillus plantarum using versatile expression vectors. Microbiology 151, 2439–2449. doi: 10.1099/mic.0.28084-0
Swint-Kruse, L., and Matthews, K. S. (2009). Allostery in the LacI/GalR family: variations on a theme. Curr. Opin. Microbiol. 12, 129–137. doi: 10.1016/j.mib.2009.01.009
Teixeira, J. S., Abdi, R., Su, S. W., Schwab, C., and Gänzle, M. G. (2013). Functional characterization of sucrose phosphorylase and scrR, a regulator of sucrose metabolism in Lactobacillus reuteri. Food Microbiol. 36, 432–439. doi: 10.1016/j.fm.2013.07.011
Teusink, B., van Enckevort, F. H., Francke, C., Wiersma, A., Wegkamp, A., Smid, E. J., et al. (2005). In silico reconstruction of the metabolic pathways of Lactobacillus plantarum: comparing predictions of nutrient requirements with those from growth experiments. Appl. Environ. Microbiol. 71, 7253–7262. doi: 10.1128/AEM.71.11.7253-7262.2005
Tiffert, Y., Supra, P., Wurm, R., Wohlleben, W., Wagner, R., and Reuther, J. (2008). The Streptomyces coelicolor GlnR regulon: identification of new GlnR targets and evidence for a central role of GlnR in nitrogen metabolism in actinomycetes. Mol. Microbiol. 67, 861–880. doi: 10.1111/j.1365-2958.2007.06092.x
Vastano, V., Perrone, F., Marasco, R., Sacco, M., and Muscariello, L. (2016). Transcriptional analysis of exopolysaccharides biosynthesis gene clusters in Lactobacillus plantarum. Arch. Microbiol. 198, 1–6. doi: 10.1007/s00203-015-1169-1
Wels, M., Overmars, L., Francke, C., Kleerebezem, M., and Siezen, R. J. (2011). Reconstruction of the regulatory network of Lactobacillus plantarum WCFS1 on basis of correlated gene expression and conserved regulatory motifs. Microb. Biotechnol. 4, 333–344. doi: 10.1111/j.1751-7915.2010.00217.x
Willenborg, J., de Greeff, A., Jarek, M., Valentin-Weigand, P., and Goethe, R. (2014). The CcpA regulon of Streptococcus suis reveals novel insights into the regulation of the streptococcal central carbon metabolism by binding of CcpA to two distinct binding motifs. Mol. Microbiol. 92, 61–83. doi: 10.1111/mmi.12537
Wu, Y., Yang, Y., Ren, C., Chen, Y., Yang, S., Yang, G., et al. (2015). Molecular modulation of pleiotropic regulator CcpA for glucose and xylose coutilization by solvent-producing Clostridium acetobutylicum. Metab. Eng. 28, 169–179. doi: 10.1016/j.ymben.2015.01.006
Yang, Y., Zhang, L., Huang, H., Yang, C., Yang, S., Gu, Y., et al. (2017). A flexible binding site architecture provides new insights into CcpA global regulation in Gram-positive bacteria. mBio 8:e02004-16. doi: 10.1128/mBio.02004-16
Yong, J. G. (2005). Genetic Characterization of Fructooligosaccharides Metabolism by Lactobacillus paracasei 1195. Lincoln, NE: University of Nebraska.
Zeng, L., Chakraborty, B., Farivar, T., and Burne, R. A. (2017). Coordinated regulation of the EII(Man) and fruRKI operons of Streptococcus mutans by global and fructose-specific pathways. Appl. Environ. Microbiol. 83:e01403-17. doi: 10.1128/AEM.01403-17
Zomer, A. L., Buist, G., Larsen, R., Kok, J., and Kuipers, O. P. (2007). Time-resolved determination of the CcpA regulon of Lactococcus lactis subsp. cremoris MG1363. J. Bacteriol. 189, 1366–1381. doi: 10.1128/JB.01013-06
Keywords: Lactobacillus plantarum, fructooligosaccharide, catabolite control protein A, carbon catabolite repression, metabolic regulation
Citation: Chen C, Lu Y, Wang L, Yu H and Tian H (2018) CcpA-Dependent Carbon Catabolite Repression Regulates Fructooligosaccharides Metabolism in Lactobacillus plantarum. Front. Microbiol. 9:1114. doi: 10.3389/fmicb.2018.01114
Received: 17 March 2018; Accepted: 11 May 2018;
Published: 29 May 2018.
Edited by:
Zhihong Sun, Inner Mongolia Agricultural University, ChinaReviewed by:
Kuan Wang, Tianjin University of Science and Technology, ChinaLianzhong Ai, University of Shanghai for Science and Technology, China
Copyright © 2018 Chen, Lu, Wang, Yu and Tian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Huaixiang Tian, dGlhbmh4QHNpdC5lZHUuY24=