Splicing is a key process that expands the coding capacity of genomes. Its kinetics remain poorly characterized, and the distribution of splicing time caused by the stochasticity of single splicing events is expected to affect regulation efficiency. We conducted a small-scale survey on 40 introns in human cells and observed that most were spliced cotranscriptionally. Consequently, we constructed a reporter system that splices cotranscriptionally and can be monitored in live cells and in real time through the use of MS2-GFP. All small nuclear ribonucleoproteins (snRNPs) are loaded on nascent pre-mRNAs, and spliceostatin A inhibits splicing but not snRNP recruitment. Intron removal occurs in minutes and is best described by a model where several successive steps are rate limiting. Each pre-mRNA molecule is predicted to require a similar time to splice, reducing kinetic noise and improving the regulation of alternative splicing. This model is relevant to other kinetically controlled processes acting on few molecules.

Alternative splicing is the main mechanism of increasing the proteome diversity coded by a limited number of genes. It is well established that different tissues or organs express different splicing variants. However, organs are composed of common major cell types, including fibroblasts, epithelial, and endothelial cells. By analyzing large-scale data sets generated by The ENCODE Project Consortium and after extensive RT-PCR validation, we demonstrate that each of the three major cell types expresses a specific splicing program independently of its organ origin. Furthermore, by analyzing splicing factor expression across samples, publicly available splicing factor binding site data sets (CLIP-seq), and exon array data sets after splicing factor depletion, we identified several splicing factors, including ESRP1 and 2, MBNL1, NOVA1, PTBP1, and RBFOX2, that contribute to establishing these cell type-specific splicing programs. All of the analyzed data sets are freely available in a user-friendly web interface named FasterDB, which describes all known splicing variants of human and mouse genes and their splicing patterns across several dozens of normal and cancer cells as well as across tissues. Information regarding splicing factors that potentially contribute to individual exon regulation is also provided via a dedicated CLIP-seq and exon array data visualization interface. To the best of our knowledge, FasterDB is the first database integrating such a variety of large-scale data sets to enable functional genomics analyses at exon-level resolution.

Genomic studies and high-throughput experiments often produce large lists of candidate genes among which only a small fraction are truly relevant to the disease, phenotype or biological process of interest. Gene prioritization tackles this problem by ranking candidate genes by profiling candidates across multiple genomic data sources and integrating this heterogeneous information into a global ranking. We describe an extended version of our gene prioritization method, Endeavour, now available for six species and integrating 75 data sources. The performance (Area Under the Curve) of Endeavour on cross-validation benchmarks using 'gold standard' gene sets varies from 88% (for human phenotypes) to 95% (for worm gene function). In addition, we have also validated our approach using a time-stamped benchmark derived from the Human Phenotype Ontology, which provides a setting close to prospective validation. With this benchmark, using 3854 novel gene-phenotype associations, we observe a performance of 82%. Altogether, our results indicate that this extended version of Endeavour efficiently prioritizes candidate genes. The Endeavour web server is freely available at https://endeavour.esat.kuleuven.be/.

Myotonic dystrophy (DM) is caused by the expression of mutant RNAs containing expanded CUG repeats that sequester muscleblind-like (MBNL) proteins, leading to alternative splicing changes. Cardiac alterations, characterized by conduction delays and arrhythmia, are the second most common cause of death in DM. Using RNA sequencing, here we identify novel splicing alterations in DM heart samples, including a switch from adult exon 6B towards fetal exon 6A in the cardiac sodium channel, SCN5A. We find that MBNL1 regulates alternative splicing of SCN5A mRNA and that the splicing variant of SCN5A produced in DM presents a reduced excitability compared with the control adult isoform. Importantly, reproducing splicing alteration of Scn5a in mice is sufficient to promote heart arrhythmia and cardiac-conduction delay, two predominant features of myotonic dystrophy. In conclusion, misregulation of the alternative splicing of SCN5A may contribute to a subset of the cardiac dysfunctions observed in myotonic dystrophy.

Human pluripotent stem cells have ushered in an exciting new era for disease modeling, drug discovery, and cell therapy development. Continued progress toward realizing the potential of human pluripotent stem cells will be facilitated by robust data sets and complementary resources that are easily accessed and interrogated by the stem cell community. In this context, we present SISTEMA, a quality-controlled curated gene expression database, built on a valuable catalog of human pluripotent stem cell lines, and their derivatives for which transcriptomic analyses have been generated using a single experimental pipeline. SISTEMA functions as a one-step resource that will assist the stem cell community to easily evaluate the expression level for genes of interest, while comparing them across different hPSC lines, cell types, pathological conditions, or after pharmacological treatments.

