Detection of metabolic adaptation in a triple-negative breast cancer animal model with [18F]choline-PET imaging as a surrogate for drug resistance

PURPOSE

Test the feasibility of an image-based method to identify taxane resistance in mouse bearing triple-negative breast cancer (TNBC) tumor xenografts.

METHODS

Xenograft tumor-bearing mice from paclitaxel-sensitive and paclitaxel-resistant TNBC cells (MDA-MD-346) were generated by orthotopic injection into female NOD-SCID mice. When tumors reached 100-150 mm3, mice were scanned using [18F]choline PET/CT. Tumors were collected and sliced for autoradiography and immunofluorescence analysis. Quantitative data was analyzed accordingly.

RESULTS

From fifteen mice scanned, five had taxane-sensitive cell line tumors of which two underwent taxol-based treatment. From the remaining 10 mice with taxane-resistant cell line tumors, four underwent taxol-based treatment. Only 13 mice had the tumor sample analyzed histologically. When normalized to the blood pool, both cell lines showed differences in metabolic uptake before and after treatment.

CONCLUSIONS

Treated and untreated taxane-sensitive and taxane-resistant cell lines have different metabolic properties that could be leveraged before the start of chemotherapy.

Chromatin insulation orchestrates matrix metalloproteinase gene cluster expression reprogramming in aggressive breast cancer tumors

BACKGROUND

Triple-negative breast cancer (TNBC) is an aggressive subtype that exhibits a high incidence of distant metastases and lacks targeted therapeutic options. Here we explored how the epigenome contributes to matrix metalloprotease (MMP) dysregulation impacting tumor invasion, which is the first step of the metastatic process.

METHODS

We combined RNA expression and chromatin interaction data to identify insulator elements potentially associated with MMP gene expression and invasion. We employed CRISPR/Cas9 to disrupt the CCCTC-Binding Factor (CTCF) binding site on an insulator element downstream of the MMP8 gene (IE8) in two TNBC cellular models. We characterized these models by combining Hi-C, ATAC-seq, and RNA-seq with functional experiments to determine invasive ability. The potential of our findings to predict the progression of ductal carcinoma in situ (DCIS), was tested in data from clinical specimens.

RESULTS

We explored the clinical relevance of an insulator element located within the Chr11q22.2 locus, downstream of the MMP8 gene (IE8). This regulatory element resulted in a topologically associating domain (TAD) boundary that isolated nine MMP genes into two anti-correlated expression clusters. This expression pattern was associated with worse relapse-free (HR = 1.57 [1.06 - 2.33]; p = 0.023) and overall (HR = 2.65 [1.31 - 5.37], p = 0.005) survival of TNBC patients. After CRISPR/Cas9-mediated disruption of IE8, cancer cells showed a switch in the MMP expression signature, specifically downregulating the pro-invasive MMP1 gene and upregulating the antitumorigenic MMP8 gene, resulting in reduced invasive ability and collagen degradation. We observed that the MMP expression pattern predicts DCIS that eventually progresses into invasive ductal carcinomas (AUC = 0.77, p < 0.01).

CONCLUSION

Our study demonstrates how the activation of an IE near the MMP8 gene determines the regional transcriptional regulation of MMP genes with opposing functional activity, ultimately influencing the invasive properties of aggressive forms of breast cancer.

3-D chromatin conformation, accessibility, and gene expression profiling of triple-negative breast cancer

OBJECTIVES

Triple-negative breast cancer (TNBC) is a highly aggressive breast cancer subtype with limited treatment options. Unlike other breast cancer subtypes, the scarcity of specific therapies and greater frequencies of distant metastases contribute to its aggressiveness. We aimed to find epigenetic changes that aid in the understanding of the dissemination process of these cancers.

DATA DESCRIPTION

Using CRISPR/Cas9, our experimental approach led us to identify and disrupt an insulator element, IE8, whose activity seemed relevant for cell invasion. The experiments were performed in two well-established TNBC cellular models, the MDA-MB-231 and the MDA-MB-436. To gain insights into the underlying molecular mechanisms of TNBC invasion ability, we generated and characterized high-resolution chromatin interaction (Hi-C) and chromatin accessibility (ATAC-seq) maps in both cell models and complemented these datasets with gene expression profiling (RNA-seq) in MDA-MB-231, the cell line that showed more significant changes in chromatin accessibility. Altogether, our data provide a comprehensive resource for understanding the spatial organization of the genome in TNBC cells, which may contribute to accelerating the discovery of TNBC-specific alterations triggering advances for this devastating disease.

Transposable Elements Are Co-opted as Oncogenic Regulatory Elements by Lineage-Specific Transcription Factors in Prostate Cancer

Transposable elements hold regulatory functions that impact cell fate determination by controlling gene expression. However, little is known about the transcriptional machinery engaged at transposable elements in pluripotent and mature versus oncogenic cell states. Through positional analysis over repetitive DNA sequences of H3K27ac chromatin immunoprecipitation sequencing data from 32 normal cell states, we report pluripotent/stem and mature cell state-specific "regulatory transposable elements." Pluripotent/stem elements are binding sites for pluripotency factors (e.g., NANOG, SOX2, OCT4). Mature cell elements are docking sites for lineage-specific transcription factors, including AR and FOXA1 in prostate epithelium. Expanding the analysis to prostate tumors, we identify a subset of regulatory transposable elements shared with pluripotent/stem cells, including Tigger3a. Using chromatin editing technology, we show how such elements promote prostate cancer growth by regulating AR transcriptional activity. Collectively, our results suggest that oncogenesis arises from lineage-specific transcription factors hijacking pluripotent/stem cell regulatory transposable elements.

SIGNIFICANCE

We show that oncogenesis relies on co-opting transposable elements from pluripotent stem cells as regulatory elements altering the recruitment of lineage-specific transcription factors. We further discover how co-option is dependent on active chromatin states with important implications for developing treatment options against drivers of oncogenesis across the repetitive DNA. This article is featured in Selected Articles from This Issue, p. 2293.

The chromatin and single-cell transcriptional landscapes of CD4 T cells in inflammatory bowel disease link risk loci with a proinflammatory Th17 cell population

Introduction

The imbalance between Th17 and regulatory T cells in inflammatory bowel diseases (IBD) promotes intestinal epithelial cell damage. In this scenario, T helper cell lineage commitment is accompanied by dynamic changes to the chromatin that facilitate or repress gene expression.

Methods

Here, we characterized the chromatin landscape and heterogeneity of intestinal and peripheral CD4 T cellsfrom IBD patients using in house ATAC-Seq and single cell RNA-Seq libraries.

Results

