Pharmacological induction of chromatin remodeling drives chemosensitization in triple-negative breast cancer

Targeted therapies have improved outcomes for certain cancer subtypes, but cytotoxic chemotherapy remains a mainstay for triple-negative breast cancer (TNBC). The epithelial-to-mesenchymal transition (EMT) is a developmental program co-opted by cancer cells that promotes metastasis and chemoresistance. There are no therapeutic strategies specifically targeting mesenchymal-like cancer cells. We report that the US Food and Drug Administration (FDA)-approved chemotherapeutic eribulin induces ZEB1-SWI/SNF-directed chromatin remodeling to reverse EMT that curtails the metastatic propensity of TNBC preclinical models. Eribulin induces mesenchymal-to-epithelial transition (MET) in primary TNBC in patients, but conventional chemotherapy does not. In the treatment-naive setting, but not after acquired resistance to other agents, eribulin sensitizes TNBC cells to subsequent treatment with other chemotherapeutics. These findings provide an epigenetic mechanism of action of eribulin, supporting its use early in the disease process for MET induction to prevent metastatic progression and chemoresistance. These findings warrant prospective clinical evaluation of the chemosensitizing effects of eribulin in the treatment-naive setting.

Characterizing control of memory CD8 T cell differentiation by BTB-ZF transcription factor Zbtb20

Members of the BTB-ZF transcription factor family regulate the immune system. Our laboratory identified that family member Zbtb20 contributes to the differentiation, recall responses, and metabolism of CD8 T cells. Here, we report a characterization of the transcriptional and epigenetic signatures controlled by Zbtb20 at single-cell resolution during the effector and memory phases of the CD8 T cell response. Without Zbtb20, transcriptional programs associated with memory CD8 T cell formation were up-regulated throughout the CD8 T response. A signature of open chromatin was associated with genes controlling T cell activation, consistent with the known impact on differentiation. In addition, memory CD8 T cells lacking Zbtb20 were characterized by open chromatin regions with overrepresentation of AP-1 transcription factor motifs and elevated RNA- and protein-level expressions of the corresponding AP-1 components. Finally, we describe motifs and genomic annotations from the DNA targets of Zbtb20 in CD8 T cells identified by cleavage under targets and release under nuclease (CUT&RUN). Together, these data establish the transcriptional and epigenetic networks contributing to the control of CD8 T cell responses by Zbtb20.

Pharmacological Induction of mesenchymal-epithelial transition chemosensitizes breast cancer cells and prevents metastatic progression

The epithelial-mesenchymal transition (EMT) is a developmental program co-opted by tumor cells that aids the initiation of the metastatic cascade. Tumor cells that undergo EMT are relatively chemoresistant, and there are currently no therapeutic avenues specifically targeting cells that have acquired mesenchymal traits. We show that treatment of mesenchymal-like triple-negative breast cancer (TNBC) cells with the microtubule-destabilizing chemotherapeutic eribulin, which is FDA-approved for the treatment of advanced breast cancer, leads to a mesenchymal-epithelial transition (MET). This MET is accompanied by loss of metastatic propensity and sensitization to subsequent treatment with other FDA-approved chemotherapeutics. We uncover a novel epigenetic mechanism of action that supports eribulin pretreatment as a path to MET induction that curtails metastatic progression and the evolution of therapy resistance.

Lineage plasticity enables low-ER luminal tumors to evolve and gain basal-like traits

Stratifying breast cancer into specific molecular or histologic subtypes aids in therapeutic decision-making and predicting outcomes; however, these subtypes may not be as distinct as previously thought. Patients with luminal-like, estrogen receptor (ER)-expressing tumors have better prognosis than patients with more aggressive, triple-negative or basal-like tumors. There is, however, a subset of luminal-like tumors that express lower levels of ER, which exhibit more basal-like features. We have found that breast tumors expressing lower levels of ER, traditionally considered to be luminal-like, represent a distinct subset of breast cancer characterized by the emergence of basal-like features. Lineage tracing of low-ER tumors in the MMTV-PyMT mouse mammary tumor model revealed that basal marker-expressing cells arose from normal luminal epithelial cells, suggesting that luminal-to-basal plasticity is responsible for the evolution and emergence of basal-like characteristics. This plasticity allows tumor cells to gain a new lumino-basal phenotype, thus leading to intratumoral lumino-basal heterogeneity. Single-cell RNA sequencing revealed SOX10 as a potential driver for this plasticity, which is known among breast tumors to be almost exclusively expressed in triple-negative breast cancer (TNBC) and was also found to be highly expressed in low-ER tumors. These findings suggest that basal-like tumors may result from the evolutionary progression of luminal tumors with low ER expression.