Despite the initial efficacy of using tyrosine kinase inhibitors of epidermal growth factor receptors (EGFR-TKIs) for treating patients with non-small cell lung cancer (NSCLC), resistance inevitably develops. Recent studies highlight a link between alternative splicing and cancer drug response. Therefore, we aimed to identify deregulated splicing events that play a role in resistance to EGFR-TKI. By using RNA sequencing, reverse-transcription PCR (RT-PCR), and RNA interference, we showed that overexpression of a splice variant of the autophagic gene ATG16-L1 that retains exon 8 and encodes the β-isoform of autophagy-related protein 16-1 (ATG16-L1 β) concurs acquired resistance to EGFR-TKI in NSCLC cells. Using matched biopsies, we found increased levels of ATG16-L1 β at the time of progression in 3 of 11 NSCLC patients treated with EGFR-TKI. Mechanistically, gefitinib-induced autophagy was impaired in resistant cells that accumulated ATG16-L1 β. Neutralization of ATG16-L1 β restored autophagy in response to gefitinib, induced apoptosis, and inhibited the growth of in ovo tumor xenografts. Conversely, overexpression of ATG16-L1 β in parental sensitive cells prevented gefitinib-induced autophagy and increased cell survival. These results support a role of defective autophagy in acquired resistance to EGFR-TKIs and identify splicing regulation of ATG16-L1 as a therapeutic vulnerability that could be explored for improving EGFR-targeted cancer therapy.

The cytidine analogues azacytidine and 5-aza-2'-deoxycytidine (decitabine) are commonly used to treat myelodysplastic syndromes, with or without a myeloproliferative component. It remains unclear whether the response to these hypomethylating agents results from a cytotoxic or an epigenetic effect. In this study, we address this question in chronic myelomonocytic leukaemia. We describe a comprehensive analysis of the mutational landscape of these tumours, combining whole-exome and whole-genome sequencing. We identify an average of 14±5 somatic mutations in coding sequences of sorted monocyte DNA and the signatures of three mutational processes. Serial sequencing demonstrates that the response to hypomethylating agents is associated with changes in DNA methylation and gene expression, without any decrease in the mutation allele burden, nor prevention of new genetic alteration occurence. Our findings indicate that cytosine analogues restore a balanced haematopoiesis without decreasing the size of the mutated clone, arguing for a predominantly epigenetic effect.

Myotonic dystrophy type 1 (DM1) is a neuromuscular disorder caused by a non-coding CTG repeat expansion in the DMPK gene. This mutation generates a toxic CUG RNA that interferes with the RNA processing of target genes in multiple tissues. Despite debilitating neurological impairment, the pathophysiological cascade of molecular and cellular events in the central nervous system (CNS) has been less extensively characterized than the molecular pathogenesis of muscle/cardiac dysfunction. Particularly, the contribution of different cell types to DM1 brain disease is not clearly understood. We first used transcriptomics to compare the impact of expanded CUG RNA on the transcriptome of primary neurons, astrocytes and oligodendrocytes derived from DMSXL mice, a transgenic model of DM1. RNA sequencing revealed more frequent expression and splicing changes in glia than neuronal cells. In particular, primary DMSXL oligodendrocytes showed the highest number of transcripts differentially expressed, while DMSXL astrocytes displayed the most severe splicing dysregulation. Interestingly, the expression and splicing defects of DMSXL glia recreated molecular signatures suggestive of impaired cell differentiation: while DMSXL oligodendrocytes failed to upregulate a subset of genes that are naturally activated during the oligodendroglia differentiation, a significant proportion of missplicing events in DMSXL oligodendrocytes and astrocytes increased the expression of RNA isoforms typical of precursor cell stages. Together these data suggest that expanded CUG RNA in glial cells affects preferentially differentiation-regulated molecular events. This hypothesis was corroborated by gene ontology (GO) analyses, which revealed an enrichment for biological processes and cellular components with critical roles during cell differentiation. Finally, we combined exon ontology with phosphoproteomics and cell imaging to explore the functional impact of CUG-associated spliceopathy on downstream protein metabolism. Changes in phosphorylation, protein isoform expression and intracellular localization in DMSXL astrocytes demonstrate the far-reaching impact of the DM1 repeat expansion on cell metabolism. Our multi-omics approaches provide insight into the mechanisms of CUG RNA toxicity in the CNS with cell type resolution, and support the priority for future research on non-neuronal mechanisms and proteomic changes in DM1 brain disease.

