RNA structure and splicing regulation

U1 snRNA interactions with deep intronic sequences regulate splicing of multiple exons of spinal muscular atrophy genes

Introduction

The U1 small nuclear RNA (snRNA) forms ribonucleoprotein particles (RNPs) such as U1 snRNP and U1-TAF15 snRNP. U1 snRNP is one of the most studied RNPs due to its critical role in pre-mRNA splicing in defining the 5' splice site (5'ss) of every exon through direct interactions with sequences at exon/intron junctions. Recent reports support the role of U1 snRNP in all steps of transcription, namely initiation, elongation, and termination. Functions of U1-TAF15 snRNP are less understood, though it associates with the transcription machinery and may modulate pre-mRNA splicing by interacting with the 5'ss and/or 5'ss-like sequences within the pre-mRNA. An anti-U1 antisense oligonucleotide (ASO) that sequesters the 5' end of U1 snRNA inhibits the functions of U1 snRNP, including transcription and splicing. However, it is not known if the inhibition of U1 snRNP influences post-transcriptional regulation of pre-mRNA splicing through deep intronic sequences.

Methods

We examined the effect of an anti-U1 ASO that sequesters the 5' end of U1 snRNA on transcription and splicing of all internal exons of the spinal muscular atrophy (SMA) genes, SMN1 and SMN2. Our study was enabled by the employment of a multi-exon-skipping detection assay (MESDA) that discriminates against prematurely terminated transcripts. We employed an SMN2 super minigene to determine if anti-U1 ASO differently affects splicing in the context of truncated introns.

Results

We observed substantial skipping of multiple internal exons of SMN1 and SMN2 triggered by anti-U1 treatment. Suggesting a role for U1 snRNP in interacting with deep intronic sequences, early exons of the SMN2 super minigene with truncated introns were resistant to anti-U1 induced skipping. Consistently, overexpression of engineered U1 snRNAs targeting the 5'ss of early SMN1 and SMN2 exons did not prevent exon skipping caused by anti-U1 treatment.

Discussion

Our results uncover a unique role of the U1 snRNA-associated RNPs in splicing regulation executed through deep intronic sequences. Findings are significant for developing novel therapies for SMA based on deep intronic targets.

Long-range RNA structures in the human transcriptome beyond evolutionarily conserved regions

RNA structure has been increasingly recognized as a critical player in the biogenesis and turnover of many transcripts classes. In eukaryotes, the prediction of RNA structure by thermodynamic modeling meets fundamental limitations due to the large sizes and complex, discontinuous organization of eukaryotic genes. Signatures of functional RNA structures can be found by detecting compensatory substitutions in homologous sequences, but a comparative approach is applicable only within conserved sequence blocks. Here, we developed a computational pipeline called PHRIC, which is not limited to conserved regions and relies on RNA contacts derived from RNA in situ conformation sequencing (RIC-seq) experiments. It extracts pairs of short RNA fragments surrounded by nested clusters of RNA contacts and predicts long, nearly perfect complementary base pairings formed between these fragments. In application to a panel of RIC-seq experiments in seven human cell lines, PHRIC predicted ~12,000 stable long-range RNA structures with equilibrium free energy below -15 kcal/mol, the vast majority of which fall outside of regions annotated as conserved among vertebrates. These structures, nevertheless, show some level of sequence conservation and remarkable compensatory substitution patterns in other clades. Furthermore, we found that introns have a higher propensity to form stable long-range RNA structures between each other, and moreover that RNA structures tend to concentrate within the same intron rather than connect adjacent introns. These results for the first time extend the application of proximity ligation assays to RNA structure prediction beyond conserved regions.

Cooperation and Competition of RNA Secondary Structure and RNA-Protein Interactions in the Regulation of Alternative Splicing

The regulation of alternative splicing in eukaryotic cells is carried out through the coordinated action of a large number of factors, including RNA-binding proteins and RNA structure. The RNA structure influences alternative splicing by blocking cis-regulatory elements, or bringing them closer or farther apart. In combination with RNA-binding proteins, it generates transcript conformations that help to achieve the necessary splicing outcome. However, the binding of regulatory proteins depends on RNA structure and, vice versa, the formation of RNA structure depends on the interaction with regulators. Therefore, RNA structure and RNA-binding proteins are inseparable components of common regulatory mechanisms. This review highlights examples of alternative splicing regulation by RNA-binding proteins, the regulation through local and long-range RNA structures, as well as how these elements work together, cooperate, and compete.

RNA in situ conformation sequencing reveals novel long-range RNA structures with impact on splicing

