Regulating PCCA gene expression by modulation of pseudoexon splicing patterns to rescue enzyme activity in propionic acidemia.

Pseudoexons are nonfunctional intronic sequences that can be activated by deep-intronic sequence variation. Activation increases pseudoexon inclusion in mRNA and interferes with normal gene expression. The PCCA c.1285-1416A>G variation activates a pseudoexon and causes the severe metabolic disorder propionic acidemia by deficiency of the propionyl-CoA carboxylase enzyme encoded by PCCA and PCCB. We characterized this pathogenic pseudoexon activation event in detail and identified hnRNP A1 to be important for normal repression. The PCCA c.1285-1416A>G variation disrupts an hnRNP A1-binding splicing silencer and simultaneously creates a splicing enhancer. We demonstrate that blocking this region of regulation with splice-switching antisense oligonucleotides restores normal splicing and rescues enzyme activity in patient fibroblasts and in a cellular model created by CRISPR gene editing. Interestingly, the PCCA pseudoexon offers an unexploited potential to upregulate gene expression because healthy tissues show relatively high inclusion levels. By blocking inclusion of the nonactivated wild-type pseudoexon, we can increase both PCCA and PCCB protein levels, which increases the activity of the heterododecameric enzyme. Surprisingly, we can increase enzyme activity from residual levels in not only patient fibroblasts harboring PCCA missense variants but also those harboring PCCB missense variants. This is a potential treatment strategy for propionic acidemia.

Genome-wide epigenetic and mRNA-expression profiling followed by CRISPR/Cas9-mediated gene-disruptions corroborate the MIR141/MIR200C-ZEB1/ZEB2-FGFR1 axis in acquired EMT-associated EGFR TKI-resistance in NSCLC cells.

Background: Epithelial-mesenchymal-transition (EMT) is an epigenetic-based mechanism contributing to the acquired treatment resistance against receptor tyrosine kinase inhibitors (TKIs) in non-small cell lung cancer (NSCLC) cells harboring epidermal growth factor receptor (EGFR)-mutations. Delineating the exact epigenetic and gene-expression alterations in EMT-associated EGFR TKI-resistance (EMT-E-TKI-R) is vital for improved diagnosis and treatment of NSCLC patients. Methods: We characterized genome-wide changes in mRNA-expression, DNA-methylation and the histone-modification H3K36me3 in EGFR-mutated NSCLC HCC827 cells in result of acquired EMT-E-TKI-R. CRISPR/Cas9 was used to functional examine key findings from the omics analyses. Results: Acquired EMT-E-TKI-R was analyzed with three omics approaches. RNA-sequencing identified 2,233 and 1,972 up- and down-regulated genes, respectively, and among these were established EMT-markers. DNA-methylation EPIC array analyses identified 14,163 and 7,999 hyper- and hypo-methylated, respectively, differential methylated positions of which several were present in EMT-markers. Finally, H3K36me3 chromatin immunoprecipitation (ChIP)-sequencing detected 2,873 and 3,836 genes with enrichment and depletion, respectively, and among these were established EMT-markers. Correlation analyses showed that EMT-E-TKI-R mRNA-expression changes correlated better with H3K36me3 changes than with DNA-methylation changes. Moreover, the omics data supported the involvement of the MIR141/MIR200C-ZEB1/ZEB2-FGFR1 signaling axis for acquired EMT-E-TKI-R. CRISPR/Cas9-mediated analyses corroborated the importance of ZEB1 in acquired EMT-E-TKI-R, MIR200C and MIR141 to be in an EMT-E-TKI-R-associated auto-regulatory loop with ZEB1, and FGFR1 to mediate cell survival in EMT-E-TKI-R. Conclusions: The current study describes the synchronous genome-wide changes in mRNA-expression, DNA-methylation, and H3K36me3 in NSCLC EMT-E-TKI-R. The omics approaches revealed potential novel diagnostic markers and treatment targets. Besides, the study consolidates the functional impact of the MIR141/MIR200C-ZEB1/ZEB2-FGFR1-signaling axis in NSCLC EMT-E-TKI-R.

