Skip to main content

Evaluation of invasive breast cancer samples using a 12-chemokine gene expression score: correlation with clinical outcomes



A unique 12-chemokine gene expression score (CS) accurately predicted the presence of tumor-localized, ectopic lymph node-like structures (TL-ELNs) and improved overall survival (OS) in primary colorectal cancer and metastatic melanoma. We analyzed the correlation between CS, clinicopathological variables, molecular data, and 366 survival in Moffitt Cancer Center’s Total Cancer Care (TCC) patients with non-metastatic breast cancer.


Affymetrix gene expression profiles were used to interrogate the CS by the principal component method. Breast tumors were classified as high or low score based on median split, and correlations between clinicopathologic variables, PAM50 molecular subtype, and ELN formation were analyzed using the TCC dataset. Differences in overall survival (OS) and recurrence-free survival (RFS) in the larger KM Plot breast cancer public datasets were compared using Kaplan-Meier curves.


We divided the Total Cancer Care (TCC) breast cancer patients into two groups of high or low CS. Mean CS was 0.24 (range, 2.2–2.1). Patients with higher CS were more likely to be white (172 vs. 159; p = 0.03), had poorly differentiated tumors (112 vs. 59; p <0.0001), ER/PR negative (41 vs. 26) and HER2 positive (36 vs. 19; p = 0.001), and contain TL-ELNs. Higher CS scores were also seen in the basal and HER2+ molecular subtypes. In the KM Plot breast cancer datasets higher CS patients demonstrated superior OS (HR = 0.73, p = 0.008) and RFS (HR 0.76, p = <0.0001), especially in basal and HER2+ patients.


High CS breast tumors tend to be higher grade, basal or HER2+, and present more frequently in Caucasians. However, this group of patients also shows the presence of TL-ELNs within the tumor microenvironment and has better survival outcomes. The CS is a novel tool that can identify breast cancer patients with tumors of a unique intratumoral immune composition and better prognosis. Whether or not the CS is a predictive response marker in breast cancer patients undergoing immunotherapy remains to be determined.


Breast cancer represents 14.0% of all new cancer cases in the United States, with 231,840 new cases and an estimated 40,290 deaths in 2015, comprising 6.8% of all cancer deaths [1]. Local and systemic treatments including surgery, chemotherapy, radiation, and endocrine therapy have all improved outcomes significantly for breast cancer patients [2]. However, the number of patients relapsing despite these treatments requires development of novel treatment modalities. One such modality that is garnering more attention recently is the use of immunotherapies [3]. Understanding how the various types of breast cancer interact with the immune system is important in informing us how to effectively utilize promising immune-oncology agents.

Although breast cancer is not perceived as a particularly immunogenic tumor when compared with melanoma and renal cell carcinoma as examples, molecular profiling of breast tumors has revealed that a subset demonstrate a high level of immunoregulatory gene activation [4]. Multiple investigators have reported that tumor-infiltrating lymphocytes and certain gene expression profiles related to immune signaling appear to have prognostic and/or predictive implications for breast cancer [especially the human epithelial growth factor receptor 2-positive (HER2+) type] [5,6,7,8,9]. These studies highlight the potential importance of the immune response in breast cancer patient outcomes. However, there are distinct types of immune cell infiltrates that can have different effects on tumor behavior. Characterization of the underlying mechanisms regulating immune infiltration in breast tumors can elucidate the key determinants for a successful host anti-tumor immune response. Secretion of chemokines within the tumor microenvironment and how certain co-morbidities like diabetes can affect the tumor chemokine milieu have gained attention as important factors that shape tumor lymphocyte infiltration [10, 11].

Chemokines act as trafficking signals for various immune cells and are important in orchestrating the spatial distribution of the immune response in a host. They also can directly affect the growth and progression of cancer cells [12]. Utilizing gene expression profiles can provide a more global assessment of immune signaling and cell populations using in silico methods such as CIBERSORT [13]. Certain chemokines have been associated with formation of a specific type of well-organized immune infiltrate known as tumor-localized, ectopic lymph node-like structures (TL-ELNs) [14]. It is hypothesized that these ELNs represent potent chemokine signaling gradients in the tumor microenvironment that attracts not only T cells but also activated B cells responding to specific tumor-associated antigens presented by co-localized dendritic cells [14].