S146L in MYC is a context-dependent activating substitution in cancer development

MYC is one of the most dysregulated oncogenes and is thought to be fundamental to tumor formation and/or maintenance in many cancer types. This dominant pro-tumor activity makes MYC an attractive target for cancer therapy. However, MYC is a transcription factor lacking enzymatic activity, and the structure of one of its two domains is unknown e.g., its transactivation domain. Consequently, few direct MYC-targeting therapies have been developed, and none have been successful in the clinic. Nevertheless, significant effort has been devoted to understanding the mechanisms of oncogenic MYC activity with the objective of uncovering novel vulnerabilities of MYC-dependent cancers. These extensive investigations have revealed in detail how MYC translocation, amplification, and other upstream perturbations contribute to MYC activity in cancer. However, missense mutations of the MYC gene have remained relatively understudied for their potential role in MYC-mediated oncogenesis. While the function of several low-frequency mutations in MYC have been described, our understanding of other equally or more frequent mutations is incomplete. Herein, we define the function of a recurrent missense mutation in MYC resulting in the substitution S146L. This mutation enhances the interaction between MYC and its cofactor TRRAP and may enhance oncogenic MYC activity in certain cellular contexts.

Phenotypic heterogeneity driven by plasticity of the intermediate EMT state governs disease progression and metastasis in breast cancer

The epithelial-to-mesenchymal transition (EMT) is frequently co-opted by cancer cells to enhance migratory and invasive cell traits. It is a key contributor to heterogeneity, chemoresistance, and metastasis in many carcinoma types, where the intermediate EMT state plays a critical tumor-initiating role. We isolate multiple distinct single-cell clones from the SUM149PT human breast cell line spanning the EMT spectrum having diverse migratory, tumor-initiating, and metastatic qualities, including three unique intermediates. Using a multiomics approach, we identify CBFβ as a key regulator of metastatic ability in the intermediate state. To quantify epithelial-mesenchymal heterogeneity within tumors, we develop an advanced multiplexed immunostaining approach using SUM149-derived orthotopic tumors and find that the EMT state and epithelial-mesenchymal heterogeneity are predictive of overall survival in a cohort of stage III breast cancer. Our model reveals previously unidentified insights into the complex EMT spectrum and its regulatory networks, as well as the contributions of epithelial-mesenchymal plasticity (EMP) in tumor heterogeneity in breast cancer.

Distinct cytosine modification profiles define epithelial-to-mesenchymal cell-state transitions

Background: Epithelial-to-mesenchymal transition (EMT) is an early step in the invasion-metastasis cascade, involving progression through intermediate cell states. Due to challenges with isolating intermediate cell states, genome-wide cytosine modifications that define transition are not completely understood. Methods: The authors measured multiple DNA cytosine modification marks and chromatin accessibility across clonal populations residing in specific EMT states. Results: Clones exhibiting more intermediate EMT phenotypes demonstrated increased 5-hydroxymethylcytosine and decreased 5-methylcytosine. Open chromatin regions containing increased 5-hydroxymethylcytosine CpG loci were enriched in EMT transcription factor motifs and were associated with Rho GTPases. Conclusion: The results indicate the importance of both distinct and shared epigenetic profiles associated with EMT processes that may be targeted to prevent EMT progression.

Abstract 2115: Multicomponent epigenetic profiling of intermediate epithelial-to-mesenchymal states in triple negative breast cancer identifies distinct regulatory signatures

