Triple-negative breast cancer (TNBC) is an aggressive subtype characterized by extensive intratumoral heterogeneity. To investigate the underlying biology, we conducted single-cell RNA-sequencing (scRNA-seq) of >1500 cells from six primary TNBC. Here, we show that intercellular heterogeneity of gene expression programs within each tumor is variable and largely correlates with clonality of inferred genomic copy number changes, suggesting that genotype drives the gene expression phenotype of individual subpopulations. Clustering of gene expression profiles identified distinct subgroups of malignant cells shared by multiple tumors, including a single subpopulation associated with multiple signatures of treatment resistance and metastasis, and characterized functionally by activation of glycosphingolipid metabolism and associated innate immunity pathways. A novel signature defining this subpopulation predicts long-term outcomes for TNBC patients in a large cohort. Collectively, this analysis reveals the functional heterogeneity and its association with genomic evolution in TNBC, and uncovers unanticipated biological principles dictating poor outcomes in this disease.

Mutational inactivation of the SWI/SNF chromatin regulator ATRX occurs frequently in gliomas, the most common primary brain tumors. Whether and how ATRX deficiency promotes oncogenesis by epigenomic dysregulation remains unclear, despite its recent implication in both genomic instability and telomere dysfunction. Here we report that Atrx loss recapitulates characteristic disease phenotypes and molecular features in putative glioma cells of origin, inducing cellular motility although also shifting differentiation state and potential toward an astrocytic rather than neuronal histiogenic profile. Moreover, Atrx deficiency drives widespread shifts in chromatin accessibility, histone composition, and transcription in a distribution almost entirely restricted to genomic sites normally bound by the protein. Finally, direct gene targets of Atrx that mediate specific Atrx-deficient phenotypes in vitro exhibit similarly selective misexpression in ATRX-mutant human gliomas. These findings demonstrate that ATRX deficiency and its epigenomic sequelae are sufficient to induce disease-defining oncogenic phenotypes in appropriate cellular and molecular contexts.

Sarcomas are malignant soft tissue and bone tumours affecting adults, adolescents and children. They represent a morphologically heterogeneous class of tumours and some entities lack defining histopathological features. Therefore, the diagnosis of sarcomas is burdened with a high inter-observer variability and misclassification rate. Here, we demonstrate classification of soft tissue and bone tumours using a machine learning classifier algorithm based on array-generated DNA methylation data. This sarcoma classifier is trained using a dataset of 1077 methylation profiles from comprehensively pre-characterized cases comprising 62 tumour methylation classes constituting a broad range of soft tissue and bone sarcoma subtypes across the entire age spectrum. The performance is validated in a cohort of 428 sarcomatous tumours, of which 322 cases were classified by the sarcoma classifier. Our results demonstrate the potential of the DNA methylation-based sarcoma classification for research and future diagnostic applications.

Plakophilin-2 (PKP2) is a component of the desmosome and known for its role in cell-cell adhesion. Mutations in human PKP2 associate with a life-threatening arrhythmogenic cardiomyopathy, often of right ventricular predominance. Here, we use a range of state-of-the-art methods and a cardiomyocyte-specific, tamoxifen-activated, PKP2 knockout mouse to demonstrate that in addition to its role in cell adhesion, PKP2 is necessary to maintain transcription of genes that control intracellular calcium cycling. Lack of PKP2 reduces expression of Ryr2 (coding for Ryanodine Receptor 2), Ank2 (coding for Ankyrin-B), Cacna1c (coding for CaV1.2) and Trdn (coding for triadin), and protein levels of calsequestrin-2 (Casq2). These factors combined lead to disruption of intracellular calcium homeostasis and isoproterenol-induced arrhythmias that are prevented by flecainide treatment. We propose a previously unrecognized arrhythmogenic mechanism related to PKP2 expression and suggest that mutations in PKP2 in humans may cause life-threatening arrhythmias even in the absence of structural disease.It is believed that mutations in desmosomal adhesion complex protein plakophilin 2 (PKP2) cause arrhythmia due to loss of cell-cell communication. Here the authors show that PKP2 controls the expression of proteins involved in calcium cycling in adult mouse hearts, and that lack of PKP2 can cause arrhythmia in a structurally normal heart.