Coppola and associates identified a unique 12-chemokine (CCL2, CCL3, CCL4, CCL5, CCL8, CCL18, CCL19, CCL21, CXCL9, CXCL10, CXCL11, and CXCL13) gene expression signature (GES) from a metagene grouping with overwhelming enrichment for immune-related and inflammation-related genes in primary colorectal cancer [15]. Messina and associates subsequently interrogated the 12-chemokine GES score (CS) across genomic arrays of 14,492 distinct solid tumors (primary and metastatic) of different histologies using the Total Cancer Care (TCC) database [16]. They found that this CS accurately predicted the presence of TL-ELNs and showed an association with improved overall survival in stage IV melanoma.

Since little was known about the effect of these aforementioned 12 chemokines on the breast tumor microenvironment, we sought to explore the relationship between the CS, presence of TL-ELNs, molecular subtype, and patient outcome in annotated breast cancer samples.


Patient inclusion criteria

A complete description of the TCC biobanking program has been previously published [17]. A retrospective review was performed on selected female patients with stage I to III breast cancer who were diagnosed between 1988 and 2012 and received primary surgery at the Moffitt Cancer Center. Patients may have received adjuvant therapies at Moffitt or at other locations. Snap-frozen tumor specimens from initially resected primary breast tumors were used for the gene expression profiles. Of the 813 unique gene expression files available, we included those from breast primary tumors only for which full clinical information was available. Patients chosen for study had available clinical and follow-up data within Moffitt Cancer Center’s electronic medical record system along with a genomic expression profile of their primary tumor. We excluded patients who had received any form of neoadjuvant therapy and with de novo metastatic disease, resulting in 366 patients in total.

mRNA microarray analysis

The 12-chemokine score for these samples was extracted from the TCC database. The 12-chemokine score for TCC was calculated, in brief, as follows. Tumors from patients treated at the Moffitt Cancer Center were arrayed on modified Affymetrix HuRSTA-2a520709 GeneChips (Affymetrix, Santa Clara, CA, USA). Chips were normalized using iterative rank-order normalization (IRON) [18]. An RNA quality-related batch effect was identified in the resulting normalized data, which was removed by training a partial least squares (PLS) model [19] to the RNA integrity number (RIN) for each sample and then subtracting the first PLS component. The final 12-chemokine score across all tumors in TCC was calculated using the first component from a principal component analysis (PCA) model based on the 12-chemokine genes. The 13 following probe sets was used; CCL18: merck-NM_002988_at, CCL19: merck-NM_006274_at, CCL2: merck-NM_002982_at, CCL21: merck-NM_002989_at, CCL3: merck-D63785_x_at, CCL4: merck-NM_002984_at, CCL5: merck-NM_002985_at, CCL8: merck-NM_005623_at, CXCL10: merck-NM_001565_at, CXCL11: merck-NM_005409_at & merck2-NM_005409_at, CXCL13: merck-NM_006419_at, CXCL9: merck-NM_002416_at.