VEGFA-targeting miR-agshRNAs combine efficacy with specificity and safety for retinal gene therapy.

Retinal gene therapy using RNA interference (RNAi) to silence targeted genes requires both efficacy and safety. Short hairpin RNAs (shRNAs) are useful for RNAi, but high expression levels and activity from the co-delivered passenger strand may cause undesirable cellular responses. Ago2-dependent shRNAs (agshRNAs) produce no passenger strand activity. To enhance efficacy and to investigate improvements in safety, we have generated VEGFA-targeting agshRNAs and microRNA (miRNA)-embedded agshRNAs (miR-agshRNAs) and inserted these RNAi effectors in Pol II/III-driven expression cassettes and lentiviral vectors (LVs). Compared with corresponding shRNAs, agshRNAs and miR-agshRNAs increased specificity and safety, while retaining a high knockdown efficacy and abolishing passenger strand activity. The agshRNAs also caused significantly smaller reductions in cell viability and reduced competition with the processing of endogenous miR21 compared with their shRNA counterparts. RNA sequencing (RNA-seq) analysis of LV-transduced ARPE19 cells revealed that expression of shRNAs in general leads to more changes in gene expression levels compared with their agshRNA counterparts and activation of immune-related pathways. In mice, subretinal delivery of LVs encoding tissue-specific miR-agshRNAs resulted in retinal pigment epithelium (RPE)-restricted expression and significant knockdown of Vegfa in transduced RPE cells. Collectively, our data suggest that agshRNAs and miR-agshRNA possess important advantages over shRNAs, thereby posing a clinically relevant approach with respect to efficacy, specificity, and safety.

DeepCLIP: predicting the effect of mutations on protein-RNA binding with deep learning.

Nucleotide variants can cause functional changes by altering protein-RNA binding in various ways that are not easy to predict. This can affect processes such as splicing, nuclear shuttling, and stability of the transcript. Therefore, correct modeling of protein-RNA binding is critical when predicting the effects of sequence variations. Many RNA-binding proteins recognize a diverse set of motifs and binding is typically also dependent on the genomic context, making this task particularly challenging. Here, we present DeepCLIP, the first method for context-aware modeling and predicting protein binding to RNA nucleic acids using exclusively sequence data as input. We show that DeepCLIP outperforms existing methods for modeling RNA-protein binding. Importantly, we demonstrate that DeepCLIP predictions correlate with the functional outcomes of nucleotide variants in independent wet lab experiments. Furthermore, we show how DeepCLIP binding profiles can be used in the design of therapeutically relevant antisense oligonucleotides, and to uncover possible position-dependent regulation in a tissue-specific manner. DeepCLIP is freely available as a stand-alone application and as a webtool at http://deepclip.compbio.sdu.dk.

Tissue-resident macrophages in omentum promote metastatic spread of ovarian cancer.

Experimental and clinical evidence suggests that tumor-associated macrophages (TAMs) play important roles in cancer progression. Here, we have characterized the ontogeny and function of TAM subsets in a mouse model of metastatic ovarian cancer that is representative for visceral peritoneal metastasis. We show that the omentum is a critical premetastatic niche for development of invasive disease in this model and define a unique subset of CD163+ Tim4+ resident omental macrophages responsible for metastatic spread of ovarian cancer cells. Transcriptomic analysis showed that resident CD163+ Tim4+ omental macrophages were phenotypically distinct and maintained their resident identity during tumor growth. Selective depletion of CD163+ Tim4+ macrophages in omentum using genetic and pharmacological tools prevented tumor progression and metastatic spread of disease. These studies describe a specific role for tissue-resident macrophages in the invasive progression of metastatic ovarian cancer. The molecular pathways of cross-talk between tissue-resident macrophages and disseminated cancer cells may represent new targets to prevent metastasis and disease recurrence.

