RNA editing by T7 RNA polymerase bypasses InDel mutations causing unexpected phenotypic changes

Nanopore-based direct sequencing of RNA transcripts with 10 different modified nucleotides reveals gaps in existing technology

RNA undergoes complex posttranscriptional processing including chemical modifications of the nucleotides. The resultant-modified nucleotides are an integral part of RNA sequences that must be considered in studying the biology of RNA and in the design of RNA therapeutics. However, the current "RNA-sequencing" methods primarily sequence complementary DNA rather than RNA itself, which means that the modifications present in RNA are not captured in the sequencing results. Emerging direct RNA-sequencing technologies, such as those offered by Oxford Nanopore, aim to address this limitation. In this study, we synthesized and used Nanopore technology to sequence RNA transcripts consisting of canonical nucleotides and 10 different modifications in various concentrations. The results show that direct RNA sequencing still has a baseline error rate of >10%, and although some modifications can be detected, many remain unidentified. Thus, there is a need to develop sequencing technologies and analysis methods that can comprehensively capture the total complexity of RNA. The RNA sequences obtained through this project are made available for benchmarking analysis methods.

Recent Advances in Novel Recombinant RNAs for Studying Post-transcriptional Gene Regulation in Drug Metabolism and Disposition

Drug-metabolizing enzymes and transporters are major determinants of the absorption, disposition, metabolism, and excretion (ADME) of drugs, and changes in ADME gene expression or function may alter the pharmacokinetics/ pharmacodynamics (PK/PD) and further influence drug safety and therapeutic outcomes. ADME gene functions are controlled by diverse factors, such as genetic polymorphism, transcriptional regulation, and coadministered medications. MicroRNAs (miRNAs) are a superfamily of regulatory small noncoding RNAs that are transcribed from the genome to regulate target gene expression at the post-transcriptional level. The roles of miRNAs in controlling ADME gene expression have been demonstrated, and such miRNAs may consequently influence cellular drug metabolism and disposition capacity. Several types of miRNA mimics and small interfering RNA (siRNA) reagents have been developed and widely used for ADME research. In this review article, we first provide a brief introduction to the mechanistic actions of miRNAs in post-transcriptional gene regulation of drug-metabolizing enzymes, transporters, and transcription factors. After summarizing conventional small RNA production methods, we highlight the latest advances in novel recombinant RNA technologies and applications of the resultant bioengineered RNA (BioRNA) agents to ADME studies. BioRNAs produced in living cells are not only powerful tools for general biological and biomedical research but also potential therapeutic agents amenable to clinical investigations.

Development of Host-Orthogonal Genetic Systems for Synthetic Biology

The construction of a host-orthogonal genetic system can not only minimize the impact of host-specific nuances on fine-tuning of gene expression, but also expand cellular functions such as in vivo continuous evolution of genes based on an error-prone DNA polymerase. It represents an emerging powerful approach for making biology easier to engineer. In this review, the recent advances are described on the design of genetic systems that can be stably inherited in the host cells and are responsible for important biological processes including DNA replication, RNA transcription, protein translation, and gene regulation. Their applications in synthetic biology are summarized and the future challenges and opportunities are discussed in developing such systems.

Mycobacterium tuberculosis and M. bovis BCG Moreau Fumarate Reductase Operons Produce Different Polypeptides That May Be Related to Non-canonical Functions

Tuberculosis is a world widespread disease, caused by Mycobacterium tuberculosis (M.tb). Although considered an obligate aerobe, this organism can resist life-limiting conditions such as microaerophily mainly due to its set of enzymes responsible for energy production and coenzyme restoration under these conditions. One of these enzymes is fumarate reductase, an heterotetrameric complex composed of a catalytic (FrdA), an iron-sulfur cluster (FrdB) and two transmembrane (FrdC and FrdD) subunits involved in anaerobic respiration and important for the maintenance of membrane potential. In this work, aiming to further characterize this enzyme function in mycobacteria, we analyzed the expression of FrdB-containing proteins in M.tb and Mycobacterium bovis Bacillus Calmette–Guérin (BCG) Moreau, the Brazilian vaccine strain against tuberculosis. We identified three isoforms in both mycobacteria, two of them corresponding to the predicted encoded polypeptides of M.tb (27 kDa) and BCG Moreau (40 kDa) frd sequences, as due to an insertion on the latter’s operon a fused FrdBC protein is expected. The third 52 kDa band can be explained by a transcriptional slippage event, typically occurring when mutation arises in a repetitive region within a coding sequence, thought to reduce its impact allowing the production of both native and variant forms. Comparative modeling of the M.tb and BCG Moreau predicted protein complexes allowed the detection of subtle overall differences, showing a high degree of structure and maybe functional resemblance among them. Axenic growth and macrophage infection assays show that the frd locus is important for proper bacterial development in both scenarios, and that both M.tb’s and BCG Moreau’s alleles can partially revert the hampered phenotype of the knockout strain. Altogether, our results show that the frdABCD operon of Mycobacteria may have evolved to possess other yet non-described functions, such as those necessary during aerobic logarithmic growth and early stage steps of infection.