A two-sided Student t test and Bonferroni correction were used to test for differences in high (n = 183) and low (n = 183) gene expression in the 366 fully annotated specimens across the 12 genes in the chemokine gene expression signature. Principal component analysis was performed and visualized using Evince™. We also studied immune gene expression levels of B and T lymphocyte attenuator (BTLA), cluster of differentiation (CD)14, CD31, CD274 [programmed death ligand 1 or (PD-L1)], CD56, CD69, CTLA4, CXCL12, fibroblast activation protein (FAP), granzyme B, T cell immunoglobulin and mucin domain containing 3 (TIM3), indoleamine 2,3 dioxygenase (IDO1), interferon gamma (IFN-γ), interleukin (IL)2, IL10, IL4, IL6, Janus-associated kinase 1 (JAK1), lymphocyte activation gene 3 (LAG3), lymphocyte expansion molecule (LEM), homologous to lymphotoxin exhibits inducible expression and competes with HSV glycoprotein D for binding to herpesvirus entry mediator, a receptor expressed on T lymphocytes (LIGHT), MHC class I polypeptide-related sequence A (MICA), neural cell adhesion molecule 1 (NCAM1), nuclear factor kappa light chain enhancer of activated B cells (NFKB), nitrogen oxide synthase 1 (NOS1), programmed death 1 (PD1), perforin 1 (PRF1), signal transducer and activator of transcription 1 (STAT1), signal transducer and activator of transcription 3 (STAT3), and vascular endothelial growth factor A (VEGFA). We performed a two tailed Student’s t test with a Bonferroni correction for multiple testing, comparing mean expression levels of the above-noted immune genes between groups with high and low 12-chemokine gene expression scores.

Pathologic analysis of tissue sections

Histological sections corresponding to 28 cases (prepared from the mirror image of the portion of tumor submitted for the mRNA microarray analysis) were retrieved from the Moffitt Cancer Center Anatomic Pathology Division’s repository as a pilot analysis to study correlation between the CS and tumor-localized, ectopic lymphoid node-like structures (TL-ELNs). Half of the specimens were from the top 10th percentile 12-chemokine gene expression scores and half were from the bottom 10th percentile. All of the specimens were 10% formalin fixed and paraffin embedded. Random representative hematoxylin and eosin (H&E)-stained sections through all selected tissue blocks were evaluated for the presence or absence of TL-ELNs. To further characterize the TL-ELNs, tissue sections were stained using the avidin-biotin complex method with retrieval under high pH. Pre-diluted monoclonal antibodies to CD3, CD4, CD8, and CD20 (Ventana Medical Systems, Tucson, AZ, USA) were used for the manual morphometric analysis of TL-ELNs by brightfield microscopy. To ensure pathologic concordance, two pathologists at our institution reviewed the tissue sections. Scores of 0 to 3 were assigned based on the following features: 0 = no lymphoid infiltrate noted in slide, 1 = 1 group lymphoid infiltrate, 2 = 2 groups of lymphoid infiltrate, and 3 = 3 or more groups of lymphoid infiltrate. Both pathologists were blinded as to the 12-chemokine gene expression scores of the individual samples. Clinical information was accessible only to the principal investigator and authorized collaborators, and all samples were anonymously coded before analysis. The Fisher’s exact test was used to test the association of 12-chemokine expression scores and H&E staining scores. The McNemar test was used to analyze the strength of agreement between the scoring methods of the two pathologists.

Patient variables and outcomes analyses

We compared clinical and pathological factors of patients with low versus high 12-chemokine gene expression scores calculated by principal component analysis (determined by median split). Correlation between clinicopathologic factors and the 12-chemokine gene expression score was tested using the chi-square test with the exact method using Monte Carlo estimation. Kaplan-Meier curves were created for both overall survival (OS) and recurrence-free survival (RFS), and log-rank tests were used to compare 12-chemokine gene expression scores and 12-chemokine gene expression scores stratified by antibody status [estrogen receptor (ER) or progesterone receptor (PR) positive], ER and PR negative (with HER2 negative or missing), and HER2 positive. Multivariable survival models were fit using Cox proportional hazards model. Final models were chosen using backward selection, with a removal alpha of 0.05. All p values were two-sided unless otherwise stated and considered statistically significant at the 0.05 level. The final multivariate survival model incorporated age, pathologic stage, and ER status based on this criteria. All statistical analyses were performed using SAS (version 9.4; SAS Institute; Cary, NC, USA). A second log rank Kaplan Meier OS and RFS analysis was done on the larger KM Plot breast cancer dataset [20] due to the small number of events within the TCC cohort using the corresponding 12-chemokine gene probeset (216598_s_at (CCL2), 205114_s_at (MGC12815), 204103_at (CCL4), 204655_at (TCP228), 214038_at (CCL8), 210072_at (CCL19), 204606_at (CCL21), 203915_at (Humig), 204533_at (CXCL10), 211122_s_at (SCYB9B), 205242_at (CXCL13)) on the available Affymetrix arrays. In all KM Plot analyses the high versus low chemokine patients were split by mean upper tertile expression of the 12 chemokine genes computed across the entire dataset with equal weighting. Intrinsic subtypes in the KM Plot dataset are based on St. Gallen criteria using estrogen receptor 1 (ESR1), HER2, and antigen KI-67 (KI67). Prediction analysis of microarray 50 (PAM50) was used to classify breast cancer molecular subtypes within the TCC dataset (11). The estrogen, progesterone, and HER2 receptor status of the samples were classified by TCC pathologists using College of American Pathologists (CAP) criteria and obtained from the TCC database. Estrogen and progesterone receptors were positive if >1% immunohistochemical (IHC) staining was noted, and HER2 status was positive if 3+ by IHC or with a positive ratio by in situ hybridization assay. GraphPad Prism software (GraphPad Software, San Diego, CA, USA) was used to compare the mean chemokine gene expression scores with a Tukey test to obtain 95% confidence intervals between the molecular subtypes and receptor statuses by immunohistochemistry.