5-Fluorouracil (5-FU) is a widely used chemotherapeutic drug in colorectal cancer. Previous studies showed that 5-FU modulates RNA metabolism and mRNA expression. In addition, it has been reported that 5-FU incorporates into the RNAs constituting the translational machinery and that 5-FU affects the amount of some mRNAs associated with ribosomes. However, the impact of 5-FU on translational regulation remains unclear. Using translatome profiling, we report that a clinically relevant dose of 5-FU induces a translational reprogramming in colorectal cancer cell lines. Comparison of mRNA distribution between polysomal and non-polysomal fractions in response to 5-FU treatment using microarray quantification identified 313 genes whose translation was selectively regulated. These regulations were mostly stimulatory (91%). Among these genes, we showed that 5-FU increases the mRNA translation of HIVEP2, which encodes a transcription factor whose translation in normal condition is known to be inhibited by mir-155. In response to 5-FU, the expression of mir-155 decreases thus stimulating the translation of HIVEP2 mRNA. Interestingly, the 5-FU-induced increase in specific mRNA translation was associated with reduction of global protein synthesis. Altogether, these findings indicate that 5-FU promotes a translational reprogramming leading to the increased translation of a subset of mRNAs that involves at least for some of them, miRNA-dependent mechanisms. This study supports a still poorly evaluated role of translational control in drug response.

Genome-wide analyses estimate that more than 90% of multi exonic human genes produce at least two transcripts through alternative splicing (AS). Various bioinformatics methods are available to analyze AS from RNAseq data. Most methods start by mapping the reads to an annotated reference genome, but some start by a de novo assembly of the reads. In this paper, we present a systematic comparison of a mapping-first approach (FARLINE) and an assembly-first approach (KISSPLICE). We applied these methods to two independent RNAseq datasets and found that the predictions of the two pipelines overlapped (70% of exon skipping events were common), but with noticeable differences. The assembly-first approach allowed to find more novel variants, including novel unannotated exons and splice sites. It also predicted AS in recently duplicated genes. The mapping-first approach allowed to find more lowly expressed splicing variants, and splice variants overlapping repeats. This work demonstrates that annotating AS with a single approach leads to missing out a large number of candidates, many of which are differentially regulated across conditions and can be validated experimentally. We therefore advocate for the combined use of both mapping-first and assembly-first approaches for the annotation and differential analysis of AS from RNAseq datasets.

Brain dysfunction in myotonic dystrophy type 1 (DM1), the prototype of toxic RNA disorders, has been mainly attributed to neuronal RNA misprocessing, while little attention has been given to non-neuronal brain cells. Here, using a transgenic mouse model of DM1 that expresses mutant RNA in various brain cell types (neurons, astroglia, and oligodendroglia), we demonstrate that astrocytes exhibit impaired ramification and polarization in vivo and defects in adhesion, spreading, and migration. RNA-dependent toxicity and phenotypes are also found in human transfected glial cells. In line with the cell phenotypes, molecular analyses reveal extensive expression and accumulation of toxic RNA in astrocytes, which result in RNA spliceopathy that is more severe than in neurons. Astrocyte missplicing affects primarily transcripts that regulate cell adhesion, cytoskeleton, and morphogenesis, and it is confirmed in human brain tissue. Our findings demonstrate that DM1 impacts astrocyte cell biology, possibly compromising their support and regulation of synaptic function.

It is widely accepted that pre-mRNA maturation, including splicing, is tightly coupled to both transcription and mRNA export, but factors linking the three processes are less understood. By analysing the estrogen-regulated expression of the c-fos mRNA that is processed during transcription, we show that the ddx5 RNA helicase, is required throughout the major nuclear steps of the expression of the c-fos gene, from transcription to mRNA export. Indeed, ddx5, whose recruitment on the c-fos gene was increased upon estrogen treatment, was required for the full transcriptional activation of the c-fos gene. In addition, ddx5 was required for c-fos co-transcriptional RNA splicing. When splicing occurred post-transcriptionally in the absence of ddx5, the c-fos mRNA was poorly exported into the cytosol because of inefficient recruitment of the TAP mRNA export receptor. Finally, ddx5 was present in the c-fos messenger ribonucleoprotein together with mRNA export factors, which further supports that ddx5 is a key operator in the c-fos 'mRNA factory'.

