Programmed death 1 ligand 1 (PD-L1) is an immune regulatory molecule that limits antitumor immune activity. Targeting of PD-L1 and other immune checkpoint proteins has shown therapeutic activity in various tumor types. The expression of PD-L1 and its correlation with response to neoadjuvant chemotherapy in breast cancer has not been studied extensively. Our goal was to assess PD-L1 expression in a cohort of breast cancer patients treated with neoadjuvant chemotherapy. Pretreatment biopsies from 105 patients with breast cancer from Yale New Haven Hospital that subsequently received neoadjuvant chemotherapy were assessed for PD-L1 protein expression by automated quantitative analysis with a rabbit monoclonal antibody (E1L3N) to the cytoplasmic domain of PD-L1. In addition, tumor-infiltrating lymphocytes (TIL) were assessed on hematoxylin and eosin slides. PD-L1 expression was observed in 30% of patients, and it was positively associated with hormone-receptor–negative and triple-negative status and high levels of TILs. Both TILs and PD-L1 measured in the epithelium or stroma predicted pathologic complete response (pCR) to neoadjuvant chemotherapy in univariate and multivariate analyses. However, because they are strongly associated, TILs and PD-L1 cannot both be included in a significant multivariate model. PD-L1 expression is prevalent in breast cancer, particularly hormone-receptor–negative and triple-negative patients, indicating a subset of patients that may benefit from immune therapy. Furthermore, PD-L1 and TILs correlate with pCR, and high PD-L1 predicts pCR in multivariate analysis. Cancer Immunol Res; 3(4); 326–32. ©2014 AACR.

Neoadjuvant chemotherapy is increasingly used in the management of stage II to III breast cancer, and pathologic complete response (pCR) is observed in 5% to 15% of estrogen-receptor (ER)–positive and 30% to 50% of triple-negative and HER2-positive patients with third-generation combination chemotherapy regimens (1, 2). High-grade and high Ki-67 expression also correlates with pCR, particularly among ER-positive cancers (3). Recent reports have found that the presence of tumor-infiltrating lymphocytes (TIL) predicts response to neoadjuvant chemotherapy (4–13). The presence of TILs may indicate immune-mediated host defense against the tumor, and TILs may contribute to and augment chemotherapy-induced cell death. The recent positive results with immune checkpoint inhibitors in melanoma and lung cancer have stimulated new interest in TILs and their relationship to tumor immunity and chemotherapy response (14, 15).

One key immune modulatory pathway is mediated by the PD-1–PD-L1 axis (programmed cell death 1 and its ligand). PD-L1 is a transmembrane protein of the B7 family of immune molecules that plays an integral role in limiting the cytotoxic immune response via interaction with programmed death-1 (PD-1) receptor (16). PD-L1 expression has been noted in a variety of cancers and reported in a variety of solid tumor types, including lung, melanoma, ovarian, colon, and breast (16–19). Its expression in tumor cells or presence in the tumor microenvironment has been correlated to the presence of TILs. Results from various preclinical studies using cell line and mouse models support the idea that inhibition of the interaction between PD-L1 and PD-1 in the tumor microenvironment may enhance antitumor immunity and promote tumor regression (16, 20–22). Various agents targeting PD-1 or PD-L1 are currently in clinical trials for a variety of solid tumor types and have demonstrated robust response rates, notably in metastatic melanoma, renal cell carcinoma, and non–small cell lung cancer (15, 23–28).

Our goal was to investigate the correlation of PD-L1, also known as B7-H1 and CD274, with TILs and pCR following neoadjuvant chemotherapy in breast cancer. We assess PD-L1 expression objectively by quantitative immunofluorescence on samples from a cohort of patients with breast cancer that received neoadjuvant chemotherapy. We describe localization and distribution of PD-L1 expression and relate this to TILs, clinical characteristics of the cancer, and response to neoadjuvant chemotherapy.

Patient cohort

The cohort used in this study consists of 94 presurgical biopsies from patients diagnosed with breast cancer between 2002 and 2010, who subsequently received neoadjuvant chemotherapy. Specimens were collected from the archives of the Department of Pathology at Yale University. A majority of patients (76.6%) received doxorubicin-based neoadjuvant chemotherapy. A more detailed characterization of the Yale Neoadjuvant Cohort has been published previously and is shown in Table 1 (3). pCR was defined as the absence of invasive carcinoma in the breast and sampled lymph nodes ypT0 ypN0.

Table 1.

Yale neoadjuvant cohort characteristics

Characteristicn (%)
Age (y) 
 <50 55 (58.5) 
 ≥50 39 (41.5) 
 Unknown 0 (0) 
Nodal status 
 Positive 48 (51.1) 
 Negative 32 (34.0) 
 Unknown 14 (14.9) 
Tumor size (cm) 
 <2 10 (10.6) 
 >2 83 (88.3) 
 Unknown 1 (1.1) 
Nuclear grade 
 1–2 48 (51.1) 
 3 43 (45.7) 
 Unknown 3 (3.2) 
Estrogen receptor 
 Negative 35 (37.2) 
 Positive 57 (60.7) 
 Unknown 2 (2.1) 