Our study included 366 TCC patients who met our inclusion criteria, with 183 patients in the low and 183 patients in the high 12-chemokine gene expression score groups. Low and high 12-chemokine gene expression score groups were compared regarding patient demographics, tumor characteristics, treatment variables, and survival status. The 12-chemokine gene expression scores ranged from -2.2 to 2.1 (median 0.24). The median age at diagnosis was 54.5 years (range, 24–90 years). The median follow-up was 66.3 months (range, 2.5–212.9 months). The chosen patient population was predominantly white (331 vs. 35), with 284 patients (77.6%) having ductal histology. There were higher numbers of moderately and poorly differentiated tumors, 148 patients (41.7%) and 171 patients (48.2%), respectively. Most of the chosen patients were ER or PR positive (212; 63.5%), followed by 67 patients (20.1%) who were negative for ER, PR, and HER2 or missing HER2 positivity. Fifty-five patients (16.5%) were HER2 positive. Patients mainly were diagnosed with stage II disease (176 patients; 51.8%). Of 366 patients, 314 (86.5%) received adjuvant systemic therapy after surgery, 71.9% were alive at time of data collection, and 83 patients (23.4%) had disease recurrence.

Patients in the high 12-chemokine gene expression group (versus low 12-chemokine gene expression group) were more likely to be Caucasian (172 vs. 159 patients; p = 0.0298), had higher rates of poorly differentiated/high grade tumors (112 vs. 59 patients; p < 0.0001), and were more likely to be ER/PR negative (41 vs. 26 patients) and HER2 positive (36 vs. 19 patients) (p = 0.001) (Table 1). When we compared 12-chemokine gene expression score with PAM50 molecular subtype, higher score correlated more with basal and HER2-positive subtypes (Fig. 1a). Based on receptor status by immunohistochemistry, higher 12-chemokine gene expression scores were associated with triple-negative breast cancer (p = 0.0007) and HER2-positive tumors (p = 0.0002) (Fig. 1b).

Table 1 Comparison of variables with 12-chemokine gene signature
Fig. 1
figure 1

Comparison of chemokine score with (a) molecular subtype and (b) immunohistochemistry (IHC). ER estrogen receptor, HER2 human epidermal growth factor receptor 2, TNBC triple-negative breast cancer

The analysis of overall and recurrence-free survival in the KM Plot dataset demonstrated that patients with high 12-chemokine gene expression tumors had superior RFS (HR = 0.85, p = 0.018) and OS (HR = 0.63, p = <0.0012). The RFS was superior in patients with high 12-chemokine gene expression in the basal subtype (HR = 0.51, p = <0.0001), HER2 subtype (HR = 0.57, p = 0.0085), and luminal B subtype (HR 0.73, p = 0.0054). There was no RFS difference noted in the luminal A subtype (Fig. 2). A similar survival analysis was performed on the smaller TCC dataset, which showed a trend toward improved RFS in the HER2 subtype (Additional file 1). The number of OS and RFS events in the TCC dataset limited the power to fully evaluate the 12-chemokine gene expression score in relation to outcomes and was primarily used to evaluate associations between the 12-chemokine gene expression scores and clinicopathologic variables.

Fig. 2
figure 2

KM Plot Kaplan-Meier survival curves of (a) RFS in all patients (b) OS in all patients (c) RFS in basal subtype (d) HER-2 subtype (e) luminal B subtype (f) luminal A subtype patients. Red line = high chemokine group, black line = low chemokine group

Of the 67 H&E tumor tissue slides analyzed histologically for immune cell infiltration, 28 were scored as 0 with absence of TL-ELNs and 39 were positive for TL-ELNs. Both the 12-chemokine gene expression score and immune cell staining score were associated with each other for both of the pathologists’ scores (p < 0.001). There were no significant differences between these two scoring methods (p = 0.052), and kappa strength of agreement of 0.6148 indicated substantial strength of agreement between the two pathologists (Tables 2 and 3). Immune cell infiltrate scores of 0 were noted in 31/34 (88.6%) and 28/28 (100%) of slides corresponding to low 12-chemokine gene expression score by each individual pathologist, respectively. Between 71 and 75% of tumors with high 12-chemokine gene expression scores were scored 1–3 for TL-ELNs on random sections evaluated for each tumor (Fig. 3). The immunohistochemistry stains of the TL-ELNs demonstrated perifollicular presence of CD3+ CD4+ and CD3+ CD8+ T cells with strong staining for CD20 centrally showing clustering of mature B cells (Fig. 4).

Table 2 Association of high and low chemokine scores and ELNs
Table 3 Scoring distribution between pathologists for ELNs
Fig. 3
figure 3

Representative examples of lymphoid aggregates. The highest 12-chemokine gene expression signature scored breast tumors (T) revealed peritumoral lymphocytic host response organized as one or more ectopic lymph node-like structures by H&E staining (arrows) (a). In contrast, the lowest 12-chemokine gene expression signature scored breast tumors were predominantly devoid of inflammatory infiltrate (b)

Fig. 4
figure 4

Analysis of primary breast tumors with H&E staining and IHC. Representative high-chemokine-scored breast slide revealed a marked peritumoral lymphocytic host response, organized as ectopic lymph node-like structures (ELN) by H&E staining (a) and by IHC (be). The immunohistochemical analysis of the lymphoid population highlighted the CD3+, CD4+, and CD8+ T cells (b, c and d) distributed in the parafollicular cortex or marginal zones and with some dispersion into the follicles while CD20 + B cells are concentrated in the center of the follicles (e). Magnification × 200

Gene expression levels of BTLA, CD274, CD69, CTLA-4, granzyme B, IDO, interferon gamma, IL10, IL2, IL6, LAG3, PD-1, PRF1, STAT1, LIGHT were all significantly higher in the 12-chemokine gene expression high group (Table 4). The enrichment of these genes indicates that the 12-chemokine gene expression score also identifies tumors with higher levels of an activated Th1-skewed cytotoxic T cell infiltrate.

Table 4 Comparison of immune gene expression between high versus low 12-gene chemokine signature scores in breast cancer


The increasing awareness surrounding the importance of the host immune response in determining breast cancer outcomes provides new opportunities to integrate this information into treatment algorithms. Efforts to systematically describe the immune response in breast cancer by entities such as the TIL working group are critical to implementing this new system in the clinic [21]. However, given the complexity of the immune response and the need to personalize immunotherapy, it is becoming prudent to use molecular markers to dissect out what immune regulatory pathways are active in a given patient’s tumor [22]. The data presented herein indicate that certain chemokine genes can identify breast tumors enriched for tumor-localized, ectopic lymph node-like structures, and potentially provide a causal mechanism for why the tumor is inflamed in this manner.

