Team Göttingen
iGEM 2018
Glyphosate on my plate?
Results
- Identification of glyphosate uptake systems
- Identification of glyphosate-regulated promoters
- A reporter-based glyphosate detection system
- Competition assay for glyphosate detection
- Engineering bacteria to disarm glyphosate
- Characterization of EPSP synthases
- References
1. Identification of glyphosate uptake systems
1.1. Interaction between glyphosate and Bacillus subtilis
To create a proper detection system using the Gram-positive model bacterium Bacillus subtilis, we first had to evaluate how this organism behaves in the presence of glyphosate. Previously, it was shown that glyphosate negatively affects growth of B. subtilis due to the inhibition of the EPSP synthase AroE (Figure 1A) (1). Moreover, it has been demonstrated that 1.8 mM of glyphosate was required to inhibit the growth rate by 50%. To re-evaluate the effect of glyphosate on growth of our B. subtilis laboratory strains 168 and SP1, we performed growth experiments in CS-Glc minimal medium that was supplemented with increasing amounts of glyphosate (Figure 1B). CS-Glc medium contains glucose and succinate as carbon sources and ammonium as the nitrogen source (see Notebook). While the B. subtilis strain 168 is auxotrophic for tryptophan due to a mutation in the trpC gene, the strain SP1 can grow in the absence of exogenous tryptophan (2).
Figure 1. (A) Glyphosate inhibits the 5-enolpyruvyl-shikimate-3-phosphate (EPSP) synthase, which converts phosphoenolpyruvate (PEP) and 3-P-shikimate into EPSP and P in B. subtilis. EPSP is a precursor for the aromatic amino acids phenylalanine (Phe), tyrosine (Tyr) and tryptophane (Trp). (B) How does glyphosate affect growth of our B. subtilis wild type strains?
As shown in Figures 2A and 2B, at a glyphosate concentration of about 1 mM the growth rate of both strains was reduced by 50% and the bacteria were not able to grow at glyphosate concentrations higher than 3 mM. In contrast to a previous study (1), we have observed that 44% fewer glyphosate is required to reduce the growth rate of B. subtilis by 50%. This discrepancy might be due to differences in the genetic makeup of the B. subtilis strains, in the medium composition, in the purity of glyphosate or due to the different cultivation conditions. However, glyphosate negatively affects growth of B. subtilis in CS-Glc minimal medium.
Figure 2. (A) Growth of the B. subtilis wild type strain 168 with increasing amounts of glyphosate (GS). (B) Growth of the B. subtilis wild type strain SP1 with increasing amounts of GS. The figure inlays show the relationship between the growth rate (μ) and the GS concentration. The bacteria were cultivated at 37°C in CS-Glc minimal medium.
1.2. Genomic adaptation of Bacillus subtilis to glyphosate
Previously, it has been reported that the growth rate of a B. subtilis culture containing 2.5 mM glyphosate strongly increases after 16 h of incubation (1). It has been suggested that this phenomenon might be either due to a physiological adaptation of the bacteria to glyphosate or due to the emergence of glyphosate-resistant mutants. To assess whether B. subtilis is able to adapt to toxic levels of glyphosate at the genome level, we propagated the wild type strain 168 on CS-Glc agar plates supplemented with increasing amounts of the herbicide (0, 10, 20 30 and 40 mM). After incubation for 24 h at 37°C, bacterial growth was only visible on the plates lacking glyphosate. When the plates were further incubated for 48 h, several glyphosate-resistant mutants appeared on the plates supplemented with 10 mM glyphosate (Figure 3A). No mutants appeared on the plates with glyphosate concentrations higher than 10 mM. Next, we repeated the plating experiment with a mutant (designated as iGEM1) that was isolated from the 10 mM glyphosate plate. After 48 h of incubation, the glyphosate-resistant strain iGEM1 formed mutants on plates containing 30 mM and 40 mM glyphosate (Fig. 2A). Again we isolated mutants from the 30 mM and the 40 mM glyphosate plates that were designated as iGEM 7 and iGEM13, respectively (Figure 3B). Our attempts to further increase the glyphosate resistance of the strains iGEM7 and iGEM13 failed. To ensure that the acquired mutations in the strains iGEM1, iGEM7 and iGEM13 are stable, the mutants were passaged three times on SP rich medium plates. The glyphosate resistance of the mutants was confirmed by streaking the bacteria on CS-Glc medium supplemented with increasing amounts of glyphosate. As in Figure 3C, all strains that were adapted to glyphosate grew on plates with 10 mM glyphosate. The strains iGEM7 and iGEM13 also formed single colonies on 20 mM glyphosate plates. Moreover, albeit weak, the strain iGEM13 showed even growth with 40 mM glyphosate.
Figure 3. Isolation and characterization of glyphosate (GS)-resistant mutants. (A) Emergence of GS-resistant B. subtilis mutants on CS-Glc minimal medium agar plates that were incubated for 48 h at 37°C. WT, B. subtilis wild type strain 168; iGEM1, GS-resistant B. subtilis strain that was isolated from plates containing 10 mM GS. (B) Genealogy of the GS-resistant B. subtilis mutants. (C) Evaluation of growth of the GS-resistant mutants on CS-Glc minimal medium plates supplemented with increasing amounts of GS. The plates were incubated for 48 h at 37°C.
Our growth experiments using CS-Glc liquid medium with increasing amounts of glyphosate also confirmed that the strains iGEM7 and iGEM13 tolerate higher amounts of the herbicide than their progenitor strain iGEM1 (Figure 4; see below). To conclude, B. subtilis rapidly adapts to toxic glyphosate concentrations at the genome level!
Figure 4. Characterization of glyphosate (GS)-resistant B. subtilis mutants. The wild type (WT) strain (A) and its derivatives iGEM1 (B), iGEM7 (C), and iGEM13 (D) were cultivated for 10 h at 37°C in CS-Glc minimal medium supplemented with increasing amounts of GS.
1.3. Identification of mutations by genome and Sanger sequencing
To identify the underlying genomic alterations causing glyphosate resistance the strains iGEM1, iGEM7 and iGEM13, we performed whole-genome sequencing with the help of the Göttingen Genomics Laboratory (Figure 5A). By analysing the genome sequences we found that the bacteria had deleted a 1240 bp-long fragment containing a part of the 5’ ends of the yhfE and gltT genes and the entire yhfF gene (Figure 5B).
Figure 5. The mutations of the glyphosate-resistant mutants and the link between glutamate and glyphosate transport. (A) The illumina sequencing reads were mapped onto the reference the B. subtilis 168 reference genome NC_000964 from GenBank (3) using the Geneious software package (Biomatters Ltd.) (4). (B) The mutations in the glyphosate (GS)-resistant B. subtilis mutants identified by genome sequencing. (C) Effect of the identified mutations on the glutamate transporters GltT and GltP in the strains iGEM1, 7 and 13. Membrane topologies of the transporters were illustrated using the web-based protein topology tool Protter (5). (D) Comparison of glutamate and glyphosate.
The 1240 bp-long deletion is also present in the strains iGEM7 and iGEM13 because they were derived from the strain iGEM1 (Figure 5B). In addition to the deletion affecting the yhfE, yhfF and gltT genes, the strains iGEM7 and iGEM13 had replaced a single nucleotide and deleted 11 nucleotides, respectively, in the gltP gene. The yhfE gene encodes a putative aminopeptidase and the function of the yhfF gene is unknown (SubtiWiki; 6). The gltT and gltP genes code for the high-affinity sodium-coupled glutamate/aspartate symporter GltT and for the low-affinity proton/glutamate symporter GltP, respectively (7,8). The fact that all strains had acquired mutations in genes encoding transporters that are involved in glutamate and aspartate uptake suggests that the proteins also transport glyphosate (Figure 5B). The 1240 bp- and 11 bp-long deletions in the gltT and gltP genes truncate the encoded proteins after 282 and 214 amino acids, respectively (Figure 5C). The G232A mutation in the gltP gene results in a replacement of a glutamate residue by a lysine residue at position 78, which is adjacent to the second transmembrane helix of the GltP transporter (Figure 5C). While the deletions certainly destroy the amino acid transporters, the E78K replacement probably reduces the transport of glyphosate via GltP. This idea is in line with the observation that the strain iGEM7 carrying the G232A mutation in the gltP gene did not grow as good as the strain iGEM13 with the deletion in gltP on plates with high amounts of glyphosate (Figure 5C). A structural comparison revealed that glutamate and glyphosate resemble each other because the molecules contain negatively charged groups encompassing positively charged groups (Figure 5D). Thus, both, GltT and GltP could be involved in glyphosate uptake. To conclude, B. subtilis rapidly adapts to toxic glyphosate concentrations by inactivating the gltT and gltP genes!
1.4. Glyphosate-resistant mutants preferentially inactivate the gltT gene
Previously, it has been shown that bacteria like Salmonella typhimurium can develop glyphosate resistance due to the acquisition of mutations in the promoter region and in the coding region of the aroA gene, elevating the EPSP synthase levels and rendering the encoded enzyme insensitive to glyphosate, respectively (9-11). Our genome sequencing analyses of the isolated iGEM1, iGEM7 and iGEM13 strains did not identify mutations that alter the synthesis or the coding sequence of the EPSP synthase. To assess whether B. subtilis can develop glyphosate resistance by the acquisition of mutations that alter the target of the herbicide, we isolated 12 additional mutants from CS-Glc plates containing 10 mM glyphosate. DNA sequencing analysis revealed that each strain had mutations within the gltT coding sequence that would truncate or lengthen GltT transporter (Figure 6). Thus, the inactivation of the gltT gene is the predominant mutational event allowing B. subtilis to grow in the presence of the weedkiller glyphosate.
Figure 6. The mutations in the gltT gene and the consequences for the encoded transporter GltT in the B. subtilis mutants that were isolated from CS-Glc medium agar plates containing 10 mM glyphosate.
1.5. Characterization of gltT, gltP and gltT gltP mutant strains
Next, we tested whether the clean deletion of either the gltT gene or the gltT gltP genes is indeed sufficient to confer high-level resistance of B. subtilis to glyphosate. For this purpose, we constructed the mutant strains BP233 (gltT), BP234 (gltP) and BP235 (gltT gltP). To assess the glyphosate resistance of the gltT, gltP and gltT gltP mutants, we cultivated the bacteria in CS-Glc minimal medium supplemented with increasing amounts of glyphosate. As shown in Figure 7A and 7C, growth of the wild type and the gltP mutant strains, respectively, was inhibited by 5 mM glyphosate. This observation indicates that GltP is a low-affinity glyphosate transporter and that the cell still contains a transport system for the herbicide. The transporter with a high-affinity for glyphosate seems to be in fact GltT because the deletion of the gltT gene conferred high-level resistance to the herbicide (Figure 7B). Moreover, the deletion of the gltP gene in the background of a gltT mutant strain further enhanced glyphosate tolerance and the growth rates of the strains BP233 (gltT) and BP235 (gltT gltP) were reduced by 50% at herbicide concentrations of 6.1 mM and 8.0 mM, respectively (Figure 7B and 7D). Thus, in comparison to the wild type strain, 6- to 8-fold higher glyphosate concentrations are needed to reduce the growth rate of the strains BP233 (gltT) and BP235 (gltT gltP) by 50% (Figure 8).
Figure 7. Inactivation of the glutamate transporters confers high-level resistance to glyphosate (GS). (A) Growth of the B. subtilis wild type (WT) strain 168 in CS-Glc minimal medium supplemented with increasing amounts of GS. (B) Growth of the ΔgltT mutant strain BP233 in CS-Glc minimal medium supplemented with increasing amounts of GS. The figure inlay shows the relationship between the growth rate (µ) and the GS concentration. (C) Growth of the ΔgltP mutant strain BP234 in CS-Glc minimal medium supplemented with increasing amounts of GS. (D) Growth of the ΔgltT ΔgltP mutant strain BP235 in CS-Glc minimal medium supplemented with increasing amounts of GS. The figure inlay shows the relationship between the growth rate (µ) and the GS concentration.
Figure 8. The relationship between the growth rate (µ) and the glyphosate (GS) concentration for the B. subtilis wild type (WT) strain 168, the ΔgltT mutant strain BP233 and the the ΔgltT ΔgltP mutant strain BP235 in CS-Glc minimal medium supplemented with increasing amounts of GS.
Figure 9. A. Growth of the strains 168 (wild type), BP233 (ΔgltT), BP234 (ΔgltP) and BP235 (ΔgltT ΔgltP) on CS-Glc minimal medium in the presence of different amounts of glyphosate (GS). B. Roundup® Alphee purchased from Westfalia. C. Growth of the strains 168, BP233, BP234 and BP235 on CS-Glc minimal medium in the presence of different amounts of Roundup® Alphee.
1.6. Complementation of the gltT mutation
We also performed a complementation experiment to provide further evidence that GltT is the major glyphosate transporter in B. subtilis. For this purpose, we fused the artificial Palf4 promoter (12) together with the gapA ribosome-binding site to the gltT gene and integrated the construct into the amyE locus of the gltT mutant (Figure 10A). Next, we propagated the strains BP233 (gltT) and BP237 (gltT Palf4-gltT) together with the wild type strain 168 on agar plates without and with glyphosate (Figure 10B). As expected, all strains grew in the absence of glyphosate. By contrast, only the gltT mutant strain BP233 grew with glyphosate. Thus, like the wild type also the complementation strain BP237 take up glyphosate via GltT. To conclude, under the tested growth conditions the high-affinity glutamate transporter GltT is the major entryway of glyphosate into B. subtilis!.
Figure 10. Complementation of the ΔgltT mutant strain. (A) The artificial Palf4 promoter (12) and the ribosome-binding site (RBS) of the B. subtilis gapA gene were fused by PCR to the gltT gene. (B) Growth of the B. subtilis wild type 168 (WT), the ΔgltT mutant BP233 and the ΔgltT Palf4-gltT complementation strain BP237 on CS-Glc minimal medium agar plates without glyphosate (0 mM GS) and with glyphosate (10 mM GS). The plates were incubated for 24 h at 37°C.
1.7. GltT is also involved in the uptake of the herbicide glufosinate
Our structural comparison of glutamate and glyphosate revealed that the two molecules resemble each other because they contain negatively charged groups encompassing positively charged groups (see Figure 5D). However, the structural similarity of glutamate and glyphosate is not very high. During a scientific discussion with our collaborator Dr. Till Ischebeck from the Department of Plant Biochemistry at the University of Göttingen, we came up with the idea that GltT may also transport the herbicide glufosinate, which is structurally more similar to glutamate than glyphosate and inhibits the glutamine synthetase (Figure 11A) (13). Glufosinate (phosphinothricin) is a naturally occurring broad-spectrum systemic herbicide produced by several species of Streptomyces soil bacteria. Plants may also metabolize bialaphos (L-Alanyl-L-alanyl-phosphinothricin), another naturally occurring herbicide, directly into glufosinate (13). The compound irreversibly inhibits glutamine synthetase, an enzyme necessary for the de novo synthesis of glutamine from ammonia and glutamate, giving it antibacterial, antifungal and herbicidal properties. Application of glufosinate to plants leads to reduced glutamine and elevated ammonia levels in tissues, halting photosynthesis, resulting in plant death. In the 1960s and early 1970s, scientists at University of Tübingen and at the Meiji Seika Kaisha Company independently discovered that species of Streptomyces bacteria produce a tripeptide they called bialaphos that inhibits bacteria; it consists of two alanine residues and a unique amino acid that is an analog of glutamate that they named "phosphinothricin". The inactivation of the gltT gene indeed allowed B. subtilis to tolerate glufosinate (Figures 11B and 11C)! Moreover, the observation that the expression of two copies of the gltT gene in strain BP236 (amyE::Palf4-gltT) increased the sensitivity of the bacteria toward glufosinate further supports the idea that GltT is involved in glufosinate uptake by B. subtilis (Figures 11B and 11C). In contrast to glyphosate, however, the simultaneous deletion of the gltT and gltP genes did not increase glufosinate tolerance of B. subtilis (Figures 11B and 11C). To conclude, under the tested growth conditions the high-affinity glutamate transporter GltT is the major entryway of glyphosate into B. subtilis and GltT is also involved in Glufosinate uptake!
Figure 11. Glufosinate is also taken up by B. subtilis via GltT. (A) Comparison of glutamate and glufosinate. (B) Disk diffusion assay to test whether glufosinate is transported via GltT and/or GltP. The strains 168 (wild type), BP233 (ΔgltT), BP234 (ΔgltP, BP235 (ΔgltT ΔgltP), BP237 (ΔgltT amyE::Palf4-gltT) and BP236 (amyE::Palf4-gltT) were cultivated in LB medium, washed once in 1 X C-salts and propagated (100 µl cell suspension, OD600 ≈ 2) on CS-Glc minimal medium agar plates. A filter disk containing 5 µl of a 0.5 M glufosinate solution was placed in the center of the plates, which were incubated for 2 days at 37°C. (C) Evaluation of the disk diffusion assay shown in (B). (D) Killing two birds with one stone: the ΔgltT mutation confers resistance to glyphosate and glufosinate!
1.8. A novel nechanism conferring resistance to glyphosate
In the past years, the underlying molecular mechanisms conferring resistance to glyphosate have been intensively studied in bacteria and plants. In many cases glyphosate resistance is directly linked to the target of the herbicide. For instance, bacteria like S. aureus are a priori resistant to glyphosate because this organism synthesizes an insensitive EPSP synthase (14-17). Moreover, increased cellular levels of the EPSP synthase, which can be achieved either by overexpression of the coding gene or by gene amplification, can confer resistance to glyphosate (Figure 12) (9, 18-22). Due to the increased cellular levels of the EPSP synthase, the glyphosate is probably titrated away and sufficient amounts of the precursor for aromatic amino acid biosynthesis can be produced. As mentioned in the design section, several glyphosate-insensitive EPSP synthase mutant variants have been isolated and engineered (11, 23-36). Often, a single amino acid exchange is sufficient to render the EPSP synthase insensitive to glyphosate. However, although glyphosate resistance is frequently linked to the target of the herbicide, resistance against the herbicide may occur by other means. Several studies have demonstrated that glyphosate can be detoxified by covalent modification. For instance, the Gat glyphosate N-acetyltransferase from Bacillus licheniformis, which was subjected to directed evolution for creating an enzyme with higher efficiency and increased specificity for the herbicide, converts glyphosate to N-acetylglyphosate, which is not herbicidal and is not an effective inhibitor of EPSP synthases (Figure 12) (37-40). Moreover, the hygromycin phosphotransferases Hph and GlpA from E. coli and Pseudomonas pseudomallei, respectively, phosphorylate glyphosate and thus confer tolerance to the herbicide (Figure 12) (41,42). Interestingly, the gat gene has also been used as a selection marker for genetic engineering of bacteria (40). The enzymes that covalently modify glyphosate have been successfully introduced into crops to increase herbicide resistance (43). Recently, it has been also demonstrated that enhanced export of glyphosate can reduce toxicity of the herbicide. For instance, the overexpression of the uncharacterized membrane proteins MFS40 and YhhS from Aspergillus oryzae and E. coli, respectively, with similarity to the major facilitator secondary transporter superfamily, enhance glyphosate tolerance of E. coli (Figure 12) (44,45). It will be interesting to test whether these transporters are suitable for engineering crops to enhance glyphosate tolerance. Proteins of unknown function can also increase glyphosate tolerance. For instance, the overexpression of the igrA gene product from the Pseudomonas sp. strain PG2982 increases glyphosate tolerance in transgenic rice (46-48). Finally, many bacteria can survive in the presence of glyphosate because they are able to degrade the herbicide (Figure 12) (49-58). To conclude, several mechanisms of glyphosate resistance have been described over the past years and novel mechanisms allowing survival with the herbicide are certainly identified in the future. We have observed that the deletion of the gltT gene in B. subtilis confers high-level glyphosate resistance to the bacteria. The B. subtilis gltT mutant strain could be a suitable host to identify glyphosate uptake systems from plants by expressing cDNA libraries and by screening for transformants that are unable to grow with glyphosate.
Figure 12. Mechanisms conferring resistance to glyphosate. GltT and GltP, B. subtilis glutamate transporters; YhhS and MFS40, major facilitator secondary transporters from E. coli and A. oryzae, respectively; GAT, B. licheniformis glyphosate N-acetyltransferase; GlpA, P. pseudomallei hygromycin phosphotransferase; PEP, phosphoenolpyruvate; 3PS, 3-phosphoshikimate; EPSP, enolpyruvylshikimate-3-phosphate; EPSPS, enolpyruvylshikimate-3-phosphate synthase. Membrane topologies of the transporters were illustrated using the web-based protein topology tool Protter (5).
2. Identification of glyphosate-regulated promoters
2.1. Effect of glyphosate on the Bacillus subtilis transcriptome
The identification of genes that are up- or down-regulated in a glyphosate-dependent manner requires the application of "omics" techniques such as proteomics or transcriptomics. RNA-sequencing (RNA-Seq) analysis is a valid methodology to identify genes that are either induced or repressed by chemical substances such as nutrients and herbicides. To assess the effect of glyphosate on the B. subtilis transcriptome, we performed an RNA-Seq analysis with the help of the Göttingen Genomics Laboratory. The effect of glyphosate on the Escherichia coli transcriptome has been analyzed in the past years (59). By contrast, the effect of the weedkiller on gene expression in B. subtilis has remained elusive. To obtain bacterial cell pellets for RNA isolation and subsequent RNA-Seq analysis, we cultivated the B. subtilis wild type strain SP1 in CS-Glc minimal medium without glyphosate and supplemented with 0.75 mM glyphosate. We selected a concentration of 0.75 mM because it was previously shown that at a glyphosate concentration of about 1 mM the growth rate of the strain SP1 was reduced by 50% (see Figure 2). By choosing a lower concentration, we wanted to ensure that the bacteria were under stress without affecting growth too severe.
Figure 13. Effect of glyphosate on the B. subtilis transcriptome. (A) The diagram shows the fold change of the B. subtilis RNAs during cultivation in presence of glyphosate. The strongest regulated operon is marked with yellow circles. (B) A list of genes whose expression was severely affected by glyphosate. The four most down- and up-regulated genes (marked in (A)).
The RNA-Seq analysis allowed us to identify genes whose expression is altered in a glyphosate-dependent manner. As illustrated in Figure 13A, the expression of only a few genes was affected by glyphosate. A 5-fold increase in gene expression was the highest increase we could observe. By contrast, we observed that the expression of some genes was 11-fold reduced. Nearly all up-regulated genes are involved in sporulation and the yciA-yciB-yczL-yciC operon whose expression is affected by glyphosate is related to zinc deprivation responses (60). Unexpectedly, we could not detect any effect of glyphosate on the expression of genes encoding proteins involved in shikimate metabolism and aromatic amino acid biosynthesis. This finding could indicate that the inhibition of the EPSP synthase and the changes of the cellular EPSP concentration does not influence the expression of the genes that are involved in aromatic amino acid biosynthesis in B. subtilis. The promoters involved in the regulation of the yciA-yciB-yczL-yciC operon were identified as PyciAB and PyciC (61).
3. A reporter-based glyphosate detection system
3.1 Utilising PyciAB and PyciC to build a reporter-based glyphosate detection system
Our RNA-Seq analysis revealed that the expression of the operon containing the yciA-yciB-yczL-yciCgenes was severely affected when B. subtilis was cultivated in the presence of glyphosate. The promoters involved in the regulation of these genes were identified as PyciAB and PyciC. Based on our RNA-Seq data, we aimed to construct a reporter-based glyphosate detection system. For this purpose, we fused the identified promoters to the xylR gene, which is a transcriptional repressor of the xylA-xylB and xynP-xynB operons in B. subtilis (Figure 14A). A B. subtilis mutant strain lacking the xylR gene was transformed with the plasmid pAC7 containing our promoter-gene constructs. We expect that upon contact of the constructed strains glyphosate, the expression of the xylR gene will be lower, which will lead to the expression of the xylA reporter gene. The expression of the xylA reporter gene can be monitored by performing a β-xylosidase assay using the colourimetric substrate 4-Nitrophenyl-β-D-xylopyranoside (PNPX). The measured β-xylosidase activity can be directly linked to the expression of the xylR gene, which controlled by the glyphosate-dependent promotor (Figure 14B). The β-xylosidase activity should be directly proportional to the glyphosate that is present in the growth medium.
Figure 14. (A) Sequence of the promotores PyciAB and PyciC including the ribosome-binding site (RBS) were fused to the xylR gene by PCR. (B) Mode of action of the glyphosate detection system.
By performing a β-xylosidase assay we could determine that this system does not react to the presence of glyphosate. This might be due to the fact, that the promotores are not regulated or the regulation is not strong enough.
4. Competition assay for glyphosate detection
4.1. Competition assay with fluorophores
Unfortunately, current analytical methods, which are based on enzyme-linked immunosorbent assays (ELISAs) and mass spectrometry (MS) are very expansive, slow, need sophisticated equipment and well-trained technical staff. Therefore, the development of a rapid and straightforward method to detect glyphosate is very attractive because the intensive use of the herbicide may lead to its accumulation in the soil and food. The basic idea of the competition assay for the detection of glyphosate is depicted in Figure 15.
Figure 15. The basic principle of the competition assay for glyphosate detection. In the absence of glyphosate (no selective pressure), both strains grow with the same rate. When glyphosate is present, the herbicide-resistant strain (expressing an orange fluorophore gene) will outcompete the wild type strain. Thus, the velocity at which the wild type strain is eliminated from the co-culture is indicative for the glyphosate present in the culture or test sample.
Two B. subtilis strains that tolerate different amounts of glyphosate are labelled with different fluorophore genes (Figure 16). The fluorescence labelling of the bacteria allows to count the amount of each strain when co-cultivated in the absence and in the presence of selective pressure that is exerted by glyphosate using fluorescence microscopy. In the absence of selective pressure, we expect that both strains will survive the co-cultivation. By contrast, in the presence of glyphosate, the strain that tolerates glyphosate will outcompete the glyphosate-sensitive B. subtilis strain. Thus, the selective pressure (≈ [glyphosate]) correlates with the amount of glyphosate-resistant cells and glyphosate-sensitive cells.
Figure 16. Integration of the fluorophore genes whose expression are controlled by an artificial σA-dependent promoter (Plij2), which is constitutively active. The constructs are be integrated into the amyE gene of the B. subtilis chromosome in single copy. RBS, ribosome-binding site; ORF, open reading frame.
Characterization of the mOrange protein in B. subtilis
The fluorophore genes are derived from the iGEM DNA distribution kit. Previous work using the BBa_E2050 biobrick was focused on the function of the mOrange protein (an mRFP (DsRed) derivative) in the yeast S. cerevisiae. The aim of our experiments was to integrate the gene encoding the mOrange protein into the genomes of different B. subtilis strains. Unfortunately, the plasmid pSB1C3 does not contain an origin of replication for B. subtilis. Therefore, we introduced the fluorophore genes into the B. subtilis plasmid pAC7. For this purpose, we amplified the fluorophore genes by PCR, ligated the digested PCR products with the plasmid pAC7 that was cut with the same enzymes and transformed competent cells of the E. coli strain DH5α. During PCR, the fluorophore genes were fused to the artificial (Plij2 promoter, which should be constitutively active (see Figure 16). The promoter was designed in silico and contains binding sites for the housekeeping sigma factor A and a perfect RBS for B. subtilis. As shown in Figure 17, all strains synthesizing mOrange formed cell pellets with an orange color.
Figure 17. A cell pellet of the B. subtilis wild type strain 168 is not orange. The cell pellets of the strains iGEM24, iGEM28 and iGEM36 synthesizing mOrange are slightly colored in orange.
In the next step, the plasmid containing the mOrange gene was used to transform the B. subtilis wild type strain 168 and the glyphosate resistant strains BP233 (ΔgltT) and BP235 (ΔgltT ΔgltP). The activity of the mOrange fluorophore protein in the B. subtilis wild type strain 168 was further analyzed by fluorescence microscopy. Therefore, the strain was cultivated overnight at 37°C in the darkness. 0.01 ml of the culture were transferred to a microscope slide that was covered with solid 1% agarose dissolved in H2O to inhibit the movement of the cells during long-time exposure. As shown in Figure 18, cells synthesizing the mOrange fluorophore protein were highly fluorescent. Thus, the mOrange fluorophore protein, which is part of the BBa_E2050 biobrick can be produced in B. subtilis.
Figure 18. The left panel shows cells of the strain iGEM24 (wild type strain producing mOrange) that were visualized using bright field microscopy. Magnification = 40 X; exposure time = 0.85 s; no filter. The right panel shows cells of the same strain visualized by fluorescence microscopy. Magnification = 40 X; exposure time 44.33 s; DAPI filter.
To create a powerful glyphosate detection system, which is based on intraspecies competition, we need to have strains in hands having similar fitness coefficients. Strains having a different fitness will not allow to generate a stable co-culture because the strain with the higher fitness will alway outcompete the other strain, even in the absence of a selective pressure. To characterize our strains and to determine their fitness we performed different growth experiments. First, we integrated fluorophore genes into the amyE gene of the B. subtilis strains 168, (wild type), ΔgltT (BP233) and ΔgltT ΔgltP (BP235). The wild type strain 168 is very sensitive to glyphosate (see Figure 2). By contrast, the ΔgltT single mutant and ΔgltT ΔgltP double mutant lacking the glyphosate transport systems tolerate high amounts of the herbicide (see Figures 7 and 8). We used the gfp, bfp and morange genes encoding the fluorophores GFP, BFP and mOrange, respectively, for strain labeling. To drive the expression of the fluorophore genes, we attached an artificial promoter to the genes during PCR (see Figure 16). The sigma factor A-dependent artificial promoter is constitutively active and contains a ribosome binding site for B. subtilis. Next, we cultivated the constructed strains in a multi-well plate reader, to determine their fitness based on the growth rates that are achieved in the medium.
Figure 19. Growth of the strains iGEM24 (mOrange), iGEM28 (ΔgltT mOrange) and iGEM36 (ΔgltT ΔgltP mOrange) in CS-Glc minimal medium. All strains synthesizing mOrange have the same growth rate.
As shown in Figure 19, all strains had similar growth rates in CS-Glc minimal medium. Thus, the expression of the mOrange gene and the absence of the glyphosate transporters GltT and GltP does not affect the fitness of the strains!
To detect glyphosate by intraspecies competition, we mixed the wild type strain synthesizing the GFP protein with the ΔgltT ΔgltP double mutant synthesizing the mOrange fluorescence protein in a 1:1 ratio. As shown in Figure 20, the amount of wild type cells decreased in the co-culture with increasing glyphosate concentrations. By contrast, the amount of cells of the strain carrying the mOrange gene increased with increasing glyphosate concentrations since the ΔgltT ΔgltP mutant tolerates the herbicide glyphosate. Thus, the intraspecies competition assay, which is based on fluorescently labelled B. subtilis strains seems to be suitable to detect glyphosate!
Figure 20. Green and orange cells represent the wild type strain (gfp) and the double transporter mutant (ΔgltT ΔgltP mOrange), respectively. The ratio indicates the relative amount of orange cells in comparison to green cells. (A) Merged orange and green channels from culture with 0 mM glyphosate after 4 h of cultivation. There are nearly equal amounts of green and orange cells. (B) Merged channels orange and green from culture with 2 mM glyphosate after 4 h of cultivation. The number of orange cells is 4.21-fold higher. (C) Merged orange and green channels from culture with 3.5 mM glyphosate after 4 h of cultivation. Nearly no green cells were visible. (D) Calculated ratio between the number of orange and green cells with increasing glyphosate concentrations.
To perform a competition assay in a greater scale, we used a platereader for our growth experiments. This device enabled us to measure different fluorophores and the optical density simultaneously. Strains were cultivated in CS-glucose (CS-Glu) medium to an OD600 of 1–1.5. A 96–well plate was prepared containing 100 µl of CS-Glc minimal medium with the desired concentrations of glyphosate. 10 µl of the cell mixture (1:1 ratio) were transferred into the wells to obtain an OD600 of 0.1 in the wells. We used the Synergy MX II multi-well plate reader from the BioTek company for our measurements. GFP was excited at 488 nm and emission was measured at 509 nm. for mOrange, the excitation wavelength was 529 nm and emission was detected at 562 nm.
Figure 21. The fluorescence is normed to the optical density, which equals the cell number. A: Intensity of GFP fluorescence decreases with inceasing glyphosate concentration(0–2 mM glyphosate). B: Intensity of mOrange fluorescence increases with increasing glyphosate concentration 0–2 mM glyphosate).
This competition assay enables us to calculate the concentration of glyphosate in samples from people that are curious, how much glyphosate is in their soil.
4.2. Competition with β–galactosidase expressing strains
Another approach using competition assays is to use enzymatic indication systems. Therefore, we used the strain BP235 (ΔgltT ΔgltP) and the strain BP193 (palf4::lacZ) (12). As previously described, BP235 is characterized by high-level glyphosate resistance, while BP193 resembles the wild type despite the lacZ gene. The promoter palf4 is an artificial one and harbors a perfect consensus sequence for B. subtilis and a ribosome binding site. This leads to constitutive expression of the lacZ gene and to blue cells if the medium is supplemented with X-Gal. Cultures of both strains in CS-Glc minimal medium without glyphosate were mixed to obtain a 1:1 ratio of both strains and incubated with increasing glyphosate concentrations (0 - 5mM) for 24 h (Figure 22).
Figure 22. Cell mixture (1:1) of strains BP193 and BP235 were incubated with increasing glyphosate concentrations (0–5 mM). Strain BP193 harbors the lacZ gene behind an artificial promoter, but does not tolerate glyphosate in contrast to strain BP235. The intensity of the blue color decreases with increasing glyphosate concentration. Strains BP193 and BP235 were also incubated without glyphosate as a control.
The intensity of the blue color drecreases with increasing glyphosate concentrations (Figure 22). The change and intensity of the blue color correlates with the glyphosate concentration and can therefore be used as detection and quantification system. Furthermore, a β–galactosidase assay can quantify the activity of the promoter palf4 and thus, the concentration of glyphosate.
The promoter-reporter gene fusion of palf4 and lacZ allows to determine the activity of the promoter. Since palf4 is constitutively active this is directly linked to the survival of the cells. Transcription results in synthesis of ß-galactosidase which in turn hydrolyses o-nitrophenyl-ß-D-galactopyranoside (ONPG), resulting in the formation of o-nitrophenol. Formed o-nitrophenol absorbs light with a wavelength of 420 nm, and can therefore be used as an indicator for the amount of cells and subsequently for the concentration of glyphosate. We determined ß-galactosidase activity for the mixture of BP193 and BP235 in presence of increasing amounts of glyphosate. In addition, we used BP193 and BP235 as controls. Cells were grown in a 1:1 ratio. ß-galactosidase activity was high in the BP193 control and at low glyphosate concentration. Increasing glyphosate concentration is directly linked to cell death of BP193 and a corresponding decrease in ß-galactosidase activity (Figure 23).
Figure 23. Cell mixture (1:1) of strains BP193 and BP235 were incubated with increasing glyphosate concentrations (0 - 5 mM). Strain BP193 harbors the lacZ gene behind an artificial promoter, but does not tolerate glyphosate in contrast to strain BP235. ß-galactosidase activity can be measured at a wavelength of 420 nm using the substrate o-nitrophenyl-ß-D-galactopyranoside. Activity decreases with increasing glyphosate concentration. BP235 does not contain the lacZ gene, hence no activity can be observed.
In an additional step, we verified that the change in coloration occurred due to the death of the sensitive strain on the one hand and the increased growth of the resistant BP235 strain on the other hand. Hence, we adopted the experimental set up and analysed growth of the two strains in the platereader. To analyse growth of BP235, we used the strain constitutively expressing GFP, which was constructed similarly to the previously used strains. For the analysis we again used 100 µL CS-Glu medium supplemented with the appropriate glyphosate concentration. The cells were added in a 1:1 ratio. GFP was excited at 488 nm and emission was measured at 509 nm.
Figure 24. Cell mixture (1:1) of strains BP193 and BP235 were incubated with increasing glyphosate concentrations (0–5 mM). Strain BP193 harbors the lacZ gene behind an artificial promoter, but does not tolerate glyphosate in contrast to strain BP235.
5. Engineering bacteria to disarm glyphosate
5.1. The glyphosate N-acetyltransferase confers resistance to weedkiller
Many years ago it has been found that glyphosate can be inactivated by acetylation yielding N-acetylglyphosate (see Figure 25) (37). Because the acetylated form of glyphosate has low affinity for the active site of EPSP synthase, it is nonherbicidal. The glyphosate N-acetyltransferase GAT from the Gram-positive soil bacterium Bacillus licheniformis transfers the acetyl group from acetyl-CoA onto the amino group of the weedkiller (37). Interestingly, B. subtilis possesses a putative glyphosate N-acetyltransferase, which is designated as YitI and shares 59% overall amino acid identity with GAT of B. licheniformis. However, the native GAT enzymes are very poor catalysts for N-acetylation of glyphosate (37). Moreover, despite extensive screening of biological amines, including amino acids, nucleotides and antibiotics, the physiological substrates for the native enzymes are unknown (39). As an alternative strategy for glyphosate resistance involving enzymatic conversion of glyphosate to N-acetylglyphosate the GAT was also subjected to directed evolution for creating an enzyme with higher efficiency and increased specificity for the herbicide (37). When introduced into plants, optimized gat genes confer robust tolerance to glyphosate (37). Therefore, we expect that bacteria like B. subtilis expressing the optimized GAT should tolerate high amounts of glyphosate even though it is taken up from the environment via the high-affinity glyphosate transporter GltT (see above). To evaluate the potential of the GAT to allow B. subtilis in the presence of glyphosate concentrations that are toxic for the wild type bacteria, we wanted to express the gat gene in B. subtilis (Figure 25).
Figure 25. (A) Sequence alignment showing the similarities of the wild type (wt) gat alleles from Bacillus licheniformis and Bacillus subtilis as well as the gat alleles R7 and R11 from B. licheniformis that were obtained by 11 rounds of gene shuffling. The red arrows indicate the 31 amino acid exchanges in the glyphosate N-acetyltransferase (Gat) variant R11. The optimized Gat variants shown an up to a 4,500-fold increase in catalytic efficiency (kcat/Km) relative to the native enzyme. The GenBank accession numbers for the sequences: B. licheniformis ST401 GAT, AX543338; R7 GAT, AY597417; R11 GAT, AY597418; B. subtilis YitI, CAA70664. (B) Structure of the Gat R7 variant (PDBid: 2JDD) in ternary complex with acetyl-CoA and the competitive inhibitor 3-phosphoglycerate (3PG). The four active site residues (Arg-21, Arg-73, Arg-111, and His-138), which are labelled in black color, contribute to a positively charged substrate-binding site that is conserved throughout the GAT subfamily. (C) Top view on the active site of Gat R7. (D) Construction of a plasmid for the expression of the gat R11 variant in B. subtilis. Expression is driven by the constitutively active σA-dependent Palf4 promoter The plasmid pAC7::gat R11 integrates into the amyE locus of the B. subtilis chromosome in single copy. lacZ, beta-galactosidase; amyE-5' and amyE-3', fragments of the B. subtilis amyE for homologous recombination; amp, ampicillin resistance gene for the selection in E. coli; kan, kanamycin resistance gene for the selection in B. subtilis; ori pBR322, origin of replication for E. coli.
For this purpose, we ordered the gat gene as a g-Block from Integrated DNA Technologies, which was free of charge because the company supports the iGEM competition. Next, we amplified the gat gene by PCR using the g-Block as template DNA. An artificial promoter and a ribosome binding site was attached during the PCR (Figure 25D). The promoter-gene fusion was cloned and introduced in single copy into the genomes of the B. subtilis wild type strain and the gltT mutant strain to further increase its glyphosate tolerance (Figure 7). To assess the potential of the GAT to disarm glyphosate in the different B. subtilis strains, we propagated the bacteria on CS-Glc minimal medium agar plates that were supplemented with different amounts of glyphosate. As shown in Figure 26, the strains synthesizing the GAT tolerated significantly more glyphosate than the wild type strain. Moreover, the inactivation of the gltT glyphosate transporter gene and the overexpression of the gat gene confers high-level glyphosate tolerance to the bacteria. Thus, when introduced into B. subtilis, the optimized gat gene confers indeed robust tolerance to glyphosate. To conclude, the engineered bacteria possess two interesting properties: (i) B. subtilis promotes plant growth by protecting plants against pathogens and (ii), the bacteria take up glyphosate from the environment and inactivate the weedkiller!
Figure 26. (A) Growth of the wild type and ΔgltT mutant and of the isogenic strains overexpressing the gat R11 and aroA* genes on CS-Glc minimal medium in the presence of different amounts of glyphosate. (B) Growth of the strains described in A on CS-Glc minimal medium in the presence of different amounts of Roundup® Alphee.
Cultivation of the different strains in liquid media allows more detailed characterization. Therefore, we analyzed growth of the bacteria in CS-Glc minimal medium supplemented with increasing amounts of glyphosate, using a multi-well plate reader. In general, the analysis supports previous results (see Figure 26). The strain containing the aroA* gene confers only a slight increase in resistance in comparison to the wild type strain and hence a significantly lower resistance in comparison the gltT mutant strain (Figure 27A, B and D). However, for the strain containing the aroA* gene and lacking the glyphosate transporter GltT a 9-fold higher glyphosate concentration is needed to reduce the growth rate by 50% (Figure 27F). The growth experiments confirm our previous results showing that the glyphosate resistance in can be strongly enhanced in B. subtitlis by expressing the gat gene or by inactivating the gltT. Moreover, a combination of both even further enhances glyphosate resistance of the bacteria. (Figure 27C and E, Figure 26).
Figure 27. (A) Growth of the B. subtilis wild type strain in presence of 0 - 3.5 mM glyphosate. (B) Growth of the B. subtilis gltT mutant strain with increasing glyphosate concentrations (0-30 mM). (C) Growth of the B. subtilis strain producing Gat with increasing glyphosate concentrations (0-30 mM. (D) Growth of the B. subtilis strain producing AroA* with increasing glyphosate concentrations (0-30 mM). (E) Growth of the B. subtilis strain producing AroA* with increasing glyphosate concentrations (0-30 mM). Growth of the B. subtilis strain lacking GltT and producing Gat with increasing glyphosate concentrations (0-30 mM). (F) Growth of the B. subtilis strain lacking GltT and producing AroA* with increasing glyphosate concentratiosn (0-30 mM). Bacteria were cultivated in liquid CS-Glc minimal medium.
Figure 28. The relationship between the growth rate (µ) and the glyphosate (GS) concentration for the B. subtilis wild type (WT) strain SP1, the ΔgltT mutant strain BP233, the B. subtilis strain producing Gat, the B. subtilis strain producing AroA*, the B. subtilis strain lacking GltT and producing AroA*, the B. subtilis strain lacking GltT and producing Gat in CS-Glc minimal medium supplemented with increasing amounts of GS.
6. Characterization of EPSP synthases
6.1. Essentiality of the B. subtilis aroE EPSP synthase gene
To our surprise, we did not find glyphosate-resistant mutants of B. subtilis that had accumulated mutations in the aroE gene rendering the encoded EPSP synthase insensitive to the herbicide. Since AroE from B. subtilis belongs to the class II EPSP synthases, which are less sensitive to glyphosate, the potential of enzyme to further evolve herbicidal insensitivity might be limited (Figure 29) (62-64).
Figure 29. The EPSP synthase fall into three different groups. While the E. coli EPSP synthase glyphosate sensitive (group I), the B. subtilis enzyme belongs to the group II of glyphosate-insensitive EPSP synthases. Not much is known about the glyphosate-dependent inhibition of the EPSP synthases that belong to group III. The phylogenetic three has been constructed using the Geneious software package (Biomatters Ltd.). Sequences of the EPSP synthases were derived from the UniProt database.
Figure 30. The EPSP synthase is essential in B. subtilis. Transformation experiment to evaluate the essentiality of the EPSP synthase in B. subtilis. The bacteria were propagated on SP plates supplemented with spectinomycin and incubated for 24 h at 37°C.
6.2. Purification and characterization of EPSP synthases
This year, our team came toghether with the iGEM team of Marburg, to investigate the interaction of the EPSP synthases AroE and AroA from B. subtilis and E. coli, respectively, with glyphosate. Glyphosate targets AroE and AroA in B. subtilis and E. coli, respectively (see Background information). Since the members of the Bange lab, which is hosting the iGEM team of Marburg, are experts in the field of Isothermal titration calorimetry (ITC), we did not hesitate to contact them for support. Our team cloned the aroE and aroA genes of B. subtilis and E. coli, respectively, and the resulting constructs are suitable for overexpression of the Strep-tagged AroE and AroA enzymes in the E. coli strain BL21. The N-terminally Strep-tagged proteins can be purified by Streptactin:Strep-tag affinity purification system from the IBA, Göttingen. The purified proteins were send to the iGEM team in Marburg (Figure 31). The ITC measurments were performed in Marburg. Unfortunately, not interaction between glyphosate and the EPSP synthases could be detected. Unfortunately, also the second ITC measurement using freshly purified AroA and AroE enzymes from E. coli and B. subtilis, respectively, did not reveal whether glyphosate may interact with the EPSP synthases. Nevertheless, we are grateful to the iGEM Team Marburg for performing the ITC measurements.
Figure 31. (A) Evaluation of the purification of the N-terminally Strep-tagged EPSP synthases from B. subtilis and E. coli Strep-AroE and Strep-AroA, respectively, by 12% SDS PAGE. The proteins were stained with Coomassie Brilliant Blue. M, unstained protein molecular weight marker Thermo Scientific; CE, crude extract; FT, flow through, W, washing steps; E, elution steps. (B) ITC measurements with Strep-AroE and Strep-AroA and glyphosate. (C) Control experiments with water and the calcium chelator EDTA.
References
- Fischer et al. (1986) J. Bacteriol. 168: 1147-1154
- Zeigler et al. (2008) J. Bacteriol. 190: 6983-6995.
- Barbe et al. (2009) Microbiology. 155: 1758-1775
- Kearse et al. (2012) Bioinformatics. 28: 1647-1649.
- Omasits et al. (2014) Bioinformatics. 30: 884-886.
- Michna et al. (2014) Nucleic Acids Res. 42: D692–D698.
- Zaprasis et al. (2015) Appl. Environ. Microbiol. 81: 250-259.
- Tolner et al. (1995) J. Bacteriol. 177: 2863-2869.
- Amrhein et al. (1980) Naturwissenschaften. 67: 356-357.
- Comai et al. (1983) Science. 221: 370-371.
- Stalker et al. (1985) J. Biol. Chem. 260: 4724-4728.
- Gundlach et al. (2017) Sci. Signal. 10: eaal3011.
- Gundlach et al. (2017) Sci. Signal. 10: eaal3011.
- Fraser & Ridley (1984) Planta. 161: 470-474.
- Priestmann et al., (2005) FEBS Lett. 579: 728-732.
- Cao et al. (2012) PLoS One. 7: e38718.
- Light et al. (2016) Biochemistry. 55: 1239-1245.
- Liu & Cao, (2018) Biotechnol. Lett. 40: 855-864.
- Rogers et al. (1983) Appl. Envrion. Microbiol. 46: 37-43.
- Gaines et al. (2011) J. Agric. Food Chem. 59: 5886-5889.
- Juglam et al. (2014) Plant Physiol. 166: 1200-1207.
- Sammons & Gaines (2014) Pest. Manag. Sci. 70: 1367-1377.
- Dillon et al. (2017) Plant Physiol. 173: 1226-1234.
- Sost & Amrhein (1990) Arch. Biochem. Biophys. 282: 433-436.
- Padgette et al. (1991) J. Biol. Chem. 266: 22364-22369.
- He et al. (2001) Biochim. Biophys. Acta 1568: 1-6.
- He et al. (2003) Biosci. Biotechnol. Biochem. 67: 1405-1409.
- Bearson et al. (2002) Plant Physiol. 129: 1265-1275.
- Eschenburg et al. (2002) Planta 216: 129-135.
- Sun et al. (2005) Appl. Environ. Microbiol. 71: 4771-4776.
- Zhou et al. (2006) Plant Physiol. 140: 184-195.
- Healy-Fried et al. (2007) J. Biol. Chem. 282: 32949-32955.
- Vande Berg et al. (2008) Pest. Manag. Sci. 64: 340-345.
- Funke et al. (2009) J. Biol. Chem. 284: 9854-9860.
- Tian et al. (2010) Appl. Environ. Microbiol. 76: 6001-6005.
- Pollegioni et al. (2011) FEBS J. 278: 2753-2766.
- Chekan et al. (2016) MedChemComm. 7: 28-36.
- Castle et al. (2004) Science. 304: 1151-1154.
- Siehl et al. (2005) Pest. Manag. Sci. 61: 235-240.
- Siehl et al. (2007) J. Biol. Chem. 282: 11446-11455.
- Norris et al. (2009) Appl. Environ. Microbiol. 75: 6062-6075.
- Rao et al. (1983) Antimicrob. Agents Chemother. 24: 689-695.
- Penaloza-Vazquez et al. (1995) Appl. Environ. Microbiol. 61: 538-543.
- Liu et al. (2015) Transgenic Res. 24: 753-763.
- Staub et al. (2012) J. Ind. Microbiol. Biotechnol. 39: 641-647.
- Tao et al. (2017) Pestic. Biochem. Physiol. 140: 65-68.
- Fitzgibbon & Braymer, (1988) Appl. Environ. Microbiol. 54: 1886-1888.
- Fitzgibbon & Braymer, (1990) Appl. Environ. Microbiol. 56: 3382-3388.
- Fartyal et al. (2018) Front. Plant. Sci. 13: 144.
- Pipke et al. (1980) Appl. Environ. Microbiol. 53: 974-978.
- Shinabarger et al. (1986) J. Bacteriol. 168: 702-707.
- Liu et al. (1991) Appl. Environ. Microbiol. 57: 1799-1804.
- Dick & Quinn, (1995) Appl. Microbiol. Biotechnol. 43: 545-550.
- Singh & Walker, (2006) FEMS Microbiol. Rev. 30: 428-471.
- Castro et al. (2007) J. Environ. Sci. Health B. 42: 883-886.
- Hove-Jensen et al. (2014) Microbiol. Mol. Biol. Rev. 78: 176-197.
- Kryuchkova et al. (2014) Microbiol. Res. 169: 99-105.
- Sviridov et al., (2011) Biochemistry (Mosc). 76: 720-725.
- Sviridov et al. (2015) Prikl. Biokhim. Mikrobiol. 51: 183-190.
- Lu et al. (2013) Mol Biosyst. 9: 522–30.
- Gaballa & Helmann (1998) J. Bacteriol. 180: 5815–5821.
- Gaballa et al. (2002) J. Bacteriol. 184: 6508–6514.
- Barry et al. (1992) WO 92/04449.
- Schönbrunn et al. (2001) Proc. Natl. Acad. Sci. USA. 98: 1376-1380.
- Light et al. (2016) Biochemistry. 66: 1239-1245.
- Koo et al. (2017) Cell Systems. 4: 291-305.
- Jordan et al. (2002) Genome Res. 12: 962-968.