Cancer-cell derived S100A11 promotes macrophage recruitment in ER+ breast cancer

Macrophages are pivotal in driving breast tumor development, progression, and resistance to treatment, particularly in estrogen receptor-positive (ER+) tumors, where they infiltrate the tumor microenvironment (TME) influenced by cancer cell-secreted factors. By analyzing single-cell RNA-sequencing data from 25 ER+ tumors, we elucidated interactions between cancer cells and macrophages, correlating macrophage density with epithelial cancer cell density. We identified that S100A11, a previously unexplored factor in macrophage-cancer crosstalk, predicts high macrophage density and poor outcomes in ER+ tumors. We found that recombinant S100A11 enhances macrophage infiltration and migration in a dose-dependent manner. Additionally, in 3D models, we showed that S100A11 expression levels in ER+ cancer cells predict macrophage infiltration patterns. Neutralizing S100A11 decreased macrophage recruitment, both in cancer cell lines and in a clinically relevant patient-derived organoid model, underscoring its role as a paracrine regulator of cancer-macrophage interactions in the protumorigenic TME. This study offers novel insights into the interplay between macrophages and cancer cells in ER+ breast tumors, highlighting S100A11 as a potential therapeutic target to modulate the macrophage-rich tumor microenvironment.

Estrogen regulation and functional role of FGFR4 in estrogen receptor positive breast cancer

Background

Resistance to endocrine therapy is a major challenge of managing estrogen receptor positive (ER+) breast cancer. We previously reported frequent overexpression of FGFR4 in endocrine resistant cell lines and breast cancers that recurred and metastasized following endocrine therapy, suggesting FGFR4 as a potential driver of endocrine resistance. In this study, we investigated the role of FGFR4 in mediating endocrine resistance and explored the therapeutic potential of targeting FGFR4 in advanced breast cancer.

Methods

A gene expression signature of FGFR4 activity was examined in ER+ breast cancer pre- and post-neoadjuvant endocrine therapy and the association between FGFR4 expression and patient survival was examined. A correlation analysis was used to uncover potential regulators of FGFR4 overexpression. To investigate if FGFR4 is necessary to drive endocrine resistance, we tested response to FGFR4 inhibition in long term estrogen deprived (LTED) cells and their paired parental cells. Doxycycline inducible FGFR4 overexpression and knockdown cell models were generated to examine if FGFR4 was sufficient to confer endocrine resistance. Finally, we examined response to FGFR4 monotherapy or combination therapy with fulvestrant in breast cancer cell lines to explore the potential of FGFR4 targeted therapy for advanced breast cancer and assessed the importance of PAM50 subtype in response to FGFR4 inhibition.

Results

A FGFR4 activity gene signature was significantly upregulated post neoadjuvant aromatase inhibitor treatment, and high FGFR4 expression predicted poorer survival in patients with ER+ breast cancer. Gene expression association analysis using TCGA, METABRIC and SCAN-B datasets uncovered ER as the most significant gene negatively correlated with FGFR4 expression. ER negatively regulates FGFR4 expression at both the mRNA and protein level across multiple ER+ breast cancer cell lines. Despite robust overexpression of FGFR4, LTED cells did not show enhanced responses to FGFR4 inhibition compared to parental cells. Similarly, FGFR4 overexpression, knockdown or hotspot mutations did not significantly alter response to endocrine treatment in ER+ cell lines, nor did FGFR4 and fulvestrant combination treatment show synergistic effects. The HER2-like subtype of breast cancer showed elevated expression of FGFR4 and an increased response to FGFR4 inhibition relative to other breast cancer subtypes.

Conclusions

Despite ER-mediated upregulation of FGFR4 post endocrine therapy, our study does not support a general role of FGFR4 in mediating endocrine resistance in ER+ breast cancer. Our data suggests that specific genomic backgrounds such as HER2 expression may be required for FGFR4 function in breast cancer and should be further explored.

Liver tropism of ER mutant breast cancer is characterized by unique molecular changes and immune infiltration

PURPOSE

Hotspot estrogen receptor alpha (ER/ESR1) mutations are recognized as the driver for both endocrine resistance and metastasis in advanced ER-positive (ER+) breast cancer, but their contributions to metastatic organ tropism remain insufficiently understood. In this study, we aim to comprehensively profile the organotropic metastatic pattern for ESR1 mutant breast cancer.

METHODS

The organ-specific metastatic pattern of ESR1 mutant breast cancer was delineated using multi-omics data from multiple publicly available cohorts of ER+ metastatic breast cancer patients. Gene mutation/copy number variation (CNV) and differential gene expression analyses were performed to identify the genomic and transcriptomic alterations uniquely associated with ESR1 mutant liver metastasis. Upstream regulator, downstream pathway, and immune infiltration analysis were conducted for subsequent mechanistic investigations.

RESULTS

ESR1 mutation-driven liver tropism was revealed by significant differences, encompassing a higher prevalence of liver metastasis in patients with ESR1 mutant breast cancer and an enrichment of mutations in liver metastatic samples. The significant enrichment of AGO2 copy number amplifications (CNAs) and multiple gene expression changes were revealed uniquely in ESR1 mutant liver metastasis. We also unveiled alterations in downstream signaling pathways and immune infiltration, particularly an enrichment of neutrophils, suggesting potential therapeutic vulnerabilities.

CONCLUSION

Our data provide a comprehensive characterization of the behaviors and mechanisms of ESR1 mutant liver metastasis, paving the way for the development of personalized therapy to target liver metastasis for patients with ESR1 mutant breast cancer.

Use of natural language understanding to facilitate surgical de-escalation of axillary staging in patients with breast cancer

Natural language understanding (NLU) may be particularly well-equipped for enhanced data capture from the electronic health record (EHR) given its examination of both content- and context-driven extraction. We developed and applied a NLU model to examine rates of pathological node positivity (pN+) and rates of lymphedema to determine if omission of routine axillary staging could be extended to younger patients with ER+/cN0 disease. We found that rates of pN+ and arm lymphedema were similar between patients 55-69yo and ≥70yo, with rates of lymphedema exceeding rates of pN+ for clinical stage T1c and smaller disease. Data from our NLU model suggest that omission of SLNB might be extended beyond Choosing Wisely recommendations, limited to those over 70 years old, to all postmenopausal women with early-stage ER+/cN0 disease. These data support the recently-reported SOUND trial results and provide additional granularity to facilitate surgical de-escalation.

Systems approach for congruence and selection of cancer models towards precision medicine

Cancer models are instrumental as a substitute for human studies and to expedite basic, translational, and clinical cancer research. For a given cancer type, a wide selection of models, such as cell lines, patient-derived xenografts, organoids and genetically modified murine models, are often available to researchers. However, how to quantify their congruence to human tumors and to select the most appropriate cancer model is a largely unsolved issue. Here, we present Congruence Analysis and Selection of CAncer Models (CASCAM), a statistical and machine learning framework for authenticating and selecting the most representative cancer models in a pathway-specific manner using transcriptomic data. CASCAM provides harmonization between human tumor and cancer model omics data, systematic congruence quantification, and pathway-based topological visualization to determine the most appropriate cancer model selection. The systems approach is presented using invasive lobular breast carcinoma (ILC) subtype and suggesting CAMA1 followed by UACC3133 as the most representative cell lines for ILC research. Two additional case studies for triple negative breast cancer (TNBC) and patient-derived xenograft/organoid (PDX/PDO) are further investigated. CASCAM is generalizable to any cancer subtype and will authenticate cancer models for faithful non-human preclinical research towards precision medicine.

MRTF promotes breast cancer cell motility through SRF-dependent upregulation of DIAPH3 expression

Dysregulated actin cytoskeleton gives rise to aberrant cell motility and metastatic spread of tumor cells. The MRTF-SRF transcriptional complex plays a key role in regulating the expressions of actin cytoskeleton-modulatory genes. In this study, we demonstrate that MRTF's interaction with SRF is critical for migration and invasion of breast cancer cells. Disruption of the MRTF-SRF interaction suppresses membrane dynamics affecting the frequency and the effectiveness of membrane protrusion during cell motility. Consistent with these phenotypic changes, we further show that MRTF promotes actin polymerization at the leading edge, a key aspect of membrane protrusion, and migration of breast cancer cells through upregulating the expression of formin-family actin nucleating/elongating protein encoding gene DIAPH3 in an SRF-dependent manner. In support of these findings, multiplexed quantitative immunohistochemistry and transcriptome analyses of clinical specimens of breast cancer further demonstrate a positive correlation between nuclear localization of MRTF with malignant traits of cancer cells as well as enrichment of MRTF/SRF gene signature in distant metastases relative to primary tumors. In conclusion, this study for the first time links the MRTF/SRF signaling axis to cell migration through the regulation of a specific actin-binding protein, and provides evidence for an association between MRTF/SRF activity and malignancy in human breast cancer.

Multi-omic characterization of ILC and ILC-like cell lines as part of ILC cell line encyclopedia (ICLE) defines new models to study potential biomarkers and explore therapeutic opportunities

