NKG2D ligand tumor expression and association with clinical outcome in early breast cancer patients: an observational study

Background Cell surface NKG2D ligands (NKG2DL) bind to the activating NKG2D receptor present on NK cells and subsets of T cells, thus playing a role in initiating an immune response. We examined tumor expression and prognostic effect of NKG2DL in breast cancer patients. Methods Our study population (n = 677) consisted of all breast cancer patients primarily treated with surgery in our center between 1985 and 1994. Formalin-fixed paraffin-embedded tumor tissue was immunohistochemically stained with antibodies directed against MIC-A/MIC-B (MIC-AB), ULBP-1, ULBP-2, ULBP-3, ULBP-4, and ULBP-5. Results NKG2DL were frequently expressed by tumors (MIC-AB, 50% of the cases; ULBP-1, 90%; ULBP-2, 99%; ULBP-3, 100%; ULBP-4, 26%; ULBP-5, 90%) and often showed co-expression: MIC-AB and ULBP-4 (p = 0.043), ULBP-1 and ULBP-5 (p = 0.006), ULBP-4 and ULBP-5 (p < 0.001). MIC-AB (p = 0.001) and ULBP-2 (p = 0.006) expression resulted in a statistically significant longer relapse free period (RFP). Combined expression of these ligands showed to be an independent prognostic parameter for RFP (p < 0.001, HR 0.41). Combined expression of all ligands showed no associations with clinical outcome. Conclusions We demonstrated for the first time that NKG2DL are frequently expressed and often co-expressed in breast cancer. Expression of MIC-AB and ULBP-2 resulted in a statistically significant beneficial outcome concerning RFP with high discriminative power. Combination of all NKG2DL showed no additive or interactive effect of ligands on each other, suggesting that similar and co-operative functioning of all NKG2DL can not be assumed. Our observations suggest that among driving forces in breast cancer outcome are immune activation on one site and tumor immune escape on the other site.


Background
Breast cancer is the most commonly diagnosed female cancer and is the leading cause of death from cancer in women in the western world [1]. Decisions regarding use of systemic therapy are mainly based on prognostic and predictive factors like lymph node status, tumor size, grade, hormone receptor and human epidermal growth factor receptor 2 (HER2) expression [2,3]. However, current prognostic and predictive factors still do not provide optimal risk-stratification. Therefore, additional prognostic and predictive information could result in an improved tailored treatment for patients with breast cancer.
There is strong evidence that the immune system plays a role in tumor growth and progression [4,5]. An effective immune response may lead to recognition of tumor cells, resulting in their eradication. However, due to their genetic unstable nature, tumor cells may arise which display properties that enables them to escape from immune recognition [4,5]. Indeed, downregulation or loss of proteins that are crucial for immune responses, like classical human leukocyte antigens (HLA) class I, or upregulation of proteins that confer resistance to immune recognition, like non-classical HLA class I, are frequently found in various types of tumors [6][7][8][9][10].
The activating receptor natural killer cell lec-tin-like receptor gene 2D (NKG2D) is a stimulatory immune receptor that is expressed on natural killer (NK) cells, NKT cells, γδ + T cells and CD8+ T cells [11]. Ligands which bind NKG2D receptors comprise major histocompatibility complex class I chain-related proteins A and B (MIC-AB) and unique long 16 (UL16) binding proteins 1-6 (ULBP1-6) [12,13]. Expression of these ligands may be induced upon infection and other inducers of cellular stress and is unusual in normal cells [14]. By binding to the NKG2D receptors on NK and T cells, the NKG2D ligands may initiate an immune response against cells expressing these ligands. Overexpression and shedding of NKG2D ligands have been reported [14]. It is, however, unclear whether these features also results in activation of an immune response or lead to overstimulation and downregulation of NKG2D on immune cells [11].
Malignant transformation of cells may be among stimuli inducing expression of NKG2D ligands as such expression has been found in various tumor types [8][9][10][15][16][17][18]. This may be a mechanism for preventing tumor growth by advancing an anti-tumor immune response. Convincing evidence has been found in in vivo studies, which have shown that in mouse models transfection with NKG2D ligands resulted in a NKG2Dmediated tumor rejection [19,20]. Other studies showed that downregulation or complete knockout of NKG2D in mice resulted in an impaired immune response against tumor cells, higher expression levels of NKG2D ligands, and an increased incidence of certain tumors [21,22].
A few studies have investigated tumor expression of NKG2D ligands and associations with clinical outcome in human breast, colorectal, and ovarian cancer [8][9][10]15,16]. Expression of MIC-A was frequently found in all tumors studied and resulted in a statistically significant favorable patient's prognosis in colorectal cancer, while it was not statistically significantly associated with outcome in breast cancer and ovarian cancer [8][9][10]16]. ULBP1-5 expression was also found to be expressed in many tumor samples of colorectal and ovarian cancer [9,10,15]. In colorectal cancer expression of ULBP5 was an independent prognostic factor for a favorable clinical outcome [9]. In contrast to these results, expression of ULBP2 and ULBP4 were found to be independent prognostic factors for a worse outcome of ovarian cancer patients [10,15]. Taken together, several studies suggest that evasion of NKG2D-mediated immune regulation plays an important role in tumor progression, but some studies contradict this suggestion. Contradictory results may be explained by assuming functional differences in immune regulation of the different ligands. Moreover, expression of NKG2D ligands may behave different among different tumor types [9]. It is known that overexpression or shedding of these ligands leads to overstimulation and downregulation of NKG2D on immune cells [10,15], thereby evading an immune response.
In breast cancer, the prognostic effect of NKG2D ligands and their mutual relationship is largely unknown. Therefore, the purpose of this study was to analyze the clinical prognostic value of MIC-AB and ULBP1-5 in a large patient cohort of early stage breast cancer.