DFI-seq identification of environment-specific gene expression in uropathogenic Escherichia coli.

BACKGROUND: During infection of the urinary tract, uropathogenic Escherichia coli (UPEC) are exposed to different environments, such as human urine and the intracellular environments of bladder epithelial cells. Each environment elicits a distinct bacterial environment-specific transcriptional response. We combined differential fluorescence induction (DFI) with next-generation sequencing, collectively termed DFI-seq, to identify differentially expressed genes in UPEC strain UTI89 during growth in human urine and bladder cells. RESULTS: DFI-seq eliminates the need for iterative cell sorting of the bacterial library and yields a genome-wide view of gene expression. By analysing the gene expression of UPEC in human urine we found that genes involved in amino acid biosynthesis were upregulated. Deletion mutants lacking genes involved in arginine biosynthesis were outcompeted by the wild type during growth in human urine and inhibited in their ability to invade or proliferate in the J82 bladder epithelial cell line. Furthermore, DFI-seq was used to identify genes involved in invasion of J82 bladder epithelial cells. 56 genes were identified to be differentially expressed of which almost 60% encoded hypothetical proteins. One such gene UTI89_C5139, displayed increased adhesion and invasion of J82 cells when deleted from UPEC strain UTI89. CONCLUSIONS: We demonstrate the usefulness of DFI-seq for identification of genes required for optimal growth of UPEC in human urine, as well as potential virulence genes upregulated during infection of bladder cell culture. DFI-seq holds potential for the study of bacterial gene expression in live-animal infection systems. By linking fitness genes, such as those genes involved in amino acid biosynthesis, to virulence, this study contributes to our understanding of UPEC pathophysiology.

Splicing factor 1 modulates dietary restriction and TORC1 pathway longevity in C. elegans.

Ageing is driven by a loss of transcriptional and protein homeostasis and is the key risk factor for multiple chronic diseases. Interventions that attenuate or reverse systemic dysfunction associated with age therefore have the potential to reduce overall disease risk in the elderly. Precursor mRNA (pre-mRNA) splicing is a fundamental link between gene expression and the proteome, and deregulation of the splicing machinery is linked to several age-related chronic illnesses. However, the role of splicing homeostasis in healthy ageing remains unclear. Here we demonstrate that pre-mRNA splicing homeostasis is a biomarker and predictor of life expectancy in Caenorhabditis elegans. Using transcriptomics and in-depth splicing analysis in young and old animals fed ad libitum or subjected to dietary restriction, we find defects in global pre-mRNA splicing with age that are reduced by dietary restriction via splicing factor 1 (SFA-1; the C. elegans homologue of SF1, also known as branchpoint binding protein, BBP). We show that SFA-1 is specifically required for lifespan extension by dietary restriction and by modulation of the TORC1 pathway components AMPK, RAGA-1 and RSKS-1/S6 kinase. We also demonstrate that overexpression of SFA-1 is sufficient to extend lifespan. Together, these data demonstrate a role for RNA splicing homeostasis in dietary restriction longevity and suggest that modulation of specific spliceosome components may prolong healthy ageing.

RNA-sequencing of a mouse-model of spinal muscular atrophy reveals tissue-wide changes in splicing of U12-dependent introns.