Stringent regulation of the interferon (IFN) signalling pathway is essential for maintaining the immune response to pathogens and tumours. The transcription factor STAT1 is a crucial mediator of this response. Here, we show that hCAF1/CNOT7 regulates class I and II IFN pathways at different crucial steps. In resting cells, hCAF1 can control STAT1 trafficking by interacting with the latent form of STAT1 in the cytoplasm. IFN treatment induces STAT1 release, suggesting that hCAF1 may shield cytoplasmic STAT1 from undesirable stimulation. Consistently, hCAF1 silencing enhances STAT1 basal promoter occupancy associated with increased expression of a subset of STAT1-regulated genes. Consequently, hCAF1 knockdown cells exhibit an increased protection against viral infection and reduced viral replication. Furthermore, hCAF1 participates in the extinction of the IFN signal, through its deadenylase activity, by speeding up the degradation of some STAT1-regulated mRNAs. Since abnormal and unbalanced JAK/STAT activation is associated with immune disorders and cancer, hCAF1 could play a major role in innate immunity and oncogenesis, contributing to tumour escape.

Nonsense-mediated mRNA Decay (NMD) degrades mutant mRNAs containing premature termination codon (PTC-mRNAs). Here we evaluate the consequence of NMD activity in colorectal cancers (CRCs) showing microsatellite instability (MSI) whose progression is associated with the accumulation of PTC-mRNAs encoding immunogenic proteins due to frameshift mutations in coding repeat sequences. Inhibition of UPF1, one of the major NMD factors, was achieved by siRNA in the HCT116 MSI CRC cell line and the resulting changes in gene expression were studied using expression microarrays. The impact of NMD activity was also investigated in primary MSI CRCs by quantifying the expression of several mRNAs relative to their mutational status and to endogenous UPF1 and UPF2 expression. Host immunity developed against MSI cancer cells was appreciated by quantifying the number of CD3epsilon-positive tumor-infiltrating lymphocytes (TILs). UPF1 silencing led to the up-regulation of 1251 genes in HCT116, among which a proportion of them (i.e. 38%) significantly higher than expected by chance contained a coding microsatellite (P<2x10(-16)). In MSI primary CRCs, UPF1 was significantly over-expressed compared to normal adjacent mucosa (P<0.002). Our data provided evidence for differential decay of PTC-mRNAs compared to wild-type that was positively correlated to UPF1 endogenous expression level (P = 0.02). A negative effect of UPF1 and UPF2 expression on the host's anti-tumor response was observed (P<0.01). Overall, our results show that NMD deeply influences MSI-driven tumorigenesis at the molecular level and indicate a functional negative impact of this system on anti-tumor immunity whose intensity has been recurrently shown to be an independent factor of favorable outcome in CRCs.

In AML, approximately one-third of expressed genes are abnormally spliced, including aberrant TET2 exon 2 expression. In a discovery cohort (n=99), TET2 exon 2 skipping (TET2E2S) was found positively associated with a significant reduction in the cumulative incidence of relapse (CIR). Age, cytogenetics, and TET2E2S were independent prognostic factors for disease-free survival (DFS), and favorable effects on outcomes predominated in cytogenetic normal (CN)-AML and younger patients. Using the same cutoff in a validation cohort of 86 CN-AML patients, TET2E2Shigh patients were found to be younger than TET2low patients without a difference in the rate of complete remission. However, TET2E2Shigh patients exhibited a significantly lower CIR (p<10-4). TET2E2S and FLT3-ITD, but not age or NPM1 mutation status were independent prognostic factors for DFS and event-free survival (EFS), while TET2E2S was the sole prognostic factor that we identified for overall survival (OS). In both the intermediate-1 and favorable ELN genetic categories, TET2E2S remained significantly associated with prolonged survival. There was no correlation between TET2E2S status and outcomes in 34 additional AML patients who were unfit for IC. Therefore our results suggest that assessments of TET2 exon 2 splicing status might improve risk stratification in CN-AML patients treated with IC.

Transcriptomic genome-wide analyses demonstrate massive variation of alternative splicing in many physiological and pathological situations. One major challenge is now to establish the biological contribution of alternative splicing variation in physiological- or pathological-associated cellular phenotypes. Toward this end, we developed a computational approach, named "Exon Ontology," based on terms corresponding to well-characterized protein features organized in an ontology tree. Exon Ontology is conceptually similar to Gene Ontology-based approaches but focuses on exon-encoded protein features instead of gene level functional annotations. Exon Ontology describes the protein features encoded by a selected list of exons and looks for potential Exon Ontology term enrichment. By applying this strategy to exons that are differentially spliced between epithelial and mesenchymal cells and after extensive experimental validation, we demonstrate that Exon Ontology provides support to discover specific protein features regulated by alternative splicing. We also show that Exon Ontology helps to unravel biological processes that depend on suites of coregulated alternative exons, as we uncovered a role of epithelial cell-enriched splicing factors in the AKT signaling pathway and of mesenchymal cell-enriched splicing factors in driving splicing events impacting on autophagy. Freely available on the web, Exon Ontology is the first computational resource that allows getting a quick insight into the protein features encoded by alternative exons and investigating whether coregulated exons contain the same biological information.

