Splicing modulators impair DNA damage response and induce killing of cohesin-mutant MDS and AML.

Splicing modulation is a promising treatment strategy pursued to date only in splicing factor-mutant cancers; however, its therapeutic potential is poorly understood outside of this context. Like splicing factors, genes encoding components of the cohesin complex are frequently mutated in cancer, including myelodysplastic syndromes (MDS) and secondary acute myeloid leukemia (AML), where they are associated with poor outcomes. Here, we showed that cohesin mutations are biomarkers of sensitivity to drugs targeting the splicing factor 3B subunit 1 (SF3B1) H3B-8800 and E-7107. We identified drug-induced alterations in splicing, and corresponding reduced gene expression, of a number of DNA repair genes, including BRCA1 and BRCA2, as the mechanism underlying this sensitivity in cell line models, primary patient samples and patient-derived xenograft (PDX) models of AML. We found that DNA damage repair genes are particularly sensitive to exon skipping induced by SF3B1 modulators due to their long length and large number of exons per transcript. Furthermore, we demonstrated that treatment of cohesin-mutant cells with SF3B1 modulators not only resulted in impaired DNA damage response and accumulation of DNA damage, but it sensitized cells to subsequent killing by poly(ADP-ribose) polymerase (PARP) inhibitors and chemotherapy and led to improved overall survival of PDX models of cohesin-mutant AML in vivo. Our findings expand the potential therapeutic benefits of SF3B1 splicing modulators to include cohesin-mutant MDS and AML.

Aberrant NOVA1 function disrupts alternative splicing in early stages of amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal disease characterized by aberrant alternative splicing (AS). Nuclear loss and cytoplasmic accumulation of the splicing factor TDP-43 in motor neurons (MN) are hallmarks of ALS at late stages of the disease. However, it is unknown if altered AS is present before TDP-43 pathology occurs. Here, we investigate altered AS and its origins in early stages of ALS using human induced pluripotent stem cell-derived motor neurons (MNs) from sporadic and familial ALS patients. We find high levels of the RNA-binding proteins NOVA1, NOVA2, and RBFOX2 in the insoluble protein fractions and observe that AS events in ALS-associated MNs are enriched for binding sites of these proteins. Our study points to an early disrupted function of NOVA1 that drives AS changes in a complex fashion, including events caused by a consistent loss of NOVA1 function. NOVA1 exhibits increased cytoplasmic protein levels in early stage MNs without TDP-43 pathology in ALS postmortem tissue. As nuclear TDP-43 protein level depletes, NOVA1 is reduced. Potential indications for a reduction of NOVA1 also came from mice over-expressing TDP-43 lacking its nuclear localization signal and iPSC-MN stressed with puromycin. This study highlights that additional RBP-RNA perturbations in ALS occur in parallel to TDP-43.

Genetic analysis of cancer drivers reveals cohesin and CTCF as suppressors of PD-L1.

Immune evasion is a significant contributor to tumor evolution, and the immunoinhibitory axis PD-1/PD-L1 is a frequent mechanism employed to escape tumor immune surveillance. To identify cancer drivers involved in immune evasion, we performed a CRISPR-Cas9 screen of tumor suppressor genes regulating the basal and interferon (IFN)-inducible cell surface levels of PD-L1. Multiple regulators of PD-L1 were identified, including IRF2, ARID2, KMT2D, and AAMP. We also identified CTCF and the cohesin complex proteins, known regulators of chromatin architecture and transcription, among the most potent negative regulators of PD-L1 cell surface expression. Additionally, loss of the cohesin subunit RAD21 was shown to up-regulate PD-L2 and MHC-I surface expression. PD-L1 and MHC-I suppression by cohesin were shown to be conserved in mammary epithelial and myeloid cells. Comprehensive examination of the transcriptional effect of STAG2 deficiency in epithelial and myeloid cells revealed an activation of strong IFN and NF-κB expression signatures. Inhibition of JAK-STAT or NF-κB pathways did not result in rescue of PD-L1 up-regulation in RAD21-deficient cells, suggesting more complex or combinatorial mechanisms at play. Discovery of the PD-L1 and IFN up-regulation in cohesin-mutant cells expands our understanding of the biology of cohesin-deficient cells as well as molecular regulation of the PD-L1 molecule.

Integrative RNA-omics Discovers GNAS Alternative Splicing as a Phenotypic Driver of Splicing Factor-Mutant Neoplasms.

