Skip to content


  • Open Access

Biosensor-guided improvements in salicylate production by recombinant Escherichia coli

Microbial Cell Factories201918:18

  • Received: 28 November 2018
  • Accepted: 20 January 2019
  • Published:



Salicylate can be biosynthesized from the common metabolic intermediate shikimate and has found applications in pharmaceuticals and in the bioplastics industry. While much metabolic engineering work focused on the shikimate pathway has led to the biosynthesis of a variety of aromatic compounds, little is known about how the relative expression levels of pathway components influence salicylate biosynthesis. Furthermore, some host strain gene deletions that improve salicylate production may be impossible to predict. Here, a salicylate-responsive transcription factor was used to optimize the expression levels of shikimate/salicylate pathway genes in recombinant E. coli, and to screen a chromosomal transposon insertion library for improved salicylate production.


A high-throughput colony screen was first developed based on a previously designed salicylate-responsive variant of the E. coli AraC regulatory protein (“AraC-SA”). Next, a combinatorial library was constructed comprising a series of ribosome binding site sequences corresponding to a range of predicted protein translation initiation rates, for each of six pathway genes (> 38,000 strain candidates). Screening for improved salicylate production allowed for the rapid identification of optimal gene expression patterns, conferring up to 123% improved production of salicylate in shake-flask culture. Finally, transposon mutagenesis and screening revealed that deletion of rnd (encoding RNase D) from the host chromosome further improved salicylate production by 27%.


These results demonstrate the effectiveness of the salicylate sensor-based screening platform to rapidly identify beneficial gene expression patterns and gene knockout targets for improving production. Such customized high-throughput tools complement other cell factory engineering strategies. This approach can be generalized for the production of other shikimate-derived compounds.


  • Metabolic engineering
  • Synthetic biology
  • Salicylate
  • AraC
  • Biosensor
  • High-throughput screening
  • Ribosome binding site


Metabolic engineering facilitates improved or novel biological conversion of low-value, renewable feedstocks into higher value chemicals. By interfacing native cell metabolism with heterologous pathways, engineered microorganisms have produced a large portfolio of compounds, ranging from bulk chemicals (e.g. ethanol and butanol), to fine chemicals (e.g. aromatics) and pharmaceutical precursors (e.g. isoprenoids and alkaloids). Advances in synthetic biology are increasingly allowing for the identification and maintenance of balanced expression levels of multiple genes within a pathway [13]. Methods of fine-tuning and even dynamically regulating gene expression at the transcriptional and/or translational levels have been developed to improve flux through engineered pathways [48]. Meanwhile, metabolite-responsive transcription factors enable detection of metabolic state, where the effector molecule triggers expression of a reporter gene for in situ phenotypic screening [1].

Phenolic compounds are popular metabolic engineering targets, due to their applications as commodity chemicals [9] and pharmaceutical precursors (e.g. anti-depressant drugs, and antitumor drugs) [10]. In microbes, the shikimate pathway is the primary biosynthetic route for synthesizing aromatic amino acids and their derivatives. The key step is a condensation of glycolytic intermediate phosphoenolpyruvate (PEP) and pentose phosphate pathway intermediate erythrose-4-phosphate (E4P), generating the precursor chorismate. The shikimate pathway has been engineered to synthesize a range of natural products including flavonoids [11] and phenylpropanoids [12]. In this study, we engineered the E. coli shikimate pathway for production of the chorismate derivative salicylate. Salicylate can serve as a precursor to a range of pharmaceuticals such as aspirin and lamivudine (an anti-HIV drug), and is a starting material to cis–cis-muconic acid and adipic acid, the building block for nylon-6,6 and polyurethane [13]. Previous metabolic engineering work led to the functional expression of a salicylate pathway in E. coli [14]. Heterologous expression of isochorismate synthase (ICS encoded by entC from E. coli) and isochorismate pyruvate lyase (IPL encoded by pchB from P. fluorescens Migula) in a phenylalanine-overproducing strain resulted in the accumulation of 8.5 mM salicylate [14]. More recently, the PEP-consuming phosphotransferase system encoded by ptsH and ptsI was replaced with the heterologous expression of galactose permease and glucokinase, and genes encoding PEP-consuming pyruvate kinases (pykA and pykF) were deleted, resulting in an engineered strain that produced 80 mM salicylate in a bench-top bioreactor under optimized conditions [15].