Progesterone receptor 
 Negative 42 (44.7) 
 Positive 50 (53.2) 
 Unknown 2 (2.1) 
HER2 
 Negative 68 (72.3) 
 Positive 25 (26.6) 
 Unknown 1 (1.1) 
Triple-negative 
 Yes 23 (24.4) 
 No 70 (74.5) 
 Unknown 1 (1.1) 
LPBC 
 Yes 8 (8.5) 
 No 84 (89.4) 
 Unknown 2 (2.1) 
Treatment 
 Adriamycin-based 72 (76.6) 
 Carboplatin + taxane + trastuzumab 10 (10.6) 
 Carboplatin + (nab)-paclitaxel + bevacizumab 6 (6.4) 
 Carboplatin + (nab)-paclitaxel + trastuzumab 3 (3.2) 
 Carboplatin + paclitaxel + etoposide 1 (1.1) 
 Cyclophosphamide + docetaxel 2 (2.1) 
Characteristicn (%)
Age (y) 
 <50 55 (58.5) 
 ≥50 39 (41.5) 
 Unknown 0 (0) 
Nodal status 
 Positive 48 (51.1) 
 Negative 32 (34.0) 
 Unknown 14 (14.9) 
Tumor size (cm) 
 <2 10 (10.6) 
 >2 83 (88.3) 
 Unknown 1 (1.1) 
Nuclear grade 
 1–2 48 (51.1) 
 3 43 (45.7) 
 Unknown 3 (3.2) 
Estrogen receptor 
 Negative 35 (37.2) 
 Positive 57 (60.7) 
 Unknown 2 (2.1) 
Progesterone receptor 
 Negative 42 (44.7) 
 Positive 50 (53.2) 
 Unknown 2 (2.1) 
HER2 
 Negative 68 (72.3) 
 Positive 25 (26.6) 
 Unknown 1 (1.1) 
Triple-negative 
 Yes 23 (24.4) 
 No 70 (74.5) 
 Unknown 1 (1.1) 
LPBC 
 Yes 8 (8.5) 
 No 84 (89.4) 
 Unknown 2 (2.1) 
Treatment 
 Adriamycin-based 72 (76.6) 
 Carboplatin + taxane + trastuzumab 10 (10.6) 
 Carboplatin + (nab)-paclitaxel + bevacizumab 6 (6.4) 
 Carboplatin + (nab)-paclitaxel + trastuzumab 3 (3.2) 
 Carboplatin + paclitaxel + etoposide 1 (1.1) 
 Cyclophosphamide + docetaxel 2 (2.1) 

Abbreviation: LBPC, lipocyte-predominant breast cancer.

Evaluation of TILs

Histopathologic analysis of TILs was performed on hematoxylin and eosin (H&E)–stained sections from the core biopsies of 94 patients from the cohort described above. Analysis was conducted by two pathologists (V. Bossuyt and C. Nixon), who were blinded to the clinical parameters and response. TILs were quantified as a percentage estimate of the stromal area adjacent to the tumor that contained lymphocytic infiltrate, as described in the literature (6). Percentages were reported in discrete increments of 10%, with 0% indicating a minimal infiltrate and 100% indicating the stroma almost exclusively consisted of TILs. Sections with 50% or greater TIL infiltrate were denoted as lymphocyte predominant breast cancer (LPBC) as discussed in the literature (4).

Quantitative immunofluorescence