Mutations in splicing factors (SF) are the predominant class of mutations in myelodysplastic syndrome (MDS), but convergent downstream disease drivers remain elusive. To identify common direct targets of missplicing by mutant U2AF1 and SRSF2, we performed RNA sequencing and enhanced version of the cross-linking and immunoprecipitation assay in human hematopoietic stem/progenitor cells derived from isogenic induced pluripotent stem cell (iPSC) models. Integrative analyses of alternative splicing and differential binding converged on a long isoform of GNAS (GNAS-L), promoted by both mutant factors. MDS population genetics, functional and biochemical analyses support that GNAS-L is a driver of MDS and encodes a hyperactive long form of the stimulatory G protein alpha subunit, Gαs-L, that activates ERK/MAPK signaling. SF-mutant MDS cells have activated ERK signaling and consequently are sensitive to MEK inhibitors. Our findings highlight an unexpected and unifying mechanism by which SRSF2 and U2AF1 mutations drive oncogenesis with potential therapeutic implications for MDS and other SF-mutant neoplasms. SIGNIFICANCE: SF mutations are disease-defining in MDS, but their critical effectors remain unknown. We discover the first direct target of convergent missplicing by mutant U2AF1 and SRSF2, a long GNAS isoform, which activates G protein and ERK/MAPK signaling, thereby driving MDS and rendering mutant cells sensitive to MEK inhibition. This article is highlighted in the In This Issue feature, p. 587.

Reintroduction of the archaic variant of NOVA1 in cortical organoids alters neurodevelopment.

The evolutionarily conserved splicing regulator neuro-oncological ventral antigen 1 (NOVA1) plays a key role in neural development and function. NOVA1 also includes a protein-coding difference between the modern human genome and Neanderthal and Denisovan genomes. To investigate the functional importance of an amino acid change in humans, we reintroduced the archaic allele into human induced pluripotent cells using genome editing and then followed their neural development through cortical organoids. This modification promoted slower development and higher surface complexity in cortical organoids with the archaic version of NOVA1 Moreover, levels of synaptic markers and synaptic protein coassociations correlated with altered electrophysiological properties in organoids expressing the archaic variant. Our results suggest that the human-specific substitution in NOVA1, which is exclusive to modern humans since divergence from Neanderthals, may have had functional consequences for our species' evolution.

5-Azacytidine Transiently Restores Dysregulated Erythroid Differentiation Gene Expression in TET2-Deficient Erythroleukemia Cells.

DNA methyltransferase inhibitors (DNMTI) like 5-Azacytidine (5-Aza) are the only disease-modifying drugs approved for the treatment of higher-risk myelodysplastic syndromes (MDS), however less than 50% of patients respond, and there are no predictors of response with clinical utility. Somatic mutations in the DNA methylation regulating gene tet-methylcytosine dioxygenase 2 (TET2) are associated with response to DNMTIs, however the mechanisms responsible for this association remain unknown. Using bisulfite padlock probes, mRNA sequencing, and hydroxymethylcytosine pull-down sequencing at several time points throughout 5-Aza treatment, we show that TET2 loss particularly influences DNA methylation (5mC) and hydroxymethylation (5hmC) patterns at erythroid gene enhancers and is associated with downregulation of erythroid gene expression in the human erythroleukemia cell line TF-1. 5-Aza disproportionately induces expression of these down-regulated genes in TET2KO cells and this effect is related to dynamic 5mC changes at erythroid gene enhancers after 5-Aza exposure. We identified differences in remethylation kinetics after 5-Aza exposure for several types of genomic regulatory elements, with distal enhancers exhibiting longer-lasting 5mC changes than other regions. This work highlights the role of 5mC and 5hmC dynamics at distal enhancers in regulating the expression of differentiation-associated gene signatures, and sheds light on how 5-Aza may be more effective in patients harboring TET2 mutations. IMPLICATIONS: TET2 loss in erythroleukemia cells induces hypermethylation and impaired expression of erythroid differentiation genes which can be specifically counteracted by 5-Azacytidine, providing a potential mechanism for the increased efficacy of 5-Aza in TET2-mutant patients with MDS. VISUAL OVERVIEW: http://mcr.aacrjournals.org/content/molcanres/19/3/451/F1.large.jpg.

Principles of RNA processing from analysis of enhanced CLIP maps for 150 RNA binding proteins.

