Automated mitotic spindle hotspot counts are highly associated with clinical outcomes in systemically untreated early-stage triple-negative breast cancer

Operable triple-negative breast cancer (TNBC) has a higher risk of recurrence and death compared to other subtypes. Tumor size and nodal status are the primary clinical factors used to guide systemic treatment, while biomarkers of proliferation have not demonstrated value. Recent studies suggest that subsets of TNBC have a favorable prognosis, even without systemic therapy. We evaluated the association of fully automated mitotic spindle hotspot (AMSH) counts with recurrence-free (RFS) and overall survival (OS) in two separate cohorts of patients with early-stage TNBC who did not receive systemic therapy. AMSH counts were obtained from areas with the highest mitotic density in digitized whole slide images processed with a convolutional neural network trained to detect mitoses. In 140 patients from the Mayo Clinic TNBC cohort, AMSH counts were significantly associated with RFS and OS in a multivariable model controlling for nodal status, tumor size, and tumor-infiltrating lymphocytes (TILs) (p < 0.0001). For every 10-point increase in AMSH counts, there was a 16% increase in the risk of an RFS event (HR 1.16, 95% CI 1.08-1.25), and a 7% increase in the risk of death (HR 1.07, 95% CI 1.00-1.14). We corroborated these findings in a separate cohort of systemically untreated TNBC patients from Radboud UMC in the Netherlands. Our findings suggest that AMSH counts offer valuable prognostic information in patients with early-stage TNBC who did not receive systemic therapy, independent of tumor size, nodal status, and TILs. If further validated, AMSH counts could help inform future systemic therapy de-escalation strategies.

Pericyte Control of Gene Expression in the Blood-Brain Barrier Endothelium: Implications for Alzheimer's Disease

Background

A strong body of evidence suggests that cerebrovascular pathologies augment the onset and progression of Alzheimer's disease (AD). One distinctive aspect of this cerebrovascular dysfunction is the degeneration of brain pericytes-often overlooked supporting cells of blood-brain barrier endothelium.

Objective

The current study investigates the influence of pericytes on gene and protein expressions in the blood-brain barrier endothelium, which is expected to facilitate the identification of pathophysiological pathways that are triggered by pericyte loss and lead to blood-brain barrier dysfunction in AD.

Methods

Bioinformatics analysis was conducted on the RNA-Seq expression counts matrix (GSE144474), which compared solo-cultured human blood-brain barrier endothelial cells against endothelial cells co-cultured with human brain pericytes in a non-contact model. We constructed a similar cell culture model to verify protein expression using western blots.

Results

The insulin resistance and ferroptosis pathways were found to be enriched. Western blots of the insulin receptor and heme oxygenase expressions were consistent with those observed in RNA-Seq data. Additionally, we observed more than 5-fold upregulation of several genes associated with neuroprotection, including insulin-like growth factor 2 and brain-derived neurotrophic factor.

Conclusions

Results suggest that pericyte influence on blood-brain barrier endothelial gene expression confers protection from insulin resistance, iron accumulation, oxidative stress, and amyloid deposition. Since these are conditions associated with AD pathophysiology, they imply mechanisms by which pericyte degeneration could contribute to disease progression.

Endoxifen downregulates AKT phosphorylation through protein kinase C beta 1 inhibition in ERα+ breast cancer

Endoxifen, a secondary tamoxifen metabolite, is a potent antiestrogen exhibiting estrogen receptor alpha (ERα) binding at nanomolar concentrations. Phase I/II clinical trials identified clinical activity of Z-endoxifen (ENDX), in endocrine-refractory metastatic breast cancer as well as ERα+ solid tumors, raising the possibility that ENDX may have a second, ERα-independent, mechanism of action. An unbiased mass spectrometry approach revealed that ENDX concentrations achieved clinically with direct ENDX administration (5 µM), but not low concentrations observed during tamoxifen treatment (<0.1 µM), profoundly altered the phosphoproteome of the aromatase expressing MCF7AC1 cells with limited impact on the total proteome. Computational analysis revealed protein kinase C beta (PKCβ) and protein kinase B alpha or AKT1 as potential kinases responsible for mediating ENDX effects on protein phosphorylation. ENDX more potently inhibited PKCβ1 kinase activity compared to other PKC isoforms, and ENDX binding to PKCβ1 was confirmed using Surface Plasma Resonance. Under conditions that activated PKC/AKT signaling, ENDX induced PKCβ1 degradation, attenuated PKCβ1-activated AKTSer473 phosphorylation, diminished AKT substrate phosphorylation, and induced apoptosis. ENDX's effects on AKT were phenocopied by siRNA-mediated PKCβ1 knockdown or treatment with the pan-AKT inhibitor, MK-2206, while overexpression of constitutively active AKT diminished ENDX-induced apoptosis. These findings, which identify PKCβ1 as an ENDX target, indicate that PKCβ1/ENDX interactions suppress AKT signaling and induce apoptosis in breast cancer.

Estrogen-regulated miRs in bone enhance osteoblast differentiation and matrix mineralization

Estrogen signaling is critical for the development and maintenance of healthy bone, and age-related decline in estrogen levels contributes to the development of post-menopausal osteoporosis. Most bones consist of a dense cortical shell and an internal mesh-like network of trabecular bone that respond differently to internal and external cues such as hormonal signaling. To date, no study has assessed the transcriptomic differences that occur specifically in cortical and trabecular bone compartments in response to hormonal changes. To investigate this, we employed a mouse model of post-menopausal osteoporosis (ovariectomy, OVX) and estrogen replacement therapy (ERT). mRNA and miR sequencing revealed distinct transcriptomic profiles between cortical and trabecular bone in the setting of OVX and ERT. Seven miRs were identified as likely contributors to the observed estrogen-mediated mRNA expression changes. Of these, four miRs were prioritized for further study and decreased predicted target gene expression in bone cells, enhanced the expression of osteoblast differentiation markers, and altered the mineralization capacity of primary osteoblasts. As such, candidate miRs and miR mimics may have therapeutic relevance for bone loss resulting from estrogen depletion without the unwanted side effects of hormone replacement therapy and therefore represent novel therapeutic approaches to combat diseases of bone loss.

Integration of multiomics data shows down regulation of mismatch repair and tubulin pathways in triple-negative chemotherapy-resistant breast tumors

BACKGROUND

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype. Patients with TNBC are primarily treated with neoadjuvant chemotherapy (NAC). The response to NAC is prognostic, with reductions in overall survival and disease-free survival rates in those patients who do not achieve a pathological complete response (pCR). Based on this premise, we hypothesized that paired analysis of primary and residual TNBC tumors following NAC could identify unique biomarkers associated with post-NAC recurrence.

METHODS AND RESULTS

We investigated 24 samples from 12 non-LAR TNBC patients with paired pre- and post-NAC data, including four patients with recurrence shortly after surgery (< 24 months) and eight who remained recurrence-free (> 48 months). These tumors were collected from a prospective NAC breast cancer study (BEAUTY) conducted at the Mayo Clinic. Differential expression analysis of pre-NAC biopsies showed minimal gene expression differences between early recurrent and nonrecurrent TNBC tumors; however, post-NAC samples demonstrated significant alterations in expression patterns in response to intervention. Topological-level differences associated with early recurrence were implicated in 251 gene sets, and an independent assessment of microarray gene expression data from the 9 paired non-LAR samples available in the NAC I-SPY1 trial confirmed 56 gene sets. Within these 56 gene sets, 113 genes were observed to be differentially expressed in the I-SPY1 and BEAUTY post-NAC studies. An independent (n = 392) breast cancer dataset with relapse-free survival (RFS) data was used to refine our gene list to a 17-gene signature. A threefold cross-validation analysis of the gene signature with the combined BEAUTY and I-SPY1 data yielded an average AUC of 0.88 for six machine-learning models. Due to the limited number of studies with pre- and post-NAC TNBC tumor data, further validation of the signature is needed.

CONCLUSION

Analysis of multiomics data from post-NAC TNBC chemoresistant tumors showed down regulation of mismatch repair and tubulin pathways. Additionally, we identified a 17-gene signature in TNBC associated with post-NAC recurrence enriched with down-regulated immune genes.

Distinct spatial immune microlandscapes are independently associated with outcomes in triple-negative breast cancer

The utility of spatial immunobiomarker quantitation in prognostication and therapeutic prediction is actively being investigated in triple-negative breast cancer (TNBC). Here, with high-plex quantitative digital spatial profiling, we map and quantitate intraepithelial and adjacent stromal tumor immune protein microenvironments in systemic treatment-naïve (female only) TNBC to assess the spatial context in immunobiomarker-based prediction of outcome. Immune protein profiles of CD45-rich and CD68-rich stromal microenvironments differ significantly. While they typically mirror adjacent, intraepithelial microenvironments, this is not uniformly true. In two TNBC cohorts, intraepithelial CD40 or HLA-DR enrichment associates with better outcomes, independently of stromal immune protein profiles or stromal TILs and other established prognostic variables. In contrast, intraepithelial or stromal microenvironment enrichment with IDO1 associates with improved survival irrespective of its spatial location. Antigen-presenting and T-cell activation states are inferred from eigenprotein scores. Such scores within the intraepithelial compartment interact with PD-L1 and IDO1 in ways that suggest prognostic and/or therapeutic potential. This characterization of the intrinsic spatial immunobiology of treatment-naïve TNBC highlights the importance of spatial microenvironments for biomarker quantitation to resolve intrinsic prognostic and predictive immune features and ultimately inform therapeutic strategies for clinically actionable immune biomarkers.

Bayesian Machine Learning Enables Identification of Transcriptional Network Disruptions Associated with Drug-Resistant Prostate Cancer

Survival rates of patients with metastatic castration-resistant prostate cancer (mCRPC) are low due to lack of response or acquired resistance to available therapies, such as abiraterone (Abi). A better understanding of the underlying molecular mechanisms is needed to identify effective targets to overcome resistance. Given the complexity of the transcriptional dynamics in cells, differential gene expression analysis of bulk transcriptomics data cannot provide sufficient detailed insights into resistance mechanisms. Incorporating network structures could overcome this limitation to provide a global and functional perspective of Abi resistance in mCRPC. Here, we developed TraRe, a computational method using sparse Bayesian models to examine phenotypically driven transcriptional mechanistic differences at three distinct levels: transcriptional networks, specific regulons, and individual transcription factors (TF). TraRe was applied to transcriptomic data from 46 patients with mCRPC with Abi-response clinical data and uncovered abrogated immune response transcriptional modules that showed strong differential regulation in Abi-responsive compared with Abi-resistant patients. These modules were replicated in an independent mCRPC study. Furthermore, key rewiring predictions and their associated TFs were experimentally validated in two prostate cancer cell lines with different Abi-resistance features. Among them, ELK3, MXD1, and MYB played a differential role in cell survival in Abi-sensitive and Abi-resistant cells. Moreover, ELK3 regulated cell migration capacity, which could have a direct impact on mCRPC. Collectively, these findings shed light on the underlying transcriptional mechanisms driving Abi response, demonstrating that TraRe is a promising tool for generating novel hypotheses based on identified transcriptional network disruptions.

SIGNIFICANCE

The computational method TraRe built on Bayesian machine learning models for investigating transcriptional network structures shows that disruption of ELK3, MXD1, and MYB signaling cascades impacts abiraterone resistance in prostate cancer.

Abstract 1689: The role of circular RNAs in triple negative breast cancer and chemotherapy resistance

Background: Breast cancer is the most frequently diagnosed cancer and the leading cause of cancer-related death among women worldwide. Triple negative breast cancer (TNBC) represents 15-20% of all breast cancers and is an aggressive subtype. Chemotherapy-based treatments remain the standard of care for TNBC. Unfortunately, chemotherapy resistance is common, and for these patients, outcomes are poor and alternative treatment strategies remain an unmet need. circRNAs are a newly identified class of noncoding RNA molecules with covalently closed circular structures. An increasing number of recent studies including ours have indicated that circRNAs play crucial roles in regulating tumor development and chemoresistance. However, the role of circRNAs in the process of chemotherapy resistance and TNBC progression is not clear.

Materials and Methods: As a first step towards identifying circRNAs that participate in the development of chemoresistance in TNBC cells, and to determine if targeting such circRNAs is a novel and efficacious therapeutic strategy, doxorubicin-resistant (Doxo-R), paclitaxel-resistant (PTX-R), and double-resistant (DP-R) cell lines were generated from MDA-MB-231. Human circRNA microarrays were utilized to profile the expression of approximately 14,000 known circRNAs in normal breast tissue, matched patient-derived xenografts (PDX) generated prior to and following neoadjuvant chemotherapy (NAC), and TNBC chemosensitive and chemoresistant cell lines. Top hits were validated using RT-PCR.

Results: circRNA microarray profiling identified 429 and 310 transcripts differentially expressed in doxorubicin and paclitaxel resistant cells, respectively, compared to parental chemosensitive cell lines (|FC| ≥ 1.5; p value &lt; 0.05). In comparison to pre-NAC derived xenografts, 1,396 circRNAs were dysregulated among post-NAC PDX models. Further, three circRNAs (hsa_circ_001388, hsa_circ_104652, and hsa_circ_061260) were upregulated in Doxo-R, PTX-R, and post-NAC PDX samples compared to their respective controls. Among these three circRNAs, hsa_circ_001388 (also known as circNSD2) was the only transcript also predicted to be translated into a novel and uncharacterized protein given the presence of a high confidence translation initiation site and IRES sequence. Ongoing studies are aimed at determining the role of circNSD2 protein in breast cancer carcinogenesis, progression, and response to standard of care chemotherapeutics and the mechanistic process by which this protein functions. Additionally, the efficacy of specifically targeting this circRNA as a novel therapeutic approach is being explored.

Conclusions: Increasing knowledge of the important functions of circRNAs underlying drug resistance will provide new opportunities for developing efficacious therapeutic strategies and prognostic/predictive biomarkers for TNBC.

