Zika virus infection in a cell culture model reflects the transcriptomic signatures in patients.

Zika virus (ZIKV), a re-emerging flavivirus, is associated with devasting developmental and neurological disease outcomes particularly in infants infected in utero. Towards understanding the molecular underpinnings of the unique ZIKV disease pathologies, numerous transcriptome-wide studies have been undertaken. Notably, these studies have overlooked the assimilation of RNA-seq analysis from ZIKV-infected patients with cell culture model systems. In this study we find that ZIKV-infection of human lung adenocarcinoma A549 cells, mirrored both the transcriptional and alternative splicing profiles from previously published RNA-seq data of peripheral blood mononuclear cells collected from pediatric patients during early acute, late acute, and convalescent phases of ZIKV infection. Our analyses show that ZIKV infection in cultured cells correlates with transcriptional changes in patients, while the overlap in alternative splicing profiles was not as extensive. Overall, our data indicate that cell culture model systems support dissection of select molecular changes detected in patients and establishes the groundwork for future studies elucidating the biological implications of alternative splicing during ZIKV infection.

Mutation of two intronic nucleotides alters RNA structure and dynamics inhibiting MBNL1 and RBFOX1 regulated splicing of the Insulin Receptor.

Alternative splicing (AS) of Exon 11 of the Insulin Receptor ( INSR ) is highly regulated and disrupted in several human disorders. To better understand INSR exon 11 AS regulation, splicing activity of an INSR exon 11 minigene reporter was measured across a gradient of the AS regulator muscleblind-like 1 protein (MBNL1). The RNA-binding protein Fox-1 (RBFOX1) was added to determine its impact on MBNL1-regulated splicing. The role of the RBFOX1 UGCAUG binding site within intron 11 was assessed across the MBNL1 gradient. Mutating the UGCAUG motif inhibited RBFOX1 regulation of exon 11 and had the unexpected effect of reducing MBNL1 regulation of this exon. Molecular dynamics simulations showed that exon 11 and the adjacent RNA adopts a dynamically stable conformation. Mutation of the RBFOX1 binding site altered RNA structure and dynamics, while a mutation that created an optimal MBNL1 binding site at the RBFOX1 site shifted the RNA back to wild type. An antisense oligonucleotide (ASO) was used to confirm the structure in this region of the pre-mRNA. This example of intronic mutations shifting pre-mRNA structure and dynamics to modulate splicing suggests RNA structure and dynamics should be taken into consideration for AS regulation and therapeutic interventions targeting pre-mRNA.

Widespread alternative splicing dysregulation occurs presymptomatically in CAG expansion spinocerebellar ataxias.

The spinocerebellar ataxias (SCAs) are a group of dominantly inherited neurodegenerative diseases, several of which are caused by CAG expansion mutations (SCAs 1, 2, 3, 6, 7 and 12) and more broadly belong to the large family of over 40 microsatellite expansion diseases. While dysregulation of alternative splicing is a well defined driver of disease pathogenesis across several microsatellite diseases, the contribution of alternative splicing in CAG expansion SCAs is poorly understood. Furthermore, despite extensive studies on differential gene expression, there remains a gap in our understanding of presymptomatic transcriptomic drivers of disease. We sought to address these knowledge gaps through a comprehensive study of 29 publicly available RNA-sequencing datasets. We identified that dysregulation of alternative splicing is widespread across CAG expansion mouse models of SCAs 1, 3 and 7. These changes were detected presymptomatically, persisted throughout disease progression, were repeat length-dependent, and were present in brain regions implicated in SCA pathogenesis including the cerebellum, pons and medulla. Across disease progression, changes in alternative splicing occurred in genes that function in pathways and processes known to be impaired in SCAs, such as ion channels, synaptic signalling, transcriptional regulation and the cytoskeleton. We validated several key alternative splicing events with known functional consequences, including Trpc3 exon 9 and Kcnma1 exon 23b, in the Atxn1154Q/2Q mouse model. Finally, we demonstrated that alternative splicing dysregulation is responsive to therapeutic intervention in CAG expansion SCAs with Atxn1 targeting antisense oligonucleotide rescuing key splicing events. Taken together, these data demonstrate that widespread presymptomatic dysregulation of alternative splicing in CAG expansion SCAs may contribute to disease onset, early neuronal dysfunction and may represent novel biomarkers across this devastating group of neurodegenerative disorders.