Epithelial-to-mesenchymal transition (EMT) is an early step in the invasion-metastasis cascade, involving cellular progression through a number of intermediate states through the transition. Due to difficulty in isolating cells with intermediate EMT phenotypes, it has been challenging to investigate their molecular characteristics. Here, we describe multicomponent epigenetic profiles of phenotypes across the EMT spectrum. Six single cell clones on different points of the EMT spectrum were isolated from SUM149 cells, a heterogeneous ER/PR negative inflammatory breast cancer line. Methylation (mC) and hydroxymethylation (hmC) of >860,000 CpGs sites were measured with tandem oxidative bisulfite treatment and Illumina EPIC arrays, global chromatin accessibility was determined using Assay for Transposase Accessible Chromatin (ATAC) seq, and gene expression was measured using RNA-seq. We observed distinct mC and hmC profiles of the EMT intermediate clones compared to clones at either the epithelial or mesenchymal extremes. Intermediate clones demonstrated 17,862 CpGs with significantly higher hmC and 7,903 CpGs with significant hypomethylation, compared to the distal clones (FDR < 0.01). CpGs with higher hmC were enriched in regulatory regions such as promoters (OR = 2.87, p < 0.001), while hypomethylated CpGs were enriched in regions within the gene bodies. High hmC included CpGs located in promoters of key epithelial genes CDH1 and OCLN. In addition, 1,653 differentially accessible chromatin regions from ATAC-seq also were identified in intermediate clones (FDR < 0.1). Differentially accessible regions in the intermediate clones were associated with extracellular matrix organization (q-val = 0.04) in Reactome Pathway analysis. Although a large proportion of differentially accessible regions where differential hmC and mC CpGs were enriched tracked to distal intergenic regions (30.8%) and introns (28.1%), we did not observe significant enrichment of differential cytosine modifications overlapping with the differentially accessible chromatin regions which indicates the importance of multicomponent measures of epigenetic states. In genes with differentially accessible regions and hydroxymethylated CpGs, intermediate clones showed higher gene expression than distal clones. In particular, RUNX3 was observed to harbor all differential chromatin accessibility, hmC and mC CpGs, in intermediate clones. There was a stronger correlation between hmC and gene expression for differentially hydroxymethylated CpGs within accessible regions (R = 0.54, p < 0.01) than within non-accessible regions (R = 0.45, p = 0.02). Our results indicate both shared and distinct epigenetic profiles on the EMT spectrum at the cytosine and chromatin level regulate EMT processes and may be targeted to prevent the progression of EMT.

