Accurate estimation of pathway activity in single cells for clustering and differential analysis.

Inferring which and how biological pathways and gene sets change is a key question in many studies that utilize single-cell RNA sequencing. Typically, these questions are addressed by quantifying the enrichment of known gene sets in lists of genes derived from global analysis. Here we offer SiPSiC, a new method to infer pathway activity in every single cell.

DNMT3B splicing dysregulation mediated by SMCHD1 loss contributes to DUX4 overexpression and FSHD pathogenesis.

Structural maintenance of chromosomes flexible hinge domain-containing 1 (SMCHD1) is a noncanonical SMC protein and an epigenetic regulator. Mutations in SMCHD1 cause facioscapulohumeral muscular dystrophy (FSHD), by overexpressing DUX4 in muscle cells. Here, we demonstrate that SMCHD1 is a key regulator of alternative splicing in various cell types.

Senescence of human pancreatic beta cells enhances functional maturation through chromatin reorganization and promotes interferon responsiveness.

Senescent cells can influence the function of tissues in which they reside, and their propensity for disease. A portion of adult human pancreatic beta cells express the senescence marker p16, yet it is unclear whether they are in a senescent state, and how this affects insulin secretion. We analyzed single-cell transcriptome datasets of adult human beta cells, and found that p16-positive cells express senescence gene signatures, as well as elevated levels of beta-cell maturation genes, consistent with enhanced functionality.

Differentiation shifts from a reversible to an irreversible heterochromatin state at the DM1 locus.

Epigenetic defects caused by hereditary or de novo mutations are implicated in various human diseases. It remains uncertain whether correcting the underlying mutation can reverse these defects in patient cells. Here we show by the analysis of myotonic dystrophy type 1 (DM1)-related locus that in mutant human embryonic stem cells (hESCs), DNA methylation and H3K9me3 enrichments are completely abolished by repeat excision (CTG2000 expansion), whereas in patient myoblasts (CTG2600 expansion), repeat deletion fails to do so.

SARS-CoV-2 infection perturbs enhancer mediated transcriptional regulation of key pathways.

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.

Isoforms of the TAL1 transcription factor have different roles in hematopoiesis and cell growth.

T-cell acute lymphoblastic leukemia (T-ALL) protein 1 (TAL1) is a central transcription factor in hematopoiesis. The timing and level of TAL1 expression orchestrate the differentiation to specialized blood cells and its overexpression is a common cause of T-ALL. Here, we studied the 2 protein isoforms of TAL1, short and long, which are generated by the use of alternative promoters as well as by alternative splicing.

Single-cell RNA sequencing identifies a paracrine interaction that may drive oncogenic notch signaling in human adenoid cystic carcinoma.

Salivary adenoid cystic carcinoma (ACC) is a rare, biologically unique biphasic tumor that consists of malignant myoepithelial and luminal cells. MYB and Notch signaling have been implicated in ACC pathophysiology, but in vivo descriptions of these two programs in human tumors and investigation into their active coordination remain incomplete. We utilize single-cell RNA sequencing to profile human head and neck ACC, including a comparison of primary ACC with a matched local recurrence.

Simultaneous Mapping of Enhancers and Enhancer Rearrangements with Paired-End H3K27ac ChIP-seq.

We describe a protocol for H3K27ac ChIP paired-end sequencing and computational analysis of rearrangements. Our approach can be used to simultaneously map enhancers and their activity and to identify structural variations at enhancers. Since changes in enhancer activity and new enhancer translocations both play a major role in tumor initiation, progression, and response to therapy, this approach holds promise to uncover some of the mechanisms behind these processes.

VEGFA's distal enhancer regulates its alternative splicing in CML.

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.

Single-cell RNA-seq reveals developmental plasticity with coexisting oncogenic states and immune evasion programs in ETP-ALL.

