Heterogeneity and transcriptional drivers of triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is a heterogeneous disease with limited treatment options. To characterize TNBC heterogeneity, we defined transcriptional, epigenetic, and metabolic subtypes and subtype-driving super-enhancers and transcription factors by combining functional and molecular profiling with computational analyses. Single-cell RNA sequencing revealed relative homogeneity of the major transcriptional subtypes (luminal, basal, and mesenchymal) within samples. We found that mesenchymal TNBCs share features with mesenchymal neuroblastoma and rhabdoid tumors and that the PRRX1 transcription factor is a key driver of these tumors. PRRX1 is sufficient for inducing mesenchymal features in basal but not in luminal TNBC cells via reprogramming super-enhancer landscapes, but it is not required for mesenchymal state maintenance or for cellular viability. Our comprehensive, large-scale, multiplatform, multiomics study of both experimental and clinical TNBC is an important resource for the scientific and clinical research communities and opens venues for future investigation.

Endocrine Therapy Synergizes with SMAC Mimetics to Potentiate Antigen Presentation and Tumor Regression in Hormone Receptor-Positive Breast Cancer.

Immunotherapies have yet to demonstrate significant efficacy in the treatment of hormone receptor-positive (HR+) breast cancer. Given that endocrine therapy (ET) is the primary approach for treating HR+ breast cancer, we investigated the effects of ET on the tumor immune microenvironment (TME) in HR+ breast cancer. Spatial proteomics of primary HR+ breast cancer samples obtained at baseline and after ET from patients enrolled in a neoadjuvant clinical trial (NCT02764541) indicated that ET upregulated ß2-microglobulin and influenced the TME in a manner that promotes enhanced immunogenicity. To gain a deeper understanding of the underlying mechanisms, the intrinsic effects of ET on cancer cells were explored, which revealed that ET plays a crucial role in facilitating the chromatin binding of RelA, a key component of the NF-κB complex. Consequently, heightened NF-κB signaling enhanced the response to interferon-gamma, leading to the upregulation of ß2-microglobulin and other antigen presentation-related genes. Further, modulation of NF-κB signaling using a SMAC mimetic in conjunction with ET augmented T-cell migration and enhanced MHC-I-specific T-cell-mediated cytotoxicity. Remarkably, the combination of ET and SMAC mimetics, which also blocks prosurvival effects of NF-κB signaling through the degradation of inhibitors of apoptosis proteins, elicited tumor regression through cell autonomous mechanisms, providing additional support for their combined use in HR+ breast cancer. SIGNIFICANCE: Adding SMAC mimetics to endocrine therapy enhances tumor regression in a cell autonomous manner while increasing tumor immunogenicity, indicating that this combination could be an effective treatment for HR+ patients with breast cancer.

Assessment of spatial transcriptomics for oncology discovery.

Tumor heterogeneity is a major challenge for oncology drug discovery and development. Understanding of the spatial tumor landscape is key to identifying new targets and impactful model systems. Here, we test the utility of spatial transcriptomics (ST) for oncology discovery by profiling 40 tissue sections and 80,024 capture spots across a diverse set of tissue types, sample formats, and RNA capture chemistries. We verify the accuracy and fidelity of ST by leveraging matched pathology analysis, which provides a ground truth for tissue section composition. We then use spatial data to demonstrate the capture of key tumor depth features, identifying hypoxia, necrosis, vasculature, and extracellular matrix variation. We also leverage spatial context to identify relative cell-type locations showing the anti-correlation of tumor and immune cells in syngeneic cancer models. Lastly, we demonstrate target identification approaches in clinical pancreatic adenocarcinoma samples, highlighting tumor intrinsic biomarkers and paracrine signaling.

A Distinct Chromatin State Drives Therapeutic Resistance in Invasive Lobular Breast Cancer.