In vivo functional investigation of oncogenes using somatic gene transfer has been successfully exploited to validate their role in tumorigenesis. For tumour suppressor genes this has proven more challenging due to technical aspects. To provide a flexible and effective method for investigating somatic loss-of-function alterations and their influence on tumorigenesis, we have established CRISPR/Cas9-mediated somatic gene disruption, allowing for in vivo targeting of TSGs. Here we demonstrate the utility of this approach by deleting single (Ptch1) or multiple genes (Trp53, Pten, Nf1) in the mouse brain, resulting in the development of medulloblastoma and glioblastoma, respectively. Using whole-genome sequencing (WGS) we characterized the medulloblastoma-driving Ptch1 deletions in detail and show that no off-targets were detected in these tumours. This method provides a fast and convenient system for validating the emerging wealth of novel candidate tumour suppressor genes and the generation of faithful animal models of human cancer.

Chromatin organization is a highly orchestrated process that influences gene expression, in part by modulating access of regulatory factors to DNA and nucleosomes. Here, we report that the chromatin accessibility regulator HMGN1, a target of recurrent DNA copy gains in leukemia, controls myeloid differentiation. HMGN1 amplification is associated with increased accessibility, expression, and histone H3K27 acetylation of loci important for hematopoietic stem cells (HSCs) and leukemia, such as HoxA cluster genes. In vivo, HMGN1 overexpression is linked to decreased quiescence and increased HSC activity in bone marrow transplantation. HMGN1 overexpression also cooperates with the AML-ETO9a fusion oncoprotein to impair myeloid differentiation and enhance leukemia stem cell (LSC) activity. Inhibition of histone acetyltransferases CBP/p300 relieves the HMGN1-associated differentiation block. These data nominate factors that modulate chromatin accessibility as regulators of HSCs and LSCs, and suggest that targeting HMGN1 or its downstream effects on histone acetylation could be therapeutically active in AML.

Pineoblastoma is a rare and highly aggressive brain cancer of childhood, histologically belonging to the spectrum of primitive neuroectodermal tumors. Patients with germline mutations in DICER1, a ribonuclease involved in microRNA processing, have increased risk of pineoblastoma, but genetic drivers of sporadic pineoblastoma remain unknown. Here, we analyzed pediatric and adult pineoblastoma samples (n = 23) using a combination of genome-wide DNA methylation profiling and whole-exome sequencing or whole-genome sequencing. Pediatric and adult pineoblastomas showed distinct methylation profiles, the latter clustering with lower-grade pineal tumors and normal pineal gland. Recurrent variants were found in genes involved in PKA- and NF-κB signaling, as well as in chromatin remodeling genes. We identified recurrent homozygous deletions of DROSHA, acting upstream of DICER1 in microRNA processing, and a novel microduplication involving chromosomal region 1q21 containing PDE4DIP (myomegalin), comprising the ancient DUF1220 protein domain. Expresion of PDE4DIP and DUF1220 proteins was present exclusively in pineoblastoma with PDE4DIP gain.

As the number of genomics datasets grows rapidly, sample mislabeling has become a high stakes issue. We present CrosscheckFingerprints (Crosscheck), a tool for quantifying sample-relatedness and detecting incorrectly paired sequencing datasets from different donors. Crosscheck outperforms similar methods and is effective even when data are sparse or from different assays. Application of Crosscheck to 8851 ENCODE ChIP-, RNA-, and DNase-seq datasets enabled us to identify and correct dozens of mislabeled samples and ambiguous metadata annotations, representing ~1% of ENCODE datasets.

The tumor microenvironment (TME) in pancreatic ductal adenocarcinoma (PDAC) is a complex ecosystem that drives tumor progression; however, in-depth single cell characterization of the PDAC TME and its role in response to therapy is lacking. Here, we perform single-cell RNA sequencing on freshly collected human PDAC samples either before or after chemotherapy. Overall, we find a heterogeneous mixture of basal and classical cancer cell subtypes, along with distinct cancer-associated fibroblast and macrophage subpopulations. Strikingly, classical and basal-like cancer cells exhibit similar transcriptional responses to chemotherapy and do not demonstrate a shift towards a basal-like transcriptional program among treated samples. We observe decreased ligand-receptor interactions in treated samples, particularly between TIGIT on CD8 + T cells and its receptor on cancer cells, and identify TIGIT as the major inhibitory checkpoint molecule of CD8 + T cells. Our results suggest that chemotherapy profoundly impacts the PDAC TME and may promote resistance to immunotherapy.