Patients and tumors
The patient population comprised all non-metastasized breast cancer patients primarily treated with surgery between 1985 and 1994 at the Leiden University Medical Center (n = 677). Patients with bilateral tumors or a prior history of cancer, other than basal cell carcinoma or cervical carcinoma in situ, were excluded. The following data were known: age, tumor morphology and differentiation grade, TNM stage, type of local and systemic therapy, recurrence and survival status, estrogen receptor (ER), progesterone receptor (PgR), and human epidermal growth factor 2 (HER2) expression (Table 1). All these parameters were determined according to current pathology standards. A tissue micro array (TMA) of available formalin-fixed paraffin-embedded (FFPE) tumors of the patient cohort has been previously constructed and described (n = 574) [23]. Approval was obtained from the Leiden University Medical Center Medical Ethics Committee. All samples were handled in a coded fashion, according to National ethical guidelines ("Code for Proper Secondary Use of Human Tissue", Dutch Federation of Medical Scientific Societies).
TMA sections of 4 μm were cut, deparaffinized and rehydrated. Endogenous peroxidase was blocked in 0.3% hydrogen-peroxide methanol for 20 minutes. Heatinduced antigen retrieval for 10 minutes at maximum power in a microwave oven was performed. Sections were incubated overnight with primary antibodies using predetermined optimal dilutions and incubations times. Sections for ULBP-2 staining were incubated with Rabbit Anti-Goat Immunoglobulins (DAKO) followed by StreptABComplex (DAKO) for 30 minutes. Sections for all other stainings were incubated with secondary antibody Envision (Dako cytomation K4001 or K4003) for 30 minutes. Stainings were visualized using DABsolution (Dako cytomation K3468), counterstained with haematoxylin, dehydrated, and finally mounted in malinol. For each type of antibody, all tissue sections were stained simultaneously to avoid inter-assay variation.

Statistical analysis
Statistical analyses were performed using the statistical package SPSS (version 16.0 for Windows, Spps Inc, Chicago, IL, USA). Cohen's kappa coefficient was used to assess inter-observer agreement in quantification. This revealed a moderate agreement for ULBP-5 (kappa = 0.410), a substantial agreement in classification for MIC-AB (kappa = 0.790) and ULBP-4 (kappa = 0.650), and an almost perfect agreement for ULBP-1 (kappa = 0.913), ULBP-2 (kappa = 0.940), and ULBP-3 (kappa = 0.869). The χ 2 test was used to evaluate associations between expression of the different NKG2D ligands. Relapse free period (RFP) was the time from date of surgery until an event (locoregional recurrence and/or a distance recurrence, whichever came first). The Kaplan-Meier method was used for survival plotting and log-rank test for comparison of survival curves. RFP is reported as cumulative incidence function, after accounting for death as competing risk [26]. Cox regression was used for univariate and multivariate analysis for RFP. Significant variables (p < 0.1) in univariate analysis were included in multivariate analysis.