The Repressor Element 1-silencing transcription factor (REST) represses a number of neuronal genes in non-neuronal cells or in undifferentiated neural progenitors. Here, we report that the DEAD box RNA helicase DDX17 controls important REST-related processes that are critical during the early phases of neuronal differentiation. First, DDX17 associates with REST, promotes its binding to the promoter of a subset of REST-targeted genes and co-regulates REST transcriptional repression activity. During neuronal differentiation, we observed a downregulation of DDX17 along with that of the REST complex that contributes to the activation of neuronal genes. Second, DDX17 and its paralog DDX5 regulate the expression of several proneural microRNAs that are known to target the REST complex during neurogenesis, including miR-26a/b that are also direct regulators of DDX17 expression. In this context, we propose a new mechanism by which RNA helicases can control the biogenesis of intronic miRNAs. We show that the processing of the miR-26a2 precursor is dependent on RNA helicases, owing to an intronic regulatory region that negatively impacts on both miRNA processing and splicing of its host intron. Our work places DDX17 in the heart of a pathway involving REST and miRNAs that allows neuronal gene repression.

TGF-β1 is involved in many aspects of tissue development and homeostasis including hematopoiesis. The TAL1 transcription factor is also an important player of this latter process and is expressed very early in the myeloid and erythroid lineages. We previously established a link between TGF-β1 signaling and TAL1 by showing that the cytokine was able to induce its proteolytic degradation by the ubiquitin proteasome pathway. In this manuscript we show that TAL1 interacts with SMAD3 that acts in the pathway downstream of TGF-β1 association with its receptor. TAL1 expression strengthens the positive or negative effect of SMAD3 on various genes. Both transcription factors activate the inhibitory SMAD7 factor through the E box motif present in its transcriptional promoter. DNA precipitation assays showed that TAL1 present in Jurkat or K562 cells binds to this SMAD binding element in a SMAD3 dependent manner. SMAD3 and TAL1 also inhibit several genes including ID1, hTERT and TGF-β1 itself. In this latter case TAL1 and SMAD3 can impair the positive effect exerted by E47. Our results indicate that TAL1 expression can modulate TGF-β1 signaling by interacting with SMAD3 and by increasing its transcriptional properties. They also suggest the existence of a negative feedback loop between TAL1 expression and TGF-β1 signaling.

Chronic NF-κB activation in inflammation and cancer has long been linked to persistent activation of NF-κB-responsive gene promoters. However, NF-κB factors also massively bind to gene bodies. Here, we demonstrate that recruitment of the NF-κB factor RELA to intragenic regions regulates alternative splicing upon NF-κB activation by the viral oncogene Tax of HTLV-1. Integrative analyses of RNA splicing and chromatin occupancy, combined with chromatin tethering assays, demonstrate that DNA-bound RELA interacts with and recruits the splicing regulator DDX17, in an NF-κB activation-dependent manner. This leads to alternative splicing of target exons due to the RNA helicase activity of DDX17. Similar results were obtained upon Tax-independent NF-κB activation, indicating that Tax likely exacerbates a physiological process where RELA provides splice target specificity. Collectively, our results demonstrate a physical and direct involvement of NF-κB in alternative splicing regulation, which significantly revisits our knowledge of HTLV-1 pathogenesis and other NF-κB-related diseases.

The NF-κB protein p65/RelA plays a pivotal role in coordinating gene expression in response to diverse stimuli, including viral infections. At the chromatin level, p65/RelA regulates gene transcription and alternative splicing through promoter enrichment and genomic exon occupancy, respectively. The intricate ways in which p65/RelA simultaneously governs these functions across various genes remain to be fully elucidated. In this study, we employed the HTLV-1 Tax oncoprotein, a potent activator of NF-κB, to investigate its influence on the three-dimensional organization of the genome, a key factor in gene regulation. We discovered that Tax restructures the 3D genomic landscape, bringing together genes based on their regulation and splicing patterns. Notably, we found that the Tax-induced gene-gene contact between the two master genes NFKBIA and RELA is associated with their respective changes in gene expression and alternative splicing. Through dCas9-mediated approaches, we demonstrated that NFKBIA-RELA interaction is required for alternative splicing regulation and is caused by an intragenic enrichment of p65/RelA on RELA. Our findings shed light on new regulatory mechanisms upon HTLV-1 Tax and underscore the integral role of p65/RelA in coordinated regulation of NF-κB-responsive genes at both transcriptional and splicing levels in the context of the 3D genome.