Tumor-infiltrating lymphocytes predict response to anthracycline-based chemotherapy in estrogen receptor-negative breast cancer
© West et al.; licensee BioMed Central Ltd. 2011
Received: 28 February 2011
Accepted: 8 December 2011
Published: 8 December 2011
Infiltration of breast tumors by tumor-infiltrating lymphocytes (TIL) has been associated with sensitivity to anthracycline-based chemotherapy. However, it is unclear whether this is true within the estrogen receptor-alpha (ER)-negative subset of breast tumors that frequently manifest high TIL levels.
The association of TIL with short-term and long-term clinical response to anthracycline-based therapy was assessed in two independent ER-negative breast cancer cohorts in which patients were categorized as TIL-high or TIL-low. We defined an eight-gene lymphocyte mRNA expression signature (including CD19, CD3D, CD48, GZMB, LCK, MS4A1, PRF1, and SELL) and used unsupervised hierarchical clustering to examine the association between TIL and short-term response to neoadjuvant chemotherapy in a previously published cohort of ER-negative tumors (n = 113). We also examined the association between TIL and long-term chemotherapeutic efficacy in a second cohort of ER-negative tumors (n = 255) with longer than 6 years of median follow-up by using tissue microarrays and immunohistochemistry (IHC) for detection of CD3, CD8, CD4, CD20, and TIA-1.
In patients with ER-negative tumors treated with neoadjuvant anthracycline-based chemotherapy, pathologic complete responses (pCRs) were achieved by 23 (74%) of 31 TIL-high patients and 25 (31%) of 80 TIL-low patients (odds ratio (OR), 6.33; 95% confidence interval (CI), 2.49 to 16.08; P < 0.0001). Multivariate logistic regression with standard clinicopathologic features demonstrated that only tumor size (P = 0.037) and TIL status (P = 0.001) were independent predictors of anthracycline response. In the second cohort, adjuvant anthracycline-based therapy was associated with increased disease-free survival (DFS) only in patients with high levels of intraepithelial CD3+ TIL (P = 0.0023). In contrast, outcomes after CMF treatment (cyclophosphamide, methotrexate, and fluorouracil) showed no association with CD3 status. In both cohorts, cytotoxic T-cells were the primary TIL subtype associated with anthracycline sensitivity. Finally, TIL significantly predicted anthracycline sensitivity for both the Her2-positive and triple-negative tumor phenotypes.
ER-negative breast cancers with high levels of TIL have heightened sensitivity to anthracycline-based chemotherapy, as assessed by the immediate response to neoadjuvant therapy and long-term outcome following adjuvant therapy. Investigations of TIL-based predictive tests to identify patients likely to benefit from anthracycline-based treatments are warranted.
The use of adjuvant systemic therapy for breast cancer has increased substantially in the last few decades, and is now a mainstay modality in clinics worldwide. Although adjuvant chemotherapy has significantly improved patient management, a fundamental limitation of this approach is that many patients fail to benefit from therapy, and physicians are frequently unable to predict the responses of individual patients to a given regimen. As such, much effort has been directed in recent years toward identifying clinical and biological predictive features to better tailor therapy to the needs of individual patients.
Anthracycline-based chemotherapy (regimens involving anthracyclines such as doxorubicin or epirubicin) has been used clinically for more than two decades and has largely supplanted first-generation regimens such as CMF . Anthracyclines are thought to exert their effects through a variety of mechanisms including intercalation of DNA, cross-linking of DNA to proteins, and generation of free radicals . Nevertheless, the precise mechanisms by which anthracyclines exert their therapeutic effects in vivo remain unclear, nor has significant progress been made in establishing predictive biomarkers to identify patients likely to derive benefit from anthracycline-based therapy [2, 3]. For example, aberrations of HER2 and TOP2A initially showed promise as anthracycline-response predictors, but have failed to demonstrate consistent clinical utility .
Attention has focused recently on the prognostic and predictive potential of antitumor immune responses, detected indirectly via gene expression signatures derived from tumor-infiltrating lymphocytes (TIL) or directly via immunohistochemical TIL staining. In breast cancer, TIL have been shown in several studies to correlate with favorable long-term prognosis, although primarily for hormone receptor-negative, Her2-positive, or high grade/highly proliferative lesions [4–10]. Denkert et al.  recently reported that TIL were associated with a favorable response to neoadjuvant anthracycline/taxane therapy in two large breast cancer cohorts, providing compelling evidence that TIL could potentially serve as predictive markers for anthracycline-based therapy. However, no such studies have been conducted specifically with the ER-negative breast cancer subset. This is important for two reasons: first, pathologic complete response to neoadjuvant chemotherapy, defined as the complete absence of invasive tumor cells in the breast and lymph nodes after treatment, occurs almost exclusively in tumors negative for ER [12–15]; second, ER-negative tumors typically feature higher levels of TIL than do ER-positive tumors [6, 11, 16–18]. Therefore, our objective in this study was to determine whether TIL correlate with sensitivity to anthracycline-based chemotherapy in ER-negative breast cancer. To do so, we examined two independent cohorts of ER-negative breast cancer cases. One cohort received anthracycline-based therapy in the neoadjuvant setting, permitting the assessment of short-term clinical responses to therapy. The second cohort included patients with long-term follow-up treated with anthracyclines in the adjuvant setting, allowing us to examine overall and disease-free survival rates after systemic therapy.
Materials and methods
Neoadjuvant chemotherapy cohort
Characteristics of all assessed patients in the EORTC and MBTB cohorts
EORTC (n= 113)
MBTB (n= 255)
Age at diagnosis
<2 cm to ≤5
Analysis of microarray data
To determine whether lymphocyte-associated gene expression correlated with pCR rates, we analyzed microarray data from the neoadjuvant (EORTC) cohort and obtained an eight-gene TIL signature by using the following steps. We first compiled a list of genes with known biological functions that are expressed at high levels by lymphocytes, based on information in the BioGPS gene portal (http://biogps.gnf.org). Expression data for these genes were compiled and z-normalized (that is, number of standard deviations from the mean). For genes with redundant probe sets, only those probes with the highest interquartile range were selected. To enrich for genes with expression patterns consistent with a lymphocyte cell of origin, expression values of all genes were Spearman correlated with those of CD247 (CD3-zeta chain) and MS4A1 (CD20), which were chosen as specific and prototypical T- and B-lymphocyte genes, respectively. Only genes that correlated to CD247 or MS4A1 with a Spearman R value of 0.7 or greater were retained, resulting in a list of 33 genes. With this gene set, cases were grouped by unsupervised hierarchical clustering using the centroid method and euclidean distance as the similarity metric (performed with Cluster 3.0). Heat maps and cluster dendograms were produced using Java TreeView.
We further refined the TIL signature by assessing differential expression of each gene in pCR versus RD cases and selecting genes with a Mann-Whitney t test P value less than 0.05. The final eight-gene signature included CD19, CD3D, CD48, GZMB, LCK, MS4A1, PRF1, and SELL. Data corresponding to this signature were used to cluster the EORTC cases, as described above. Two cases that failed to cluster with any of the primary centroids were considered outliers and discarded from further analysis, leaving a final assessed cohort size of 111 cases. Univariate associations of TIL status and pathologic parameters with pCR rates were calculated using two-sided Fisher's exact tests in Prism 5.0 (GraphPad, La Jolla, CA). Multivariate logistic regression models for prediction of pCR were constructed in SPSS Statistics 14 (SPSS, Chicago, IL).
Adjuvant chemotherapy cohort
To assess the relationship between TIL and long-term outcome after adjuvant chemotherapy, we analyzed TIL in primary tumor biopsies or resections taken before adjuvant therapy from a cohort of patients with ER-negative invasive ductal breast carcinomas registered by the Manitoba Breast Tumor Bank (MBTB). Cases were accrued between the years 1988 and 2000 and had a median follow-up time of 83 months. Clinical characteristics of the study cohort are summarized in Table 1 and Additional file 1 (Supplementary Table S1). All selected cases were determined to be ER-negative by a single central provincial clinical laboratory and subsequently managed as ER-negative cases at a single provincial cancer center, where adjuvant systemic therapy was administered to 61% of the cohort. Of these patients, 58% were treated with CMF (cyclophosphamide, methotrexate, fluorouracil), and 37% were treated with one of several anthracycline-based regimens (for example, AC, CAF, CEF, where A is Adriamycin/doxorubicin, and E is epirubicin). Disease-free survival (DFS) was defined as the time from surgery to the first instance of disease recurrence or disease-specific death, and overall survival as the time to death from all causes. The MBTB has approval from the Research Ethics Board, Faculty of Medicine, University of Manitoba, to collect, store, and distribute anonymous cases from its archive under the Canadian Tri-Council Policy Statement waiver of consent. The current study was conducted with approval from the University of British Columbia/British Columbia Cancer Agency Research Ethics Board.
Tissue microarray construction
An initial cohort of 255 ER-negative cases was selected on the basis of (a) ER-negative status defined by ligand-binding analysis of <10 fmol/mg protein; (b) invasive ductal carcinoma components occupying more than 20% of the tumor section, and (c) no prior therapy. To construct a tissue microarray (TMA), all cases were rereviewed to confirm and select areas for coring of corresponding blocks. Duplicate tissue cores (0.6 mm diameter) were taken from central cellular areas of each tumor with a tissue arrayer instrument (Beecher Instruments, Silver Spring, MD). Prior use and exhaustion of some tissue cores reduced the interpretable cohort size to approximately 160, depending on the TIL marker analyzed.
Immunohistochemistry and TMA scoring
Immunohistochemistry (IHC) was performed for CD3, CD8, CD4, CD20, and TIA-1 on deparaffinized sections from TMAs by using a Ventana Discovery XT autostainer (Ventana, Tucson, AZ). Ventana's standard CC1 protocol was used for antigen retrieval. Primary antibodies used are as follows (clone, supplier, animal source, concentration): CD3 (RM-9107, Lab Vision, rabbit, 1/150); CD8 (RM-9116, Lab Vision, rabbit, 1/100); CD4 (MS-1528, Lab Vision, mouse, 1/10); CD20 (polyclonal, catalogue no. RB-9013, Lab Vision, rabbit, 1/250); TIA-1 (TIA-1, Abcam, mouse, 1/50). TMA sections were incubated with primary antibodies for 60 minutes at room temperature followed by the appropriate cross-adsorbed, biotinylated secondary antibody (Jackson Immunoresearch, West Grove, PA) for 32 minutes. Antibodies were detected using the DABMap kit (Ventana). Slides were counterstained with hematoxylin (Ventana) and coverslipped manually with Cytoseal-60 (Richard Allan, Kalamazoo, MI).
IHC scoring was undertaken by using a microscope eyepiece grid to standardize the assessed area. In brief, duplicate cores of each immunostained tumor were reviewed at low magnification, and the core exhibiting a tumor/stroma ratio closest to 50:50 and the highest density of positive cells was selected. This core was then assessed at higher magnification (×20 objective) with a grid overlaid on the center of the core. Under a ×20 objective magnification, this grid defined an area of 0.56 mm2. The number of positive intraepithelial lymphocytes was quantified within the area of the grid (intraepithelial localization was defined as lymphocytes within tumor cell nests or in direct contact with tumor cells, consistent with the method used by Denkert et al. ). To account for variation in epithelial-stromal proportions between different samples, intraepithelial TIL levels were calculated by dividing the number of observed intraepithelial TIL by the fraction of grid area occupied by epithelium.
Statistical analysis of IHC data
All analyses were performed using Prism 5.0. Unless otherwise specified, median values for TIL markers were used as predetermined cut-points to define high versus low cases. Associations of TIL with pathologic features were evaluated with Fisher's exact test. Survival outcomes were assessed via Kaplan-Meier methods and compared using Log-rank tests. All statistical tests were two-sided, with significance established at P values less than 0.05.
Additional validation datasets
We assessed publically available gene-expression data from ER-negative breast cancers within three additional cohorts (accessible via the Gene Expression Omnibus website), the characteristics of which are as follows: GSE21974, a cohort of 32 breast cancer patients (14 ER-negative) treated with neoadjuvant epirubicin, cyclophosphamide, and docetaxel (not previously published); GSE19615 , a cohort of 115 breast cancers, 36 of which were ER-negative and treated with adjuvant anthracycline-based therapy (primarily AC and AC plus taxol); and GSE18864 , a cohort of 28 triple-negative breast cancers (24 with complete data) treated with neoadjuvant cisplatin. We tested the eight-gene TIL signature in these cohorts by extracting expression data for the eight genes and processing it as described above. In each cohort, patients were divided into two groups based on clinical outcome. For GSE21974, this was based on pathologic complete response or residual disease; for GSE19615, patients were divided based on the presence or absence of distant recurrence after 36 months of follow-up; and for GSE18864, patients were divided into good and poor response groups based on Miller-Payne scores of 3, 4, or 5 versus scores of 0, 1, or 2, respectively. For each response group, the average expression level of each gene was calculated and plotted as a box-and-whiskers plot; Mann-Whitney t tests were used to assess differences in overall TIL-signature expression.
TIL-enriched tumors are highly sensitive to neoadjuvant anthracycline-based therapy
Two-sided Fisher's exact tests for association of clinical variables and TIL-signature status with probability of pathological complete response (pCR) in the EORTC cohort
Total (n= 111)
FEC treated (n= 54)
TET treated (n= 57)
OR (95% CI)
OR (95% CI)
OR (95% CI)
Size T3 vs T1/T2
Grade 3 vs 1/2
Node 1 vs 0a
Age ≥50 vs <50
Her2 A vs Nb
TET vs FEC
TIL high vs low
Logistic regression model for prediction of pathologic complete response (pCR) in the EORTC cohort (n = 99)
OR (95% CI)
Size T3 vs T1/T2
Grade 3 vs 1/2
Node 1 vs 0a
Age ≥50 vs <50
Her2 A vs Nb
FEC vs TET
TIL high vs low
The eight genes of our TIL signature were selected because they are expressed primarily by T and B lymphocytes. As shown in Figure 1a, the genes clustered into separate B-cell (CD19 and MS4A1) and T-cell groups (CD3D, LCK, CD48, SELL, GZMB, PRF1). To determine the relative importance of each gene in the signature, we individually assessed their associations with outcome by comparing pCR rates in cases with expression levels in the upper quartile with those of cases at or below this cut-point (Figure 1c). Intriguingly, both B-cell genes failed to individually predict pCR (P > 0.1), whereas all T-cell genes reached statistical significance. GZMB, a marker of cytotoxic T lymphocytes (CTL), was the strongest single predictor of pCR (P = 0.0008), performing equivalently to the mean expression level of all eight genes.
We further assessed the association of the eight-gene TIL signature with response to anthracycline-based therapy in two small independent cohorts of ER-negative breast cancer with publically available gene-expression and clinical outcomes data (see Methods for details and Additional file 1: Supplementary Figure S1). One cohort consisted of both pre-chemo biopsies and post-chemo surgical specimens; in this group, the samples from complete clinical responders (both before and after chemotherapy) featured significantly higher levels of TIL-signature expression than did those with residual disease (P < 0.001). The second cohort consisted of patients treated with anthracycline-based therapy in the adjuvant setting. When patients were categorized based on the presence or absence of distant metastasis after 36 months of follow-up, TIL-signature expression was significantly elevated in disease-free cases (P < 0.001).
To determine the association of the eight-gene TIL signature with outcome in a distinct treatment setting, we performed the same analysis on a cohort of patients treated with neoadjuvant cisplatin therapy (Figure S1). In contrast with the anthracycline-based cohorts described thus far, TIL-signature expression was not elevated in patients with favorable responses to cisplatin therapy, suggesting that the eight-gene signature may be associated with clinical response in a regimen-specific manner.
Tumor-infiltrating T-cells are associated with disease-free survival in patients treated with adjuvant anthracycline-based therapy
Two-sided Fisher's exact tests for association of CD3 statusa and clinical parameters in the MBTB cohort
Age at diagnosis
CMF chemotherapy was administered to approximately one third of the patients in this cohort. Unlike anthracycline-based therapy, CMF effectively increased survival rates for both CD3-high (HR, 0.37; 95% CI, 0.18 to 0.76; P = 0.0064) and CD3-low patients (HR, 0.34; 95% CI, 0.17 to 0.67; P = 0.0017; Figure 2c). Because no significant clinical differences were found between anthracycline- and CMF-treated patients in this cohort, this implies that CD3 status is predictive of treatment response specifically for anthracycline-based regimens. Although CD3-high, anthracycline-treated patients had a higher 5-year DFS rate than did the CMF-treated patients (86.7% and 70.8%, respectively), this did not reach statistical significance.
Cytotoxic T-cells are the primary TIL subtype associated with anthracycline sensitivity
TIL predict anthracycline sensitivity in distinct molecular subtypes of ER-negative breast cancer
Her2 overexpression is frequently observed in ER-negative breast cancer and is a well-known marker of poor prognosis. While trastuzumab use has markedly improved management of Her2-positive breast tumors, standard chemotherapeutic approaches continue to play a role in the treatment of these cancers. Although our second cohort did not contain enough Her2-positive cases treated with anthracyclines to permit rigorous survival analysis, we assessed clinical responses in the Her2-positive subset of the neoadjuvant (EORTC) cohort. Of seven cases classified as TIL-high, six achieved a pCR (85.7%), compared with only eight of 23 TIL-low cases (34.8%; Fisher's exact test: P = 0.0309; OR, 11.25; 95% CI, 1.145 to 110.5).
Unlike ER-positive and Her2-positive breast tumors, which are candidates for endocrine and trastuzumab therapy, respectively, no targeted therapies for triple-negative breast tumors (TN; those lacking ER, PR, and Her2 overexpression) currently exist, making these a priority of translational research efforts. We therefore investigated whether our findings also applied to TN tumors. Among cases with documented Her2, ER, and PR status, TN tumors constituted 71% of the neoadjuvant (EORTC) cohort and 61% of the adjuvant (MBTB) cohort. In the neoadjuvant cohort, the predictive capacity of TIL status for anthracycline responsiveness in the TN subset was comparable to that of the entire dataset; the pCR rates of the TIL-high and -low groups were 73.9% and 30.9%, respectively (Fisher's exact test: P = 0.0009; OR, 6.333; 95% CI, 2.124 to 18.89). Neither CD19 nor MS4A1 alone correlated with pCR, whereas GZMB did so significantly (P = 0.0042). Similarly, in TN tumors of the MBTB adjuvant cohort, anthracycline-based therapy was associated with improved survival in the CD3-high group (log-rank HR, 0.25; 95% CI, 0.09 to 0.66; P = 0.0056), but not the CD3-low group (log-rank HR, 0.65; 95% CI, 0.24 to 1.79; P = 0.4054; data not shown). Consistent with the neoadjuvant cohort, CD8+ (P = 0.0390) and TIA-1+ (P = 0.0040) TIL were associated with improved outcomes after anthracycline-based therapy, but not CD4+ or CD20+ TIL. Therefore, TIL (particularly CTL) are associated with sensitivity to anthracycline-based therapy in the clinically challenging triple-negative subset of breast cancer.
We have demonstrated that a subgroup of ER-negative breast cancer defined by high lymphocyte gene expression has a remarkably high rate (74%) of pathologic complete response to neoadjuvant anthracycline-based therapy. We have further shown that the clinical benefit of adjuvant anthracycline-based therapy may be restricted to patients with high numbers of tumor-infiltrating T-cells, especially those of the CD8+ and TIA-1+ subsets. In contrast, TIL show no association with outcome after classic CMF therapy.
Our findings are consistent with several prior reports. In 2008, Hornychová et al.  observed high levels of intraepithelial CD3+ TIL in pretreatment biopsies of patients who achieved pCR after neoadjuvant therapy with doxorubicin and paclitaxel. Ladoire et al.  analyzed pre- and post-treatment samples from 56 patients treated with neoadjuvant anthracycline-based therapy and found that high ratios of CD8+ to Foxp3+ TIL (putative regulatory T-cells (Tregs)) in surgical specimens were correlated with a high rate of pCR. Consistent with our own findings, they also reported that, relative to pretreatment levels, TIA-1+ and GzmB+ TIL (CTL) were significantly increased in the posttreatment samples of patients who achieved pCR, but not in those who had residual disease. Finally, Denkert et al.  reported a significant relationship between TIL (identified by a combination of H&E assessment, IHC for CD3 and CD20, and expression analysis of several TIL genes by PCR) and pathologic response to neoadjuvant anthracycline/taxane therapy in a large group of 1,058 patients, approximately one fourth of whom were ER-negative. As such, based on our data and prior reports, TIL are consistently associated with improved outcomes after anthracycline-based therapy. Intriguingly, in 2001, Demaria et al.  observed in a small group of breast cancer patients (n = 25) treated with neoadjuvant paclitaxel therapy that levels of CD3+ TIL increased in the tumor bed after treatment, and that this correlated with therapeutic response. Although our data imply that TIL are not related to therapeutic responses in patients treated with cisplatin or CMF, the Demaria study nevertheless suggests that the correlation of TIL with treatment response may not be restricted entirely to anthracycline-based therapy.
It is important to note that TIL have been correlated with improved survival in several breast cancer studies to date, the association being observed primarily in high-grade, ER-negative, or Her2-positive lesions [4–10]. Given that many of the patients evaluated in these studies were not systemically treated, TIL may also be a feature of breast cancers with a naturally favorable clinical outcome, which could have implications for the interpretation of the results in our study. This may also be the case for a variety of other cancer types, including malignancies of the colon , ovary , lung , and bladder . In the systemically untreated patients of our long-term follow-up cohort, however, those with high levels of CD3+ or CD8+ TIL had only a modest, nonsignificant survival advantage over those with low TIL levels. It is therefore unlikely that the inherent association between TIL and improved prognosis was a significant underlying feature of our observations with regard to chemotherapy response.
The mechanism underlying the relationship between TIL and anthracycline sensitivity in humans is not currently understood. It is possible that TIL may not be causally related to therapeutic outcomes, and are instead biomarkers of an unknown feature that burdens tumors with anthracycline sensitivity. However, it is equally plausible that TIL directly mediate treatment responses. Indeed, this is supported in preclinical models, discussed here at length.
Four primary explanations for the relationship between anthracyclines and immunity have arisen based on observations in both the laboratory and the clinic [36, 37]. One possible scenario is that anthracyclines kill or suppress Tregs in the tumor microenvironment, thus relieving inhibition of the anti-tumor immune response. For example, in the aforementioned study by Ladoire et al. , Foxp3+ TIL were significantly reduced in the tumor bed after successful neoadjuvant anthracycline-based therapy, whereas their levels were unchanged in patients with residual disease. However, interpretation of this result is difficult because patients in the Ladoire study received a regimen (FEC) containing cyclophosphamide, a drug known to deplete Tregs selectively at low doses (as well as conventional T-cells at high doses) [38–41]. Similarly, many of the drug combinations used to treat patients in our study also included cyclophosphamide. Further studies designed to elucidate the specific effects of anthracyclines on intratumoral Treg populations will be required to determine the validity of this concept.
A second possible mechanism is derived from observations that transient lymphopenia (induced therapeutically via low-dose total-body irradiation or by myelosuppressive drugs such as cyclophosphamide) can enhance immunotherapy in the context of adoptive transfer or vaccination. Through a variety of mechanisms, therapeutic lymphopenia appears to trigger homeostatic processes that favor the proliferation and functionality of anti-tumor effector T-cells, while also ameliorating tumor-induced immunosuppression . However, considering the myelosuppressive properties of cyclophosphamide, this idea does not account for our observation that TIL are associated with outcome after anthracycline-based but not CMF chemotherapy.
A third hypothesis is that chemotherapy can sensitize tumor cells to T-cell-mediated cytotoxicity. For example, Yang and Haluska  found that treatment of human melanoma cell lines with either 5-fluorouracil or dacarbazine improved the efficacy of perforin/granzyme-mediated killing by antigen-specific CTL. Similar results have been reported with respect to paclitaxel, cisplatin, and doxorubicin . Cyclophosphamide has also been shown to sensitize mesothelioma cells to CTL killing, although through a TRAIL (TNF-related apoptosis-inducing ligand)-mediated effect . Finally, doxorubicin has been shown to induce Fas expression in breast cancer cell lines, leading to increased susceptibility to Fas ligand-mediated apoptosis . These studies collectively suggest an attractive model involving cytotoxic synergy between chemotherapy and the immune system. Nevertheless, as with the myelosuppression model discussed above, this mechanism would support a relationship between TIL and CMF therapy, which is not observed in our study.
The fourth hypothesis is that anthracyclines, unlike most chemotherapeutic drugs, induce immunogenic tumor cell death and thereby function indirectly as immunostimulatory agents . Casares et al.  treated murine CT26 colon carcinoma cells with doxorubicin and inoculated the dying cells into syngeneic hosts, which resulted in protection from concurrent or subsequent challenge with live CT26 cells. Notably, this effect was dependent on dendritic cells and CD8+ T-cells and was specific for the CT26 tumor line, indicating an immunologic mechanism. In contrast, cells stimulated to undergo apoptosis by using mitomycin-C, a non-anthracycline DNA cross-linking agent, did not elicit protective immunity. Obeid et al.  expanded on this work by demonstrating that anthracyclines trigger preapoptotic shuttling of calreticulin to the dying cell surface, where it acts as an opsonin to stimulate phagocytosis and processing by dendritic cells. Compared with a variety of apoptosis-inducing drugs with various mechanisms of action (including staurosporine, tunicamycin, mitomycin-C, and etoposide), anthracyclines induced greater calreticulin surface presentation on target cells, which correlated with the ability to elicit protective immunity. Indeed, blockade or knockdown of calreticulin in anthracycline-treated cells abolished their immunogenicity in mice . Finally, Apetoh et al.  demonstrated that doxorubicin additionally causes the extracellular release of HMGB1 (high-mobility-group box 1), an endogenous ligand of Toll-like receptor 4 (TLR4). HMGB1 stimulation of dendritic cells was necessary for processing and cross-presentation of antigens from dying tumor cells. The authors further showed that breast cancer patients bearing loss-of-function TLR4 alleles had significantly worse prognosis when treated with anthracycline-based chemotherapy . Taken together, mounting evidence indicates that anthracyclines have the rare ability to elicit specific immunity against target cells via stimulation of dendritic cells in a calreticulin and HMGB1/TLR-4-dependent fashion. This may explain the observation that anthracycline-based therapy appears most effective in breast cancer patients with high levels of endogenous tumor immunity, as well as the finding of Ladoire et al.  that CTL accumulate in the tumor beds of patients treated successfully with neoadjuvant anthracycline-based therapy. One unanswered question is whether anthracycline-based treatments can induce new tumor antigen-specific T-cells (de novo immunity) or, alternatively, augment the expansion or functional activity of preexisting tumor-reactive T-cells.
An important confounding feature of our study is the heterogeneous and multivalent chemotherapy used in the assessed cohorts. Although we have emphasized anthracyclines as the unifying components of the treatments we assessed, it is difficult to determine the extent to which anthracyclines are specifically responsible for our observations. For example, many of the patients in this study were treated with regimens involving taxanes, which may also have clinically relevant effects on host immunity . Furthermore, the potentially competing or cooperative effects on tumor immunity that agents such as anthracyclines, cyclophosphamide, and taxanes may possess when administered sequentially or simultaneously are not currently understood. A second caveat worthy of note pertains to the use of gene-expression arrays to infer the presence of TIL. The use of IHC to detect TIL is advantageous for several reasons, including the ability to directly quantify cells that express a given marker, as well as to accurately determine their specific localization within a tissue. This information cannot be acquired by using conventional methods of gene-expression analysis. Thus, although the data presented for the eight-gene TIL signature reflect a relative abundance of TIL, they lack the desirable single-cell resolution that is achievable through IHC.
A further limitation of this study is that the adjuvant cohort did not contain enough anthracycline-treated cases to permit rigorous multivariate analysis (although this was performed in the neoadjuvant cohort). This was largely because ER-negative tumors constitute only ~30% of breast cancer cases, and CMF was the dominant systemic regimen for breast cancer in Canada at the time that this cohort was treated. Our additional validation cohorts are also small for two primary reasons: first, because many publically available datasets are derived from fine-needle aspirates (as opposed to core or open-biopsy specimens), and so may not accurately reflect overall lymphocyte infiltration; and second, because most cohorts contain both ER-positive and ER-negative tumors, with relatively few cases of the latter. Our observations thus warrant validation in large, independent cohorts with long-term follow-up. Ultimately, prospective trials of therapeutic outcome prediction will be required to fully assess the clinical utility of the concepts presented herein.
Our findings indicate that TIL, which are easily detectable with standard IHC, may prove useful for identifying patients who are likely to benefit from anthracycline-based therapy. Moreover, our findings suggest that patients with low TIL levels may benefit most from alternative regimens such as CMF. In addition, our study provides support for the idea of combining anthracycline-based chemotherapy with immunotherapeutic strategies such as cancer vaccines or adoptive T-cell therapy to fully engage the host immune response against breast cancer.
TIL that express cytotoxic markers are strongly associated with favorable outcome after anthracycline-based treatment of ER-negative breast cancer, including the clinically challenging triple-negative subset. Further studies to develop TIL-based predictive and therapeutic strategies for breast cancer are warranted.
cytotoxic T lymphocyte
estrogen receptor alpha
human epidermal growth factor receptor 2
pathologic complete response
T-cell-restricted intracellular antigen 1
regulatory T cell.
We thank Rebecca Barnes, Sindy Babinszky (BC Cancer Agency), and Michelle Parisien (MBTB) for technical assistance and members of the Deeley Research Centre for helpful discussion. This study was supported by funding from the Canadian Institutes of Health Research (CIHR, grant MOP-64349) and the BC Cancer Foundation. The Manitoba Breast Tumor Bank, a member of the Canadian Tumor Repository Network, is supported by CIHR grant PRG80155. NRW is supported by a US DOD Breast Cancer Research Program predoctoral traineeship award (W81XWH-08-1-0781).
- Gianni L, Norton L, Wolmark N, Suter TM, Bonadonna G, Hortobagyi GN: Role of anthracyclines in the treatment of early breast cancer. J Clin Oncol. 2009, 27: 4798-4808. 10.1200/JCO.2008.21.4791.PubMedView ArticleGoogle Scholar
- Munro AF, Cameron DA, Bartlett JM: Targeting anthracyclines in early breast cancer: new candidate predictive biomarkers emerge. Oncogene. 2010, 29: 5231-5240. 10.1038/onc.2010.286.PubMedView ArticleGoogle Scholar
- Moretti E, Oakman C, Di Leo A: Predicting anthracycline benefit: have we made any progress?. Curr Opin Oncol. 2009, 21: 507-515. 10.1097/CCO.0b013e328331a501.PubMedView ArticleGoogle Scholar
- Alexe G, Dalgin GS, Scanfeld D, Tamayo P, Mesirov JP, DeLisi C, Harris L, Barnard N, Martel M, Levine AJ, Ganesan S, Bhanot G: High expression of lymphocyte-associated genes in node-negative HER2+ breast cancers correlates with lower recurrence rates. Cancer Res. 2007, 67: 10669-10676. 10.1158/0008-5472.CAN-07-0539.PubMedView ArticleGoogle Scholar
- Arnould L, Gelly M, Penault-Llorca F, Benoit L, Bonnetain F, Migeon C, Cabaret V, Fermeaux V, Bertheau P, Garnier J, Jeannin JF, Coudert B: Trastuzumab-based treatment of HER2-positive breast cancer: an antibody-dependent cellular cytotoxicity mechanism?. Br J Cancer. 2006, 94: 259-267. 10.1038/sj.bjc.6602930.PubMedPubMed CentralView ArticleGoogle Scholar
- Calabrò A, Beissbarth T, Kuner R, Stojanov M, Benner A, Asslaber M, Ploner F, Zatloukal K, Samonigg H, Poustka A, Sültmann H: Effects of infiltrating lymphocytes and estrogen receptor on gene expression and prognosis in breast cancer. Breast Cancer Res Treat. 2009, 116: 69-77. 10.1007/s10549-008-0105-3.PubMedView ArticleGoogle Scholar
- Desmedt C, Haibe-Kains B, Wirapati P, Buyse M, Larsimont D, Bontempi G, Delorenzi M, Piccart M, Sotiriou C: Biological processes associated with breast cancer clinical outcome depend on the molecular subtypes. Clin Cancer Res. 2008, 14: 5158-5165. 10.1158/1078-0432.CCR-07-4756.PubMedView ArticleGoogle Scholar
- Schmidt M, Bohm D, von Torne C, Steiner E, Puhl A, Pilch H, Lehr HA, Hengstler JG, Kolbl H, Gehrmann M: The humoral immune system has a key prognostic impact in node-negative breast cancer. Cancer Res. 2008, 68: 5405-5413. 10.1158/0008-5472.CAN-07-5206.PubMedView ArticleGoogle Scholar
- Teschendorff AE, Miremadi A, Pinder SE, Ellis IO, Caldas C: An immune response gene expression module identifies a good prognosis subtype in estrogen receptor negative breast cancer. Genome Biol. 2007, 8: R157-10.1186/gb-2007-8-8-r157.PubMedPubMed CentralView ArticleGoogle Scholar
- Rody A, Holtrich U, Pusztai L, Liedtke C, Gaetje R, Ruckhaeberle E, Solbach C, Hanker L, Ahr A, Metzler D, Engels K, Karn T, Kaufmann M: T-cell metagene predicts a favorable prognosis in estrogen receptor-negative and HER2-positive breast cancers. Breast Cancer Res. 2009, 11: R15-10.1186/bcr2234.PubMedPubMed CentralView ArticleGoogle Scholar
- Denkert C, Loibl S, Noske A, Roller M, Muller BM, Komor M, Budczies J, Darb-Esfahani S, Kronenwett R, Hanusch C, von Törne C, Weichert W, Engels K, Solbach C, Schrader I, Dietel M, von Minckwitz G: Tumor-associated lymphocytes as an independent predictor of response to neoadjuvant chemotherapy in breast cancer. J Clin Oncol. 2010, 28: 105-113. 10.1200/JCO.2009.23.7370.PubMedView ArticleGoogle Scholar
- Colleoni M, Bagnardi V, Rotmensz N, Gelber RD, Viale G, Pruneri G, Veronesi P, Torrisi R, Cardillo A, Montagna E, Campagnoli E, Luini A, Intra M, Galimberti V, Scarano E, Peruzzotti G, Goldhirsch A: Increasing steroid hormone receptors expression defines breast cancer subtypes nonresponsive to preoperative chemotherapy. Breast Cancer Res Treat. 2009, 116: 359-369. 10.1007/s10549-008-0223-y.PubMedView ArticleGoogle Scholar
- Colleoni M, Viale G, Zahrieh D, Pruneri G, Gentilini O, Veronesi P, Gelber RD, Curigliano G, Torrisi R, Luini A, Intra M, Galimberti V, Renne G, Nolè F, Peruzzotti G, Goldhirsch A: Chemotherapy is more effective in patients with breast cancer not expressing steroid hormone receptors: a study of preoperative treatment. Clin Cancer Res. 2004, 10: 6622-6628. 10.1158/1078-0432.CCR-04-0380.PubMedView ArticleGoogle Scholar
- Guarneri V, Broglio K, Kau SW, Cristofanilli M, Buzdar AU, Valero V, Buchholz T, Meric F, Middleton L, Hortobagyi GN, Gonzalez-Angulo AM: Prognostic value of pathologic complete response after primary chemotherapy in relation to hormone receptor status and other factors. J Clin Oncol. 2006, 24: 1037-1044. 10.1200/JCO.2005.02.6914.PubMedView ArticleGoogle Scholar
- Ring AE, Smith IE, Ashley S, Fulford LG, Lakhani SR: Oestrogen receptor status, pathological complete response and prognosis in patients receiving neoadjuvant chemotherapy for early breast cancer. Br J Cancer. 2004, 91: 2012-2017. 10.1038/sj.bjc.6602235.PubMedPubMed CentralView ArticleGoogle Scholar
- Bates GJ, Fox SB, Han C, Leek RD, Garcia JF, Harris AL, Banham AH: Quantification of regulatory T cells enables the identification of high-risk breast cancer patients and those at risk of late relapse. J Clin Oncol. 2006, 24: 5373-5380. 10.1200/JCO.2006.05.9584.PubMedView ArticleGoogle Scholar
- Chavey C, Bibeau F, Gourgou-Bourgade S, Burlinchon S, Boissiere F, Laune D, Roques S, Lazennec G: Oestrogen receptor negative breast cancers exhibit high cytokine content. Breast Cancer Res. 2007, 9: R15-10.1186/bcr1648.PubMedPubMed CentralView ArticleGoogle Scholar
- Reis-Filho JS, Tutt AN: Triple negative tumours: a critical review. Histopathology. 2008, 52: 108-118.PubMedView ArticleGoogle Scholar
- Bonnefoi H, Potti A, Delorenzi M, Mauriac L, Campone M, Tubiana-Hulin M, Petit T, Rouanet P, Jassem J, Blot E, Becette V, Farmer P, André S, Acharya CR, Mukherjee S, Cameron D, Bergh J, Nevins JR, Iggo RD: Validation of gene signatures that predict the response of breast cancer to neoadjuvant chemotherapy: a substudy of the EORTC 10994/BIG 00-01 clinical trial. Lancet Oncol. 2007, 8: 1071-1078. 10.1016/S1470-2045(07)70345-5.PubMedView ArticleGoogle Scholar
- Li Y, Zou L, Li Q, Haibe-Kains B, Tian R, Li Y, Desmedt C, Sotiriou C, Szallasi Z, Iglehart JD, Richardson AL, Wang ZC: Amplification of LAPTM4B and YWHAZ contributes to chemotherapy resistance and recurrence of breast cancer. Nat Med. 2010, 16: 214-218. 10.1038/nm.2090.PubMedPubMed CentralView ArticleGoogle Scholar
- Silver DP, Richardson AL, Eklund AC, Wang ZC, Szallasi Z, Li Q, Juul N, Leong CO, Calogrias D, Buraimoh A, Fatima A, Gelman RS, Ryan PD, Tung NM, De Nicolo A, Ganesan S, Miron A, Colin C, Sgroi DC, Ellisen LW, Winer EP, Garber JE: Efficacy of neoadjuvant cisplatin in triple-negative breast cancer. J Clin Oncol. 2010, 28: 1145-1153. 10.1200/JCO.2009.22.4725.PubMedPubMed CentralView ArticleGoogle Scholar
- Alvaro-Naranjo T, Lejeune M, Salvado-Usach MT, Bosch-Princep R, Reverter-Branchat G, Jaen-Martinez J, Pons-Ferre LE: Tumor-infiltrating cells as a prognostic factor in Hodgkin's lymphoma: a quantitative tissue microarray study in a large retrospective cohort of 267 patients. Leuk Lymphoma. 2005, 46: 1581-1591. 10.1080/10428190500220654.PubMedView ArticleGoogle Scholar
- Anderson P, Nagler-Anderson C, O'Brien C, Levine H, Watkins S, Slayter HS, Blue ML, Schlossman SF: A monoclonal antibody reactive with a 15-kDa cytoplasmic granule-associated protein defines a subpopulation of CD8+ T lymphocytes. J Immunol. 1990, 144: 574-582.PubMedGoogle Scholar
- Dalmazzo LF, Jacomo RH, Marinato AF, Figueiredo-Pontes LL, Cunha RL, Garcia AB, Rego EM, Falcao RP: The presence of CD56/CD16 in T-cell acute lymphoblastic leukaemia correlates with the expression of cytotoxic molecules and is associated with worse response to treatment. Br J Haematol. 2009, 144: 223-229. 10.1111/j.1365-2141.2008.07457.x.PubMedView ArticleGoogle Scholar
- Felgar RE, Salhany KE, Macon WR, Pietra GG, Kinney MC: The expression of TIA-1+ cytolytic-type granules and other cytolytic lymphocyte-associated markers in CD30+ anaplastic large cell lymphomas (ALCL): correlation with morphology, immunophenotype, ultrastructure, and clinical features. Hum Pathol. 1999, 30: 228-236. 10.1016/S0046-8177(99)90281-2.PubMedView ArticleGoogle Scholar
- Milne K, Barnes RO, Girardin A, Mawer MA, Nesslinger NJ, Ng A, Nielsen JS, Sahota R, Tran E, Webb JR, Wong MQ, Wick DA, Wray A, McMurtrie E, Köbel M, Kalloger SE, Gilks CB, Watson PH, Nelson BH: Tumor-infiltrating T cells correlate with NY-ESO-1-specific autoantibodies in ovarian cancer. PLoS One. 2008, 3: e3409-10.1371/journal.pone.0003409.PubMedPubMed CentralView ArticleGoogle Scholar
- Tian Q, Streuli M, Saito H, Schlossman SF, Anderson P: A polyadenylate binding protein localized to the granules of cytolytic lymphocytes induces DNA fragmentation in target cells. Cell. 1991, 67: 629-639. 10.1016/0092-8674(91)90536-8.PubMedView ArticleGoogle Scholar
- Vermeer MH, van Doorn R, Dukers D, Bekkenk MW, Meijer CJ, Willemze R: CD8+ T cells in cutaneous T-cell lymphoma: expression of cytotoxic proteins, Fas ligand, and killing inhibitory receptors and their relationship with clinical behavior. J Clin Oncol. 2001, 19: 4322-4329.PubMedGoogle Scholar
- Hornychova H, Melichar B, Tomsova M, Mergancova J, Urminska H, Ryska A: Tumor-infiltrating lymphocytes predict response to neoadjuvant chemotherapy in patients with breast carcinoma. Cancer Invest. 2008, 26: 1024-1031. 10.1080/07357900802098165.PubMedView ArticleGoogle Scholar
- Ladoire S, Arnould L, Apetoh L, Coudert B, Martin F, Chauffert B, Fumoleau P, Ghiringhelli F: Pathologic complete response to neoadjuvant chemotherapy of breast carcinoma is associated with the disappearance of tumor-infiltrating foxp3+ regulatory T cells. Clin Cancer Res. 2008, 14: 2413-2420. 10.1158/1078-0432.CCR-07-4491.PubMedView ArticleGoogle Scholar
- Demaria S, Volm MD, Shapiro RL, Yee HT, Oratz R, Formenti SC, Muggia F, Symmans WF: Development of tumor-infiltrating lymphocytes in breast cancer after neoadjuvant paclitaxel chemotherapy. Clin Cancer Res. 2001, 7: 3025-3030.PubMedGoogle Scholar
- Galon J, Costes A, Sanchez-Cabo F, Kirilovsky A, Mlecnik B, Lagorce-Pagès C, Tosolini M, Camus M, Berger A, Wind P, Zinzindohoué F, Bruneval P, Cugnenc PH, Trajanoski Z, Fridman WH, Pagès F: Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Science. 2006, 313: 1960-1964. 10.1126/science.1129139.PubMedView ArticleGoogle Scholar
- Nelson BH: The impact of T-cell immunity on ovarian cancer outcomes. Immunol Rev. 2008, 222: 101-116. 10.1111/j.1600-065X.2008.00614.x.PubMedView ArticleGoogle Scholar
- Al-Shibli KI, Donnem T, Al-Saad S, Persson M, Bremnes RM, Busund LT: Prognostic effect of epithelial and stromal lymphocyte infiltration in non-small cell lung cancer. Clin Cancer Res. 2008, 14: 5220-5227. 10.1158/1078-0432.CCR-08-0133.PubMedView ArticleGoogle Scholar
- Sharma P, Shen Y, Wen S, Yamada S, Jungbluth AA, Gnjatic S, Bajorin DF, Reuter VE, Herr H, Old LJ, Sato E: CD8 tumor-infiltrating lymphocytes are predictive of survival in muscle-invasive urothelial carcinoma. Proc Natl Acad Sci USA. 2007, 104: 3967-3972. 10.1073/pnas.0611618104.PubMedPubMed CentralView ArticleGoogle Scholar
- Andre F, Berrada N, Desmedt C: Implication of tumor microenvironment in the resistance to chemotherapy in breast cancer patients. Curr Opin Oncol. 2010, 22: 547-551. 10.1097/CCO.0b013e32833fb384.PubMedView ArticleGoogle Scholar
- Emens LA: Chemoimmunotherapy. Cancer J. 2010, 16: 295-303. 10.1097/PPO.0b013e3181eb5066.PubMedPubMed CentralView ArticleGoogle Scholar
- Ercolini AM, Ladle BH, Manning EA, Pfannenstiel LW, Armstrong TD, Machiels JP, Bieler JG, Emens LA, Reilly RT, Jaffee EM: Recruitment of latent pools of high-avidity CD8(+) T cells to the antitumor immune response. J Exp Med. 2005, 201: 1591-1602. 10.1084/jem.20042167.PubMedPubMed CentralView ArticleGoogle Scholar
- Ghiringhelli F, Larmonier N, Schmitt E, Parcellier A, Cathelin D, Garrido C, Chauffert B, Solary E, Bonnotte B, Martin F: CD4+CD25+ regulatory T cells suppress tumor immunity but are sensitive to cyclophosphamide which allows immunotherapy of established tumors to be curative. Eur J Immunol. 2004, 34: 336-344. 10.1002/eji.200324181.PubMedView ArticleGoogle Scholar
- Lutsiak ME, Semnani RT, De Pascalis R, Kashmiri SV, Schlom J, Sabzevari H: Inhibition of CD4(+)25+ T regulatory cell function implicated in enhanced immune response by low-dose cyclophosphamide. Blood. 2005, 105: 2862-2868. 10.1182/blood-2004-06-2410.PubMedView ArticleGoogle Scholar
- Taieb J, Chaput N, Schartz N, Roux S, Novault S, Menard C, Ghiringhelli F, Terme M, Carpentier AF, Darrasse-Jeze G, Lemonnier F, Zitvogel L: Chemoimmunotherapy of tumors: cyclophosphamide synergizes with exosome based vaccines. J Immunol. 2006, 176: 2722-2729.PubMedView ArticleGoogle Scholar
- Zitvogel L, Apetoh L, Ghiringhelli F, Kroemer G: Immunological aspects of cancer chemotherapy. Nat Rev Immunol. 2008, 8: 59-73. 10.1038/nri2216.PubMedView ArticleGoogle Scholar
- Yang S, Haluska FG: Treatment of melanoma with 5-fluorouracil or dacarbazine in vitro sensitizes cells to antigen-specific CTL lysis through perforin/granzyme- and Fas-mediated pathways. J Immunol. 2004, 172: 4599-4608.PubMedView ArticleGoogle Scholar
- Ramakrishnan R, Assudani D, Nagaraj S, Hunter T, Cho HI, Antonia S, Altiok S, Celis E, Gabrilovich DI: Chemotherapy enhances tumor cell susceptibility to CTL-mediated killing during cancer immunotherapy in mice. J Clin Invest. 2010, 120: 1111-1124. 10.1172/JCI40269.PubMedPubMed CentralView ArticleGoogle Scholar
- van der Most RG, Currie AJ, Cleaver AL, Salmons J, Nowak AK, Mahendran S, Larma I, Prosser A, Robinson BW, Smyth MJ, Scalzo AA, Degli-Esposti MA, Lake RA: Cyclophosphamide chemotherapy sensitizes tumor cells to TRAIL-dependent CD8 T cell-mediated immune attack resulting in suppression of tumor growth. PLoS One. 2009, 4: e6982-10.1371/journal.pone.0006982.PubMedPubMed CentralView ArticleGoogle Scholar
- Crescenzi E, Pacifico F, Lavorgna A, De Palma R, D'Aiuto E, Palumbo G, Formisano S, Leonardi A: NF-kappaB-dependent cytokine secretion controls Fas expression on chemotherapy-induced premature senescent tumor cells. Oncogene. 2011, 30: 2707-2717. 10.1038/onc.2011.1.PubMedView ArticleGoogle Scholar
- Casares N, Pequignot MO, Tesniere A, Ghiringhelli F, Roux S, Chaput N, Schmitt E, Hamai A, Hervas-Stubbs S, Obeid M, Coutant F, Métivier D, Pichard E, Aucouturier P, Pierron G, Garrido C, Zitvogel L, Kroemer G: Caspase-dependent immunogenicity of doxorubicin-induced tumor cell death. J Exp Med. 2005, 202: 1691-1701. 10.1084/jem.20050915.PubMedPubMed CentralView ArticleGoogle Scholar
- Obeid M, Tesniere A, Ghiringhelli F, Fimia GM, Apetoh L, Perfettini JL, Castedo M, Mignot G, Panaretakis T, Casares N, Métivier D, Larochette N, van Endert P, Ciccosanti F, Piacentini M, Zitvogel L, Kroemer G: Calreticulin exposure dictates the immunogenicity of cancer cell death. Nat Med. 2007, 13: 54-61. 10.1038/nm1523.PubMedView ArticleGoogle Scholar
- Apetoh L, Ghiringhelli F, Tesniere A, Obeid M, Ortiz C, Criollo A, Mignot G, Maiuri MC, Ullrich E, Saulnier P, Yang H, Amigorena S, Ryffel B, Barrat FJ, Saftig P, Levi F, Lidereau R, Nogues C, Mira JP, Chompret A, Joulin V, Clavel-Chapelon F, Bourhis J, André F, Delaloge S, Tursz T, Kroemer G, Zitvogel L: Toll-like receptor 4-dependent contribution of the immune system to anticancer chemotherapy and radiotherapy. Nat Med. 2007, 13: 1050-1059. 10.1038/nm1622.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.