Transcriptional frameshifts contribute to protein allergenicity

Transcription infidelity (TI) is a mechanism that increases RNA and protein diversity. We found that single-base omissions (i.e., gaps) occurred at significantly higher rates in the RNA of highly allergenic legumes. Transcripts from peanut, soybean, sesame, and mite allergens contained a higher density of gaps than those of nonallergens. Allergen transcripts translate into proteins with a cationic carboxy terminus depleted in hydrophobic residues. In mice, recombinant TI variants of the peanut allergen Ara h 2, but not the canonical allergen itself, induced, without adjuvant, the production of anaphylactogenic specific IgE (sIgE), binding to linear epitopes on both canonical and TI segments of the TI variants. The removal of cationic proteins from bovine lactoserum markedly reduced its capacity to induce sIgE. In peanut-allergic children, the sIgE reactivity was directed toward both canonical and TI segments of Ara h 2 variants. We discovered 2 peanut allergens, which we believe to be previously unreported, because of their RNA-DNA divergence gap patterns and TI peptide amino acid composition. Finally, we showed that the sIgE of children with IgE-negative milk allergy targeted cationic proteins in lactoserum. We propose that it is not the canonical allergens, but their TI variants, that initiate sIgE isotype switching, while both canonical and TI variants elicit clinical allergic reactions.

Novel approaches for efficient in vivo fermentation production of noncoding RNAs

Genome-derived noncoding RNAs (ncRNAs), including microRNAs (miRNAs), small interfering RNAs (siRNAs), and long noncoding RNAs (lncRNAs), play an essential role in the control of target gene expression underlying various cellular processes, and dysregulation of ncRNAs is involved in the pathogenesis and progression of various diseases in virtually all species including humans. Understanding ncRNA biology has opened new avenues to develop novel RNA-based therapeutics. Presently, ncRNA research and drug development is dominated by the use of ncRNA mimics that are synthesized chemically in vitro and supplemented with extensive and various types of artificial modifications and thus may not necessarily recapitulate the properties of natural RNAs generated and folded in living cells in vivo. Therefore, there are growing interests in developing novel technologies for in vivo production of RNA molecules. The two most recent major breakthroughs in achieving an efficient, large-scale, and cost-effective fermentation production of recombinant or bioengineered RNAs (e.g., tens of milligrams from 1 L of bacterial culture) are (1) using stable RNA carriers and (2) direct overexpression in RNase III-deficient bacteria, while other approaches offer a low yield (e.g., nano- to microgram scales per liter). In this article, we highlight these novel microbial fermentation-based technologies that have shifted the paradigm to the production of true biological ncRNA molecules for research and development.

Non-programmed transcriptional frameshifting is common and highly RNA polymerase type-dependent

Background

The viral or host systems for a gene expression assume repeatability of the process and high quality of the protein product. Since level and fidelity of transcription primarily determines the overall efficiency, all factors contributing to their decrease should be identified and optimized. Among many observed processes, non-programmed insertion/deletion (indel) of nucleotide during transcription (slippage) occurring at homopolymeric A/T sequences within a gene can considerably impact its expression. To date, no comparative study of the most utilized Escherichia coli and T7 bacteriophage RNA polymerases (RNAP) propensity for this type of erroneous mRNA synthesis has been reported. To address this issue we evaluated the influence of shift-prone A/T sequences by assessing indel-dependent phenotypic changes. RNAP-specific expression profile was examined using two of the most potent promoters, P_araBAD of E. coli and φ10 of phage T7.

Results

Here we report on the first systematic study on requirements for efficient transcriptional slippage by T7 phage and cellular RNAPs considering three parameters: homopolymer length, template type, and frameshift directionality preferences. Using a series of out-of-frame gfp reporter genes fused to a variety of A/T homopolymeric sequences we show that T7 RNAP has an exceptional potential for generating frameshifts and is capable of slipping on as few as three adenine or four thymidine residues in a row, in a flanking sequence-dependent manner. In contrast, bacterial RNAP exhibits a relatively low ability to baypass indel mutations and requires a run of at least 7 tymidine and even more adenine residues. This difference comes from involvement of various intrinsic proofreading properties. Our studies demonstrate distinct preference towards a specific homopolymer in slippage induction. Whereas insertion slippage performed by T7 RNAP (but not deletion) occurs tendentiously on poly(A) rather than on poly(T) runs, strong bias towards poly(T) for the host RNAP is observed.