Patient and tumor characteristics
Median age of patients was 57 years (range: 23-96 years). Median follow-up of patients alive was 19 years (range: 14-23 years). Clinicopathological and treatment characteristics are shown in table 1.

Association of NKG2D ligands with clinicopathological parameters
High expression of NKG2D ligands was generally associated with favorable clinicopathological parameters (table 1 and 2): statistically significant associations were found between high expression of MIC-AB and lower tumor grade (p = 0.012); high expression of ULBP-1 and higher tumor grade (p < 0.001), smaller tumor size (p < 0.001) and more lymph node positive tumors (p = 0.049); high expression of ULBP-2 and younger age (p = 0.022), lower tumor grade (p < 0.001), smaller tumor size (p = 0.005) and more lymph node negative tumors (p = 0.046); high expression of ULBP-3 and higher tumor grade (p = 0.001); high expression of ULBP-4 and smaller tumor size (p = 0.001); high expression of ULBP-5 and more PgR negative tumor status (p = 0.016).  Table 3).

Associations with outcome of NKG2D ligands
In order to analyze the frequencies and prognostic effect of number of co-expressed and amount of coexpression of NKG2D ligands, two new variables were constructed. First, the total number of the different NKG2D ligands that were expressed. For that purpose, the number of NKG2D ligands with high expression was counted. So for each tumor, this resulted in a minimal and maximal possible score of respectively 0 and 6. Second, the total amount of NKG2D ligand expression. For that purpose, the intensity of staining (ranging from 0 to 3) of NKG2D ligands was added, obtaining a total NKG2D ligand intensity score. So for each tumor, this resulted in a minimal and maximal possible score of respectively 0 and 18. The median number of NKG2D ligands with high expression was 1 (range 0-6). For statistical reasons (too small patient groups) in outcome analyses, the groups with 3, 4, 5 and 6 numbers of different NKG2D ligands highly expressed were combined as one single group: ≥ 3 ligands of high expression.
No associations were seen for the number of NKG2D ligands with high expression for RFP outcome analyses (log rank p-value: 0.967); patients with tumors with a low number of NKG2D ligands with high expression resulted in a similar RFP compared to a high number of NKG2D ligands with high expression ( Figure 4A).
For outcome analyses, NKG2D scores 14-16 were combined and classified as ≥ 14, since these subgroups separately contained only one patient. No association was seen for amount of total NKG2D ligand expression and RFP (log rank p-value: 0.721); high total NKG2D ligand expression resulted in some cases in a worse RFP (e.g. score 11) while in others it resulted in a favorable RFP (e.g. score 13) and vice versa, low total NKG2D ligand expression resulted for some patients in a favorable RFP (e.g. score 6) and for other patients in a worse RFP (e.g. score 4) ( Figure 4B).

