- Research note
- Open Access
Attempted use of PACE for riboswitch discovery generates three new translational theophylline riboswitch side products
BMC Research Notesvolume 11, Article number: 861 (2018)
The purpose of this project was to use an in vivo method to discover riboswitches that are activated by new ligands. We employed phage-assisted continuous evolution (PACE) to evolve new riboswitches in vivo. We started with one translational riboswitch and one transcriptional riboswitch, both of which were activated by theophylline. We used xanthine as the new target ligand during positive selection followed by negative selection using theophylline. The goal was to generate very large M13 phage populations that contained unknown mutations, some of which would result in new aptamer specificity. We discovered side products of three new theophylline translational riboswitches with different levels of protein production.
We used next generation sequencing to identify M13 phage that carried riboswitch mutations. We cloned and characterized the most abundant riboswitch mutants and discovered three variants that produce different levels of translational output while retaining their theophylline specificity. Although we were unable to demonstrate evolution of new riboswitch ligand specificity using PACE, we recommend careful design of recombinant M13 phage to avoid evolution of “cheaters” that short circuit the intended selection pressure.
An important goal of microbial metabolic engineering is the development of improved methods for the integration of host and orthogonal metabolism [1,2,3,4,5]. Our research group developed a generalized method called programmed evolution for optimizing orthogonal metabolic engineering . The key to programmed evolution is a fitness module containing a riboswitch (RS) that is regulated by the desired metabolic product. The extent to which programmed evolution can be applied is limited by the availability of RSs. RSs discovered in vitro often do not predictably control gene expression in vivo . Rational approaches to RS discovery exist but they do not facilitate exhaustive sequence space searches . We hypothesized that phage-assisted continuous evolution (PACE) would have advantages for RS discovery [9,10,11]. PACE relies on sustained phage populations and induced mutations coupled with positive and negative selection of M13 phage grown in chemostats.
The experimental design we developed was intended to be a proof-of-concept for using PACE to discover RSs with new ligand specificity (Fig. 1). Our goal was to evolve xanthine-specific transcriptional and translational RSs from existing RSs that bind theophylline. Our attempt to evolve new xanthine RSs from theophylline RSs did not succeed. The RSs accumulated many mutations after positive selection but negative selection did not occur. This report provides details of our experiments and their limitations so that others might build on our work to develop PACE for RS discovery. We also report side products in the form of thee new theophylline-specific translational RS with different translational efficiencies and induction ratios.
Our adaptation of PACE placed a RS between the gene three (gIII) promoter (PIII) and the T7 RNA polymerase (RNAP) gene so that the RS was the target of directed evolution (Additional file 1: Figure S1 in PACE Additional Information). We chose RS-C from Topp et al. as the translational RS, and RS10 from Wachsmuth et al. as the transcriptional riboswitch [12, 13]. RS-C blocks production of T7 RNA polymerase unless theophylline binds to its aptamer and alters its folding so that the RBS is available for the initiation of translation. Binding of theophylline to RS10 disrupts the formation of a stem and loop structure required for the termination of transcription. We conducted three simultaneous PACE experiments in our chemostat apparatus (see details in Additional file 2: PACE Chemostat Manual 2016). One PACE chemostat contained phages that only carried RS-C, a second chemostat contained only phages with RS10, and a third contained both RSs. The objective for positive selection was to evolve a phage variant that contained a xanthine-activated RS. The experiment began with a 24 h gradual increase in positive selection for RSs that respond to xanthine instead of theophylline. The inclusion of anhydrotetracycline (aTc) in the experiment is for T7 RNAP-driven protein three (pIII) production initially until xanthine-responsive RSs arise and accumulate, after which aTc levels were gradually reduced (see Fig. 1 and Additional file 1: Appendix A in PACE Additional Information). During positive selection, we used arabinose to induce a host cell mutagenesis gene expression cassette that increases the rate of DNA mutations among phage progeny. We successfully measured green fluorescent protein (GFP) fluorescence as an indirect measure of pIII production, and used polymerase chain reaction (PCR) to confirm the persistence of M13 genomes in the chemostats (see Additional file 1: Methods in PACE Additional Information).
The objective of negative selection was to eliminate phages with RSs that are constitutively activated and do not require ligand binding to achieve an ON state and to eliminate phages with RSs that still have theophylline responsiveness. After a short transitional period to purge xanthine from the experiment, we began 16 h of negative selection with a new host strain that lacked the mutagenesis cassette to prevent evolution away from the desired RS outcome. We provided a steady supply of theophylline but gradually reduced the amount of aTc (see Fig. 1 and Additional file 1: Appendix A in PACE Additional Information). Throughout negative selection, we were able to confirm the presence of M13 DNA by PCR, but we failed to detect red fluorescent protein (RFP) as an indirect measure of the gIII neg expression (see Fig. 1 and Additional file 1: Figure S1 in PACE Additional Information). Despite our inability to confirm negative selection, we isolated phage from all three evolving chemostat populations and used PCR to amplify the RS regions of M13.
To determine the level of sequence variation among RSs carried by phage populations after PACE, we used paired-end MiSeq to sequence 3.9 million RSs from the three chemostats. Our sequence analysis revealed a high frequency of two types of RSs that were present at the ends of our three PACE experiments (see Additional file 1: Appendix B in PACE Additional Information). From the two PACE lagoons with RS10, we isolated RS10Δ7, which has a seven nucleotide deletion that disrupts the stem and loop structure required for transcriptional termination. We verified that RS10Δ7 is constitutively activated, which explains how it survived positive selection, but not how it escaped negative selection (see Additional file 1: Figure S2 in PACE Additional Information). In the RS-C experiments, we found multiple RSs that should not have persisted through either positive or negative selection (see Additional file 1: Table S2 in PACE Additional Information). The failure of negative selection was consistent with our inability to detect RFP.
The inability of our PACE experiments to exert positive selection for xanthine riboswitches and negative selection against theophylline-responsive RSs might be explained by an inherent experimental design flaw. There are more ways to evolve a RS that remains in an ON state than to evolve a RS with new ligand specificity. Our evolving populations accumulated “cheaters” during positive and negative selection. Our suggestion for others who attempt PACE for RS evolution is to avoid an experimental strategy in which disabling the DNA of interest is adaptive for the M13 phage under selection. A better strategy would enable simultaneous negative selection against RS OFF states and positive selection for ON states that are dependent on a desired ligand specificity. Positive and negative selection might need to be run longer since the RS target sequence is small compared to the entire M13 phage genome. It might take longer to accumulate sufficient mutations to alter ligand binding.
Despite our inability to evolve xanthine-specific RS, we discovered side products in the form of three new theophylline-specific RSs with three levels of ON state translational output (Fig. 2). These new RSs retained their original ligand specificity and shared similar mutations in and near their RBSs (Fig. 3). Each of the three new RS carried the same single base change plus unique insertions of 4 bases adjacent to the RBS. To explain our experimentally observed differences in translational output, we used mFold to predict secondary structures of the new RSs and compared them to original RS-C. Each of the three new RSs adopts a stem and loop structure of similar stability to the one found in RS-C, which supports the measurement of their OFF states. To assess ON states, we analyzed the ability of each of the RS RBSs to base pair with the 3′ end of 16S rRNA. The new RS J100377 can form 6 base pairs, whereas the other three RSs form only 5. These results are consistent with our observation of higher expression from the J100377 RS than the other RSs.
We acknowledge several limitations to our contribution to the synthetic biology literature of three new theophylline-specific translational RSs and of a protocol for using PACE to evolve RSs. Positive selection during our PACE experiments worked to some degree, but was confounded by the fitness advantage of a “cheater” with an easily-produced deletion that enabled constitutive pIII production. Negative selection failed during our PACE experiments, perhaps because of undetected mutations in gIII neg. We acknowledge that it might not be possible to evolve theophylline RSs into xanthine RSs with PACE. This shortcoming could be addressed by using other ligands or other RS starting points. The paper from which we chose RS-C includes several RSs that show higher levels of reporter gene expression than J100377, J100381, and J100382, and one with a larger induction ratio . Translational RSs with similar characteristics can also be found in other publications [14,15,16]. We have not measured the theophylline dose–response or explored their modularity with regard to various promoter strengths and reporter genes.
phage-assisted continuous evolution
promoter to gene three (gIII)
green fluorescent protein
RS10 with a 7 base pair deletion
polymerase chain reaction
red fluorescent protein
Sinha J, Topp S, Gallivan JP. From SELEX to cell dual selections for synthetic riboswitches. Methods Enzymol. 2011;497:207–20.
Lynch SA, Topp S, Gallivan JP. High-throughput screens to discover synthetic riboswitches. Methods Mol Biol Clifton NJ. 2009;540:321–33.
Meyer S, Chappell J, Sankar S, Chew R, Lucks JB. Improving fold activation of small transcription activating RNAs (STARs) with rational RNA engineering strategies. Biotechnol Bioeng. 2016;113:216–25.
Boerneke MA, Hermann T. Ligand-responsive RNA mechanical switches. RNA Biol. 2015;12:780–6.
Robinson CJ, Vincent HA, Wu M-C, Lowe PT, Dunstan MS, Leys D, et al. Modular riboswitch toolsets for synthetic genetic control in diverse bacterial species. J Am Chem Soc. 2014;136(30):10615–24.
Eckdahl TT, Campbell AM, Heyer LJ, Poet JL, Blauch DN, Snyder NL, et al. Programmed evolution for optimization of orthogonal metabolic output in bacteria. PLoS ONE. 2015;10(2):e0118322.
Lakhin AV, Tarantul VZ, Gening LV. Aptamers: problems, solutions and prospects. Acta Nat. 2013;5(4):34–43.
Townshed B, Kennedy AB, Xiang JS, Smolke CD. High-throughput cellular RNA device engineering. Nat Methods. 2015;12(10):989–94.
Hallberg ZF, Su Y, Kitto RZ, Hammond MC. Engineering and in vivo applications of riboswitches. Annu Rev Biochem. 2017;20(86):515–39.
Denby CM, Li RA, Vu VT, Costello Z, Lin W, Chan LJG, et al. Industrial brewing yeast engineered for the production of primary flavor determinants in hopped beer. Nat Commun. 2018;9(1):965.
Galanie S, Thodey K, Trenchard IJ, Interrante MF, Smolke CD. Complete biosynthesis of opioids in yeast. Science. 2015;349(6252):7.
Topp S, Reynoso CMK, Seeliger JC, Goldlust IS, Desai SK, Murat D, et al. Synthetic riboswitches that induce gene expression in diverse bacterial species. Appl Environ Microbiol. 2010;76(23):7881–4.
Wachsmuth M, Findeiss S, Weissheimer N, Stadler PF, Morl M. De novo design of a synthetic riboswitch that regulates transcription termination. Nucleic Acids Res. 2013;41(4):2541–51.
Mishler DM, Gallivan JP. A family of synthetic riboswitches adopts a kinetic trapping mechanism. Nucleic Acids Res. 2014;42(10):6753–61.
Felletti M, Hartig JS. Ligand-dependent ribozymes: ligand-dependent ribozymes. Wiley Interdiscip Rev RNA. 2017;8(2):e1395.
Page K, Shaffer J, Lin S, Zhang M, Liu JM. Engineering riboswitches in vivo using dual genetic selection and fluorescence-activated cell sorting. ACS Synth Biol. 2018;7(9):2000–6.
ZMS tested all candidate riboswitches and helped write the manuscript. SSB isolated and sequenced the RS variants from the evolving populations. IGC wrote the python script that identified the most abundant RS variants. SRB, AB, IGC, AJD, LD, FCE, SH, JCH, AMH, MH, HJ, OK, CCM, DyM, HS, AT, AU, CW, contributed to the building of DNA constructs, fabrication and testing of the chemostats and performing PACE experiments. LJH, JLP, MB, AE, DeM, HTM, AJS, DS, AT, and DZ developed mathematical models and web tools for use with PACE experiments. LJH, JLP, TTE and AMC mentored the undergraduates and contributed to all aspects of the project. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
Sequences for the three new riboswitches are available in the Registry for Standard Biological Parts (http://parts.igem.org/Main_Page) as Part Numbers J100377, J100381 and J100382.
Consent for publication
Ethics approval and consent to participate
NSF RUI Grant MCB-1613203 to Davidson College and MCB-1613281 to Missouri Western State University, Davidson College Martin Genomics Program. These funders played no role in the design of the study nor collection, analysis, and interpretation of data nor in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.