Whole-tissue sections were baked overnight at 60°C and then soaked in xylene twice for 20 minutes each. Slides were rehydrated in two 1-minute washes in 100% ethanol followed by one wash in 70% ethanol and finally rinsed in streaming tap water for 5 minutes. Antigen retrieval was performed in sodium citrate buffer, pH 6, in the PT module from LabVision. Endogenous peroxidases were blocked by 30-minute incubation in 2.5% hydrogen peroxide in methanol. Subsequent steps were carried out on the LabVision 720 Autostainer (Thermo-Scientific). Nonspecific antigens were blocked by a 30-minute incubation in 0.3% BSA in TBST. Primary PD-L1 (E1L3N) rabbit monoclonal antibody (Cell Signaling Technology; clone E1L3N; see Supplementary Fig. S2 for antibody validation) was prepared to a working concentration of 3.5 μg/mL combined with 1:100 pan-cytokeratin antibody (Dako; Cat#Z062201-2) in 0.3% BSA in TBST and transferred to 4°C overnight. Primary antibodies were followed by incubation with Alexa 546–conjugated goat anti-rabbit secondary antibody (Life Technologies; Cat#A-11010) diluted 1:100 in mouse EnVision reagent (Dako; Cat#K400111-2) for 1 hour. Signal was amplified with Cy5-Tyramide (Perkin Elmer; Cat#SAT705A001EA) for 10 minutes, and then slides were mounted with ProlongGold + DAPI (Life Technologies; Cat#P36931).

Immunofluorescence was quantified using automated quantitative analysis (AQUA) on all regions of tissue on each slide. Briefly, fluorescent images of DAPI, Cy3 (Alexa 546-cytokeratin), and Cy5 (PD-L1) for each field of view were collected. The number of fields of view assessed per case ranged from 5 to 93. The average number of fields of view was 32. Image analysis was carried out using the AQUAnalysis software (Genoptix), which is generated for each compartment by dividing the sum of target pixel intensities by the area of the compartment in which the target is measured (29, 30). The stromal compartment was created by subtracting the epithelial tumor mask from a DAPI mask.

Statistical analysis

T tests were used to determine the correlations between continuous quantitative scores of PD-L1 expression and clinicopathologic factors as well as pCR. The χ2 tests were used to determine correlation of binary PD-L1 expression with pCR. Logistic regression was used for univariate and multivariate analyses. All statistical tests mentioned above were carried out using StatView (SAS Institute Inc.). Joinpoint software was used to dichotomize continuous PD-L1 AQUA scores. Briefly, average quantitative scores and the SD for each patient with greater than 4 fields of view of tissue were imported into Joinpoint software, which identifies trends in the population distribution, enabling an objective method of splitting a population in two (31).

Of the 94 cases collected for the neoadjuvant cohort, 14.9% (14/94) and 19.1% (18/94) were eliminated from consideration for PD-L1 staining in the epithelium and stroma, respectively, due to insufficient measureable tissue. Our criteria for the minimum amount of tissue for evaluation are 4 fields of view per tissue section with a minimum of 3% area within the field of view of the epithelial or stromal compartment. The number of fields of view analyzed per biopsy ranged from 5 to 93, with an average of 34. We could assess epithelial PD-L1 expression in 80 cases and stromal PD-L1 expression in 76 cases. Examples of epithelial and stromal PD-L1 expression can be seen in Supplementary Fig. S1C–S1F. Heat maps of AQUA scores generated on one whole-tissue section in the epithelial and stromal compartments are shown in Supplementary Fig. S1A and S1B, demonstrating a higher level of PD-L1 expression in the epithelium for the example given. The AQUA scores were reflective in each case of the predominant localization of PD-L1 expression, and Fig. 1A shows an example of the distribution of AQUA scores within a tissue section of epithelial-predominant expression, whereas the distribution of a case with predominantly stromal PD-L1 expression is illustrated in Fig. 1B. The distribution of PD-L1 expression in the epithelial and stromal compartments is similar (Fig. 1C and D), although epithelial expression showed higher signal, reflected by higher AQUA scores, than stromal expression (Fig. 1C and D). Although the distributions in Fig. 1C and D represent averages of scores from a number of fields of view from each patient, we note the heterogeneity of PD-L1 expression within each biopsy. Figure 1A and B illustrate that the level of PD-L1 expression can vary up to 4× in different areas of the same biopsy. When measured as a continuous quantitative score, high PD-L1 expression in epithelial cells or stroma is significantly associated with hormone-receptor–negative and triple-negative breast cancers. Using a threshold adopted from the literature as described in Materials and Methods, 8.5% (8/94) of cases were LPBC. We also find a significant positive association of PD-L1 with LPBC (Table 2).

Figure 1.

Epithelial and stromal PD-L1 expression correlates with pCR. A and B, distributions of epithelial and stromal PD-L1 expression with quantitative scores for each field of view in two representative cases. A frequency distribution of the epithelial (C) and stromal (D) PD-L1 quantitative scores. Box plots with continuous PD-L1 scores in epithelium (E) and stroma (F) on the y axis and pCR or no pCR on the x axis.

Figure 1.

Epithelial and stromal PD-L1 expression correlates with pCR. A and B, distributions of epithelial and stromal PD-L1 expression with quantitative scores for each field of view in two representative cases. A frequency distribution of the epithelial (C) and stromal (D) PD-L1 quantitative scores. Box plots with continuous PD-L1 scores in epithelium (E) and stroma (F) on the y axis and pCR or no pCR on the x axis.

Close modal
Table 2.

Epithelial and stromal PD-L1 expression association with Yale neoadjuvant cohort characteristics

P value
CharacteristicTumor maskStroma
Age (y) 
 <50 0.3036 0.4228 
 ≥50   
Nodal status 
 Positive 0.3062 0.3811 
 Negative   
Tumor size (cm) 
 <2 0.4708 0.8819 
 >2   
Estrogen receptor 
 Negative 0.0378a 0.0219a 
 Positive   
Progesterone receptor 
 Negative 0.0099a 0.0009a 
 Positive   
HER2 
 Negative 0.7381 0.7826 
 Positive   
Triple-negative 
 Yes 0.0026b 0.0021b 
 No   
LPBC 
 Yes <0.0001b <0.0001b 
 No   
P value
CharacteristicTumor maskStroma
Age (y) 
 <50 0.3036 0.4228 
 ≥50   
Nodal status 
 Positive 0.3062 0.3811 
 Negative   
Tumor size (cm) 
 <2 0.4708 0.8819 
 >2   
Estrogen receptor 
 Negative 0.0378a 0.0219a 
 Positive   
Progesterone receptor 
 Negative 0.0099a 0.0009a 
 Positive   
HER2 
 Negative 0.7381 0.7826 
 Positive   
Triple-negative 
 Yes 0.0026b 0.0021b 
 No   
LPBC 
 Yes <0.0001b <0.0001b 
 No   

aNegative association with PD-L1.

bPositive association with PD-L1.

PD-L1 in the epithelium and stroma correlates with pCR when measured as a continuous quantitative score (Fig. 1E and F; epithelial P = 0.0189; stromal P = 0.0050). Analysis of PD-L1 expression as a continuous quantitative score in subsets of patients that were LPBC-positive, hormone receptor–positive, HER2-amplified, and triple-negative revealed that PD-L1 in the epithelium and stroma correlates with pCR only in hormone-receptor–positive and HER2-amplified breast cancers, though analyses in LPBC and triple-negative subsets may be underpowered (Supplementary Fig. S3).

To determine a statistically rigorous cut point of PD-L1 using the continuous quantitative data, Joinpoint software was used (31). As shown in Supplementary Fig. S4, Joinpoint identified three differential points in the distribution of both epithelial (Supplementary Fig. S4A) and stromal (Supplementary Fig. S4B) continuous PD-L1 scores in the cohort, one of which located to the approximate visual threshold of PD-L1 positivity. We dichotomized PD-L1 at Joinpoint#2 for both epithelial and stromal scores, making approximately 30% of the cohort positive for PD-L1 expression. When dichotomized by Joinpoint software in this way, PD-L1 correlates with pCR (Supplementary Fig. S5; epithelial χ2P = 0.0595; stromal χ2P = 0.0499).

Examples of H&E images from cases scored as non-LPBC, with TIL infiltrate less than 50%, and LPBC, with TIL infiltrate ≥50%, are shown in Fig. 2A and B. Figure 2A represents a case scored as non-LPBC, with the two pathologists scoring <5 and 20% TIL component, respectively. Figure 2C is an image of PD-L1 staining, showing little to no reactivity in the same TIL-low case. Figure 2B represents a case scored as LPBC, with the two pathologists scoring 70% and 80% TIL component. PD-L1 staining in the same case shows robust expression in the stroma (Fig. 2D). Further, PD-L1 as measured on the entire cohort as a continuous quantitative score in the epithelium and stroma positively correlates with high TIL component (Fig. 3A and B; epithelial P < 0.0001; stromal P = 0.0001).

Figure 2.

Examples of PD-L1 expression in breast cancers with low and high levels of TILs. Representative H&E and PD-L1 in cases scored non-LPBC (A and C) and LPBC (B and D). Immunofluorescence images: blue, DAPI; green, pan-cytokeratin; red, PD-L1.

Figure 2.

Examples of PD-L1 expression in breast cancers with low and high levels of TILs. Representative H&E and PD-L1 in cases scored non-LPBC (A and C) and LPBC (B and D). Immunofluorescence images: blue, DAPI; green, pan-cytokeratin; red, PD-L1.

Close modal
Figure 3.

Quantitative assessment of epithelial and stromal PD-L1 expression correlates with TILs. Box plots with continuous PD-L1 scores in epithelium (A) and stroma (B) on the y axis and LPBC or non-LPBC on the x axis.

Figure 3.

Quantitative assessment of epithelial and stromal PD-L1 expression correlates with TILs. Box plots with continuous PD-L1 scores in epithelium (A) and stroma (B) on the y axis and LPBC or non-LPBC on the x axis.

Close modal

Univariate analyses using logistic regression identified node status and LPBC as predictors of pCR to neoadjuvant chemotherapy (Table 3). Multivariate analyses, including age, nodal status, tumor size, molecular classification (hormone therapy sensitive, HER2 positive, or triple-negative), nuclear grade, Ki-67 AQUA score, TILs, and PD-L1 expression, revealed that both epithelial and stromal PD-L1 expressions are nearing significance in predicting pCR in this multivariate analysis (Table 4). The P value becomes significant for both if Ki-67 is excluded from the analysis (Supplementary Table S1).

Table 3.

Univariate analysis of likelihood of pCR on Yale neoadjuvant cohort

Univariate
VariableORa (95% CI)P value
Age (y) 
 <50  
 ≥50 0.923 (0.377–2.259) 0.8609 
Nodal status 
 Positive 0.126 (0.042–0.379) 0.0002 
 Negative  
Tumor size (cm) 
 <2  
 >2 0.383 (0.101–1.449) 0.1575 
Nuclear grade 
 1–2  
 3 2.485 (0.983–6.278) 0.0543 
Estrogen receptor 
 Negative  
 Positive 0.488 (0.197–1.209) 0.1211 
Progesterone receptor 
 Negative  
 Positive 0.415 (0.167–1.030) 0.0579 
HER2 
 Negative  
 Positive 2.357 (0.901–6.167) 0.0806 
Triple-negative 
 No  
 Yes 2.222 (0.831–5.941) 0.1115 
LPBC 
 No  
 Yes 4.697 (1.036–21.302) 0.0449 
Epithelial PD-L1 
 Low  
 High 2.667 (0.945–7.528) 0.0639 
Stromal PD-L1 
 Low  
 High 2.804 (0.982–8.006) 0.0541 
Univariate
VariableORa (95% CI)P value
Age (y) 
 <50  
 ≥50 0.923 (0.377–2.259) 0.8609 
Nodal status 
 Positive 0.126 (0.042–0.379) 0.0002 
 Negative  
Tumor size (cm) 
 <2  
 >2 0.383 (0.101–1.449) 0.1575 
Nuclear grade 
 1–2  
 3 2.485 (0.983–6.278) 0.0543 
Estrogen receptor 
 Negative  
 Positive 0.488 (0.197–1.209) 0.1211 
Progesterone receptor 
 Negative  
 Positive 0.415 (0.167–1.030) 0.0579 
HER2 
 Negative  
 Positive 2.357 (0.901–6.167) 0.0806 
Triple-negative 
 No  
 Yes 2.222 (0.831–5.941) 0.1115 
LPBC 
 No  
 Yes 4.697 (1.036–21.302) 0.0449 
Epithelial PD-L1 
 Low  
 High 2.667 (0.945–7.528) 0.0639 
Stromal PD-L1 
 Low  
 High 2.804 (0.982–8.006) 0.0541 

NOTE: Bold P values are highlighted as being less than 0.05.

aOR indicates likelihood of pCR.

Table 4.

Multivariate analysis of likelihood of pCR on Yale neoadjuvant cohort

Multivariate (TIL)Multivariate (epithelial PD-L1)Multivariate (stromal PD-L1)
VariableORa (95% CI)P valueORa (95% CI)P valueORa (95% CI)P value
Age (y) 
 <50    
 ≥50 0.482 (0.086–2.711) 0.4077 0.164 (0.011–2.329) 0.1815 0.147 (0.009–2.290) 0.1713 
Nodal status 
 Negative    
 Positive 0.026 (0.003–0.219) 0.0008 0.009 (0.0003–0.254) 0.0057 0.009 (0.0003–0.284) 0.0072 
Tumor size (cm) 
 <2    
 >2 0.110 (0.014–0.841) 0.0334 0.054 (0.004–0.676) 0.0236 0.072 (0.006–0.915) 0.0425 
Nuclear grade 
 1–2    
 3 4.780 (0.766–29.822) 0.0939 2.920 (0.225–37.903) 0.4126 3.105 (0.238–40.520) 0.3872 
Molecular category 
 HR+/HER2    
 HER2+ 2.908 (0.390–21.704) 0.2979 3.624 (0.172–76.513) 0.4079 4.156 (0.206–83.696) 0.3524 
 Triple-negative 0.614 (0.061–6.167) 0.6783 0.705 (0.046–10.724) 0.8013 0.864 (0.057–13.138) 0.9160 
Ki-67 
 Low    
 High 6.186 (0.838–45.651) 0.0739 7.945 (0.530–119.073) 0.1335 8.550 (0.570–128.29) 0.1204 
LPBC 
 No      
 Yes 6.019 (0.035–1029.2) 0.4938     
Epithelial PD-L1 
 Low      
 High   11.120 (0.870–142.159) 0.0639   
Stromal PD-L1 
 Low      
 High     11.267 (0.779–162.95) 0.0756 
Multivariate (TIL)Multivariate (epithelial PD-L1)Multivariate (stromal PD-L1)
VariableORa (95% CI)P valueORa (95% CI)P valueORa (95% CI)P value
Age (y) 
 <50    
 ≥50 0.482 (0.086–2.711) 0.4077 0.164 (0.011–2.329) 0.1815 0.147 (0.009–2.290) 0.1713 
Nodal status 
 Negative    
 Positive 0.026 (0.003–0.219) 0.0008 0.009 (0.0003–0.254) 0.0057 0.009 (0.0003–0.284) 0.0072 
Tumor size (cm) 
 <2    
 >2 0.110 (0.014–0.841) 0.0334 0.054 (0.004–0.676) 0.0236 0.072 (0.006–0.915) 0.0425 
Nuclear grade 
 1–2    
 3 4.780 (0.766–29.822) 0.0939 2.920 (0.225–37.903) 0.4126 3.105 (0.238–40.520) 0.3872 
Molecular category 
 HR+/HER2    
 HER2+ 2.908 (0.390–21.704) 0.2979 3.624 (0.172–76.513) 0.4079 4.156 (0.206–83.696) 0.3524 
 Triple-negative 0.614 (0.061–6.167) 0.6783 0.705 (0.046–10.724) 0.8013 0.864 (0.057–13.138) 0.9160 
Ki-67 
 Low    
 High 6.186 (0.838–45.651) 0.0739 7.945 (0.530–119.073) 0.1335 8.550 (0.570–128.29) 0.1204 
LPBC 
 No      
 Yes 6.019 (0.035–1029.2) 0.4938     
Epithelial PD-L1 
 Low      
 High   11.120 (0.870–142.159) 0.0639   
Stromal PD-L1 
 Low      
 High     11.267 (0.779–162.95) 0.0756 

NOTE: Bold P values are highlighted as being less than 0.05.

Abbreviation: HR, hormone receptor.

aOR indicates likelihood of pCR.

PD-L1 expression in the epithelium or stroma as a continuous quantitative score or dichotomized into negative and positive predicts pCR. Dichotomized scores of epithelial or stromal expression are not independent of one another, though individually both are predictive of pCR in a multivariate model, including age, nodal status, tumor size, hormonal receptor status, HER2, and triple-negative status. PD-L1 expression in the epithelium and stroma is associated with ER-negative and progesterone receptor (PR)–negative, triple-negative, and LPBC breast cancers.

The evaluation of PD-L1 expression is challenging due to heterogeneity in expression and nonreproducibility of antibody reagents (32). Only a few studies have described PD-L1 protein expression in patients with breast cancer, finding an association of PD-L1 with proliferative markers and FoxP3 T-regulatory cells (18, 33). We found that roughly 30% of breast cancers express PD-L1 in the epithelium and/or stroma. The heterogeneous nature of PD-L1 expression in breast cancer, being both epithelial and/or stromal as well as present in only select fields of view, supports previous studies that suggest its expression is limited to specific regions of the tumor such as the invasive front (17). Previous studies have also described localization in both epithelial cells and on specific cells in the stroma. Our work confirms expression in both compartments, but we found no differences in epithelial or stromal expression with respect to prediction of pCR or association with various clinical features.

In addition to predicting response to neoadjuvant chemotherapy, we investigated whether PD-L1 expression may also be a biomarker for predicting response to immune therapies targeting the PD-1–PD-L1 pathway. Indeed, a landmark phase I trial of an anti–PD-1 antibody on various solid tumor types showed that the only patients with objective response to anti–PD-1 therapy were those whose tumors expressed PD-L1 (15). A limitation thus far in the assessment of PD-L1 protein expression has been a lack of specific and reproducible antibodies for use on formalin-fixed paraffin-embedded tissue. Of the eight commercially available antibodies, only three, including the E1L3N clone from Cell Signaling Technology, used in this study passed our quality control (data not shown). In previous studies, we have examined PD-L1 expression using clone 5H1 through a collaboration with Lieping Chen (Department of Immunobiology, Yale University; ref. 32). The antibody used in this study to detect PD-L1 shows similar, but not identical staining compared with clone 5H1 (see Supplementary Fig. S2C–S2E). This clone is commercially available, uses a more standardized staining protocol, and yields reproducible results.

Although we believe this study illustrates the value of PD-L1 as a potential prognostic marker, there are a number of limitations. This is a retrospective study comprising a modest sample size of patients from a single institution. Although the majority of patients (76.6%) received anthracycline-based neoadjuvant therapy, the remainder received variations of taxanes and/or carboplatin, limiting any interpretation of treatment-specific results. Another limitation is that only a single monoclonal antibody was used to assess PD-L1 expression. Although this is considered acceptable in a publication, efforts are under way to evaluate other validated antibodies. We have measured PD-L1 mRNA expression on tissue microarrays containing a large number of breast cancer samples, although not yet on this neoadjuvant cohort (34). Furthermore, the potential information contained within the dynamic nature of PD-L1 expression (heterogeneity of localization as well as intensity of expression) may be oversimplified by our methods of analysis. Finally, this cohort is too recent to provide a meaningful prognostic evaluation of PD-L1 expression. Although TILs have been reported to be prognostic and PD-L1 expression correlates with the presence of TILs, the prognostic value of PD-L1 will be assessed in future work.

The close correlation of PD-L1 with TILs and the ease with which PD-L1 can be induced by expression of inflammatory cytokine IFNγ suggest that PD-L1 may act as a surrogate marker for an antitumor immune response, albeit one that is being downregulated. The presence of both TILs and PD-L1 in the tumor microenvironment could indicate an adaptive immune resistance to endogenous antitumor activity, suggesting that patients with both of these components would benefit from immunotherapy (17).

In summary, we demonstrate a reproducible assay for evaluating PD-L1 protein expression on formalin-fixed, paraffin-embedded tissue sections that predicts response to neoadjuvant chemotherapy, independent of localization and treatment. PD-L1 expression also correlates with the presence of TILs. However, the value of this marker will be its use in patients treated with PD-L1 axis–directed therapies. In the future, we look forward to using these same reagents and methods to evaluate treated patients.

D.L. Rimm reports receiving commercial research grants from Genoptix, Kolltan, and Gilead. He is also a consultant/advisory board member for Novartis, BMS, Genoptix, and ACD, Inc. No potential conflicts of interest were disclosed by the other authors.

Conception and design: H. Wimberly, K. Schalper, V. Bossuyt, D.L. Rimm

Development of methodology: H. Wimberly, J.R. Brown, K. Schalper, V. Bossuyt, D.L. Rimm

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): H. Wimberly, J.R. Brown, M.R. Silver, C. Nixon, V. Bossuyt, D.R. Lannin

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): H. Wimberly, J.R. Brown, V. Bossuyt, D.L. Rimm