Here we describe work to further improve salicylate production in E. coli by optimizing gene expression of the key enzymes in the shikimate pathway, with the help of a high-throughput sensor-reporter system. We previously engineered the E. coli transcription factor AraC to respond to and hence report on the presence of a number of non-native effector molecules, including triacetic acid lactone (TAL), mevalonate, and phenolic compounds [1620]. In this study, a salicylate sensor-reporter system based on variant “AraC-SA” was rewired to control the chromosomal expression of β-galactosidase (LacZ), which cleaves X-gal and generates blue color.

As depicted in Fig. 1, the designed reporter system was applied to optimize the expression level of key enzymes in the shikimate pathway. A series of ribosome binding site (RBS) sequences corresponding to a spectrum of predicted translation rates were designed for each of six genes involved in the heterologous salicylate biosynthesis pathway (Fig. 1a). Screening the combinatorial expression library resulted in isolation of several clones with enhanced salicylate production (up to 123% improvement). Subsequently, screening a random transposon insertion library identified a knock-out gene target that further improved salicylate titer by 27%.
Fig. 1
Fig. 1

Sensor-based pathway optimization for improving salicylate production. a Synthetic pathway for producing salicylate from glycerol. Six genes highlighted in blue encode the key enzymes involved in salicylate biosynthesis: ppsA encodes PEP synthetase, which converts pyruvate to PEP. tktA encodes a transketolase that increases supply of E4P. DAHP synthase encoded by aroGfbr catalyzes the aldol condensation to generate DAHP. Shikimate kinase encoded by aroL produces chorismate. Isochorismate synthase (encoded by entC) and isochorismate pyruvate lyase (encoded by pchB) enable the synthesis of salicylate form chorismate. b AraC-based salicylate biosensor was used to screen for clones showing enhanced salicylate production


Construction of AraC-based sensor system for agar plate-based high-throughput screening

AraC variant “AraC-SA” (P8G, T24L, H80G, Y82I and H93R) responds to salicylate, as described previously [20]. To assess the specificity of this sensor, we compared the response (fold activation of GFP expression under the control of promoter PBAD) to salicylate with that of several other similar molecules (substituted benzoates) (Fig. 2b). A concentration of 5 mM was used to compare compounds, as this amount is relevant to screening salicylate production, without major growth inhibition. Salicylate induces GFP expression over 200-fold, while there is essentially no response to the other compounds tested. The complete AraC-SA salicylate dose response is provided as Additional file 1.
Fig. 2
Fig. 2