Most invasive lobular breast cancers (ILC) are of the luminal A subtype and are strongly hormone receptor-positive. Yet, ILC is relatively resistant to tamoxifen and associated with inferior long-term outcomes compared with invasive ductal cancers (IDC). In this study, we sought to gain mechanistic insights into these clinical findings that are not explained by the genetic landscape of ILC and to identify strategies to improve patient outcomes. A comprehensive analysis of the epigenome of ILC in preclinical models and clinical samples showed that, compared with IDC, ILC harbored a distinct chromatin state linked to gained recruitment of FOXA1, a lineage-defining pioneer transcription factor. This resulted in an ILC-unique FOXA1-estrogen receptor (ER) axis that promoted the transcription of genes associated with tumor progression and poor outcomes. The ILC-unique FOXA1-ER axis led to retained ER chromatin binding after tamoxifen treatment, which facilitated tamoxifen resistance while remaining strongly dependent on ER signaling. Mechanistically, gained FOXA1 binding was associated with the autoinduction of FOXA1 in ILC through an ILC-unique FOXA1 binding site. Targeted silencing of this regulatory site resulted in the disruption of the feed-forward loop and growth inhibition in ILC. In summary, ILC is characterized by a unique chromatin state and FOXA1-ER axis that is associated with tumor progression, offering a novel mechanism of tamoxifen resistance. These results underscore the importance of conducting clinical trials dedicated to patients with ILC in order to optimize treatments in this breast cancer subtype. SIGNIFICANCE: A unique FOXA1-ER axis in invasive lobular breast cancer promotes disease progression and tamoxifen resistance, highlighting a potential therapeutic avenue for clinical investigations dedicated to this disease. See related commentary by Blawski and Toska, p. 3668.

FGFR-inhibitor-mediated dismissal of SWI/SNF complexes from YAP-dependent enhancers induces adaptive therapeutic resistance.

How cancer cells adapt to evade the therapeutic effects of drugs targeting oncogenic drivers is poorly understood. Here we report an epigenetic mechanism leading to the adaptive resistance of triple-negative breast cancer (TNBC) to fibroblast growth factor receptor (FGFR) inhibitors. Prolonged FGFR inhibition suppresses the function of BRG1-dependent chromatin remodelling, leading to an epigenetic state that derepresses YAP-associated enhancers. These chromatin changes induce the expression of several amino acid transporters, resulting in increased intracellular levels of specific amino acids that reactivate mTORC1. Consistent with this mechanism, addition of mTORC1 or YAP inhibitors to FGFR blockade synergistically attenuated the growth of TNBC patient-derived xenograft models. Collectively, these findings reveal a feedback loop involving an epigenetic state transition and metabolic reprogramming that leads to adaptive therapeutic resistance and provides potential therapeutic strategies to overcome this mechanism of resistance.

Predicting master transcription factors from pan-cancer expression data.

Critical developmental "master transcription factors" (MTFs) can be subverted during tumorigenesis to control oncogenic transcriptional programs. Current approaches to identifying MTFs rely on ChIP-seq data, which is unavailable for many cancers. We developed the CaCTS (Cancer Core Transcription factor Specificity) algorithm to prioritize candidate MTFs using pan-cancer RNA sequencing data. CaCTS identified candidate MTFs across 34 tumor types and 140 subtypes including predictions for cancer types/subtypes for which MTFs are unknown, including e.g. PAX8, SOX17, and MECOM as candidates in ovarian cancer (OvCa). In OvCa cells, consistent with known MTF properties, these factors are required for viability, lie proximal to superenhancers, co-occupy regulatory elements globally, co-bind loci encoding OvCa biomarkers, and are sensitive to pharmacologic inhibition of transcription. Our predictions of MTFs, especially for tumor types with limited understanding of transcriptional drivers, pave the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Subtype heterogeneity and epigenetic convergence in neuroendocrine prostate cancer.

