Skip to content

Advertisement

  • Research note
  • Open Access

Attempted use of PACE for riboswitch discovery generates three new translational theophylline riboswitch side products

  • 1,
  • 1,
  • 2,
  • 3,
  • 1,
  • 1,
  • 2,
  • 2,
  • 1,
  • 3,
  • 2,
  • 1,
  • 2,
  • 1,
  • 1,
  • 1,
  • 2,
  • 1, 4,
  • 3,
  • 3,
  • 1, 4,
  • 3,
  • 3,
  • 4,
  • 2,
  • 2,
  • 3,
  • 4,
  • 3,
  • 2 and
  • 1Email authorView ORCID ID profile
BMC Research Notes201811:861

https://doi.org/10.1186/s13104-018-3965-6

  • Received: 24 October 2018
  • Accepted: 29 November 2018
  • Published:

Abstract

Objective

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.

Results

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.

Keywords

  • Phage-assisted continuous evolution (PACE)
  • Riboswitch
  • Synthetic biology
  • Metabolic engineering
  • Theophylline
  • Xanthine
  • M13

Introduction

An important goal of microbial metabolic engineering is the development of improved methods for the integration of host and orthogonal metabolism [15]. Our research group developed a generalized method called programmed evolution for optimizing orthogonal metabolic engineering [6]. 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 [7]. Rational approaches to RS discovery exist but they do not facilitate exhaustive sequence space searches [8]. We hypothesized that phage-assisted continuous evolution (PACE) would have advantages for RS discovery [911]. 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.
Fig. 1
Fig. 1

PACE experimental design for the discovery of new RSs. a During positive selection, anhydrotetracycline (aTc) concentration is reduced while the new target ligand xanthine concentration remained constant. Mutagenesis cassette is induced using arabinose. The chemostat replenishes growth media and uninfected E. coli. b Negative selection employed a dominant gIII negative allele (gIII neg) induced by a steady concentration of theophylline to select against theophylline-specific riboswitches while gradually reducing aTc to maintain phage with xanthine-specific riboswitches

Main text

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.
Fig. 2
Fig. 2

Directed evolution of three new riboswitches. a Each new riboswitch was cloned into rClone Red and grown in triplicate over 10 h with one of three treatments as indicated. Each population was measured for RFP and cell density three times to produce an average end point value ± s.e.m. b Theophylline-induction values from three independent populations from a were normalized by comparing each construct to the fold induction of RS-C in the absence of theophylline. Error bars represent s.e.m

Fig. 3
Fig. 3

Predicted 2D structures of RS-C and three new riboswitches. a The original RS-C sequence (left) was submitted to mFold to produce a predicted 2D structure and quantification of stability (ΔG); the boxed region is the predicted RBS. The three new riboswitch sequences were also evaluated by mFold with RBS denoted by boxed nucleotides. Yellow highlighting shows inserted bases and the blue A within a box indicates a base substitution conserved in all three new riboswitches. b RNA base pairing (dashed lines) between ribosome (red) and the complementary anti-parallel  riboswitch sequence (black) for the original RS-C as well as the three new RSs. The number of base pairs between rRNA and RS sequence is shown on the far right

Limitations

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 [12]. Translational RSs with similar characteristics can also be found in other publications [1416]. We have not measured the theophylline dose–response or explored their modularity with regard to various promoter strengths and reporter genes.

Notes

Abbreviations

RS: 

riboswitch

PACE: 

phage-assisted continuous evolution

PIII: 

promoter to gene three (gIII)

gIII: 

gene three

pIII: 

protein three

RNAP: 

RNA polymerase

aTc: 

anhydrotetracycline

GFP: 

green fluorescent protein

RS10Δ7: 

RS10 with a 7 base pair deletion

PCR: 

polymerase chain reaction

RFP: 

red fluorescent protein

Declarations

Authors’ contributions

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.

Competing interests

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

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

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.

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 (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Department of Biology, Davidson College, Davidson, USA
(2)
Department of Biology, Missouri Western State University, Saint Joseph, USA
(3)
Department of Computer Science, Math, and Physics, Missouri Western State University, Saint Joseph, USA
(4)
Department of Math and Computer Science, Davidson College, Davidson, USA

References

  1. Sinha J, Topp S, Gallivan JP. From SELEX to cell dual selections for synthetic riboswitches. Methods Enzymol. 2011;497:207–20.View ArticlePubMedGoogle Scholar
  2. Lynch SA, Topp S, Gallivan JP. High-throughput screens to discover synthetic riboswitches. Methods Mol Biol Clifton NJ. 2009;540:321–33.View ArticleGoogle Scholar
  3. 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.View ArticlePubMedGoogle Scholar
  4. Boerneke MA, Hermann T. Ligand-responsive RNA mechanical switches. RNA Biol. 2015;12:780–6.View ArticlePubMedPubMed CentralGoogle Scholar
  5. 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.View ArticlePubMedGoogle Scholar
  6. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Lakhin AV, Tarantul VZ, Gening LV. Aptamers: problems, solutions and prospects. Acta Nat. 2013;5(4):34–43.Google Scholar
  8. Townshed B, Kennedy AB, Xiang JS, Smolke CD. High-throughput cellular RNA device engineering. Nat Methods. 2015;12(10):989–94.View ArticleGoogle Scholar
  9. Hallberg ZF, Su Y, Kitto RZ, Hammond MC. Engineering and in vivo applications of riboswitches. Annu Rev Biochem. 2017;20(86):515–39.View ArticleGoogle Scholar
  10. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  11. Galanie S, Thodey K, Trenchard IJ, Interrante MF, Smolke CD. Complete biosynthesis of opioids in yeast. Science. 2015;349(6252):7.View ArticleGoogle Scholar
  12. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  13. 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.View ArticlePubMedGoogle Scholar
  14. Mishler DM, Gallivan JP. A family of synthetic riboswitches adopts a kinetic trapping mechanism. Nucleic Acids Res. 2014;42(10):6753–61.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Felletti M, Hartig JS. Ligand-dependent ribozymes: ligand-dependent ribozymes. Wiley Interdiscip Rev RNA. 2017;8(2):e1395.View ArticleGoogle Scholar
  16. 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.View ArticlePubMedGoogle Scholar

Copyright

© The Author(s) 2018

Advertisement