NKG2A Is a Therapeutic Vulnerability in Immunotherapy Resistant MHC-I Heterogeneous Triple-Negative Breast Cancer

Despite the success of immune checkpoint inhibition (ICI) in treating cancer, patients with triple-negative breast cancer (TNBC) often develop resistance to therapy, and the underlying mechanisms are unclear. MHC-I expression is essential for antigen presentation and T-cell-directed immunotherapy responses. This study demonstrates that TNBC patients display intratumor heterogeneity in regional MHC-I expression. In murine models, loss of MHC-I negates antitumor immunity and ICI response, whereas intratumor MHC-I heterogeneity leads to increased infiltration of natural killer (NK) cells in an IFNγ-dependent manner. Using spatial technologies, MHC-I heterogeneity is associated with clinical resistance to anti-programmed death (PD) L1 therapy and increased NK:T-cell ratios in human breast tumors. MHC-I heterogeneous tumors require NKG2A to suppress NK-cell function. Combining anti-NKG2A and anti-PD-L1 therapies restores complete response in heterogeneous MHC-I murine models, dependent on the presence of activated, tumor-infiltrating NK and CD8+ T cells. These results suggest that similar strategies may enhance patient benefit in clinical trials.

SIGNIFICANCE

Clinical resistance to immunotherapy is common in breast cancer, and many patients will likely require combination therapy to maximize immunotherapeutic benefit. This study demonstrates that heterogeneous MHC-I expression drives resistance to anti-PD-L1 therapy and exposes NKG2A on NK cells as a target to overcome resistance. This article is featured in Selected Articles from This Issue, p. 201.

Signaling Pathway Alterations Driven by BRCA1 and BRCA2 Germline Mutations are Sufficient to Initiate Breast Tumorigenesis by the PIK3CAH1047R Oncogene

Single-cell transcriptomics studies have begun to identify breast epithelial cell and stromal cell specific transcriptome differences between BRCA1/2 mutation carriers and non-carriers. We generated a single-cell transcriptome atlas of breast tissues from BRCA1, BRCA2 mutation carriers and compared this single-cell atlas of mutation carriers with our previously described single-cell breast atlas of healthy non-carriers. We observed that BRCA1 but not BRCA2 mutations altered the ratio between basal (basal-myoepithelial), luminal progenitor (luminal adaptive secretory precursor, LASP), and mature luminal (luminal hormone sensing) cells in breast tissues. A unique subcluster of cells within LASP cells is underrepresented in case of BRCA1 and BRCA2 mutation carriers compared with non-carriers. Both BRCA1 and BRCA2 mutations specifically altered transcriptomes in epithelial cells which are an integral part of NFκB, LARP1, and MYC signaling. Signaling pathway alterations in epithelial cells unique to BRCA1 mutations included STAT3, BRD4, SMARCA4, HIF2A/EPAS1, and Inhibin A signaling. BRCA2 mutations were associated with upregulation of IL6, PDK1, FOXO3, and TNFSF11 signaling. These signaling pathway alterations are sufficient to alter sensitivity of BRCA1/BRCA2-mutant breast epithelial cells to transformation as epithelial cells from BRCA1 mutation carriers overexpressing hTERT + PIK3CAH1047R generated adenocarcinomas, whereas similarly modified mutant BRCA2 cells generated basal carcinomas in NSG mice. Thus, our studies provide a high-resolution transcriptome atlas of breast epithelial cells of BRCA1 and BRCA2 mutation carriers and reveal their susceptibility to PIK3CA mutation-driven transformation.

SIGNIFICANCE

This study provides a single-cell atlas of breast tissues of BRCA1/2 mutation carriers and demonstrates that aberrant signaling due to BRCA1/2 mutations is sufficient to initiate breast cancer by mutant PIK3CA.

Peripheral blood TCR clonotype diversity as an age-associated marker of breast cancer progression