Conclusions

Intrinsic RNAP slippage properties involve trade-offs between accuracy, speed and processivity of transcription. Viral T7 RNAP manifests far greater inclinations to the transcriptional slippage than E. coli RNAP. This possibly plays an important role in driving bacteriophage adaptation and therefore could be considered as beneficial. However, from biotechnological and experimental viewpoint, this might create some problems, and strongly argues for employing bacterial expression systems, stocked with proofreading mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1034-4) contains supplementary material, which is available to authorized users.

Evaluation of GFP reporter utility for analysis of transcriptional slippage during gene expression

Background

Epimutations arising from transcriptional slippage seem to have more important role in regulating gene expression than earlier though. Since the level and the fidelity of transcription primarily determine the overall efficiency of gene expression, all factors contributing to their decrease should be identified and optimized.

Results

To examine the influence of A/T homopolymeric sequences on introduction of erroneous nucleotides by slippage mechanism green fluorescence protein (GFP) reporter was chosen. The in- or out-of-frame gfp gene was fused to upstream fragment with variable number of adenine or thymine stretches resulting in several hybrid GFP proteins with diverse amino acids at N-terminus. Here, by using T7 phage expression system we showed that the intensity of GFP fluorescence mainly depends on the number of the retained natural amino acids. While the lack of serine (S₂) residue results in negligible effects, the lack of serine and lysine (S₂K₃) contributed to a significant reduction in fluorescence by 2.7-fold for polyA-based in-frame controls and twofold for polyTs. What is more, N-terminal tails amino acid composition was rather of secondary importance, since the whole-cell fluorescence differed in a range of 9–18% between corresponding polyA- and polyT-based constructs.

Conclusions

Here we present experimental evidence for utility of GFP reporter for accurate estimation of A/T homopolymeric sequence contribution in transcriptional slippage induction. We showed that the intensity of GFP hybrid fluorescence mainly depends on the number of retained natural amino acids, thus fluorescence raw data need to be referred to appropriate positive control. Moreover, only in case of GFP hybrids with relatively short N-terminal tags the fluorescence level solely reflects production yield, what further indicates the impact of an individual slippage sequence. Our results demonstrate that in contrast to the E. coli enzyme, T7 RNA polymerase exhibits extremely high propensity to slippage even on runs as short as 3 adenine or 4 thymine residues.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-0999-3) contains supplementary material, which is available to authorized users.

RNA structure inference through chemical mapping after accidental or intentional mutations

Despite the critical roles RNA structures play in regulating gene expression, sequencing-based methods for experimentally determining RNA base pairs have remained inaccurate. Here, we describe a multidimensional chemical-mapping method called "mutate-and-map read out through next-generation sequencing" (M2-seq) that takes advantage of sparsely mutated nucleotides to induce structural perturbations at partner nucleotides and then detects these events through dimethyl sulfate (DMS) probing and mutational profiling. In special cases, fortuitous errors introduced during DNA template preparation and RNA transcription are sufficient to give M2-seq helix signatures; these signals were previously overlooked or mistaken for correlated double-DMS events. When mutations are enhanced through error-prone PCR, in vitro M2-seq experimentally resolves 33 of 68 helices in diverse structured RNAs including ribozyme domains, riboswitch aptamers, and viral RNA domains with a single false positive. These inferences do not require energy minimization algorithms and can be made by either direct visual inspection or by a neural-network-inspired algorithm called M2-net. Measurements on the P4-P6 domain of the Tetrahymena group I ribozyme embedded in Xenopus egg extract demonstrate the ability of M2-seq to detect RNA helices in a complex biological environment.

Interaction Analysis of T7 RNA Polymerase with Heparin and Its Low Molecular Weight Derivatives - An In Silico Approach

The single subunit T7 RNA polymerase (T7RNAP) is a model enzyme for studying the transcription process and for various biochemical and biophysical studies. Heparin is a commonly used inhibitor against T7RNAP and other RNA polymerases. However, exact interaction between heparin and T7RNAP is still not completely understood. In this work, we analyzed the binding pattern of heparin by docking heparin and few of its low molecular weight derivatives to T7RNAP, which helps in better understanding of T7RNAP inhibition mechanism. The efficiency of the compounds was calculated by docking the selected compounds and post-docking molecular mechanics/generalized Born surface area analysis. Evaluation of the simulation trajectories and binding free energies of the complexes after simulation showed enoxaparin to be the best among low molecular weight heparins. Binding free energy analysis revealed that van der Waals interactions and polar solvation energy provided the substantial driving force for the binding process. Furthermore, per-residue free energy decomposition analysis revealed that the residues Asp 471, Asp 506, Asp 537, Tyr 571, Met 635, Asp 653, Pro 780, and Asp 812 are important for heparin interaction. Apart from these residues, most favorable contribution in all the three complexes came from Asp 506, Tyr 571, Met 635, Glu 652, and Asp 653, which can be essential for binding of heparin-like structures with T7RNAP. The results obtained from this study will be valuable for the future rational design of novel and potent inhibitors against T7RNAP and related proteins.

