Recent evidence points to a role of chromatin in regulation of alternative pre-mRNA splicing (AS). In order to identify novel chromatin regulators of AS, we screened an RNAi library of chromatin proteins using a cell-based high-throughput in vivo assay. We identified a set of chromatin proteins that regulate AS. Using simultaneous genome-wide expression and AS analysis, we demonstrate distinct and non-overlapping functions of these chromatin modifiers on transcription and AS. Detailed mechanistic characterization of one dual function chromatin modifier, the H3K9 methyltransferase EHMT2 (G9a), identified VEGFA as a major chromatin-mediated AS target. Silencing of EHMT2, or its heterodimer partner EHMT1, affects AS by promoting exclusion of VEGFA exon 6a, but does not alter total VEGFA mRNA levels. The epigenetic regulatory mechanism of AS by EHMT2 involves an adaptor system consisting of the chromatin modulator HP1γ, which binds methylated H3K9 and recruits splicing regulator SRSF1. The epigenetic regulation of VEGFA is physiologically relevant since EHMT2 is transcriptionally induced in response to hypoxia and triggers concomitant changes in AS of VEGFA. These results characterize a novel epigenetic regulatory mechanism of AS and they demonstrate separate roles of epigenetic modifiers in transcription and alternative splicing.

SMARCB1 (also known as SNF5, INI1, and BAF47), a core subunit of the SWI/SNF (BAF) chromatin-remodeling complex, is inactivated in nearly all pediatric rhabdoid tumors. These aggressive cancers are among the most genomically stable, suggesting an epigenetic mechanism by which SMARCB1 loss drives transformation. Here we show that, despite having indistinguishable mutational landscapes, human rhabdoid tumors exhibit distinct enhancer H3K27ac signatures, which identify remnants of differentiation programs. We show that SMARCB1 is required for the integrity of SWI/SNF complexes and that its loss alters enhancer targeting-markedly impairing SWI/SNF binding to typical enhancers, particularly those required for differentiation, while maintaining SWI/SNF binding at super-enhancers. We show that these retained super-enhancers are essential for rhabdoid tumor survival, including some that are shared by all subtypes, such as SPRY1, and other lineage-specific super-enhancers, such as SOX2 in brain-derived rhabdoid tumors. Taken together, our findings identify a new chromatin-based epigenetic mechanism underlying the tumor-suppressive activity of SMARCB1.

The analysis of exonic DNA from prostate cancers has identified recurrently mutated genes, but the spectrum of genome-wide alterations has not been profiled extensively in this disease. We sequenced the genomes of 57 prostate tumors and matched normal tissues to characterize somatic alterations and to study how they accumulate during oncogenesis and progression. By modeling the genesis of genomic rearrangements, we identified abundant DNA translocations and deletions that arise in a highly interdependent manner. This phenomenon, which we term "chromoplexy," frequently accounts for the dysregulation of prostate cancer genes and appears to disrupt multiple cancer genes coordinately. Our modeling suggests that chromoplexy may induce considerable genomic derangement over relatively few events in prostate cancer and other neoplasms, supporting a model of punctuated cancer evolution. By characterizing the clonal hierarchy of genomic lesions in prostate tumors, we charted a path of oncogenic events along which chromoplexy may drive prostate carcinogenesis.