Lineage plasticity and stemness have been invoked as causes of therapy resistance in cancer, because these flexible states allow cancer cells to dedifferentiate and alter their dependencies. We investigated such resistance mechanisms in relapsed/refractory early T-cell progenitor acute lymphoblastic leukemia (ETP-ALL) carrying activating NOTCH1 mutations via full-length single-cell RNA sequencing (scRNA-seq) of malignant and microenvironmental cells. We identified 2 highly distinct stem-like states that critically differed with regard to cell cycle and oncogenic signaling.

Large-Scale Topological Changes Restrain Malignant Progression in Colorectal Cancer.

Widespread changes to DNA methylation and chromatin are well documented in cancer, but the fate of higher-order chromosomal structure remains obscure. Here we integrated topological maps for colon tumors and normal colons with epigenetic, transcriptional, and imaging data to characterize alterations to chromatin loops, topologically associated domains, and large-scale compartments. We found that spatial partitioning of the open and closed genome compartments is profoundly compromised in tumors.

Altered chromosomal topology drives oncogenic programs in SDH-deficient GISTs.

Epigenetic aberrations are widespread in cancer, yet the underlying mechanisms and causality remain poorly understood. A subset of gastrointestinal stromal tumours (GISTs) lack canonical kinase mutations but instead have succinate dehydrogenase (SDH) deficiency and global DNA hyper-methylation. Here, we associate this hyper-methylation with changes in genome topology that activate oncogenic programs.

Publisher Correction: Enhancer signatures stratify and predict outcomes of non-functional pancreatic neuroendocrine tumors.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Targeting FGFR overcomes EMT-mediated resistance in EGFR mutant non-small cell lung cancer.

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.

Enhancer signatures stratify and predict outcomes of non-functional pancreatic neuroendocrine tumors.

Most pancreatic neuroendocrine tumors (PNETs) do not produce excess hormones and are therefore considered 'non-functional'. As clinical behaviors vary widely and distant metastases are eventually lethal, 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.

Epithelial-to-Mesenchymal Transition Antagonizes Response to Targeted Therapies in Lung Cancer by Suppressing BIM.

Epithelial-to-mesenchymal transition (EMT) confers resistance to a number of targeted therapies and chemotherapies. However, it has been unclear why EMT promotes resistance, thereby impairing progress to overcome it. We have developed several models of EMT-mediated resistance to EGFR inhibitors (EGFRi) in -mutant lung cancers to evaluate a novel mechanism of EMT-mediated resistance.

Regulation of the glucocorticoid receptor via a BET-dependent enhancer drives antiandrogen resistance in prostate cancer.

In prostate cancer, resistance to the antiandrogen enzalutamide (Enz) can occur through bypass of androgen receptor (AR) blockade by the glucocorticoid receptor (GR). In contrast to fixed genomic alterations, here we show that GR-mediated antiandrogen resistance is adaptive and reversible due to regulation of GR expression by a tissue-specific enhancer. GR expression is silenced in prostate cancer by a combination of AR binding and EZH2-mediated repression at the GR locus, but is restored in advanced prostate cancers upon reversion of both repressive signals.

SMARCB1-mediated SWI/SNF complex function is essential for enhancer regulation.

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.

Genetic interrogation of circulating multiple myeloma cells at single-cell resolution.

Multiple myeloma (MM) remains an incurable disease, with a treatment-refractory state eventually developing in all patients. Constant clonal evolution and genetic heterogeneity of MM are a likely explanation for the emergence of drug-resistant disease. Monitoring of MM genomic evolution on therapy by serial bone marrow biopsy is unfortunately impractical because it involves an invasive and painful procedure.

An oncogenic MYB feedback loop drives alternate cell fates in adenoid cystic carcinoma.

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.

Insulator dysfunction and oncogene activation in IDH mutant gliomas.

Gain-of-function IDH mutations are initiating events that define major clinical and prognostic classes of gliomas. Mutant IDH protein produces a new onco-metabolite, 2-hydroxyglutarate, which interferes with iron-dependent hydroxylases, including the TET family of 5'-methylcytosine hydroxylases. TET enzymes catalyse a key step in the removal of DNA methylation.