Individual transcriptomic response to strength training for patients with myotonic dystrophy type 1.

Myotonic dystrophy type 1 (DM1), the most common form of adult-onset muscular dystrophy, is caused by a CTG expansion resulting in significant transcriptomic dysregulation that leads to muscle weakness and wasting. While strength training is clinically beneficial in DM1, molecular effects had not been studied. To determine whether training rescued transcriptomic defects, RNA-Seq was performed on vastus lateralis samples from 9 male patients with DM1 before and after a 12-week strength-training program and 6 male controls who did not undergo training. Differential gene expression and alternative splicing analysis were correlated with the one-repetition maximum strength evaluation method (leg extension, leg press, hip abduction, and squat). While training program-induced improvements in splicing were similar among most individuals, rescued splicing events varied considerably between individuals. Gene expression improvements were highly varied between individuals, and the percentage of differentially expressed genes rescued after training were strongly correlated with strength improvements. Evaluating transcriptome changes individually revealed responses to the training not evident from grouped analysis, likely due to disease heterogeneity and individual exercise response differences. Our analyses indicate that transcriptomic changes are associated with clinical outcomes in patients with DM1 undergoing training and that these changes are often specific to the individual and should be analyzed accordingly.

Quercetin selectively reduces expanded repeat RNA levels in models of myotonic dystrophy.

Myotonic dystrophy is a multisystemic neuromuscular disease caused by either a CTG repeat expansion in DMPK (DM1) or a CCTG repeat expansion in CNBP (DM2). Transcription of the expanded alleles produces toxic gain-of-function RNA that sequester the MBNL family of alternative splicing regulators into ribonuclear foci, leading to pathogenic mis-splicing. There are currently no approved treatments that target the root cause of disease which is the production of the toxic expansion RNA molecules. In this study, using our previously established HeLa DM1 repeat selective screening platform, we identified the natural product quercetin as a selective modulator of toxic RNA levels. Quercetin treatment selectively reduced toxic RNA levels and rescued MBNL dependent mis-splicing in DM1 and DM2 patient derived cell lines and in the HSALR transgenic DM1 mouse model where rescue of myotonia was also observed. Based on our data and its safety profile for use in humans, we have identified quercetin as a priority disease-targeting therapeutic lead for clinical evaluation for the treatment of DM1 and DM2.

Disease-associated inosine misincorporation into RNA hinders translation.

Failure to prevent accumulation of the non-canonical nucleotide inosine triphosphate (ITP) by inosine triphosphate pyrophosphatase (ITPase) during nucleotide synthesis results in misincorporation of inosine into RNA and can cause severe and fatal developmental anomalies in humans. While the biochemical activity of ITPase is well understood, the pathogenic basis of ITPase deficiency and the molecular and cellular consequences of ITP misincorporation into RNA remain cryptic. Here, we demonstrate that excess ITP in the nucleotide pool during in vitro transcription results in T7 polymerase-mediated inosine misincorporation in luciferase RNA. In vitro translation of inosine-containing luciferase RNA reduces resulting luciferase activity, which is only partly explained by reduced abundance of the luciferase protein produced. Using Oxford Nanopore Direct RNA sequencing, we reveal inosine misincorporation to be stochastic but biased largely towards misincorporation in place of guanosine, with evidence for misincorporation also in place of cytidine, adenosine and uridine. Inosine misincorporation into RNA is also detected in Itpa-null mouse embryonic heart tissue as an increase in relative variants compared with the wild type using Illumina RNA sequencing. By generating CRISPR/Cas9 rat H9c2 Itpa-null cardiomyoblast cells, we validate a translation defect in cells that accumulate inosine within endogenous RNA. Furthermore, we observe hindered cellular translation of transfected luciferase RNA containing misincorporated inosine in both wild-type and Itpa-null cells. We therefore conclude that inosine misincorporation into RNA perturbs translation, thus providing mechanistic insight linking ITPase deficiency, inosine accumulation and pathogenesis.