Evolved resistance to tyrosine kinase inhibitor (TKI)-targeted therapies remains a major clinical challenge. In epidermal growth factor receptor (EGFR) mutant non-small-cell lung cancer (NSCLC), failure of EGFR TKIs can result from both genetic and epigenetic mechanisms of acquired drug resistance. Widespread reports of histologic and gene expression changes consistent with an epithelial-to-mesenchymal transition (EMT) have been associated with initially surviving drug-tolerant persister cells, which can seed bona fide genetic mechanisms of resistance to EGFR TKIs. While therapeutic approaches targeting fully resistant cells, such as those harboring an EGFRT790M mutation, have been developed, a clinical strategy for preventing the emergence of persister cells remains elusive. Using mesenchymal cell lines derived from biopsies of patients who progressed on EGFR TKI as surrogates for persister populations, we performed whole-genome CRISPR screening and identified fibroblast growth factor receptor 1 (FGFR1) as the top target promoting survival of mesenchymal EGFR mutant cancers. Although numerous previous reports of FGFR signaling contributing to EGFR TKI resistance in vitro exist, the data have not yet been sufficiently compelling to instigate a clinical trial testing this hypothesis, nor has the role of FGFR in promoting the survival of persister cells been elucidated. In this study, we find that combining EGFR and FGFR inhibitors inhibited the survival and expansion of EGFR mutant drug-tolerant cells over long time periods, preventing the development of fully resistant cancers in multiple vitro models and in vivo. These results suggest that dual EGFR and FGFR blockade may be a promising clinical strategy for both preventing and overcoming EMT-associated acquired drug resistance and provide motivation for the clinical study of combined EGFR and FGFR inhibition in EGFR-mutated NSCLCs.

The soluble isoform of leptin receptor (sOb-R), secreted by the liver, regulates leptin bioavailability and bioactivity. Its reduced levels in diet-induced obesity (DIO) contribute to hyperleptinemia and leptin resistance, effects that are regulated by the endocannabinoid (eCB)/CB1R system. Here we show that pharmacological activation/blockade and genetic overexpression/deletion of hepatic CB1R modulates sOb-R levels and hepatic leptin resistance. Interestingly, peripheral CB1R blockade failed to reverse DIO-induced reduction of sOb-R levels, increased fat mass and dyslipidemia, and hepatic steatosis in mice lacking C/EBP homologous protein (CHOP), whereas direct activation of CB1R in wild-type hepatocytes reduced sOb-R levels in a CHOP-dependent manner. Moreover, CHOP stimulation increased sOb-R expression and release via a direct regulation of its promoter, while CHOP deletion reduced leptin sensitivity. Our findings highlight a novel molecular aspect by which the hepatic eCB/CB1R system is involved in the development of hepatic leptin resistance and in the regulation of sOb-R levels via CHOP.

Enhancer demethylation in leukemia has been shown to lead to overexpression of genes which promote cancer characteristics. The vascular endothelial growth factor A (VEGFA) enhancer, located 157 Kb downstream of its promoter, is demethylated in chronic myeloid leukemia (CML). VEGFA has several alternative splicing isoforms with different roles in cancer progression. Since transcription and splicing are coupled, we wondered whether VEGFA enhancer activity can also regulate the gene's alternative splicing to contribute to the pathology of CML. Our results show that mutating the VEGFA +157 enhancer promotes exclusion of exons 6a and 7 and activating the enhancer by tethering a chromatin activator has the opposite effect. In line with these results, CML patients present with high expression of +157 eRNA and inclusion of VEGFA exons 6a and 7. In addition, our results show that the positive regulator of RNAPII transcription elongation, CCNT2, binds VEGFA's promoter and enhancer, and its silencing promotes exclusion of exons 6a and 7 as it slows down RNAPII elongation rate. Thus our results suggest that VEGFA's +157 enhancer regulates its alternative splicing by increasing RNAPII elongation rate via CCNT2. Our work demonstrates for the first time a connection between an endogenous enhancer and alternative splicing regulation of its target gene.

Despite extensive studies on the effects of SARS-CoV-2 infection, there is still a lack of understanding of the downstream epigenetic and regulatory alterations in infected cells. In this study, we investigated changes in enhancer acetylation in epithelial lung cells infected with SARS-CoV-2 and their influence on transcriptional regulation and pathway activity. To achieve this, we integrated and reanalyzed data of enhancer acetylation, ex-vivo infection and single cell RNA-seq data from human patients. Our findings revealed coordinated changes in enhancers and transcriptional networks. We found that infected cells lose the WT1 transcription factor and demonstrate disruption of WT1-bound enhancers and of their associated target genes. Downstream targets of WT1 are involved in the regulation of the Wnt signaling and the mitogen-activated protein kinase cascade, which indeed exhibit increased activation levels. These findings may provide a potential explanation for the development of pulmonary fibrosis, a lethal complication of COVID-19. Moreover, we revealed over-acetylated enhancers associated with upregulated genes involved in cell adhesion, which could contribute to cell-cell infection of SARS-CoV-2. Furthermore, we demonstrated that enhancers may play a role in the activation of pro-inflammatory cytokines and contribute to excessive inflammation in the lungs, a typical complication of COVID-19. Overall, our analysis provided novel insights into the cell-autonomous dysregulation of enhancer regulation caused by SARS-CoV-2 infection, a step on the path to a deeper molecular understanding of the disease.