Spinal Muscular Atrophy (SMA) is a neuromuscular disorder caused by insufficient levels of the Survival of Motor Neuron (SMN) protein. SMN is expressed ubiquitously and functions in RNA processing pathways that include trafficking of mRNA and assembly of snRNP complexes. Importantly, SMA severity is correlated with decreased snRNP assembly activity. In particular, the minor spliceosomal snRNPs are affected, and some U12-dependent introns have been reported to be aberrantly spliced in patient cells and animal models. SMA is characterized by loss of motor neurons, but the underlying mechanism is largely unknown. It is likely that aberrant splicing of genes expressed in motor neurons is involved in SMA pathogenesis, but increasing evidence indicates that pathologies also exist in other tissues. We present here a comprehensive RNA-seq study that covers multiple tissues in an SMA mouse model. We show elevated U12-intron retention in all examined tissues from SMA mice, and that U12-dependent intron retention is induced upon siRNA knock-down of SMN in HeLa cells. Furthermore, we show that retention of U12-dependent introns is mitigated by ASO treatment of SMA mice and that many transcriptional changes are reversed. Finally, we report on missplicing of several Ca2+ channel genes that may explain disrupted Ca2+ homeostasis in SMA and activation of Cdk5.

The prevalent deep intronic c. 639+919 G>A GLA mutation causes pseudoexon activation and Fabry disease by abolishing the binding of hnRNPA1 and hnRNP A2/B1 to a splicing silencer.

Fabry disease is an X-linked recessive inborn disorder of the glycosphingolipid metabolism, caused by total or partial deficiency of the lysosomal α-galactosidase A enzyme due to mutations in the GLA gene. The prevalent c.639+919 G>A mutation in GLA leads to pathogenic insertion of a 57bp pseudoexon sequence from intron 4, which is responsible for the cardiac variant phenotype. In this study we investigate the splicing regulatory mechanism leading to GLA pseudoexon activation. Splicing analysis of GLA minigenes revealed that pseudoexon activation is influenced by cell-type. We demonstrate that the wild-type sequence harbors an hnRNP A1 and hnRNP A2/B1-binding exonic splicing silencer (ESS) overlapping the 5'splice site (5'ss) that prevents pseudoexon inclusion. The c.639+919 G>A mutation disrupts this ESS allowing U1 snRNP recognition of the 5'ss. We show that the wild-type GLA 5'ss motif with the ESS is also able to inhibit inclusion of an unrelated pseudoexon in the FGB gene, and that also in the FGB context inactivation of the ESS by the c.639+919 G>A mutation causes pseudoexon activation, underscoring the universal nature of the ESS. Finally, we demonstrate that splice switching oligonucleotide (SSO) mediated blocking of the pseudoexon 3'ss and 5'ss effectively restores normal GLA splicing. This indicates that SSO based splicing correction may be a therapeutic alternative in the treatment of Fabry disease.

Splice-shifting oligonucleotide (SSO) mediated blocking of an exonic splicing enhancer (ESE) created by the prevalent c.903+469T>C MTRR mutation corrects splicing and restores enzyme activity in patient cells.

The prevalent c.903+469T>C mutation in MTRR causes the cblE type of homocystinuria by strengthening an SRSF1 binding site in an ESE leading to activation of a pseudoexon. We hypothesized that other splicing regulatory elements (SREs) are also critical for MTRR pseudoexon inclusion. We demonstrate that the MTRR pseudoexon is on the verge of being recognized and is therefore vulnerable to several point mutations that disrupt a fine-tuned balance between the different SREs. Normally, pseudoexon inclusion is suppressed by a hnRNP A1 binding exonic splicing silencer (ESS). When the c.903+469T>C mutation is present two ESEs abrogate the activity of the ESS and promote pseudoexon inclusion. Blocking the 3'splice site or the ESEs by SSOs is effective in restoring normal splicing of minigenes and endogenous MTRR transcripts in patient cells. By employing an SSO complementary to both ESEs, we were able to rescue MTRR enzymatic activity in patient cells to approximately 50% of that in controls. We show that several point mutations, individually, can activate a pseudoexon, illustrating that this mechanism can occur more frequently than previously expected. Moreover, we demonstrate that SSO blocking of critical ESEs is a promising strategy to treat the increasing number of activated pseudoexons.

Absence of an intron splicing silencer in porcine Smn1 intron 7 confers immunity to the exon skipping mutation in human SMN2.