Citation Format: Min Kyung Lee, Meredith S. Brown, Owen M. Wilkins, Diwakar R. Pattabiraman, Brock C. Christensen. Multicomponent epigenetic profiling of intermediate epithelial-to-mesenchymal states in triple negative breast cancer identifies distinct regulatory signatures [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2115.

Kynurenine-Induced Aryl Hydrocarbon Receptor Signaling in Mice Causes Body Mass Gain, Liver Steatosis, and Hyperglycemia

OBJECTIVE

The aryl hydrocarbon receptor (AHR) plays a key role in obesity. In vitro studies revealed that the tryptophan metabolite kynurenine (Kyn) activates AHR signaling in cultured hepatocytes. The objective of this study was to determine whether Kyn activated the AHR in mice to induce obesity.

METHODS

Mice were fed a low-fat diet and the same diet supplemented with Kyn. Body mass, liver status, and the expression of identified relevant genes were determined.

RESULTS

Kyn caused mice to gain significant body mass, develop fatty liver and hyperglycemia, and increase expression levels of cytochrome P450 1B1 and stearoyl-CoA desaturase 1. The hyperglycemia was accompanied with decreased insulin levels, which may have been due to the repression of genes involved in insulin secretion. Kyn plasma concentrations and BMI were measured in female patients, and a significant association was observed between Kyn and age in patients with obesity but not in patients who were lean.

CONCLUSIONS

Results show that (1) Kyn or a metabolite thereof is a ligand responsible for inducing AHR-based obesity, fatty liver, and hyperglycemia in mice; (2) plasma Kyn levels increase with age in women with obesity but not in lean women; and (3) an activated AHR is necessary but not sufficient to attain obesity, a status that also requires fat in the diet.

Enrichment of CpG island shore region hypermethylation in epigenetic breast field cancerization

While changes in DNA methylation are known to occur early in breast carcinogenesis and the landscape of breast tumour DNA methylation is profoundly altered compared with normal tissue, there have been limited efforts to identify DNA methylation field cancerization effects in histologically normal breast tissue adjacent to tumour. Matched tumour, histologically normal tissue of the ipsilateral breast (ipsilateral-normal), and histologically normal tissue of the contralateral breast (contralateral-normal) were obtained from nine women undergoing bilateral mastectomy. Laser capture microdissection was used to select epithelial cells from normal tissue, and neoplastic cells from tumour for genome-scale measures of DNA methylation with the Illumina HumanMethylationEPIC array. We identified substantially more CpG loci that were differentially methylated between contralateral-normal and tumour (63,271 CpG loci q < 0.01), than between ipsilateral-normal and tumour (38,346 CpG loci q < 0.01). We identified differential methylation in ipsilateral-normal relative to contralateral-normal tissue (9,562 CpG loci p < 0.01). In this comparison, hypomethylated loci were significantly enriched for breast cancer-relevant transcription factor binding sites including those for ESR1, FoxA1, and GATA3 and hypermethylated loci were significantly enriched for CpG island shore regions. In addition, progression of shore hypermethylation was observed in tumours compared to matched ipsilateral normal tissue, and these alterations tracked to several well-established tumour suppressor genes. Our results indicate an epigenetic field effect in surrounding histologically normal tissue. This work offers an opportunity to focus investigations of early DNA methylation alterations in breast carcinogenesis and potentially develop epigenetic biomarkers of disease risk.

ABBREVIATIONS

DCIS: ductal carcinoma in situ; GO: gene ontology; OR: odds ratio; CI: confidence interval; TFBS: transcription factor binding site; LOLA: Locus Overlap Analysis.

Genome-wide characterization of cytosine-specific 5-hydroxymethylation in normal breast tissue

Despite recent evidence that 5-hydroxymethylcytosine (5hmC) possesses roles in gene regulation distinct from 5-methylcytosine (5mC), relatively little is known regarding the functions of 5hmC in mammalian tissues. To address this issue, we utilized an approach combining both paired bisulfite (BS) and oxidative bisulfite (oxBS) DNA treatment, to resolve genome-wide patterns of 5hmC and 5mC in normal breast tissue from disease-free women. Although less abundant than 5mC, 5hmC was differentially distributed, and consistently enriched among breast-specific enhancers and transcriptionally active chromatin. In contrast, regulatory regions associated with transcriptional inactivity, such as heterochromatin and repressed Polycomb regions, were relatively depleted of 5hmC. Gene regions containing abundant 5hmC were significantly associated with lactate oxidation, immune cell function, and prolactin signaling pathways. Furthermore, genes containing abundant 5hmC were enriched among those actively transcribed in normal breast tissue. Finally, in independent data sets, normal breast tissue 5hmC was significantly enriched among CpG loci demonstrated to have altered methylation in pre-invasive breast cancer and invasive breast tumors. Primarily, our findings identify genomic loci containing abundant 5hmC in breast tissues and provide a genome-wide map of nucleotide-level 5hmC in normal breast tissue. Additionally, these data suggest 5hmC may participate in gene regulatory programs that are dysregulated during breast-related carcinogenesis.

MicroRNA-Related Genetic Variants Associated with Survival of Head and Neck Squamous Cell Carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is commonly diagnosed at an advanced stage, and prognosis for such patients is poor. There remains a gap in our understanding of genetic variants related with HNSCC prognosis. miRNA-related single nucleotide polymorphisms (miR-SNPs) are a class of genetic variants with gene-regulatory potential.

METHODS

We used a genome-scale approach and independent patient populations in a two-stage approach to test 40,286 common miR-SNPs for association with HNSCC survival in the discovery population (n = 847), and selected the strongest associations for replication in validation phase cases (n = 1,236). Furthermore, we leveraged miRNA interaction databases and miRNA expression data from The Cancer Genome Atlas, to provide functional insight for the identified and replicated associations.

RESULTS

Joint population analyses identified novel miR-SNPs associated with overall survival in oral and laryngeal cancers. rs1816158, located within long noncoding RNA MIR100HG, was associated with overall survival in oral cavity cancer (HR, 1.56; 95% confidence interval (CI), 1.21-2.00). In addition, expression of MIR100HG-embedded miRNA, miR-100, was significantly associated with overall survival in an independent cohort of HNSCC cases (HR, 1.25; 95% CI, 1.06-1.49). A SNP in the 3'UTR of SH3BP4 (rs56161233) that overlaps predicted miRNA-binding sites and is predicted to disrupt several miRNA-mRNA interactions was associated with overall survival of laryngeal cancer (HR, 2.57; 95% CI, 1.71-3.86).

CONCLUSIONS

This work reveals novel miR-SNPs associated with HNSCC survival, and utilizes miRNA-mRNA interaction and expression data to provide functional support for these associations.

IMPACT

These findings extend our understanding of how genetic variation contributes to HNSCC survival, and may contribute to future prognostic models for improved risk stratification.

Genome-scale identification of microRNA-related SNPs associated with risk of head and neck squamous cell carcinoma

Polymorphisms in microRNAs and their target sites can disrupt microRNA-dependent gene regulation, and have been associated with cancer susceptibility. However, genome-scale analyses of microRNA-related genetic variation in cancer are lacking. We tested the associations of ~40 000 common [minor allele frequency (MAF) ≥5%], microRNA-related single nucleotide polymorphisms (miR-SNPs), with risk of head and neck squamous cell carcinoma (HNSCC) in a discovery population, and validated selected loci in an independent population among a total of 2198 cases and 2180 controls. Joint analyses across the discovery and validation populations revealed six novel miR-SNP associations with risk of HNSCC. An upstream variant of MIR548H4 (rs7834169), replicated its association with overall HNSCC risk as well as risk of oral cavity cancer. Four other variants were specifically associated with oral cavity cancer risk (rs16914640, rs1134367, rs7306991 and rs1373756). 3'UTR variant of HADH, rs221347 and rs4975616, located within known cancer risk locus 5p15.33, were specific to risk of laryngeal cancer. High confidence predicted microRNA binding sites were identified for CLEC2D, LOC37443, KDM8 and HADH overlapping rs16914640, rs7306991, rs1134367 and rs221347, respectively. Furthermore, we identified several microRNA interactions with KDM8 and HADH predicted to be disrupted by genetic variation at rs1134367 and rs221347. These results suggest microRNA-related genetic variation may contribute to the genetic susceptibility of HNSCC, and that more powerful evaluation of this class of genetic variation and their relationship with cancer risk is warranted.

Inhibition of the Proteasome β2 Site Sensitizes Triple-Negative Breast Cancer Cells to β5 Inhibitors and Suppresses Nrf1 Activation

The proteasome inhibitors carfilzomib (Cfz) and bortezomib (Btz) are used successfully to treat multiple myeloma, but have not shown clinical efficacy in solid tumors. Here we show that clinically achievable inhibition of the β5 site of the proteasome by Cfz and Btz does not result in loss of viability of triple-negative breast cancer cell lines. We use site-specific inhibitors and CRISPR-mediated genetic inactivation of β1 and β2 to demonstrate that inhibiting a second site of the proteasome, particularly the β2 site, sensitizes cell lines to Btz and Cfz in vitro and in vivo. Inhibiting both β5 and β2 suppresses production of the soluble, active form of the transcription factor Nrf1 and prevents the recovery of proteasome activity through induction of new proteasomes. These findings provide a strong rationale for the development of dual β5 and β2 inhibitors for the treatment of solid tumors.

Receptor tyrosine kinase ERBB4 mediates acquired resistance to ERBB2 inhibitors in breast cancer cells

Approximately 25% of breast cancers overexpress and depend on the receptor tyrosine kinase ERBB2, one of 4 ERBB family members. Targeted therapies directed against ERBB2 have been developed and used clinically, but many patients continue to develop resistance to such therapies. Although much effort has been focused on elucidating the mechanisms of acquired resistance to ERBB2-targeted therapies, the involvement of ERBB4 remains elusive and controversial. We demonstrate that genetic ablation of ERBB4, but not ERBB1-3, led to apoptosis in lapatinib-resistant cells, suggesting that the efficacy of pan-ERBB inhibitors was, at least in part, mediated by the inhibition of ERBB4. Moreover, ERBB4 was upregulated at the protein level in ERBB2+ breast cancer cell lines selected for acquired lapatinib resistance in vitro and in MMTV-Neu mice following prolonged lapatinib treatment. Knockdown of ERBB4 caused a decrease in AKT phosphorylation in resistant cells but not in sensitive cells, suggesting that ERBB4 activated the PI3K/AKT pathway in lapatinib-resistant cells. Importantly, ERBB4 knockdown triggered apoptosis not only in lapatinib-resistant cells but also in trastuzumab-resistant cells. Our results suggest that although ERBB4 is dispensable for naïve ERBB2+ breast cancer cells, it may play a key role in the survival of ERBB2+ cancer cells after they develop resistance to ERBB2 inhibitors, lapatinib and trastuzumab.