Discussion
The importance of interaction between tumor development and the immune system for cancer outcome is highlighted by an overwhelming number of studies, performed in vitro, in vivo and using patient cohorts. Recent studies have shown that NKG2D ligands may play an important role in cancer immunosurveillance and cancer immunoediting [8][9][10][15][16][17][18][19][20][21][22]. In this study, we examined the impact of tumor expression of NKG2D ligands on the prognosis of breast cancer patients. The data of our study indicate that NKG2D ligands are frequently high expressed in breast tumors and that this expression influences prognosis of patients. We were able to statistically prove that high expression levels of MIC-AB and ULBP-2 resulted in a RFP benefit. Combining expression of MIC-AB and ULBP-2 resulted in a very accurate stratification of patients for prognosis  concerning RFP. The prognostic potential of this combined variable was comparable to that of lymph node status: patients with low tumor expression of both ligands had an almost 2.5 times increased risk of developing relapses compared to patients with high tumor expression of both ligands.
NKG2D ligands are expressed on the cell surface in response to stress or malignant transformation [11]. Our study confirms that breast cancer tumor cells show frequent and high expression of NKG2D ligands, as has been found in other studies for various types of tumors such as ovarian cancer, colorectal cancer and breast cancer [8][9][10]15,16]. Though all studies show consistently frequent expression of NKG2D ligands, very diverse prognostic effects have been described for these types of ligands in different cancer types [8][9][10]15,16]. This may be explained by functional differences in immune regulation for varying expression levels of different ligands in different environments. Expression of NKG2D ligands may induce an immune response through binding to the NKG2D receptor, present on NK cells and a subset of T cells [11]. Therefore, selective outgrowth of malignant cells that do not express these NKG2D ligands may be a mechanism of tumor immune escape. On the other hand, overexpression of NKG2D ligands could lead to overstimulation and thereby insensibility or anergy of immune cells, which would result in evasion of immune attack by tumors overexpressing NKG2D ligands [11]. Adding to this hypothesis, it has been reported that NKG2D ligands on the cell membrane may be cleaved and produce soluble molecules. This shedding of NKG2D ligands could systemically downregulate NKG2D receptor expression and thereby result in an impaired anti-tumor reactivity of NK and T cells [11,27]. Taken together, the mechanisms by which NKG2D ligands mediate immune function or dysfunction may be diverse in different tumors and differ according to circumstances. The contradictory results on the prognostic effect of NKG2D ligands found between different studies on different tumors may be reflected by the functional and mechanistic implications of interaction between NKG2D and its ligands. In ovarian cancer expression of NKG2D ligands resulted in a worse patient outcome, probably due to chronic MICAB:ULBP2 both low MICAB:ULBP2 either low MICAB:ULBP2 both high overexpression and shedding of these ligands, leading to overstimulation and downregulation of the NKG2D receptor of NK and T cells and, therefore, an impaired immune response [10,15]. Supporting the hypothesis that elevated expression of NKG2D ligands results in immune escape in ovarian cancer, one study found elevated levels of MIC-AB and ULBP-2 to be positively correlated to less intra-tumor epithelial CD57+ cells. The results found in breast cancer in the present study are contradictory to the results found in ovarian cancer, but similar to those found in colorectal cancer [9,16]. The results in our study and colorectal cancer are supported by the theory that expression of NKG2D ligands results in activation of immune cells which is reflected in a patient beneficial outcome for high ligand expression [9,16]. We found frequent and high expression of ligands in our study and statistically significant associations between expression levels of these ligands, indicating their cooperation with each other. Adding to the hypothesis that low expression of these ligands is a result of selective pressure by the immune system that results in cancer immune evasion or immunoediting, low expression of MIC-AB and ULBP-2 were prognostic factors for an unfavorable RFP of patients. When expression of MIC-AB and ULBP-2 were combined they showed to add to each others prognostic effect which is in line with the results found in previous studies [9,10] and suggests that NKG2D ligands operate together and in a similar manner.
Since the exact functioning of all NKG2D ligands and their cooperative function is largely unknown, we performed outcome analyses with two different variables that represented combined number of highly coexpressed ligands and amount of co-expression of all ligands. The results of these analyses revealed no patterns of any cooperative functioning between all ligands, as both variables showed no consistent and significant relationship with clinical outcome of disease. This suggests that the original hypothesis of all NKG2D ligands having a similar functioning and additive effect on each other's functioning in activating or evading an immune response, may be too simplistic. Considering our results and those as found in literature, altogether, each NKG2D ligand analysed separately does not show equal effects on clinical outcome, and different ligands show varying prognostic effects in different tumors. Specific combinations of ligands (e.g. MIC-AB and ULBP-2 in our study, ULBP2 and ULBP4 in ovarian cancer [10]) do show additive effects or statistical interactions on prognostic value. However, as highlighted by our combined analyses, a simple additive effect of all NKG2D ligands, by considering a similar or cooperative functioning of all these ligands, can not be assumed. This indicates the complexity of NKG2D ligands functioning and emphasizes the importance of further research on the precise mechanisms of actions of NKG2D ligands, separately, in combination with each other, and under different circumstances.

Conclusions
We have shown in this study, for the first time, that breast tumors may express all of the known NKG2D ligands and that expression of MIC-AB and ULBP-2 results in a favorable outcome concerning RFP. A variable combining MIC-AB and ULBP-2 expression has shown to be a prognostic parameter independently of known clinicopathological parameters and with high discriminative power. Our results suggest that NKG2D ligands play a crucial role in tumor immunoediting in breast cancer and provide further evidence that tumorimmune interactions play an important role in breast cancer. In addition, by NKG2D ligand combined analyses we highlight the importance of further studies on unraveling the precise separate functioning of these ligands. Submit your next manuscript to BioMed Central and take full advantage of: