DNA barcoded competitive clone-initiating cell analysis reveals novel features of metastatic growth in a cancer xenograft model

A problematic feature of many human cancers is a lack of understanding of mechanisms controlling organ-specific patterns of metastasis, despite recent progress in identifying many mutations and transcriptional programs shown to confer this potential. To address this gap, we developed a methodology that enables different aspects of the metastatic process to be comprehensively characterized at a clonal resolution. Our approach exploits the application of a computational pipeline to analyze and visualize clonal data obtained from transplant experiments in which a cellular DNA barcoding strategy is used to distinguish the separate clonal contributions of two or more competing cell populations. To illustrate the power of this methodology, we demonstrate its ability to discriminate the metastatic behavior in immunodeficient mice of a well-established human metastatic cancer cell line and its co-transplanted LRRC15 knockdown derivative. We also show how the use of machine learning to quantify clone-initiating cell (CIC) numbers and their subsequent metastatic progeny generated in different sites can reveal previously unknown relationships between different cellular genotypes and their initial sites of implantation with their subsequent respective dissemination patterns. These findings underscore the potential of such combined genomic and computational methodologies to identify new clonally-relevant drivers of site-specific patterns of metastasis.

Luminal androgen receptor breast cancer subtype and investigation of the microenvironment and neoadjuvant chemotherapy response

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype with low overall survival rates and high molecular heterogeneity; therefore, few targeted therapies are available. The luminal androgen receptor (LAR) is the most consistently identified TNBC subtype, but the clinical utility has yet to be established. Here, we constructed a novel genomic classifier, LAR-Sig, that distinguishes the LAR subtype from other TNBC subtypes and provide evidence that it is a clinically distinct disease. A meta-analysis of seven TNBC datasets (n = 1086 samples) from neoadjuvant clinical trials demonstrated that LAR patients have significantly reduced response (pCR) rates than non-LAR TNBC patients (odds ratio = 2.11, 95% CI: 1.33, 2.89). Moreover, deconvolution of the tumor microenvironment confirmed an enrichment of luminal epithelium corresponding with a decrease in basal and myoepithelium in LAR TNBC tumors. Increased immunosuppression in LAR patients may lead to a decreased presence of cycling T-cells and plasma cells. While, an increased presence of myofibroblast-like cancer-associated cells may impede drug delivery and treatment. In summary, the lower levels of tumor infiltrating lymphocytes (TILs), reduced immune activity in the micro-environment, and lower pCR rates after NAC, suggest that new therapeutic strategies for the LAR TNBC subtype need to be developed.

Abstract P3-05-06: Prognostic value of the neutrophil-to-lymphocyte ratio (NLR) and its relation to stromal tumor infiltrating lymphocytes (sTILs) in triple negative breast cancer (TNBC)

Background: While TNBC remains the most aggressive type of breast cancer (BC), substantial heterogeneity in biology and outcomes exists among TNBC subtypes. Historically, risk stratification of TNBC has been based on anatomic factors such as tumor size, nodal involvement and presence of distant metastases. However, these features alone fail to accurately predict outcomes. Tumor immune infiltration (sTILs) and distribution of immune cell subsets in the perip heral blood (NLR) have emerged as variables reported to be associated with outcomes in TNBC. We sought to evaluate whether NLR and sTILs provided independent prognostic information in TNBC.

Methods: From a cohort of 9,982 women who underwent BC surgery at Mayo Clinic, Rochester, MN between Jan 1985 and Dec 2012, we identified 605 centrally-confirmed TNBC tumors. Patients (pts) with prior BC, bilateral BC, non-invasive disease, stage IV, neoadjuvant therapy, endocrine therapy, or adenoid cystic histology were excluded. For eligible tumors, clinical and pathologic variables were evaluated, including peripheral blood NLR and central assessment of sTILs per the 2014 International TILs Working Group recommendations. We calculated the Pearson correlation coefficient (PCC) between NLR and sTILs and constructed Cox Proportional Hazards Models to evaluate their association with invasive-disease free (IDFS) and overall survival (OS). NLR and sTILs were both analyzed as continuous variables.

Results: Most pts had T1-2 (95%) and N0-1 disease (86%). Median OS follow-up was 10.6yrs. Median IDFS was 12yrs (95%CI 10.2-16.7) and median OS was 18.8yrs (95%CI 15.6-20.8). NLR and sTILs were available in 408 and 599 pts, respectively. The median NLR and sTIL content were 2.29 (0.14-10.50) and 20% (0-90%), respectively. NLR and sTILs were poorly correlated (PCC 0.0237). On univariate analysis (UVA), a higher NLR was associated with worse IDFS (HR 1.13; 95%CI 1.02-1.26, p=0.02) and OS (HR 1.17; 95%CI 1.04-1.31, p=0.01). Each 1% increment in sTILs was associated with improved IDFS (HR 0.99; 95%CI 0.98-0.99, p<0.001) and OS (HR 0.99, 95%CI 0.98-1.00, p<0.001). Among pts with high sTILs (≥20%), a higher NLR remained significantly associated with worse IDFS (HR 1.21; 95%CI 1.05-1.38, p=0.007) and OS (HR 1.25; 95%CI 1.09-1.44, p=0.001). In contrast, among pts with low sTILs (<20%), NLR was not associated with IDFS (HR 1.07; 95%CI 0.89-1.28, p=0.49) or OS (HR 1.07; 95%CI 0.88-1.30, p=0.49). The interaction test between NLR and sTILs did not reach statistical significance. A multivariate analysis (MVA; including age, menopausal status, histologic subtype, grade, tumor size, nodal stage, Ki-67, NLR, sTILs, adjuvant chemotherapy, type of surgery and adjuvant radiation) showed that sTILs remained independently associated with IDFS (HR 0.99, 95%CI 0.97-1.0, p=0.019) and OS (HR 0.99, 95% CI 0.97-1.0, p=0.044), whereas NLR did not.

Conclusions: A lower NLR and a higher sTIL content were each associated with improved IDFS and OS among pts with nonmetastatic TNBC on UVA. However, when evaluated on a MVA, only sTILs remained independently associated with IDFS and OS. Our data suggest that the effect of sTILs on outcomes may not be modified by the NLR.

Citation Format: Leon-Ferre RA, Polley M-Y, Liu H, Gilbert J, Cafourek V, Hillman D, Negron V, Boughey JC, Liu MC, Ingle JN, Kalari K, Couch FJ, Visscher DW, Goetz MP. Prognostic value of the neutrophil-to-lymphocyte ratio (NLR) and its relation to stromal tumor infiltrating lymphocytes (sTILs) in triple negative breast cancer (TNBC) [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P3-05-06.

Abstract P1-06-07: Mayo clinic TNBC outcome calculator: A clinical calculator to predict disease relapse and survival in women with triple-negative breast cancer

Purpose: Triple negative breast cancer (TNBC) is an aggressive breast cancer subtype with substantial risks of disease recurrence. While cytotoxic chemotherapy is commonly administered and reduces recurrence, disease outcomes vary considerably and few prognostic tools are available for risk stratification for TNBC patients. We constructed and validated clinical calculators for invasive-disease free survival (IDFS) and overall survival (OS) for TNBC and compared their performance against AJCC-based models which include only tumor size and nodal status.

Methods: From a surgical cohort of 9,982 patients who underwent breast cancer surgery at Mayo Clinic between January 1985 and December 2012, 605 centrally reviewed TNBC patients were identified and used to construct Cox models for IDFS and OS. Patients treated with neoadjuvant chemotherapy were excluded. Variables considered included age, menopausal status, tumor size, nodal status, Nottingham grade, type of breast surgery (mastectomy vs. lumpectomy), adjuvant radiation therapy, adjuvant chemotherapy, Ki67, stromal tumor infiltrating lymphocytes (sTILs), and neutrophil-to-lymphocyte ratio (NLR). Missing values were imputed using single imputation with all variables (including outcomes) included in the imputation model. Backward step-down procedure was used for model selections. The final models were internally validated for calibration and discrimination using bootstrapping methods and compared with AJCC-based models.

Results: For both IDFS and OS, higher sTIL's, less extensive nodal involvement, use of adjuvant chemotherapy, and lower NLR were significant predictors of improved clinical outcomes. Premenopausal status and younger age were additionally predictive of improved IDFS and OS, respectively. Models for IDFS and OS have good calibration and are associated with bias-corrected C-indices of 0.68 and 0.71, respectively, as compared with C-indices of 0.59 and 0.62 for AJCC-based models.

Conclusions: Our data indicate that a clinical calculator that includes sTIL's, NLR, menopausal status, age, nodal involvement as well as chemotherapy use can provide significantly greater prediction of clinical risk than tumor size and nodal status alone. These tools may be used to identify TNBC patients at elevated risk of disease relapse and to aid physician's communication with patients regarding their long-term disease outlook and planning treatment strategies. External validation is required to further evaluate broader applicability of this tool, which was developed utilizing a single-institutional experience.

Citation Format: Polley M-YC, Leon-Ferre RA, Liu H, Gilbert J, Cafourek V, Hillman DW, Negron V, Boughey JC, Liu MC, Ingle JN, Kalari K, Couch F, Visscher DW, Goetz MP. Mayo clinic TNBC outcome calculator: A clinical calculator to predict disease relapse and survival in women with triple-negative breast cancer [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P1-06-07.

Abstract PD1-04: CSMD1 SNPs selectively affect anastrozole response in postmenopausal breast cancer patients

BACKGROUND: Based on prospective clinical trials, there is no evidence for differences in efficacy between the 3 aromatase inhibitors (AIs) anastrozole, exemestane, and letrozole. The purpose of this study was to identify germline genetic variants associated with response to AIs and to help identify novel mechanisms associated with drug disease efficacy.

METHODS: A genome-wide association study (GWAS) was performed for 624 patients (Steroids 2015;99:32-38) to identify SNPs associated with estrogen level change in women with estrogen receptor (ER) positive breast cancer treated with anastrozole. Replication of associated SNPs was performed in a GWAS from the MA.27 trial that compared adjuvant anastrozole and exemestane treatment of post-menopausal women with ER+ breast cancer. Functional studies were subsequently performed to determine SNP effects and underlying mechanisms.

RESULTS: Our initial GWAS identified SNPs within CSMD1 that were associated with changes in estrogen levels during anastrozole therapy. An additional SNP in CSMD1 was also associated with breast cancer events in CCTG MA.27. Functionally, we showed that CSMD1 regulates CYP19 expression in a SNP-, and in an anastrozole- dependent fashion. These phenomena were not observed for either letrozole or exemestane. In MA.27, an anastrozole- specific effect was also seen with the minor allele having a protective effect on time to distant metastasis (HR=0.49, p=0.00259), but this was not the case for exemestane (HR=0.71, p=0.111). Our in vitro functional studies indicated that overexpression of CSMD1 sensitized anastrozole or letrozole resistant cells to anastrozole but not to the other two AIs. The SNP in CSMD1 that was associated with increased CSMD1 and CYP19 expression levels increased anastrozole sensitivity, but not letrozole or exemestane in lymphoblastoid cell lines (LCLs) homozygous for either WT or variant CSMD1 SNP genotypes. Based on these observations, we explored whether anastrozole has additional mechanisms beyond its function as a CYP19 inhibitor. Utilizing an estrogen response element (ERE) luciferase reporter assay in a CYP19 CRISPR knockout breast cancer T47D cell line and a surface plasmon resonance (SPR) assay, we found that anastrozole can also function as an ERα agonist, and can bind to, and result in, proteasome dependent ERα degradation, especially in the presence of E2. Treatment of these CYP19 CRISPR knockout cells with anastrozole in the presence of increasing concentrations of E2 results in greater sensitivity compared with anastrozole alone, while the addition of E2, as expected, does not improve letrozole or exemestane sensitivity. These same observations were also seen in letrozole and anastrazole resistant cells.

CONCLUSIONS: Our findings suggest that anastrozole might be more effective than letrozole or exemestane in patients with the CSMD1 SNP. Furthermore, anastrozole can function as an ERα agonist, binding to ERα and resulting in its degradation, especially in the presence of E2. These findings should help to make it possible to develop precision endocrine therapies for women who are candidates for AIs.

Citation Format: Cairns J, Ingle J, Dudenkov T, Kalari K, Buzdar A, Kubo M, Robson M, Ellis M, Goss P, Shepherd L, Goetz M, Weinshilboum R, Wang L. CSMD1 SNPs selectively affect anastrozole response in postmenopausal breast cancer patients [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr PD1-04.

TSPAN5, ERICH3 and selective serotonin reuptake inhibitors in major depressive disorder: pharmacometabolomics-informed pharmacogenomics

Millions of patients suffer from major depressive disorder (MDD), but many do not respond to selective serotonin reuptake inhibitor (SSRI) therapy. We used a pharmacometabolomics-informed pharmacogenomics research strategy to identify genes associated with metabolites that were related to SSRI response. Specifically, 306 MDD patients were treated with citalopram or escitalopram and blood was drawn at baseline, 4 and 8 weeks for blood drug levels, genome-wide single nucleotide polymorphism (SNP) genotyping and metabolomic analyses. SSRI treatment decreased plasma serotonin concentrations (P<0.0001). Baseline and plasma serotonin concentration changes were associated with clinical outcomes (P<0.05). Therefore, baseline and serotonin concentration changes were used as phenotypes for genome-wide association studies (GWAS). GWAS for baseline plasma serotonin concentrations revealed a genome-wide significant (P=7.84E-09) SNP cluster on chromosome four 5' of TSPAN5 and a cluster across ERICH3 on chromosome one (P=9.28E-08) that were also observed during GWAS for change in serotonin at 4 (P=5.6E-08 and P=7.54E-07, respectively) and 8 weeks (P=1.25E-06 and P=3.99E-07, respectively). The SNPs on chromosome four were expression quantitative trait loci for TSPAN5. Knockdown (KD) and overexpression (OE) of TSPAN5 in a neuroblastoma cell line significantly altered the expression of serotonin pathway genes (TPH1, TPH2, DDC and MAOA). Chromosome one SNPs included two ERICH3 nonsynonymous SNPs that resulted in accelerated proteasome-mediated degradation. In addition, ERICH3 and TSPAN5 KD and OE altered media serotonin concentrations. Application of a pharmacometabolomics-informed pharmacogenomic research strategy, followed by functional validation, indicated that TSPAN5 and ERICH3 are associated with plasma serotonin concentrations and may have a role in SSRI treatment outcomes.

Abstract P5-04-01: ERβ elicits tumor suppressive effects in triple negative breast cancer through the induction of cystatins and suppression of TGFβ signaling

Background: Triple negative breast cancer (TNBC) accounts for approximately 20% of all breast cancer diagnoses. Clinical management of TNBC is limited to surgery, chemotherapy and radiation due to lack of estrogen receptor alpha and HER2 expression. Recently, we have shown that approximately 40% of TNBCs express estrogen receptor beta (ERβ) and have begun to explore the possibility that this receptor could be utilized as a novel therapeutic target for this disease.

Methods: To examine the biological functions of ERβ in TNBC, novel ERβ expressing TN cell lines (MDA-MB-231 and Hs578T) were developed. In vitro experiments were employed to determine alterations in the global gene expression profiles, biological pathways, proliferation rates, and cell cycle progression following estrogen or ERβ-specific agonist treatment. Cell line xenografts were also established in athymic ovariectomized nude mice to examine tumoral responses to ERβ targeting agents and to investigate gene and protein expression patterns as well as potential serum biomarkers indicative of therapeutic response. Additionally, using the resources of the Mayo Clinic Breast Cancer Genome Guided Therapy Study (BEAUTY), we have identified, and begun to analyze, ERβ+ and ERβ- patient derived xenografts (PDX) established from women with TNBC.

Results: Our studies have revealed that both estrogen and multiple ERβ-specific agonists elicit significant anti-proliferative effects in ERβ+ TNBC cells primarily through a G1/S phase cell cycle arrest. These anti-proliferative effects appear to be mediated by cystatins, a family of small secreted cysteine protease inhibitors which are highly induced following estrogen and ERβ-specific agonist treatment. Conditioned media isolated from estrogen or ERβ-specific agonist treated cells decreased the proliferation rates of multiple non-ERβ expressing cell lines; effects that were completely reversed when cystatins were depleted from the media. In addition, we have shown that activation of ERβ, and the subsequent induction of cystatin gene expression, leads to suppression of canonical TGFβ signaling through multiple mechanisms including suppression of TGFβR2 expression, induction of Smad7 expression and blockade of TGFβ ligand-mediated activation of this pathway both in vitro and in vivo. Finally, ERβ+ TNBC PDXs exhibit significantly decreased tumor growth rates in estrogen-treated mice compared to ERβ- TN breast tumors.

Conclusions: Our in vitro and in vivo data show that estrogen and ERβ-specific agonists elicit anti-cancer effects in ERβ+ TNBC. These effects appear to be mediated, in part, by cystatins through their inhibitory effects on canonical TGFB signaling, a pathway known to drive TNBC progression. Importantly, these data lay the foundation for studies aimed at examining the ability to therapeutically target ERβ in TNBC patients.

Citation Format: Reese J, Bruinsma E, Subramaniam M, Suman V, Pitel K, Kalari K, Yu J, Wang L, Goetz M, Ingle J, Hawse J. ERβ elicits tumor suppressive effects in triple negative breast cancer through the induction of cystatins and suppression of TGFβ signaling. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P5-04-01.