Writing, review, and/or revision of the manuscript: H. Wimberly, J.R. Brown, K. Schalper, H. Haack, M.R. Silver, C. Nixon, V. Bossuyt, L. Pusztai, D.R. Lannin, D.L. Rimm

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): H. Wimberly, H. Haack, D.R. Lannin, D.L. Rimm

Study supervision: D.L. Rimm

This work was supported by the Breast Cancer Research Foundation.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Rastogi
P
,
Anderson
SJ
,
Bear
HD
,
Geyer
CE
,
Kahlenberg
MS
,
Robidoux
A
, et al
Preoperative chemotherapy: updates of National Surgical Adjuvant Breast and Bowel Project Protocols B-18 and B-27
.
J Clin Oncol
2008
;
26
:
778
85
.
2.
Guarneri
V
,
Broglio
K
,
Kau
SW
,
Cristofanilli
M
,
Buzdar
AU
,
Valero
V
, et al
Prognostic value of pathologic complete response after primary chemotherapy in relation to hormone receptor status and other factors
.
J Clin Oncol
2006
;
24
:
1037
44
.
3.
Brown
JR
,
DiGiovanna
MP
,
Killelea
B
,
Lannin
DR
,
Rimm
DL
. 
Quantitative assessment Ki-67 score for prediction of response to neoadjuvant chemotherapy in breast cancer
.
Lab Invest
2014
;
94
:
98
106
.
4.
Loi
S
,
Sirtaine
N
,
Piette
F
,
Salgado
R
,
Viale
G
,
Van Eenoo
F
, et al
Prognostic and predictive value of tumor-infiltrating lymphocytes in a phase III randomized adjuvant breast cancer trial in node-positive breast cancer comparing the addition of docetaxel to doxorubicin with doxorubicin-based chemotherapy: BIG 02–98
.
J Clin Oncol
2013
;
31
:
860
7
.
5.
Issa-Nummer
Y
,
Darb-Esfahani
S
,
Loibl
S
,
Kunz
G
,
Nekljudova
V
,
Schrader
I
, et al
Prospective validation of immunological infiltrate for prediction of response to neoadjuvant chemotherapy in HER2-negative breast cancer - A substudy of the Neoadjuvant GeparQuinto Trial
.
PLoS ONE
2013
;
8
:
e79775
.
6.
Denkert
C
,
Loibl
S
,
Noske
A
,
Roller
M
,
Muller
BM
,
Komor
M
, et al
Tumor-associated lymphocytes as an independent predictor of response to neoadjuvant chemotherapy in breast cancer
.
J Clin Oncol
2010
;
28
:
105
13
.
7.
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
31
.
8.
Gianni
L
,
Zambetti
M
,
Clark
K
,
Baker
J
,
Cronin
M
,
Wu
J
, et al
Gene expression profiles in paraffin-embedded core biopsy tissue predict response to chemotherapy in women with locally advanced breast cancer
.
J Clin Oncol
2005
;
23
:
7265
77
.
9.
Ono
M
,
Tsuda
H
,
Shimizu
C
,
Yamamoto
S
,
Shibata
T
,
Yamamoto
H
, et al
Tumor-infiltrating lymphocytes are correlated with response to neoadjuvant chemotherapy in triple-negative breast cancer
.
Breast Cancer Res Treat
2012
;
132
:
793
805
.
10.
Yamaguchi
R
,
Tanaka
M
,
Yano
A
,
Tse
GM
,
Yamaguchi
M
,
Koura
K
, et al
Tumor-infiltrating lymphocytes are important pathologic predictors for neoadjuvant chemotherapy in patients with breast cancer
.
Hum Pathol
2012
;
43
:
1688
94
.
11.
Mahmoud
SM
,
Paish
EC
,
Powe
DG
,
Macmillan
RD
,
Grainge
MJ
,
Lee
AH
, et al
Tumor-infiltrating CD8+ lymphocytes predict clinical outcome in breast cancer
.
J Clin Oncol
2011
;
29
:
1949
55
.
12.
Seo
AN
,
Lee
HJ
,
Kim
EJ
,
Kim
HJ
,
Jang
MH
,
Lee
HE
, et al
Tumour-infiltrating CD8+ lymphocytes as an independent predictive factor for pathological complete response to primary systemic therapy in breast cancer
.
Br J Cancer
2013
;
109
:
2705
13
.
13.
Lee
HJ
,
Seo
JY
,
Ahn
JH
,
Ahn
SH
,
Gong
G
. 
Tumor-associated lymphocytes predict response to neoadjuvant chemotherapy in breast cancer patients
.
J Breast Cancer
2013
;
16
:
32
9
.
14.
Brahmer
JR
,
Tykodi
SS
,
Chow
LQ
,
Hwu
WJ
,
Topalian
SL
,
Hwu
P
, et al
Safety and activity of anti-PD-L1 antibody in patients with advanced cancer
.
N Engl J Med
2012
;
366
:
2455
65
.
15.
Topalian
SL
,
Hodi
FS
,
Brahmer
JR
,
Gettinger
SN
,
Smith
DC
,
McDermott
DF
, et al
Safety, activity, and immune correlates of anti-PD-1 antibody in cancer
.
N Engl J Med
2012
;
366
:
2443
54
.
16.
Dong
HD
,
Strome
SE
,
Salomao
DR
,
Tamura
H
,
Hirano
F
,
Flies
DB
, et al
Tumor-associated B7-H1 promotes T-cell apoptosis: a potential mechanism of immune evasion
.
Nat Med
2002
;
8
:
793
800
.
17.
Taube
JM
,
Anders
RA
,
Young
GD
,
Xu
HY
,
Sharma
R
,
McMiller
TL
, et al
Colocalization of inflammatory response with B7-H1 expression in human melanocytic lesions supports an adaptive resistance mechanism of immune escape
.
Sci Transl Med
2012
;
4
:
127ra37
.
18.
Ghebeh
H
,
Tulbah
A
,
Mohammed
S
,
Eikum
N
,
Bin Amer
SM
,
Al-Tweigeri
T
, et al
Expression of B7-H1 in breast cancer patients is strongly associated with high proliferative Ki-67-expressing tumor cells
.
Int J Cancer
2007
;
121
:
751
8
.
19.
Hamanishi
J
,
Mandai
M
,
Iwasaki
M
,
Okazaki
T
,
Tanaka
Y
,
Yamaguchi
K
, et al
Programmed cell death 1 ligand 1 and tumor-infiltrating CD8+ T lymphocytes are prognostic factors of human ovarian cancer
.
Proc Natl Acad Sci U S A
2007
;
104
:
3360
5
.
20.
Iwai
Y
,
Ishida
M
,
Tanaka
Y
,
Okazaki
T
,
Honjo
T
,
Minato
N
. 
Involvement of PD-L1 on tumor cells in the escape from host immune system and tumor immunotherapy by PD-L1 blockade
.
Proc Natl Acad Sci U S A
2002
;
99
:
12293
7
.
21.
Hirano
F
,
Kaneko
K
,
Tamura
H
,
Dong
HD
,
Wang
SD
,
Ichikawa
M
, et al
Blockade of B7-H1 and PD-1 by monoclonal antibodies potentiates cancer therapeutic immunity
.
Cancer Res
2005
;
65
:
1089
96
.
22.
Azuma
T
,
Yao
S
,
Zhu
G
,
Flies
AS
,
Flies
SJ
,
Chen
L
. 
B7-H1 is a ubiquitous antiapoptotic receptor on cancer cells
.
Blood
2008
;
111
:
3635
43
.
23.
Topalian
SL
,
Hodi
FS
,
Brahmer
JR
,
Gettinger
SN
,
Smith
DC
,
McDermott
DF
, et al
Safety, activity, and immune correlates of anti-PD-1 antibody in cancer
.
N Engl J Med
2012
;
366
:
2443
54
.
24.
Topalian
SL
,
Sznol
M
,
McDermott
DF
,
Kluger
HM
,
Carvajal
RD
,
Sharfman
WH
, et al
Survival, durable tumor remission, and long-term safety in patients with advanced melanoma receiving nivolumab
.
J Clin Oncol
2014
;
32
:
1020
30
.
25.
Brahmer
JR
,
Tykodi
SS
,
Chow
LQ
,
Hwu
WJ
,
Topalian
SL
,
Hwu
P
, et al
Safety and activity of anti-PD-L1 antibody in patients with advanced cancer
.
N Engl J Med
2012
;
366
:
2455
65
.
26.
Herbst
RS
,
Soria
JC
,
Kowanetz
M
,
Fine
GD
,
Hamid
O
,
Gordon
MS
, et al
Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients
.
Nature
2014
;
515
:
563
7
.
27.
Taube
JM
,
Klein
A
,
Brahmer
JR
,
Xu
H
,
Pan
X
,
Kim
JH
, et al
Association of PD-1, PD-1 ligands, and other features of the tumor immune microenvironment with response to anti-PD-1 therapy
.
Clin Cancer Res
2014
;
20
:
5064
74
.
28.
Hamid
O
,
Robert
C
,
Daud
A
,
Hodi
FS
,
Hwu
WJ
,
Kefford
R
, et al
Safety and tumor responses with lambrolizumab (anti-PD-1) in melanoma
.
N Engl J Med
2013
;
369
:
134
44
.
29.
Camp
RL
,
Chung
GG
,
Rimm
DL
. 
Automated subcellular localization and quantification of protein expression in tissue microarrays
.
Nat Med
2002
;
8
:
1323
7
.
30.
Moeder
CB
,
Giltnane
JM
,
Moulis
SP
,
Rimm
DL
. 
Quantitative, fluorescence-based in-situ assessment of protein expression
.
Methods Mol Biol
2009
;
520
:
163
75
.
31.
Kim
HJ
,
Fay
MP
,
Feuer
EJ
,
Midthune
DN
. 
Permutation tests for joinpoint regression with applications to cancer rates
.
Stat Med
2000
;
19
:
335
51
.
32.
Velcheti
V
,
Schalper
KA
,
Carvajal
DE
,
Anagnostou
VK
,
Syrigos
KN
,
Sznol
M
, et al
Programmed death ligand-1 expression in non-small cell lung cancer
.
Lab Invest
2014
;
94
:
107
16
.
33.
Ghebeh
H
,
Barhoush
E
,
Tulbah
A
,
Elkum
N
,
Al-Tweigeri
T
,
Dermime
S
. 
FOXP3(+) Tregs and B7-HI+/PD-I+T lymphocytes co-infiltrate the tumor tissues of high-risk breast cancer patients: implication for immunotherapy
.
BMC Cancer
2008
;
8
:
57
.
34.
Schalper
KA
,
Velcheti
V
,
Carvajal
D
,
Wimberly
H
,
Brown
J
,
Pusztai
L
, et al
In situ tumor PD-L1 mRNA expression is associated with increased TILs and better outcome in breast carcinomas
.
Clin Cancer Res
2014
;
20
:
2773
82
.