Development of a salicylate-responsive sensor-reporter system. a Heterologous pathway for producing salicylate from the native E. coli central metabolite chorismate. b Fold activation of GFP expression under control of AraC-SA/PBAD, in the presence of 5 mM salicylate (SA) and similar compounds (l-ara: l-arabinose (inducer for wild-type AraC); BA: benzoic acid; 2-MOBA: 2-methoxybenzoic acid; 3-HBA: 3-hydroxybenzoic acid; 4-HBA: 4-hydroxybenzoic acid; and 2MBA: 2-methylbenzoic acid. Reported data represents mean ± std. dev., from four biological replicates. c Response of AraC sensor in E. coli strain QH4 carrying plasmid pGF29-SA (neither entC nor pchB is expressed), plasmid pFG29-SA-entC (entC but not pchB are expressed), and plasmid pFG29-SA-EP (enabling salicylate production). d Colorimetric colony assay using the designed salicylate sensor system (SQ22 expressing AraC-SA growing on LB-agar plates containing X-GAL and different concentrations of salicylate)

To confirm that AraC-SA does not also respond to compounds upstream in the salicylate biosynthetic pathway (namely chorismate and isochorismate), strains expressing the full or partial salicylate pathway with the AraC-SA/GFP module were constructed (Fig. 2c). Only the strain expressing the full pathway produced salicylate (~ 4.5 mM after 14 h post-induction), resulting in 15-fold increased GFP fluorescence. Strains without the complete pathway exhibited only basal fluorescence signal. The AraC-SA sensor/reporter thus responds specifically to salicylate, at least in the context of salicylate biosynthesis in recombinant E. coli.

In using our sensor/reporter system to screen libraries for mutants with enhanced salicylate production, a colony-based screen was developed [21]. This approach provides a simple means of spatially isolating individual mutants and preventing “cheaters” from responding to salicylate produced and secreted by other mutants. To that end, a single copy of lacZ under control of PBAD was integrated into the salicylate production host strain (QH4) [14], resulting in reporter strain SQ22 for library screening. As shown in Fig. 2d, in the presence of X-gal (5-bromo-4-cholo-3-indolyl β-d-galactopyranoside), the colony’s blue color intensity resulting from LacZ-catalyzed X-gal hydrolysis correlates with the salicylate concentration. Further validation of this colony screen for isolating clones with enhanced salicylate production is provided in the (Additional file 2).

Salicylate biosynthetic pathway RBS library construction and screening

Given the known metabolic burden imposed by gene overexpression, and potential toxicities or enzyme inhibition resulting from accumulation of pathway intermediates [22, 23], we reasoned that one approach to improving salicylate production in recombinant E. coli is to fine-tune the relative expression levels of the six key biosynthetic pathway genes (described above). The RBS calculator [7] was used to design a series of RBS sequences resulting in a range of predicted translation initiation strengths for the six targeted genes (entC, pchB, aroL, ppsA, tktA, aroGfbr) (refer to Additional file 3 for details). These six genes with varied RBS sequences were assembled combinatorially into plasmid pPCC1244, leading to pathway library “QSAlib3” composed of theoretically 38,880 unique members (Fig. 3a). Strain SQ22 containing a chromosomal PBAD-lacZ reporter was transformed with the QSAlib3 library and spread onto large, X-Gal-containing plates for screening. Strain SQ22 is derived from strain QH4 [14] having pheA and tyrA gene deletions (Table 1), which prevents conversion of chorismate to phenylalanine and tyrosine. Yeast extract (5 g/L) was therefore added to our modified minimal medium to support growth.
Fig. 3
Fig. 3

Implementation of sensor-reporter screening to optimize expression of salicylate pathway genes. a Schematic of RBS library plasmid constructs. Six RBS sequences ranging from weakest (“gene-1”) to strongest (“gene-6”) calculated translation initiation rates were placed upstream of each pathway gene on plasmid pPCC1253, also containing the AraC-SA sensor. b Cell density and salicylate titers measured from the top clones identified from screening. c List of RBS combinations for each gene from the top six selected clones. Reported data represents mean ± std. dev., from four independent experiments. The corresponding predicted TIRs are shown in an analogous table in the (Additional file 3). Note that TIRs calculated for a specific gene are relative numbers for predicting the targeted protein relative expression level, and the actual protein expression level is affected by other factors as well. Hence it should not be considered constructive to compare calculated TIRs for different genes

Table 1

Plasmids and strains used in this study


Relevant features



entC, pchB, and luc genes under control of Plac; AmpR



aroL ppsA, tktA, aroGfbr under control of Plac; KanR



CRIM plasmid with PBAD-lacZ;;AprR

[21, 29]


araC-SA gene under the control of Ptac; AprR

In this study


entC, pchB genes under control of Ptac; AmpR

In this study


araC-SA, entC, pchB genes under control of Ptac; AprR

In this study


aroL ppsA, tktA, aroGfbr genes under control of Ptac; AprR

In this study


araC-SA, entC, pchB genes under control of Ptac; aroL ppsA, tktA, aroGfbr genes under control of Ptac; AprR

In this study


araC-SA, entC, pchB genes under control of Ptac; ppsA, tktA, aroGfbr genes under control of Ptac; AprR

In this study


From vector pJA1. Transposable element containing a chloramphenicol resistance gene flanked by IS10 inverted repeat sequences; R6 K replication origin; tnp gene encoding a mutant Tn10 transposase with 100-fold lower frequency of insertion at hot spots; AmpR

[21, 31]


araC gene under the control of Ptac, tac; AprR



araC-SA gene under the control of Ptac; gfpuv under control of PBAD promoter, AprR




BW27786 (∆araFGH, ∆araBAD,lacZ), with araC deleted



ATCC31884 with pheA and tyrA disrupted



QH4 with PBAD-lacZ reporter integrated into chromosome at HK022 site

In this study


SQ18 with LacZ in the lac operon disrupted

In this study


QH4 with rnd gene disrupted

In this study


SQ22 with rnd gene disrupted

In this study

Approximately 100,000 colonies were screened by eye. Rescreening of clones representing the 75 darkest blue colonies was performed by quantifying salicylate production in liquid culture (using strain SQ22 harboring the selected mutant plasmids). Among these, the top six mutant plasmids based on salicylate titer were sequenced (Fig. 3c) and salicylate production using these constructs was further characterized with strain QH4 (Fig. 3b). Two clones produced about 13.8 ± 1.2 mM salicylate, which is > 120% improved relative to strain QH4 overexpressing all six target genes, each with the same “original” RBS sequence (AGGAGA) (plasmid pPCC1253). The time courses of salicylate production and cell density for QH4 + pPCC1253 and QH4 + pQSA-50 (provided in Additional file 4) show that salicylate production occurs primarily after growth has already slowed and culture densities are near maximum.

As shown in Fig. 3c, the RBS sequences appearing in the top six selected clones encompass the full range of weak to strong, for all six genes. This gene expression library is represented on polycistronic operons (Fig. 3a), so translational coupling [24] may at least partly explain the general lack of patterns to RBS strengths for the selected clones. While there are no prominent expression patterns for these genes, we note that most clones contain a relatively strong RBS upstream of tktA, and relatively weak RBS sequences upstream of aroL. SDS-PAGE analysis of cell lysates (using Coomassie blue staining) confirms these observations, in that strong bands are present for transketolase (though less prominent for the weaker RBSs in pQSA-78 and pQSA-94), and the two relatively stronger RBS sequences for aroL (pQSA-43 and pQSA-94) show higher aroL expression than the others (Additional file 5). We tested whether further lowering aroL expression by removing this gene from the plasmid (construct pPCC1253-aroL), and relying on only chromosomal expression of this gene, would further improve salicylate production. This however significantly reduced salicylate production (Fig. 3b). These results demonstrate the challenges of rationally tuning pathway gene expression levels, and the value in using combinatorial approaches to identify improvements.

Screening a transposon insertion library for further improvements in salicylate production

We next sought to mutate the host genome and screen for further improvements in salicylate production. A suicide plasmid vector containing a transposable DNA fragment and IPTG-inducible transposase was used to generate the transposon insertion library in strain SQ22 harboring plasmid pQSA-50. The SQ22-transposon insertion library was screened using our AraC-PBAD-lacZ sensor-reporter system. After screening ~ 70,000 colonies on plates, 50 selected clones were re-screened in liquid culture for HPLC analysis. Among these, clone QSA#228 consistently produced ~ 16.5 mM salicylate after 48 h induction, which is ~ 27% improved over parent strain SQ22 harboring mutant plasmid pQSA-50 (Fig. 4). Curing the originally selected strain of plasmid, and re-transforming with fresh plasmid pQSA-50 resulted in the same improvement in salicylate titer. Self-formed adaptor PCR (SEFA-PCR) analysis [25] located the transposable element within the rnd gene, encoding RNase D. RNase D is an exonuclease involved in the 3′ ribonucleotide processing of tRNA [26]. Direct deletion of rnd from the genomes of host strains SQ22 and QH4 resulted in similar improvements in salicylate production. It is not known how changes in protein expression levels due to the rnd mutation are beneficial to salicylate biosynthesis.
Fig. 4
Fig. 4

Transposon insertion library screening for improved salicylate production led to the identification of an rnd mutation in the host. Red column: Salicylate titer after 48 h fermentation of different host strains expressing plasmid pQSA-50; Black column: cell density of culture after 48 h fermentation. Reported data represents mean ± std. dev., from four independent experiments


Salicylate has been widely used as a precursor in the pharmaceuticals and bioplastics industry. Metabolic engineers have successfully synthesized salicylate from recombinant E. coli by extending the shikimate pathway [14, 15]. Here we further improved salicylate production in E. coli by optimizing expression of key enzymes in the shikimate pathway, with the help of a high-throughput sensor-reporter system. In addition, screening a transposon library with the same reporter system revealed a gene knockout target that further contributes to salicylate biosynthesis. The resulting salicylate titer of 17 mM is the highest reported from shake-flask culture.

It is worth noting that all six of the highest salicylate-producing mutants presented in Fig. 3b (results from RBS library screening) show similar final titers, suggesting a potential “upper limit” to production, at least for this library and under these culture conditions. At least 40 mM glycerol remained in all cultures tested (data not shown), suggesting carbon supply does not limit production. The inhibitory effects of salicylate on growth of strain QH4 are notable even at 5 mM, and more severe as salicylate concentration increases (refer to Additional file 6). However, as noted above, salicylate titers do not exceed 5 mM until culture densities are already near maximum, and maximum culture density is not reduced as a result of increased salicylate production (Fig. 3b and Additional file 4). Nonetheless, salicylate toxicity may contribute to an upper limit on production. Further strain engineering to improve salicylate tolerance and/or export may therefore be promising approaches to further advance the commercial potential of microbial salicylate production.



Restriction enzymes, Phusion High-Fidelity DNA polymerase, and T4 DNA ligase were purchased from New England Biolabs (Ipswich, MA). Oligonucleotides were synthesized by Integrated DNA Technologies (Coralville, IA). DNA sequencing was performed at SeqWright (Houston, TX) or Genewiz (South Plainfield, NJ). All chemicals were purchased from Sigma-Aldrich (St. Louis, MO). Molecular biology techniques for DNA manipulation were performed according to standard protocols [27] and all cultures were grown in lysogeny broth (LB). Antibiotics and IPTG were prepared as a 1000× stock solution in purified water and sterile filtered with EMD Millipore Millex-GP syringe driven filters (EMD Millipore, Cat. No. SLGP033RS). The modified NBS medium contains MOPS (50 mM, pH 7.4) and per liter: 10 g glycerol; 2.5 g glucose; 5 g yeast extract; 3.4 g KH2PO4; 5.2 g K2HPO4; 3.3 g (NH4)2HPO4, 0.25 g MgSO4·7 H2O, 15 mg CaCl2·2 H2O, 0.5 mg thiamine, and 1 mL of trace metal stock (described by Chin and Cirino [28]). The concentrations of antibiotics used for maintaining the plasmids are as follows: apramycin 50 µg/ml, chloramphenicol 25 µg/ml, ampicillin 100 µg/ml.


Plasmids and strains used in this study are listed in Table 1. All primers used in this study are listed in Additional file 7. Primers pPCC1244-gib-for and pPCC1244-gib-rev were used to amplify lacI-Ptac-araC-SA DNA fragment from pFG29-SA, and the PCR product was assembled with PciI- and XmnI-digested pFG1, resulting in plasmid pPCC1244.

To create plasmid pPCC1250, pZE-EP was double digested with Bsu36I and SphI, and the digested DNA fragment was subjected to a DNA-blunting reaction. The blunt-ended DNA fragment was self-ligated, resulting in plasmid pPCC1250. A DNA fragment containing entC and pchB was amplified from pPCC1250 using primers pPCC1251-gib-for and pPCC1251-gib-rev2, and assembled with XmnI-digested pPCC1244, resulting in plasmid pPCC1251. A DNA fragment containing aroL ppsA, tktA, aroGfbr genes was amplified from pCS-APTA using primers pPCC1252-gib-for and pPCC1251-gib-rev2, and the amplified DNA fragment was assembled with NdeI- and XmnI-digested pPCC1244, resulting in plasmid pPCC1252. Then, a DNA fragment of Ptac-araC-SA-entC-pchB was amplified from pPCC1251 using primers pPCC1253-NheI-for and pPCC1253-NotI-rev, and digested with NheI and NotI. The digested DNA fragment was ligated into pPCC1252 digested with the same restriction enzymes, resulting in plasmid pPCC1253. To remove aroL from plasmid pPCC1253, primers pPCC1251-gib-for and pPCC1253-Ptac-rev were used to amplify the entC-pchB DNA fragment from pPCC1253, and primers pPCC1253-ppsA-for and pPCC1253-tktA-rev were used to amplify the ppsA-tktA section from pPCC1253. The two amplified PCR products were assembled into KpnI/SphI double-digested pPCC1253, resulting in plasmid pPCC1253-aroL.

Construction of pathway RBS library

The RBS calculator [7] was used to design six RBS sequences having different translation initiation rates (TIRs) for each gene (see Additional file 1 for the RBS sequences and calculated TIRs). For each gene, the six designed upstream primers containing different RBS sequences were mixed equimolar, resulting in primer mixtures entC-RBSs-for, pchB-RBSs-for, etc., and each primer mixture was used to construct the RBS library. Primers pFG29_araC_GS_fwd_1 and AraC-gib-rev were used to amplify Ptac-AraC-SA DNA fragment from pPCC1244. Primers entC-RBSs-for and entC-RBS-rev were used to amplify entC from pZE-EP, and primers pchB-RBSs-for and pPCC1251-gib-rev2 were used to amplify pchB from pZE-EP. Plasmid pPCC1244 was double digested by BstAPI and XmnI, and the linearized vector was assembled with the three PCR products, resulting in QSALib1, which contains the RBS libraries for entC and pchB genes.

Primers pFG29_araC_GS_fwd_1 and Ptac-gib-rev were used to amplify a DNA fragment containing promoter Ptac. The aroL gene was amplified from pCS-APTA by primers aroL-RBSs-for and AroL-RBS-rev. Similarly, ppsA was amplified from pCS-APTA using primers ppsA-RBSs-for and ppsA-RBS-rev. Overlap extension PCR (OE-PCR) was performed to assemble these three PCR products using primers pFG29_araC_GS_fwd_1 and QSAlib2-OE123-rev, resulting in DNA fragment QSAlib2-f123. Next, primers tktA-RBSs-for and tktA-RBS-rev were used to amplify tktA from plasmid pCS-APTA, and primers AroG-RBSs-for and pPCC1251-gib-rev2 were used to amplify aroG. tktA and aroG were assembled by OE-PCR using primers QSAlib2-OE45-for and QSAlib2-OE45-rev, resulting in DNA fragment QSAlib2-f45. pPCC1252 was then double digested with BstAPI/BamHI, QSAlib2-f123 was double digested with BstAPI/SpeI, and QSAlib2-f45 was double digested with SpeI/BamHI. Ligation of these three digest fragments resulted in library QSAlib2.

Finally, primers pPCC1253-NheI-for and pPCC1253-NotI-rev were used to amplify Ptac-araC-SA-RBS-entC-RBS-pchB from QSAlib1. The PCR product was digested with NheI and NotI, and ligated with QSAlib3 digested with the same enzymes, resulting in QSAlib3. QSAlib3 contains the RBS libraries for all six genes. Sanger sequencing of QSAlib3 confirmed proper library construction.

Strain construction

The strains used in this study are listed in Table 1. Plasmid pPCC1155-5 was integrated into the chromosome of QH4 at site HK022 as described [29]. Removal of FRT-flanked apramycin resistance cassette resulted in strain SQ18. A Phage λ Red disruption method [30] was used to delete lacZ from the lac operon of strain SQ18, resulting in strain SQ22. Deletion of rnd in strains QH4 and SQ22 was similarly performed, resulting in strains QH4∆rnd and SQ22∆rnd, respectively.

Sensor-reporter fluorescence assays

Essentially as described [20], 500 μl LB + apramycin in 2-ml-well, 96-well plate was inoculated with strain HF19 harboring pFG29-SA. These starter cultures were incubated for 6 h at 37 °C, 900 rpm, then diluted to OD595 = 0.05 in 500 μl “biosensor medium” containing different concentrations of the compound of interest. After 6 h, the cultures were pelleted and washed with an equal volume of phosphate-buffered saline, before measuring OD595 and fluorescence (400 nm excitation, 510 nm emission) using plate readers.

Salicylate production in baffled flasks

A colony of the salicylate-producing strain was used to inoculate 3 ml LB + apramycin, and grown in a test tube for 8 h at 37 °C and 250 rpm. This seed culture was then diluted to OD595 = 0.05 into 25 ml modified NBS medium containing apramycin and 250 mM IPTG, in 125 ml baffled flasks. The flasks were shaken at 37 °C and 250 rpm for 48 h, at which time OD595 values were measured and salicylate concentrations were analyzed by HPLC.

Screening the RBS library for improved salicylate producers

Strain SQ22 was transformed with QSAlib3. The outgrowth was transferred into LB + apramycin, and grown at 37 °C, 250 rpm for 12 h. The resulting culture was diluted and spread onto large plates containing modified NBS-agar with IPTG (250 µM), X-Gal (40 µg/ml), and apramycin. After 24 h of incubation, the top 5 blue colonies from each screening plate were picked and streaked onto fresh LB plates supplied with apramycin. The resulting colonies were tested for salicylate production in liquid culture.

Construction and screening of SQ22 transposon insertion library

Strain SQ22 harboring the highest-producing salicylate plasmid (pQSA-50) was transformed with 1 µg of plasmid pPCC507, and the outgrowth was grown in 1 ml SOB medium supplied with 20 µM IPTG at 37 °C for 1 h. The outgrowth was transferred into 500 ml LB + apramycin and 12.5 µg/ml chloramphenicol, and grown at 37 °C for 12 h. The resulting culture was diluted and plated on modified NBS-agar plates containing IPTG (250 µM), X-Gal (40 µg/ml), apramycin, and 12.5 μg/ml chloramphenicol. Totally 70,000 colonies were screened on 10 screening plates. After 24 h of incubation, the top 5 blue colonies from each screening plate were picked and streaked onto fresh LB plates containing apramycin and 12.5 µg/ml chloramphenicol. The resulting colonies were tested for salicylate production in liquid culture.

Salicylate quantification by HPLC

500 µl of cell culture was centrifuged at 17,900×g, and the supernatant was filtered through a 0.45 µm filter. The salicylate concentration in the filtrate was determined by reverse-phase HPLC using a C18 column on a Shimadzu LC-20AD HPLC system (Kyoto, Japan) equipped with a UV monitor. The elution profile was as follows: Solvent A, 1% (v/v) acetic acid in water; solvent B, 1% (v/v) acetic acid in acetonitrile; gradient: 5–95% B (0–15 min), 95–5% B (15–17 min), 5% B (17–20 min). The column temperature was set to 50 °C. Salicylate eluted around 11.2 min at a flow rate of 1 ml/min. Elution absorbance at 310 nm was monitored and peak areas were converted to sample concentrations based on calibration with pure salicylate.





benzoic acid


dihydroxyacetone phosphate












3-hydroxybenzoic acid


4-hydroxybenzoic acid




2-methoxybenzoic acid


2-methylbenzoic acid






ribosome binding site








translation initiation rate




Authors’ contributions

SQ and YL conducted the experiments and analyzed the data. SQ and PCC designed the experiments and prepared the manuscript. All authors read and approved the final manuscript.


This project was supported by the National Science Foundation (Award Number: CBET-1511425). We thank Dr. Yajun Yan (University of Georgia) for providing strain QH4 and plasmids pZE-EP and pCS-APTA.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its additional files.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


Indicated in the acknowledgments.

Publisher’s Note

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

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Department of Chemical & Biomolecular Engineering, University of Houston, S222 Engineering Building 1, Houston, TX 77204-4004, USA


  1. Qian S, Cirino PC. Using metabolite-responsive gene regulators to improve microbial biosynthesis. Curr Opin Chem Eng. 2016;14:93–102.View ArticleGoogle Scholar
  2. Xu P. Production of chemicals using dynamic control of metabolic fluxes. Curr Opin Biotechnol. 2018;53:12–9.View ArticleGoogle Scholar
  3. Wang Z, Cirino PC. New and improved tools and methods for enhanced biosynthesis of natural products in microorganisms. Curr Opin Biotechnol. 2016;42:159–68.View ArticleGoogle Scholar
  4. Juminaga D, Baidoo EEK, Redding-Johanson AM, Batth TS, Burd H, Mukhopadhyay A, et al. Modular engineering of l-tyrosine production in Escherichia coli. Appl Environ Microbiol. 2012;78:89–98.View ArticleGoogle Scholar
  5. Anthony JR, Anthony LC, Nowroozi F, Kwon G, Newman JD, Keasling JD. Optimization of the mevalonate-based isoprenoid biosynthetic pathway in Escherichia coli for production of the anti-malarial drug precursor amorpha-4,11-diene. Metab Eng. 2009;11:13–9.View ArticleGoogle Scholar
  6. Bokinsky G, Peralta-Yahya PP, George A, Holmes BM, Steen EJ, Dietrich J, et al. Synthesis of three advanced biofuels from ionic liquid-pretreated switchgrass using engineered Escherichia coli. Proc Natl Acad Sci USA. 2011;108:19949–54.View ArticleGoogle Scholar
  7. Salis HM, Mirsky EA, Voigt CA. Automated design of synthetic ribosome binding sites to control protein expression. Nat Biotechnol. 2009;27:946–50.View ArticleGoogle Scholar
  8. Zelcbuch L, Antonovsky N, Bar-Even A, Levin-Karp A, Barenholz U, Dayagi M, et al. Spanning high-dimensional expression space using ribosome-binding site combinatorics. Nucleic Acids Res. 2013;41:98–98.View ArticleGoogle Scholar
  9. Gosset G. Production of aromatic compounds in bacteria. Curr Opin Biotechnol. 2009;20:651–8.View ArticleGoogle Scholar
  10. Jiang M, Zhang H. Engineering the shikimate pathway for biosynthesis of molecules with pharmaceutical activities in E. coli. Curr Opin Biotechnol. 2016;42:1–6.View ArticleGoogle Scholar
  11. Xu P, Ranganathan S, Fowler ZL, Maranas CD, Koffas MAG. Genome-scale metabolic network modeling results in minimal interventions that cooperatively force carbon flux towards malonyl-CoA. Metab Eng. 2011;13:578–87.View ArticleGoogle Scholar
  12. Lin Y, Sun X, Yuan Q, Yan Y. Combinatorial biosynthesis of plant-specific coumarins in bacteria. Metab Eng. 2013;18:69–77.View ArticleGoogle Scholar
  13. Sun X, Lin Y, Huang Q, Yuan Q, Yan Y. A novel muconic acid biosynthesis approach by shunting tryptophan biosynthesis via anthranilate. Appl Environ Microbiol. 2013;79:4024–30.View ArticleGoogle Scholar
  14. Lin Y, Sun X, Yuan Q, Yan Y. Extending shikimate pathway for the production of muconic acid and its precursor salicylic acid in Escherichia coli. Metab Eng. 2014;23:62–9.View ArticleGoogle Scholar
  15. Noda S, Shirai T, Oyama S, Kondo A. Metabolic design of a platform Escherichia coli strain producing various chorismate derivatives. Metab Eng. 2016;33:119–29.View ArticleGoogle Scholar
  16. Tang S-Y, Fazelinia H, Cirino PC. AraC regulatory protein mutants with altered effector specificity. J Am Chem Soc. 2008;130:5267–71.View ArticleGoogle Scholar
  17. Tang SY, Cirino PC. Design and application of a mevalonate-responsive regulatory protein. Angew Chem Int Ed. 2011;50:1084–6.View ArticleGoogle Scholar
  18. Tang SY, Qian S, Akinterinwa O, Frei CS, Gredell JA, Cirino PC. Screening for enhanced triacetic acid lactone production by recombinant Escherichia coli expressing a designed triacetic acid lactone reporter. J Am Chem Soc. 2013;135:10099–103.View ArticleGoogle Scholar
  19. Frei CS, Wang Z, Qian S, Deutsch S, Sutter M, Cirino PC. Analysis of amino acid substitutions in AraC variants that respond to triacetic acid lactone. Protein Sci. 2016;25:804–14.View ArticleGoogle Scholar
  20. Frei CS, Qian S, Cirino PC. New engineered phenolic biosensors based on the AraC regulatory protein. Protein Eng Des Sel. 2018;31:213–20.View ArticleGoogle Scholar
  21. Li Y, Qian S, Dunn R, Cirino PC. Engineering Escherichia coli to increase triacetic acid lactone (TAL) production using an optimized TAL sensor-reporter system. J Ind Microbiol Biotechnol. 2018;45:789–93.View ArticleGoogle Scholar
  22. Dahl RH, Zhang F, Alonso-Gutierrez J, Baidoo E, Batth TS, Redding-Johanson AM, et al. Engineering dynamic pathway regulation using stress-response promoters. Nat Biotechnol. 2013;31:1039–46.View ArticleGoogle Scholar
  23. Zhang F, Carothers JM, Keasling JD. Design of a dynamic sensor-regulator system for production of chemicals and fuels derived from fatty acids. Nat Biotechnol. 2012;30:354–9.View ArticleGoogle Scholar
  24. Levin-Karp A, Barenholz U, Bareia T, Dayagi M, Zelcbuch L, Antonovsky N, et al. Quantifying translational coupling in E. coli synthetic operons using RBS modulation and fluorescent reporters. ACS Synth Biol. 2013;2:327–36.View ArticleGoogle Scholar
  25. Wang S, He J, Cui Z, Li S. Self-formed adaptor PCR: a simple and efficient method for chromosome walking. Appl Environ Microbiol. 2007;73:5048–51.View ArticleGoogle Scholar
  26. Reuven NB, Deutscher MP. Multiple exoribonucleases are required for the 3′ processing of Escherichia coli tRNA precursors in vivo. FASEB J. 1993;7:143–8.View ArticleGoogle Scholar
  27. Sambrook J, Russell DW. Molecular cloning : a laboratory manual. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 2001.Google Scholar
  28. Chin JW, Cirino PC. Improved NADPH supply for xylitol production by engineered Escherichia coli with glycolytic mutations. Biotechnol Prog. 2009;27:333–41.View ArticleGoogle Scholar
  29. Haldimann A, Wanner BL. Conditional-replication, integration, excision, and retrieval plasmid-host systems for gene structure-function studies of bacteria. J Bacteriol. 2001;183:6384–93.View ArticleGoogle Scholar
  30. Datsenko KA, Wanner BL. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci. 2000;97:6640–5.View ArticleGoogle Scholar
  31. Badarinarayana V, Estep PW III, Shendure J, Edwards J, Tavazoie S, Lam F, et al. Selection analysis of insertional mutants using subgenic-resolution arrays. Nat Biotechnol. 2001;19:1060–5.View ArticleGoogle Scholar


© The Author(s) 2019