Invasive lobular carcinoma (ILC), the most common histological "special type", accounts for ∼10-15% of all BC diagnoses, is characterized by unique features such as E-cadherin loss/deficiency, lower grade, hormone receptor positivity, larger diffuse tumors, and specific metastatic patterns. Despite ILC being acknowledged as a disease with distinct biology that necessitates specialized and precision medicine treatments, the further exploration of its molecular alterations with the goal of discovering new treatments has been hindered due to the scarcity of well-characterized cell line models for studying this disease. To address this, we generated the ILC Cell Line Encyclopedia (ICLE), providing a comprehensive multi-omic characterization of ILC and ILC-like cell lines. Using consensus multi-omic subtyping, we confirmed luminal status of previously established ILC cell lines and uncovered additional ILC/ILC-like cell lines with luminal features for modeling ILC disease. Furthermore, most of these luminal ILC/ILC-like cell lines also showed RNA and copy number similarity to ILC patient tumors. Similarly, ILC/ILC-like cell lines also retained molecular alterations in key ILC genes at similar frequency to both primary and metastatic ILC tumors. Importantly, ILC/ILC-like cell lines recapitulated the CDH1 alteration landscape of ILC patient tumors including enrichment of truncating mutations in and biallelic inactivation of CDH1 gene. Using whole-genome optical mapping, we uncovered novel genomic-rearrangements including novel structural variations in CDH1 and functional gene fusions and characterized breast cancer specific patterns of chromothripsis in chromosomes 8, 11 and 17. In addition, we systematically analyzed aberrant DNAm events and integrative analysis with RNA expression revealed epigenetic activation of TFAP2B - an emerging biomarker of lobular disease that is preferentially expressed in lobular disease. Finally, towards the goal of identifying novel druggable vulnerabilities in ILC, we analyzed publicly available RNAi loss of function breast cancer cell line datasets and revealed numerous putative vulnerabilities cytoskeletal components, focal adhesion and PI3K/AKT pathway in ILC/ILC-like vs NST cell lines. In summary, we addressed the lack of suitable models to study E-cadherin deficient breast cancers by first collecting both established and putative ILC models, then characterizing them comprehensively to show their molecular similarity to patient tumors along with uncovering their novel multi-omic features as well as highlighting putative novel druggable vulnerabilities. Not only we expand the array of suitable E-cadherin deficient cell lines available for modelling human-ILC disease but also employ them for studying epigenetic activation of a putative lobular biomarker as well as identifying potential druggable vulnerabilities for this disease towards enabling precision medicine research for human-ILC.

MRTF activity in breast cancer cells promotes osteoclastogenesis through a paracrine action of CTGF

Bone is a frequent site for breast cancer metastasis. Conditioning of local tumor microenvironment through crosstalk between tumor cells and bone resident cells in the metastatic niche is a major driving force for bone colonization of cancer cells. This study demonstrates that Myocardin-related transcription factor (MRTF - a major cofactor for the transcription factor serum-response factor, SRF) activity in breast cancer cells is required for their ability to promote osteoclast differentiation of bone marrow-derived monocytes and colonize in bone. MRTF depletion in breast cancer cells affects a wide range of cell-secreted osteoclast-regulatory factors including connective tissue growth factor (CTGF), a prominent bone metastasis-associated gene that exhibit strong positive association in expression with MRTF activity in human breast cancer. Rescue experiments demonstrate that CTGF is an important paracrine mediator of pro-osteoclastogenic action of MRTF in breast cancer cells. Both SRF-dependent and -independent (SAP-domain directed) functions of MRTF are required for its ability to regulate CTGF expression and osteoclast differentiation. In conclusion, this study uncovers a novel MRTF-directed tumor-extrinsic mechanism of bone colonization of cancer cells and suggest that MRTF inhibition could be a novel strategy to suppress osteoclast activity and skeletal involvement in metastatic breast cancer.

iGenSig-Rx: an integral genomic signature based white-box tool for modeling cancer therapeutic responses using multi-omics data

Multi-omics sequencing is expected to become clinically routine within the next decade and transform clinical care. However, there is a paucity of viable and interpretable genome-wide modeling methods that can facilitate rational selection of patients for tailored intervention. Here we develop an integral genomic signature-based method called iGenSig-Rx as a white-box tool for modeling therapeutic response based on clinical trial datasets with improved cross-dataset applicability and tolerance to sequencing bias. This method leverages high-dimensional redundant genomic features to address the challenges of cross-dataset modeling, a concept similar to the use of redundant steel rods to reinforce the pillars of a building. Using genomic datasets for HER2 targeted therapies, the iGenSig-Rx model demonstrates stable predictive power across four independent clinical trials. More importantly, the iGenSig-Rx model offers the level of transparency much needed for clinical application, allowing for clear explanations as to how the predictions are produced, how the features contribute to the prediction, and what are the key underlying pathways. We expect that iGenSig-Rx as a class of biologically interpretable multi-omics modeling methods will have broad applications in big-data based precision oncology. The R package is available: https://github.com/wangxlab/iGenSig-Rx. NOTE: the Github website will be released upon publication and the R package is available for review through google drive: https://drive.google.com/drive/folders/1KgecmUoon9-h2Dg1rPCyEGFPOp28Ols3?usp=sharing.

Effects of socioeconomic status and race on survival and treatment in metastatic breast cancer

Race and socioeconomic factors affect outcomes in breast cancer. We aimed to assess the effect of race and neighborhood socioeconomic status (SES) on overall survival and treatment patterns in patients with metastatic breast cancer (MBC). This is a retrospective cohort study involving patients (N = 1246) with distant breast cancer metastases diagnosed at UPMC Magee Women's Breast Cancer Clinic from 2000-2017. Overall survival and treatment patterns were compared between races (Blacks and whites) and SES groups (defined using Area Deprivation Index). Low SES, but not tumor characteristics, was associated with Black race (P < 0.001) in the study population. Low SES (Median [Interquartile Range, IQR] survival 2.3[2.2-2.5] years vs high SES 2.7[2.5-3.1] years, P = 0.01) and Black race (Median [IQR] survival 1.8[1.3-2.3] years, vs white 2.5[2.3-2.7] years P = 0.008) separately were associated with worse overall survival in patients with MBC. In the Cox Proportional Hazard model with SES, race, age, subtype, number of metastases, visceral metastasis, and year of diagnosis as covariates, low SES (Hazard ratio 1.19[1.04-1.35], P = 0.01), but not Black race (Hazard ratio 1.19[0.96-1.49], P = 0.12), independently predicted overall survival in MBC. Moreover, patients from low SES neighborhoods and Black race received fewer lines of chemotherapy than high SES and whites. In conclusion, low neighborhood SES is associated with worse outcomes in patients with MBC. Poor outcomes in Black patients with MBC, at least in part is driven by socioeconomic factors. Future studies should delineate the interplay between neighborhood SES, race, and their effects on tumor biology in MBC.

Small molecule nitroalkenes inhibit RAD51-mediated homologous recombination and amplify triple-negative breast cancer cell killing by DNA-directed therapies

Nitro fatty acids (NO2-FAs) are endogenously generated lipid signaling mediators from metabolic and inflammatory reactions between conjugated diene fatty acids and nitric oxide or nitrite-derived reactive species. NO2-FAs undergo reversible Michael addition with hyperreactive protein cysteine thiolates to induce posttranslational protein modifications that can impact protein function. Herein, we report a novel mechanism of action of natural and non-natural nitroalkenes structurally similar to (E) 10-nitro-octadec-9-enoic acid (CP-6), recently de-risked by preclinical Investigational New Drug-enabling studies and Phase 1 and Phase 2 clinical trials and found to induce DNA damage in a TNBC xenograft by inhibiting homologous-recombination (HR)-mediated repair of DNA double-strand breaks (DSB). CP-6 specifically targets Cys319, essential in RAD51-controlled HR-mediated DNA DSB repair in cells. A nitroalkene library screen identified two structurally different nitroalkenes, a non-natural fatty acid [(E) 8-nitro-nonadec-7-enoic acid (CP-8)] and a dicarboxylate ester [dimethyl (E)nitro-oct-4-enedioate (CP-23)] superior to CP-6 in TNBC cells killing, synergism with three different inhibitors of the poly ADP-ribose polymerase (PARP) and γ-IR. CP-8 and CP-23 effectively inhibited γ-IR-induced RAD51 foci formation and HR in a GFP-reported assay but did not affect benign human epithelial cells or cell cycle phases. In vivo, CP-8 and CP-23's efficacies diverged as only CP-8 showed promising anticancer activities alone and combined with the PARP inhibitor talazoparib in an HR-proficient TNBC mouse model. As preliminary preclinical toxicology analysis also suggests CP-8 as safe, our data endorse CP-8 as a novel anticancer molecule for treating cancers sensitive to homologous recombination-mediated DNA repair inhibitors.

ESR1 fusion proteins invoke breast cancer subtype-dependent enrichment of ligand independent pro-oncogenic signatures and phenotypes

Breast cancer is a leading cause of female mortality and despite advancements in diagnostics and personalized therapeutics, metastatic disease largely remains incurable due to drug resistance. Fortunately, identification of mechanisms of therapeutic resistance have rapidly transformed our understanding of cancer evasion and is enabling targeted treatment regimens. When the druggable estrogen receptor (ER, ESR1 ), expressed in two-thirds of all breast cancer, is exposed to endocrine therapy, there is risk of somatic mutation development in approximately 30% of cases and subsequent treatment resistance. A more recently discovered mechanism of ER mediated endocrine resistance is the expression of ER fusion proteins. ER fusions, which retain the protein's DNA binding domain, harbor ESR1 exons 1-6 fused to an in-frame gene partner resulting in loss of the 3' ER ligand binding domain (LBD). In this report we demonstrate that in no-special type (NST) and invasive lobular carcinoma (ILC) cell line models, ER fusion proteins exhibit robust hyperactivation of canonical ER signaling pathways independent of the ligand estradiol or anti-endocrine therapies such as Fulvestrant and Tamoxifen. We employ cell line models stably overexpressing ER fusion proteins with concurrent endogenous ER knockdown to minimize the influence of endogenous wildtype ER. Cell lines exhibited shared transcriptomic enrichment in pathways known to be drivers of metastatic disease, notably the MYC pathway. The heterogeneous 3' fusion partners, particularly transcription factors SOX9 and YAP1 , evoked varying degrees of transcriptomic and cistromic activity that translated into unique phenotypic readouts. Herein we report that cell line activity is subtype-, fusion-, and assay-specific suggesting that the loss of the LBD, the 3' fusion partner, and the cellular landscape all influence fusion activity. Therefore, it will be critical to generate additional data on frequency of the ER fusions, in the context of the clinicopathological features of the tumor.

Significance

ER fusion proteins exhibit diverse mechanisms of endocrine resistance in breast cancer cell lines representing the no special type (NST) and invasive lobular cancer (ILC) subtypes. Our emphasize upon both the shared and unique cellular adaptations imparted by ER fusions offers the foundation for further translational research and clinical decision making.