Citation Format: Xiyin Wang, Michael J. Emch, Xiaojia Tang, Jia Yu, Krishna R. Kalari, Liewei Wang, Matthew P. Goetz, John R. Hawse. The role of circular RNAs in triple negative breast cancer and chemotherapy resistance [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 1689.

Dynamic assessment of serum chromogranin A and treatment response with abiraterone acetate in metastatic castration-resistant prostate cancer

OBJECTIVE

Elevated serum chromogranin A (CGA) is associated with intrinsic or treatment-related neuroendocrine differentiation (NED) in men with metastatic castration-resistant prostate cancer (mCRPC). Fluctuations in serum CGA during treatment of mCRPC have had conflicting results. We analyzed the impact of (i) rising serum CGA and (ii) baseline CGA/PSA ratio during treatment to identify associations with abiraterone acetate (AA) therapy.

METHODS

Between June 2013 and August 2015, 92 men with mCRPC were enrolled in a prospective trial with uniform serum CGA processing performed before initiating abiraterone acetate/prednisone (AA/P) and serially after 12 weeks of AA/P treatments. Serum CGA was measured using a homogenous automated immunofluorescent assay. Patients receiving proton pump inhibitors or with abnormal renal function were excluded due to possible false elevations of serum CGA (n = 21 excluded), therefore 71 patients were analyzed. All patients underwent a composite response assessment at 12-weeks. Kaplan-Meier estimates and Cox Regression models were used to calculate the association with time-to-treatment failure analyses and overall survival.

RESULTS

An increase in chromogranin was associated with a lower risk of treatment failure (hazard ratio [HR]: 0.52, p = 0.0181). The median CGA/PSA ratio was 7.8 (2.6-16.0) and an elevated pretreatment CGA/PSA ratio above the median was associated with a lower risk of treatment failure (HR: 0.54 p value = 0.0185). An increase in CGA was not found to be associated with OS (HR: 0.71, 95% CI: 0.42-1.21, p = 0.207). An elevated baseline CGA/PSA ratio was not associated with OS (HR: 0.62, 95% CI: 0.37-1.03, p = 0.062). An increase in PSA after 12 weeks of treatment was associated with an increased risk of treatment failure (HR: 4.14, CI: 2.21-7.73, p = < 0.0001) and worse OS (HR: 2.93, CI: 1.57-4.45, p = < 0.0001).

CONCLUSIONS

We show that an increasing chromogranin on AA/P and an elevated baseline CGA/PSA in patients with mCRPC were associated with a favorable response to AA/P with no changes in survival. There may be limited clinical utility in serum CGA testing to evaluate for lethal NED as AA/P did not induce lethal NED in this cohort. This highlights that not all patients with an increasing CGA have a worse OS.

Abstract P2-11-34: Mitotic spindle hotspot counting using deep learning networks is highly associated with clinical outcomes in patients with early-stage triple-negative breast cancer who did not receive systemic therapy

Background: Triple-negative breast cancers (TNBC) exhibit high rates of recurrence and mortality. However, recent studies suggest that a subset of patients (pts) with early-stage TNBC enriched in tumor-infiltrating lymphocytes (TILs) have excellent clinical outcomes even in the absence of systemic therapy. Additional histological biomarkers that could identify pts for future systemic therapy escalation/de-escalation strategies are of great interest. TNBC are frequently highly proliferative with abundant mitoses. However, classic markers of proliferation (manual mitosis counting and Ki-67) appear to offer no prognostic value. Here, we evaluated the prognostic effects of automated mitotic spindle hotspot (AMSH) counting on RFS in independent cohorts of systemically untreated early-stage TNBC.

Methods: AMSH counting was conducted with a state-of-the-art deep learning algorithm trained on the detection of mitoses within 2 mm2 areas with the highest mitotic density (i.e. hotspots) in digital H&E images. Details of the development, training and validation of the algorithm were published previously [1] in a cohort of unselected TNBC. We obtained AMSH counts in a centrally confirmed TNBC cohort from Mayo Clinic [2] and focused our analysis on pts who received locoregional therapy but no systemic therapy. Using a fractional polynomial analysis with a multivariable proportional hazards regression model, we confirmed the assumption of linearity in the log hazard for the continuous variable AMSH and evaluated whether AMSH counts were prognostic of RFS. We corroborated our findings in an independent cohort of systemically untreated TNBC pts from the Radboud University Medical Center in the Netherlands (Radboud Cohort). Results are reported at a median follow-up of 8.1 and 6.7 years for the Mayo and Netherlands cohorts, respectively.

Results: Among 182 pts with who did not receive systemic therapy in the Mayo Cohort, 140 (77%) with available AMSH counts were included. The mean age was 61 (range: 31-94), 71% were postmenopausal, 67% had tumors ≤ 2cm, and 83% were node-negative. As expected, most tumors were Nottingham grade 3 (84%) and had a high Ki-67 proliferation index (54% with Ki-67 &gt;30%). Most tumors (73%) had stromal TILs ≤ 30%. The median AMSH count was 18 (IQR: 8, 42). AMSH counts were linearly associated with grade and tumor size, with the proportion of pts with grade 3 tumors and size &gt; 2 cm increasing as the AMSH counts increased (p=0.007 and p=0.059, respectively). In a multivariate model controlling for nodal status, tumor size, and stromal TILs, AMSH counts were independently associated with RFS (p&lt; 0.0001). For every 10-point increase in the AMSH count, we observed a 17% increase in the risk of experiencing an RFS event (HR 1.17, 95% CI 1.08-1.26). We corroborated our findings in the Radboud Cohort (n=126). The mean age was 68 (range: 40-96), and 81% were node-negative. While the median AMSH count was 36 (IQR: 16-63), higher than in the Mayo Cohort (p=0.004), the prognostic impact was similar, with a significant association between AMSH count and RFS (p=0.028) in a multivariate model corrected for nodal status, tumor size, and stromal TILs. For every 10-point increase in the AMSH count in the Netherlands cohort, we observed a 9% increase in the risk of experiencing an RFS event (HR 1.09, 95% CI 1.01-1.17). RFS rates according to AMSH counts for both cohorts are shown in the Table.

Conclusions: AMSH counting is a new proliferation biomarker that provides prognostic value independent of nodal status, tumor size, and stromal TILs in systemically untreated early-stage TNBC. Plans are underway to evaluate AMSH counts in additional cohorts of systemically untreated TNBC, and in other disease settings such as prior to neoadjuvant systemic therapy. If validated, this biomarker should be prospectively evaluated as a potential selection biomarker in clinical trials of systemic therapy de-escalation.

References:

1. PMID: 29994086

2. PMID: 28913760

Table RFS according to AMSH counts in the Mayo and Radboud Cohorts

Citation Format: Roberto A. Leon-Ferre, Jodi M. Carter, David Zahrieh, Jason P. Sinnwell, Roberto Salgado, Vera Suman, David Hillman, Judy C. Boughey, Krishna R. Kalari, Fergus J. Couch, James N. Ingle, Maschenka Balkenkohl, Francesco Ciompi, Jeroen van der Laak, Matthew P. Goetz. Mitotic spindle hotspot counting using deep learning networks is highly associated with clinical outcomes in patients with early-stage triple-negative breast cancer who did not receive systemic therapy [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 P2-11-34.

Molecular Profile Changes in Patients with Castrate-Resistant Prostate Cancer Pre- and Post-Abiraterone/Prednisone Treatment

We identified resistance mechanisms to abiraterone acetate/prednisone (AA/P) in patients with metastatic castration-resistant prostate cancer (mCRPC) in the Prostate Cancer Medically Optimized Genome-Enhanced Therapy (PROMOTE) study.

We analyzed whole-exome sequencing (WES) and RNA-sequencing data from 83 patients with metastatic biopsies before (V1) and after 12 weeks of AA/P treatment (V2). Resistance was determined by time to treatment change (TTTC).

At V2, 18 and 11 of 58 patients had either short-term (median 3.6 months; range 1.4-4.5) or long-term (median 29 months; range 23.5-41.7) responses, respectively. Nonresponders had low expression of TGFBR3 and increased activation of the Wnt pathway, cell cycle, upregulation of AR variants, both pre- and posttreatment, with further deletion of AR inhibitor CDK11B posttreatment. Deletion of androgen processing genes, HSD17B11, CYP19A1 were observed in nonresponders posttreatment. Genes involved in cell cycle, DNA repair, Wnt-signaling, and Aurora kinase pathways were differentially expressed between the responder and non-responder at V2. Activation of Wnt signaling in nonresponder and deactivation of MYC or its target genes in responders was detected via SCN loss, somatic mutations, and transcriptomics. Upregulation of genes in the AURKA pathway are consistent with the activation of MYC regulated genes in nonresponders. Several genes in the AKT1 axis had increased mutation rate in nonresponders. We also found evidence of resistance via PDCD1 overexpression in responders.

IMPLICATIONS

Finally, we identified candidates drugs to reverse AA/P resistance: topoisomerase inhibitors and drugs targeting the cell cycle via the MYC/AURKA/AURKB/TOP2A and/or PI3K_AKT_MTOR pathways.

Identification of Src Family Kinases as Potential Therapeutic Targets for Chemotherapy-Resistant Triple Negative Breast Cancer

Neoadjuvant chemotherapy (NAC) remains the cornerstone of the treatment for triple negative breast cancer (TNBC), with the goal of complete eradication of disease. However, for patients with residual disease after NAC, recurrence and mortality rates are high and the identification of novel therapeutic targets is urgently needed. We quantified tyrosine phosphorylation (pTyr)-mediated signaling networks in chemotherapy sensitive (CS) and resistant (CR) TNBC patient-derived xenografts (PDX), to gain novel therapeutic insights. The antitumor activity of SFK inhibition was examined in vivo. Treated tumors were further subjected to phosphoproteomic and RNAseq analysis, to identify the mechanism of actions of the drug. We identified Src Family Kinases (SFKs) as potential therapeutic targets in CR TNBC PDXs. Treatment with dasatinib, an FDA approved SFK inhibitor, led to inhibition of tumor growth in vivo. Further analysis of post-treatment PDXs revealed multiple mechanisms of actions of the drug, confirming the multi-target inhibition of dasatinib. Analysis of pTyr in tumor specimens suggested a low prevalence of SFK-driven tumors, which may provide insight into prior clinical trial results demonstrating a lack of dasatinib antitumor activity in unselected breast cancer patients. Taken together, these results underscore the importance of pTyr characterization of tumors, in identifying new targets, as well as stratifying patients based on their activated signaling networks for therapeutic options. Our data provide a strong rationale for studying SFK inhibitors in biomarker-selected SFK-driven TNBC.

Biomarkers for Predicting Abiraterone Treatment Outcome and Selecting Alternative Therapies in Castration-Resistant Prostate Cancer

Approximately one‐third of patients with metastatic castration‐resistant prostate cancer (CRPC) exhibited primary abiraterone resistance. To identify alternative treatment for abiraterone nonresponders, we performed drug discovery analyses using the L1000 database using differentially expressed genes identified in tumor biopsies and patient‐derived xenograft (PDX) tumors between abiraterone responders and nonresponders enrolled in PROMOTE trial. This approach identified 3 drugs, including topoisomerase II (TOP2) inhibitor mitoxantrone, CDK4/6 inhibitor palbociclib, and pan‐CDK inhibitor PHA‐793887. These drugs significantly suppressed the growth of abiraterone‐resistant cell lines and PDX models. Moreover, we identified 11 genes targeted by all 3 drugs that were associated with worse outcomes in both the PROMOTE and Stand Up To Cancer cohorts. This 11‐gene panel might also function as biomarkers to select the 3 alternative therapies for this subgroup of patients with CRPC, warranting further clinical investigation.

Toward Individualized Prediction of Response to Methotrexate in Early Rheumatoid Arthritis: A Pharmacogenomics-Driven Machine Learning Approach

OBJECTIVE

To test the ability of machine learning (ML) approaches with clinical and genomic biomarkers to predict methotrexate treatment response in patients with early rheumatoid arthritis (RA).

METHODS

Demographic, clinical, and genomic data from 643 patients of European ancestry with early RA (mean age 54 years; 70% female) subdivided into a training (n = 336) and validation cohort (n = 307) were used. The genomic data comprised 160 single-nucleotide polymorphisms (SNPs) previously associated with RA or methotrexate metabolism. Response to methotrexate monotherapy was defined as good or moderate by the European Alliance of Associations for Rheumatology (EULAR) response criteria at the 3-month follow-up. Supervised ML methods were trained with 5 repeats and 10-fold cross-validation using the training cohort. Prediction performance was validated in the independent validation cohort.

RESULTS

Supervised ML methods combining age, sex, smoking, rheumatoid factor, baseline Disease Activity Score in 28 joints (DAS28) scores and 160 SNPs predicted EULAR response at 3 months with the area under the receiver operating curve of 0.84 (P = 0.05) in the training cohort and achieved a prediction accuracy of 76% (P = 0.05) in the validation cohort (sensitivity 72%, specificity 77%). Intergenic SNPs rs12446816, rs13385025, rs113798271, and ATIC (rs2372536) had variable importance above 60.0 and along with baseline DAS28 scores were among the top predictors of methotrexate response.

CONCLUSION

Pharmacogenomic biomarkers combined with baseline DAS28 scores can be useful in predicting response to methotrexate in patients with early RA. Applying ML to predict treatment response holds promise for guiding effective RA treatment choices, including timely escalation of RA therapies.

Implementation of preemptive DNA sequence-based pharmacogenomics testing across a large academic medical center: The Mayo-Baylor RIGHT 10K Study

PURPOSE

The Mayo-Baylor RIGHT 10K Study enabled preemptive, sequence-based pharmacogenomics (PGx)-driven drug prescribing practices in routine clinical care within a large cohort. We also generated the tools and resources necessary for clinical PGx implementation and identified challenges that need to be overcome. Furthermore, we measured the frequency of both common genetic variation for which clinical guidelines already exist and rare variation that could be detected by DNA sequencing, rather than genotyping.

METHODS

Targeted oligonucleotide-capture sequencing of 77 pharmacogenes was performed using DNA from 10,077 consented Mayo Clinic Biobank volunteers. The resulting predicted drug response-related phenotypes for 13 genes, including CYP2D6 and HLA, affecting 21 drug-gene pairs, were deposited preemptively in the Mayo electronic health record.

RESULTS

For the 13 pharmacogenes of interest, the genomes of 79% of participants carried clinically actionable variants in 3 or more genes, and DNA sequencing identified an average of 3.3 additional conservatively predicted deleterious variants that would not have been evident using genotyping.

CONCLUSION

Implementation of preemptive rather than reactive and sequence-based rather than genotype-based PGx prescribing revealed nearly universal patient applicability and required integrated institution-wide resources to fully realize individualized drug therapy and to show more efficient use of health care resources.

Estrogen receptor beta repurposes EZH2 to suppress oncogenic NFκB/p65 signaling in triple negative breast cancer

Triple Negative Breast Cancer (TNBC) accounts for 15-20% of all breast cancer cases, yet is responsible for a disproportionately high percentage of breast cancer mortalities. Thus, there is an urgent need to identify novel biomarkers and therapeutic targets based on the molecular events driving TNBC pathobiology. Estrogen receptor beta (ERβ) is known to elicit anti-cancer effects in TNBC, however its mechanisms of action remain elusive. Here, we report the expression profiles of ERβ and its association with clinicopathological features and patient outcomes in the largest cohort of TNBC to date. In this cohort, ERβ was expressed in approximately 18% of TNBCs, and expression of ERβ was associated with favorable clinicopathological features, but correlated with different overall survival outcomes according to menopausal status. Mechanistically, ERβ formed a co-repressor complex involving enhancer of zeste homologue 2/polycomb repressive complex 2 (EZH2/PRC2) that functioned to suppress oncogenic NFκB/RELA (p65) activity. Importantly, p65 was shown to be required for formation of this complex and for ERβ-mediated suppression of TNBC. Our findings indicate that ERβ+ tumors exhibit different characteristics compared to ERβ- tumors and demonstrate that ERβ functions as a molecular switch for EZH2, repurposing it for tumor suppressive activities and repression of oncogenic p65 signaling.

Abstract P5-01-07: Patient-derived xenografts (PDXs) generated from hormone receptor-positive breast cancer (BC) before and after cyclin-dependent kinase 4 and 6 inhibitors (CDK4/6i) treatment: Initial findings from the PROMISE study

Background: Hormone receptor-positive (HR+) breast cancer (BC) accounts for ~60-70% of all BC cases. CDK4/6i combined with endocrine therapy (ET) is the standard of care for first- and second-line therapy of patients with HR+ metastatic breast cancer (MBC) based on substantial improvements in progression free survival (PFS) and overall survival compared to ET alone. Unfortunately, most patients (pts) develop resistance to CDK4/6i. The mechanisms underlying CDK4/6i resistance have not been fully elucidated. Robust in vivo CDK4/6i resistant animal models from pts with de novo and acquired CDK 4/6i resistance are lacking. Patient-derived xenografts (PDXs) are important preclinical tools in oncology to advance our understanding of cancer biology and to evaluate drug response phenotypes. Unfortunately, the establishment of sustainable PDXs from pts with HR+ BC is much more difficult than from HR-negative BC. In order to obtain clinically relevant models that can be utilized to better understand mechanisms of primary and acquired resistance to ET and CDK4/6i, we report the development of PDXs from a prospective study in which percutaneous biopsies were obtained from women receiving palbociclib for the treatment of HR+ advanced or MBC (PROMISE; NCT03281902). Methods: PROMISE (Prospective Study to Evaluate the Role of Tumor Sequencing in Women Receiving Palbociclib for Advanced Hormone Receptor (HR)-Positive Breast Cancer) is a prospective trial which enrolled pts with HR+, HER2-negative MBC who planned to start treatment with palbociclib and letrozole (1st line [FL]) or palbociclib and fulvestrant (2nd line [SL]). For pts enrolled at Mayo Clinic Rochester (MCR), serial tumor biopsies were obtained prior to the start of treatment and at disease progression to establish PDXs. Because baseline take rate was low, we amended the protocol to attempt PDX generation after the 2nd cycle (C2) of therapy. Tissues were received within an hour of excision and were orthotopically implanted in the mammary fat pad of 6~8-week-old nonobese diabetic (NOD)/SCID/IL-2g-receptor null (NSG) mice. Tumor growth was defined as tumor volume ≥ 50 mm3. Tumor take rate was defined as percent of patients with development of at least one stably transplantable (passed at least for two generations) xenograft that was pathologically confirmed as BC. Results: Of the 50 pts enrolled at MCR, 37 were FL pts (2 ineligible) and 13 were SL pts (1 ineligible). Fresh tumor samples obtained pre-treatment from 43 eligible pts (32 FL), after C2 from 17 eligible pts (13 FL), and at disease progression from 10 eligible pts (5 FL) were implanted. Tumor growth was seen in at least one mouse from a pre-treatment specimen in 1 (3.1%) FL pt (confirmed BC) and 4 (36.4%) SL pts (2 confirmed BC); from end of C2 specimens in 2 (15.4%) FL pts and 2 (50%) SL pts (1 confirmed BC); and at progression in 1 (20%) FL pt and 2 (40%) SL pts (1 confirmed BC). The take rate from pre-treatment samples among FL pts was 1/32 (3.1%; 95% CI: 0.08-16.2%). Take rates in other groups are pending histological review. Discussion: In pts receiving palbociclib for HR+ MBC, the take rate using a pre-treatment percutaneous biopsy was low (3.1%) in the FL setting. However, higher growth rates were seen in the SL setting. The observation of higher growth rates at the end of C2 and at progression may relate to the activation of G1/S that occurs with withdrawal of palbociclib, and should be further studied. PDX remains a valuable tool for the evaluation of tumor biology and the development of new therapeutics, especially for those models established from pts with palbociclib and ET resistance. Proteogenomic and sequencing studies are ongoing to fully characterize these PDXs.

Citation Format: Jia Yu, Vera J. Suman, Jun He, Jason P. Sinnwell, Ann M. Moyer, Bo Qin, Yayun Gu, Huan Zhang, Jodi M. Carter, Krishna R. Kalari, Brendan P. McMenomy, Minetta C. Liu, Tufia C. Haddad, Kathryn J. Ruddy, Fergus J. Couch, Alvaro Moreno-Aspitia, Donald W. Northfelt, Richard Weinshilboum, Ciara O'Sullivan, Matthew P. Goetz, Liewei Wang. Patient-derived xenografts (PDXs) generated from hormone receptor-positive breast cancer (BC) before and after cyclin-dependent kinase 4 and 6 inhibitors (CDK4/6i) treatment: Initial findings from the PROMISE study [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-01-07.

Abstract OT2-19-05: A phase I/II trial of abemaciclib and T-DM1 in women and men with HER2-positive advanced or metastatic breast cancer that has progressed on treatment with a taxane, trastuzumab and pertuzumab (THP) (ACCRU-BR-1801)

Background: Although major advances have been made in the treatment of HER2+ metastatic breast cancer (MBC), the goal of care remains largely palliative, therefore better treatments are needed. Given encouraging preclinical and clinical data, the combination of cyclin dependent 4/6 kinase inhibitors and HER2-directed therapy is further being evaluated in this trial. Trial Design: In this phase I/II multicenter trial, we will determine the maximum tolerated dose (MTD) of abemaciclib (Abem) combined with T-DM1. Three Abem dose levels will be examined, 50 mg, 100 mg and 150 mg. The phase 2 portion of this trial will examine whether PFS is increased with addition of Abem to T-DM1 in two pt cohorts - those with ER+ HER2+ MBC and those with ER-HER2+ MBC. For phase 2, a pre-registration biopsy is required to confirm ER and HER2-status and to determine levels of tumor infiltrating lymphocytes; vimentin (an epithelial-mesenchymal transition marker); and CD8 and FOXP3 expression. Blood samples will be collected pretreatment , at 6 weeks, and at progression for all pts. Eligibility Criteria: Phase I&II: All pts must have HER2+ MBC per ASCO-CAP guidelines and prior treatment with a taxane, trastuzumab and pertuzumab. For the phase Iportion, pts can have measurable or non-measurable disease with no restriction on the number of prior lines of therapy. In the phase II portion, pts must have measurable disease with ≤1-2 prior lines of chemotherapy alone, ≤1 HER2-directed therapy alone, and/or chemotherapy with HER2-directed therapies. There is no limit on prior endocrine therapy. Specific Aims: The primary objective of the phase II trial is to assess whether addition of Abem to T-DM1 increases PFS in one or both patient cohorts. Secondary objectives include an assessment of toxicity, objective response rates and overall survival. Correlative studies will assess the association between baseline TIL levels, vimentin expression, and CD8/FOXP3 expression with PFS. Changes in peripheral blood mononuclear cells, CTCs, ctDNA and serum thymidine kinase 1 during the course of treatment will be examined to determine if there is a link with PFS outcomes overall and separately for each cohort. Pharmacogenomic studies will determine if pts with the FCGR3A-158 polymorphism derive less benefit from T-DM1 and have inferior PFS outcomes compared with pts who do not have this polymorphism. Statistical Methods:Phase I: Standard 3+3 design, with dose limiting toxicities as per protocol. Phase II: For each pt cohort, a stratified randomization scheme will be used to assign pts to treatment with liver mets as a stratification factor.

For each pt cohort, a stratified log rank test will be used to assess whether PFS is increased with the addition of Abem to T-DM1. A non-binding futility analysis will be applied in each cohort after 58 events in the ER+ HER2+ MBC study cohort and 48 events in the ER- HER2+ MBC study cohort. Present Accrual: 0; target accrual: minimal 120 pts., maximal 140 pts

CohortOne-sidedalphaPowerAccrual Period (accrual rate)Follow-up after close of enrollmentPFS with T-DM1PFS with abema and T-DM1Number of eligible patientsER+/HER2+0.100.912 months (5-6 pts per month)12 months12 weeks24 weeks64 (32 per arm)ER-/HER2+0.100.8512 months (3-4 pts per month)12 months6 weeks12 weeks50 (25 pts per arm)

Citation Format: Ciara C O Sullivan, Jun He, Vera J Suman, Krishna R Kalari, Roberto A Leon-Ferre, Jose C Villasboas-Bisneto, Pratima Chalasani, Demet Gokalp Yasar, Daniel M Anderson, Philip J Stella, Anthony J Jaslowski, Susan H Tannenbaum, Angela Saverimuthu, Donald Northfelt, Alvaro Moreno-Aspitia, Jodi M Carter, Minetta C Liu, Liewei Wang, Zhenkun Lou, Matthew P Goetz. A phase I/II trial of abemaciclib and T-DM1 in women and men with HER2-positive advanced or metastatic breast cancer that has progressed on treatment with a taxane, trastuzumab and pertuzumab (THP) (ACCRU-BR-1801) [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 OT2-19-05.

Abstract P3-08-02: The frequency and somatic mutation landscape of Fibroblast growth factor receptor (FGFR) alterations in breast cancer

Background: FGFR dysregulation is observed in multiple cancers and targeting FGFR is an emerging therapeutic strategy with FDA approved treatments in bladder and cholangiocarcinoma. Here we examined the prevalence of FGFR mutations, fusions, and high-level amplifications in breast cancer, stratified by receptor subtype and local/metastatic status, in both Foundation Medicine (FM) and institutional Mayo Clinic (MC) cohorts. Methods: For the FM cohort, comprehensive genomic profiling (CGP) examining at least 324 genes for all classes of alterations, including FGFR1-4 was carried out for 32,048 breast cancers during the course of routine clinical care in a Clinical Laboratory Improvement Amendments (CLIA)-certified lab (Foundation Medicine Inc., Cambridge, MA, USA). Tumor mutational burden (TMB) was determined on up to 1.1 Mb, microsatellite instability high (MSI-High) was determined on up to 114 loci and predicted ancestry from &gt;10,000 SNPs. Estrogen receptor (ER) and HER2 status were available for a subset of FM samples. Additionally, 131 patients with metastatic breast cancer from a subset of patients at three Mayo Clinic sites (MC cohort) with clinical characteristics and cancer-panel DNA sequencing data from a CLIA-certified lab (Tempus, Chicago, IL) were included. Results: In the FM cohort, the prevalence of FGFR1-4 high-level amplification (CN≥10) was 10.1%, while mutations (1.5%) and fusions (0.72%) were rare. Most amplifications occurred in FGFR1 (9.2%); most fusions and mutations occurred in FGFR2 (0.46%, 0.77%). FGFR alteration prevalence was highest in ER+/HER2- subtype (14.4%) and lowest in HER2+ disease (7.7%). FGFR alterations were more common in IDC (11.7%) than ILC (7.7%), p&lt;3E-08. FGFR alterations were more prevalent in the metastatic setting relative to breast-biopsied disease (13.6% v 10.1%; OR = 1.4; p=2E-17), especially in the HER2+ (OR =1.9, p=0.004) and ER-/HER2- (OR = 1.9, p = 0.05) disease; no enrichment was seen in the ER+/HER2- metastases (OR =1.0, p = 1). FGFR amplifications were observed at a higher prevalence in patients with predicted East Asian ancestry, relative to patients with European ancestry (12.1% v 10.0%; p = 0.03). Overall, the most common activating mutations in FGFR were FGFR2 N549K (n=85), FGFR1 N546K (n=78), FGFR4 V510M (n=28), FGFR2 K659E (n=28), FGFR4 V510L (n=20), and FGFR2 Y375C (n=15). The most common recurrent fusions were FGFR3:TACC3 (n=36), FGFR2:TACC2 (n=17), FGFR1:TACC1 (n=9), FGFR1:BAG4 (n=6), and FGFR2:ATE1 (n=5). In patients with FGFR amplifications, the most frequently co-occurring alterations were ZNF703 (78.4%), TP53 (51.5%), CCND1 (36.1%), FGF3/4/19 (32.9 - 34.4%), PIK3CA (30.7%), MYC (29.6%), ESR1 (17.2%), EMSY (16.3%), and PTEN (10.6%). Significant co-occurrence was observed for a number of genes including FGF3/4/19, CDK4, and CDK8 (all OR&gt;2, p&lt;1E-07); mutual exclusivity was observed with PIK3R1, BRCA1, and BRCA2 (all OR &lt;0.5, p&lt;4E-13), among other genes. In the 131 metastatic tumors from MC, the prevalence of FGFR1-4 high-level amplifications was 19.8% [FGFR1 (12.4%), FGFR2 (7.4%), and FGFR3 (0.8%)]. The prevalence of high-level FGFR amplifications did not differ by clinical subtypes: HR-/HER2- (7/31), HR+/HER2- (15/79), and HER2+ (2/11), p=0.68. Conclusions: High-level FGFR amplifications are observed in &gt;11% of breast cancers, especially the ER+/HER2- subtype, while mutations/fusions are rare. These data support the ongoing studies evaluating targeted therapies for FGFR amplified ER + breast cancer. Correlations with clinical information (MC cohort) and associations with actionable alterations are ongoing and may inform potential combination strategies.

Citation Format: Karthik V Giridhar, Ethan S Sokol, Peter T Vedell, Jason P Sinnwell, Aakash Desai, Tufia C Haddad, Ciara C O’Sullivan, Roberto A Leon-Ferre, Siddhartha Yadav, Kostandinos Sideras, Brenda Ernst, Minetta C Liu, Abe Eyman Casey, Xiaojia Tang, Zoe Fleischmann, Karthikeyan Murugesan, Krishna R Kalari, Matthew P Goetz. The frequency and somatic mutation landscape of Fibroblast growth factor receptor (FGFR) alterations in breast cancer [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 P3-08-02.

Abstract P1-04-01: Digital spatial profiling of immune-related proteins in luminal androgen receptor (LAR) vs non-LAR triple-negative breast cancer (TNBC)

Background: The importance of the antitumor immune response in TNBC is well established. TNBC with higher TILs are less likely to recur and more responsive to systemic therapy. Likewise, PD-L1+ TNBC are more likely to benefit from chemoimmunotherapy. However, TNBC is highly heterogeneous. Of the TNBC molecular subtypes, LAR TNBC is less sensitive to systemic therapy, has lower TILs and lower rates of PD-L1 positivity. The role of other immune related proteins in LAR TNBC is not well established. Here, we evaluated differentially expressed (DE) immune related proteins in the stromal and intratumoral compartments of LAR vs non-LAR TNBC tumors. Methods: We used the Nanostring GeoMX DSP platform to quantitate 58 proteins within spatially distinct intraepithelial, cytokeratin (CK)-positive tumor segments and adjacent CK-negative/nuclei-positive stromal segments in 248 TNBC tumors included in a tissue microarray generated from a cohort of pts with centrally confirmed TNBC who underwent breast surgery without prior neoadjuvant therapy. A subset (n=111) underwent bulk tumor RNA sequencing and were classified as LAR or non-LAR TNBC. DE proteins were identified using a negative binomial generalized linear model (SNR&gt;2, p&lt;0.05). A targeted set of DE proteins was dichotomized at the 80th percentile. Results: Of 111 TNBC tumors, 17 (15%) were LAR and 94 (85%) non-LAR. Compared to non-LAR TNBC, pts with LAR TNBC were older (age ≥50: 82% vs 52%, p&lt;0.01), with tumors that were more often of apocrine histology (35% vs 0%, p &lt;0.01), grade 1-2 (24% vs 1%, p&lt;0.01), and had lower Ki67 (Ki67 ≤15: 24% vs 11%, p=0.06). Most tumors were T1-2 (94% vs 93%, p=0.82) and N0 (53% vs 62%, p=0.09), respectively. As expected, expression of most immune-related proteins was higher in the stromal vs the intratumoral compartment for both LAR and non-LAR TNBC. When focusing on the stromal compartment, expression of multiple immune related proteins was significantly lower in LAR compared to non-LAR TNBC, including the pan-leukocyte marker CD45 (log-2 fold change [log2FC]: 0.552, p=0.05), the macrophage marker CD14 (log2FC: 0.834, p=0.06), CD44 (lof2FC: 0.637, p=0.07), and the immune checkpoint proteins IDO1 (log2FC: 0.914, p=0.04), VISTA (log2FC: 0.471, p=0.07), ICOS (log2FC: 0.444, p=0.08), and STING (log2FC: 0.544, p=0.09). Proteins with expression levels too low for comparisons included PD-L1, LAG3, FOXP3 and BCL-2. When focusing on the intratumoral compartment, expression of most immune-related proteins was very low in both LAR and non-LAR TNBC. Like in the stromal compartment, CD45 expression was lower in LAR TNBC (log2FC: 0.78, p=0.02). Expression of the immune checkpoint B7-H3 was lower in LAR TNBC (log2FC: 0.737, p=0.02), while expression of the T cell marker CD127 was higher (log2FC: -0.528, p=0.34). With regards to relevant non-immune markers, expression of Ki67 was lower in LAR TNBC (log2FC: 0.5498, p=0.05), consistent with the clinical assay. Conclusion: In this ultra high-plex spatial analysis, we provide first insights into the differential expression at the protein level of several targetable immune checkpoint molecules in LAR vs non-LAR TNBC. The lower expression of several immune related proteins in LAR TNBC is consistent with the hypothesis that LAR TNBC exhibits a “cold” immune microenvironment compared to other TNBC subtypes, potentially rendering itself less susceptible to immunotherapy-based strategies. These data support the need to consider TNBC molecular subtypes in future evaluations of immune-based therapeutic approaches. Funding: This work was supported by NIH grant P50CA116201 to RLF, JMC, KRK, FJC, DZ, JNI, and MPG; BCRF grant 19-161 to EAT and NCATS grant CTSA KL2 TR002379 to RLF. The contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIH

Citation Format: Roberto A Leon-Ferre, Jodi M. Carter, David M. Zahrieh, David W. Hillman, Saranya Chumsri, Yaohua Ma, Jennifer M. Kachergus, Xue Wang, Judy C. Boughey, Minetta C. Liu, James N. Ingle, Krishna R. Kalari, Jose C. Villasboas Bisneto, Fergus J. Couch, E. Aubrey Thompson, Matthew P. Goetz. Digital spatial profiling of immune-related proteins in luminal androgen receptor (LAR) vs non-LAR triple-negative breast cancer (TNBC) [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 P1-04-01.

Abstract P5-13-22: Serum thymidine kinase 1 activity (TKa) levels and progression-free survival (PFS) in patients (pts) with hormone receptor positive (HR+) HER2-negative metastatic breast cancer (MBC) on palbociclib (Pb) and endocrine therapy (ET)

Background: Cyclin dependent 4/6 kinase inhibitors (CDK4/6i) and endocrine therapy (ET) have improved progression-free survival (PFS) and overall survival in HR+ MBC, but progression of disease ultimately occurs. Apart from HR+ status, there are no clinically available biomarkers that enable oncologists to determine prognosis and predict response to CDK4/6i. An emerging biomarker is serum thymidine kinase 1 (TK1), a secreted marker of proliferation that is prognostic in pts with HR+ HER2- MBC. High levels of TKa are associated with inferior PFS, whereas pts with low TKa levels pretreatment, or TKa levels that decrease on ET and a CDK4/6i, have superior PFS. Notably, TKa levels rebound ≥ 5 days off Pb, with resumption of cell cycling. PROMISE (NCT0281902) is a prospective study that enrolled women with HR+ MBC starting Pb + letrozole (L) in the 1st line [FL] or Pb + fulvestrant in the 2nd line [SL] setting. The trial includes a comprehensive “omic” assessment of blood, tumor, urine and the fecal microbiome to identify novel genomic variants and pathways associated with an early decline in TKa (measured after 2 months or end of cycle [C]2) and PFS. Here, we report the association between i) pre-treatment TKa (pre-TKa) levels and PFS (i.e. from registration to the 1st disease event) and ii) TKa levels at the end of C2 (C2-TKa) and PFS-2 (i.e. from the start of C3 to the 1st disease event).Methods: TKa testing was performed using the DiviTum assay (Biovica). TKa+ disease was defined as ≥ 200 Du/L and TKa- disease as below limit of detection to 200 Du/L. Log-rank test and univariate Cox modeling were used to assess the association between pre-TKa levels and PFS and between end of C2-TKa levels and PFS-2. The database was locked on June 28, 2021. Results: Of 68 pts enrolled, 4 were ineligible and pre-TKa data was unavailable for 4. Of the remaining 60 pts (45 FL, 15 SL), the percentage of pts with pre-TKa+ disease was 33.3% in FL (15/45, 95% CI: 20.0-49.0%), and 46.7% (7/15, 95% CI: 21.4-71.9%) in the SL. The median follow-up time for pts on study was 24 months (range: 2-42 months). There were 22 disease events (13 in FL, 9 in SL). In the FL setting, PFS was significantly shorter for preTKa+ pts compared to preTKa- pts (HR: 4.15, 95% CI:1.35-12.74; p=0.007), but not for SL pts (HR: 1.11, 95% CI: 0.30-4.18, p=0.875). End of C2 TKa data was obtained for pts while on Pb (n=5), or after stopping Pb as follows: 1-4 days (n=9), 5-8 days (n=28) and 9-36 days (n=11). PFS-2 was not associated with C2-TKa in the FL (p=0.834) or SL (p=0.454) settings. An analysis of TKa levels by metastatic site will be presented at the meeting.Conclusions: A secreted biomarker of proliferation (TK1) obtained prior to initiating CDK4/6i and ET for the treatment of HR+ MBC is associated with PFS in pts receiving 1st line Pb + L, but not in those receiving 2nd line Pb + fulvestrant. While the end of C2 TKa levels were not associated with PFS, the interpretability of these data are limited, given treatment delays (0-36 days) prior to starting C3 that may result in TKa rebound. Future studies evaluating the predictive nature of TKa and Pb response should focus on earlier timepoints while on drug.

Citation Format: Ciara C O'Sullivan, Jun He, Jason Sinnwell, Vera J Suman, Krishna R Kalari, Peter T Vedell, Ann M Moyer, Xiaojia Tang, Kevin J Thompson, Abe Eyman Casey, Alvaro Moreno-Aspitia, Donald W Northfelt, Minetta C Liu, Tufia C Haddad, Saranya Chumsri, Brendan McMenomy, Prema Peethambaram, Kathryn J Ruddy, Karthik V Giridhar, Roberto A Leon-Ferre, Mattias Bergqvist, Adrian Nordmark, Richard M Weinshilboum, Liewei Wang, Matthew P Goetz. Serum thymidine kinase 1 activity (TKa) levels and progression-free survival (PFS) in patients (pts) with hormone receptor positive (HR+) HER2-negative metastatic breast cancer (MBC) on palbociclib (Pb) and endocrine therapy (ET) [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-13-22.

Abstract P4-01-03: Multiomics data reveal novel biomarkers for CDK4/6 resistance

Background: Cyclin-dependent 4/6 kinase inhibitors (CDK4/6i) and endocrine therapy (ET) have improved progression-free survival (PFS) and overall survival in hormone-receptor-positive (HR+) metastatic breast cancer (MBC), but endocrine resistance is a major challenge. PROMISE [NCT0281902; n=63] is a multicenter study that enrolled women with HR+ HER2- MBC commencing palbociclib (Pb) with letrozole (1st line [1L]) or fulvestrant (2nd line [2L]), and was designed to perform a comprehensive “omic” assessment of prospectively collected biospecimens (pre-treatment (M1), at 2 months (M2), and at disease progression). The goal is to identify novel genomic variants and pathways associated with resistance to CDK4/6i and ET and PFS outcomes. Here we report the association between the proteomic, metabolomics, and lipidomics data generated from pre-Pb and 2-month serum samples and PFS. Methods: Untargeted mass spectrometry data was generated from Metabolon, assaying 1308 metabolites and 831 lipids. Additionally, 1436 proteins were assayed on the Olink platform. Cox proportional hazard models were used to evaluate the univariate hazard ratio (HR) for all features with respect to PFS. The analyses were performed on samples from 45 patients (N=33 1Lwith 9 progression events and 12 2L with 8 progression events), obtained from M1 and M2 timepoints on Pb + ET. Enrichment analysis p-values are calculated using Fisher’s exact test. Results: Proteomics: In the M1 timepoint, 93 and 43 proteins were associated with PFS in the 1L and 2L settings, respectively; inflammation genes were enriched among the 1L setting (p= 0.034); 33 proteins presented HRs ranging between 0.026 and 0.56. The FABP9 protein (HR of 1.98, 95% CI 1.02-3.83) was associated with worse PFS. Conversely, inflammation genes were not observed to be enriched in 2L. In the M2 timepoint, we observed 60 and 21 proteins significantly associated with PFS, but no biological function was enriched in 1L and 2L. Metabolites: In the M1 timepoint, metabolism of the sulfur-containing amino acids (methionine, cysteine, SAM and taurine) were enriched in the 1L setting (p= 0.035, HR range 0.15-0.33); and the branched-chain amino acids (leucine, isoleucine, and valine) were significantly associated with PFS in the 2L setting (p= 0.028, HR range 0.013-0.33). At the M2 timepoint, the amino acids were no longer enriched, but fatty acid metabolism was significantly enriched for both 1L and 2L (p= 0.048 and 0.067, respectively). Pathways involving lipids, amino acids, and xenobiotics were enriched in metabolites related to PFS (p &lt;0.05) for both treatment lines at M1 and M2. Lipidomics: In the M1 timepoint, 10 and 19 lipids were associated with PFS for 1L and 2L, respectively. The most notable lipid associated with worse PFS in the 1L was an 18 carbon phosphatidylinositol, PI(18:1/18:2), (HR 7.34 (CI 1.27-42.50); 8 triglycerides were associated with improved PFS (HR range 0.39 and 0.55). In 2L, the 19 lipids associated with PFS included 12 phosphatidylcholines (enrichment p = 5.6X10-8). In the M2 timepoint, 15 and 8 lipids were significantly associated with PFS for 1L and 2L. An enrichment of phosphatidylinositols was observed in 1L (p= 1.2X10-5); none were observed in the 2L.Future Directions: Networks are being constructed using the proximity scores of the proteins, lipids, and metabolites associated with PFS in M1 and M2 for 1L and 2L. Network similarities and analyses will be conducted.Conclusion: Distinct multi-omic changes identified in serum samples obtained from PROMISE participants M1 and M2 on Pb correlate with disease progression in both 1L and 2L settings. Additionally, validation studies will determine the significance of these findings.

Citation Format: Krishna R. Kalari, Kevin J. Thompson, Jason Sinnwell, Xiaojia Tang, Vera J. Suman, Jun He, Seul Kee Byeon, Akhilesh Pandey, Abe Eyman Casey, Peter T. Vedell, Ann M. Moyer, Alvaro Moreno-Aspitia, Donald W. Northfelt, Minetta C. Liu, Tufia C. Haddad, Saranya Chumsri, Prema Peethambaram, Kathryn J. Ruddy, Karthik V. Giridhar, Roberto A. Leon-Ferre, Richard M. Weinshilboum, Liewei Wang, Ciara C. O’ Sullivan, Matthew P. Goetz. Multiomics data reveal novel biomarkers for CDK4/6 resistance [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 P4-01-03.

Abstract P4-01-05: Multi-omics data shows downregulation of mismatch repair, purin and tublin pathways in AR-negative triple-negative chemotherapy-resistant tumors

Introduction. The standard treatment for non-metastatic triple-negative breast cancer (TNBC) is neoadjuvant chemotherapy (NAC) and nearly 50% exhibit pathological complete response (pCR). However, patients with residual disease after NAC are at increased risk for recurrence and death. Prior studies examining the transcriptome of TNBC pre/post-NAC have examined a limited number of genes (&lt;500) in heterogeneous subgroups of TNBC (e.g. LAR and non-LAR). We explored the transcriptome of androgen-receptor (AR) negative (non-LAR) TNBC subtype both pre/post NAC to identify pathways associated with NAC response. Methods. Tumors obtained pre/post NAC from TNBC patients enrolled in the Breast Cancer Genome Guided therapy study (BEAUTY) underwent RNA sequencing and reverse-phase protein array (RPPA). EdgeR was applied for differentially expressed (DE) analysis and regression methods for RPPA. Digital deconvolution method (CIBERSORTx) and TNBC single-cell data were used to obtain cell types. Pathway analysis was carried out using 2972 gene sets and gene set variation analysis (GSVA). Functional enrichment analysis was conducted with significant genes. Results. Of the 44 TNBC patients, 32 patients were excluded from the analysis cohort due to: LAR tumor (6 pts.), non-LAR tumor with pCR (23 pts.), and cell type issues with RNA-seq data (3 pt.). Paired RNA-Seq data were available for 12 TNBC patients (4 with progression &lt;2 years [EP]) and 8 who were progression-free &gt; 4 years [NP]) and paired RPPA data were available for 9 of these 12 patients. Differentially expressed genes, proteins and cell types between EP and NP in post-NAC. We identified 489 genes differentially expressed (DE) between EP and NP (logFC=|2|, FDR &lt; 0.05). Analysis of cytobands from these 489 genes showed an enrichment of genes on chromosome 6p22.1-2 and 17q25.3 regions (enrichment ratio &gt;5; p-value &lt;10E-4). Critical genes identified in the AR- network (p-value &lt; 10E-3) were IL1RN, SLAMF9, KRT81, BHLHE22, B3GALT5, PCP4, TREM1, AQP9, NRTN, and COL2A1.In addition, preliminary results from RPPA data of post-NAC tumors showed astrocytic phosphoprotein (PEA-15), involved in apoptosis, proliferation, glucose metabolism, as well as cell proliferation and Y box binding (YB1) proteins (involved in metastases), were more DE in EP than NP (p &lt; 0.05). CIBERSORTx was applied to estimate the proportions of different cell types in post-NAC tumors. Cancer-associated fibroblasts iCAFs were low and myCAFs are high in EP vs NP. It is known that the cross-talk between CAFs and tumor cells may induce tumor resistance to chemotherapy. Differentially expressed pathways in post and pre-NAC EP tumors. Using genome-wide expression data from the paired 12 tumors and the GSVA method, we obtained individual pathway scores for 2972 pathways. One hundred ninety pathways were downregulated and 61 pathways were upregulated (p-value &lt;= 0.05) in the post-NAC residual disease of EP relative to NP. We further examined these 190 pathways in the paired EPs and found 71% of those pathways were upregulated in the pre-NAC. These 190 downregulated pathways were enriched with FOXO, TGF-beta, PI3k, FGFR1, insulin and others. The 61 upregulated pathways in post-NAC EP tumors were enriched with mismatch repair, purine, tubulin, telomere, polymerase and gap-junction related pathways; 77% of those 61 pathways were downregulated in pre-NAC. Conclusions. Using a comprehensive “omics” approach, we have identified novel cancer and drug response pathways associated with recurrence in AR-TNBC disease. Further work to evaluate these as markers of outcome and potential drug targets is warranted.

Citation Format: Krishna R Kalari, Vera J Suman, Xiaojia Tang, Jason P Sinnwell, Kevin J Thompson, Peter T Vedell, Jodi M Carter, Sarah A McLaughlin, Alvaro Moreno Aspitia, Donald W Northfelt, Richard J Gray, Richard Weinshilboum, Liewei Wang, Judy C Boughey, Matthew Goetz. Multi-omics data shows downregulation of mismatch repair, purin and tublin pathways in AR-negative triple-negative chemotherapy-resistant tumors [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 P4-01-05.

Anastrozole Regulates Fatty Acid Synthase in Breast Cancer

Our previous matched case-control study of postmenopausal women with resected early-stage breast cancer revealed that only anastrozole, but not exemestane or letrozole, showed a significant association between the 6-month estrogen concentrations and risk of breast cancer. Anastrozole, but not exemestane or letrozole, is a ligand for estrogen receptor α. The mechanisms of endocrine resistance are heterogenous and with the new mechanism of anastrozole, we have found that treatment of anastrozole maintains fatty acid synthase (FASN) protein level by limiting the ubiquitin-mediated FASN degradation, leading to increased breast cancer cell growth. Mechanistically, anastrozole decreases the guided entry of tail-anchored proteins factor 4 (GET4) expression, resulting in decreased BCL2-associated athanogene cochaperone 6 (BAG6) complex activity, which in turn, prevents RNF126-mediated degradation of FASN. Increased FASN protein level can induce a negative feedback loop mediated by the MAPK pathway. High levels of FASN are associated with poor outcome only in patients with anastrozole-treated breast cancer, but not in patients treated with exemestane or letrozole. Repressing FASN causes regression of breast cancer cell growth. The anastrozole-FASN signaling pathway is eminently targetable in endocrine-resistant breast cancer.

Age-Dependent Changes in the Plasma and Brain Pharmacokinetics of Amyloid-β Peptides and Insulin

BACKGROUND

Age is the most common risk factor for Alzheimer's disease (AD), a neurodegenerative disorder characterized by the hallmarks of toxic amyloid-β (Aβ) plaques and hyperphosphorylated tau tangles. Moreover, sub-physiological brain insulin levels have emerged as a pathological manifestation of AD.

OBJECTIVE

Identify age-related changes in the plasma disposition and blood-brain barrier (BBB) trafficking of Aβ peptides and insulin in mice.

METHODS

Upon systemic injection of 125I-Aβ40, 125I-Aβ42, or 125I-insulin, the plasma pharmacokinetics and brain influx were assessed in wild-type (WT) or AD transgenic (APP/PS1) mice at various ages. Additionally, publicly available single-cell RNA-Seq data [GSE129788] was employed to investigate pathways regulating BBB transport in WT mice at different ages.

RESULTS

The brain influx of 125I-Aβ40, estimated as the permeability-surface area product, decreased with age, accompanied by an increase in plasma AUC. In contrast, the brain influx of 125I-Aβ42 increased with age, accompanied by a decrease in plasma AUC. The age-dependent changes observed in WT mice were accelerated in APP/PS1 mice. As seen with 125I-Aβ40, the brain influx of 125I-insulin decreased with age in WT mice, accompanied by an increase in plasma AUC. This finding was further supported by dynamic single-photon emission computed tomography (SPECT/CT) imaging studies. RAGE and PI3K/AKT signaling pathways at the BBB, which are implicated in Aβ and insulin transcytosis, respectively, were upregulated with age in WT mice, indicating BBB insulin resistance.

CONCLUSION

Aging differentially affects the plasma pharmacokinetics and brain influx of Aβ isoforms and insulin in a manner that could potentially augment AD risk.

Characteristics and Spatially Defined Immune (micro)landscapes of Early-stage PD-L1-positive Triple-negative Breast Cancer

PURPOSE

Programmed death ligand 1 [PD-(L)1]-targeted therapies have shown modest survival benefit in triple-negative breast cancer (TNBC). PD-L1+ microenvironments in TNBC are not well characterized and may inform combinatorial immune therapies. Herein, we characterized clinicopathologic features, RNA-based immune signatures, and spatially defined protein-based tumor-immune microenvironments (TIME) in early-stage PD-L1+ and PD-L1- TNBC.

EXPERIMENTAL DESIGN

From a large cohort of chemotherapy-naïve TNBC, clinicopathologic features, deconvoluted RNA immune signatures, and intraepithelial and stromal TIME (Nanostring GeoMX) were identified in subsets of PD-L1+ and PD-L1- TNBC, as defined by FDA-approved PD-L1 companion assays.

RESULTS

228 of 499 (46%) TNBC were PD-L1+ (SP142: ≥1% immune cells-positive). Using PD-L1 22C3, 46% had combined positive score (CPS) ≥ 1 and 16% had CPS ≥10. PD-L1+ TNBC were higher grade with higher tumor-infiltrating lymphocytes (TIL; P < 0.05). PD-L1 was not associated with improved survival following adjustment for TILs and other variables. RNA profiles of PD-L1+ TNBC had increased dendritic cell, macrophage, and T/B cell subset features; and decreased myeloid-derived suppressor cells. PD-L1+ stromal and intraepithelial TIMEs were highly enriched in IDO-1, HLA-DR, CD40, and CD163 compared with PD-L1-TIME, with spatially specific alterations in CTLA-4, Stimulator of Interferon Genes (STING), and fibronectin. Macrophage- and antigen presentation-related proteins correlated most strongly with PD-L1 protein.

CONCLUSIONS

In this early-stage TNBC cohort, nearly 50% were PD-L1+ (SP142 companion assay) while 16% were PD-L1+ with the 22C3 companion assay. PD-L1+ TNBC had specific myeloid-derived and lymphoid features. Spatially defined PD-L1+ TIME were enriched in several clinically actionable immune proteins. These data may inform future studies on combinatorial immunotherapies for patients with PD-L1+ TNBC.See related commentary by Symmans, p. 5446.

Interaction Between SNP Genotype and Efficacy of Anastrozole and Exemestane in Early-Stage Breast Cancer

Aromatase inhibitors (AIs) are the treatment of choice for hormone receptor-positive early breast cancer in postmenopausal women. None of the third-generation AIs are superior to the others in terms of efficacy. We attempted to identify genetic factors that could differentiate between the effectiveness of adjuvant anastrozole and exemestane by examining single-nucleotide polymorphism (SNP)-treatment interaction in 4,465 patients. A group of SNPs were found to be differentially associated between anastrozole and exemestane regarding outcomes. However, they showed no association with outcome in the combined analysis. We followed up common SNPs near LY75 and GPR160 that could differentiate anastrozole from exemestane efficacy. LY75 and GPR160 participate in epithelial-to-mesenchymal transition and growth pathways, in both cases with SNP-dependent variation in regulation. Collectively, these studies identified SNPs that differentiate the efficacy of anastrozole and exemestane and they suggest additional genetic biomarkers for possible use in selecting an AI for a given patient.

Multi-Omics Analyses Show Disease, Diet, and Transcriptome Interactions With the Virome

BACKGROUND & AIMS

The gut virome includes eukaryotic viruses and bacteriophages that can shape the gut bacterial community and elicit host responses. The virome can be implicated in diseases, such as irritable bowel syndrome (IBS), where gut bacteria play an important role in pathogenesis. We provide a comprehensive and longitudinal characterization of the virome, including DNA and RNA viruses and paired multi-omics data in a cohort of healthy subjects and patients with IBS.

METHODS

We selected 2 consecutive stool samples per subject from a longitudinal study cohort and performed metagenomic sequencing on DNA and RNA viruses after enriching for viral-like particles. Viral sequence abundance was evaluated over time, as well as in the context of diet, bacterial composition and function, metabolite levels, colonic gene expression, host genetics, and IBS subsets.

RESULTS

We found that the gut virome was temporally stable and correlated with the colonic transcriptome. We identified IBS-subset-specific changes in phage populations; Microviridae, Myoviridae, and Podoviridae species were elevated in diarrhea-predominant IBS, and other Microviridae and Myoviridae species were elevated in constipation-predominant IBS compared to healthy controls. We identified correlations between subsets of the virome and bacterial composition (unclassifiable "dark matter" and phages) and diet (eukaryotic viruses).

CONCLUSIONS

We found that the gut virome is stable over time but varies among subsets of patients with IBS. It can be affected by diet and potentially influences host function via interactions with gut bacteria and/or altering host gene expression.

FOXA1 overexpression suppresses interferon signaling and immune response in cancer

Androgen receptor-positive prostate cancer (PCa) and estrogen receptor-positive luminal breast cancer (BCa) are generally less responsive to immunotherapy compared with certain tumor types such as melanoma. However, the underlying mechanisms are not fully elucidated. In this study, we found that FOXA1 overexpression inversely correlated with interferon (IFN) signature and antigen presentation gene expression in PCa and BCa patients. FOXA1 bound the STAT2 DNA-binding domain and suppressed STAT2 DNA-binding activity, IFN signaling gene expression, and cancer immune response independently of the transactivation activity of FOXA1 and its mutations detected in PCa and BCa. Increased FOXA1 expression promoted cancer immuno- and chemotherapy resistance in mice and PCa and BCa patients. These findings were also validated in bladder cancer expressing high levels of FOXA1. FOXA1 overexpression could be a prognostic factor to predict therapy resistance and a viable target to sensitize luminal PCa, BCa, and bladder cancer to immuno- and chemotherapy.

Establishment and characterization of immortalized human breast cancer cell lines from breast cancer patient-derived xenografts (PDX)

The application of patient-derived xenografts (PDX) in drug screening and testing is a costly and time-consuming endeavor. While cell lines permit extensive mechanistic studies, many human breast cancer cell lines lack patient characteristics and clinical treatment information. Establishing cell lines that retain patient's genetic and drug response information would enable greater drug screening and mechanistic studies. Therefore, we utilized breast cancer PDX from the Mayo Breast Cancer Genome Guided Therapy Study (BEAUTY) to establish two immortalized, genomically unique breast cancer cell lines. Through extensive genetic and therapeutic testing, the cell lines were found to retain the same clinical subtype, major somatic alterations, and drug response phenotypes as their corresponding PDX and patient tumor. Our findings demonstrate PDX can be utilized to develop immortalized breast cancer cell lines and provide a valuable tool for understanding the molecular mechanism of drug resistance and exploring novel treatment strategies.

Patient-Derived Xenograft Engraftment and Breast Cancer Outcomes in a Prospective Neoadjuvant Study (BEAUTY)

PURPOSE

Patient-derived xenografts (PDX) are a research tool for studying cancer biology and drug response phenotypes. While engraftment rates are higher for tumors with more aggressive characteristics, it is uncertain whether engraftment is prognostic for cancer recurrence.

PATIENTS AND METHODS

In a prospective study of patients with breast cancer treated with neoadjuvant chemotherapy (NAC) with taxane ± trastuzumab followed by anthracycline-based chemotherapy, we report the association between breast cancer events and PDX engraftment using tumors derived from treatment naïve (pre-NAC biopsies from 113 patients) and treatment resistant (post-NAC at surgery from 34 patients). Gray test was used to assess whether the cumulative incidence of a breast cancer event differs with respect to either pre-NAC PDX engraftment or post-NAC PDX engraftment.

RESULTS

With a median follow-up of 5.7 years, the cumulative incidence of breast cancer relapse did not differ significantly according to pre-NAC PDX engraftment (5-year rate: 13.6% vs. 13.4%; P = 0.89). However, the incidence of a breast event was greater for patients with post-NAC PDX engraftment (5-year rate: 50.0% vs. 19.6%), but this did not achieve significance (P = 0.11).

CONCLUSIONS

In treatment-naïve breast cancer receiving standard NAC, PDX engraftment was not prognostic for breast cancer recurrence. Further study is needed to establish whether PDX engraftment in the treatment-resistant setting is prognostic for cancer recurrence.

Long-term outcomes of patients with node-negative (N0), triple-negative breast cancer (TNBC) who did not receive adjuvant chemotherapy according to stromal TILs (sTILs)

548

Background: sTILs are a well-established prognostic and predictive biomarker in patients with operable TNBC receiving pre or postoperative systemic therapy. We1 and others2,3 have also shown that sTILs are prognostic in patients who did not receive adjuvant chemotherapy. Here, we detail the outcomes of systemically untreated patients with N0 TNBC according to sTIL score. We focused on the N0 subset as a group of patients who may be candidates for future prospective therapy de-escalation trials. Methods: From a clinically annotated cohort of 605 patients with centrally confirmed TNBC (ER/PR < 1% and HER2 negative) with long-term outcomes data, we identified 182 patients treated with locoregional therapy only (breast surgery +/- radiation therapy and no chemotherapy). The clinicopathological characteristics of this cohort have previously been published1. In this analysis, we report the 5- and 10-year invasive disease-free survival (iDFS) and overall survival (OS) rates of patients with N0 TNBC according to sTIL levels. IDFS and OS were defined as per the STEEP classification and estimated using the Kaplan–Meier method. Comparisons of the survival distributions between groups were assessed by the log-rank test. sTILs were assessed as a continuous parameter according to the International TIL Working Group guidelines. For comparisons of outcomes between groups, tumors were classified as lymphocyte-predominant TNBC (defined as containing ≥50% sTILs) vs non-lymphocyte-predominant ( < 50% sTILs). Results: Of 182 systemically untreated patients, 149 (82%) were N0 and most (78%) were post-menopausal. T stage distribution was T1: 68%, T2: 28%, T3/4: 4%. Among N0 patients, 31 (21%) had lymphocyte-predominant TNBC, and in this group the 5-year iDFS and OS were 89% (95% CI 76-100) and 96% (95% CI 89-100), while the 10-year iDFS and OS were 89% (95% CI 76-100) and 87% (95% CI 73-100), respectively. In contrast, outcomes for patients with non-lymphocyte predominant TNBC were significantly worse. For this group, 5-year iDFS and OS were 62% (95% CI 53-73) and 78% (95% CI 71-86) while the 10-year iDFS and OS were 45% (95% CI 36-58) and 66% (95% CI 68-76), respectively ( log-rank p = 0.02 for iDFS and log-rank p = 0.03 for OS). Conclusions: sTIL quantification identifies a subset of patients with early-stage N0 TNBC with an exceedingly good prognosis, even in the absence of adjuvant chemotherapy. These data provide support for the evaluation of sTILs as part of prospective investigation of systemic therapy de-escalation strategies in N0 TNBC. References:1Leon-Ferre et al, Breast Cancer Res Treat (2018) 167:89-99 2Park et al, Ann Oncol (2019) 12:1941-1949 3De Jong et al, ESMO 2020

SLCO1B1: Application and Limitations of Deep Mutational Scanning for Genomic Missense Variant Function

SLCO1B1 (solute carrier organic anion transporter family member 1B1) is an important transmembrane hepatic uptake transporter. Genetic variants in the SLCO1B1 gene have been associated with altered protein folding, resulting in protein degradation and decreased transporter activity. Next-generation sequencing (NGS) of pharmacogenes is being applied increasingly to associate variation in drug response with genetic sequence variants. However, it is difficult to link variants of unknown significance with functional phenotypes using "one-at-a-time" functional systems. Deep mutational scanning (DMS) using a "landing pad cell-based system" is a high-throughput technique designed to analyze hundreds of gene open reading frame (ORF) missense variants in a parallel and scalable fashion. We have applied DMS to analyze 137 missense variants in the SLCO1B1 ORF obtained from the Exome Aggregation Consortium project. ORFs containing these variants were fused to green fluorescent protein and were integrated into "landing pad" cells. Florescence-activated cell sorting was performed to separate the cells into four groups based on fluorescence readout indicating protein expression at the single cell level. NGS was then performed and SLCO1B1 variant frequencies were used to determine protein abundance. We found that six variants not previously characterized functionally displayed less than 25% and another 12 displayed approximately 50% of wild-type protein expression. These results were then functionally validated by transporter studies. Severely damaging variants identified by DMS may have clinical relevance for SLCO1B1-dependent drug transport, but we need to exercise caution since the relatively small number of severely damaging variants identified raise questions with regard to the application of DMS to intrinsic membrane proteins such as organic anion transporter protein 1B1. SIGNIFICANCE STATEMENT: The functional implications of a large numbers of open reading frame (ORF) "variants of unknown significance" (VUS) in transporter genes have not been characterized. This study applied deep mutational scanning to determine the functional effects of VUS that have been observed in the ORF of SLCO1B1(s olute carrier organic anion transporter family member 1B1). Several severely damaging variants were identified, studied, and validated. These observations have implications for both the application of deep mutational scanning to intrinsic membrane proteins and for the clinical effect of drugs and endogenous compounds transported by SLCO1B1.

COVID-19 Transmission, Current Treatment, and Future Therapeutic Strategies

At the stroke of the New Year 2020, COVID-19, a zoonotic disease that would turn into a global pandemic, was identified in the Chinese city of Wuhan. Although unique in its transmission and virulence, COVID-19 is similar to zoonotic diseases, including other SARS variants (e.g., SARS-CoV) and MERS, in exhibiting severe flu-like symptoms and acute respiratory distress. Even at the molecular level, many parallels have been identified between SARS and COVID-19 so much so that the COVID-19 virus has been named SARS-CoV-2. These similarities have provided several opportunities to treat COVID-19 patients using clinical approaches that were proven to be effective against SARS. Importantly, the identification of similarities in how SARS-CoV and SARS-CoV-2 access the host, replicate, and trigger life-threatening pathological conditions have revealed opportunities to repurpose drugs that were proven to be effective against SARS. In this article, we first provided an overview of COVID-19 etiology vis-à-vis other zoonotic diseases, particularly SARS and MERS. Then, we summarized the characteristics of droplets/aerosols emitted by COVID-19 patients and how they aid in the transmission of the virus among people. Moreover, we discussed the molecular mechanisms that enable SARS-CoV-2 to access the host and become more contagious than other betacoronaviruses such as SARS-CoV. Further, we outlined various approaches that are currently being employed to diagnose and symptomatically treat COVID-19 in the clinic. Finally, we reviewed various approaches and technologies employed to develop vaccines against COVID-19 and summarized the attempts to repurpose various classes of drugs and novel therapeutic approaches.

Abstract PD8-04: Antitumor activity of Z-endoxifen (ENDX) is mediated via PKCβ1-dependent ERα loss and cell cycle arrest in ERα-positive breast cancer

Background: The tamoxifen (TAM) metabolite, ENDX, demonstrated promising antitumor activity in endocrine resistant breast cancer (BC) in both phase I and phase II settings. Furthermore, ENDX resulted in superior in vivo antitumor activity compared to TAM and letrozole in aromatase-expressing aromatase inhibitor-sensitive and resistant MCF7AC1 models. Recently, we identified protein kinase C beta 1 (PKCβ1), which regulates cell proliferation and tumorigenic transformation, as a novel target of ENDX. ENDX-bound PKCβ1 at concentrations achieved in phase I/II ENDX studies (100-300 nM). In contrast, TAM binding to PKCβ1 occurred at concentrations 7-10 folds higher (2 μM) than achievable with TAM 20 mg/day dosing. However, the clinical relevance of targeting PKCβ1 kinase activity is unclear, since drugs that target PKCβ1 (enzastaurin) have been ineffective in BC and other solid tumors. Therefore, we sought to understand how ENDX altered PKCβ1 and to further compare and contrast ENDXs effects to that of PKCβ1 kinase inhibition in ERα+ BC. Methods: The effects of PKCβ1 silencing and ENDX treatment on gene expression was analyzed by RNAseq in MCF7AC1 cells. The impact of PKCβ1-silencing on cell cycle was evaluated by flow cytometry. Protein expression of cell cycle regulators in PKCβ1 and ENDX-treated MCF7AC1 and T47D cells were compared to TAM and enzastaurin in vitro and to letrozole, TAM or control in vivo. The effects of PKCβ1 and drugs on growth were analyzed by cell proliferation assays. PRKCB gene amplification was assessed in primary tumors using TCGA data and in metastatic tumors using whole-exome sequencing data from patients enrolled in the PROMISE study (NCT 03281902). Results: RNAseq analysis revealed E2F targets and G2M checkpoints as the top hallmark genesets significantly downregulated in both PKCβ1-silenced and ENDX-treated MCF7AC1 cells. Flow cytometry demonstrated that PKCβ1 silencing increased G1 and reduced S phases of the cell cycle. Western blot analyses of PKCβ1-silenced MCF7AC1 and T47D cells displayed reduced protein levels of the cell cycle regulators Cyclin D1, Retinoblastoma (Rb), phospho-RbS807/811, CDK4, Chk1 and E2F1 that regulate G1/S transition. While short term ENDX (48 hours) treatment did not alter PKCβ1 levels, prolonged in vitro ENDX treatment profoundly reduced PKCβ1 protein levels and the aforementioned cell cycle regulators, faithfully replicating PKCβ1 silencing effects. In contrast, enzastaurin had no impact on proliferation or cell cycle proteins in either model. Consistent with this finding, ENDX, but not TAM or letrozole, reduced protein levels of ERα and cell cycle regulators in vivo. Overexpression of PKCβ1 induced TAM, but not ENDX, resistance and had little impact on responsiveness to enzastaurin. While PRKCB gene amplification was uncommon in newly diagnosed ERα+/HER2- BC (5%, TCGA), PRKCB was amplified in 40% of metastatic ERα+/HER2- BC (PROMISE study). Conclusion: We have confirmed the relevance of a new ENDX target, PKCβ1, in ERα+/HER2- BC. While targeting PKCβ1 kinase activity elicited no anticancer effects in ERα+ cells, PKCβ1 downregulation, either by siRNA or ENDX, resulted in profound ERα turnover, reduced protein levels of essential cell cycle mediators and profoundly inhibited cell proliferation. Furthermore, PKCβ1 protein expression is associated with TAM, but not ENDX, resistance, a finding whose clinical relevance is further magnified by identification of PRKCB amplification in metastatic ERα+ BC, confirming its potential importance in progression. Efforts are currently underway to elucidate the mechanistic basis for ENDX-induced PKCβ1 and ERα degradation and the contribution of these effects to the superior antitumor activity of ENDX in ERα+ BC.

Citation Format: Swaathi Jayaraman, Mary J Kuffel, Krishna R Kalari, Kevin J Thompson, Xiaojia Tang, Vera J Suman, Elizabeth S Bruinsma, Ciara C O'Sullivan, Liewei Wang, Richard Weinshilboum, James N Ingle, John R Hawse, Matthew P Goetz. Antitumor activity of Z-endoxifen (ENDX) is mediated via PKCβ1-dependent ERα loss and cell cycle arrest in ERα-positive breast cancer [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PD8-04.

Abstract PS6-02: Spatially defined immune-related proteins and outcome in triple negative breast cancer in the FinXX trial and Mayo Clinic cohort

Background: Growing data established the pivotal role of preexisting immune response in triple negative breast cancer (TNBC). Conventionally, preexisting immune response can be evaluated by quantifying tumor infiltrating lymphocytes mainly in the stroma or gene expression analysis from the whole tumor section. Due to technical challenges with these conventional methods, limited data regarding specific subtypes and spatial distribution of these immune infiltrates are currently available. Methods: NanoString IO360 gene expression analysis and Digital Spatial Profiling (DSP) were used. DSP was used to quantify 29 immune-related proteins in stromal and tumor-enriched segments from 44 TNBC samples from the FinXX trial (NCT00114816) and 335 samples from the Mayo Clinic (MC) cohort of centrally reviewed TNBC (Leon-Ferre BCRT 2018). In FinXX trial, 22 patients with recurrence and 22 patients without recurrence were included. In MC cohort, 217/335 patients received adjuvant chemotherapy while 118 patients had surgery only without adjuvant chemotherapy. Regions were segmented based on pancytokeratin staining. The general linear model was used for statistical analysis of differential expression with recurrence free survival (RFS) as a categorical variable (recur yes or no). Kaplan-Meier (KM) estimates and Cox regression models were also used for analysis. Results: In the FinXX trial, there were 12 out of 29 proteins in tumor epithelial segments (intraepithelial) which were significantly expressed at higher levels among patients who were free of recurrence. These proteins include Beta-2 microglobulin, CD11c, CD20, CD40, CD56, CD8, Granzyme B, HLA-DR, ICOS, PD-L1, PD-L2, and TGFB1. In contrast, merely 5 out of 29 proteins in stromal segments were significantly differentially expressed in these 2 groups of patients. Granzyme B, IDO1, PD-L1, and PD-L2 in stroma were significantly higher and SMA was significantly lower in patients without recurrence. Using Cox regression models, intraepithelial CD56, CD40, and HLA-DR were significantly associated with outcome. When comparing between highest and lowest intraepithelial protein expression by tertile, intraepithelial CD56 (HR 0.12, 95%CI 0.03-0.39, p &lt; 0.001), CD40 (HR 0.13, 95%CI 0.04-0.46, p = 0.002), and HLA-DR (HR 0.24, 95%CI 0.06-0.89, p = 0.032) were significantly associated with improved outcome. However, expression of these same proteins in stroma was not associated with outcome. Using KM estimates, intraepithelial CD56 (p &lt; 0.0001), CD40 (p = 0.0006), and HLA-DR (p = 0.013) were also significantly associated with improved outcome. Nonetheless, RNA expression of these proteins by IO360 from whole tumor sections were not significantly associated with outcome (CD56 p = 0.27, CD40 p = 0.21, HLA-DR p = 0.48). Similar findings with DSP were observed in MC TNBC cohort. Comparing between the highest and lowest quartiles, there were significantly fewer patients who developed recurrence with high protein expression of intraepithelial CD56 (p &lt; 0.001), CD40 (p = 0.002), and HLA-DR (p = 0.006). Conclusions: Using an in-depth analysis with spatially defined context, we identify that there were numerically more intraepithelial immune-related proteins associated with outcome compared to proteins in stroma. Specifically, intraepithelial CD56, CD40, and HLA-DR were significantly associated with improved outcome in both FinXX and MC TNBC cohorts. However, neither expression of these proteins in stroma nor RNA expression from whole tumor were associated with outcome. Our study highlights the impact of spatial biology and the importance of evaluating each potential biomarker in a spatially defined manner. Support: W81XWH-15-1-0292, BCRF 19-161, P50CA116201-9, P50CA015083

Citation Format: Saranya Chumsri, Jodi M. Carter, Yaohua Ma, Douglas Hinerfeld, Heather Ann Brauer, Sarah Warren, Torsten O. Nielsen, Karama Asleh, Heikki Joensuu, Edith A. Perez, Roberto A. Leon-Ferre, David W. Hillman, Judy C. Boughey, Minetta C. Liu, James N. Ingle, Krishna R. Kalari, Fergus J. Couch, Keith L. Knutson, Matthew P. Goetz, E. A. Thompson. Spatially defined immune-related proteins and outcome in triple negative breast cancer in the FinXX trial and Mayo Clinic cohort [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS6-02.

Abstract PD7-05: Neoadjuvant chemotherapy selectively alters spatially-defined immune landscapes in clinical luminal B HR+/HER2- breast cancers: Analysis of the breast cancer genome guided therapy study (BEAUTY)

Introduction: The immune microenvironment of high-risk HR+ breast cancer (BC) is poorly understood. BEAUTY is a prospective neoadjuvant chemotherapy (NAC) study of stage I-III BC patients treated with neoadjuvant weekly taxane followed by anthracycline-based chemotherapy. Among clinical luminal B BC (St. Gallen criteria) from BEAUTY, we used high-plex digital spatial profiling (DSP) to characterize the intra-epithelial tumor and stromal immune landscapes to 1) assess the impact of NAC on these landscapes and 2) identify immune biomarkers predictive of response to NAC. Methods: The tissue set included FFPE sections of resected tumors from 35 patients (median 51y; range: 21-71y) with clinically-defined luminal B BC (ER &gt; 10%/grade 2/Ki67 ≥ 15% or ER &gt; 10%/grade 3), and 16 paired pre-NAC biopsies. Nanostring GeoMX DSP was used to quantitate 58 immune and BC biomarker proteins in intra-epithelial, cytokeratin-positive tumor segments and adjacent stromal (cytokeratin-negative/SYTO13 (nuclear stain)-positive) segments. Data were normalized using the geometric mean of two negative controls (RbIgG and MsIgG1). A general linear model with negative binomial identified differentially-expressed (DE) proteins in pre/post-NAC tumors. Based on DE protein data and biologic function, a targeted protein subset (N=19) was evaluated in pre-NAC biopsies for associations between protein abundance and residual cancer burden (RCB) class (0-II vs. III) (Wilcoxon-rank sum test: p-value &lt; 0.025 considered significant) or RCB ‘breast only’ index (Spearman correlation). Results: Comparing tumor segments in pre-NAC specimens, intra-epithelial segments were predictably enriched in cytokeratin, ER, PR and Ki-67, but also CD127 and NY-ESO-1; whereas stromal segments were enriched in proteins associated with T cell subsets (e.g. CD3, CD4, CD8), macrophages (CD68 and CD163), antigen presentation/dendritic cells (CD11c, HLA-DR) and immune checkpoint proteins (PD-L1, PD-L2, B7-H3, TIM-3, VISTA) among others. While NAC did not alter the spatial distribution of proteins (intra-epithelial vs. stromal segments), it markedly attenuated the immune landscape, with decreased abundance of functionally-diverse immune proteins (e.g. CD45, CD3, CD4, CD127, granzyme B, CTLA4; STING, B7-H3, CD11c, and CD68, log2FC: 0.27-1.52 p &lt; 0.05), including low abundance proteins PD-1, PD-L1, PD-L2, CD20, and OX40L (log 2FC: 0.15-1.16, p &lt; 0.05). Both PR and Ki-67 decreased in post-NAC tumors whereas ER was not significantly altered (p &lt; 0.05). CD8, CD14, CD163, HLA-DR, IDO-1 and TIM-3 were not significantly altered by NAC. In the pre-NAC biopsies (which had subsequent residual tumor burdens of RCB class 0 or I (n=1 each), II (n=6), and III (n=8), and median RCB ‘breast only’ index of 3.39 (range: 0.00-38.02), 19 proteins were used for RCB analysis (CD3, CD4, CD8, CD14, CD34, CD44, CD45, CD68, CD127, CD163, CTLA4, granzyme B, STING, B7-H3, fibronectin, Ki-67, HLA-DR, SMA and TIM3). Among them, stromal CD127 was significantly higher in RCB class III than RCB class 0-II (p=0.021) and RCB ‘breast only’ index was positively correlated with intra-epithelial granzyme B (rho = 0.61; p=0.012), and negatively correlated with intra-epithelial Ki-67 (rho= -0.71; p=0.0022). Conclusion: In this series of clinical luminal B BC, NAC markedly attenuated the tumoral immune landscapes with a small set of “NAC-resistant” immune proteins. Among a targeted set, stromal CD127 was significantly higher in RCB III, and RCB breast-only index was positively correlated with intra-epithelial granzyme B. These data provide insight into the impact of NAC on HR+ BC, and identify potential immune biomarkers to predict response to neoadjuvant chemotherapy.

Citation Format: Jodi M Carter, Judy C Boughey, Jun He, Vera J Suman, Xue Wang, Jennifer M Kachergus, Krishna R Kalari, Liewei Wang, Richard Weinshilboum, Ann M Moyer, Sarah A McLaughlin, Alvaro Moreno-Aspitia, Donald W Northfelt, Richard J Gray, James N Ingle, E. Aubrey Thompson, Matthew P Goetz. Neoadjuvant chemotherapy selectively alters spatially-defined immune landscapes in clinical luminal B HR+/HER2- breast cancers: Analysis of the breast cancer genome guided therapy study (BEAUTY) [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PD7-05.

Abstract PS5-24: Novel genomic variants and pathways associated with baseline serum thymidine kinase 1 levels in HR-positive HER2-negative MBC patients commencing palbociclib and letrozole

Background: Cyclin dependent 4/6 kinase inhibitors (CDK4/6i) and endocrine therapy (ET) have improved progression-free survival (PFS) and overall survival in hormone-receptor (HR)-positive metastatic breast cancer (MBC), but progression of disease is inevitable. Serum thymidine kinase-1 (TK1) is a secreted marker of proliferation that is prognostic in patients (pts) with HR-positive, HER2-negative MBC and may be predictive of ET and CDK 4/6i response. PROMISE (NCT0281902) is a prospective study enrolling women with HR-positive MBC starting palbociclib (Pb) + letrozole (L) (1st line) or Pb + fulvestrant (2nd line). We undertook a comprehensive “omic” assessment of blood, tumor, urine and the fecal microbiome in order to identify novel genomic variants and pathways associated with an early decline in TK1 (measured after 2 months) and PFS. Additionally, patient derived xenografts/organoids were generated at baseline and progression to test new therapeutic approaches to overcome resistance to CDK4/6i and ET. We report the initial association between the baseline genomic landscape and baseline TK1 levels. Methods: FFPE tumor biopsies were obtained for DNA/RNA sequencing (TempusTM) and blood samples for TK1 (Divitum® assay, Biovica) were collected pretreatment (pre-Pb) and after 2 cycles of Pb + ET (post-Pb2). Both whole-exome (exome capture) sequencing (WES) and RNA-Seq used the Integrated DNA Technologies xGen Exome Research Panel v1.0 capture kit. TK1+ disease was defined as &gt; 200 Du/L and TK1- disease as below limit of detection up to 200 Du/L. We tested the association between genomic and transcriptomic characteristics with baseline TK1 data in pretreatment samples where both WES and RNA-seq and TK1 was available. The data were analyzed using bioinformatics pipelines for somatic and germline mutations and copy number alterations. The current analysis focuses on baseline 1st-line pre-Pb omics data in conjunction with baseline TK1 levels. The database was locked for analysis on 5/29/2020. Results: Thirty-three pts (median age: 59 yrs.) were evaluable, with paired samples for TK1 in 32. Six pts had TK1+ disease pre-Pb and post-Pb2. Twenty-two pts had TK1- disease pre-Pb and post-Pb2. Four pts had a decrease in TK1 after 2 cycles of treatment that altered the classification from TK1+ to TK1-. Both baseline RNA seq and serum TK1 (n=16) were available for 4 TK1+ and 12 TK1- pts. In this group, 476 genes were differentially regulated (398 upregulated; 78 downregulated). Pathway analysis demonstrated enrichment in complement and coagulation cascade pathway, PPAR signaling pathway, and metabolism-related pathways related to up-regulation of CYP and UGT gene families. Further testing for the association of WES data with baseline TK1+ (n=8) and TK1- (n=16) disease demonstrated somatic copy number variations on chromosomes 6, 11, 12 and 15. CDK4 copy number gains were observed in 3/8 TK1+ pts and 0/16 TK1- pts. We also observed that somatic mutations (LOH, copy number and/or SNV/INDELs) were more prevalent in the TK1+ compared to the TK1- pre-Pb group in several cancer-associated genes (FAS [p=0.06] PTEN, PIK3CB, NAB2, SOX9 and FAT1 [p=0.08], TP53, and MAP2K4 [p=0.22]). Conversely, we also noted that 6/7 pts with GATA3 mutations had TK1- disease (p=0.23). Conclusions: Using a comprehensive “omics” approach, our data suggest that a secreted biomarker of proliferation (TK1) obtained prior to initiating CDK4/6i and ET for the first line treatment of HR+ MBC is associated with established and novel genes and pathways associated with prognosis of pts receiving ET and CDK 4/6i. Analysis of on-treatment (after 2 cycles) tumor RNA seq and its association with change in TK1 as well as data from the 2nd-line cohort will be presented at the meeting.

Citation Format: Ciara C O Sullivan, Krishna R Kalari, Vera J Suman, Peter T Vedell, Ann Moyer, Abraham D Eyman Casey, Jason Sinnwell, Xiaojia Tang, Kevin Thompson, Alvaro Moreno-Aspitia, Donald W Northfelt, Minetta C Liu, Tufia C Haddad, Saranya Chumsri, Prema Peethambaram, Kathryn J Ruddy, Karthik V Giridhar, Roberto A Leon-Ferre, Adrian Nordmark, Mattias Bergqvist, Brendan P McMenomy, Richard M Weinshilboum, Liewei Wang, Matthew P Goetz. Novel genomic variants and pathways associated with baseline serum thymidine kinase 1 levels in HR-positive HER2-negative MBC patients commencing palbociclib and letrozole [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS5-24.

Abstract PS16-01: Intra-epithelial tumor immune landscapes are associated with clinical outcomes in early-stage triple-negative breast cancer

Introduction: Stromal tumor-infiltrating lymphocytes (sTILs) have established prognostic and predictive significance in triple-negative breast cancer (TNBC). However, the roles of other immune cells in TNBC are less well-established. We performed high-plex quantitative spatial profiling in a cohort of early-stage TNBC to 1) apply spatial context to tumoral immune landscapes and 2) identify immune proteins associated with clinical outcomes, independently of TILs and other established prognostic clinicopathologic variables, in patients (pts) treated with or without adjuvant chemotherapy (CTX). Methods: The Mayo TNBC cohort comprises pts with centrally-verified, CTX-naive tumors resected from 1985-2012. Using a cohort-based TMA, with Nanostring GeoMX DSP, we quantitated 58 proteins within spatially-distinct intra-epithelial, cytokeratin-positive tumor segments and adjacent cytokeratin-negative/nuclei-positive stromal segments. Differentially-expressed (DE) proteins were identified using a negative binomial generalized linear model (SNR&gt;2, p&lt; 0.05) and a target DE protein set was dichotomized (80th percentile). After adjusting for prognostic clinicopathologic variables, proteins associated with recurrence-free survival (RFS, defined as time from surgery to either local, regional, and distant recurrence, or death by any cause) were identified by performing variable selection using the Akaike Information Criterion (AIC) obtained from fitting all possible Cox proportional hazards regression models (performed separately for intra-epithelial/stromal segments, and in groups +/- adjuvant CTX. Results: From the TNBC TMA, DSP data (N=250 tumors) included 169 pts who received adjuvant CTX+ and 81 who did not (CTX-). Overall, 85/250 developed recurrent disease. In the CTX+ group, intra-epithelial tumor segments from pts without recurrent disease were enriched in 10 immune proteins, including CD8, markers involved in antigen presentation/dendritic cells (CD11c, CD40, HLA-DR) or NK cells (CD56) (FC: 1.4-2.1, p&lt;0.05); CD14 was increased in stroma (FC: 1.5, p&lt;0.05). In contrast, in the CTX- group, both the intra-epithelial tumor and stromal segments from pts without recurrences were enriched in immune proteins (N= 12 and 15 respectively; FC 1.6-5.5, p&lt; 0.05) most markedly CD40, IDO1 and HLA-DR (FC: 3.2-5.5, p&lt; 0.05). Overall, CD3, CD4, CD27, CD44, and ICOS among others were enriched only in the CTX- group; CD14 and CD56 were enriched only in the CTX+ group. Based on these spatial data, biologic function and DSP data from another set of TNBC (FinXX trial), CD11c, CD14, CD27, CD40, CD56, and IDO1 were selected for RFS analysis. After applying our model selection criterion and adjusting for pt age at surgery, tumor size, lymph node status, and sTILs, intra-epithelial CD56 was independently associated with improved RFS in the CTX+ group (HR: 0.31[0.12, 0.81]). In the CTX- group, intra-epithelial CD11c was independently associated with improved RFS (0.10 [0.01, 0.81]). Conclusion: In this early-stage TNBC cohort, spatially-distinct tumor immune landscapes were associated with RFS but differed according to receipt of CTX after surgical resection. In the patients who received CTX, the intra-epithelial compartment, rather than stromal compartment, was immune-enriched in pts without recurrences. Among a targeted protein set, intra-epithelial CD56 remained associated with improved outcomes, independent of sTILs and other clinicopathologic features. In the CTX- group, spatial landscapes were more balanced, and intra-epithelial CD11c was independently associated with improved outcomes. These data provide insight into the spatial context of intrinsic immune landscapes in TNBC, and identify candidate prognostic immune biomarkers which may inform therapeutic strategies.

Citation Format: Jodi M Carter, Saranya Chumsri, David W Hillman, David M Zahrieh, Yaohua Ma, Xue Wang, Jennifer M Kachergus, Judy C Boughey, Minetta C Liu, Krishna R Kalari, JC Villasboas, Roberto A Leon Ferre, Fergus J Couch, Matthew P Goetz, E. Aubrey Thompson. Intra-epithelial tumor immune landscapes are associated with clinical outcomes in early-stage triple-negative breast cancer [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS16-01.

Single-nucleotide polymorphism biomarkers of adjuvant anastrozole-induced estrogen suppression in early breast cancer

OBJECTIVES

Based on our previous findings that postmenopausal women with estrone (E1) and estradiol (E2) concentrations at or above 1.3 pg/ml and 0.5 pg/ml, respectively, after 6 months of adjuvant anastrozole therapy had a three-fold risk of recurrence, we aimed to identify a single-nucleotide polymorphism (SNP)-based model that would predict elevated E1 and E2 and then validate it in an independent dataset.

PATIENTS AND METHODS

The test set consisted of 322 women from the M3 study and the validation set consisted of 152 patients from MA.27. All patients were treated with adjuvant anastrozole, had on-anastrozole E1 and E2 concentrations and genome-wide genotyping.

RESULTS

SNPs were identified from the M3 genome-wide association study. The best model to predict the E1-E2 phenotype with high balanced accuracy was a support vector machine model using clinical factors plus 46 SNPs. We did not have an independent cohort that is similar to the M3 study with clinical, E1-E2 phenotypes and genotype data to test our model. Hence, we chose a nested matched case-control cohort (MA.27 study) for testing. Our E1-E2 model was not validated but we found the MA.27 validation cohort was both clinically and genomically different.

CONCLUSIONS

We identified a SNP-based model that had excellent performance characteristics for predicting the phenotype of elevated E1 and E2 in women treated with anastrozole. This model was not validated in an independent dataset but that dataset was clinically and genomically substantially different. The model will need validation in a prospective study.

CDC25B partners with PP2A to induce AMPK activation and tumor suppression in triple negative breast cancer

Cell division cycle 25 (CDC25) dual specificity phosphatases positively regulate the cell cycle by activating cyclin-dependent kinase/cyclin complexes. Here, we demonstrate that in addition to its role in cell cycle regulation, CDC25B functions as a regulator of protein phosphatase 2A (PP2A), a major cellular Ser/Thr phosphatase, through its direct interaction with PP2A catalytic subunit. Importantly, CDC25B alters the regulation of AMP-activated protein kinase signaling (AMPK) by PP2A, increasing AMPK activity by inhibiting PP2A to dephosphorylate AMPK. CDC25B depletion leads to metformin resistance by inhibiting metformin-induced AMPK activation. Furthermore, dual inhibition of CDC25B and PP2A further inhibits growth of 3D organoids isolated from patient derived xenograft model of breast cancer compared to CDC25B inhibition alone. Our study identifies CDC25B as a regulator of PP2A, and uncovers a mechanism of controlling the activity of a key energy metabolism marker, AMPK.

Longitudinal Multi-omics Reveals Subset-Specific Mechanisms Underlying Irritable Bowel Syndrome

The gut microbiome has been implicated in multiple human chronic gastrointestinal (GI) disorders. Determining its mechanistic role in disease has been difficult due to apparent disconnects between animal and human studies and lack of an integrated multi-omics view of disease-specific physiological changes. We integrated longitudinal multi-omics data from the gut microbiome, metabolome, host epigenome, and transcriptome in the context of irritable bowel syndrome (IBS) host physiology. We identified IBS subtype-specific and symptom-related variation in microbial composition and function. A subset of identified changes in microbial metabolites correspond to host physiological mechanisms that are relevant to IBS. By integrating multiple data layers, we identified purine metabolism as a novel host-microbial metabolic pathway in IBS with translational potential. Our study highlights the importance of longitudinal sampling and integrating complementary multi-omics data to identify functional mechanisms that can serve as therapeutic targets in a comprehensive treatment strategy for chronic GI diseases. VIDEO ABSTRACT.

Pharmacogenomics of aromatase inhibitors in postmenopausal breast cancer and additional mechanisms of anastrozole action

Aromatase inhibitors (AIs) reduce breast cancer recurrence and prolong survival, but up to 30% of patients exhibit recurrence. Using a genome-wide association study of patients entered on MA.27, a phase III randomized trial of anastrozole versus exemestane, we identified a single nucleotide polymorphism (SNP) in CUB And Sushi multiple domains 1 (CSMD1) associated with breast cancer–free interval, with the variant allele associated with fewer distant recurrences. Mechanistically, CSMD1 regulates CYP19 expression in an SNP- and drug-dependent fashion, and this regulation is different among 3 AIs: anastrozole, exemestane, and letrozole. Overexpression of CSMD1 sensitized AI-resistant cells to anastrozole but not to the other 2 AIs. The SNP in CSMD1 that was associated with increased CSMD1 and CYP19 expression levels increased anastrozole sensitivity, but not letrozole or exemestane sensitivity. Anastrozole degrades estrogen receptor α (ERα), especially in the presence of estradiol (E2). ER⁺ breast cancer organoids and AI- or fulvestrant-resistant breast cancer cells were more sensitive to anastrozole plus E2 than to AI alone. Our findings suggest that the CSMD1 SNP might help to predict AI response, and anastrozole plus E2 serves as a potential new therapeutic strategy for patients with AI- or fulvestrant-resistant breast cancers.

A germline variation within the CSMD1 gene predicts aromatase inhibitor response in breast cancer.

Abstract 1928: Study of copy number amplification and chemo-response in triple negative breast cancer

Triple negative breast cancer (TNBC) accounts for 15-20% of breast cancer cases in the United States. With currently no targeted therapies for TNBC, identification of novel biomarkers to improve chemo-response is vital. Recently, copy number amplifications (CNAs) have been implicated in treatment response in several diseases including breast cancer. Our clinical studies have associated a 5 Mb CNA in TNBC patients with chemo-response, which may elucidate novel biomarkers and treatment strategies for TNBC. Analysis of RNAseq and Whole Exome Sequencing data from 42 TNBC patient tumors enrolled in the BEAUTY Study, identified a 5 Mb CNA on chromosome 8q which associated with pathological complete response following taxane and anthracycline chemotherapy. By performing siRNA knockdown cytotoxicity studies in TNBC cell lines for genes within this amplified region, we identified RDH10 as a potential novel biomarker for increasing chemo-sensitivity. RDH10 is the rate limiting step in the synthesis of all-trans retinoic acid (ATRA). Previous studies have demonstrated ATRA to be a potent inhibitor of PIN1, a key regulator of several oncogenic signaling pathways, amplified in TNBC. We discovered overexpression of RDH10 to increase endogenous levels of ATRA, contributing to the degradation of PIN1 and taxane sensitivity. A rescue experiment with PIN1 was able to restore the chemo-resistance phenotype. Similarly, we found treating TNBC cells with ATRA phenocopied RDH10's effect on PIN1 and improved chemo-sensitivity. Furthermore, analyses of chemo-responsive PDX tumors containing the 5 Mb CNA, overexpressed RDH10 and exhibited lower levels of PIN1 when compared to non-responders. In conclusion, RDH10 amplification increases endogenous levels of ATRA and degrades PIN1 in TNBC. ATRA-mediated degradation of PIN1 sensitizes TNBC to chemotherapy, suggesting RDH10 as a novel biomarker for chemo-response and ATRA as a novel treatment to improve chemo-sensitivity.

Citation Format: Jordan Mark Grainger, Jia Yu, Krishna R. Kalari, Peter Vedell, Kevin Thompson, Matthew P. Goetz, Vera J. Suman, Judy C. Boughey, Liewei Wang. Study of copy number amplification and chemo-response in triple negative breast cancer [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1928.

Abstract LB-103: The interaction between SNP genotype and aromatase inhibitor treatment response in early breast cancer

Purpose: To determine genetic variants that could differentiate aromatase inhibitors (AIs) efficacy in early stage breast cancer. Experimental Design: We performed a stratified cox-proportional hazards analysis utilizing stratification factors and other covariates examining Breast Cancer Free Interval (BCFI) SNP-treatment interaction in Caucasian patients entered on the MA.27 trial. This analysis involved 4465 patients (2226 on anastrozole and 2239 on exemestane arm) including 245 (121 on anastrozole and 124 events on exemestane arm) with a breast event. Preclinical laboratory studies included luciferase activity, chromatin immunoprecipitation (ChIP) assay, and cell migration assays. Results: We identified 887 SNPs with a p-value &lt;1E-4 that could differentiate anastrozole from exemestane efficacy. We next used GTEx databases to determine whether these SNPs might be cis-eQTLs with nearby genes, and found that 95 SNPs were eQTL with 14 genes. Functional validation of SNP effect in these 14 genes on response to anastrozole and exemestane revealed that 3 SNPs showed genotype-dependent differences between anastrozole and exemestane. Two of the three SNPs, rs1877193 and rs6735923 located upstream of LY75 gene, were associated with higher LY75 gene expression. Our drug-SNP interaction GWAS showed that both SNPs were associated with better BCFI for exemestane compared to anastrozole (HR= 0.447, 0.458). The SNP rs62293499 located downstream of the GPR160 gene was also associated with better BCFI for exemestane compared to anastrozole (HR=0.433). Interestingly, all 3 SNPs were associated with worse outcome (shorter BCFI) in anastrozole treated patients (HR=1.39~1.58), but with longer BCFI in the exemestane treated patients (HR=0.58~0.71) based on our GWAS results in MA.27 trial. However, all 3 SNPs showed no association with BCFI if the two treatment arms, anastrozole and exemestane, were combined (HR=0.92~1.08). Consistently, cells with Ly75 or GPR160 SNP variant genotypes were more sensitive to exemestane compared to anastrozole. Mechanistically, the 3 SNPs regulated estrogen receptor-dependant LY75 and GPR160 expression. LY75 suppression induced epithelial-to-mesenchymal transition (EMT) in breast cancer cell lines, accompanied by increased migratory capacity in vitro. LY75 knockdown also resulted in predominant downregulation of functional pathways such as cell proliferation, while pathways associated with mesenchymal stimulation were generally increased. GPR160 knockdown also resulted in downregulation of cell proliferation pathways. Conclusions: This SNP genotype and AI treatment interaction clinical study revealed unique genetic variants that differentiate anastrozole and exemestane efficacy. The signals were lost in the GWAS analysis when anastrozole and exemestane were combined. Preclinical laboratory studies revealed novel functions of LY75 and GPR160 in breast cancer. These findings represent potential steps towards individualized AI therapy.

Citation Format: Junmei Cairns, James N. Ingle, Krishna R. Kalari, Lois E. Shepherd, Matthew J. Ellis, Paul E. Goss, Poulami Barman, Erin E. Carlson, Matthew P. Goetz, Richard M. Weinshilboum, Liewei Wang. The interaction between SNP genotype and aromatase inhibitor treatment response in early breast cancer [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr LB-103.

Cell-level somatic mutation detection from single-cell RNA sequencing

MOTIVATION

Both single-cell RNA sequencing (scRNA-seq) and DNA sequencing (scDNA-seq) have been applied for cell-level genomic profiling. For mutation profiling, the latter seems more natural. However, the task is highly challenging due to the limited input materials from only two copies of DNA molecules, while whole-genome amplification generates biases and other technical noises. ScRNA-seq starts with a higher input amount, so generally has better data quality. There exists various methods for mutation detection from DNA sequencing, it is not clear whether these methods work for scRNA-seq data.

RESULTS

Mutation detection methods developed for either bulk-cell sequencing data or scDNA-seq data do not work well for the scRNA-seq data, as they produce substantial numbers of false positives. We develop a novel and robust statistical method-called SCmut-to identify specific cells that harbor mutations discovered in bulk-cell data. Statistically SCmut controls the false positives using the 2D local false discovery rate method. We apply SCmut to several scRNA-seq datasets. In scRNA-seq breast cancer datasets SCmut identifies a number of highly confident cell-level mutations that are recurrent in many cells and consistent in different samples. In a scRNA-seq glioblastoma dataset, we discover a recurrent cell-level mutation in the PDGFRA gene that is highly correlated with a well-known in-frame deletion in the gene. To conclude, this study contributes a novel method to discover cell-level mutation information from scRNA-seq that can facilitate investigation of cell-to-cell heterogeneity.

AVAILABILITY AND IMPLEMENTATION

The source codes and bioinformatics pipeline of SCmut are available at https://github.com/nghiavtr/SCmut.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CYP2C9 and CYP2C19: Deep Mutational Scanning and Functional Characterization of Genomic Missense Variants

Single nucleotide variants in the open reading frames (ORFs) of pharmacogenes are important causes of interindividual variability in drug response. The functional characterization of variants of unknown significance within ORFs remains a major challenge for pharmacogenomics. Deep mutational scanning (DMS) is a high-throughput technique that makes it possible to analyze the functional effect of hundreds of variants in a parallel and scalable fashion. We adapted a "landing pad" DMS system to study the function of missense variants in the ORFs of cytochrome P450 family 2 subfamily C member 9 (CYP2C9) and cytochrome P450 family 2 subfamily C member 19 (CYP2C19). We studied 230 observed missense variants in the CYP2C9 and CYP2C19 ORFs and found that 19 of 109 CYP2C9 and 36 of 121 CYP2C19 variants displayed less than ~ 25% of the wild-type protein expression, a level that may have clinical relevance. Our results support DMS as an efficient method for the identification of damaging ORF variants that might have potential clinical pharmacogenomic application.

Concordance between predicted HLA type using next generation sequencing data generated for non-HLA purposes and clinical HLA type

We explored the feasibility of obtaining accurate HLA type using pre-existing NGS data not generated for HLA purposes. 83 exomes and 500 targeted NGS pharmacogenomic panels were analyzed using Omixon HLA Explore, OptiType, and/or HLA-Genotyper software. Results were compared against clinical HLA genotyping. 765 (94.2%) Omixon and 769 (94.7%) HLA-Genotyper of 812 germline allele calls across class I/II loci and 402 (99.5%) of 404 OptiType class I calls were concordant to the second field (i.e. HLA-A*02:01). An additional 19 (2.3%) Omixon, 39 (4.8%) HLA-Genotyper, and 2 (0.5%) OptiType allele calls were first field concordant (i.e. HLA-A*02). Using Omixon, four alleles (0.4%) were discordant and 24 (3.0%) failed to call, while 4 alleles (0.4%) were discordant using HLA-Genotyper. Tumor exomes were also evaluated and were 85.4%, 91.6%, and 100% concordant (Omixon and HLA-Genotyper with 96 alleles tested, and Optitype with 48 class I alleles, respectively). The 15 exomes and 500 pharmacogenomic panels were 100% concordant for each pharmacogenomic allele tested. This work has broad implications spanning future clinical care (pharmacogenomics, tumor response to immunotherapy, autoimmunity, etc.) and research applications.

Frequency of MicroRNA Response Elements Identifies Pathologically Relevant Signaling Pathways in Triple-Negative Breast Cancer

Complex interactions between mRNAs and microRNAs influence cellular functions. The mRNA-microRNA interactions also determine the post-transcriptional availability of mRNAs and unbound microRNAs. MicroRNAs binds to one or more microRNA response elements (MREs) located on the 3′UTR of mRNAs. In this study, we leveraged MREs and their frequencies in cancer and matched normal tissues to obtain insights into disease-specific interactions between mRNAs and microRNAs. We developed a bioinformatics method “ReMIx” that utilizes RNA sequencing (RNA-Seq) data to quantify MRE frequencies across the transcriptome. We applied ReMIx to triple-negative (TN) breast cancer tumor-normal adjacent pairs and identified MREs specific to TN tumors. ReMIx identified candidate mRNAs and microRNAs in the MAPK signaling cascade. Further analysis of MAPK gene regulatory networks revealed microRNA partners that influence and modulate MAPK signaling. In conclusion, we demonstrate a novel method of using MREs in the identification of functionally relevant mRNA-microRNA interactions in TN breast cancer.

•

Bioinformatics method ReMIx identify differential microRNA response rlements (MRE)

•

Tumor-specific MREs frequency observed in triple-negative breast cancer (TNBC)

•

MRE analysis identify MAPK signaling genes as therapeutic target for TNBC

•

MREs frequency can be used to identify pathologically relevant pathways