BACKGROUND: A critical step in uncovering rules of RNA processing is to study the in vivo regulatory networks of RNA binding proteins (RBPs). Crosslinking and immunoprecipitation (CLIP) methods enable mapping RBP targets transcriptome-wide, but methodological differences present challenges to large-scale analysis across datasets. The development of enhanced CLIP (eCLIP) enabled the mapping of targets for 150 RBPs in K562 and HepG2, creating a unique resource of RBP interactomes profiled with a standardized methodology in the same cell types. RESULTS: Our analysis of 223 eCLIP datasets reveals a range of binding modalities, including highly resolved positioning around splicing signals and mRNA untranslated regions that associate with distinct RBP functions. Quantification of enrichment for repetitive and abundant multicopy elements reveals 70% of RBPs have enrichment for non-mRNA element classes, enables identification of novel ribosomal RNA processing factors and sites, and suggests that association with retrotransposable elements reflects multiple RBP mechanisms of action. Analysis of spliceosomal RBPs indicates that eCLIP resolves AQR association after intronic lariat formation, enabling identification of branch points with single-nucleotide resolution, and provides genome-wide validation for a branch point-based scanning model for 3' splice site recognition. Finally, we show that eCLIP peak co-occurrences across RBPs enable the discovery of novel co-interacting RBPs. CONCLUSIONS: This work reveals novel insights into RNA biology by integrated analysis of eCLIP profiling of 150 RBPs with distinct functions. Further, our quantification of both mRNA and other element association will enable further research to identify novel roles of RBPs in regulating RNA processing.

Disruption in A-to-I Editing Levels Affects C. elegans Development More Than a Complete Lack of Editing.

A-to-I RNA editing, catalyzed by ADAR proteins, is widespread in eukaryotic transcriptomes. Studies showed that, in C. elegans, ADR-2 can actively deaminate dsRNA, whereas ADR-1 cannot. Therefore, we set out to study the effect of each of the ADAR genes on the RNA editing process. We performed comprehensive phenotypic, transcriptomics, proteomics, and RNA binding screens on worms mutated in a single ADAR gene. We found that ADR-1 mutants exhibit more-severe phenotypes than ADR-2, and some of them are a result of non-editing functions of ADR-1. We also show that ADR-1 significantly binds edited genes and regulates mRNA expression, whereas the effect on protein levels is minor. In addition, ADR-1 primarily promotes editing by ADR-2 at the L4 stage of development. Our results suggest that ADR-1 has a significant role in the RNA editing process and in altering editing levels that affect RNA expression; loss of ADR-1 results in severe phenotypes.

The C. elegans neural editome reveals an ADAR target mRNA required for proper chemotaxis.

ADAR proteins alter gene expression both by catalyzing adenosine (A) to inosine (I) RNA editing and binding to regulatory elements in target RNAs. Loss of ADARs affects neuronal function in all animals studied to date. Caenorhabditis elegans lacking ADARs exhibit reduced chemotaxis, but the targets responsible for this phenotype remain unknown. To identify critical neural ADAR targets in C. elegans, we performed an unbiased assessment of the effects of ADR-2, the only A-to-I editing enzyme in C. elegans, on the neural transcriptome. Development and implementation of publicly available software, SAILOR, identified 7361 A-to-I editing events across the neural transcriptome. Intersecting the neural editome with adr-2 associated gene expression changes, revealed an edited mRNA, clec-41, whose neural expression is dependent on deamination. Restoring clec-41 expression in adr-2 deficient neural cells rescued the chemotaxis defect, providing the first evidence that neuronal phenotypes of ADAR mutants can be caused by altered gene expression.

Noncoding regions of C. elegans mRNA undergo selective adenosine to inosine deamination and contain a small number of editing sites per transcript.

ADARs (Adenosine deaminases that act on RNA) "edit" RNA by converting adenosines to inosines within double-stranded regions. The primary targets of ADARs are long duplexes present within noncoding regions of mRNAs, such as introns and 3' untranslated regions (UTRs). Because adenosine and inosine have different base-pairing properties, editing within these regions can alter splicing and recognition by small RNAs. However, despite numerous studies identifying multiple editing sites in these genomic regions, little is known about the extent to which editing sites co-occur on individual transcripts or the functional output of these combinatorial editing events. To begin to address these questions, we performed an ultra-deep sequencing analysis of 4 Caenorhabditis elegans 3' UTRs that are known ADAR targets. Synchronous editing events were determined for the long duplexes in vivo. Furthermore, the validity of each editing event was confirmed by sequencing the same regions of mRNA from worms that lack A-to-I editing. This analysis identified a large number of editing sites that can occur within each 3' UTR, but interestingly, each individual transcript contained only a small fraction of these A-to-I editing events. In addition, editing patterns were not random, indicating that an editing event can affect the efficiency of editing at subsequent adenosines. Furthermore, we identified specific sites that can be both positively and negatively correlated with additional sites leading to mutually exclusive editing patterns. These results suggest that editing in noncoding regions is selective and hyper-editing of cellular RNAs is rare.