WCRC-25: A novel luminal Invasive Lobular Carcinoma cell line model

Breast cancer is categorized by the molecular and histologic presentation of the tumor, with the major histologic subtypes being No Special Type (NST) and Invasive Lobular Carcinoma (ILC). ILC are characterized by growth in a single file discohesive manner with stromal infiltration attributed to their hallmark pathognomonic loss of E-cadherin ( CDH1 ). Few ILC cell line models are available to researchers. Here we report the successful establishment and characterization of a novel ILC cell line, WCRC-25, from a metastatic pleural effusion from a postmenopausal Caucasian woman with metastatic ILC. WCRC-25 is an ER-negative luminal epithelial ILC cell line with both luminal and Her2-like features. It exhibits anchorage independent growth and haptotactic migration towards Collagen I. Sequencing revealed a CDH1 Q706* truncating mutation, together with mutations in FOXA1, CTCF, BRCA2 and TP53 , which were also seen in a series of metastatic lesions from the patient. Copy number analyses revealed amplification and deletion of genes frequently altered in ILC while optical genome mapping revealed novel structural rearrangements. RNA-seq analysis comparing the primary tumor, metastases and the cell line revealed signatures for cell cycle progression and receptor tyrosine kinase signaling. To assess targetability, we treated WCRC-25 with AZD5363 and Alpelisib confirming WCRC-25 as susceptible to PI3K/AKT signaling inhibition as predicted by our RNA sequencing analysis. In conclusion, we report WCRC-25 as a novel ILC cell line with promise as a valuable research tool to advance our understanding of ILC and its therapeutic vulnerabilities.

Financial support

The work was in part supported by a Susan G Komen Leadership Grant to SO (SAC160073) and NCI R01 CA252378 (SO/AVL). AVL and SO are Komen Scholars, Hillman Foundation Fellows and supported by BCRF. This project used the UPMC Hillman Cancer Center and Tissue and Research Pathology/Pitt Biospecimen Core shared resource which is supported in part by award P30CA047904. This research was also supported in part by the University of Pittsburgh Center for Research Computing, RRID:SCR_022735, through the resources provided. Specifically, this work used the HTC cluster, which is supported by NIH award number S10OD028483. Finally, partial support was provided by the Magee-Womens Research Institute and Foundation, The Shear Family Foundation, and The Metastatic Breast Cancer Network.

The EstroGene Database Reveals Diverse Temporal, Context-Dependent, and Bidirectional Estrogen Receptor Regulomes in Breast Cancer

As one of the most successful cancer therapeutic targets, estrogen receptor-α (ER/ESR1) has been extensively studied over the past few decades. Sequencing technological advances have enabled genome-wide analysis of ER action. However, comparison of individual studies is limited by different experimental designs, and few meta-analyses are available. Here, we established the EstroGene database through unified processing of data from 246 experiments including 136 transcriptomic, cistromic, and epigenetic datasets focusing on estradiol (E2)-triggered ER activation across 19 breast cancer cell lines. A user-friendly browser (https://estrogene.org/) was generated for multiomic data visualization involving gene inquiry under user-defined experimental conditions and statistical thresholds. Notably, annotation of metadata associated with public datasets revealed a considerable lack of experimental details. Comparison of independent RNA-seq or ER ChIP-seq data with the same design showed large variability and only strong effects could be consistently detected. Temporal estrogen response metasignatures were defined, and the association of E2 response rate with temporal transcriptional factors, chromatin accessibility, and heterogeneity of ER expression was evaluated. Unexpectedly, harmonizing 146 E2-induced transcriptomic datasets uncovered a subset of genes harboring bidirectional E2 regulation, which was linked to unique transcriptional factors and highly associated with immune surveillance in the clinical setting. Furthermore, the context dependent E2 response programs were characterized in MCF7 and T47D cell lines, the two most frequently used models in the EstroGene database. Collectively, the EstroGene database provides an informative and practical resource to the cancer research community to uniformly evaluate key reproducible features of ER regulomes and unravels modes of ER signaling.

SIGNIFICANCE

A resource database integrating 246 publicly available ER profiling datasets facilitates meta-analyses and identifies estrogen response temporal signatures, a bidirectional program, and model-specific biases.

Small molecule nitroalkenes inhibit RAD51-mediated homologous recombination and amplify triple-negative breast cancer cell killing by DNA-directed therapies

Nitro fatty acids (NO 2 -FAs) are endogenously generated lipid signaling mediators from metabolic and inflammatory reactions between conjugated diene fatty acids and nitric oxide or nitrite-derived reactive species. NO 2 -FAs undergo reversible Michael addition with hyperreactive protein cysteine thiolates to induce posttranslational protein modifications that can impact protein function. Herein, we report a novel mechanism of action of natural and non-natural nitroalkenes structurally similar to ( E ) 10-nitro-octadec-9-enoic acid (CP-6), recently de-risked by preclinical Investigational New Drug-enabling studies and Phase 1 and Phase 2 clinical trials and found to induce DNA damage in a TNBC xenograft by inhibiting homologous-recombination (HR)-mediated repair of DNA double-strand breaks (DSB). CP-6 specifically targets Cys319, essential in RAD51-controlled HR-mediated DNA DSB repair in cells. A nitroalkene library screen identified two structurally different nitroalkenes, a non-natural fatty acid [( E ) 8-nitro- nonadec-7-enoic acid (CP-8)] and a dicarboxylate ester [dimethyl ( E )nitro-oct-4-enedioate (CP- 23)] superior to CP-6 in TNBC cells killing, synergism with three different inhibitors of the poly ADP-ribose polymerase (PARP) and γ-IR. CP-8 and CP-23 effectively inhibited γ-IR-induced RAD51 foci formation and HR in a GFP-reported assay but did not affect benign human epithelial cells or cell cycle phases. In vivo, CP-8 and CP-23's efficacies diverged as only CP-8 showed promising anticancer activities alone and combined with the PARP inhibitor talazoparib in an HR-proficient TNBC mouse model. As preliminary preclinical toxicology analysis also suggests CP-8 as safe, our data endorse CP-8 as a novel anticancer molecule for treating cancers sensitive to homologous recombination-mediated DNA repair inhibitors.

Clinicopathologic and genomic features of lobular like invasive mammary carcinoma: is it a distinct entity?

This study describes "lobular-like invasive mammary carcinomas" (LLIMCas), a group of low- to intermediate-grade invasive mammary carcinomas with discohesive, diffusely infiltrative cells showing retained circumferential membranous immunoreactivity for both E-cadherin and p120. We analyzed the clinical-pathologic features of 166 LLIMCas compared to 104 classical invasive lobular carcinomas (ILCs) and 100 grade 1 and 2 invasive ductal carcinomas (IDCs). Tumor size and pT stage of LLIMCas were intermediate between IDCs and ILCs, and yet often underestimated on imaging and showed frequent positive margins on the first resection. Despite histomorphologic similarities to classical ILC, the discohesion in LLIMCa was independent of E-cadherin/p120 immunophenotypic alteration. An exploratory, hypothesis-generating analysis of the genomic features of 14 randomly selected LLIMCas and classical ILCs (7 from each category) was performed utilizing an FDA-authorized targeted capture sequencing assay (MSK-IMPACT). None of the seven LLIMCas harbored CDH1 loss-of-function mutations, and none of the CDH1 alterations detected in two of the LLIMCas was pathogenic. In contrast, all seven ILCs harbored CDH1 loss-of-function mutations coupled with the loss of heterozygosity of the CDH1 wild-type allele. Four of the six evaluable LLIMCas were positive for CDH1 promoter methylation, which may partially explain the single-cell infiltrative morphology seen in LLIMCa. Further studies are warranted to better define the molecular basis of the discohesive cellular morphology in LLIMCa. Until more data becomes available, identifying LLIMCas and distinguishing them from typical IDCs and ILCs would be justified. In patients with LLIMCas, preoperative MRI should be entertained to guide surgical management.

p38MAPKα Stromal Reprogramming Sensitizes Metastatic Breast Cancer to Immunotherapy

Metastatic breast cancer is an intractable disease that responds poorly to immunotherapy. We show that p38MAPKα inhibition (p38i) limits tumor growth by reprogramming the metastatic tumor microenvironment in a CD4+ T cell-, IFNγ-, and macrophage-dependent manner. To identify targets that further increased p38i efficacy, we utilized a stromal labeling approach and single-cell RNA sequencing. Thus, we combined p38i and an OX40 agonist that synergistically reduced metastatic growth and increased overall survival. Intriguingly, patients with a p38i metastatic stromal signature had better overall survival that was further improved by the presence of an increased mutational load, leading us to ask if our approach would be effective in antigenic breast cancer. The combination of p38i, anti-OX40, and cytotoxic T-cell engagement cured mice of metastatic disease and produced long-term immunologic memory. Our findings demonstrate that a detailed understanding of the stromal compartment can be used to design effective antimetastatic therapies.

SIGNIFICANCE

Immunotherapy is rarely effective in breast cancer. We dissected the metastatic tumor stroma, which revealed a novel therapeutic approach that targets the stromal p38MAPK pathway and creates an opportunity to unleash an immunologic response. Our work underscores the importance of understanding the tumor stromal compartment in therapeutic design. This article is highlighted in the In This Issue feature, p. 1275.

FOXA1 Reprogramming Dictates Retinoid X Receptor Response in ESR1-Mutant Breast Cancer

Estrogen receptor alpha (ER/ESR1) mutations occur in 30% to 40% of endocrine resistant ER-positive (ER+) breast cancer. Forkhead box A1 (FOXA1) is a key pioneer factor mediating ER-chromatin interactions and endocrine response in ER+ breast cancer, but its role in ESR1-mutant breast cancer remains unclear. Our previous FOXA1 chromatin immunoprecipitation sequencing (ChIP-seq) identified a large portion of redistributed binding sites in T47D genome-edited Y537S and D538G ESR1-mutant cells. Here, we further integrated FOXA1 genomic binding profile with the isogenic ER cistrome, accessible genome, and transcriptome data of T47D cell model. FOXA1 redistribution was significantly associated with transcriptomic alterations caused by ESR1 mutations. Furthermore, in ESR1-mutant cells, FOXA1-binding sites less frequently overlapped with ER, and differential gene expression was less associated with the canonical FOXA1-ER axis. Motif analysis revealed a unique enrichment of retinoid X receptor (RXR) motifs in FOXA1-binding sites of ESR1-mutant cells. Consistently, ESR1-mutant cells were more sensitive to growth stimulation with the RXR agonist LG268. The mutant-specific response was dependent on two RXR isoforms, RXR-α and RXR-β, with a stronger dependency on the latter. In addition, T3, the agonist of thyroid receptor (TR) also showed a similar growth-promoting effect in ESR1-mutant cells. Importantly, RXR antagonist HX531 blocked growth of ESR1-mutant cells and a patient-derived xenograft (PDX)-derived organoid with an ESR1 D538G mutation. Collectively, our data support the evidence for a stronger RXR response associated with FOXA1 reprograming in ESR1-mutant cells, suggesting development of therapeutic strategies targeting RXR pathways in breast tumors with ESR1 mutation.

IMPLICATIONS

It provides comprehensive characterization of the role of FOXA1 in ESR1-mutant breast cancer and potential therapeutic strategy through blocking RXR activation.

Immune landscape in invasive ductal and lobular breast cancer reveals a divergent macrophage-driven microenvironment

T cell-centric immunotherapies have shown modest clinical benefit thus far for estrogen receptor-positive (ER+) breast cancer. Despite accounting for 70% of all breast cancers, relatively little is known about the immunobiology of ER+ breast cancer in women with invasive ductal carcinoma (IDC) and invasive lobular carcinoma (ILC). To investigate this, we performed phenotypic, transcriptional and functional analyses for a cohort of treatment-naive IDC (n = 94) and ILC (n = 87) tumors. We show that macrophages, and not T cells, are the predominant immune cells infiltrating the tumor bed and the most transcriptionally diverse cell subset between IDC and ILC. Analysis of cellular neighborhoods revealed an interplay between macrophages and T cells associated with longer disease-free survival in IDC but not ILC. Our datasets provide a rich resource for further interrogation into immune cell dynamics in ER+ IDC and ILC and highlight macrophages as a potential target for ER+ breast cancer.

Mixed invasive ductal lobular carcinoma is clinically and pathologically more similar to invasive lobular than ductal carcinoma

BACKGROUND

Mixed invasive ductal lobular carcinoma (mDLC) remains a poorly understood subtype of breast cancer composed of coexisting ductal and lobular components.

METHODS

We sought to describe clinicopathologic characteristics and determine whether mDLC is clinically more similar to invasive ductal carcinoma (IDC) or invasive lobular carcinoma (ILC), using data from patients seen at the University of Pittsburgh Medical Center.

RESULTS

We observed a higher concordance in clinicopathologic characteristics between mDLC and ILC, compared to IDC. There is a trend for higher rates of successful breast-conserving surgery after neoadjuvant chemotherapy in patients with mDLC compared to patients with ILC, in which it is known to be lower than in those with IDC. Metastatic patterns of mDLC demonstrate a propensity to develop in sites characteristic of both IDC and ILC. A meta-analysis evaluating mDLC showed shared features with both ILC and IDC with significantly more ER-positive and fewer high grades in mDLC compared to IDC, although mDLCs were significantly smaller and included fewer late-stage tumours compared to ILC.

CONCLUSIONS

These findings support clinicopathologic characteristics of mDLC driven by individual ductal vs lobular components and given the dominance of lobular pathology, mDLC features are often more similar to ILC than IDC. This study exemplifies the complexity of mixed disease.

RET signaling in breast cancer therapeutic resistance and metastasis

RET, a single-pass receptor tyrosine kinase encoded on human chromosome 10, is well known to the field of developmental biology for its role in the ontogenesis of the central and enteric nervous systems and the kidney. In adults, RET alterations have been characterized as drivers of non-small cell lung cancer and multiple neuroendocrine neoplasms. In breast cancer, RET signaling networks have been shown to influence diverse functions including tumor development, metastasis, and therapeutic resistance. While RET is known to drive the development and progression of multiple solid tumors, therapeutic agents selectively targeting RET are relatively new, though multiple multi-kinase inhibitors have shown promise as RET inhibitors in the past; further, RET has been historically neglected as a potential therapeutic co-target in endocrine-refractory breast cancers despite mounting evidence for a key pathologic role and repeated description of a bi-directional relationship with the estrogen receptor, the principal driver of most breast tumors. Additionally, the recent discovery of RET enrichment in breast cancer brain metastases suggests a role for RET inhibition specific to advanced disease. This review assesses the status of research on RET in breast cancer and evaluates the therapeutic potential of RET-selective kinase inhibitors across major breast cancer subtypes.

Abstract P6-14-02: Hope for OTHERS – An organ donation program for metastatic breast cancer research

Background: Previous studies have shown that rapid autopsies (RA) provide a unique opportunity for tissue collection from patients who succumb to the disease. Because cancer patients are unable to donate their organs to other people, this program provides the patient an opportunity to leave a legacy by donating their body to research. These donations are vital for advancing breast cancer research. The UPMC/Pitt RA group revamped an existing program in 2018 through the formation of a larger multidisciplinary team that includes breast cancer laboratory and clinical researchers, pathologists, nurses, bioinformaticians, and tissue bankers. Because recruitment to the RA program was a challenge, we recently added patient advocates to the team to provide their essential perspective, and a dedicated research coordinator who serves as an ambassador for the program. Methods: Autopsy is performed by the Autopsy and Forensic Pathology Center of Excellence/Decedent Affairs Service of UPMC. Samples are banked in the Pitt Biospecimen Core (PBC), in addition to immediate processing including preparing of samples for sequencing and growing of organoids in the laboratory. Immunohistochemical (IHC) analysis is performed by UPMC/Magee Pathology. Results: The research coordinator quickly became an integral part of the program and closely interacts with care providers, patients and their families, pathologists on call, and manages interactions with transport services. Five breast cancer advocates have been instrumental in advising on additional changes to the program. The advocates attend regular team meetings and have formulated patient considerations for the the RA program, including appropriate and sensitive recruitment of patients, the role of physicians in decision making by the patient, registration for more than one RA program, potential issues with transporting a body across state lines and more. The advocates also developed the name for the program - “Hope for OTHERS” with Others standing for “Our Tissue Helping Enhance Research & Science”. As of June 2022, the team has completed 26 autopsies, and an additional 20 patients have consented to the program. The completed autopsies include patients with breast tumors representing different molecular and histological classes, ethnicities, and genders. The average disease-free survival and overall survival of patients that underwent autopsy was 81.6 and 127.8 months, respectively. Most patients passed outside the hospital (86%), with 62% in home hospice and 24% in inpatient hospice. Average time between death and start and end of autopsy was 4.56 hrs and 7.09 hours, respectively. The most common metastatic sites from which specimens are collected are liver, lung and lymph nodes. Per patient we collect on average specimens from 4 different organs. In addition to the metastatic lesions, we have access to primary tumor tissue and liquid biopsies obtained during the breast cancer disease progression for 44% and 73% of the patients, respectively. For a subset of the patients, tissue has been grown as patient-derived organoids or xenograft models. Preliminary IHC and sequencing analysis has provided insight into inter- and intra-patient and intra-tumor heterogeneity. Further molecular studies are ongoing. Conclusion In summary, over the last 5+ years, we have established a successful post-mortem tissue collection program, by addressing a series of barriers through the formation and work of a multi-disciplinary well-coordinated team. We are currently expanding our omics studies using state-of-the-art technologies to improve understanding how intra- and inter-tumor heterogeneity play a role in the clinical course of advanced breast cancer, to increase diversity of the patients enrolled in the RA program, and to support the successful implementation of other RA programs nationwide and worldwide.

Citation Format: Steffi Oesterreich, Lori Miller, Margaret Q. Rosenzweig, Tanner L. Bartholow, Megan Yates, Ashuvinee Elangovan, Laura Savariau, Allison N. Casey, Nolan Priedigkeit, Kai Ding, Abdalla Wedn, Jie Bin Liu, Daniel D. Brown, Tara Hyder, Geoffrey Pecar, Neil Carleton, Humberto Trejo Bittar, Daniel Geisler, Oscar Lopez-Nunez, Amanda M. Clark, Alan Wells, Partha Roy, Shannon Puhalla, Naomi Howard, Christine Needles, Susan Trent, Stephanie Walker, Christine Hodgdon, Rohit Bhargava, Jennifer M. Atkinson, Adrian V. Lee. Hope for OTHERS – An organ donation program for metastatic breast cancer research [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P6-14-02.

Transcriptomic congruence analysis for evaluating model organisms

Model organisms are instrumental substitutes for human studies to expedite basic, translational, and clinical research. Despite their indispensable role in mechanistic investigation and drug development, molecular congruence of animal models to humans has long been questioned and debated. Little effort has been made for an objective quantification and mechanistic exploration of a model organism's resemblance to humans in terms of molecular response under disease or drug treatment. We hereby propose a framework, namely Congruence Analysis for Model Organisms (CAMO), for transcriptomic response analysis by developing threshold-free differential expression analysis, quantitative concordance/discordance scores incorporating data variabilities, pathway-centric downstream investigation, knowledge retrieval by text mining, and topological gene module detection for hypothesis generation. Instead of a genome-wide vague and dichotomous answer of "poorly" or "greatly" mimicking humans, CAMO assists researchers to numerically quantify congruence, to dissect true cross-species differences from unwanted biological or cohort variabilities, and to visually identify molecular mechanisms and pathway subnetworks that are best or least mimicked by model organisms, which altogether provides foundations for hypothesis generation and subsequent translational decisions.

EstroGene database reveals diverse temporal, context-dependent and directional estrogen receptor regulomes in breast cancer

As one of the most successful cancer therapeutic targets, estrogen receptor-α (ER/ESR1) has been extensively studied in decade-long. Sequencing technological advances have enabled genome-wide analysis of ER action. However, reproducibility is limited by different experimental design. Here, we established the EstroGene database through centralizing 246 experiments from 136 transcriptomic, cistromic and epigenetic datasets focusing on estradiol-treated ER activation across 19 breast cancer cell lines. We generated a user-friendly browser ( https://estrogene.org/ ) for data visualization and gene inquiry under user-defined experimental conditions and statistical thresholds. Notably, documentation-based meta-analysis revealed a considerable lack of experimental details. Comparison of independent RNA-seq or ER ChIP-seq data with the same design showed large variability and only strong effects could be consistently detected. We defined temporal estrogen response metasignatures and showed the association with specific transcriptional factors, chromatin accessibility and ER heterogeneity. Unexpectedly, harmonizing 146 transcriptomic analyses uncovered a subset of E2-bidirectionally regulated genes, which linked to immune surveillance in the clinical setting. Furthermore, we defined context dependent E2 response programs in MCF7 and T47D cell lines, the two most frequently used models in the field. Collectively, the EstroGene database provides an informative resource to the cancer research community and reveals a diverse mode of ER signaling.

The Great Immune Escape: Understanding the Divergent Immune Response in Breast Cancer Subtypes

Breast cancer, the most common type of cancer affecting women, encompasses a collection of histologic (mainly ductal and lobular) and molecular subtypes exhibiting diverse clinical presentation, disease trajectories, treatment options, and outcomes. Immunotherapy has revolutionized treatment for some solid tumors but has shown limited promise for breast cancers. In this review, we summarize recent advances in our understanding of the complex interactions between tumor and immune cells in subtypes of breast cancer at the cellular and microenvironmental levels. We aim to provide a perspective on opportunities for future immunotherapy agents tailored to specific features of each subtype of breast cancer.

SIGNIFICANCE

Although there are currently over 200 ongoing clinical trials testing immunotherapeutics, such as immune-checkpoint blockade agents, these are largely restricted to the triple-negative and HER2+ subtypes and primarily focus on T cells. With the rapid expansion of new in vitro, in vivo, and clinical data, it is critical to identify and highlight the challenges and opportunities unique for each breast cancer subtype to drive the next generation of treatments that harness the immune system.

Multiomics in primary and metastatic breast tumors from the AURORA US network finds microenvironment and epigenetic drivers of metastasis

The AURORA US Metastasis Project was established with the goal to identify molecular features associated with metastasis. We assayed 55 females with metastatic breast cancer (51 primary cancers and 102 metastases) by RNA sequencing, tumor/germline DNA exome and low-pass whole-genome sequencing and global DNA methylation microarrays. Expression subtype changes were observed in ~30% of samples and were coincident with DNA clonality shifts, especially involving HER2. Downregulation of estrogen receptor (ER)-mediated cell-cell adhesion genes through DNA methylation mechanisms was observed in metastases. Microenvironment differences varied according to tumor subtype; the ER+/luminal subtype had lower fibroblast and endothelial content, while triple-negative breast cancer/basal metastases showed a decrease in B and T cells. In 17% of metastases, DNA hypermethylation and/or focal deletions were identified near HLA-A and were associated with reduced expression and lower immune cell infiltrates, especially in brain and liver metastases. These findings could have implications for treating individuals with metastatic breast cancer with immune- and HER2-targeting therapies.

Clinicopathological Features and Outcomes Comparing Patients With Invasive Ductal and Lobular Breast Cancer

BACKGROUND

There is increasing interest in better understanding the biology and clinical presentation of invasive lobular cancer (ILC), which is the most common special histological subtype of breast cancer. Limited large contemporary data sets are available allowing comparison of clinicopathologic features between ILC and invasive ductal cancer (IDC).

METHODS

The Great Lakes Breast Cancer Consortium was formed to compare clinical behavior of ILC (n = 3617) and IDC (n = 30 045) from 33 662 patients treated between 1990 and 2017 at 3 large clinical centers. We used Kaplan-Meier analysis, Cox proportional hazards modeling, and propensity score matching to evaluate treatment differences and outcomes. All statistical testing used 2-sided P values.

RESULTS

Compared with IDC, patients with ILC were more frequently diagnosed at later stages and with more lymph node involvement (corrected P < .001). Estrogen receptor-positive ILCs were of lower grade (grade 1 and 2: 90% in ILC vs 72% in IDC) but larger in size (T3 and 4: 14.3% in ILC vs 3.4% in IDC) (corrected P < .001), and since 1990, the mean ILC size detected at diagnosis increased yearly. Patients with estrogen receptor (ER)-positive ILC underwent statistically significantly more mastectomies compared with ER-positive IDC (57% vs 46%). Using Kaplan-Meier analysis, patients with ER-positive ILC had statistically significantly worse disease-free survival and overall survival than ER-positive IDC although 6 times more IDCs were classified as high risk by OncotypeDx Breast Recurrence Score assay.

CONCLUSIONS

This large, retrospective, collaborative analysis with 3 clinical centers identified meaningful differences in clinicopathological features between ILC and IDC, providing further evidence that these are 2 different entities requiring different clinical management.

Loss of E-cadherin Induces IGF1R Activation and Reveals a Targetable Pathway in Invasive Lobular Breast Carcinoma

No special-type breast cancer [NST; commonly known as invasive ductal carcinoma (IDC)] and invasive lobular carcinoma (ILC) are the two major histological subtypes of breast cancer with significant differences in clinicopathological and molecular characteristics. The defining pathognomonic feature of ILC is loss of cellular adhesion protein, E-cadherin (CDH1). We have previously shown that E-cadherin functions as a negative regulator of the IGF1R and propose that E-cadherin loss in ILC sensitizes cells to growth factor signaling that thus alters their sensitivity to growth factor-signaling inhibitors and their downstream activators. To investigate this potential therapeutic vulnerability, we generated CRISPR-mediated CDH1 knockout (CDH1 KO) IDC cell lines (MCF7, T47D, and ZR75.1) to uncover the mechanism by which loss of E-cadherin results in IGF pathway activation. CDH1 KO cells demonstrated enhanced invasion and migration that was further elevated in response to IGF1, serum and collagen I. CDH1 KO cells exhibited increased sensitivity to IGF resulting in elevated downstream signaling. Despite minimal differences in membranous IGF1R levels between wild-type (WT) and CDH1 KO cells, significantly higher ligand-receptor interaction was observed in the CDH1 KO cells, potentially conferring enhanced downstream signaling activation. Critically, increased sensitivity to IGF1R, PI3K, Akt, and MEK inhibitors was observed in CDH1 KO cells and ILC patient-derived organoids.

IMPLICATIONS

Overall, this suggests that these targets require further exploration in ILC treatment and that CDH1 loss may be exploited as a biomarker of response for patient stratification.

Abstract 881: Functional characterization and clinical prevalence of ESR1 fusions in advanced endocrine resistant breast cancer

Breast cancer is the most frequently diagnosed women’s cancer with one in eight women being diagnosed during her lifetime. Even within the era of novel therapeutics and improved diagnostic screening, mortality rates remain stable, with an estimated 43,600 deaths occurring in 2021 alone. Malignant transformation of breast cancer to advanced disease is the mainstay driver of mortality and is largely attributed to treatment resistance. Fortunately, identification of therapeutic resistance mechanisms has rapidly transformed our understanding of cancer evasion and is enabling personalized treatment guidance and regimens. Notably, mutations and amplifications of the estrogen receptor (ER, ESR1) remain a prominent source of endocrine therapy resistance. Two-thirds of all breast cancer are positive for ER, making it an attractive druggable target that has revolutionized patient care. Unfortunately, however, 30% of tumors exposed to endocrine therapy subsequently develop resistance. Recently identified ER fusion proteins, discovered in metastatic endocrine resistant disease, result in substitution of the C-terminal ER ligand binding domain with a variety of fused in-frame gene partners due to genomic translocations occurring at ESR1 intron 6. Our lab has demonstrated that ER fusion proteins drive robust ER hyperactivation and metastatic potential independent of ER ligand estradiol and in the presence of endocrine therapies such as Fulvestrant. Overexpression of ER fusion proteins, particularly ESR1-SOX9 and ESR1-YAP1, in invasive ductal and invasive lobular carcinoma models with endogenous ER knockdown, resulted in enhanced 3D growth, increased migration by wound scratch assay, enhanced colony formation as a proxy of survival and ER pathway hyperactivation. Targeted endogenous ER knockdown utilizing short hairpin RNA enhances fusion driven phenotypes that are unique to fusion partner proteins and cellular context. ER fusion prevalence in advanced breast cancer remains a challenge due to the promiscuous fusion partners disrupting traditional sequencing techniques. Our group has developed a novel method of ESR1 fusion detection, EnRich, which exploits the ESR1 fusion breakpoints. Importantly, we have assessed EnRich fusion discovery with the use of circulating tumor DNA (ctDNA) from non-invasive liquid biopsies. Our pilot study of 15 blood ctDNA samples detected three different ESR1 fusion gene products which have been further validated by PCR. EnRich will help us identify and monitor treatment resistance of known and novel ESR1 fusions in patient samples irrespective of treatment regime as well as in longitudinal samples collected from patients undergoing endocrine therapy. The in vitro resistant properties of fusion proteins emphasize the necessity to better detect these mutations in patients with advanced disease and determine more appropriate treatment measures.

Citation Format: Megan E. Yates, Tiantong Liu, Jagmohan Hooda, Jennifer M. Atkinson, Adrian V. Lee, Steffi Oesterreich. Functional characterization and clinical prevalence of ESR1 fusions in advanced endocrine resistant breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 881.

Abstract 5031: Improved deconvolution of combined bulk and single-cell RNA-sequencing data

Identifying how cell types and their abundances evolve during tumor progression is critical to understanding the mechanisms and identifying predictors of metastasis. Single-cell RNA sequencing (scRNA-seq) has been especially promising in resolving heterogeneity of expression programs at the single cell level but is not always available, for example for large cohort studies or longitudinal analysis of archived samples. In such cases, cell subpopulations must be inferred by deconvolution, a process that can infer single-cell genomic data from bulk data but has limited ability to resolve fine clonal structure. We extend our previous bulk genomic deconvolution tool, Robust and Accurate Deconvolution (RAD), to establish a new method, scRAD, that can use reference scRNA-seq to interpret sample collections for which only bulk RNA-seq is available for some samples, e.g., clonally resolving archived primary (PRM) tissues using scRNA-seq from metastases (METs). We preprocess scRNA-seq data to accurately represent gene expression profiles, yielding a signature matrix S then extend our RAD method via a regularization term to deconvolve bulk data while maximizing consistency with S. We validate our method on semi-synthetic data derived from human PRM breast cancer cases and bone and ovary METs, showing that scRAD improves inference of single-cell gene expression profiles and their frequencies relative to the prior RAD with random initialization or initialization using the single-cell matrix S (Table 1). We then apply scRAD to a collection of paired PRM and MET tumors to quantify progression changes in common cell types. One-sided Kaplan-Meier analysis shows that tumors inferred to increase the mast cell fraction from PRM to MET exhibit lower overall survival (p&lt;0.05), consistent with the role of mast cells in metastatic growth and propagation. Tumors that show increased macrophage cell fraction from PRM to MET show improved overall survival (p&lt;0.04), consistent with the role of immune infiltration in survival.

mean square error (MSE) of gene expression and mixture fraction inference on semi-simulated data Method RAD with random initialization RAD with random initialization RAD with random initialization RAD initialized with S RAD initialized with S RAD initialized with S scRAD scRAD scRAD Sample number 2 4 8 2 4 8 2 4 8 Gene Expression MSE 0.41 0.28 0.37 0.31 0.27 0.37 0.22 0.15 0.13 Mixture fraction MSE 0.58 0.82 0.59 0.39 0.83 0.58 0.19 0.20 0.22

Citation Format: Haoyun Lei, Xiaoyan A. Guo, Yifeng Tao, Kai Ding, Xuecong Fu, Steffi Oesterreich, Adrian V. Lee, Russell Schwartz. Improved deconvolution of combined bulk and single-cell RNA-sequencing data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5031.

Abstract 81: In depth single cell profiling of a case of breast cancer bone metastases with associated organoid models reveal a precision medicine approach to treatment

Occurring in 65-80% of metastatic breast cancer (BC), bone metastasis (BoM) is the major cause of BC related mortality. Lack of understanding of BoM evolution and heterogeneity at patient specific levels and precise application of targeted therapies are major challenges of managing BoM. In this study, we described in depth histologic and molecular characterization of a case of invasive lobular breast cancer (ILC) bilateral metastasis to bone, and identified and tested potential targets for treating BoM. H&E/IHC staining, whole exome sequencing (WES), and RNAseq were performed on FFPE primary tumor and pelvis and tibia BoMs collected from our study case. Organoids were derived from the two BoMs and single cell RNA sequencing (scRNAseq) undertaken on the organoids and their originating tumors. H&E/IHC demonstrated evolution of the disease from an ER+ primary ILC to ER- BoM with mixed lobular and ductal carcinoma features. WES revealed two druggable mutations including PIK3CA (E545K) and BRCA1 (D1813H/H399P). RNAseq revealed upregulation of TGF-β, Wnt/beta-catenin and PI3K pathways, epithelial to mesenchymal transition (EMT) and angiogenesis in BoM compared to primary tumor, representing promising targets for BoM. scRNAseq revealed 5 major cell types including epithelial, fibroblasts, immune, osteoclasts and endothelial cells, and pronounced intracellular heterogeneity. Six epithelial clusters were identified, featuring high TNF-α signaling, high partial EMT signatures regulated by PRRX1/2, TWIST1/2, and FOXS1, high proliferation, and endocrine resistance signatures. In fibroblasts, 9 clusters were identified representing ECM remodeling, angiogenesis, osteoclast-like, MSC, IFN response and myofibroblasts. Immune cells majorly composed of monocytes/macrophages, CD4+, CD8+ and Treg T cells, and NK cell. WES and scRNAseq analysis demonstrated that organoids preserved mutational landscape and cellular heterogeneity of matched BoMs. Consistent with the BRCA1 and PIK3CA mutations, organoids were responsive to a PARP (Talazoparib: IC50 1.3uM) and PI3K (Alpelisib: IC50 4-9uM) inhibitors. In summary, we have identified potential therapeutic targets from understanding evolution and heterogeneity of BC BoM, and evaluated these in patient-specific organoids, thereby providing insights for the design of a precision medicine based clinical treatment strategy.

Citation Format: Kai Ding, Fangyuan Chen, Nolan Priedigkeit, Daniel D. Brown, Tanya Heim, Rebecca Watters, Kurt Weiss, Peter C. Lucas, Jennifer M. Atkinson, Steffi Oesterreich, Adrian V. Lee. In depth single cell profiling of a case of breast cancer bone metastases with associated organoid models reveal a precision medicine approach to treatment [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 81.

Abstract 2690: Loss of E-cadherin induces IGF1R activation revealing a targetable pathway in invasive lobular breast carcinoma

Invasive Ductal Carcinoma (IDC) and Invasive Lobular Carcinoma (ILC) are two major subtypes of breast cancer with significant differences in their histological and molecular underpinnings. ILC has a unique loss of E-cadherin (CDH1) which we have previously demonstrated as a negative regulator of the Insulin-like Growth Factor 1 receptor (IGF1R) through a comprehensive analysis of cell line models and tumor samples. We propose that loss of E-cadherin in ILC sensitizes cells to growth factor signaling and thus alters their sensitivity to growth factor signaling inhibitors. To investigate this, we used CRISPR to generate CDH1 knockout (KO) IDC cell lines (MCF7, T47D, ZR75.1) to uncover the mechanism by which E-cadherin loss activates the IGF pathway while also assessing its targetability. CDH1 KO cells exhibited anchorage independent growth in suspension culture and altered p120 catenin localization as observed in ILC tumors. Through in vitro studies, we show increased signaling sensitivity to IGF/insulin ligands and enhanced signaling duration in CDH1 KO cells. In addition, we observed a higher migratory potential of CDH1 KO cells compared to wild type (WT) cells, which was further enhanced as a chemotactic response to IGF1 or serum. Further, this phenotype could be reversed with an IGF1R inhibitor, BMS-754807. We additionally identified an increase in Collagen I haptotaxis in the CDH1 KO cells, which was also translated into a novel invasive phenotype towards serum in the T47D CDH1 KO cells. Despite no significant differences in membranous IGF1R levels between WT and CDH1 KO cells, higher ligand-receptor interaction was observed with CDH1 KO cells, demonstrating an increased ligand-receptor complex formation upon stimulation. Our results suggest that loss of CDH1 results in an increase in IGF1 receptor availability for ligand binding which in turn allows for an enhanced downstream signaling activation. Interestingly, a physical repression of E-cadherin on IGF1R could not be demonstrated, suggesting spatial changes on the membrane following E-cadherin loss may control ligand binding. Critically, increased sensitivity to IGF1R, PI3K, AKT and MEK inhibitors was observed in CDH1 KO cells suggesting that these targets should be further explored in ILC and that CDH1 loss may be exploited as a biomarker of response, or for patient stratification to inhibitors targeting these pathways.

Citation Format: Ashuvinee Elangovan, Laura Savariau, Megan E. Yates, Jagmohan Hooda, Alison M. Nagle, Steffi Oesterreich, Jennifer M. Atkinson, Adrian V. Lee. Loss of E-cadherin induces IGF1R activation revealing a targetable pathway in invasive lobular breast carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2690.

Abstract 956: Role of E-cadherin in progression of invasive breast lobular carcinoma

Invasive lobular carcinoma (ILC) is the second most frequently diagnosed histologic subtype of invasive breast cancer following invasive ductal carcinoma (IDC) and accounts for 15% of all cases. The hallmark of ILC is the genetic loss of E-cadherin (CDH1) causing the disruption of adherens junctions and resulting in discohesive, linear growth of ILC cells in tissues. Beyond its effects on this unique histology, there is limited data on the role of E-cadherin loss in ILC metastasis. Therefore, we generated three ILC cell lines (MDA-MB-134, SUM44PE and BCK4) with inducible E-cadherin overexpression. Using immunoblotting and immunofluorescence microscopy we confirmed E-cadherin expression, re-expression of other adherens junction proteins including α and β-catenin, and re-localization of cytoplasmic p120-catenin to the membrane. Successful cell adhesion to E-cadherin coated plates suggested formation of functional junctions. E-cadherin expression had no effect on 2D growth but did diminish cell growth in ultra-low attachment (ULA) conditions consistent with its previously described role in anoikis. We confirmed increased cell death of MM134 with E-cadherin overexpression compared to control cells when grown in ULA. Overexpression of E-cadherin failed to rescue poor migratory and invasive ability of ILC cell lines as measured by transwell assays. However, E-cadherin expression in SUM44PE lead to decreased haptotaxis to collagen I, and altered morphology from loose grape-like structure to tight sphere-like formation when embedded in collagen I.In vivo analyses revealed reduced primary tumor formation in E-cadherin-overexpression SUM44PE cells after mammary fat pad injection. Additional studies aimed at understanding E-cadherin’s effect on in vivo metastasis are ongoing.Collectively, our studies using novel cell line models will lead to an improved understanding of the hallmark loss of E-cadherin and disease mechanisms in ILC, which we hope will ultimately enable the development of more effective therapies and improve the outcome of patients with this understudied histological subtype of breast cancer.

Citation Format: Laura Savariau, Nilgun Tasdemir, Ashuvinee Elangovan, Kai Ding, Pooja Sree Tallapaneni, Jagmohan Hooda, Jennifer M. Atkinson, Adrian V. Lee, Steffi Oesterreich. Role of E-cadherin in progression of invasive breast lobular carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 956.

Isoforms of Neuropilin-2 Denote Unique Tumor-Associated Macrophages in Breast Cancer

Tumor-associated macrophages (TAMs) exert profound influence over breast cancer progression, promoting immunosuppression, angiogenesis, and metastasis. Neuropilin-2 (NRP2), consisting of the NRP2a and NRP2b isoforms, is a co-receptor for heparin-binding growth factors including VEGF-C and Class 3 Semaphorins. Selective upregulation in response to environmental stimuli and independent signaling pathways endow the NRP2 isoforms with unique functionality, with NRP2b promoting increased Akt signaling via receptor tyrosine kinases including VEGFRs, MET, and PDGFR. Although NRP2 has been shown to regulate macrophage/TAM biology, the role of the individual NRP2a/NRP2b isoforms in TAMs has yet to be evaluated. Using transcriptional profiling and spectral flow cytometry, we show that NRP2 isoform expression was significantly higher in TAMs from murine mammary tumors. NRP2a/NRP2b levels in human breast cancer metastasis were dependent upon the anatomic location of the tumor and significantly correlated with TAM infiltration in both primary and metastatic breast cancers. We define distinct phenotypes of NRP2 isoform-expressing TAMs in mouse models of breast cancer and within malignant pleural effusions from breast cancer patients which were exclusive of neuropilin-1 expression. Genetic depletion of either NRP2 isoform in macrophages resulted in a dramatic reduction of LPS-induced IL-10 production, defects in phagosomal processing of apoptotic breast cancer cells, and increase in cancer cell migration following co-culture. By contrast, depletion of NRP2b, but not NRP2a, inhibited production of IL-6. These results suggest that NRP2 isoforms regulate both shared and unique functionality in macrophages and are associated with distinct TAM subsets in breast cancer.

Effect of socioeconomic status as measured by Neighborhood Deprivation Index on survival in metastatic breast cancer

1013

Background: Socioeconomic status (SES) and race are major determinants of health outcomes in the United States. We aim to assess the effect of SES as measured by the Neighborhood Deprivation Index (NDI) and race on outcomes in metastatic breast cancer patients at our center. Methods: The NDI scores for patients with metastatic breast cancer who were treated at our center between 2000 and 2017 were obtained from the Neighborhood Atlas using their Zip-Code (N = 1246). The SES groups were defined as low deprivation with an NDI score in the bottom tertile and high deprivation with NDI in the top or middle tertiles. Baseline characteristics were compared between the SES groups with Bonferroni correction. Univariate and multivariate survival analysis were performed using the R packages “survival” and “survminer”. Results: Race was the only baseline characteristic that was significantly different between the SES groups, the high deprivation group had a higher proportion of African Americans (10.5%) than the low deprivation group (3.7%, P = 9.3e-05). In univariate Kaplan-Meier survival analysis, both SES and race had significant effect on overall survival such that the high deprivation group had worse survival than low deprivation (Log Rank P = 0.01) and African Americans had worse survival than Caucasians (P = 0.008). In multivariate Cox proportional hazard model, SES, but not race, had a significant effect on overall survival (hazard ratio for high deprivation was 1.19 [95% Confidence interval 1.04 - 1.37], P = 0.01; Table). Progression-free survival on first-line chemotherapy was not different between the SES groups or racial groups in both univariate and multivariate analysis. Conclusions: The current study shows that patients from the high deprivation group (i.e., low SES), have worse survival in metastatic breast cancer. Race was no longer a significant predictor of survival when SES was accounted for in the analysis. This possibly suggests that poor outcomes in the African American population is explained by the association between low SES and African American race. Based on these results, there is an urgent need for healthcare investments in the low SES neighborhoods. [Table: see text]

Mutual exclusivity of ESR1 and TP53 mutations in endocrine resistant metastatic breast cancer

Both TP53 and ESR1 mutations occur frequently in estrogen receptor positive (ER+) metastatic breast cancers (MBC) and their distinct roles in breast cancer tumorigenesis and progression are well appreciated. Recent clinical studies discovered mutual exclusivity between TP53 and ESR1 mutations in metastatic breast cancers; however, mechanisms underlying this intriguing clinical observation remain largely understudied and unknown. Here, we explored the interplay between TP53 and ESR1 mutations using publicly available clinical and experimental data sets. We first confirmed the robust mutational exclusivity using six independent cohorts with 1,056 ER+ MBC samples and found that the exclusivity broadly applies to all ER+ breast tumors regardless of their clinical and distinct mutational features. ESR1 mutant tumors do not exhibit differential p53 pathway activity, whereas we identified attenuated ER activity and expression in TP53 mutant tumors, driven by a p53-associated E2 response gene signature. Further, 81% of these p53-associated E2 response genes are either direct targets of wild-type (WT) p53-regulated transactivation or are mutant p53-associated microRNAs, representing bimodal mechanisms of ER suppression. Lastly, we analyzed the very rare cases with co-occurrences of TP53 and ESR1 mutations and found that their simultaneous presence was also associated with reduced ER activity. In addition, tumors with dual mutations showed higher levels of total and PD-L1 positive macrophages. In summary, our study utilized multiple publicly available sources to explore the mechanism underlying the mutual exclusivity between ESR1 and TP53 mutations, providing further insights and testable hypotheses of the molecular interplay between these two pivotal genes in ER+ MBC.

ESR1 mutant breast cancers show elevated basal cytokeratins and immune activation

Estrogen receptor alpha (ER/ESR1) is frequently mutated in endocrine resistant ER-positive (ER+) breast cancer and linked to ligand-independent growth and metastasis. Despite the distinct clinical features of ESR1 mutations, their role in intrinsic subtype switching remains largely unknown. Here we find that ESR1 mutant cells and clinical samples show a significant enrichment of basal subtype markers, and six basal cytokeratins (BCKs) are the most enriched genes. Induction of BCKs is independent of ER binding and instead associated with chromatin reprogramming centered around a progesterone receptor-orchestrated insulated neighborhood. BCK-high ER+ primary breast tumors exhibit a number of enriched immune pathways, shared with ESR1 mutant tumors. S100A8 and S100A9 are among the most induced immune mediators and involve in tumor-stroma paracrine crosstalk inferred by single-cell RNA-seq from metastatic tumors. Collectively, these observations demonstrate that ESR1 mutant tumors gain basal features associated with increased immune activation, encouraging additional studies of immune therapeutic vulnerabilities.

Hotspot ESR1 Mutations Are Multimodal and Contextual Modulators of Breast Cancer Metastasis

Constitutively active estrogen receptor α (ER/ESR1) mutations have been identified in approximately one-third of ER+ metastatic breast cancers. Although these mutations are known as mediators of endocrine resistance, their potential role in promoting metastatic disease has not yet been mechanistically addressed. In this study, we show the presence of ESR1 mutations exclusively in distant but not local recurrences in five independent breast cancer cohorts. In concordance with transcriptomic profiling of ESR1-mutant tumors, genome-edited ESR1 Y537S and D538G-mutant cell models exhibited a reprogrammed cell adhesive gene network via alterations in desmosome/gap junction genes and the TIMP3/MMP axis, which functionally conferred enhanced cell-cell contacts while decreasing cell-extracellular matrix adhesion. In vivo studies showed ESR1-mutant cells were associated with larger multicellular circulating tumor cell (CTC) clusters with increased compactness compared with ESR1 wild-type CTCs. These preclinical findings translated to clinical observations, where CTC clusters were enriched in patients with ESR1-mutated metastatic breast cancer. Conversely, context-dependent migratory phenotypes revealed cotargeting of Wnt and ER as a vulnerability in a D538G cell model. Mechanistically, mutant ESR1 exhibited noncanonical regulation of several metastatic pathways, including secondary transcriptional regulation and de novo FOXA1-driven chromatin remodeling. Collectively, these data provide evidence for ESR1 mutation-modulated metastasis and suggest future therapeutic strategies for targeting ESR1-mutant breast cancer.

SIGNIFICANCE

Context- and allele-dependent transcriptome and cistrome reprogramming in mutant ESR1 cell models elicit diverse metastatic phenotypes related to cell adhesion and migration, which can be pharmacologically targeted in metastatic breast cancer.

Abstract PD1-08: Esr1 mutant breast cancers show elevated basal cytokeratins and immune activation

Estrogen receptor alpha (ER/ESR1) is mutated in 30-40% of endocrine resistant ER-positive (ER+) breast cancer. ESR1 mutations cause ligand-independent growth and increased metastasis in vivo and in vitro. Despite the distinct clinical features and changes in therapeutic response associated with ESR1 mutations, there are no data about their potential role in intrinsic subtype switching.Applying five different luminal and basal gene set pairs derived from cell lines and tumors, ESR1 mutant cell models and clinical samples showed a significant enrichment of basal subtype markers. Among them, the six basal cytokeratins (BCKs) were the most enriched genes (KRT5/6A/6B/14/16/17) uniquely in ESR1 mutant cells but not other endocrine resistant cell models. BCKs were observed to heterogeneously express in a minor cell subpopulation in ESR1 mutant cell models and clinical specimens. ER ChIP-seq showed the mutant-specific induction of BCKs was independent of ER binding and instead selectively expressed in clones with low ER expression. In contrast, BCKs are associated with chromatin reprogramming centered around a progesterone receptor-orchestrated 154 kb insulated neighborhood at the KRT14/16/17 genomic region. Stronger CTCF binding was detected at the bases of chromatin loop in ESR1 mutant cells. Knockdown of progesterone receptor but not glucocorticoid receptor drastically blocked the induction of KRT14/16/17 in ESR1 mutant cells. Unexpectedly, high BCK expression in ER+ primary breast cancer is associated with good prognosis, and these tumors show enriched activation of a number of immune pathways, a distinctive feature shared with ESR1 mutant tumors. While the BCK-associated immune activation is not related to tumor mutation burdens, S100A8 and S100A9 were identified as the most highly induced immune mediators shared between high-BCKs ER+ and ESR1 mutant tumors, which was further validated in the plasma samples of a cohort of 18 patients with ER+ metastases (11 WT vs 7 mutant). Finally, single-cell RNA-seq analysis in an ER+ bone metastasis case inferred the involvement of S100A8 and S100A9 in paracrine crosstalk between epithelial and stromal cells, particularly macrophages and fibroblasts through TLR4 signaling. Collectively, these observations demonstrate that ESR1 mutant tumors gain basal features with induction of basal cytokeratins via epigenetic mechanisms in rare subpopulation of cells. This is associated with increased immune activation, encouraging additional studies of immune therapeutic vulnerabilities in ESR1 mutant tumors.

Citation Format: Zheqi Li, Yang Wu, Olivia Mcginn, Amir Bahreini, Nolan M. Priedigkeit, Kai Ding, Sayali Onkar, Carol A. Sartorius, Lori Miller, Margaret Rosenzweig, Nikhil Wagle, Jennifer K. Richer, William J. Muller, Laki Buluwela, Simak Ali, Dario A.A. Vignali, Yusi Fang, Li Zhu, George C. Tseng, Jason Gertz, Jennifer M. Atkinson, Adrian V. Lee, Steffi Oesterreich. Esr1 mutant breast cancers show elevated basal cytokeratins and immune activation [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr PD1-08.

Abstract P5-12-03: Wcrc-25: A novel luminal invasive lobular carcinoma cell line model

Breast cancer is categorized by the molecular and histologic presentation of the tumor; with the major histologic subtypes observed in patients being Invasive Ductal Carcinoma (IDC) and Invasive Lobular Carcinoma (ILC). ILC are characterized by growth in a single file discohesive manner with stromal infiltration attributed to their hallmark loss of E-cadherin. The lower incidence of ILC, relative to IDC, resulted in this subtype being understudied in the field. The Cancer Cell Line Encyclopedia (CCLE) currently only lists 2 ILC cell lines; emphasizing the need to expand the inventory of models available to researchers. To that end, here we report the successful establishment and characterization of a novel ILC cell line from a metastatic pleural aspirate from a postmenopausal Caucasian woman with ER+ metastatic ILC. The patient underwent a bilateral mastectomy, several pleural fluid aspirations and multiple treatment changes as her disease progressed until she succumbed to the disease. Various cell culture methods were attempted with the pleural aspirate and a standard practice of DMEM supplemented with 10% FBS in normoxic conditions was deemed effective. A cell line was developed, WCRC-25 which is ER-. WCRC-25 appears to be E2 unresponsive in vitro and is morphologically larger than other ILC cells currently in use, namely, MDA-MB-134, BCK4 and IPH-926, while exhibiting better adherence to 2D plates and a less rounded morphology. Sequencing revealed a CDH1 Q706 truncating mutation in the cells, which was also observed in the patient’s cfDNA with its enrichment increasing over time, supporting elevated disease burden. RNA-seq analysis comparing the primary tumor to metastases and the cell line reveal a signature for cell cycle progression and Akt signaling. To assess targetability, we subjected WCRC-25 to AZD5363 and Alpelisib which confirmed WCRC-25 as susceptible to PI3K/Akt signaling inhibition. Finally, we generated mouse xenograft models of WCRC-25 through mammary fat pad injection and observed small, but palpable tumors and metastasis to several sites. In conclusion, we report WCRC-25 as a novel ILC cell line with much promise as an invaluable research tool to advance our understanding of ILC and its therapeutic vulnerabilities.

Citation Format: Ashuvinee Elangovan, Emily A Bossart, Ahmed Basudan, Nilgun Tasdemir, Osama Shiraz Shah, Kai Ding, Carolin Meier, Tanya Heim, Lori Miller, Tiantong Liu, Shannon L Puhalla, Grzegorz Gurda, Peter C Lucas, Priscilla F McAuliffe, Jennifer M Atkinson, Adrian V Lee, Steffi Oesterreich. Wcrc-25: A novel luminal invasive lobular carcinoma cell line model [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P5-12-03.

Fusion-Associated Carcinomas of the Breast: Diagnostic, Prognostic, and Therapeutic Significance

Recurrent gene fusions comprise a class of viable genetic targets in solid tumors that have culminated several recent break-through cancer therapies. Their role in breast cancer, however, remains largely underappreciated due to the complexity of genomic rearrangements in breast malignancy. Just recently, we and others have identified several recurrent gene fusions in breast cancer with important clinical and biological implications. Examples of the most significant recurrent gene fusions to date include 1) ESR1-CCDC170 gene fusions in luminal B and endocrine resistant breast cancer that exert oncogenic function via modulating the HER2/HER3/SRC complex, 2) ESR1 exon 6 fusions in metastatic disease that drive estrogen-independent ER transcriptional activity, 3) BCL2L14-ETV6 fusions in a more aggressive form of the triple negative subtype that prime epithelial-mesenchymal transition and endow paclitaxel resistance, 4) the ETV6-NTRK3 fusion in secretory breast carcinoma that constitutively activates NTRK3 kinase, 5) the oncogenic MYB-NFIB fusion as a genetic driver underpinning adenoid cystic carcinomas of the breast that activates MYB pathway, and 6) the NOTCH/MAST kinase gene fusions that activate NOTCH and MAST signaling. Importantly, these fusions are enriched in more aggressive and lethal breast cancer presentations and appear to confer therapeutic resistance. Thus, these gene fusions could be utilized as genetic biomarkers to identify patients that require more intensive treatment and surveillance. In addition, kinase fusions are currently being evaluated in breast cancer clinical trials and on-going mechanistic investigation is exposing therapeutic vulnerabilities in patients with fusion positive disease. This article is protected by copyright. All rights reserved.