Molecular characterization of myotonic dystrophy fibroblast cell lines for use in small molecule screening.

Myotonic dystrophy type 1 (DM1) and type 2 (DM2) are common forms of adult onset muscular dystrophy. Pathogenesis in both diseases is largely driven by production of toxic-expanded repeat RNAs that sequester MBNL RNA-binding proteins, causing mis-splicing. Given this shared pathogenesis, we hypothesized that diamidines, small molecules that rescue mis-splicing in DM1 models, could also rescue mis-splicing in DM2 models. While several DM1 cell models exist, few are available for DM2 limiting research and therapeutic development. Here, we characterize DM1 and DM2 patient-derived fibroblasts for use in small molecule screens and therapeutic studies. We identify mis-splicing events unique to DM2 fibroblasts and common events shared with DM1 fibroblasts. We show that diamidines can partially rescue molecular phenotypes in both DM1 and DM2 fibroblasts. This study demonstrates the potential of fibroblasts as models for DM1 and DM2, which will help meet an important need for well-characterized DM2 cell models.

CCG•CGG interruptions in high-penetrance SCA8 families increase RAN translation and protein toxicity.

Spinocerebellar ataxia type 8 (SCA8), a dominantly inherited neurodegenerative disorder caused by a CTG•CAG expansion, is unusual because most individuals that carry the mutation do not develop ataxia. To understand the variable penetrance of SCA8, we studied the molecular differences between highly penetrant families and more common sporadic cases (82%) using a large cohort of SCA8 families (n = 77). We show that repeat expansion mutations from individuals with multiple affected family members have CCG•CGG interruptions at a higher frequency than sporadic SCA8 cases and that the number of CCG•CGG interruptions correlates with age at onset. At the molecular level, CCG•CGG interruptions increase RNA hairpin stability, and in cell culture experiments, increase p-eIF2α and polyAla and polySer RAN protein levels. Additionally, CCG•CGG interruptions, which encode arginine interruptions in the polyGln frame, increase toxicity of the resulting proteins. In summary, SCA8 CCG•CGG interruptions increase polyAla and polySer RAN protein levels, polyGln protein toxicity, and disease penetrance and provide novel insight into the molecular differences between SCA8 families with high vs. low disease penetrance.

Zebrafish mbnl mutants model physical and molecular phenotypes of myotonic dystrophy.

The muscleblind RNA-binding proteins (MBNL1, MBNL2 and MBNL3) are highly conserved across vertebrates and are important regulators of RNA alternative splicing. Loss of MBNL protein function through sequestration by CUG or CCUG RNA repeats is largely responsible for the phenotypes of the human genetic disorder myotonic dystrophy (DM). We generated the first stable zebrafish (Danio rerio) models of DM-associated MBNL loss of function through mutation of the three zebrafish mbnl genes. In contrast to mouse models, zebrafish double and triple homozygous mbnl mutants were viable to adulthood. Zebrafish mbnl mutants displayed disease-relevant physical phenotypes including decreased body size and impaired movement. They also exhibited widespread alternative splicing changes, including the misregulation of many DM-relevant exons. Physical and molecular phenotypes were more severe in compound mbnl mutants than in single mbnl mutants, suggesting partially redundant functions of Mbnl proteins. The high fecundity and larval optical transparency of this complete series of zebrafish mbnl mutants will make them useful for studying DM-related phenotypes and how individual Mbnl proteins contribute to them, and for testing potential therapeutics. This article has an associated First Person interview with the first author of the paper.

Repeat length increases disease penetrance and severity in C9orf72 ALS/FTD BAC transgenic mice.