Our study demonstrates that a 12-chemokine gene expression signature can identify a group of breast cancers with more favorable long-term outcomes. This is despite the fact that this group also contains greater number of tumors with traditionally adverse pathologic factors such as higher grade, ER negativity, and HER2 overexpression. In contrast to other immune infiltrate scoring methods, the chemokine score can provide a mechanistic explanation for why a particular tumor is forming TL-ELNs and exhibiting higher levels of activated T cell infiltrates. Another advantage of this approach is that chemokine scores can be obtained from limited core biopsies (i.e., prior to neoadjuvant therapy) while whole tissue sections would be required to histologically evaluate for the presence/absence of TL-ELNs in a tumor. An important question is what tumor-specific molecular features are conducive to the emergence of the high chemokine score phenotype. Future analyses should focus on analyzing other datasets combining RNA sequencing data that can provide information on mutational load and specific mutations or epitopes associated with a high chemokine gene expression score.

In our study, the TL-ELNs have the appearance of typical peripheral lymph nodes and are constructed of the necessary immune components, with CD3+, CD4+, and CD8+ T cells appearing in the parafollicular cortex or marginal zones and with some dispersion into the follicle and CD20+ B cells concentrated in the center of the follicle. The formation of these TL-ELNs is likely a different process compared to the lesser organized, dispersed infiltration of stromal tumor-infiltrating lymphocytes. Approximately 20% of invasive breast cancers contain perivascular TL-ELNs [20]. In particular, these infiltrates were associated with medullary breast cancers in one analysis, possibly accounting, in part, for its favorable prognosis [22].

Our study sheds light on the role of chemokine gene signaling in the tumor microenvironment and the formation of TL-ELNs, which potentially provides novel therapeutic opportunities. These may include manipulating TL-ELN-negative tumors to become TL-ELN-positive ones or isolating antibodies from the reactive B cell clones resident within TL-ELNs that potentially target tumor-associated antigens. Investigation of tumor chemokine gene expression scores in groups of breast cancer patients treated with checkpoint inhibitors and comparing its association with programmed death ligand 1 staining and clinical response is another possibility. In this respect, the chemokine score may prove useful to select patients for checkpoint blockade therapy. For chemokine-score-low breast tumors, increasing levels of key chemokines may ultimately prime those patients to respond more effectively to subsequent immunotherapies. Using immune gene expression signatures to personalize immunotherapy approaches could be critical in the future to maximizing clinical benefit in breast cancer patients.


The 12-gene chemokine score evaluated in our study was associated with ectopic lymph node formation in breast tumors, increased gene expression of immune signaling pathways, and improved outcomes. The chemokine score should be further explored as a prognostic factor and predictive marker for emerging immunotherapy approaches in breast cancer patients.



B and T lymphocyte attenuator


Cluster of differentiation


Chemokine score


Cytotoxic T lymphocyte associated protein 4


Ectopic lymph nodes


Estrogen receptor


Estrogen receptor 1


Fibroblast-activation protein


Gene expression signature


Hematoxylin and eosin


Human epithelial growth factor receptor 2


Indoleamine 2,3 dioxygenase


Interferon gamma






Janus-associated kinase 1


Antigen KI-67


Lymphocyte activation gene 3


Lymphocyte expansion molecule


Homologous to lymphotoxin exhibits inducible expression and competes with HSV glycoprotein D for binding to herpesvirus entry mediator, a receptor expressed on T lymphocytes


MHC class I polypeptide-related sequence A


Neural cell adhesion molecule 1


Nuclear factor kappa light chain enhancer of activated B cells


Nitrogen oxide synthase 1


Overall survival


Prediction analysis of microarray 50


Programmed death 1


Programmed death ligand 1


Progesterone receptor


Perforin 1


Recurrence-free survival


Signal transducer and activator of transcription 1


Signal transducer and activator of transcription 3


Total Cancer Care


T cell immunoglobulin and mucin domain containing 3