Osteosarcoma is an aggressive bone tumor that primarily affects children and adolescents. This malignancy is highly aggressive, associated with poor clinical outcomes, and primarily metastasizes to the lungs. Due to its rarity and biological heterogeneity, limited studies on its molecular basis exist, hindering the development of effective therapies. The WW domain-containing oxidoreductase (WWOX) is frequently altered in human osteosarcoma. Combined deletion of Wwox and Trp53 using Osterix1-Cre transgenic mice has been shown to accelerate osteosarcoma development. In this study, we generated a traceable osteosarcoma mouse model harboring the deletion of Trp53 alone (single-knockout) or combined deletion of Wwox/Trp53 (double-knockout) and expressing a tdTomato reporter. By tracking Tomato expression at different time points, we detected the early presence of tdTomato-positive cells in the bone marrow mesenchymal stem cells of non-osteosarcoma-bearing mice (young BM). We found that double-knockout young BM cells, but not single-knockout young BM cells, exhibited tumorigenic traits both in vitro and in vivo. Molecular and cellular characterization of these double-knockout young BM cells revealed their resemblance to osteosarcoma tumor cells. Interestingly, one of the observed significant transcriptomic changes in double-knockout young BM cells was the upregulation of Myc and its target genes compared to single-knockout young BM cells. Intriguingly, Myc-chromatin immunoprecipitation sequencing revealed its increased enrichment on Myc targets, which were upregulated in double-knockout young BM cells. Restoration of WWOX in double-knockout young BM cells reduced Myc protein levels. As a prototype target, we demonstrated the upregulation of MCM7, a known Myc target, in double-knockout young BM relative to single-knockout young BM cells. Inhibition of MCM7 expression using simvastatin resulted in reduced proliferation and tumor cell growth of double-knockout young BM cells. Our findings reveal BM mesenchymal stem cells as a platform to study osteosarcoma and Myc and its targets as WWOX effectors and early molecular events during osteosarcomagenesis.

Most pancreatic neuroendocrine tumors (PNETs) do not produce excess hormones and are therefore considered 'non-functional'1-3. As clinical behaviors vary widely and distant metastases are eventually lethal2,4, biological classifications might guide treatment. Using enhancer maps to infer gene regulatory programs, we find that non-functional PNETs fall into two major subtypes, with epigenomes and transcriptomes that partially resemble islet α- and β-cells. Transcription factors ARX and PDX1 specify these normal cells, respectively5,6, and 84% of 142 non-functional PNETs expressed one or the other factor, occasionally both. Among 103 cases, distant relapses occurred almost exclusively in patients with ARX+PDX1- tumors and, within this subtype, in cases with alternative lengthening of telomeres. These markedly different outcomes belied similar clinical presentations and histology and, in one cohort, occurred irrespective of MEN1 mutation. This robust molecular stratification provides insight into cell lineage correlates of non-functional PNETs, accurately predicts disease course and can inform postoperative clinical decisions.

Prostate cancer is the second most common cause of male cancer deaths in the United States. However, the full range of prostate cancer genomic alterations is incompletely characterized. Here we present the complete sequence of seven primary human prostate cancers and their paired normal counterparts. Several tumours contained complex chains of balanced (that is, 'copy-neutral') rearrangements that occurred within or adjacent to known cancer genes. Rearrangement breakpoints were enriched near open chromatin, androgen receptor and ERG DNA binding sites in the setting of the ETS gene fusion TMPRSS2-ERG, but inversely correlated with these regions in tumours lacking ETS fusions. This observation suggests a link between chromatin or transcriptional regulation and the genesis of genomic aberrations. Three tumours contained rearrangements that disrupted CADM2, and four harboured events disrupting either PTEN (unbalanced events), a prostate tumour suppressor, or MAGI2 (balanced events), a PTEN interacting protein not previously implicated in prostate tumorigenesis. Thus, genomic rearrangements may arise from transcriptional or chromatin aberrancies and engage prostate tumorigenic mechanisms.