Over recent years, long-range RNA structure has emerged as a factor that is fundamental to alternative splicing regulation. An increasing number of human disorders are now being associated with splicing defects; hence it is essential to develop methods that assess long-range RNA structure experimentally. RNA in situ conformation sequencing (RIC-seq) is a method that recapitulates RNA structure within physiological RNA-protein complexes. In this work, we juxtapose pairs of conserved complementary regions (PCCRs) that were predicted in silico with the results of RIC-seq experiments conducted in seven human cell lines. We show statistically that RIC-seq support of PCCRs correlates with their properties, such as equilibrium free energy, presence of compensatory substitutions, and occurrence of A-to-I RNA editing sites and forked eCLIP peaks. Exons enclosed in PCCRs that are supported by RIC-seq tend to have weaker splice sites and lower inclusion rates, which is indicative of post-transcriptional splicing regulation mediated by RNA structure. Based on these findings, we prioritize PCCRs according to their RIC-seq support and show, using antisense nucleotides and minigene mutagenesis, that PCCRs in two disease-associated human genes, PHF20L1 and CASK, and also PCCRs in their murine orthologs, impact alternative splicing. In sum, we demonstrate how RIC-seq experiments can be used to discover functional long-range RNA structures, and particularly those that regulate alternative splicing.

Diverse targets of SMN2-directed splicing-modulating small molecule therapeutics for spinal muscular atrophy

Designing an RNA-interacting molecule that displays high therapeutic efficacy while retaining specificity within a broad concentration range remains a challenging task. Risdiplam is an FDA-approved small molecule for the treatment of spinal muscular atrophy (SMA), the leading genetic cause of infant mortality. Branaplam is another small molecule which has undergone clinical trials. The therapeutic merit of both compounds is based on their ability to restore body-wide inclusion of Survival Motor Neuron 2 (SMN2) exon 7 upon oral administration. Here we compare the transcriptome-wide off-target effects of these compounds in SMA patient cells. We captured concentration-dependent compound-specific changes, including aberrant expression of genes associated with DNA replication, cell cycle, RNA metabolism, cell signaling and metabolic pathways. Both compounds triggered massive perturbations of splicing events, inducing off-target exon inclusion, exon skipping, intron retention, intron removal and alternative splice site usage. Our results of minigenes expressed in HeLa cells provide mechanistic insights into how these molecules targeted towards a single gene produce different off-target effects. We show the advantages of combined treatments with low doses of risdiplam and branaplam. Our findings are instructive for devising better dosing regimens as well as for developing the next generation of small molecule therapeutics aimed at splicing modulation.

Analysis of the impact of DGAT1 p.M435L and p.K232A variants on pre-mRNA splicing in a full-length gene assay

DGAT1 is playing a major role in fat metabolism and triacylglyceride synthesis. Only two DGAT1 loss-of-function variants altering milk production traits in cattle have been reported to date, namely p.M435L and p.K232A. The p.M435L variant is a rare alteration and has been associated with skipping of exon 16 which results in a non-functional truncated protein, and the p.K232A-containing haplotype has been associated with modifications of the splicing rate of several DGAT1 introns. In particular, the direct causality of the p.K232A variant in decreasing the splicing rate of the intron 7 junction was validated using a minigene assay in MAC-T cells. As both these DGAT1 variants were shown to be spliceogenic, we developed a full-length gene assay (FLGA) to re-analyse p.M435L and p.K232A variants in HEK293T and MAC-T cells. Qualitative RT-PCR analysis of cells transfected with the full-length DGAT1 expression construct carrying the p.M435L variant highlighted complete skipping of exon 16. The same analysis performed using the construct carrying the p.K232A variant showed moderate differences compared to the wild-type construct, suggesting a possible effect of this variant on the splicing of intron 7. Finally, quantitative RT-PCR analyses of cells transfected with the p.K232A-carrying construct did not show any significant modification on the splicing rate of introns 1, 2 and 7. In conclusion, the DGAT1 FLGA confirmed the p.M435L impact previously observed in vivo, but invalidated the hypothesis whereby the p.K232A variant strongly decreased the splicing rate of intron 7.

ShapeSorter: a fully probabilistic method for detecting conserved RNA structure features supported by SHAPE evidence

There is an increased interest in the determination of RNA structures in vivo as it is now possible to probe them in a high-throughput manner, e.g. using SHAPE protocols. By now, there exist a range of computational methods that integrate experimental SHAPE-probing evidence into computational RNA secondary structure prediction. The state-of-the-art in this field is currently provided by computational methods that employ the minimum-free energy strategy for prediction RNA secondary structures with SHAPE-probing evidence. These methods, however, rely on the assumption that transcripts in vivo fold into the thermodynamically most stable configuration and ignore evolutionary evidence for conserved RNA structure features. We here present a new computational method, ShapeSorter, that predicts RNA structure features without employing the thermodynamic strategy. Instead, ShapeSorter employs a fully probabilistic framework to identify RNA structure features that are supported by evolutionary and SHAPE-probing evidence. Our method can capture RNA structure heterogeneity, pseudo-knotted RNA structures as well as transient and mutually exclusive RNA structure features. Moreover, it estimates P-values for the predicted RNA structure features which allows for easy filtering and ranking. We investigate the merits of our method in a comprehensive performance benchmarking and conclude that ShapeSorter has a significantly superior performance for predicting base-pairs than the existing state-of-the-art methods.

High Concentration of an ISS-N1-Targeting Antisense Oligonucleotide Causes Massive Perturbation of the Transcriptome