Immune escape is a prerequisite for tumor growth. We previously described a decline in intratumor activated cytotoxic T cells and T cell receptor (TCR) clonotype diversity in invasive breast carcinomas compared to ductal carcinoma in situ (DCIS), implying a central role of decreasing T cell responses in tumor progression. To determine potential associations between peripheral immunity and breast tumor progression, here, we assessed the peripheral blood TCR clonotype of 485 breast cancer patients diagnosed with either DCIS or de novo stage IV disease at younger (<45) or older (≥45) age. TCR clonotype diversity was significantly lower in older compared to younger breast cancer patients regardless of tumor stage at diagnosis. In the younger age group, TCR-α clonotype diversity was lower in patients diagnosed with de novo stage IV breast cancer compared to those diagnosed with DCIS. In the older age group, DCIS patients with higher TCR-α clonotype diversity were more likely to have a recurrence compared to those with lower diversity. Whole blood transcriptome profiles were distinct depending on the TCR-α Chao1 diversity score. There were more CD8+ T cells and a more active immune environment in DCIS tumors of young patients with higher peripheral blood TCR-α Chao1 diversity than in those with lower diversity. These results provide insights into the role that host immunity plays in breast cancer development across different age groups.

Combination Therapies to Improve the Efficacy of Immunotherapy in Triple-negative Breast Cancer

Immune checkpoint inhibition combined with chemotherapy is currently approved as first-line treatment for patients with advanced PD-L1-positive triple-negative breast cancer (TNBC). However, a significant proportion of metastatic TNBC is PD-L1-negative and, in this population, chemotherapy alone largely remains the standard-of-care and novel therapeutic strategies are needed to improve clinical outcomes. Here, we describe a triple combination of anti-PD-L1 immune checkpoint blockade, epigenetic modulation thorough bromodomain and extra-terminal (BET) bromodomain inhibition (BBDI), and chemotherapy with paclitaxel that effectively inhibits both primary and metastatic tumor growth in two different syngeneic murine models of TNBC. Detailed cellular and molecular profiling of tumors from single and combination treatment arms revealed increased T- and B-cell infiltration and macrophage reprogramming from MHCIIlow to a MHCIIhigh phenotype in mice treated with triple combination. Triple combination also had a major impact on gene expression and chromatin profiles shifting cells to a more immunogenic and senescent state. Our results provide strong preclinical evidence to justify clinical testing of BBDI, paclitaxel, and immune checkpoint blockade combination.

Kinome Reprogramming Is a Targetable Vulnerability in ESR1 Fusion-Driven Breast Cancer

Transcriptionally active ESR1 fusions (ESR1-TAF) are a potent cause of breast cancer endocrine therapy (ET) resistance. ESR1-TAFs are not directly druggable because the C-terminal estrogen/anti-estrogen-binding domain is replaced with translocated in-frame partner gene sequences that confer constitutive transactivation. To discover alternative treatments, a mass spectrometry (MS)-based kinase inhibitor pulldown assay (KIPA) was deployed to identify druggable kinases that are upregulated by diverse ESR1-TAFs. Subsequent explorations of drug sensitivity validated RET kinase as a common therapeutic vulnerability despite remarkable ESR1-TAF C-terminal sequence and structural diversity. Organoids and xenografts from a pan-ET-resistant patient-derived xenograft model that harbors the ESR1-e6>YAP1 TAF were concordantly inhibited by the selective RET inhibitor pralsetinib to a similar extent as the CDK4/6 inhibitor palbociclib. Together, these findings provide preclinical rationale for clinical evaluation of RET inhibition for the treatment of ESR1-TAF-driven ET-resistant breast cancer.

SIGNIFICANCE

Kinome analysis of ESR1 translocated and mutated breast tumors using drug bead-based mass spectrometry followed by drug-sensitivity studies nominates RET as a therapeutic target. See related commentary by Wu and Subbiah, p. 3159.

Dynamic Glycoprotein Hyposialylation Promotes Chemotherapy Evasion and Metastatic Seeding of Quiescent Circulating Tumor Cell Clusters in Breast Cancer

Most circulating tumor cells (CTC) are detected as single cells, whereas a small proportion of CTCs in multicellular clusters with stemness properties possess 20- to 100-times higher metastatic propensity than the single cells. Here we report that CTC dynamics in both singles and clusters in response to therapies predict overall survival for breast cancer. Chemotherapy-evasive CTC clusters are relatively quiescent with a specific loss of ST6GAL1-catalyzed α2,6-sialylation in glycoproteins. Dynamic hyposialylation in CTCs or deficiency of ST6GAL1 promotes cluster formation for metastatic seeding and enables cellular quiescence to evade paclitaxel treatment in breast cancer. Glycoproteomic analysis reveals newly identified protein substrates of ST6GAL1, such as adhesion or stemness markers PODXL, ICAM1, ECE1, ALCAM1, CD97, and CD44, contributing to CTC clustering (aggregation) and metastatic seeding. As a proof of concept, neutralizing antibodies against one newly identified contributor, PODXL, inhibit CTC cluster formation and lung metastasis associated with paclitaxel treatment for triple-negative breast cancer.

