Disease-associated human genetic variation through the lens of precursor and mature RNA structure

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

An intronic RNA element modulates Factor VIII exon-16 splicing

Pathogenic variants in the human Factor VIII (F8) gene cause Hemophilia A (HA). Here, we investigated the impact of 97 HA-causing single-nucleotide variants on the splicing of 11 exons from F8. For the majority of F8 exons, splicing was insensitive to the presence of HA-causing variants. However, splicing of several exons, including exon-16, was impacted by variants predicted to alter exonic splicing regulatory sequences. Using exon-16 as a model, we investigated the structure-function relationship of HA-causing variants on splicing. Intriguingly, RNA chemical probing analyses revealed a three-way junction structure at the 3'-end of intron-15 (TWJ-3-15) capable of sequestering the polypyrimidine tract. We discovered antisense oligonucleotides (ASOs) targeting TWJ-3-15 partially rescue splicing-deficient exon-16 variants by increasing accessibility of the polypyrimidine tract. The apical stem loop region of TWJ-3-15 also contains two hnRNPA1-dependent intronic splicing silencers (ISSs). ASOs blocking these ISSs also partially rescued splicing. When used in combination, ASOs targeting both the ISSs and the region sequestering the polypyrimidine tract, fully rescue pre-mRNA splicing of multiple HA-linked variants of exon-16. Together, our data reveal a putative RNA structure that sensitizes F8 exon-16 to aberrant splicing.

The Effect of Pseudoknot Base Pairing on Cotranscriptional Structural Switching of the Fluoride Riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of E. coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37 °C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65 °C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

Structural analysis of MALAT1 long noncoding RNA in cells and in evolution

Although not canonically polyadenylated, the long noncoding RNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) is stabilized by a highly conserved 76-nt triple helix structure on its 3' end. The entire MALAT1 transcript is over 8000 nt long in humans. The strongest structural conservation signal in MALAT1 (as measured by covariation of base pairs) is in the triple helix structure. Primary sequence analysis of covariation alone does not reveal the degree of structural conservation of the entire full-length transcript, however. Furthermore, RNA structure is often context dependent; RNA binding proteins that are differentially expressed in different cell types may alter structure. We investigate here the in-cell and cell-free structures of the full-length human and green monkey (Chlorocebus sabaeus) MALAT1 transcripts in multiple tissue-derived cell lines using SHAPE chemical probing. Our data reveal levels of uniform structural conservation in different cell lines, in cells and cell-free, and even between species, despite significant differences in primary sequence. The uniformity of the structural conservation across the entire transcript suggests that, despite seeing covariation signals only in the triple helix junction of the lncRNA, the rest of the transcript's structure is remarkably conserved, at least in primates and across multiple cell types and conditions.

Rescue of blood coagulation Factor VIII exon-16 mis-splicing by antisense oligonucleotides

The human Factor VIII ( F8 ) protein is essential for the blood coagulation cascade and specific F8 mutations cause the rare bleeding disorder Hemophilia A (HA). Here, we investigated the impact of HA-causing single-nucleotide mutations on F8 pre-mRNA splicing. We found that 14/97 (∼14.4%) coding sequence mutations tested in our study induced exon skipping. Splicing patterns of 4/11 (∼36.4%) F8 exons tested were especially sensitive to the presence of common disease-causing mutations. RNA-chemical probing analyses revealed a three-way junction structure at the 3' end of intron 15 (TWJ-3-15). TWJ-3-15 sequesters the polypyrimidine tract, a key determinant of 3' splice site strength. Using exon-16 of the F8 gene as a model, we designed specific antisense oligonucleotides (ASOs) that target TWJ-3-15 and identified three that promote the splicing of F8 exon-16. Interaction of TWJ-3-15 with ASOs increases accessibility of the polypyrimidine tract and inhibits the binding of hnRNPA1-dependent splicing silencing factors. Moreover, ASOs targeting TWJ-3-15 rescue diverse splicing-sensitive HA-causing mutations, most of which are distal to the 3' splice site being impacted. The TWJ-3-15 structure and its effect on mRNA splicing provide a model for HA etiology in patients harboring specific F8 mutations and provide a framework for precision RNA-based HA therapies.

How does precursor RNA structure influence RNA processing and gene expression?

RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.

Predicting functional riboSNitches in the context of alternative splicing

RNAs are the major regulators of gene expression, and their secondary structures play crucial roles at different levels. RiboSNitches are disease-associated SNPs that cause changes in the pre-mRNA secondary structural ensemble. Several riboSNitches have been detected in the 5' and 3' untranslated regions and lncRNA. Although cases of secondary structural elements playing a regulatory role in alternative splicing are known, regions specific to splicing events, such as splice junctions have not received much attention. We tested splice-site mutations for their efficiency in disrupting the secondary structure and hypothesized that these could play a crucial role in alternative splicing. Multiple riboSNitch prediction methods were applied to obtain overlapping results that are potentially more reliable. Putative riboSNitches were identified from aberrant 5' and 3' splice site mutations, cancer-causing somatic mutations, and genes that harbor the regulatory RNA secondary structural elements. Our workflow for predicting riboSNitches associated with alternative splicing is novel and paves the way for subsequent experimental validation.

A novel algorithm for ranking RNA structure candidates

RNA research is advancing at an ever increasing pace. The newest and most state-of-the-art instruments and techniques have made possible the discoveries of new RNAs, and they have carried the field to new frontiers of disease research, vaccine development, therapeutics, and architectonics. Like proteins, RNAs show a marked relationship between structure and function. A deeper grasp of RNAs requires a finer understanding of their elaborate structures. In pursuit of this, cutting-edge experimental and computational structure-probing techniques output several candidate geometries for a given RNA, each of which is perfectly aligned with experimentally determined parameters. Identifying which structure is the most accurate, however, remains a major obstacle. In recent years, several algorithms have been developed for ranking candidate RNA structures in order from most to least probable, though their levels of accuracy and transparency leave room for improvement. Most recently, advances in both areas are demonstrated by rsRNASP, a novel algorithm proposed by Tan et al. rsRNASP is a residue-separation-based statistical potential for three-dimensional structure evaluation, and it outperforms the leading algorithms in the field.