We show that chromatin accessibility profiles of CD4 T cells from inflamed intestinal biopsies relate to genes associated with a network of inflammatory processes. After integrating the chromatin profiles of tissue-derived CD4 T cells and in-vitro polarized CD4 T cell subpopulations, we found that the chromatin accessibility changes of CD4 T cells were associated with a higher predominance of pathogenic Th17 cells (pTh17 cells) in inflamed biopsies. In addition, IBD risk loci in CD4 T cells were colocalized with accessible chromatin changes near pTh17-related genes, as shown in intronic STAT3 and IL23R regions enriched in areas of active intestinal inflammation. Moreover, single cell RNA-Seq analysis revealed a population of pTh17 cells that co-expresses Th1 and cytotoxic transcriptional programs associated with IBD severity.

Discussion

Altogether, we show that cytotoxic pTh17 cells were specifically associated with IBD genetic variants and linked to intestinal inflammation of IBD patients.

Human papillomavirus integration transforms chromatin to drive oncogenesis

BACKGROUND

Human papillomavirus (HPV) drives almost all cervical cancers and up to 70% of head and neck cancers. Frequent integration into the host genome occurs predominantly in tumorigenic types of HPV. We hypothesize that changes in chromatin state at the location of integration can result in changes in gene expression that contribute to the tumorigenicity of HPV.

RESULTS

We find that viral integration events often occur along with changes in chromatin state and expression of genes near the integration site. We investigate whether introduction of new transcription factor binding sites due to HPV integration could invoke these changes. Some regions within the HPV genome, particularly the position of a conserved CTCF binding site, show enriched chromatin accessibility signal. ChIP-seq reveals that the conserved CTCF binding site within the HPV genome binds CTCF in 4 HPV^+ cancer cell lines. Significant changes in CTCF binding pattern and increases in chromatin accessibility occur exclusively within 100 kbp of HPV integration sites. The chromatin changes co-occur with out-sized changes in transcription and alternative splicing of local genes. Analysis of The Cancer Genome Atlas (TCGA) HPV⁺ tumors indicates that HPV integration upregulates genes which have significantly higher essentiality scores compared to randomly selected upregulated genes from the same tumors.

CONCLUSIONS

Our results suggest that introduction of a new CTCF binding site due to HPV integration reorganizes chromatin state and upregulates genes essential for tumor viability in some HPV⁺ tumors. These findings emphasize a newly recognized role of HPV integration in oncogenesis.

Chronic hypoxia favours adoption to a castration-resistant cell state in prostate cancer

Predicting and treating recurrence in intermediate-risk prostate cancer patients remains a challenge despite having identified genomic instability [1] and hypoxia [2, 3] as risk factors. This underlies challenges in assigning the functional impact of these risk factors to mechanisms promoting prostate cancer progression. Here we show chronic hypoxia (CH), as observed in prostate tumours [4], leads to the adoption of an androgen-independent state in prostate cancer cells. Specifically, CH results in prostate cancer cells adopting transcriptional and metabolic alterations typical of castration-resistant prostate cancer cells. These changes include the increased expression of transmembrane transporters for the methionine cycle and related pathways leading to increased abundance of metabolites and expression of enzymes related to glycolysis. Targeting of the Glucose Transporter 1 (GLUT1) identified a dependency on glycolysis in androgen-independent cells. Overall, we identified a therapeutically targetable weakness in chronic hypoxia and androgen-independent prostate cancer. These findings may offer additional strategies for treatment development against hypoxic prostate cancer.

EPCO-38. TYPE B ULTRA LONG-RANGE INTERACTIONS IN PFAS (TULIPS) ARE RECURRENT EPIGENOMIC FEATURES OF PFA EPENDYMOMA

Posterior Fossa Group A (PFA) ependymomas are pediatric brain tumors with extremely poor survival outcomes. As protein-coding mutations in PFA are exceedingly rare, the underlying etiology of these tumors remains elusive. Elevated CpG island methylation and depletion of H3K27me3 have been described in PFA, leading to the hypothesis that PFA may be driven by a dysregulated epigenetic state. In this study, we sought to determine how three-dimensional (3D) genome features (such as DNA loops, domains, and compartments) differ between pediatric brain tumors. We performed Hi-C sequencing on a collection of 64 patient specimens and patient-derived primary cultures that collectively span multiple subgroups of ependymoma, medulloblastoma, high-grade glioma, and non-neoplastic brain. For certain samples, we further performed RNA-seq, histone modification ChIP-seq, or whole-genome bisulfite sequencing to allow multiomic data integration. Overall, the 3D genome organization of PFA samples appeared distinct from other tumor types. We identified and defined TULIPs: a subset of type B compartments, separated by genomic distances greater than 10 Mbp, that exhibit a striking fivefold increase in reciprocal interaction strength. These TULIPs recurred at the same genomic positions across the vast majority of PFA samples with minimal representation among other tumor or non-tumor samples. TULIPs displayed enrichment for heterochromatic features such as H3K9me3 and late replication timing and were depleted of euchromatic features such as H3K27ac and protein-coding genes. By using immuno-fluorescence for H3K9me3 and oligo-FISH to label TULIP regions, we demonstrated that TULIP regions are more compact in PFA than other tumors. Finally, by applying inhibitors of H3K9 lysine methylation to PFA cultures we showed that TULIPs become more diffuse and cell viability is reduced. Altogether, this work defines TULIPs as highly recurrent epigenetic features of PFA tumors.

A noncoding single-nucleotide polymorphism at 8q24 drives IDH1-mutant glioma formation

Establishing causal links between inherited polymorphisms and cancer risk is challenging. Here, we focus on the single-nucleotide polymorphism rs55705857, which confers a sixfold greater risk of isocitrate dehydrogenase (IDH)-mutant low-grade glioma (LGG). We reveal that rs55705857 itself is the causal variant and is associated with molecular pathways that drive LGG. Mechanistically, we show that rs55705857 resides within a brain-specific enhancer, where the risk allele disrupts OCT2/4 binding, allowing increased interaction with the Myc promoter and increased Myc expression. Mutating the orthologous mouse rs55705857 locus accelerated tumor development in an Idh1R132H-driven LGG mouse model from 472 to 172 days and increased penetrance from 30% to 75%. Our work reveals mechanisms of the heritable predisposition to lethal glioma in ~40% of LGG patients.

Chromatin Variants Reveal the Genetic Determinants of Oncogenesis in Breast Cancer

Breast cancer presents as multiple distinct disease entities. Each tumor harbors diverse cell populations defining a phenotypic heterogeneity that impinges on our ability to treat patients. To date, efforts mainly focused on genetic variants to find drivers of inter- and intratumor phenotypic heterogeneity. However, these efforts have failed to fully capture the genetic basis of breast cancer. Through recent technological and analytical approaches, the genetic basis of phenotypes can now be decoded by characterizing chromatin variants. These variants correspond to polymorphisms in chromatin states at DNA sequences that serve a distinct role across cell populations. Here, we review the function and causes of chromatin variants as they relate to breast cancer inter- and intratumor heterogeneity and how they can guide the development of treatment alternatives to fulfill the goal of precision cancer medicine.

