Mutations causing aberrant splicing are frequently implicated in human diseases including cancer. Here, we establish a high-throughput screen of randomly mutated minigenes to decode the cis-regulatory landscape that determines alternative splicing of exon 11 in the proto-oncogene MST1R (RON). Mathematical modelling of splicing kinetics enables us to identify more than 1000 mutations affecting RON exon 11 skipping, which corresponds to the pathological isoform RON∆165. Importantly, the effects correlate with RON alternative splicing in cancer patients bearing the same mutations. Moreover, we highlight heterogeneous nuclear ribonucleoprotein H (HNRNPH) as a key regulator of RON splicing in healthy tissues and cancer. Using iCLIP and synergy analysis, we pinpoint the functionally most relevant HNRNPH binding sites and demonstrate how cooperative HNRNPH binding facilitates a splicing switch of RON exon 11. Our results thereby offer insights into splicing regulation and the impact of mutations on alternative splicing in cancer.

Background Alternative polyadenylation (APA) refers to the regulated selection of polyadenylation sites (PASs) in transcripts, which determines the length of their 3′ untranslated regions (3′UTRs). We have recently shown that SRSF3 and SRSF7, two closely related SR proteins, connect APA with mRNA export. The mechanism underlying APA regulation by SRSF3 and SRSF7 remained unknown. Results Here we combine iCLIP and 3′-end sequencing and find that SRSF3 and SRSF7 bind upstream of proximal PASs (pPASs), but they exert opposite effects on 3′UTR length. SRSF7 enhances pPAS usage in a concentration-dependent but splicing-independent manner by recruiting the cleavage factor FIP1, generating short 3′UTRs. Protein domains unique to SRSF7, which are absent from SRSF3, contribute to FIP1 recruitment. In contrast, SRSF3 promotes distal PAS (dPAS) usage and hence long 3′UTRs directly by counteracting SRSF7, but also indirectly by maintaining high levels of cleavage factor Im (CFIm) via alternative splicing. Upon SRSF3 depletion, CFIm levels decrease and 3′UTRs are shortened. The indirect SRSF3 targets are particularly sensitive to low CFIm levels, because here CFIm serves a dual function; it enhances dPAS and inhibits pPAS usage by binding immediately downstream and assembling unproductive cleavage complexes, which together promotes long 3′UTRs. Conclusions We demonstrate that SRSF3 and SRSF7 are direct modulators of pPAS usage and show how small differences in the domain architecture of SR proteins can confer opposite effects on pPAS regulation.

Chromosomal translocations frequently promote carcinogenesis by producing gain-of-function fusion proteins. Recent studies have identified highly recurrent chromosomal translocations in patients with endometrial stromal sarcomas (ESSs) and ossifying fibromyxoid tumors (OFMTs), leading to an in-frame fusion of PHF1 (PCL1) to six different subunits of the NuA4/TIP60 complex. While NuA4/TIP60 is a coactivator that acetylates chromatin and loads the H2A.Z histone variant, PHF1 is part of the Polycomb repressive complex 2 (PRC2) linked to transcriptional repression of key developmental genes through methylation of histone H3 on lysine 27. In this study, we characterize the fusion protein produced by the EPC1-PHF1 translocation. The chimeric protein assembles a megacomplex harboring both NuA4/TIP60 and PRC2 activities and leads to mislocalization of chromatin marks in the genome, in particular over an entire topologically associating domain including part of the HOXD cluster. This is linked to aberrant gene expression—most notably increased expression of PRC2 target genes. Furthermore, we show that JAZF1—implicated with a PRC2 component in the most frequent translocation in ESSs, JAZF1-SUZ12—is a potent transcription activator that physically associates with NuA4/TIP60, its fusion creating outcomes similar to those of EPC1-PHF1. Importantly, the specific increased expression of PRC2 targets/HOX genes was also confirmed with ESS patient samples. Altogether, these results indicate that most chromosomal translocations linked to these sarcomas use the same molecular oncogenic mechanism through a physical merge of NuA4/TIP60 and PRC2 complexes, leading to mislocalization of histone marks and aberrant Polycomb target gene expression.

Alternative splicing is a key step in eukaryotic gene expression that allows for the production of multiple transcript and protein isoforms from the same gene. Even though splicing is perturbed in many diseases, we currently lack insights into regulatory mechanisms promoting its precision and efficiency. We analyze high-throughput mutagenesis data obtained for an alternatively spliced exon in the proto-oncogene RON and determine the functional units that control this splicing event. Using mathematical modeling of distinct splicing mechanisms, we show that alternative splicing is based in RON on a so-called ''exon definition'' mechanism. Here, the recognition of the adjacent exons by the spliceosome is required for removal of an intron. We use our model to analyze the differences between the exon and intron definition scenarios and find that exon definition prevents the accumulation of deleterious, partially spliced retention products during alternative splicing regulation. Furthermore, it modularizes splicing control, as multiple regulatory inputs are integrated into a common net input, irrespective of the location and nature of the corresponding cis-regulatory elements in the pre-messenger RNA. Our analysis suggests that exon definition promotes robust and reliable splicing outcomes in RON splicing.