Spinal Muscular Atrophy is caused by homozygous loss of SMN1. All patients retain at least one copy of SMN2 which produces an identical protein but at lower levels due to a silent mutation in exon 7 which results in predominant exclusion of the exon. Therapies targeting the splicing of SMN2 exon 7 have been in development for several years, and their efficacy has been measured using either in vitro cellular assays or in vivo small animal models such as mice. In this study we evaluated the potential for constructing a mini-pig animal model by introducing minimal changes in the endogenous porcine Smn1 gene to maintain the native genomic structure and regulation. We found that while a Smn2-like mutation can be introduced in the porcine Smn1 gene and can diminish the function of the ESE, it would not recapitulate the splicing pattern seen in human SMN2 due to absence of a functional ISS immediately downstream of exon 7. We investigated the ISS region and show here that the porcine ISS is inactive due to disruption of a proximal hnRNP A1 binding site, while a distal hnRNP A1 binding site remains functional but is unable to maintain the functionality of the ISS as a whole.

A synonymous polymorphic variation in ACADM exon 11 affects splicing efficiency and may affect fatty acid oxidation.

In recent studies combining genome-wide association and tandem-MS based metabolic profiling, a single-nucleotide polymorphism (SNP), rs211718C>T, located far upstream of the MCAD gene (ACADM) was found to be associated with serum concentrations of medium-chain acylcarnitines indicating improved beta-oxidation of medium-chain fatty acids. We examined the functional basis for this association and identified linkage between rs211718 and the intragenic synonymous polymorphic variant c.1161A>G in ACADM exon 11 (rs1061337). Employing minigene studies we show that the c.1161A allele is associated with exon 11 missplicing, and that the c.1161G allele corrects this missplicing. This may result in production of more full length MCAD protein from the c.1161G allele. Our analysis suggests that the improved splicing of the c.1161G allele is due to changes in the relative binding of splicing regulatory proteins SRSF1 and hnRNP A1. Using publicly available pre-aligned RNA-seq data, we find that the ACADM c.1161G allele is expressed at significantly higher levels than the c.1161A allele across different tissues. This supports that c.1161A>G is a functional SNP, which leads to higher MCAD expression, perhaps due to improved splicing. This study is a proof of principle that synonymous SNPs are not neutral. By changing the binding sites for splicing regulatory proteins they can have significant effects on pre-mRNA splicing and thus protein function. In addition, this study shows that for a sequence variation to have an effect, it might need to change the balance in the relative binding of positive and negative splicing factors.

CUGBP1 and MBNL1 preferentially bind to 3' UTRs and facilitate mRNA decay.

CUGBP1 and MBNL1 are developmentally regulated RNA-binding proteins that are causally associated with myotonic dystrophy type 1. We globally determined the in vivo RNA-binding sites of CUGBP1 and MBNL1. Interestingly, CUGBP1 and MBNL1 are both preferentially bound to 39 UTRs. Analysis of CUGBP1- and MBNL1-bound 39 UTRs demonstrated that both factors mediate accelerated mRNA decay and temporal profiles of expression arrays supported this. Role of CUGBP1 on accelerated mRNA decay has been previously reported, but the similar function of MBNL1 has not been reported to date. It is well established that CUGBP1 and MBNL1 regulate alternative splicing. Screening by exon array and validation by RT-PCR revealed position dependence of CUGBP1- and MBNL1-binding sites on the resulting alternative splicing pattern. This study suggests that regulation of CUGBP1 and MBNL1 is essential for accurate control of destabilization of a broad spectrum of mRNAs as well as of alternative splicing events.

SMN2 exon 7 splicing is inhibited by binding of hnRNP A1 to a common ESS motif that spans the 3' splice site.

