Metabolism-dependent secondary effect of anti-MAPK cancer therapy on DNA repair.

Amino acid bioavailability impacts mRNA translation in a codon-dependent manner. Here, we report that the anti-cancer MAPK inhibitors (MAPKi) decrease the intracellular concentration of aspartate and glutamate in melanoma cells. This coincides with the accumulation of ribosomes on codons corresponding to these amino acids and triggers the translation-dependent degradation of mRNAs encoding aspartate- and glutamate-rich proteins, involved in DNA metabolism such as DNA replication and repair. Consequently, cells that survive MAPKi degrade aspartate and glutamate likely to generate energy, which simultaneously decreases their requirement for amino acids due to the downregulation of aspartate- and glutamate-rich proteins involved in cell proliferation. Concomitantly, the downregulation of aspartate- and glutamate-rich proteins involved in DNA repair increases DNA damage loads. Thus, DNA repair defects, and therefore mutations, are at least in part a secondary effect of the metabolic adaptation of cells exposed to MAPKi.

Translation-dependent and -independent mRNA decay occur through mutually exclusive pathways defined by ribosome density during T cell activation.

mRNA translation and decay are tightly interconnected processes both in the context of mRNA quality-control pathways and for the degradation of functional mRNAs. Cotranslational mRNA degradation through codon usage, ribosome collisions, and the recruitment of specific proteins to ribosomes is an important determinant of mRNA turnover. However, the extent to which translation-dependent mRNA decay (TDD) and translation-independent mRNA decay (TID) pathways participate in the degradation of mRNAs has not been studied yet. Here we describe a comprehensive analysis of basal and signal-induced TDD and TID in mouse primary CD4+ T cells. Our results indicate that most cellular transcripts are decayed to some extent in a translation-dependent manner. Our analysis further identifies the length of untranslated regions, the density of ribosomes, and GC3 content as important determinants of TDD magnitude. Consistently, all transcripts that undergo changes in ribosome density within their coding sequence upon T cell activation display a corresponding change in their TDD level. Moreover, we reveal a dynamic modulation in the relationship between GC3 content and TDD upon T cell activation, with a reversal in the impact of GC3- and AU3-rich codons. Altogether, our data show a strong and dynamic interconnection between mRNA translation and decay in mammalian primary cells.

Gene-to-gene coordinated regulation of transcription and alternative splicing by 3D chromatin remodeling upon NF-κB activation.

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.

The NF-κB pathway is essential for coordinating gene expression in response to various stimuli, including viral infections. Most studies have focused on the role of NF-κB in transcriptional regulation. In the present study, the impact of the potent NF-κB activator HTLV-1 Tax oncoprotein on the three-dimensional organization of the genome was investigated. Tax-mediated NF-κB activation was found to restructure the 3D genomic landscape in cells and to bring genes together in multigene complexes that are coordinately regulated either transcriptionally or through alternative splicing by NF-κB. Induced coordinate changes in transcription and alternative splicing included the two master genes of NF-κB pathway NFKBIA and RELA. The findings have significant implications for understanding cell fate determination and disease development associated with HTLV-1 infection, as well as chronic NF-κB activation in various human inflammatory diseases and cancer.

SplicingLore: a web resource for studying the regulation of cassette exons by human splicing factors.

One challenge faced by scientists from the alternative RNA splicing field is to decode the cooperative or antagonistic effects of splicing factors (SFs) to understand and eventually predict splicing outcomes on a genome-wide scale. In this manuscript, we introduce SplicingLore, an open-access database and web resource that help to fill this gap in a straightforward manner. The database contains a collection of RNA-sequencing-derived lists of alternative exons regulated by a total of 75 different SFs. All datasets were processed in a standardized manner, ensuring valid comparisons and correlation analyses. The user can easily retrieve a factor-specific set of differentially included exons from the database or provide a list of exons and search which SF(s) control(s) their inclusion. Our simple workflow is fast and easy to run, and it ensures a reliable calculation of correlation scores between the tested datasets. As a proof of concept, we predicted and experimentally validated a novel functional cooperation between the RNA helicases DDX17 and DDX5 and the heterogeneous nuclear ribonucleoprotein C (HNRNPC) protein. SplicingLore is available at https://splicinglore.ens-lyon.fr/. Database URL:  https://splicinglore.ens-lyon.fr/.

RNA helicase-dependent gene looping impacts messenger RNA processing.

DDX5 and DDX17 are DEAD-box RNA helicase paralogs which regulate several aspects of gene expression, especially transcription and splicing, through incompletely understood mechanisms. A transcriptome analysis of DDX5/DDX17-depleted human cells confirmed the large impact of these RNA helicases on splicing and revealed a widespread deregulation of 3' end processing. In silico analyses and experiments in cultured cells showed the binding and functional contribution of the genome organizing factor CTCF to chromatin sites at or near a subset of DDX5/DDX17-dependent exons that are characterized by a high GC content and a high density of RNA Polymerase II. We propose the existence of an RNA helicase-dependent relationship between CTCF and the dynamics of transcription across DNA and/or RNA structured regions, that contributes to the processing of internal and terminal exons. Moreover, local DDX5/DDX17-dependent chromatin loops spatially connect RNA helicase-regulated exons with their cognate promoter, and we provide the first direct evidence that de novo gene looping modifies alternative splicing and polyadenylation. Overall our findings uncover the impact of DDX5/DDX17-dependent chromatin folding on pre-messenger RNA processing.