Approximately 30% of early-stage lung adenocarcinoma patients present with disease progression after successful surgical resection. Despite efforts of mapping the genetic landscape, there has been limited success in discovering predictive biomarkers of disease outcomes. Here we performed a systematic multi-omic assessment of 143 tumors and matched tumor-adjacent, histologically-normal lung tissue with long-term patient follow-up. Through histologic, mutational, and transcriptomic profiling of tumor and adjacent-normal tissue, we identified an inflammatory gene signature in tumor-adjacent tissue as the strongest clinical predictor of disease progression. Single-cell transcriptomic analysis demonstrated the progression-associated inflammatory signature was expressed in both immune and non-immune cells, and cell type-specific profiling in monocytes further improved outcome predictions. Additional analyses of tumor-adjacent transcriptomic data from The Cancer Genome Atlas validated the association of the inflammatory signature with worse outcomes across cancers. Collectively, our study suggests that molecular profiling of tumor-adjacent tissue can identify patients at high risk for disease progression.

Cancer evolution is fueled by epigenetic as well as genetic diversity. In chronic lymphocytic leukemia (CLL), intra-tumoral DNA methylation (DNAme) heterogeneity empowers evolution. Here, to comprehensively study the epigenetic dimension of cancer evolution, we integrate DNAme analysis with histone modification mapping and single cell analyses of RNA expression and DNAme in 22 primary CLL and 13 healthy donor B lymphocyte samples. Our data reveal corrupted coherence across different layers of the CLL epigenome. This manifests in decreased mutual information across epigenetic modifications and gene expression attributed to cell-to-cell heterogeneity. Disrupted epigenetic-transcriptional coordination in CLL is also reflected in the dysregulation of the transcriptional output as a function of the combinatorial chromatin states, including incomplete Polycomb-mediated gene silencing. Notably, we observe unexpected co-mapping of typically mutually exclusive activating and repressing histone modifications, suggestive of intra-tumoral epigenetic diversity. Thus, CLL epigenetic diversification leads to decreased coordination across layers of epigenetic information, likely reflecting an admixture of cells with diverging cellular identities.

The human genome is folded into regulatory units termed 'topologically-associated domains' (TADs). Genome-wide studies support a global role for the insulator protein CTCF in mediating chromosomal looping and the topological constraint of TAD boundaries. However, the impact of individual insulators on enhancer-gene interactions and transcription remains poorly understood. Here, we investigate epigenome editing strategies for perturbing individual CTCF insulators and evaluating consequent effects on genome topology and transcription. We show that fusions of catalytically-inactive Cas9 (dCas9) to transcriptional repressors (dCas9-KRAB) and DNA methyltransferases (dCas9-DNMT3A, dCas9-DNMT3A3L) can selectively displace CTCF from specific insulators, but only when precisely targeted to the cognate motif. We further demonstrate that stable, partially-heritable insulator disruption can be achieved through combinatorial hit-and-run epigenome editing. Finally, we apply these strategies to simulate an insulator loss mechanism implicated in brain tumorigenesis. Our study provides strategies for stably modifying genome organization and gene activity without altering the underlying DNA sequence.

Chromatin regulators are frequently mutated in human cancer and are attractive drug targets. They include diverse proteins that share functional domains and assemble into related multi-subunit complexes. To investigate functional relationships among these regulators, here we apply combinatorial CRISPR knockouts (KOs) to test over 35,000 gene-gene pairings in leukemia cells, using a library of over 300,000 constructs. Top pairs that demonstrate either compensatory non-lethal interactions or synergistic lethality enrich for paralogs and targets that occupy the same protein complex. The screen highlights protein complex dependencies not apparent in single KO screens, for example MCM histone exchange, the nucleosome remodeling and deacetylase (NuRD) complex, and HBO1 (KAT7) complex. We explore two approaches to NuRD complex inactivation. Paralog and non-paralog combinations of the KAT7 complex emerge as synergistic lethal and specifically nominate the ING5 PHD domain as a potential therapeutic target when paired with other KAT7 complex member losses. These findings highlight the power of combinatorial screening to provide mechanistic insight and identify therapeutic targets within redundant networks.