Neuroendocrine carcinomas (NEC) are tumors expressing markers of neuronal differentiation that can arise at different anatomic sites but have strong histological and clinical similarities. Here we report the chromatin landscapes of a range of human NECs and show convergence to the activation of a common epigenetic program. With a particular focus on treatment emergent neuroendocrine prostate cancer (NEPC), we analyze cell lines, patient-derived xenograft (PDX) models and human clinical samples to show the existence of two distinct NEPC subtypes based on the expression of the neuronal transcription factors ASCL1 and NEUROD1. While in cell lines and PDX models these subtypes are mutually exclusive, single-cell analysis of human clinical samples exhibits a more complex tumor structure with subtypes coexisting as separate sub-populations within the same tumor. These tumor sub-populations differ genetically and epigenetically contributing to intra- and inter-tumoral heterogeneity in human metastases. Overall, our results provide a deeper understanding of the shared clinicopathological characteristics shown by NECs. Furthermore, the intratumoral heterogeneity of human NEPCs suggests the requirement of simultaneous targeting of coexisting tumor populations as a therapeutic strategy.

CDK4/6 inhibition reprograms the breast cancer enhancer landscape by stimulating AP-1 transcriptional activity.

Pharmacologic inhibitors of cyclin-dependent kinases 4 and 6 (CDK4/6) were designed to induce cancer cell cycle arrest. Recent studies have suggested that these agents also exert other effects, influencing cancer cell immunogenicity, apoptotic responses, and differentiation. Using cell-based and mouse models of breast cancer together with clinical specimens, we show that CDK4/6 inhibitors induce remodeling of cancer cell chromatin characterized by widespread enhancer activation, and that this explains many of these effects. The newly activated enhancers include classical super-enhancers that drive luminal differentiation and apoptotic evasion, as well as a set of enhancers overlying endogenous retroviral elements that is enriched for proximity to interferon-driven genes. Mechanistically, CDK4/6 inhibition increases the level of several Activator Protein-1 (AP-1) transcription factor proteins, which are in turn implicated in the activity of many of the new enhancers. Our findings offer insights into CDK4/6 pathway biology and should inform the future development of CDK4/6 inhibitors.

A hotspot mutation in transcription factor IKZF3 drives B cell neoplasia via transcriptional dysregulation.

Hotspot mutation of IKZF3 (IKZF3-L162R) has been identified as a putative driver of chronic lymphocytic leukemia (CLL), but its function remains unknown. Here, we demonstrate its driving role in CLL through a B cell-restricted conditional knockin mouse model. Mutant Ikzf3 alters DNA binding specificity and target selection, leading to hyperactivation of B cell receptor (BCR) signaling, overexpression of nuclear factor κB (NF-κB) target genes, and development of CLL-like disease in elderly mice with a penetrance of ~40%. Human CLL carrying either IKZF3 mutation or high IKZF3 expression was associated with overexpression of BCR/NF-κB pathway members and reduced sensitivity to BCR signaling inhibition by ibrutinib. Our results thus highlight IKZF3 oncogenic function in CLL via transcriptional dysregulation and demonstrate that this pro-survival function can be achieved by either somatic mutation or overexpression of this CLL driver. This emphasizes the need for combinatorial approaches to overcome IKZF3-mediated BCR inhibitor resistance.

Therapeutically Increasing MHC-I Expression Potentiates Immune Checkpoint Blockade.