Cancer-associated chromatin variants uncover the oncogenic role of transposable elements

The vast array of cell states found across human tissue arises from chromatin variants, which correspond to segments of the genome, known as DNA elements, adopting a different chromatin state over cell state transitions. Oncogenesis stems from alterations to the chromatin states over DNA elements that result in cancer-associated chromatin variants. Here, we review how cancer-associated chromatin variants call attention to repetitive DNA elements, and guide the functional characterization of transposable elements to decode their role in oncogenesis. We further discuss prevailing opportunities in the study of repetitive DNA elements to move towards the 'complete cancer genome' goal for precision medicine in oncology.

PRMT inhibition induces a viral mimicry response in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype with the worst prognosis and few effective therapies. Here we identified MS023, an inhibitor of type I protein arginine methyltransferases (PRMTs), which has antitumor growth activity in TNBC. Pathway analysis of TNBC cell lines indicates that the activation of interferon responses before and after MS023 treatment is a functional biomarker and determinant of response, and these observations extend to a panel of human-derived organoids. Inhibition of type I PRMT triggers an interferon response through the antiviral defense pathway with the induction of double-stranded RNA, which is derived, at least in part, from inverted repeat Alu elements. Together, our results represent a shift in understanding the antitumor mechanism of type I PRMT inhibitors and provide a rationale and biomarker approach for the clinical development of type I PRMT inhibitors.

Type I PRMT inhibition elicits potent antitumor activity associated with increased interferon response and intron-retained dsRNA accumulation, suggesting its potential combination with immune checkpoint inhibitors for cancer treatment.

PRMT5 regulates ATF4 transcript splicing and oxidative stress response

Protein methyltransferase 5 (PRMT5) symmetrically dimethylates arginine residues leading to regulation of transcription and splicing programs. Although PRMT5 has emerged as an attractive oncology target, the molecular determinants of PRMT5 dependency in cancer remain incompletely understood. Our transcriptomic analysis identified PRMT5 regulation of the activating transcription factor 4 (ATF4) pathway in acute myelogenous leukemia (AML). PRMT5 inhibition resulted in the expression of unstable, intron-retaining ATF4 mRNA that is detained in the nucleus. Concurrently, the decrease in the spliced cytoplasmic transcript of ATF4 led to lower levels of ATF4 protein and downregulation of ATF4 target genes. Upon loss of functional PRMT5, cells with low ATF4 displayed increased oxidative stress, growth arrest, and cellular senescence. Interestingly, leukemia cells with EVI1 oncogene overexpression demonstrated dependence on PRMT5 function. EVI1 and ATF4 regulated gene signatures were inversely correlated. We show that EVI1-high AML cells have reduced ATF4 levels, elevated baseline reactive oxygen species and increased sensitivity to PRMT5 inhibition. Thus, EVI1-high cells demonstrate dependence on PRMT5 function and regulation of oxidative stress response. Overall, our findings identify the PRMT5-ATF4 axis to be safeguarding the cellular redox balance that is especially important in high oxidative stress states, such as those that occur with EVI1 overexpression.

ecDNAs personify cancer gangsters

Alterations to gene regulatory plexuses typify oncogenesis, and two recent studies from Hung, Yost, Xie et al. (Hung et al., 2021) and Yi et al. (2021) collectively reveal the competitive advantage of extrachromosomal DNAs (ecDNAs) to set their own rules to control transcription by clustering into ecDNA hubs through intermolecular interactions that forge hub-specific regulatory plexuses driving high-oncogene expression.

Mutations in Noncoding Cis-Regulatory Elements Reveal Cancer Driver Cistromes in Luminal Breast Cancer

Whole-genome sequencing of primary breast tumors enabled the identification of cancer driver genes and noncoding cancer driver plexuses from somatic mutations. However, differentiating driver from passenger events among noncoding genetic variants remains a challenge. Herein, we reveal cancer-driver cis-regulatory elements linked to transcription factors previously shown to be involved in development of luminal breast cancers by defining a tumor-enriched catalogue of approximately 100,000 unique cis-regulatory elements from 26 primary luminal estrogen receptor (ER)⁺ progesterone receptor (PR)⁺ breast tumors. Integrating this catalog with somatic mutations from 350 publicly available breast tumor whole genomes, we uncovered cancer driver cistromes, defined as the sum of binding sites for a transcription factor, for ten transcription factors in luminal breast cancer such as FOXA1 and ER, nine of which are essential for growth in breast cancer with four exclusive to the luminal subtype. Collectively, we present a strategy to find cancer driver cistromes relying on quantifying the enrichment of noncoding mutations over cis-regulatory elements concatenated into a functional unit. IMPLICATIONS: Mapping the accessible chromatin of luminal breast cancer led to discovery of an accumulation of mutations within cistromes of transcription factors essential to luminal breast cancer. This demonstrates coopting of regulatory networks to drive cancer and provides a framework to derive insight into the noncoding space of cancer.

Alu insertion variants alter gene transcript levels

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

Copy-scAT: Deconvoluting single-cell chromatin accessibility of genetic subclones in cancer

Single-cell epigenomic assays have tremendous potential to illuminate mechanisms of transcriptional control in functionally diverse cancer cell populations. However, application of these techniques to clinical tumor specimens has been hampered by the current inability to distinguish malignant from nonmalignant cells, which potently confounds data analysis and interpretation. Here, we describe Copy-scAT, an R package that uses single-cell epigenomic data to infer copy number variants (CNVs) that define cancer cells. Copy-scAT enables studies of subclonal chromatin dynamics in complex tumors like glioblastoma. By deploying Copy-scAT, we uncovered potent influences of genetics on chromatin accessibility profiles in individual subclones. Consequently, some genetic subclones were predisposed to acquire stem-like or more differentiated molecular phenotypes, reminiscent of developmental paradigms. Copy-scAT is ideal for studies of the relationships between genetics and epigenetics in malignancies with high levels of intratumoral heterogeneity and to investigate how cancer cells interface with their microenvironment.

Reorganization of the 3D genome pinpoints non-coding drivers of primary prostate tumors

Prostate cancer is a heterogeneous disease whose progression is linked to genome instability. However, the impact of this instability on the non-coding genome and its three-dimensional organization to aid progression is unclear. Using primary benign and tumor tissue, we find a high concordance in higher order three-dimensional genome organization. This concordance argues for constraints to the topology of prostate tumor genomes. Nonetheless, we identified changes in focal chromatin interactions, typical of loops bridging non-coding cis-regulatory elements, and showed how structural variants can induce these changes to guide cis-regulatory element hijacking. Such events resulted in opposing differential expression of genes found at antipodes of rearrangements. Collectively, these results argue that changes to focal chromatin interactions, as opposed to higher order genome organization, allow for aberrant gene regulation and are repeatedly mediated by structural variants in primary prostate cancer.