The P1N-PISPO trans-Frame Gene of Sweet Potato Feathery Mottle Potyvirus Is Produced during Virus Infection and Functions as an RNA Silencing Suppressor

The positive-sense RNA genome of Sweet potato feathery mottle virus (SPFMV) (genus Potyvirus, family Potyviridae) contains a large open reading frame (ORF) of 3,494 codons translatable as a polyprotein and two embedded shorter ORFs in the -1 frame: PISPO, of 230 codons, and PIPO, of 66 codons, located in the P1 and P3 regions, respectively. PISPO is specific to some sweet potato-infecting potyviruses, while PIPO is present in all potyvirids. In SPFMV these two extra ORFs are preceded by conserved G2A6 motifs. We have shown recently that a polymerase slippage mechanism at these sites could produce transcripts bringing these ORFs in frame with the upstream polyprotein, thus leading to P1N-PISPO and P3N-PIPO products (B. Rodamilans, A. Valli, A. Mingot, D. San Leon, D. B. Baulcombe, J. J. Lopez-Moya, and J.A. Garcia, J Virol 89:6965-6967, 2015, doi:10.1128/JVI.00337-15). Here, we demonstrate by liquid chromatography coupled to mass spectrometry that both P1 and P1N-PISPO are produced during viral infection and coexist in SPFMV-infected Ipomoea batatas plants. Interestingly, transient expression of SPFMV gene products coagroinfiltrated with a reporter gene in Nicotiana benthamiana revealed that P1N-PISPO acts as an RNA silencing suppressor, a role normally associated with HCPro in other potyviruses. Moreover, mutation of WG/GW motifs present in P1N-PISPO abolished its silencing suppression activity, suggesting that the function might require interaction with Argonaute components of the silencing machinery, as was shown for other viral suppressors. Altogether, our results reveal a further layer of complexity of the RNA silencing suppression activity within the Potyviridae family.

IMPORTANCE

Gene products of potyviruses include P1, HCPro, P3, 6K1, CI, 6K2, VPg/NIaPro, NIb, and CP, all derived from the proteolytic processing of a large polyprotein, and an additional P3N-PIPO product, with the PIPO segment encoded in a different frame within the P3 cistron. In sweet potato feathery mottle virus (SPFMV), another out-of-frame element (PISPO) was predicted within the P1 region. We have shown recently that a polymerase slippage mechanism can generate the transcript variants with extra nucleotides that could be translated into P1N-PISPO and P3N-PIPO. Now, we demonstrate by mass spectrometry analysis that P1N-PISPO is indeed produced in SPFMV-infected plants, in addition to P1. Interestingly, while in other potyviruses the suppressor of RNA silencing is HCPro, we show here that P1N-PISPO exhibited this activity in SPFMV, revealing how the complexity of the gene content could contribute to supply this essential function in members of the Potyviridae family.

Bioengineering of noncoding RNAs for research agents and therapeutics

The discovery of functional small noncoding RNAs (ncRNAs), such as microRNAs and small interfering RNAs, in the control of human cellular processes has opened new avenues to develop RNA-based therapies for various diseases including viral infections and cancers. However, studying ncRNA functions and developing RNA-based therapeutics relies on access to large quantities of affordable ncRNA agents. Currently, synthetic RNAs account for the major source of agents for RNA research and development, yet carry artificial modifications on the ribose ring and phosphate backbone in sharp contrast to posttranscriptional modifications present on the nucleobases or unmodified natural RNA molecules produced within cells. Therefore, large efforts have been made in recent years to develop recombinant RNA techniques to cost-effectively produce biological RNA agents that may better capture the structure, function, and safety properties of natural RNAs. In this article, we summarize and compare current in vitro and in vivo methods for the production of RNA agents including chemical synthesis, in vitro transcription, and bioengineering approaches. We highlight the latest recombinant RNA approaches using transfer RNA (tRNA), ribosomal RNA (rRNA), and optimal ncRNA scaffold (OnRS), and discuss the applications of bioengineered ncRNA agents (BERAs) that should facilitate RNA research and development.