Vascular endothelial growth factor A


  1. SEER Cancer Statistics Review. 1975-2013. Accessed 12 Apr 2016.

  2. Gradishar WJ, Anderson BO, Balassanian R, Blair SL, Burstein HJ, Cyr A, Elias AD, Farrar WB, Forero A, Giordano SH, et al. Breast cancer, Version 1.2016. J Natl Compr Cancer Netw. 2015;13(12):1475–85.

    Article  CAS  Google Scholar 

  3. Emens LA, Tuohy VK, Stanton SE, Clarke E. Immunotherapy for breast cancer: is it feasible? Immunotherapy. 2015;7(11):1135–43.

    Article  CAS  Google Scholar 

  4. Ascierto ML, Kmieciak M, Idowu MO, Manjili R, Zhao Y, Grimes M, Dumur C, Wang E, Ramakrishnan V, Wang XY, et al. A signature of immune function genes associated with recurrence-free survival in breast cancer patients. Breast Cancer Res Treat. 2012;131(3):871–80.

    Article  CAS  PubMed  Google Scholar 

  5. Bedognetti D, Wang E, Marincola FM. Meta-analysis and metagenes: CXCL-13-driven signature as a robust marker of intratumoral immune response and predictor of breast cancer chemotherapeutic outcome. Oncoimmunology. 2014;3:e28727.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Sota Y, Naoi Y, Tsunashima R, Kagara N, Shimazu K, Maruyama N, Shimomura A, Shimoda M, Kishi K, Baba Y, et al. Construction of novel immune-related signature for prediction of pathological complete response to neoadjuvant chemotherapy in human breast cancer. Ann Oncol. 2014;25(1):100–6.

    Article  CAS  PubMed  Google Scholar 

  7. Perez EA, Thompson EA, Ballman KV, Anderson SK, Asmann YW, Kalari KR, Eckel-Passow JE, Dueck AC, Tenner KS, Jen J, et al. Genomic analysis reveals that immune function genes are strongly linked to clinical outcome in the North Central Cancer Treatment Group n9831 Adjuvant Trastuzumab Trial. J Clin Oncol. 2015;33(7):701–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Salgado R, Denkert C, Campbell C, Savas P, Nuciforo P, Aura C, de Azambuja E, Eidtmann H, Ellis CE, Baselga J, et al. Tumor-infiltrating lymphocytes and associations with pathological complete response and event-free survival in HER2-positive early-stage breast cancer treated with lapatinib and trastuzumab: a secondary analysis of the NeoALTTO Trial. JAMA Oncol. 2015;1(4):448–54.

    Article  PubMed  Google Scholar 

  9. Loi S, Michiels S, Salgado R, Sirtaine N, Jose V, Fumagalli D, Kellokumpu-Lehtinen PL, Bono P, Kataja V, Desmedt C, et al. Tumor infiltrating lymphocytes are prognostic in triple negative breast cancer and predictive for trastuzumab benefit in early breast cancer: results from the FinHER trial. Ann Oncol. 2014;25(8):1544–50.

    Article  CAS  PubMed  Google Scholar 

  10. Palacios-Arreola MI, Nava-Castro KE, Castro JI, Garcia-Zepeda E, Carrero JC, Morales-Montor J. The role of chemokines in breast cancer pathology and its possible use as therapeutic targets. J Immunol Res. 2014;2014:849720.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Wintrob ZA, Hammel JP, Nimako GK, Fayazi ZS, Gaile DP, Davis EE, Forrest A, Ceacareanu AC. Data report on inflammatory C-C chemokines among insulin-using women with diabetes mellitus and breast cancer. Data Brief. 2017;11:446–58.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Chow MT, Luster AD. Chemokines in cancer. Cancer Immunol Res. 2014;2(12):1125–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Coppola D, Mule JJ. Ectopic lymph nodes within human solid tumors. J Clin Oncol. 2008;26(27):4369–70.

    Article  PubMed  Google Scholar 

  15. Coppola D, Nebozhyn M, Khalil F, Dai H, Yeatman T, Loboda A, Mule JJ. Unique ectopic lymph node-like structures present in human primary colorectal carcinoma are identified by immune gene array profiling. Am J Pathol. 2011;179(1):37–45.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Messina JL, Fenstermacher DA, Eschrich S, Qu X, Berglund AE, Lloyd MC, Schell MJ, Sondak VK, Weber JS, Mule JJ. 12-Chemokine gene signature identifies lymph node-like structures in melanoma: potential for patient selection for immunotherapy? Sci Rep. 2012;2:765.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Fenstermacher DA, Wenham RM, Rollison DE, Dalton WS. Implementing personalized medicine in a cancer center. Cancer J. 2011;17(6):528–36.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Welsh EA, Eschrich SA, Berglund AE, Fenstermacher DA. Iterative rank-order normalization of gene expression microarray data. BMC Bioinformatics. 2013;14:153.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Wold S, Ruhe A, Wold H, Dunn IWJ. The collinearity problem in linear regression. The partial least squares (PLS) approach to generalized inverses. SIAM J Sci Stat Comput. 1984;5(3):735–43.

    Article  Google Scholar 

  20. Gyorffy B, Lanczky A, Eklund AC, Denkert C, Budczies J, Li Q, Szallasi Z. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat. 2010;123(3):725–31.

    Article  PubMed  Google Scholar 

  21. Salgado R, Denkert C, Demaria S, Sirtaine N, Klauschen F, Pruneri G, Wienert S, Van den Eynden G, Baehner FL, Penault-Llorca F, et al. The evaluation of tumor-infiltrating lymphocytes (TILs) in breast cancer: recommendations by an International TILs Working Group 2014. Ann Oncol. 2015;26(2):259–71.

    Article  CAS  PubMed  Google Scholar 

  22. Mellman I, Hubbard-Lucey VM, Tontonoz MJ, Kalos MD, Chen DS, Allison JP, Drake CG, Levitsky H, Lonberg N, van der Burg SH, et al. De-risking immunotherapy: report of a consensus workshop of the Cancer Immunotherapy Consortium of the Cancer Research Institute. Cancer Immunol Res. 2016;4(4):279–88.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank Rasa Hamilton (Moffitt Cancer Center) for editorial assistance.


