Systematic Quantification of Sequence and Structural Determinants Controlling mRNA stability in Bacterial Operons

ABSTRACT mRNA degradation is a central process that affects all gene expression levels, and yet the determinants that control mRNA decay rates remain poorly characterized. Here, we applied a synthetic biology, learn-by-design approach to elucidate the sequence and structural determinants that control mRNA stability in bacterial operons. We designed, constructed, and characterized 82 operons, systematically varying RNAse binding site characteristics, translation initiation rates, and transcriptional terminator efficiencies in the 5’ UTR, intergenic, and 3’ UTR regions, and measuring their mRNA levels using RT-qPCR assays. We show that introducing long single-stranded RNA into 5’ UTRs reduced mRNA levels by up to 9.4-fold and that lowering translation rates reduced mRNA levels by up to 11.8-fold. We also found that RNAse binding sites in intergenic regions had much lower effects on mRNA levels. Surprisingly, changing transcriptional termination efficiency or introducing long single-stranded RNA into 3’ UTRs had no effect on upstream mRNA levels. From these measurements, we developed and validated biophysical models of ribosome protection and RNAse activity with excellent quantitative correspondence. We also formulated design rules to rationally control a mRNA’s stability, facilitating the automated design of engineered genetic systems with desired functionalities.