Mapping molecular subtype specific alterations in breast cancer brain metastases identifies clinically relevant vulnerabilities

The molecular events and transcriptional plasticity driving brain metastasis in clinically relevant breast tumor subtypes has not been determined. Here we comprehensively dissect genomic, transcriptomic and clinical data in patient-matched longitudinal tumor samples, and unravel distinct transcriptional programs enriched in brain metastasis. We report on subtype specific hub genes and functional processes, central to disease-affected networks in brain metastasis. Importantly, in luminal brain metastases we identify homologous recombination deficiency operative in transcriptomic and genomic data with recurrent breast mutational signatures A, F and K, associated with mismatch repair defects, TP53 mutations and homologous recombination deficiency (HRD) respectively. Utilizing PARP inhibition in patient-derived brain metastatic tumor explants we functionally validate HRD as a key vulnerability. Here, we demonstrate a functionally relevant HRD evident at genomic and transcriptomic levels pointing to genomic instability in breast cancer brain metastasis which is of potential translational significance.

Personalising therapy for early-stage oestrogen receptor-positive breast cancer in older women

Age is one of the most important risk factors for the development of breast cancer. Nearly a third of all breast cancer cases occur in older women (aged ≥70 years), with most cases being oestrogen receptor-positive (ER+). Such tumours are often indolent and unlikely to be the ultimate cause of death for older women, particularly when considering other comorbidities. This Review focuses on unique clinical considerations for screening, detection, and treatment regimens for older women who develop ER+ breast cancers-specifically, we focus on recent trends for de-implementation of screening, staging, surgery, and adjuvant therapies along the continuum of care. Additionally, we also review emerging basic and translational research that will further uncover the unique underlying biology of these tumours, which develop in the context of systemic age-related inflammation and changing hormone profiles. With prevailing trends of clinical de-implementation, new insights into mechanistic biology might provide an opportunity for precision medicine approaches to treat patients with well tolerated, low-toxicity agents to extend patients' lives with a higher quality of life, prevent tumour recurrences, and reduce cancer-related burdens.