This work has been supported in part by the Cancer Informatics Core Facility, Tissue Core Facility, and Total Cancer Care Initiative and Collaborative Data Services Core at the H. Lee Moffitt Cancer Center & Research Institute, an NCI-designated Comprehensive Cancer Center (P30-CA076292). This work was supported, in part, by gifts from benefactors Mrs. Julie Wooley and Mrs. Jackie Lewis and family, as well as the V Foundation.

Availability of data and materials

The dataset supporting the conclusions of this article are included within the article (and its additional file(s)).

Authors’ contributions

SP participated in manuscript authorship, data collection, and data analysis/interpretation. VR participated in data collection. LC, CC, and AB participated in manuscript authorship, data collection, biostatistical/bioinformatics and data analysis/interpretation. DC and FK participated in specimen collection, and pathology analysis/interpretation. JM participated in manuscript authorship, and data analysis/interpretation. HS participated in study design, manuscript authorship, and data collection/analysis/interpretation. All authors reviewed and approved the final manuscript.

Competing interests

James Mulé is a member of the following External Advisory Boards: Celgene Corp., Antigen Express, Inc., OncoPep, Inc., Lion Biotechnologies, Inc., and Vault Pharma, Inc. Dr. Hatem Soliman serves on Advisory Boards for Celgene, and Eli Lilly.

Consent for publication

Not applicable.

Ethics approval and consent to participate

All tumor and patient information were retrieved from the Moffitt Total Cancer Care (TCC) Biorepository from prospectively consented patients in accordance with policies of the University of South Florida Institutional Review Board.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Hatem H. Soliman.

Additional file

Additional file 1:

Outcomes (OS, RFS for entire cohort and RFS for HER2+) plus multivariate regression analysis of overall survival within the TCC dataset. (DOCX 119 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Prabhakaran, S., Rizk, V.T., Ma, Z. et al. Evaluation of invasive breast cancer samples using a 12-chemokine gene expression score: correlation with clinical outcomes. Breast Cancer Res 19, 71 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: