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.

Non-microRNA binding competitively inhibits LIN28 regulation.

RNA binding protein (RBP) expression is finite. For RBPs that are vastly outnumbered by their potential target sites, a simple competition for binding can set the magnitude of post-transcriptional control. Here, we show that LIN28, best known for its direct regulation of let-7 miRNA biogenesis, is also indirectly regulated by its widespread binding of non-miRNA transcripts. Approximately 99% of LIN28 binding sites are found on non-miRNA transcripts, like protein coding and ribosomal RNAs. These sites are bound specifically and strongly, but they do not appear to mediate direct post-transcriptional regulation. Instead, non-miRNA sites act to sequester LIN28 protein and effectively change its functional availability, thus impeding the regulation of let-7 in cells. Together, these data show that the binding properties of the transcriptome broadly influence the ability of an RBP to mediate changes in RNA metabolism and gene expression.

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.

A role for alternative splicing in circadian control of exocytosis and glucose homeostasis.

The circadian clock is encoded by a negative transcriptional feedback loop that coordinates physiology and behavior through molecular programs that remain incompletely understood. Here, we reveal rhythmic genome-wide alternative splicing (AS) of pre-mRNAs encoding regulators of peptidergic secretion within pancreatic ß cells that are perturbed in Clock-/- and Bmal1-/- ß-cell lines. We show that the RNA-binding protein THRAP3 (thyroid hormone receptor-associated protein 3) regulates circadian clock-dependent AS by binding to exons at coding sequences flanking exons that are more frequently skipped in clock mutant ß cells, including transcripts encoding Cask (calcium/calmodulin-dependent serine protein kinase) and Madd (MAP kinase-activating death domain). Depletion of THRAP3 restores expression of the long isoforms of Cask and Madd, and mimicking exon skipping in these transcripts through antisense oligonucleotide delivery in wild-type islets reduces glucose-stimulated insulin secretion. Finally, we identify shared networks of alternatively spliced exocytic genes from islets of rodent models of diet-induced obesity that significantly overlap with clock mutants. Our results establish a role for pre-mRNA alternative splicing in ß-cell function across the sleep/wake cycle.

Pooled CRISPR screens with imaging on microraft arrays reveals stress granule-regulatory factors.

Genetic screens using pooled CRISPR-based approaches are scalable and inexpensive, but restricted to standard readouts, including survival, proliferation and sortable markers. However, many biologically relevant cell states involve cellular and subcellular changes that are only accessible by microscopic visualization, and are currently impossible to screen with pooled methods. Here we combine pooled CRISPR-Cas9 screening with microraft array technology and high-content imaging to screen image-based phenotypes (CRaft-ID; CRISPR-based microRaft followed by guide RNA identification). By isolating microrafts that contain genetic clones harboring individual guide RNAs (gRNA), we identify RNA-binding proteins (RBPs) that influence the formation of stress granules, the punctate protein-RNA assemblies that form during stress. To automate hit identification, we developed a machine-learning model trained on nuclear morphology to remove unhealthy cells or imaging artifacts. In doing so, we identified and validated previously uncharacterized RBPs that modulate stress granule abundance, highlighting the applicability of our approach to facilitate image-based pooled CRISPR screens.

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.

A Transcriptome-wide Translational Program Defined by LIN28B Expression Level.

LIN28 RNA binding proteins are dynamically expressed throughout mammalian development and during disease. However, it remains unclear how changes in LIN28 expression define patterns of post-transcriptional gene regulation. Here we show that LIN28 expression level is a key variable that sets the magnitude of protein translation. By systematically varying LIN28B protein levels in human cells, we discovered a dose-dependent divergence in transcriptome-wide ribosome occupancy that enabled the formation of two discrete translational subpopulations composed of nearly all expressed genes. This bifurcation in gene expression was mediated by a redistribution in Argonaute association, from let-7 to non-let-7 microRNA families, resulting in a global shift in cellular miRNA activity. Post-transcriptional effects were scaled across the physiological LIN28 expression range. Together, these data highlight the central importance of RBP expression level and its ability to encode regulation.