Spinal Muscular Atrophy is caused by homozygous loss of SMN1 with phenotypic modulation by SMN2. SMN2 expresses only limited amounts of full-length transcript due to skipping of exon 7 caused by disruption of an SF2/ASF binding ESE. Additionally, hnRNP A1 has been reported to inhibit inclusion of SMN2 exon 7. We previously reported high similarity between the sequence spanning the 3' ss of SMN1 and SMN2 exon 7 and an hnRNP A1 binding ESS, which regulates MCAD exon 5 splicing. We show here that this 3' ss motif indeed functions as a crucial hnRNP A1 binding ESS, which inhibits inclusion of SMN1/2 exon 7 and is antagonized by the SMN1 ESE, but not by the inactive SMN2 sequence. Pull-down experiments revealed a specific interaction between hnRNP A1 and the 3' ss AG-dinucleotide, which could be disrupted by mutations shown to improve splicing in reporter minigenes. Genomic analyses revealed that in the human genome, 3' ss matching the SMN1/2 ESS motif region are much less abundant than 3' ss with a disrupted ESS motif. This indicates that this ESS may be a general splicing inhibitory motif, which binds hnRNP A1 and inhibits exon inclusion by binding to 3' ss harboring this ESS motif.

The deep intronic c.903+469T>C mutation in the MTRR gene creates an SF2/ASF binding exonic splicing enhancer, which leads to pseudoexon activation and causes the cblE type of homocystinuria.

Deep intronic mutations are often ignored as possible causes of human diseases. A deep intronic mutation in the MTRR gene, c.903+469T>C, is the most frequent mutation causing the cblE type of homocystinuria. It is well known to be associated with pre-mRNA mis-splicing, resulting in pseudoexon inclusion; however, the pathological mechanism remains unknown. We used minigenes to demonstrate that this mutation is the direct cause of MTRR pseudoexon inclusion, and that the pseudoexon is normally not recognized due to a suboptimal 5' splice site. Within the pseudoexon we identified an exonic splicing enhancer (ESE), which is activated by the mutation. Cotransfection and siRNA experiments showed that pseudoexon inclusion depends on the cellular amounts of SF2/ASF and in vitro RNA-binding assays showed dramatically increased SF2/ASF binding to the mutant MTRR ESE. The mutant MTRR ESE sequence is identical to an ESE of the alternatively spliced MST1R proto-oncogene, which suggests that this ESE could be frequently involved in splicing regulation. Our study conclusively demonstrates that an intronic single nucleotide change is sufficient to cause pseudoexon activation via creation of a functional ESE, which binds a specific splicing factor. We suggest that this mechanism may cause genetic disease much more frequently than previously reported.

Seemingly neutral polymorphic variants may confer immunity to splicing-inactivating mutations: a synonymous SNP in exon 5 of MCAD protects from deleterious mutations in a flanking exonic splicing enhancer.

The idea that point mutations in exons may affect splicing is intriguing and adds an additional layer of complexity when evaluating their possible effects. Even in the best-studied examples, the molecular mechanisms are not fully understood. Here, we use patient cells, model minigenes, and in vitro assays to show that a missense mutation in exon 5 of the medium-chain acyl-CoA dehydrogenase (MCAD) gene primarily causes exon skipping by inactivating a crucial exonic splicing enhancer (ESE), thus leading to loss of a functional protein and to MCAD deficiency. This ESE functions by antagonizing a juxtaposed exonic splicing silencer (ESS) and is necessary to define a suboptimal 3' splice site. Remarkably, a synonymous polymorphic variation in MCAD exon 5 inactivates the ESS, and, although this has no effect on splicing by itself, it makes splicing immune to deleterious mutations in the ESE. Furthermore, the region of MCAD exon 5 that harbors these elements is nearly identical to the exon 7 region of the survival of motor neuron (SMN) genes that contains the deleterious silent mutation in SMN2, indicating a very similar and finely tuned interplay between regulatory elements in these two genes. Our findings illustrate a mechanism for dramatic context-dependent effects of single-nucleotide polymorphisms on gene-expression regulation and show that it is essential that potential deleterious effects of mutations on splicing be evaluated in the context of the relevant haplotype.