Eliminating the bacterial cloning step has been a major factor in the vastly improved efficiency of massively parallel sequencing approaches. However, this also has made it a technical challenge to produce the modern equivalent of the Fosmid- or BAC-end sequences that were crucial for assembling and analyzing complex genomes during the Sanger-based sequencing era. To close this technology gap, we developed Fosill, a method for converting Fosmids to Illumina-compatible jumping libraries. We constructed Fosmid libraries in vectors with Illumina primer sequences and specific nicking sites flanking the cloning site. Our family of pFosill vectors allows multiplex Fosmid cloning of end-tagged genomic fragments without physical size selection and is compatible with standard and multiplex paired-end Illumina sequencing. To excise the bulk of each cloned insert, we introduced two nicks in the vector, translated them into the inserts, and cleaved them. Recircularization of the vector via coligation of insert termini followed by inverse PCR generates a jumping library for paired-end sequencing with 101-base reads. The yield of unique Fosmid-sized jumps is sufficiently high, and the background of short, incorrectly spaced and chimeric artifacts sufficiently low, to enable applications such as mapping of structural variation and scaffolding of de novo assemblies. We demonstrate the power of Fosill to map genome rearrangements in a cancer cell line and identified three fusion genes that were corroborated by RNA-seq data. Our Fosill-powered assembly of the mouse genome has an N50 scaffold length of 17.0 Mb, rivaling the connectivity (16.9 Mb) of the Sanger-sequencing based draft assembly.

Translocation events are frequent in cancer and may create chimeric fusions or 'regulatory rearrangements' that drive oncogene overexpression. Here we identify super-enhancer translocations that drive overexpression of the oncogenic transcription factor MYB as a recurrent theme in adenoid cystic carcinoma (ACC). Whole-genome sequencing data and chromatin maps highlight distinct chromosomal rearrangements that juxtapose super-enhancers to the MYB locus. Chromosome conformation capture confirms that the translocated enhancers interact with the MYB promoter. Remarkably, MYB protein binds to the translocated enhancers, creating a positive feedback loop that sustains its expression. MYB also binds enhancers that drive different regulatory programs in alternate cell lineages in ACC, cooperating with TP63 in myoepithelial cells and a Notch program in luminal epithelial cells. Bromodomain inhibitors slow tumor growth in ACC primagraft models in vivo. Thus, our study identifies super-enhancer translocations that drive MYB expression and provides insight into downstream MYB functions in alternate ACC lineages.

Multiple myeloma is an incurable malignancy of plasma cells, and its pathogenesis is poorly understood. Here we report the massively parallel sequencing of 38 tumour genomes and their comparison to matched normal DNAs. Several new and unexpected oncogenic mechanisms were suggested by the pattern of somatic mutation across the data set. These include the mutation of genes involved in protein translation (seen in nearly half of the patients), genes involved in histone methylation, and genes involved in blood coagulation. In addition, a broader than anticipated role of NF-κB signalling was indicated by mutations in 11 members of the NF-κB pathway. Of potential immediate clinical relevance, activating mutations of the kinase BRAF were observed in 4% of patients, suggesting the evaluation of BRAF inhibitors in multiple myeloma clinical trials. These results indicate that cancer genome sequencing of large collections of samples will yield new insights into cancer not anticipated by existing knowledge.