Immune checkpoint blockade (ICB) therapy revolutionized cancer treatment, but many patients with impaired MHC-I expression remain refractory. Here, we combined FACS-based genome-wide CRISPR screens with a data-mining approach to identify drugs that can upregulate MHC-I without inducing PD-L1. CRISPR screening identified TRAF3, a suppressor of the NFκB pathway, as a negative regulator of MHC-I but not PD-L1. The Traf3-knockout gene expression signature is associated with better survival in ICB-naïve patients with cancer and better ICB response. We then screened for drugs with similar transcriptional effects as this signature and identified Second Mitochondria-derived Activator of Caspase (SMAC) mimetics. We experimentally validated that the SMAC mimetic birinapant upregulates MHC-I, sensitizes cancer cells to T cell-dependent killing, and adds to ICB efficacy. Our findings provide preclinical rationale for treating tumors expressing low MHC-I expression with SMAC mimetics to enhance sensitivity to immunotherapy. The approach used in this study can be generalized to identify other drugs that enhance immunotherapy efficacy. SIGNIFICANCE: MHC-I loss or downregulation in cancer cells is a major mechanism of resistance to T cell-based immunotherapies. Our study reveals that birinapant may be used for patients with low baseline MHC-I to enhance ICB response. This represents promising immunotherapy opportunities given the biosafety profile of birinapant from multiple clinical trials.This article is highlighted in the In This Issue feature, p. 1307.

ERG-Mediated Coregulator Complex Formation Maintains Androgen Receptor Signaling in Prostate Cancer.

The TMPRSS2-ERG fusion is the most common genomic rearrangement in human prostate cancer. However, in established adenocarcinoma, it is unknown how the ERG oncogene promotes a cancerous phenotype and maintains downstream androgen receptor (AR) signaling pathways. In this study, we utilized a murine prostate organoid system to explore the effects of ERG on tumorigenesis and determined the mechanism underlying prostate cancer dependence on ERG. Prostate organoids lacking PTEN and overexpressing ERG (Pten-/- R26-ERG) faithfully recapitulated distinct stages of prostate cancer disease progression. In this model, deletion of ERG significantly dampened AR-dependent gene expression. While ERG was able to reprogram the AR cistrome in the process of prostate carcinogenesis, ERG knockout in established prostate cancer organoids did not drastically alter AR binding, H3K27ac enhancer, or open chromatin profiles at these reprogrammed sites. Proteomic analysis of DNA-bound AR complexes demonstrated that ERG deletion causes a loss of recruitment of critical AR coregulators and basal transcriptional machinery, including NCOA3 and RNA polymerase II, but does not alter AR binding itself. Together, these data reveal a novel mechanism of ERG oncogene addiction in prostate cancer, whereby ERG facilitates AR signaling by maintaining coregulator complexes at AR bound sites across the genome. SIGNIFICANCE: These findings exploit murine organoid models to uncover the mechanism of ERG-mediated tumorigenesis and subsequent oncogenic dependencies in prostate cancer.

Prostate cancer reactivates developmental epigenomic programs during metastatic progression.

Epigenetic processes govern prostate cancer (PCa) biology, as evidenced by the dependency of PCa cells on the androgen receptor (AR), a prostate master transcription factor. We generated 268 epigenomic datasets spanning two state transitions-from normal prostate epithelium to localized PCa to metastases-in specimens derived from human tissue. We discovered that reprogrammed AR sites in metastatic PCa are not created de novo; rather, they are prepopulated by the transcription factors FOXA1 and HOXB13 in normal prostate epithelium. Reprogrammed regulatory elements commissioned in metastatic disease hijack latent developmental programs, accessing sites that are implicated in prostate organogenesis. Analysis of reactivated regulatory elements enabled the identification and functional validation of previously unknown metastasis-specific enhancers at HOXB13, FOXA1 and NKX3-1. Finally, we observed that prostate lineage-specific regulatory elements were strongly associated with PCa risk heritability and somatic mutation density. Examining prostate biology through an epigenomic lens is fundamental for understanding the mechanisms underlying tumor progression.

FiTAc-seq: fixed-tissue ChIP-seq for H3K27ac profiling and super-enhancer analysis of FFPE tissues.