Myotonic dystrophy RNA toxicity alters morphology, adhesion and migration of mouse and human astrocytes.

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.

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.

Modulation of alternative splicing during early infection of human primary B lymphocytes with Epstein-Barr virus (EBV): a novel function for the viral EBNA-LP protein.

Epstein-Barr virus (EBV) is a human herpesvirus associated with human cancers worldwide. Ex vivo, the virus efficiently infects resting human B lymphocytes and induces their continuous proliferation. This process is accompanied by a global reprogramming of cellular gene transcription. However, very little is known on the impact of EBV infection on the regulation of alternative splicing, a pivotal mechanism that plays an essential role in cell fate determination and is often deregulated in cancer. In this study, we have developed a systematic time-resolved analysis of cellular mRNA splice variant expression during EBV infection of resting B lymphocytes. Our results reveal that major modifications of alternative splice variant expression appear as early as day 1 post-infection and suggest that splicing regulation provides-besides transcription-an additional mechanism of gene expression regulation at the onset of B cell activation and proliferation. We also report a role for the viral proteins, EBNA2 and EBNA-LP, in the modulation of specific alternative splicing events and reveal a previously unknown function for EBNA-LP-together with the RBM4 splicing factor-in the alternative splicing regulation of two important modulators of cell proliferation and apoptosis respectively, NUMB and BCL-X.

SISTEMA: A large and standardized collection of transcriptome data sets for human pluripotent stem cell research.

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.

The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype.

The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.

Integrative Cell Type-Specific Multi-Omics Approaches Reveal Impaired Programs of Glial Cell Differentiation in Mouse Culture Models of DM1.

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.

Intragenic recruitment of NF-κB drives splicing modifications upon activation by the oncogene Tax of HTLV-1.

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.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces.

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

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.

Characterizing the interplay between gene nucleotide composition bias and splicing.

BACKGROUND: Nucleotide composition bias plays an important role in the 1D and 3D organization of the human genome. Here, we investigate the potential interplay between nucleotide composition bias and the regulation of exon recognition during splicing. RESULTS: By analyzing dozens of RNA-seq datasets, we identify two groups of splicing factors that activate either about 3200 GC-rich exons or about 4000 AT-rich exons. We show that splicing factor-dependent GC-rich exons have predicted RNA secondary structures at 5' ss and are dependent on U1 snRNP-associated proteins. In contrast, splicing factor-dependent AT-rich exons have a large number of decoy branch points, SF1- or U2AF2-binding sites and are dependent on U2 snRNP-associated proteins. Nucleotide composition bias also influences local chromatin organization, with consequences for exon recognition during splicing. Interestingly, the GC content of exons correlates with that of their hosting genes, isochores, and topologically associated domains. CONCLUSIONS: We propose that regional nucleotide composition bias over several dozens of kilobase pairs leaves a local footprint at the exon level and induces constraints during splicing that can be alleviated by local chromatin organization at the DNA level and recruitment of specific splicing factors at the RNA level. Therefore, nucleotide composition bias establishes a direct link between genome organization and local regulatory processes, like alternative splicing.

Interplay between coding and exonic splicing regulatory sequences.

The inclusion of exons during the splicing process depends on the binding of splicing factors to short low-complexity regulatory sequences. The relationship between exonic splicing regulatory sequences and coding sequences is still poorly understood. We demonstrate that exons that are coregulated by any given splicing factor share a similar nucleotide composition bias and preferentially code for amino acids with similar physicochemical properties because of the nonrandomness of the genetic code. Indeed, amino acids sharing similar physicochemical properties correspond to codons that have the same nucleotide composition bias. In particular, we uncover that the TRA2A and TRA2B splicing factors that bind to adenine-rich motifs promote the inclusion of adenine-rich exons coding preferentially for hydrophilic amino acids that correspond to adenine-rich codons. SRSF2 that binds guanine/cytosine-rich motifs promotes the inclusion of GC-rich exons coding preferentially for small amino acids, whereas SRSF3 that binds cytosine-rich motifs promotes the inclusion of exons coding preferentially for uncharged amino acids, like serine and threonine that can be phosphorylated. Finally, coregulated exons encoding amino acids with similar physicochemical properties correspond to specific protein features. In conclusion, the regulation of an exon by a splicing factor that relies on the affinity of this factor for specific nucleotide(s) is tightly interconnected with the exon-encoded physicochemical properties. We therefore uncover an unanticipated bidirectional interplay between the splicing regulatory process and its biological functional outcome.

Pluripotent Stem Cell-Based Drug Screening Reveals Cardiac Glycosides as Modulators of Myotonic Dystrophy Type 1.

There is currently no treatment for myotonic dystrophy type 1 (DM1), the most frequent myopathy of genetic origin. This progressive neuromuscular disease is caused by nuclear-retained RNAs containing expanded CUG repeats. These toxic RNAs alter the activities of RNA splicing factors, resulting in alternative splicing misregulation. By combining human mutated pluripotent stem cells and phenotypic drug screening, we revealed that cardiac glycosides act as modulators for both upstream nuclear aggregations of DMPK mRNAs and several downstream alternative mRNA splicing defects. However, these occurred at different drug concentration ranges. Similar biological effects were recorded in a DM1 mouse model. At the mechanistic level, we demonstrated that this effect was calcium dependent and was synergic with inhibition of the ERK pathway. These results further underscore the value of stem-cell-based assays for drug discovery in monogenic diseases.

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.

Complementarity of assembly-first and mapping-first approaches for alternative splicing annotation and differential analysis from RNAseq data.

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.