C9orf72 ALS/FTD patients show remarkable clinical heterogeneity, but the complex biology of the repeat expansion mutation has limited our understanding of the disease. BAC transgenic mice were used to better understand the molecular mechanisms and repeat length effects of C9orf72 ALS/FTD. Genetic analyses of these mice demonstrate that the BAC transgene and not integration site effects cause ALS/FTD phenotypes. Transcriptomic changes in cell proliferation, inflammation and neuronal pathways are found late in disease and alternative splicing changes provide early molecular markers that worsen with disease progression. Isogenic sublines of mice with 800, 500 or 50 G4C2 repeats generated from the single-copy C9-500 line show longer repeats result in earlier onset, increased disease penetrance and increased levels of RNA foci and dipeptide RAN protein aggregates. These data demonstrate G4C2 repeat length is an important driver of disease and identify alternative splicing changes as early biomarkers of C9orf72 ALS/FTD.

RNA structure probing to characterize RNA-protein interactions on a low abundance pre-mRNA in living cells.

In vivo RNA structure analysis has become a powerful tool in molecular biology, largely due to the coupling of an increasingly diverse set of chemical approaches with high-throughput sequencing. This has resulted in a transition from single target to transcriptome-wide approaches. However, these methods require sequencing depths that preclude studying low abundance targets, which are not sufficiently captured in transcriptome-wide approaches. Here we present a ligation-free method to enrich for low abundance RNA sequences, which improves the diversity of molecules analyzed and results in improved analysis. In addition, this method is compatible with any choice of chemical adduct or read-out approach. We utilized this approach to study an autoregulated event in the pre-mRNA of the splicing factor, muscleblind-like splicing regulator 1 (MBNL1).

Asian Zika Virus Isolate Significantly Changes the Transcriptional Profile and Alternative RNA Splicing Events in a Neuroblastoma Cell Line.

The alternative splicing of pre-mRNAs expands a single genetic blueprint to encode multiple, functionally diverse protein isoforms. Viruses have previously been shown to interact with, depend on, and alter host splicing machinery. The consequences, however, incited by viral infection on the global alternative slicing (AS) landscape are under-appreciated. Here, we investigated the transcriptional and alternative splicing profile of neuronal cells infected with a contemporary Puerto Rican Zika virus (ZIKVPR) isolate, an isolate of the prototypical Ugandan ZIKV (ZIKVMR), and dengue virus 2 (DENV2). Our analyses revealed that ZIKVPR induced significantly more differential changes in expressed genes compared to ZIKVMR or DENV2, despite all three viruses showing equivalent infectivity and viral RNA levels. Consistent with the transcriptional profile, ZIKVPR induced a higher number of alternative splicing events compared to ZIKVMR or DENV2, and gene ontology analyses highlighted alternative splicing changes in genes associated with mRNA splicing. In summary, we show that ZIKV affects cellular RNA homeostasis not only at the transcriptional levels but also through the alternative splicing of cellular transcripts. These findings could provide new molecular insights into the neuropathologies associated with this virus.

Drug Screen Tugs at Common Thread for Repeat Disorders.

Repeat-associated non-ATG (RAN) translation is emerging as a driver of pathogenesis in microsatellite expansion disorders. Green and colleagues recently identified several candidate RAN translation inhibitors from a high-throughput small-molecule screen for fragile X tremor ataxia syndrome. Their study establishes a path forward for identifying inhibitors of RAN translation for multiple disorders.

The potential of engineered eukaryotic RNA-binding proteins as molecular tools and therapeutics.

Eukaroytic RNA-binding proteins (RBPs) recognize and process RNAs through recognition of their sequence motifs via RNA-binding domains (RBDs). RBPs usually consist of one or more RBDs and can include additional functional domains that modify or cleave RNA. Engineered RBPs have been used to answer basic biology questions, control gene expression, locate viral RNA in vivo, as well as many other tasks. Given the growing number of diseases associated with RNA and RBPs, engineered RBPs also have the potential to serve as therapeutics. This review provides an in depth description of recent advances in engineered RBPs and discusses opportunities and challenges in the field. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Methods > RNA Nanotechnology RNA in Disease and Development > RNA in Disease.