Fixed-tissue ChIP-seq for H3K27 acetylation (H3K27ac) profiling (FiTAc-seq) is an epigenetic method for profiling active enhancers and promoters in formalin-fixed, paraffin-embedded (FFPE) tissues. We previously developed a modified ChIP-seq protocol (FiT-seq) for chromatin profiling in FFPE. FiT-seq produces high-quality chromatin profiles particularly for methylated histone marks but is not optimized for H3K27ac profiling. FiTAc-seq is a modified protocol that replaces the proteinase K digestion applied in FiT-seq with extended heating at 65 °C in a higher concentration of detergent and a minimized sonication step, to produce robust genome-wide H3K27ac maps from clinical samples. FiTAc-seq generates high-quality enhancer landscapes and super-enhancer (SE) annotation in numerous archived FFPE samples from distinct tumor types. This approach will be of great interest for both basic and clinical researchers. The entire protocol from FFPE blocks to sequence-ready library can be accomplished within 4 d.

Synthetic Lethal and Resistance Interactions with BET Bromodomain Inhibitors in Triple-Negative Breast Cancer.

BET bromodomain inhibitors (BBDIs) are candidate therapeutic agents for triple-negative breast cancer (TNBC) and other cancer types, but inherent and acquired resistance to BBDIs limits their potential clinical use. Using CRISPR and small-molecule inhibitor screens combined with comprehensive molecular profiling of BBDI response and resistance, we identified synthetic lethal interactions with BBDIs and genes that, when deleted, confer resistance. We observed synergy with regulators of cell cycle progression, YAP, AXL, and SRC signaling, and chemotherapeutic agents. We also uncovered functional similarities and differences among BRD2, BRD4, and BRD7. Although deletion of BRD2 enhances sensitivity to BBDIs, BRD7 loss leads to gain of TEAD-YAP chromatin binding and luminal features associated with BBDI resistance. Single-cell RNA-seq, ATAC-seq, and cellular barcoding analysis of BBDI responses in sensitive and resistant cell lines highlight significant heterogeneity among samples and demonstrate that BBDI resistance can be pre-existing or acquired.

Androgen Receptor Is a Non-canonical Inhibitor of Wild-Type and Mutant Estrogen Receptors in Hormone Receptor-Positive Breast Cancers.

Sustained treatment of estrogen receptor (ER)-positive breast cancer with ER-targeting drugs results in ER mutations and refractory unresponsive cancers. Androgen receptor (AR), which is expressed in 80%-95% of ER-positive breast cancers, could serve as an alternate therapeutic target. Although AR agonists were used in the past to treat breast cancer, their use is currently infrequent due to virilizing side effects. Discovery of tissue-selective AR modulators (SARMs) has renewed interest in using AR agonists to treat breast cancer. Using translational models, we show that AR agonist and SARM, but not antagonist, inhibit the proliferation and growth of ER-positive breast cancer cells, patient-derived tissues, and patient-derived xenografts (PDX). Ligand-activated AR inhibits wild-type and mutant ER activity by reprogramming the ER and FOXA1 cistrome and rendering tumor growth inhibition. These findings suggest that ligand-activated AR may function as a non-canonical inhibitor of ER and that AR agonists may offer a safe and effective treatment for ER-positive breast cancer.

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.

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

ERG overexpression plus SLC45A3 (prostein) and PTEN expression loss: Strong association of the triple hit phenotype with an aggressive pathway of prostate cancer progression.

TMPRSS2 and SLC45A3 rearrangements may coexist in the same tumor. ERG rearrangements and PTEN loss are concomitant events in prostate cancer (PrCa), and can cooperate in progression. We have reported that mRNA expression of TMPRSS2-ERG and SLC45A3-ERG rearrangements plus PTEN loss define an aggressive tumor subset. The aim of this study has been to validate these results by immunohistochemistry in a large cohort of tumors. ERG, SLC45A3 and PTEN immunostaining and their association with pathological features and PSA progression-free survival were analyzed in 220 PrCa (PSMAR-Biobank, Barcelona, Spain). ERG protein expression was found in 46.8% and SLC45A3 and PTEN loss in 30% and 34% tumors, respectively. Single ERG positive immunostaining was associated with GS = 6 tumors (p = 0.016), double ERG+/PTEN loss with GS = 7 (p = 0.008) and Grade Group 2 (GG) or GG3 cases (p = 0.042), ERG+/SLC45A3 loss/PTEN loss ("triple hit") with GS ≥ 8 (p < 0.0001) and GG4 or GG5 tumors (p = 0.0003). None of GS = 6 nor = GG1 cases showed this combination. In the GS ≥ 8 group, ERG+ (p = 0.002), PTEN loss (p = 0.009) and "triple hit" (p = 0.003) were associated with Gleason pattern 3 component, and single SLC45A3 loss (p = 0.036) with GS ≥ 8 without pattern 3. The number of aberrant events and the triple hit were strongly associated with shorter PSA progression-free survival. In GS = 6 PrCa, single ERG+ was also associated with progression. ERG+ identifies a distinct pathway of PrCa. Additional assessment of PTEN and SLC45A3 adds relevant prognostic information. The triple hit phenotype (ERG+/SLC45A3 loss/PTEN loss) is associated with progression and could be used for patient stratification, treatment and follow-up.

Enhancer-Mediated Oncogenic Function of the Menin Tumor Suppressor in Breast Cancer.

While the multiple endocrine neoplasia type 1 (MEN1) gene functions as a tumor suppressor in a variety of cancer types, we explored its oncogenic role in breast tumorigenesis. The MEN1 gene product menin is involved in H3K4 trimethylation and co-activates transcription. We integrated ChIP-seq and RNA-seq data to identify menin target genes. Our analysis revealed that menin-dependent target gene promoters display looping to distal enhancers that are bound by menin, FOXA1 and GATA3. In this fashion, MEN1 co-regulates a proliferative breast cancer-specific gene expression program in ER+ cells. In primary mammary cells, MEN1 exerts an anti-proliferative function by regulating a distinct expression signature. Our findings clarify the cell-type-specific functions of MEN1 and inform the development of menin-directed treatments for breast cancer.

FOXO1 down-regulation is associated with worse outcome in bladder cancer and adds significant prognostic information to p53 overexpression.

Nuclear FOXOs mediate cell cycle arrest and promote apoptosis. FOXOs and p53 could have similar effects as tumor suppressor genes. In spite of extensive literature, little is known about the role of FOXO1 and its relationship with p53 status in bladder cancer. Expression of FOXO1 and p53 were analyzed by immunohistochemistry in 162 urothelial carcinomas (UC). Decreased FOXO1 expression, p53 overexpression and the combination FOXO1 down-regulation/p53 overexpression were strongly associated with high grade (P=.030; P=.017; P=.004, respectively), high stage (P=.0001; P<.0001; P<.0001, respectively) or both (P=.0004; P<.0001; P<.0001, respectively). In the overall series of cases, p53 overexpression was associated with tumor progression (hazard ratio [HR]=3.18, 95% confidence interval [CI] 1.19-8.48, P=.02), but this association was even stronger if having any alteration in any of the 2 genes was considered (HR=3.51, 95% CI 1.34-9.21, P=.01). Having both FOXO1 down-regulation and p53 overexpression was associated with disease recurrence (HR=2.75, 95% CI 1.06-7.13, P=.03). In the analysis of the different subgroups, having any alteration in any of the 2 genes was associated with progression in low-grade (P=.005) and pTa (P=.006) tumors. Finally, the combined FOXO1 down-regulation/p53 overexpression was associated with disease recurrence specifically in high-grade (P=.04) and in pT1 stage tumors (P=.007). Adding FOXO1 expression to the immunohistochemical analysis of p53 can provide relevant prognostic information on progression and recurrence of bladder cancer. It may be particularly informative on the risk of progression in the more indolent and on the risk of recurrence in the more aggressive tumors.