SIGNIFICANCE

This study discovers that dynamic loss of terminal sialylation in glycoproteins of CTC clusters contributes to the fate of cellular dormancy, advantageous evasion to chemotherapy, and enhanced metastatic seeding. It identifies PODXL as a glycoprotein substrate of ST6GAL1 and a candidate target to counter chemoevasion-associated metastasis of quiescent tumor cells. This article is featured in Selected Articles from This Issue, p. 1949.

Assessing shoulder muscle stretch reflexes following breast cancer treatment and postmastectomy breast reconstruction

Muscle stiffness is altered following postmastectomy breast reconstruction and breast cancer treatment. The exact mechanisms underlying these alterations are unknown; however, muscle stretch reflexes may play a role. This work examined short- (SLR) and long-latency (LLR) shoulder muscle stretch reflexes in breast cancer survivors. Forty-nine patients who had undergone postmastectomy breast reconstruction, 17 who had undergone chemoradiation, and 18 healthy, age-matched controls were enrolled. Muscle activity was recorded from the clavicular and sternocostal regions of the pectoralis major and anterior, middle, and posterior deltoids during vertical ab/adduction or horizontal flex/extension perturbations while participants maintained minimal torques. SLR and LLR were quantified for each muscle. Our major finding was that following postmastectomy breast reconstruction, SLR and LLR are impaired in the clavicular region of the pectoralis major. Individuals who had chemoradiation had impaired stretch reflexes in the clavicular and sternocostal region of the pectoralis major, anterior, middle, and posterior deltoid. These findings indicate that breast cancer treatments alter the regulation of shoulder muscle stretch reflexes and may be associated with surgical or nonsurgical damage to the pectoral fascia, muscle spindles, and/or sensory Ia afferents.NEW & NOTEWORTHY Shoulder muscle stretch reflexes may be impacted following postmastectomy breast reconstruction and chemoradiation. Here, we examined short- and long-latency shoulder muscle stretch reflexes in two experiments following common breast reconstruction procedures and chemoradiation. We show impairments in pectoralis major stretch reflexes following postmastectomy breast reconstruction and pectoralis major and deltoid muscle stretch reflexes following chemoradiation. These findings indicate that breast cancer treatments alter the regulation of shoulder muscle stretch reflexes.

Metastasis from the tumor interior and necrotic core formation are regulated by breast cancer-derived angiopoietin-like 7

Necrosis in the tumor interior is a common feature of aggressive cancers that is associated with poor clinical prognosis and the development of metastasis. How the necrotic core promotes metastasis remains unclear. Here, we report that emergence of necrosis inside the tumor is correlated temporally with increased tumor dissemination in a rat breast cancer model and in human breast cancer patients. By performing spatially focused transcriptional profiling, we identified angiopoietin-like 7 (Angptl7) as a tumor-specific factor localized to the perinecrotic zone. Functional studies showed that Angptl7 loss normalizes central necrosis, perinecrotic dilated vessels, metastasis, and reduces circulating tumor cell counts to nearly zero. Mechanistically, Angptl7 promotes vascular permeability and supports vascular remodeling in the perinecrotic zone. Taken together, these findings show that breast tumors actively produce factors controlling central necrosis formation and metastatic dissemination from the tumor core.

Landscape of Genetic Alterations Underlying Hallmark Signature Changes in Cancer Reveals TP53 Aneuploidy-driven Metabolic Reprogramming

The hallmark signatures based on gene expression capture core cancer processes. Through a pan-cancer analysis, we describe the overview of hallmark signatures across tumor types/subtypes and reveal significant relationships between these signatures and genetic alterations. TP53 mutation exerts diverse changes, including increased proliferation and glycolysis, which are closely mimicked by widespread copy-number alterations. Hallmark signature and copy-number clustering identify a cluster of squamous tumors and basal-like breast and bladder cancers with elevated proliferation signatures, frequent TP53 mutation, and high aneuploidy. In these basal-like/squamous TP53-mutated tumors, a specific and consistent spectrum of copy-number alterations is preferentially selected prior to whole-genome duplication. Within Trp53-null breast cancer mouse models, these copy-number alterations spontaneously occur and recapitulate the hallmark signature changes observed in the human condition. Together, our analysis reveals intertumor and intratumor heterogeneity of the hallmark signatures, uncovering an oncogenic program induced by TP53 mutation and select aneuploidy events to drive a worsened prognosis.

Significance

Our data demonstrate that TP53 mutation and a resultant selected pattern of aneuploidies cause an aggressive transcriptional program including upregulation of glycolysis signature with prognostic implications. Importantly, basal-like breast cancer demonstrates genetic and/or phenotypic changes closely related to squamous tumors including 5q deletion that reveal alterations that could offer therapeutic options across tumor types regardless of tissue of origin.

Incorporating RNA-based Risk Scores for Genomic Instability to Predict Breast Cancer Recurrence and Immunogenicity in a Diverse Population

Markers of genomic instability, including TP53 status and homologous recombination deficiency (HRD), are candidate biomarkers of immunogenicity and immune-mediated survival, but little is known about the distribution of these markers in large, population-based cohorts of racially diverse patients with breast cancer. In prior clinical trials, DNA-based approaches have been emphasized, but recent data suggest that RNA-based assessment can capture pathway differences conveniently and may be streamlined with other RNA-based genomic risk scores. Thus, we used RNA expression to study genomic instability (HRD and TP53 pathways) in context of the breast cancer immune microenvironment in three datasets (total n = 4,892), including 1,942 samples from the Carolina Breast Cancer Study, a population-based study that oversampled Black (n = 1,026) and younger women (n = 1,032). Across all studies, 36.9% of estrogen receptor (ER)-positive and 92.6% of ER-negative breast cancer had presence of at least one genomic instability signature. TP53 and HRD status were significantly associated with immune expression in both ER-positive and ER-negative breast cancer. RNA-based genomic instability signatures were associated with higher PD-L1, CD8 T-cell marker, and global and multimarker immune cell expression. Among tumors with genomic instability signatures, adaptive immune response was associated with improved recurrence-free survival regardless of ER status, highlighting genomic instability as a candidate marker for predicting immunotherapy response. Leveraging a convenient, integrated RNA-based approach, this analysis shows that genomic instability interacts with immune response, an important target in breast cancer overall and in Black women who experience higher frequency of TP53 and HR deficiency.

Significance

Despite promising advances in breast cancer immunotherapy, predictive biomarkers that are valid across diverse populations and breast cancer subtypes are needed. Genomic instability signatures can be coordinated with other RNA-based scores to define immunogenic breast cancers and may have value in stratifying immunotherapy trial participants.

RNA-Based Classification of Homologous Recombination Deficiency in Racially Diverse Patients with Breast Cancer

BACKGROUND

Aberrant expression of DNA repair pathways such as homologous recombination (HR) can lead to DNA repair imbalance, genomic instability, and altered chemotherapy response. DNA repair imbalance may predict prognosis, but variation in DNA repair in diverse cohorts of breast cancer patients is understudied.

METHODS

To identify RNA-based patterns of DNA repair expression, we performed unsupervised clustering on 51 DNA repair-related genes in the Cancer Genome Atlas Breast Cancer [TCGA BRCA (n = 1,094)] and Carolina Breast Cancer Study [CBCS (n = 1,461)]. Using published DNA-based HR deficiency (HRD) scores (high-HRD ≥ 42) from TCGA, we trained an RNA-based supervised classifier. Unsupervised and supervised HRD classifiers were evaluated in association with demographics, tumor characteristics, and clinical outcomes.

RESULTS

: Unsupervised clustering on DNA repair genes identified four clusters of breast tumors, with one group having high expression of HR genes. Approximately 39.7% of CBCS and 29.3% of TCGA breast tumors had this unsupervised high-HRD (U-HRD) profile. A supervised HRD classifier (S-HRD) trained on TCGA had 84% sensitivity and 73% specificity to detect HRD-high samples. Both U-HRD and S-HRD tumors in CBCS had higher frequency of TP53 mutant-like status (45% and 41% enrichment) and basal-like subtype (63% and 58% enrichment). S-HRD high was more common among black patients. Among chemotherapy-treated participants, recurrence was associated with S-HRD high (HR: 2.38, 95% confidence interval = 1.50-3.78).

CONCLUSIONS

HRD is associated with poor prognosis and enriched in the tumors of black women.

IMPACT

RNA-level indicators of HRD are predictive of breast cancer outcomes in diverse populations.

Development of an automated liquid biopsy assay for methylated markers in advanced breast cancer

Current molecular liquid biopsy assays to detect recurrence or monitor response to treatment require sophisticated technology, highly trained personnel, and a turnaround time of weeks. We describe the development and technical validation of an automated Liquid Biopsy for Breast Cancer Methylation (LBx-BCM) prototype, a DNA methylation detection cartridge assay that is simple to perform and quantitatively detects nine methylated markers within 4.5 h. LBx-BCM demonstrated high interassay reproducibility when analyzing exogenous methylated DNA (75-300 DNA copies) spiked into plasma (Coefficient of Variation, CV = 7.1 - 10.9%) and serum (CV = 19.1 - 36.1%). It also demonstrated high interuser reproducibility (Spearman r = 0.887, P < 0.0001) when samples of metastatic breast cancer (MBC, N = 11) and normal control (N = 4) were evaluated independently by two users. Analyses of interplatform reproducibility indicated very high concordance between LBx-BCM and the reference assay, cMethDNA, among 66 paired plasma samples (MBC N = 40, controls N = 26; Spearman r = 0.891; 95% CI = 0.825 - 0.933, P< 0.0001). LBx-BCM achieved a ROC AUC = 0.909 (95% CI = 0.836 - 0.982), 83% sensitivity and 92% specificity; cMethDNA achieved a ROC AUC = 0.896 (95% CI = 0.817 - 0.974), 83% sensitivity and 92% specificity in test set samples. The automated LBx-BCM cartridge prototype is fast, with performance levels equivalent to the highly sensitive, manual cMethDNA method. Future prospective clinical studies will evaluate LBx-BCM detection sensitivity and its ability to monitor therapeutic response during treatment for advanced breast cancer.

Peripheral Blood Monocyte Abundance Predicts Outcomes in Patients with Breast Cancer

Biomarkers of response are needed in breast cancer to stratify patients to appropriate therapies and avoid unnecessary toxicity. We used peripheral blood gene expression and cell type abundance to identify biomarkers of response and recurrence in neoadjuvant chemotherapy treated breast cancer patients. We identified a signature of interferon and complement response that was higher in the blood of patients with pathologic complete response. This signature was preferentially expressed by monocytes in single cell RNA sequencing. Monocytes are routinely measured clinically, enabling examination of clinically measured monocytes in multiple independent cohorts. We found that peripheral monocytes were higher in patients with good outcomes in four cohorts of breast cancer patients. Blood gene expression and cell type abundance biomarkers may be useful for prognostication in breast cancer.

Significance

Biomarkers are needed in breast cancer to identify patients at risk for recurrence. Blood is an attractive site for biomarker identification due to the relative ease of longitudinal sampling. Our study suggests that blood-based gene expression and cell type abundance biomarkers may have clinical utility in breast cancer.

Tumor suppressor PLK2 may serve as a biomarker in triple-negative breast cancer for improved response to PLK1 therapeutics

Polo-like kinase (PLK) family members play important roles in cell cycle regulation. The founding member PLK1 is oncogenic and preclinically validated as a cancer therapeutic target. Paradoxically, frequent loss of chromosome 5q11-35 which includes PLK2 is observed in basal-like breast cancer. In this study, we found that PLK2 was tumor suppressive in breast cancer, preferentially in basal-like and triple-negative breast cancer (TNBC) subtypes. Knockdown of PLK1 rescued phenotypes induced by PLK2-loss both in vitro and in vivo. We also demonstrated that PLK2 directly interacted with PLK1 at prometaphase through the kinase but not the polo-box domains of PLK2, suggesting PLK2 functioned at least partially through the interaction with PLK1. Furthermore, an improved treatment response was seen in both Plk2-deleted/low mouse preclinical and PDX TNBC models using the PLK1 inhibitor volasertib alone or in combination with carboplatin. Re-expression of PLK2 in an inducible PLK2-null mouse model reduced the therapeutic efficacy of volasertib. In summary, this study delineates the effects of chromosome 5q loss in TNBC that includes PLK2, the relationship between PLK2 and PLK1, and how this may render PLK2-deleted/low tumors more sensitive to PLK1 inhibition in combination with chemotherapy.

Systems-level Analysis Reveals Multiple Modulators of Epithelial-mesenchymal Transition and Identifies DNAJB4 and CD81 as Novel Metastasis Inducers in Breast Cancer

Epithelial-mesenchymal transition (EMT) is driven by complex signaling events that induce dramatic biochemical and morphological changes whereby epithelial cells are converted into cancer cells. However, the underlying molecular mechanisms remain elusive. Here, we used mass spectrometry based quantitative proteomics approach to systematically analyze the post-translational biochemical changes that drive differentiation of human mammary epithelial (HMLE) cells into mesenchymal. We identified 314 proteins out of more than 6,000 unique proteins and 871 phosphopeptides out of more than 7,000 unique phosphopeptides as differentially regulated. We found that phosphoproteome is more unstable and prone to changes during EMT compared with the proteome and multiple alterations at proteome level are not thoroughly represented by transcriptional data highlighting the necessity of proteome level analysis. We discovered cell state specific signaling pathways, such as Hippo, sphingolipid signaling, and unfolded protein response (UPR) by modeling the networks of regulated proteins and potential kinase-substrate groups. We identified two novel factors for EMT whose expression increased on EMT induction: DnaJ heat shock protein family (Hsp40) member B4 (DNAJB4) and cluster of differentiation 81 (CD81). Suppression of DNAJB4 or CD81 in mesenchymal breast cancer cells resulted in decreased cell migration in vitro and led to reduced primary tumor growth, extravasation, and lung metastasis in vivo Overall, we performed the global proteomic and phosphoproteomic analyses of EMT, identified and validated new mRNA and/or protein level modulators of EMT. This work also provides a unique platform and resource for future studies focusing on metastasis and drug resistance.

Insights into Impact of DNA Copy Number Alteration and Methylation on the Proteogenomic Landscape of Human Ovarian Cancer via a Multi-omics Integrative Analysis

In this work, we propose iProFun, an integrative analysis tool to screen for proteogenomic functional traits perturbed by DNA copy number alterations (CNAs) and DNA methylations. The goal is to characterize functional consequences of DNA copy number and methylation alterations in tumors and to facilitate screening for cancer drivers contributing to tumor initiation and progression. Specifically, we consider three functional molecular quantitative traits: mRNA expression levels, global protein abundances, and phosphoprotein abundances. We aim to identify those genes whose CNAs and/or DNA methylations have cis-associations with either some or all three types of molecular traits. Compared with analyzing each molecular trait separately, the joint modeling of multi-omics data enjoys several benefits: iProFun experienced enhanced power for detecting significant cis-associations shared across different omics data types, and it also achieved better accuracy in inferring cis-associations unique to certain type(s) of molecular trait(s). For example, unique associations of CNAs/methylations to global/phospho protein abundances may imply posttranslational regulations.We applied iProFun to ovarian high-grade serous carcinoma tumor data from The Cancer Genome Atlas and Clinical Proteomic Tumor Analysis Consortium and identified CNAs and methylations of 500 and 121 genes, respectively, affecting the cis-functional molecular quantitative traits of the corresponding genes. We observed substantial power gain via the joint analysis of iProFun. For example, iProFun identified 117 genes whose CNAs were associated with phosphoprotein abundances by leveraging mRNA expression levels and global protein abundances. By comparison, analyses based on phosphoprotein data alone identified none. A network analysis of these 117 genes revealed the known oncogene AKT1 as a key hub node interacting with many of the rest. In addition, iProFun identified one gene, BIN2, whose DNA methylation has cis-associations with its mRNA expression, global protein, and phosphoprotein abundances. These and other genes identified by iProFun could serve as potential drug targets for ovarian cancer.