Breast cancer is the second leading cause of death in women above 60 years in the US. Screening mammography is recommended for women above 50 years; however, 22% of breast cancer cases are diagnosed in women below this age. We set out to develop a test based on the detection of cell-free RNA from saliva. To this end, we sequenced RNA from a pool of ten women. The 1254 transcripts identified were enriched for genes with an annotation of alternative pre-mRNA splicing. Pre-mRNA splicing is a tightly regulated process and its misregulation in cancer cells promotes the formation of cancer-driving isoforms. For these reasons, we chose to focus on splicing factors as biomarkers for the early detection of breast cancer. We found that the level of the splicing factors is unique to each woman and consistent in the same woman at different time points. Next, we extracted RNA from 36 healthy subjects and 31 breast cancer patients. Recording the mRNA level of seven splicing factors in these samples demonstrated that the combination of all these factors is different in the two groups (p value = 0.005). Our results demonstrate a differential abundance of splicing factor mRNA in the saliva of breast cancer patients.

Epigenetic aberrations are widespread in cancer, yet the underlying mechanisms and causality remain poorly understood1-3. A subset of gastrointestinal stromal tumours (GISTs) lack canonical kinase mutations but instead have succinate dehydrogenase (SDH) deficiency and global DNA hyper-methylation4,5. Here, we associate this hyper-methylation with changes in genome topology that activate oncogenic programs. To investigate epigenetic alterations systematically, we mapped DNA methylation, CTCF insulators, enhancers, and chromosome topology in KIT-mutant, PDGFRA-mutant and SDH-deficient GISTs. Although these respective subtypes shared similar enhancer landscapes, we identified hundreds of putative insulators where DNA methylation replaced CTCF binding in SDH-deficient GISTs. We focused on a disrupted insulator that normally partitions a core GIST super-enhancer from the FGF4 oncogene. Recurrent loss of this insulator alters locus topology in SDH-deficient GISTs, allowing aberrant physical interaction between enhancer and oncogene. CRISPR-mediated excision of the corresponding CTCF motifs in an SDH-intact GIST model disrupted the boundary between enhancer and oncogene, and strongly upregulated FGF4 expression. We also identified a second recurrent insulator loss event near the KIT oncogene, which is also highly expressed across SDH-deficient GISTs. Finally, we established a patient-derived xenograft (PDX) from an SDH-deficient GIST that faithfully maintains the epigenetics of the parental tumour, including hypermethylation and insulator defects. This PDX model is highly sensitive to FGF receptor (FGFR) inhibition, and more so to combined FGFR and KIT inhibition, validating the functional significance of the underlying epigenetic lesions. Our study reveals how epigenetic alterations can drive oncogenic programs in the absence of canonical kinase mutations, with implications for mechanistic targeting of aberrant pathways in cancers.

Genome-wide profiling of histone modifications can provide systematic insight into the regulatory elements and programs engaged in a given cell type. However, conventional chromatin immunoprecipitation and sequencing (ChIP-seq) does not capture quantitative information on histone modification levels, requires large amounts of starting material, and involves tedious processing of each individual sample. Here, we address these limitations with a technology that leverages DNA barcoding to profile chromatin quantitatively and in multiplexed format. We concurrently map relative levels of multiple histone modifications across multiple samples, each comprising as few as a thousand cells. We demonstrate the technology by monitoring dynamic changes following inhibition of p300, EZH2, or KDM5, by linking altered epigenetic landscapes to chromatin regulator mutations, and by mapping active and repressive marks in purified human hematopoietic stem cells. Hence, this technology enables quantitative studies of chromatin state dynamics across rare cell types, genotypes, environmental conditions, and drug treatments.

Mutations in the serine/threonine kinase BRAF are found in more than 60% of melanomas. The most prevalent melanoma mutation is BRAF(V600E), which constitutively activates downstream MAPK signalling. Vemurafenib is a potent RAF kinase inhibitor with remarkable clinical activity in BRAF(V600E)-positive melanoma tumours. However, patients rapidly develop resistance to vemurafenib treatment. One resistance mechanism is the emergence of BRAF alternative splicing isoforms leading to elimination of the RAS-binding domain. Here we identify interference with pre-mRNA splicing as a mechanism to combat vemurafenib resistance. We find that small-molecule pre-mRNA splicing modulators reduce BRAF3-9 production and limit in-vitro cell growth of vemurafenib-resistant cells. In xenograft models, interference with pre-mRNA splicing prevents tumour formation and slows growth of vemurafenib-resistant tumours. Our results identify an intronic mutation as the molecular basis for a RNA splicing-mediated RAF inhibitor resistance mechanism and we identify pre-mRNA splicing interference as a potential therapeutic strategy for drug resistance in BRAF melanoma.