OUP accepted manuscript

Motivation

Identifying cell types and their abundances and how these evolve during tumor progression is critical to understanding the mechanisms of metastasis and identifying predictors of metastatic potential that can guide the development of new diagnostics or therapeutics. Single-cell RNA sequencing (scRNA-seq) has been especially promising in resolving heterogeneity of expression programs at the single-cell level, but is not always feasible, e.g. for large cohort studies or longitudinal analysis of archived samples. In such cases, clonal subpopulations may still be inferred via genomic deconvolution, but deconvolution methods have limited ability to resolve fine clonal structure and may require reference cell type profiles that are missing or imprecise. Prior methods can eliminate the need for reference profiles but show unstable performance when few bulk samples are available.

Results

In this work, we develop a new method using reference scRNA-seq to interpret sample collections for which only bulk RNA-seq is available for some samples, e.g. clonally resolving archived primary tissues using scRNA-seq from metastases. By integrating such information in a Quadratic Programming framework, our method can recover more accurate cell types and corresponding cell type abundances in bulk samples. Application to a breast tumor bone metastases dataset confirms the power of scRNA-seq data to improve cell type inference and quantification in same-patient bulk samples.

Availability and implementation

Source code is available on Github at https://github.com/CMUSchwartzLab/RADs.

The Development and Implementation of an Autopsy/ Tissue Donation for Breast Cancer Research

There is growing interest in tissue procurement for cancer research through autopsy. Establishing an autopsy/tissue donation programme for breast cancer research within an academic medical centre in the United States requires consideration, planning, multi-departmental collaboration and labour-intensive maintenance. It is the purpose of this paper to outline the necessary considerations in implementing and maintaining a tissue donation and autopsy programme within a breast cancer centre at a comprehensive cancer centre. Considerations of programme planning include: patient engagement, the recruitment of patients and families into the programme, the role and scope of work of the clinical coordinator, regulatory issues and the coordination with both pathology and the research team at time of death and autopsy/tissue donation. All aspects of the tissue donation/rapid autopsy programme development and implementation are discussed and illustrated through case study. An Autopsy/ Tissue Donation for breast cancer research can be successfully developed and implemented.