Deficiency of the minor spliceosome component U4atac snRNA secondarily results in ciliary defects in human and zebrafish.

In the human genome, about 750 genes contain one intron excised by the minor spliceosome. This spliceosome comprises its own set of snRNAs, among which U4atac. Its noncoding gene, RNU4ATAC, has been found mutated in Taybi-Linder (TALS/microcephalic osteodysplastic primordial dwarfism type 1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes. These rare developmental disorders, whose physiopathological mechanisms remain unsolved, associate ante- and post-natal growth retardation, microcephaly, skeletal dysplasia, intellectual disability, retinal dystrophy, and immunodeficiency. Here, we report bi-allelic RNU4ATAC mutations in five patients presenting with traits suggestive of the Joubert syndrome (JBTS), a well-characterized ciliopathy. These patients also present with traits typical of TALS/RFMN/LWS, thus widening the clinical spectrum of RNU4ATAC-associated disorders and indicating ciliary dysfunction as a mechanism downstream of minor splicing defects. Intriguingly, all five patients carry the n.16G>A mutation, in the Stem II domain, either at the homozygous or compound heterozygous state. A gene ontology term enrichment analysis on minor intron-containing genes reveals that the cilium assembly process is over-represented, with no less than 86 cilium-related genes containing at least one minor intron, among which there are 23 ciliopathy-related genes. The link between RNU4ATAC mutations and ciliopathy traits is supported by alterations of primary cilium function in TALS and JBTS-like patient fibroblasts, as well as by u4atac zebrafish model, which exhibits ciliopathy-related phenotypes and ciliary defects. These phenotypes could be rescued by WT but not by pathogenic variants-carrying human U4atac. Altogether, our data indicate that alteration of cilium biogenesis is part of the physiopathological mechanisms of TALS/RFMN/LWS, secondarily to defects of minor intron splicing.

Altered splicing of ATG16-L1 mediates acquired resistance to tyrosine kinase inhibitors of EGFR by blocking autophagy in non-small cell lung cancer.

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.

Clinical interpretation of variants identified in RNU4ATAC, a non-coding spliceosomal gene.

Biallelic variants in RNU4ATAC, a non-coding gene transcribed into the minor spliceosome component U4atac snRNA, are responsible for three rare recessive developmental diseases, namely Taybi-Linder/MOPD1, Roifman and Lowry-Wood syndromes. Next-generation sequencing of clinically heterogeneous cohorts (children with either a suspected genetic disorder or a congenital microcephaly) recently identified mutations in this gene, illustrating how profoundly these technologies are modifying genetic testing and assessment. As RNU4ATAC has a single non-coding exon, the bioinformatic prediction algorithms assessing the effect of sequence variants on splicing or protein function are irrelevant, which makes variant interpretation challenging to molecular diagnostic laboratories. In order to facilitate and improve clinical diagnostic assessment and genetic counseling, we present i) an update of the previously reported RNU4ATAC mutations and an analysis of the genetic variations affecting this gene using the Genome Aggregation Database (gnomAD) resource; ii) the pathogenicity prediction performances of scores computed based on an RNA structure prediction tool and of those produced by the Combined Annotation Dependent Depletion tool for the 285 RNU4ATAC variants identified in patients or in large-scale sequencing projects; iii) a method, based on a cellular assay, that allows to measure the effect of RNU4ATAC variants on splicing efficiency of a minor (U12-type) reporter intron. Lastly, the concordance of bioinformatic predictions and cellular assay results was investigated.

ZRANB2 and SYF2-mediated splicing programs converging on ECT2 are involved in breast cancer cell resistance to doxorubicin.

Besides analyses of specific alternative splicing (AS) variants, little is known about AS regulatory pathways and programs involved in anticancer drug resistance. Doxorubicin is widely used in breast cancer chemotherapy. Here, we identified 1723 AS events and 41 splicing factors regulated in a breast cancer cell model of acquired resistance to doxorubicin. An RNAi screen on splicing factors identified the little studied ZRANB2 and SYF2, whose depletion partially reversed doxorubicin resistance. By RNAi and RNA-seq in resistant cells, we found that the AS programs controlled by ZRANB2 and SYF2 were enriched in resistance-associated AS events, and converged on the ECT2 splice variant including exon 5 (ECT2-Ex5+). Both ZRANB2 and SYF2 were found associated with ECT2 pre-messenger RNA, and ECT2-Ex5+ isoform depletion reduced doxorubicin resistance. Following doxorubicin treatment, resistant cells accumulated in S phase, which partially depended on ZRANB2, SYF2 and the ECT2-Ex5+ isoform. Finally, doxorubicin combination with an oligonucleotide inhibiting ECT2-Ex5 inclusion reduced doxorubicin-resistant tumor growth in mouse xenografts, and high ECT2-Ex5 inclusion levels were associated with bad prognosis in breast cancer treated with chemotherapy. Altogether, our data identify AS programs controlled by ZRANB2 and SYF2 and converging on ECT2, that participate to breast cancer cell resistance to doxorubicin.

New insights into minor splicing-a transcriptomic analysis of cells derived from TALS patients.

Minor intron splicing plays a central role in human embryonic development and survival. Indeed, biallelic mutations in RNU4ATAC, transcribed into the minor spliceosomal U4atac snRNA, are responsible for three rare autosomal recessive multimalformation disorders named Taybi-Linder (TALS/MOPD1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes, which associate numerous overlapping signs of varying severity. Although RNA-seq experiments have been conducted on a few RFMN patient cells, none have been performed in TALS, and more generally no in-depth transcriptomic analysis of the ∼700 human genes containing a minor (U12-type) intron had been published as yet. We thus sequenced RNA from cells derived from five skin, three amniotic fluid, and one blood biosamples obtained from seven unrelated TALS cases and from age- and sex-matched controls. This allowed us to describe for the first time the mRNA expression and splicing profile of genes containing U12-type introns, in the context of a functional minor spliceosome. Concerning RNU4ATAC-mutated patients, we show that as expected, they display distinct U12-type intron splicing profiles compared to controls, but that rather unexpectedly mRNA expression levels are mostly unchanged. Furthermore, although U12-type intron missplicing concerns most of the expressed U12 genes, the level of U12-type intron retention is surprisingly low in fibroblasts and amniocytes, and much more pronounced in blood cells. Interestingly, we found several occurrences of introns that can be spliced using either U2, U12, or a combination of both types of splice site consensus sequences, with a shift towards splicing using preferentially U2 sites in TALS patients' cells compared to controls.

The RNA helicase DDX17 controls the transcriptional activity of REST and the expression of proneural microRNAs in neuronal differentiation.

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.

Identification of protein features encoded by alternative exons using Exon Ontology.

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.