A CTG repeat-selective chemical screen identifies microtubule inhibitors as selective modulators of toxic CUG RNA levels.

A CTG repeat expansion in the DMPK gene is the causative mutation of myotonic dystrophy type 1 (DM1). Transcription of the expanded CTG repeat produces toxic gain-of-function CUG RNA, leading to disease symptoms. A screening platform that targets production or stability of the toxic CUG RNA in a selective manner has the potential to provide new biological and therapeutic insights. A DM1 HeLa cell model was generated that stably expresses a toxic r(CUG)480 and an analogous r(CUG)0 control from DMPK and was used to measure the ratio-metric level of r(CUG)480 versus r(CUG)0. This DM1 HeLa model recapitulates pathogenic hallmarks of DM1, including CUG ribonuclear foci and missplicing of pre-mRNA targets of the muscleblind (MBNL) alternative splicing factors. Repeat-selective screening using this cell line led to the unexpected identification of multiple microtubule inhibitors as hits that selectively reduce r(CUG)480 levels and partially rescue MBNL-dependent missplicing. These results were validated by using the Food and Drug Administration-approved clinical microtubule inhibitor colchicine in DM1 mouse and primary patient cell models. The mechanism of action was found to involve selective reduced transcription of the CTG expansion that we hypothesize to involve the LINC (linker of nucleoskeleton and cytoskeleton) complex. The unanticipated identification of microtubule inhibitors as selective modulators of toxic CUG RNA opens research directions for this form of muscular dystrophy and may shed light on the biology of CTG repeat expansion and inform therapeutic avenues. This approach has the potential to identify modulators of expanded repeat-containing gene expression for over 30 microsatellite expansion disorders.

Combination Treatment of Erythromycin and Furamidine Provides Additive and Synergistic Rescue of Mis-Splicing in Myotonic Dystrophy Type 1 Models.

Myotonic dystrophy type 1 (DM1) is a multi-systemic disease that presents with clinical symptoms including myotonia, cardiac dysfunction and cognitive impairment. DM1 is caused by a CTG expansion in the 3' UTR of the DMPK gene. The transcribed expanded CUG repeat RNA sequester the muscleblind-like (MBNL) and up-regulate the CUG-BP Elav-like (CELF) families of RNA-binding proteins leading to global mis-regulation of RNA processing and altered gene expression. Currently, there are no disease-targeting treatments for DM1. Given the multi-step pathogenic mechanism, combination therapies targeting different aspects of the disease mechanism may be a viable therapeutic approach. Here, as proof-of-concept, we studied a combination of two previously characterized small molecules, erythromycin and furamidine, in two DM1 models. In DM1 patient-derived myotubes, rescue of mis-splicing was observed with little to no cell toxicity. In a DM1 mouse model, a combination of erythromycin and the prodrug of furamidine (pafuramidine), administered orally, displayed both additive and synergistic mis-splicing rescue. Gene expression was only modestly affected and over 40 % of the genes showing significant expression changes were rescued back toward WT expression levels. Further, the combination treatment partially rescued the myotonia phenotype in the DM1 mouse. This combination treatment showed a high degree of mis-splicing rescue coupled with low off-target gene expression changes. These results indicate that combination therapies are a promising therapeutic approach for DM1.

Mitigating RNA Toxicity in Myotonic Dystrophy using Small Molecules.

This review, one in a series on myotonic dystrophy (DM), is focused on the development and potential use of small molecules as therapeutics for DM. The complex mechanisms and pathogenesis of DM are covered in the associated reviews. Here, we examine the various small molecule approaches taken to target the DNA, RNA, and proteins that contribute to disease onset and progression in myotonic dystrophy type 1 (DM1) and 2 (DM2).

Repeat-associated RNA structure and aberrant splicing.