Pyridoxine‐dependent epilepsy (PDE) is a relatively rare subgroup of epileptic disorders. They generally present in infancy as an early onset epileptic encephalopathy or seizures, refractory to standard treatments, with rapid and variable responses to vitamin B6 treatment. Whole exome sequencing of three unrelated families identified homozygous pathogenic mutation c.370_373del, p.Asp124fs in PLPBP gene in five persons. Haplotype analysis showed a single shared profile for the affected persons and their parents, leading to a hypothesis about founder effect of the mutation in Saguenay‐Lac‐St‐Jean region of French Canadians. All affected probands also shared one single mitochondrial haplotype T2b3 and two rare variations in the mitochondrial genome m.801A>G and m.5166A>G suggesting that a single individual female introduced PLPBP mutation c.370_373del, p.Asp124fs in Quebec. The mutation p.Asp124fs causes a severe disease phenotype with delayed myelination and cortical/subcortical brain atrophy. The most noteworthy radiological finding in this Quebec founder mutation is the presence of the temporal cysts that can be used as a marker of the disease. Also, both patients, who are alive, had a history of prenatal supplements taken by their mothers as antiemetic medication with high doses of pyridoxine. In the context of suspected PDE in patients with neonatal refractory seizures, treatment with pyridoxine and/or Pyridoxal‐5‐phophate has to be started immediately and continued until the results of genetic analysis received. Even with early appropriate treatment, neurological outcome of our patient is still poor.

Nemaline myopathy is a rare disorder affecting the muscle sarcomere. Mutations in nebulin gene (NEB) are known to be responsible for about 50% of nemaline myopathy cases. Nebulin is a giant protein which is formed integrally with the sarcomeric thin filament. This complex gene is under extensive alternative splicing giving rise to multiple isoforms. In this study, we report a 6-year-old boy presenting with general muscular weaknesses. Identification of rod-shaped structures in the patient' biopsy raised doubt about the presence of a nemaline myopathy. Next-generation sequencing was used to identify a causative mutation for the patient syndrome. A homozygous deep intronic substitution was found in the intron 144 of the NEB. The variant was predicted by in silico tools to create a new donor splice site. Molecular analysis has shown that the mutation could alter splicing events of the nebulin gene leading to a significant decrease of isoforms level. This change in the expression level of nebulin could give rise to functional consequences in the sarcomere. These results are consistent with the phenotypes observed in the patient. Such a discovery of variants in this gene will allow a better understanding of the involvement of nebulin in neuromuscular diseases and help find new treatments for the nemaline myopathy.

Rare diseases impact the lives of 300 million people in the world. Rapid advances in bioinformatics and genomic technologies have enabled the discovery of causes of 20–30% of rare diseases. However, most rare diseases have remained as unsolved enigmas to date. Newer tools and availability of high throughput sequencing data have enabled the reanalysis of previously undiagnosed patients. In this review, we have systematically compiled the latest developments in the discovery of the genetic causes of rare diseases using machine learning methods. Importantly, we have detailed methods available to reanalyze existing whole exome sequencing data of unsolved rare diseases. We have identified different reanalysis methodologies to solve problems associated with sequence alterations/mutations, variation re-annotation, protein stability, splice isoform malfunctions and oligogenic analysis. In addition, we give an overview of new developments in the field of rare disease research using whole genome sequencing data and other omics.

Alternative splicing is a key step in eukaryotic gene expression that allows the production of multiple protein isoforms from the same gene. Even though splicing is perturbed in many diseases, we currently lack insights into regulatory mechanisms promoting its precision and efficiency. Using mechanistic mathematical modeling, we show that alternative splicing control is facilitated if spliceosomes recognize exons as functional units ('exon definition'). We find that exon definition is crucial to prevent the accumulation of partially spliced retention products during alternative splicing regulation. Furthermore, it modularizes splicing control, as multiple regulatory inputs are integrated into a common net input, irrespective of the location and nature of the corresponding cis-regulatory elements in the pre-mRNA. These predictions of our model are qualitatively and quantitatively supported by high-throughput mutagenesis data obtained for an alternatively spliced exon in the proto-oncogene RON (MST1R). Our analysis suggests that exon definition has evolved as the dominant splice-regulatory mechanism in higher organisms to promote robust and reliable splicing outcomes.One Sentence Summary: Exon definition is required for precise alternative splicing control and prevents accumulation of undesired retention isoforms.