Immune Resistance Interrogation Study (IRIS): A prospective comprehensive multi-omic analysis in patients with intrinsic and acquired resistance to immunotherapy

TPS2679

Background: Immune checkpoint inhibitors (ICI) have demonstrated efficacy in a wide variety of cancers. Nevertheless, only a small proportion of patients derive a durable benefit. Mechanisms underlying primary and acquired resistance are still incompletely understood. They comprise tumor-intrinsic factors such as genomic and transcriptomic changes; upregulation of immunosuppressive subsets; T cell exhaustion; and promotion of an immune-tolerant tumor microenvironment. The collection of tumor biopsy at disease progression (PD) is challenging both in clinical and research settings as this often occurs at the time of treatment discontinuation. However, the analysis of these samples can lead to novel strategies to prevent or reverse immune resistance. Thus, the current approach to begin a profiling study with patients at the time of PD on ICI enables access and interrogation of such samples. Methods: IRIS is a prospective, investigator-initiated trial at the Princess Margaret Cancer Centre that aims to extensively characterize the genomic, transcriptomic, epigenetic and immunophenotypic profiles of tumors with primary versus acquired resistance to ICI-based therapy. Primary resistance is defined as PD at the first on-treatment imaging, whereas acquired resistance is defined as PD occurring after an initial partial or complete response or following disease stability lasting ≥6 months. Additional objectives include the evaluation of radiomic parameters on standard radiological imaging, investigation of fecal microbiome, generation of patient-derived organoids and facilitation of data and sample sharing with the research community. The planned samples size is 100 patients. A one-time fresh tumor biopsy, blood and stool samples and archival tissue (when available) are collected at the time of PD on ICI (baseline) from all the participants. Longitudinal blood samples are obtained every 2-3 months (around the time of tumor imaging) until PD in patients receiving a subsequent treatment. Subjects who are not amenable for therapy undergo blood collections at the time of further PD. Molecular characterization of tumor samples includes: DNA/RNA sequencing, Assay of Transposase Accessible Chromatin (ATAC)-sequencing, Cellular Indexing of Transcriptomes and Epitopes (CITE)-sequencing, multiplexed immunohistochemistry and flow cytometry. Results of NGS performed on the first biopsy core are returned to patient and physician. Key eligibility criteria include diagnosis of solid tumor, progression to ICI as the most recent line of treatment and disease amenable to core needle biopsy. The IRIS trial, activated in October 2020, is currently open to enrollment. As of January 2021, 21 patients have been enrolled and a total of 92 tissue cores, 42 blood and 20 stool samples have been collected. Clinical trial information: NCT04243720.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

Abstract SP087: Epigenetic therapy against high-risk breast cancer

Standard-of-care therapy in high-risk stage II-III ER+/HER2-negative and Triple-Negative Breast Cancer (TNBC) includes neoadjuvant chemotherapy (NAC) prior to surgery, which is then followed by adjuvant endocrine (for ER+) +/- salvage chemotherapy (for TNBC). These treatments are delivered with limited stratification by disease histology and, while many patients are cured, more than one third of those in whom invasive disease remains in the breast at the time of surgery experience early relapses of metastatic disease within the first few years. The poor survival rate stems from incomplete response to NAC and the lack of alternative tailored therapy against the chemo-resistant disease that recurs. Hence, there is an urgent need for new therapies to treat recurrence and improve outcome. Tumor progression following standard-of-care therapy arises from pre-existing resistant and/or adapting persister cancer cells committing to an expansion phase. These differ from treatment sensitive cells in the epigenetic landscape of their genome, regulating the on/off state of genes for instance. This underlies epigenetic variants, akin to genetic variants, that define vulnerabilities to new treatment strategies. Here, we demonstrate the value of epigenetic therapy in preclinical settings against high-risk breast cancer, putting forward new treatment options to be tested in clinical trials.