Over 30 hereditary disorders attributed to the expansion of microsatellite repeats have been identified. Despite variant nucleotide content, number of consecutive repeats, and different locations in the genome, many of these diseases have pathogenic RNA gain-of-function mechanisms. The repeat-containing RNAs can form structures in vitro predicted to contribute to the disease through assembly of intracellular RNA aggregates termed foci. The expanded repeat RNAs within these foci sequester RNA binding proteins (RBPs) with important roles in the regulation of RNA metabolism, most notably alternative splicing (AS). These deleterious interactions lead to downstream alterations in transcriptome-wide AS directly linked with disease symptoms. This review summarizes existing knowledge about the association between the repeat RNA structures and RBPs as well as the resulting aberrant AS patterns, specifically in the context of myotonic dystrophy. The connection between toxic, structured RNAs and dysregulation of AS in other repeat expansion diseases is also discussed. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Transcriptome alterations in myotonic dystrophy skeletal muscle and heart.

Myotonic dystrophy (dystrophia myotonica, DM) is a multi-systemic disease caused by expanded CTG or CCTG microsatellite repeats. Characterized by symptoms in muscle, heart and central nervous system, among others, it is one of the most variable diseases known. A major pathogenic event in DM is the sequestration of muscleblind-like proteins by CUG or CCUG repeat-containing RNAs transcribed from expanded repeats, and differences in the extent of MBNL sequestration dependent on repeat length and expression level may account for some portion of the variability. However, many other cellular pathways are reported to be perturbed in DM, and the severity of specific disease symptoms varies among individuals. To help understand this variability and facilitate research into DM, we generated 120 RNASeq transcriptomes from skeletal and heart muscle derived from healthy and DM1 biopsies and autopsies. A limited number of DM2 and Duchenne muscular dystrophy samples were also sequenced. We analyzed splicing and gene expression, identified tissue-specific changes in RNA processing and uncovered transcriptome changes strongly correlating with muscle strength. We created a web resource at http://DMseq.org that hosts raw and processed transcriptome data and provides a lightweight, responsive interface that enables browsing of processed data across the genome.

Furamidine Rescues Myotonic Dystrophy Type I Associated Mis-Splicing through Multiple Mechanisms.

Myotonic dystrophy type 1 (DM1) is an autosomal dominant, CTG•CAG microsatellite expansion disease. Expanded CUG repeat RNA sequester the muscleblind-like (MBNL) family of RNA-binding proteins, thereby disrupting their normal cellular function which leads to global mis-regulation of RNA processing. Previously, the small molecule furamidine was shown to reduce CUG foci and rescue mis-splicing in a DM1 HeLa cell model and to rescue mis-splicing in the HSALR DM1 mouse model, but furamidine's mechanism of action was not explored. Here we use a combination of biochemical, cell toxicity, and genomic studies in DM1 patient-derived myotubes and the HSALR DM1 mouse model to investigate furamidine's mechanism of action. Mis-splicing rescue was observed in DM1 myotubes and the HSALR DM1 mouse with furamidine treatment. Interestingly, while furamidine was found to bind CTG•CAG repeat DNA with nanomolar affinity, a reduction in expanded CUG repeat transcript levels was observed in the HSALR DM1 mouse but not DM1 patient-derived myotubes. Further investigation in these cells revealed that furamidine treatment at nanomolar concentrations led to up-regulation of MBNL1 and MBNL2 protein levels and a reduction of ribonuclear foci. Additionally, furamidine was shown to bind CUG RNA with nanomolar affinity and disrupted the MBNL1 -CUG RNA complex in vitro at micromolar concentrations. Furamidine's likely promiscuous interactions in vitro and in vivo appear to affect multiple pathways in the DM1 mechanism to rescue mis-splicing, yet surprisingly furamidine was shown globally to rescue many mis-splicing events with only modest off-target effects on gene expression in the HSALR DM1 mouse model. Importantly, over 20% of the differentially expressed genes were shown to be returned, to varying degrees, to wild-type expression levels.