Matrin 3 (MATR3) is a highly conserved, inner nuclear matrix protein with two zinc finger domains and two RNA recognition motifs (RRM), whose function is largely unknown. Recently we found MATR3 to be phosphorylated by the protein kinase ATM, which activates the cellular response to double strand breaks in the DNA. Here, we show that MATR3 interacts in an RNA-dependent manner with several proteins with established roles in RNA processing, and maintains its interaction with RNA via its RRM2 domain. Deep sequencing of the bound RNA (RIP-seq) identified several small noncoding RNA species. Using microarray analysis to explore MATR3's role in transcription, we identified 77 transcripts whose amounts depended on the presence of MATR3. We validated this finding with nine transcripts which were also bound to the MATR3 complex. Finally, we demonstrated the importance of MATR3 for maintaining the stability of several of these mRNA species and conclude that it has a role in mRNA stabilization. The data suggest that the cellular level of MATR3, known to be highly regulated, modulates the stability of a group of gene transcripts.

Whole-genome sequencing using massively parallel sequencing technologies enables accurate detection of somatic rearrangements in cancer. Pinpointing large numbers of rearrangement breakpoints to base-pair resolution allows analysis of rearrangement microhomology and genomic location for every sample. Here we analyze 95 tumor genome sequences from breast, head and neck, colorectal, and prostate carcinomas, and from melanoma, multiple myeloma, and chronic lymphocytic leukemia. We discover three genomic factors that are significantly correlated with the distribution of rearrangements: replication time, transcription rate, and GC content. The correlation is complex, and different patterns are observed between tumor types, within tumor types, and even between different types of rearrangements. Mutations in the APC gene correlate with and, hence, potentially contribute to DNA breakage in late-replicating, low %GC, untranscribed regions of the genome. We show that somatic rearrangements display less microhomology than germline rearrangements, and that breakpoint loci are correlated with local hypermutability with a particular enrichment for transversions.

Jumonji-domain-containing protein 6 (JMJD6) is a Fe(II) and 2-oxogluterate (2OG) dependent oxygenase involved in gene regulation through post-translationally modifying nuclear proteins. It is highly expressed in many cancer types and linked to tumor progression and metastasis. Four alternatively-spliced jmjd6 transcripts were annotated. Here, we focus on the two most abundantly expressed ones, which we call jmjd6-2 and jmjd6-Ex5. TCGA SpliceSeq data revealed a significant decrease of jmjd6-Ex5 transcripts in patients and postmortem tissue of several tumors. The two protein isoforms are distinguished by their C-terminal sequences, which include a serine-rich region (polyS-domain) in JMJD6-2 that is not present in JMJD6-Ex5. Immunoprecipitation followed by LC-MS/MS for JMJD6-Ex5 shows that different sets of proteins interact with JMJD6-2 and JMJD6-Ex5 with only a few overlaps. In particular, we found TFIIF-associating CTD phosphatase (FCP1), proteins of the survival of motor neurons (SMN) complex, heterogeneous nuclear ribonucleoproteins (hnRNPs) and upstream binding factor (UBF) to interact with JMJD6-Ex5. Like JMJD6-2, both UBF and FCP1 comprise a polyS-domain. The polyS domain of JMJD6-2 might block the interaction with polyS-domains of other proteins. In contrast, JMJD6-2 interacts with many SR-like proteins with arginine/serine-rich (RS)-domains, including several splicing factors. In an HIV-based splicing reporter assay, co-expression of JMJD6-2 inhibited exon inclusion, whereas JMJD6-Ex5 did not have any effect. Furthermore, the silencing of jmjd6 by siRNAs favored jmjd6-Ex5 transcripts, suggesting that JMJD6 controls splicing of its own pre-mRNA. The distinct molecular properties of JMJD6-2 and JMJD6-Ex5 open a lead into the functional implications of the variations of their relative abundance in tumors.