Transcriptome-pathology correlation identifies interplay between TDP-43 and the expression of its kinase CK1E in sporadic ALS.

Sporadic amyotrophic lateral sclerosis (sALS) is the most common form of ALS, however, the molecular mechanisms underlying cellular damage and motor neuron degeneration remain elusive. To identify molecular signatures of sALS we performed genome-wide expression profiling in laser capture microdissection-enriched surviving motor neurons (MNs) from lumbar spinal cords of sALS patients with rostral onset and caudal progression. After correcting for immunological background, we discover a highly specific gene expression signature for sALS that is associated with phosphorylated TDP-43 (pTDP-43) pathology. Transcriptome-pathology correlation identified casein kinase 1ε (CSNK1E) mRNA as tightly correlated to levels of pTDP-43 in sALS patients. Enhanced crosslinking and immunoprecipitation in human sALS patient- and healthy control-derived frontal cortex, revealed that TDP-43 binds directly to and regulates the expression of CSNK1E mRNA. Additionally, we were able to show that pTDP-43 itself binds RNA. CK1E, the protein product of CSNK1E, in turn interacts with TDP-43 and promotes cytoplasmic accumulation of pTDP-43 in human stem-cell-derived MNs. Pathological TDP-43 phosphorylation is therefore, reciprocally regulated by CK1E activity and TDP-43 RNA binding. Our framework of transcriptome-pathology correlations identifies candidate genes with relevance to novel mechanisms of neurodegeneration.

Advances and challenges in the detection of transcriptome-wide protein-RNA interactions.

RNA binding proteins (RBPs) play key roles in determining cellular behavior by manipulating the processing of target RNAs. Robust methods are required to detect the numerous binding sites of RBPs across the transcriptome. RNA-immunoprecipitation followed by sequencing (RIP-seq) and crosslinking followed by immunoprecipitation and sequencing (CLIP-seq) are state-of-the-art methods used to identify the RNA targets and specific binding sites of RBPs. Historically, CLIP methods have been confounded with challenges such as the requirement for tens of millions of cells per experiment, low RNA yields resulting in libraries that contain a high number of polymerase chain reaction duplicated reads, and technical inconveniences such as radioactive labeling of RNAs. However, recent improvements in the recovery of bound RNAs and the efficiency of converting isolated RNAs into a library for sequencing have enhanced our ability to perform the experiment at scale, from less starting material than has previously been possible, and resulting in high quality datasets for the confident identification of protein binding sites. These, along with additional improvements to protein capture, removal of nonspecific signals, and methods to isolate noncanonical RBP targets have revolutionized the study of RNA processing regulation, and reveal a promising future for mapping the human protein-RNA regulatory network. WIREs RNA 2018, 9:e1436. doi: 10.1002/wrna.1436 This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Methods > RNA Analyses in Cells.

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.

RNA-binding protein CPEB1 remodels host and viral RNA landscapes.

Host and virus interactions occurring at the post-transcriptional level are critical for infection but remain poorly understood. Here, we performed comprehensive transcriptome-wide analyses revealing that human cytomegalovirus (HCMV) infection results in widespread alternative splicing (AS), shortening of 3' untranslated regions (3' UTRs) and lengthening of poly(A)-tails in host gene transcripts. We found that the host RNA-binding protein CPEB1 was highly induced after infection, and ectopic expression of CPEB1 in noninfected cells recapitulated infection-related post-transcriptional changes. CPEB1 was also required for poly(A)-tail lengthening of viral RNAs important for productive infection. Strikingly, depletion of CPEB1 reversed infection-related cytopathology and post-transcriptional changes, and decreased productive HCMV titers. Host RNA processing was also altered in herpes simplex virus-2 (HSV-2)-infected cells, thereby indicating that this phenomenon might be a common occurrence during herpesvirus infections. We anticipate that our work may serve as a starting point for therapeutic targeting of host RNA-binding proteins in herpesvirus infections.

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.