Citation Format: M Lupien. Epigenetic therapy against high-risk breast cancer [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr SP087.

PRMT5 inhibition disrupts splicing and stemness in glioblastoma

Glioblastoma (GBM) is a deadly cancer in which cancer stem cells (CSCs) sustain tumor growth and contribute to therapeutic resistance. Protein arginine methyltransferase 5 (PRMT5) has recently emerged as a promising target in GBM. Using two orthogonal-acting inhibitors of PRMT5 (GSK591 or LLY-283), we show that pharmacological inhibition of PRMT5 suppresses the growth of a cohort of 46 patient-derived GBM stem cell cultures, with the proneural subtype showing greater sensitivity. We show that PRMT5 inhibition causes widespread disruption of splicing across the transcriptome, particularly affecting cell cycle gene products. We identify a GBM splicing signature that correlates with the degree of response to PRMT5 inhibition. Importantly, we demonstrate that LLY-283 is brain-penetrant and significantly prolongs the survival of mice with orthotopic patient-derived xenografts. Collectively, our findings provide a rationale for the clinical development of brain penetrant PRMT5 inhibitors as treatment for GBM.

Gradient of Developmental and Injury Response transcriptional states defines functional vulnerabilities underpinning glioblastoma heterogeneity

Glioblastomas harbor diverse cell populations, including rare glioblastoma stem cells (GSCs) that drive tumorigenesis. To characterize functional diversity within this population, we performed single-cell RNA sequencing on >69,000 GSCs cultured from the tumors of 26 patients. We observed a high degree of inter- and intra-GSC transcriptional heterogeneity that could not be fully explained by DNA somatic alterations. Instead, we found that GSCs mapped along a transcriptional gradient spanning two cellular states reminiscent of normal neural development and inflammatory wound response. Genome-wide CRISPR-Cas9 dropout screens independently recapitulated this observation, with each state characterized by unique essential genes. Further single-cell RNA sequencing of >56,000 malignant cells from primary tumors found that the majority organize along an orthogonal astrocyte maturation gradient yet retain expression of founder GSC transcriptional programs. We propose that glioblastomas grow out of a fundamental GSC-based neural wound response transcriptional program, which is a promising target for new therapy development.

Single-cell chromatin accessibility profiling of glioblastoma identifies an invasive cancer stem cell population associated with lower survival

Chromatin accessibility discriminates stem from mature cell populations, enabling the identification of primitive stem-like cells in primary tumors, such as glioblastoma (GBM) where self-renewing cells driving cancer progression and recurrence are prime targets for therapeutic intervention. We show, using single-cell chromatin accessibility, that primary human GBMs harbor a heterogeneous self-renewing population whose diversity is captured in patient-derived glioblastoma stem cells (GSCs). In-depth characterization of chromatin accessibility in GSCs identifies three GSC states: Reactive, Constructive, and Invasive, each governed by uniquely essential transcription factors and present within GBMs in varying proportions. Orthotopic xenografts reveal that GSC states associate with survival, and identify an invasive GSC signature predictive of low patient survival, in line with the higher invasive properties of Invasive state GSCs compared to Reactive and Constructive GSCs as shown by in vitro and in vivo assays. Our chromatin-driven characterization of GSC states improves prognostic precision and identifies dependencies to guide combination therapies.

The Transition from Quiescent to Activated States in Human Hematopoietic Stem Cells Is Governed by Dynamic 3D Genome Reorganization

Lifelong blood production requires long-term hematopoietic stem cells (LT-HSCs), marked by stemness states involving quiescence and self-renewal, to transition into activated short-term HSCs (ST-HSCs) with reduced stemness. As few transcriptional changes underlie this transition, we used single-cell and bulk assay for transposase-accessible chromatin sequencing (ATAC-seq) on human HSCs and hematopoietic stem and progenitor cell (HSPC) subsets to uncover chromatin accessibility signatures, one including LT-HSCs (LT/HSPC signature) and another excluding LT-HSCs (activated HSPC [Act/HSPC] signature). These signatures inversely correlated during early hematopoietic commitment and differentiation. The Act/HSPC signature contains CCCTC-binding factor (CTCF) binding sites mediating 351 chromatin interactions engaged in ST-HSCs, but not LT-HSCs, enclosing multiple stemness pathway genes active in LT-HSCs and repressed in ST-HSCs. CTCF silencing derepressed stemness genes, restraining quiescent LT-HSCs from transitioning to activated ST-HSCs. Hence, 3D chromatin interactions centrally mediated by CTCF endow a gatekeeper function that governs the earliest fate transitions HSCs make by coordinating disparate stemness pathways linked to quiescence and self-renewal.

GLUT1 inhibition blocks growth of RB1-positive triple negative breast cancer

Triple negative breast cancer (TNBC) is a deadly form of breast cancer due to the development of resistance to chemotherapy affecting over 30% of patients. New therapeutics and companion biomarkers are urgently needed. Recognizing the elevated expression of glucose transporter 1 (GLUT1, encoded by SLC2A1) and associated metabolic dependencies in TNBC, we investigated the vulnerability of TNBC cell lines and patient-derived samples to GLUT1 inhibition. We report that genetic or pharmacological inhibition of GLUT1 with BAY-876 impairs the growth of a subset of TNBC cells displaying high glycolytic and lower oxidative phosphorylation (OXPHOS) rates. Pathway enrichment analysis of gene expression data suggests that the functionality of the E2F pathway may reflect to some extent OXPHOS activity. Furthermore, the protein levels of retinoblastoma tumor suppressor (RB1) strongly correlate with the degree of sensitivity to GLUT1 inhibition in TNBC, where RB1-negative cells are insensitive to GLUT1 inhibition. Collectively, our results highlight a strong and targetable RB1-GLUT1 metabolic axis in TNBC and warrant clinical evaluation of GLUT1 inhibition in TNBC patients stratified according to RB1 protein expression levels.

Abstract 3438: Hypoxia inhibits BMP4-induced astrocytic differentiation in glioma stem cells

The tumor microenvironment is an important contributor to malignancy in glioblastoma (GBM) by influencing stemness properties in glioma stem cells (GSCs) and their capacity for differentiation into various cell types. Tumor hypoxia is a common feature of the GBM microenvironment and can influence cell state through a variety of processes including the inhibition of oxygen-dependent lysine histone demethylases (KDMs) and subsequent epigenetic changes required for differentiation or maintenance of the stem cell state. However, the precise mechanisms underlying hypoxic restriction of differentiation to specific neural lineages is unclear. Here, we show that hypoxia suppresses astrocytic differentiation induced by bone morphogenetic protein 4 (BMP4) in patient-derived GSCs. In 8/10 patient derived GSC lines, BMP4 induced expression of astrocytic markers, such as glial fibrillary acidic protein, and this induction was repressed by hypoxia in 7/10 lines. Interestingly, significant anti-proliferative effects of BMP4 induced differentiation were only observed in 3/10 lines and was rescued in only one GSC line by hypoxia, suggesting that these processes are uncoupled. Although hypoxia may directly limit the activity of oxygen-dependent KDMs important in this process, levels of S-2-hydroxyglutarate (S-2HG), an endogenous KDM inhibitor, were also increased in 5/8 GSC lines under hypoxia. Exogenous treatment with S-2HG or a small molecule KDM6 family inhibitor, phenocopied the repressive effects of hypoxia on BMP4 induced differentiation. In patient tumors, hypoxic areas surrounding pseudopalisading necrosis showed high expression of hypoxia markers, and these hypoxia markers were anti-correlated with markers of astrocytic differentiation. Together, these results demonstrate that the hypoxic tumor microenvironment in GBM, through multiple mechanisms limiting oxygen-dependent KDM activity, promote maintenance of the malignant stemness state by inhibiting astrocytic differentiation. This study implicates tumour hypoxia as a key regulator of stemness and a therapeutic target to promote differentiation in GBM. Treatment strategies that target hypoxic cells are urgently needed to overcome these challenges.

Citation Format: Ronald S. Wu, Elisabeth Liedtke, Sheila Mansouri, Samantha Brazier, Paul Guilhamon, Nicole I. Park, Eric Chen, Mathieu Lupien, Daniel D. De Carvalho, Peter G. Dirks, Gelareh Zadeh, Bradley G. Wouters. Hypoxia inhibits BMP4-induced astrocytic differentiation in glioma stem cells [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3438.

Abstract 5724: Revelation of shared transcriptional gradients in glioblastoma tumors and cultured glioma stem cells

The brain tumor Glioblastoma (GBM) is a virtual death sentence for anyone who is diagnosed, with a median survival time of 12-15 months with standard treatments1 and a 5 year survival rate of under 10% 2. There is a strong body of evidence supporting the existence of stem like cells, termed glioma stem cells (GSCs), that can repopulate the tumor after removal and therapy application3-6. Thus, GSCs present a tantalizing target for potential therapies against GBM. However, studies of GSCs have shown heterogeneity at the level of the transcriptome and drug response5,7, suggesting that there is significant biological variation that translates into differential sensitivity to various drugs. A full characterization of biological heterogeneity may aid in the search for targeted therapies.

Here, we profile the transcriptomes of 72 patient derived GSC cultures, and obtain scRNA-seq on 29 cultures from 26 patients (&gt; 69,000 cells) plus 5 GBM tumors (&gt; 14,000 cells). With this data, we find two anticorrelated transcriptional programs in the GSC cultures, one associated with immune or injury response related pathways and the other with neural developmental pathways. We then compare the GSC cultures to patient tumors in the scRNA-seq data, and find that a portion of GBM tumor cells are similar to GSC cultures. We find that a gradient between GSC and astrocyte programs separates cells with stemness properties from those that are not stem-like, and that within stem-like tumor cells and non stem-like tumor cells a gradient between the neural developmental and immune related programs exists as was seen for the cultured GSCs. In GSCs, we also find epigenetic variation in DNA methylation associated with the developmental and immune related programs. Overall these data suggest variation between two biological programs manifests at the level of gene expression and epigenetic regulation in GSCs in tumors.

1. Stupp, R. et al. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N. Engl. J. Med. 352, 987-996 (2005).

2. Brennan, C. W. et al. The somatic genomic landscape of glioblastoma. Cell 155, 462-477 (2013).

3. Singh, S. K. et al. Identification of human brain tumour initiating cells. Nature 432, 396-401 (2004).

4. Chen, J. et al. A restricted cell population propagates glioblastoma growth after chemotherapy. Nature 488, 522-526 (2012).

5. Lan, X. et al. Fate mapping of human glioblastoma reveals an invariant stem cell hierarchy. Nature 549, 227-232 (2017).

6. Patel, A. P. et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 344, 1396-1401 (2014).

7. Meyer, M. et al. Single cell-derived clonal analysis of human glioblastoma links functional and genomic heterogeneity. Proc. Natl. Acad. Sci. U. S. A. 112, 851-856 (2015).

Citation Format: Owen K. Whitley, Laura M. Richards, Florence Cavalli, Paul Guilhamon, Fiona Coutinho, Michelle Kushida, H. Artee Luchman, Samuel Weiss, Mathieu Lupien, Peter Dirks, Trevor Pugh, Gary Bader. Revelation of shared transcriptional gradients in glioblastoma tumors and cultured glioma stem cells [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5724.

Metabolic Regulation of the Epigenome Drives Lethal Infantile Ependymoma

Posterior fossa A (PFA) ependymomas are lethal malignancies of the hindbrain in infants and toddlers. Lacking highly recurrent somatic mutations, PFA ependymomas are proposed to be epigenetically driven tumors for which model systems are lacking. Here we demonstrate that PFA ependymomas are maintained under hypoxia, associated with restricted availability of specific metabolites to diminish histone methylation, and increase histone demethylation and acetylation at histone 3 lysine 27 (H3K27). PFA ependymomas initiate from a cell lineage in the first trimester of human development that resides in restricted oxygen. Unlike other ependymomas, transient exposure of PFA cells to ambient oxygen induces irreversible cellular toxicity. PFA tumors exhibit a low basal level of H3K27me3, and, paradoxically, inhibition of H3K27 methylation specifically disrupts PFA tumor growth. Targeting metabolism and/or the epigenome presents a unique opportunity for rational therapy for infants with PFA ependymoma.

Disrupting Mitochondrial Copper Distribution Inhibits Leukemic Stem Cell Self-Renewal

Leukemic stem cells (LSCs) rely on oxidative metabolism and are differentially sensitive to targeting mitochondrial pathways, which spares normal hematopoietic cells. A subset of mitochondrial proteins is folded in the intermembrane space via the mitochondrial intermembrane assembly (MIA) pathway. We found increased mRNA expression of MIA pathway substrates in acute myeloid leukemia (AML) stem cells. Therefore, we evaluated the effects of inhibiting this pathway in AML. Genetic and chemical inhibition of ALR reduces AML growth and viability, disrupts LSC self-renewal, and induces their differentiation. ALR inhibition preferentially decreases its substrate COX17, a mitochondrial copper chaperone, and knockdown of COX17 phenocopies ALR loss. Inhibiting ALR and COX17 increases mitochondrial copper levels which in turn inhibit S-adenosylhomocysteine hydrolase (SAHH) and lower levels of S-adenosylmethionine (SAM), DNA methylation, and chromatin accessibility to lower LSC viability. These results provide insight into mechanisms through which mitochondrial copper controls epigenetic status and viability of LSCs.

Telomere dysfunction cooperates with epigenetic alterations to impair murine embryonic stem cell fate commitment

The precise relationship between epigenetic alterations and telomere dysfunction is still an extant question. Previously, we showed that eroded telomeres lead to differentiation instability in murine embryonic stem cells (mESCs) via DNA hypomethylation at pluripotency-factor promoters. Here, we uncovered that telomerase reverse transcriptase null (Tert^-/-) mESCs exhibit genome-wide alterations in chromatin accessibility and gene expression during differentiation. These changes were accompanied by an increase of H3K27me3 globally, an altered chromatin landscape at the Pou5f1/Oct4 promoter, and a refractory response to differentiation cues. Inhibition of the Polycomb Repressive Complex 2 (PRC2), an H3K27 tri-methyltransferase, exacerbated the impairment in differentiation and pluripotency gene repression in Tert^-/-mESCs but not wild-type mESCs, whereas inhibition of H3K27me3 demethylation led to a partial rescue of the Tert^-/- phenotype. These data reveal a new interdependent relationship between H3K27me3 and telomere integrity in stem cell lineage commitment that may have implications in aging and cancer.

Author Correction: Human somatic cell mutagenesis creates genetically tractable sarcomas

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

Noncoding mutations target cis-regulatory elements of the FOXA1 plexus in prostate cancer

Prostate cancer is the second most commonly diagnosed malignancy among men worldwide. Recurrently mutated in primary and metastatic prostate tumors, FOXA1 encodes a pioneer transcription factor involved in disease onset and progression through both androgen receptor-dependent and androgen receptor-independent mechanisms. Despite its oncogenic properties however, the regulation of FOXA1 expression remains unknown. Here, we identify a set of six cis-regulatory elements in the FOXA1 regulatory plexus harboring somatic single-nucleotide variants in primary prostate tumors. We find that deletion and repression of these cis-regulatory elements significantly decreases FOXA1 expression and prostate cancer cell growth. Six of the ten single-nucleotide variants mapping to FOXA1 regulatory plexus significantly alter the transactivation potential of cis-regulatory elements by modulating the binding of transcription factors. Collectively, our results identify cis-regulatory elements within the FOXA1 plexus mutated in primary prostate tumors as potential targets for therapeutic intervention.

FOXA1 pioneer transcription factor is recurrently mutated in primary and metastatic prostate tumors. Here, authors identify a set of six cis-regulatory elements in the FOXA1 regulatory plexus harboring somatic SNVs in primary prostate tumors and characterize their role in regulating FOXA1 expression and prostate cancer cell growth.

Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

Abstract B46: Targeting pancreatic cancer organoids with dual BET and CBP/P300 inhibitor NEO2734

Pancreatic ductal adenocarcinoma (PDA) has a dismal prognosis due to largely ineffective therapeutics available for this disease. Recently, inhibitors of bromodomain and extraterminal domain (BET) inhibitors have been shown to have therapeutic efficacy against various malignancies, including PDA. Here we present data on sensitivity of PDA patient-derived organoids (PDO) to a novel dual BET and CBP/P300 inhibitor, NEO2734 (Epigene Therapeutics Inc.). Our cell viability screen on a panel of 30 PDOs derived from PDA patient and xenograft tumor tissue revealed differential drug sensitivities to NEO2734 inhibitor. Antitumor effect of NEO2734 on the most sensitive PPTO.2 model was further confirmed in vivo using patient-derived xenograft PPTO.2 model. To identify potential biomarkers associated with differential sensitivities of PDA PDOs to NEO2734 inhibitor, we generated gene expression profiles from 13 PDOs derived from primary patient tumors. Using the probabilistic concordance index (C-index), 13 genes (Bonferroni corrected P&lt;1.0E-10) were identified as potential therapeutic targets. Among the biomarkers identified, we found three associations with known BET inhibitor target genes, including CA9 (C-index= 0.38, PBonferroni= 2.7E-19), which is known to be downregulated by BET inhibitor and coexpressed with the target gene BRDT. Furthermore, the biomarker SLC7A2 (0.28, 3.4E-37) has a known genetic interaction with BRD4, and ABCB6 (0.31, 1.5E-27) is predicted to be coexpressed with CREBBP. Overall, this work reveals that dual BET and CBP/P300 inhibitor, NEO2734, could be a novel therapeutic option for a subset of PDA patients.

Citation Format: Nikolina Radulovich, Laura Tamblyn, Heewon Seo, Benjamin Haibe-Kains, Mathieu Lupien, Francis Gilles, Ming S. Tsao. Targeting pancreatic cancer organoids with dual BET and CBP/P300 inhibitor NEO2734 [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer: Advances in Science and Clinical Care; 2019 Sept 6-9; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2019;79(24 Suppl):Abstract nr B46.

Cistrome Partitioning Reveals Convergence of Somatic Mutations and Risk Variants on Master Transcription Regulators in Primary Prostate Tumors

Thousands of noncoding somatic single-nucleotide variants (SNVs) of unknown function are reported in tumors. Partitioning the genome according to cistromes reveals the enrichment of somatic SNVs in prostate tumors as opposed to adjacent normal tissue cistromes of master transcription regulators, including AR, FOXA1, and HOXB13. This parallels enrichment of prostate cancer genetic predispositions over these transcription regulators' tumor cistromes, exemplified at the 8q24 locus harboring both risk variants and somatic SNVs in cis-regulatory elements upregulating MYC expression. However, Massively Parallel Reporter Assays reveal that few SNVs can alter the transactivation potential of individual cis-regulatory elements. Instead, similar to inherited risk variants, SNVs accumulate in cistromes of master transcription regulators required for prostate cancer development.

Functional Enhancers Shape Extrachromosomal Oncogene Amplifications

Non-coding regions amplified beyond oncogene borders have largely been ignored. Using a computational approach, we find signatures of significant co-amplification of non-coding DNA beyond the boundaries of amplified oncogenes across five cancer types. In glioblastoma, EGFR is preferentially co-amplified with its two endogenous enhancer elements active in the cell type of origin. These regulatory elements, their contacts, and their contribution to cell fitness are preserved on high-level circular extrachromosomal DNA amplifications. Interrogating the locus with a CRISPR interference screening approach reveals a diversity of additional elements that impact cell fitness. The pattern of fitness dependencies mirrors the rearrangement of regulatory elements and accompanying rewiring of the chromatin topology on the extrachromosomal amplicon. Our studies indicate that oncogene amplifications are shaped by regulatory dependencies in the non-coding genome.

Identifying clusters of cis-regulatory elements underpinning TAD structures and lineage-specific regulatory networks

Cellular identity relies on cell-type-specific gene expression controlled at the transcriptional level by cis-regulatory elements (CREs). CREs are unevenly distributed across the genome, giving rise to individual CREs and clusters of CREs (COREs). Technical and biological features hinder CORE identification. We addressed these issues by developing an unsupervised machine learning approach termed clustering of genomic regions analysis method (CREAM). CREAM automates CORE detection from chromatin accessibility profiles that are enriched in CREs strongly bound by master transcription regulators, proximal to highly expressed and essential genes, and discriminating cell identity. Although COREs share similarities with super-enhancers, we highlight differences in terms of the genomic distribution and structure of these cis-regulatory units. We further show the enhanced value of COREs over super-enhancers to identify master transcription regulators, highly expressed and essential genes defining cell identity. COREs enrich at topologically associated domain (TAD) boundaries. They are also preferentially bound by the chromatin looping factors CTCF and cohesin, in contrast to super-enhancers, forming clusters of CTCF and cohesin binding regions and defining homotypic clusters of transcription regulator binding regions (HCTs). Finally, we show the clinical utility of CREAM to identify COREs across chromatin accessibility profiles to stratify more than 400 tumor samples according to their cancer type and to delineate cancer type-specific active biological pathways. Collectively, our results support the utility of CREAM to delineate COREs underlying, with greater accuracy than individual CREs or super-enhancers, the cell-type-specific biological underpinning across a wide range of normal and cancer cell types.

Genome-wide germline correlates of the epigenetic landscape of prostate cancer

Oncogenesis is driven by germline, environmental and stochastic factors. It is unknown how these interact to produce the molecular phenotypes of tumors. We therefore quantified the influence of germline polymorphisms on the somatic epigenome of 589 localized prostate tumors. Predisposition risk loci influence a tumor's epigenome, uncovering a mechanism for cancer susceptibility. We identified and validated 1,178 loci associated with altered methylation in tumoral but not nonmalignant tissue. These tumor methylation quantitative trait loci influence chromatin structure, as well as RNA and protein abundance. One prominent tumor methylation quantitative trait locus is associated with AKT1 expression and is predictive of relapse after definitive local therapy in both discovery and validation cohorts. These data reveal intricate crosstalk between the germ line and the epigenome of primary tumors, which may help identify germline biomarkers of aggressive disease to aid patient triage and optimize the use of more invasive or expensive diagnostic assays.

High-resolution structural genomics reveals new therapeutic vulnerabilities in glioblastoma

We investigated the role of 3D genome architecture in instructing functional properties of glioblastoma stem cells (GSCs) by generating sub-5-kb resolution 3D genome maps by in situ Hi-C. Contact maps at sub-5-kb resolution allow identification of individual DNA loops, domain organization, and large-scale genome compartmentalization. We observed differences in looping architectures among GSCs from different patients, suggesting that 3D genome architecture is a further layer of inter-patient heterogeneity for glioblastoma. Integration of DNA contact maps with chromatin and transcriptional profiles identified specific mechanisms of gene regulation, including the convergence of multiple super enhancers to individual stemness genes within individual cells. We show that the number of loops contacting a gene correlates with elevated transcription. These results indicate that stemness genes are hubs of interaction between multiple regulatory regions, likely to ensure their sustained expression. Regions of open chromatin common among the GSCs tested were poised for expression of immune-related genes, including CD276. We demonstrate that this gene is co-expressed with stemness genes in GSCs and that CD276 can be targeted with an antibody-drug conjugate to eliminate self-renewing cells. Our results demonstrate that integrated structural genomics data sets can be employed to rationally identify therapeutic vulnerabilities in self-renewing cells.

Pervasive H3K27 Acetylation Leads to ERV Expression and a Therapeutic Vulnerability in H3K27M Gliomas

High-grade gliomas defined by histone 3 K27M driver mutations exhibit global loss of H3K27 trimethylation and reciprocal gain of H3K27 acetylation, respectively shaping repressive and active chromatin landscapes. We generated tumor-derived isogenic models bearing this mutation and show that it leads to pervasive H3K27ac deposition across the genome. In turn, active enhancers and promoters are not created de novo and instead reflect the epigenomic landscape of the cell of origin. H3K27ac is enriched at repeat elements, resulting in their increased expression, which in turn can be further amplified by DNA demethylation and histone deacetylase inhibitors providing an exquisite therapeutic vulnerability. These agents may therefore modulate anti-tumor immune responses as a therapeutic modality for this untreatable disease.

Targeting bivalency de-represses Indian Hedgehog and inhibits self-renewal of colorectal cancer-initiating cells

In embryonic stem cells, promoters of key lineage-specific differentiation genes are found in a bivalent state, having both activating H3K4me3 and repressive H3K27me3 histone marks, making them poised for transcription upon loss of H3K27me3. Whether cancer-initiating cells (C-ICs) have similar epigenetic mechanisms that prevent lineage commitment is unknown. Here we show that colorectal C-ICs (CC-ICs) are maintained in a stem-like state through a bivalent epigenetic mechanism. Disruption of the bivalent state through inhibition of the H3K27 methyltransferase EZH2, resulted in decreased self-renewal of patient-derived C-ICs. Epigenomic analyses revealed that the promoter of Indian Hedgehog (IHH), a canonical driver of normal colonocyte differentiation, exists in a bivalent chromatin state. Inhibition of EZH2 resulted in de-repression of IHH, decreased self-renewal, and increased sensitivity to chemotherapy in vivo. Our results reveal an epigenetic block to differentiation in CC-ICs and demonstrate the potential for epigenetic differentiation therapy of a solid tumour through EZH2 inhibition.

The Mitochondrial Transacylase, Tafazzin, Regulates for AML Stemness by Modulating Intracellular Levels of Phospholipids

Tafazzin (TAZ) is a mitochondrial transacylase that remodels the mitochondrial cardiolipin into its mature form. Through a CRISPR screen, we identified TAZ as necessary for the growth and viability of acute myeloid leukemia (AML) cells. Genetic inhibition of TAZ reduced stemness and increased differentiation of AML cells both in vitro and in vivo. In contrast, knockdown of TAZ did not impair normal hematopoiesis under basal conditions. Mechanistically, inhibition of TAZ decreased levels of cardiolipin but also altered global levels of intracellular phospholipids, including phosphatidylserine, which controlled AML stemness and differentiation by modulating toll-like receptor (TLR) signaling.

The Proteogenomic Landscape of Curable Prostate Cancer

DNA sequencing has identified recurrent mutations that drive the aggressiveness of prostate cancers. Surprisingly, the influence of genomic, epigenomic, and transcriptomic dysregulation on the tumor proteome remains poorly understood. We profiled the genomes, epigenomes, transcriptomes, and proteomes of 76 localized, intermediate-risk prostate cancers. We discovered that the genomic subtypes of prostate cancer converge on five proteomic subtypes, with distinct clinical trajectories. ETS fusions, the most common alteration in prostate tumors, affect different genes and pathways in the proteome and transcriptome. Globally, mRNA abundance changes explain only ∼10% of protein abundance variability. As a result, prognostic biomarkers combining genomic or epigenomic features with proteomic ones significantly outperform biomarkers comprised of a single data type.

Abstract ES10-2: Mining noncoding mutations for drivers of development in ER alpha-positive breast cancer

Large-scale analysis of whole genome sequencing of breast tumors has identified thousands of genetic alterations, such as single-nucleotide variants (SNVs), outside of gene coding sequences. Delineating the functional impact of these noncoding SNVs is challenging, limiting their inclusion in precision genomics medicine pipelines. Pan-cancer attempts to detect evidence for positive selection within individual CREs identified promoters of a limited number of genes such as TERT. By considering the enrichment of noncoding mutations in cistromes as opposed to individual CREs, we identify the enrichment of noncoding SNVs in cistromes of breast cancer driver transcription factors, including FOXA1. This parallels enrichment of breast cancer genetic predispositions over these transcription factor cistromes. Furthermore, somatic SNVs accumulate in the regulatory plexus of some of these driver transcription factors providing evidence for a multi-targeted approach to disrupt their activity in breast cancer. In conclusion, our results support a functional impact of somatic SNVs converging on the cistrome of driver transcription factors in breast cancer. Altogether, we show how to interrogate noncoding SNVs to delineate the oncogenic determinants of breast cancer to be considered for therapeutic intervention.

Citation Format: Lupien M. Mining noncoding mutations for drivers of development in ER alpha-positive breast cancer [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr ES10-2.

SMuRF: a novel tool to identify regulatory elements enriched for somatic point mutations

BACKGROUND

Single Nucleotide Variants (SNVs), including somatic point mutations and Single Nucleotide Polymorphisms (SNPs), in noncoding cis-regulatory elements (CREs) can affect gene regulation and lead to disease development. Several approaches have been developed to identify highly mutated regions, but these do not take into account the specific genomic context, and thus likelihood of mutation, of CREs.

RESULTS

Here, we present SMuRF (Significantly Mutated Region Finder), a user-friendly command-line tool to identify these significantly mutated regions from user-defined genomic intervals and SNVs. We demonstrate this using publicly available datasets in which SMuRF identifies 72 significantly mutated CREs in liver cancer, including known mutated gene promoters as well as previously unreported regions.

CONCLUSIONS

SMuRF is a helpful tool to allow the simple identification of significantly mutated regulatory elements. It is open-source and freely available on GitHub ( https://github.com/LupienLab/SMURF ).

Characterization of inv(3) cell line OCI-AML-20 with stroma-dependent CD34 expression

Acute myeloid leukemia (AML) is a complex, heterogeneous disease with variable outcomes following curative intent chemotherapy. AML with inv(3) is a genetic subgroup characterized by a very low response rate to current induction type chemotherapy and thus has among the worst long-term survivorship of the AMLs. Here, we describe OCI-AML-20, a new AML cell line with inv(3) and deletion of chromosome 7; the latter is a common co-occurrence in inv(3) AML. In OCI-AML-20, CD34 expression is maintained and required for repopulation in vitro and in vivo. CD34 expression in OCI-AML-20 shows dependence on the co-culture with stromal cells. Transcriptome analysis indicates that the OCI-AML-20 clusters with other AML patient data sets that have poor prognosis, as well as other AML cell lines, including another inv(3) line, MUTZ-3. OCI-AML-20 is a new cell line resource for studying the biology of inv(3) AML that can be used to identify potential therapies for this poor outcome disease.