Intronic splicing silencer N1 (ISS-N1) located within Survival Motor Neuron 2 (SMN2) intron 7 is the target of a therapeutic antisense oligonucleotide (ASO), nusinersen (Spinraza), which is currently being used for the treatment of spinal muscular atrophy (SMA), a leading genetic disease associated with infant mortality. The discovery of ISS-N1 as a promising therapeutic target was enabled in part by Anti-N1, a 20-mer ASO that restored SMN2 exon 7 inclusion by annealing to ISS-N1. Here, we analyzed the transcriptome of SMA patient cells treated with 100 nM of Anti-N1 for 30 h. Such concentrations are routinely used to demonstrate the efficacy of an ASO. While 100 nM of Anti-N1 substantially stimulated SMN2 exon 7 inclusion, it also caused massive perturbations in the transcriptome and triggered widespread aberrant splicing, affecting expression of essential genes associated with multiple cellular processes such as transcription, splicing, translation, cell signaling, cell cycle, macromolecular trafficking, cytoskeletal dynamics, and innate immunity. We validated our findings with quantitative and semiquantitative PCR of 39 candidate genes associated with diverse pathways. We also showed a substantial reduction in off-target effects with shorter ISS-N1-targeting ASOs. Our findings are significant for implementing better ASO design and dosing regimens of ASO-based drugs.

Conserved long-range base pairings are associated with pre-mRNA processing of human genes

The ability of nucleic acids to form double-stranded structures is essential for all living systems on Earth. Current knowledge on functional RNA structures is focused on locally-occurring base pairs. However, crosslinking and proximity ligation experiments demonstrated that long-range RNA structures are highly abundant. Here, we present the most complete to-date catalog of conserved complementary regions (PCCRs) in human protein-coding genes. PCCRs tend to occur within introns, suppress intervening exons, and obstruct cryptic and inactive splice sites. Double-stranded structure of PCCRs is supported by decreased icSHAPE nucleotide accessibility, high abundance of RNA editing sites, and frequent occurrence of forked eCLIP peaks. Introns with PCCRs show a distinct splicing pattern in response to RNAPII slowdown suggesting that splicing is widely affected by co-transcriptional RNA folding. The enrichment of 3'-ends within PCCRs raises the intriguing hypothesis that coupling between RNA folding and splicing could mediate co-transcriptional suppression of premature pre-mRNA cleavage and polyadenylation.

The First Orally Deliverable Small Molecule for the Treatment of Spinal Muscular Atrophy

Spinal muscular atrophy (SMA) is 1 of the leading causes of infant mortality. SMA is mostly caused by low levels of Survival Motor Neuron (SMN) protein due to deletion of or mutation in the SMN1 gene. Its nearly identical copy, SMN2, fails to compensate for the loss of SMN1 due to predominant skipping of exon 7. Correction of SMN2 exon 7 splicing by an antisense oligonucleotide (ASO), nusinersen (Spinraza™), that targets the intronic splicing silencer N1 (ISS-N1) became the first approved therapy for SMA. Restoration of SMN levels using gene therapy was the next. Very recently, an orally deliverable small molecule, risdiplam (Evrysdi™), became the third approved therapy for SMA. Here we discuss how these therapies are positioned to meet the needs of the broad phenotypic spectrum of SMA patients.

A survey of transcripts generated by spinal muscular atrophy genes

Human Survival Motor Neuron (SMN) genes code for SMN, an essential multifunctional protein. Complete loss of SMN is embryonic lethal, while low levels of SMN lead to spinal muscular atrophy (SMA), a major genetic disease of children and infants. Reduced levels of SMN are associated with the abnormal development of heart, lung, muscle, gastro-intestinal system and testis. The SMN loci have been shown to generate a vast repertoire of transcripts, including linear, back- and trans-spliced RNAs as well as antisense long noncoding RNAs. However, functions of the majority of these transcripts remain unknown. Here we review the nature of RNAs generated from the SMN loci and discuss their potential functions in cellular metabolism.

RNA in spinal muscular atrophy: therapeutic implications of targeting

INTRODUCTION

Spinal muscular atrophy (SMA) is caused by low levels of the Survival Motor Neuron (SMN) protein due to deletions of or mutations in the SMN1 gene. Humans carry another nearly identical gene, SMN2, which mostly produces a truncated and less stable protein SMNΔ7 due to predominant skipping of exon 7. Elevation of SMN upon correction of SMN2 exon 7 splicing and gene therapy have been proven to be the effective treatment strategies for SMA.

AREAS COVERED

This review summarizes existing and potential SMA therapies that are based on RNA targeting.We also discuss the mechanistic basis of RNA-targeting molecules.

EXPERT OPINION

The discovery of intronic splicing silencer N1 (ISS-N1) was the first major step towards developing the currently approved antisense-oligonucleotide (ASO)-directed therapy (SpinrazaTM) based on the correction of exon 7 splicing of the endogenous SMN2pre-mRNA. Recently, gene therapy (Zolgensma) has become the second approved treatment for SMA. Small compounds (currently in clinical trials) capable of restoring SMN2 exon 7 inclusion further expand the class of the RNA targeting molecules for SMA therapy. Endogenous RNA targets, such as long non-coding RNAs, circular RNAs, microRNAs and ribonucleoproteins, could be potentially exploited for developing additional SMA therapies.