Expression of HLA-DR genes in gliomas: correlation with clinicopathological features and prognosis

Abstract

Background

Human leukocyte antigen (HLA)-DR is a classical major histocompatibility complex (MHC) class II molecule encoded by five genes: HLA-DRA, HLA-DRB1, HLA-DRB3, HLA-DRB4 and HLA-DRB5. The current study aimed to investigate the role of these genes in gliomas by analyzing microarray data.

Methods

We enrolled 305 patients with histologically confirmed gliomas, and performed microarray data analysis along with studying their clinical characteristics. A new variable, termed HLA-DR score, was defined to explain the expression information of all five HLA-DR genes by factor analysis. HLA-DR scores in each grade of glioma and normal brain tissue were compared using one-way ANOVA. Lastly, correlations of HLA-DR scores with progression-free survival (PFS) and overall survival (OS) were analyzed with Kaplan-Meier and Cox analysis.

Results

Our study indicated that an increased HLA-DR score, i.e. overexpression of HLA-DR genes, was correlated with a more aggressive glioma tumor grade (p < 0.001, One-way ANOVA). Moreover, the HLA-DR score was significantly higher in astrocytic tumors than oligodendroglial tumors (−0.718 ± 3.177 versus −2.975 ± 2.662, t-test) in low-grade gliomas (LGGs). Kaplan-Meier analysis of both PFS (p = 0.046, log-rank test; p = 0.021, Breslow test) and OS (p = 0.029, Breslow test) showed significant differences in the clinical outcomes between LGG patients with high versus low HLA-DR scores. Finally, the HLA-DR score was further identified to be an independent prognostic factor of clinical outcomes by multivariate analysis (p = 0.042 and p = 0.025, for PFS and OS, respectively) in LGG patients.

Conclusion

Expression of HLA-DR genes can be used to predict the tumor grade in gliomas, and the histological subtype in LGG. Furthermore, they are also an independent predictor for LGG patient survival.

Background

Glioma, which accounts for approximately 81% of all primary malignant brain tumors, is the most common and aggressive form of primary brain cancer in adults [1, 2]. According to the World Health Organization (WHO) classification [3], gliomas can be classified into four grades based on histopathological criteria: WHO grade I are often benign and curable with gross total resection; WHO grade II and grade III gliomas often display more aggressive biological behaviors and eventually progress to high grade malignancies; and WHO grade IV tumors, also referred to as glioblastomas (GBMs), are the most frequent and lethal form. Despite current treatments, the median survival of patients with GBMs remains dismal [4]. Additionally, gliomas can also be classified as either astrocytic tumors or oligodendroglial tumors, based on their histopathological appearance [3]. In recent decades, there has been a substantial increase in the general understanding of glioma development and progression. However, there are still many issues which need further investigation, including the role of the host immune system in glioma development/biology.

Tumor-related immune surveillance evasion, and host immune suppression, have been widely accepted as typical characteristics of human cancers, but how the host immune system affects tumor development and progression has remained a challenge in the field of immunology [5, 6]. In this context, Schreiber et al. have first proposed the concept of cancer immunoediting, and suggested that the immune system played a dual role in cancer development and progression, i.e. in both the suppression and promotion of cancer [7, 8]. Gliomas is a disease carries profound genomic alterations, and a variety of tumor-specific genetic changes have been identified as prognostic indicators and/or therapeutic targets in glioma [9,10,11]. However, there are few papers regarding the roles of immune-related genetic changes in gliomas [12].

The human major histocompatibility complex (MHC), i.e. human leukocyte antigen (HLA), comprises of a set of genes which play an essential role in immune modulation. It is a 4000 kb gene complex on chromosome 6 [13], and can be divided into three groups based on the DNA sequence: HLA class I, HLA class II, and HLA class III. Aberrant expression of HLA genes has been identified in many human cancers. For instance, overexpression of HLA-G (subtype of HLA class I) has been demonstrated in various human malignancies, such as melanoma, breast cancer, and ovarian cancer [14,15,16,17]. Moreover, correlation of HLA-G overexpression with poor prognosis in low-grade gliomas (LGGs) has also been demonstrated previously by our group [18].

HLA-DR is a classical MHC class II molecule, with a αβ heterodimer anchored in the membrane. Its genetics are complex, and it is encoded by several genes with different functions. The α-chain is encoded by the HLA-DRA gene, while the β-chain is mainly encoded by 4 genes, including HLA-DRB1, HLA-DRB3, HLA-DRB4 and HLA-DRB5 [19]. Recently, a study reported that glioma patients with high HLA-DR expression had a significantly lower survival rate [20]. This finding was unanticipated and piqued our interest. As HLA-DR has been known as the antigen which is most responsible for graft loss during the first 6 months, whereas HLA-B is more important during the first 2 years, and HLA-A is associated with long-term survival [21]. Overexpression of these molecules is supposed to increase tumor immunogenicity, and thus slow or even prevent tumor growth. However, this study suggested otherwise and showed that overexpression of HLA-DR led to decreased survival in glioma patients [20]. Additionally, this study had some limitations, such as a small cohort of 60 patients, and the survival analysis was performed in patients with all grades of gliomas grouped together. In our view, the conclusions of this study may not accurately reflect the differences of HLA-DR distribution in different glioma grades. Thus, in this study we aimed to assess the mRNA expression levels of five HLA-DR genes using a large cohort, and specifically explore their potential correlation with the clinical outcome in gliomas of each grade.

Methods

Study population

A total of 305 patients who had undergone neurosurgical operations at Beijing Tiantan Hospital were enrolled. All patients had histological confirmation of gliomas along with mRNA microarray analysis. Histopathologic evaluation was performed independently by two experienced neuropathologists [3]. The clinical information of all patients was collected from the Chinese Glioma Genome Atlas database (http://www.cgga.org.cn). Among 305 patients, only six were lost in follow-up evaluation, and no patient deaths were recorded due to other diseases or unexpected events. The average follow-up length was 42.8 ± 35.5 months (median, 27.4 months; range, 0.7-111.4 months). In addition, tissue from five normal brains were surgically resected from trauma patients for use as negative controls. The present study was approved by the Ethics Committee of Beijing, Tiantan Hospital.

RNA extraction and whole genome gene profiling

RNA isolation and microarray analysis was carried out as previously described [22]. Briefly, tissue samples from patients were immediately snap frozen in liquid nitrogen following resection, and stored at −80 °C. Later, the frozen tumor samples were processed, and assessed based on cell morphology following hematoxylin and eosin-staining. Samples containing more than 20% of normal cells were excluded from RNA extraction to reduce the effects of contamination. Total RNA was isolated using MirVana miRNA Isolation kit (Thermo Fisher Scientific, Waltham, USA). Subsequently, quantification of exacted total RNA was performed with a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, USA).

The microarray data analysis was performed using an Agilent Whole Human Genome Array (Agilent, California, USA). First, an Agilent 2100 Bioanalyzer was used to determine total RNA integrity, then cDNA and biotinylated cRNA were isolated from the total RNA and hybridized onto microarrays. Next, the Agilent G2565BA Microarray Scanner System and Agilent Feature Extraction Software were used for data acquisition. The normalization of probe intensity was performed with Gene Spring GX 11.0 software.

Statistical analysis

Statistical analyses were performed using the IBM SPSS software, version 19.0. Continuous data were expressed as mean ± standard deviation (SD), and a p value less than 0.05 was considered statistically significant. First, the mRNA expression of five HLA-DR genes in gliomas of each grade and normal brain tissue were compared using one-way analysis of variance (ANOVA). Next, to evaluate the consistency of HLA-DR gene expression trends, a Pearson correlation matrix was created, and the Kaiser-Meyer-Olkin (KMO, criterion: > 0.5) and Bartlett’s test for sphericity (criterion: p < 0.05) were used to test the adequacy of the correlation matrix. To eliminate the interference of multicollinearity, a second factor analysis was performed to define the new variable, which can represent the expression characteristics of the entire HLA-DR gene group. Patients with different tumor grades were classified into two subgroups for further analysis, based on their values of the defined variable (cut off at 50% of the entire group). Univariate analysis was carried out with the chi-square test for dichotomous variables. The log-rank test and Breslow test were both used to explore correlations of the defined variable with progression-free survival (PFS) and overall survival (OS) in each glioma grade. PFS was defined as time from surgical resection until tumor recurrence (diagnosed by magnetic resonance (MR) imaging) or date of last follow-up, while OS was defined as the time from surgical resection until death or date of last follow-up. Finally, a Cox proportional hazards model was used to generate the multivariate model that described the correlation of prognostic factors with clinical outcomes.

Results

Patient characteristics

Among the 305 glioma patients enrolled in the current study, 126 were characterized as LGGs, 51 as anaplastic gliomas (AGs), and 128 as GBMs, according to the 2007 WHO classification. Table 1 shows the baseline characteristics of the cohort.

Table 1 Baseline characteristics of glioma patients (n = 305)

Expression of five HLA-DR genes in gliomas

First, the mRNA expression levels of HLA-DR genes between normal brain tissue and gliomas of each grade were assessed using one-way ANOVA. Significant differences were observed in the mRNA expression of all five HLA-DR genes among normal brain tissues and different grades of gliomas (Fig. 1, p < 0.001 for all five HLA-DR genes, One-way ANOVA). Moreover, overexpression of those genes was associated with higher tumor grade (Fig. 1, normal brain tissue < LGGs < AGs < GBMs, mean). The expression of all five HLA-DR genes in normal brain tissue have been summarized in Additional file 1: Table S1.

Fig. 1

Mean mRNA expression values of a HLA-DRA; b HLA-DRB1; c HLA-DRB3; d HLA-DRB4; e HLA-DRB5 in all 310 samples of normal brain tissue and gliomas (5 normal brain tissues, 126 low-grade gliomas, 51 anaplastic gliomas and 128 glioblastomas). Statistics were performed with one-way ANOVA analysis

Factor analysis and HLA-DR score calculation

Next, the consistency of HLA-DR gene expression in all 310 samples was confirmed by Pearson correlation analysis and the correlation matrix shown in Table 2. This correlation matrix in combination with a KMO index of 0.867 and a p-value less than 0.001 in the Bartlett test was appropriate for factor analysis, which was performed to avoid multicollinearity issues in regression analysis. The results are summarized in Table 3, and the Gravel Map (Fig. 2) indicated that the cumulative contribution rate of the first principal component eigenvalue is 94.952% for the total variance.

Table 2 Correlation matrix of HLA-DR genes

Table 3 Total variance explained^a

Fig. 2

Gravel Map of Factor analysis

Under normal circumstances, when the cumulative contribution eigenvalue of the first k principal component is over 85%, there is enough information of all the original factors. Therefore, the information provided by the first component in our analysis was equivalent and representative of the entire HLA-DR gene group, and we defined it as the HLA-DR score. According to the component score coefficient matrix, the factor score of each HLA-DR gene was calculated (Table 4). Finally, the detailed method to determine the HLA-DR score was as follows:

$$ {\displaystyle \begin{array}{l}\mathrm{HLA}\hbox{-} \mathrm{DR}\ \mathrm{score}=\left(0.445\times \mathrm{HLA}\hbox{-} \mathrm{DR}\mathrm{A}\right)+\left(0.455\times \mathrm{HLA}\hbox{-} \mathrm{DR}\mathrm{B}1\right)\\ {}\kern13em +\left(0.440\times \mathrm{HLA}\hbox{-} \mathrm{DR}\mathrm{B}3\right)+\left(0.439\times \mathrm{HLA}\hbox{-} \mathrm{DR}\mathrm{B}4\right)+\left(0.456\times \mathrm{HLA}\hbox{-} \mathrm{DR}\mathrm{B}5\right).\end{array}} $$

Table 4 Rotated component matrix and component score coefficient matrix for five HLA-DR genes

Analysis of HLA-DR score in different grades of gliomas

The HLA-DR scores were compared in gliomas from different grades and normal brain tissue (Fig. 3). There was a significant difference (p < 0.001, One-way ANOVA) in HLA-DR scores between test groups, and higher HLA-DR scores were associated with increased tumor grade. Next, patients with different tumor grades were further divided into two subgroups according to their HLA-DR score. A chi-square test was subsequently used to analyze the associations between the HLA-DR scores and clinicopathological features (Table 5). A significant association was observed between the HLA-DR score and histological subtype of LGG (p < 0.001, chi-square test). Subsequently, a Student’s t-test was used to further analyze these associations, and we observed that in LGGs, the HLA-DR score in astrocytic tumors was significantly higher compared to oligodendroglial tumors (−0.718 ± 3.177 versus −2.975 ± 2.662, p < 0.001, t-test; Fig. 4).

Fig. 3

HLA-DR scores in all 310 samples of normal brain tissue and gliomas (5 normal brain tissues, 126 low-grade gliomas, 51 anaplastic gliomas and 128 glioblastomas). Statistics were performed with one-way ANOVA analysis

Table 5 Association between HLA-DR score and clinical characteristics^a

Fig. 4

Comparison of HLA-DR scores between low-grade glioma (LGG) tumor samples with different histological subtypes. The astrocytic tumors displayed an average HLA-DR score of −0.718 ± 3.177, while oligodendroglial tumors displayed an average of −2.975 ± 2.662. The difference was significant between the two subtypes (p < 0.001, t-test)

Correlation of HLA-DR score with glioma patient survival

Finally, the correlations of the HLA-DR score with clinical outcomes were examined using Kaplan-Meier analysis. The mean follow-up time was 42.8 ± 35.5 months (median, 27.4 months; range, 0.7-111.4 months) and 176 patients died (60.8% males) during the follow-up time. For PFS, significant differences (p = 0.046, log-rank test; p = 0.021, Breslow test; Fig. 5a) were found in LGG patients when high and low HLA-DR scores were compared. However, for OS, the Breslow test showed significance between the two subgroups (p = 0.029, Fig. 5b), while the log-rank test did not confirm this difference (p = 0.053, Fig. 5b). Based on this, we further analyzed the correlations of the HLA-DR score with clinical outcomes in astrocytic and oligodendroglial tumors. For patients with oligodendroglial tumors, a low HLA-DR score was associated with significantly improved PFS (p = 0.023, log-rank test; p = 0.024, Breslow test; Fig. 5c) and OS (p = 0.020, log-rank test; p = 0.021, Breslow test; Fig. 5d). For patients with astrocytic tumors, a low HLA-DR score also seemed to indicate better survival, but these data were not statistically significant (Fig. 5e and f). Overall, LGG patients with low HLA-DR scores appeared to have a better clinical outcome. In contrast, no significant difference was found in either PFS or OS between AG and GBM patients with high and low HLA-DR scores (Additional file 1: Figure S1).

Fig. 5

Kaplan-Meier survival analysis. a Comparison of the PFS between high and low HLA-DR score groups in patients with LGGs (p = 0.046, log-rank test; p = 0.021, Breslow test). b Comparison of the OS between the two groups in patients with LGGs (p = 0.053, log-rank test; p = 0.029, Breslow test). c Comparison of the PFS between the two groups in patients with LGG oligodendroglial tumors (p = 0.023, log-rank test; p = 0.024, Breslow test). d Comparison of the OS between the two groups in patients with LGG oligodendroglial tumors (p = 0.020, log-rank test; p = 0.021, Breslow test). e Comparison of the PFS between the two groups in patients with LGG astrocytic tumors (p = 0.377, log-rank test; p = 0.078, Breslow test). f Comparison of the OS between the two groups in patients with LGG astrocytic tumors (p = 0.205, log-rank test; p = 0.092, Breslow test). PFS Progression-free survival; OS overall survival; LGG low-grade glioma

Furthermore, we also performed multivariate progression analysis to evaluate the independent value of the HLA-DR score and other variables for predicting PFS and OS in patients with LGGs. The HLA-DR score was enrolled into the Cox model as a continuous variable, and we observed the HLA-DR score to be an independent predictor of patient survival for LGGs (p = 0.042 and p = 0.025 for PFS and OS, respectively; Table 6).

Table 6 Multivariate predictors of PFS and OS for patients with LGGs^a

Discussion

HLA-DR is the most abundant form of MHC class II antigens and it can inhibit tumor growth through two possible mechanisms. In the first immunological mechanism, HLA-DR contributes to the recognition of tumor-associated antigens by CD4⁺ T cells, which subsequently produce specific cytokines, such as interleukins and interferon-γ, which eventually results in tumor suppression [23]. The second mechanism is non-immunological, where HLA-DR is directly involved in the regulation of antitumor action by interferon-γ [24]. Overexpression of HLA-DR has been demonstrated in a variety of human malignancies including thyroid carcinoma, gastric carcinoma, colorectal cancer and cervical cancer [25,26,27,28]. However, in other tumors like breast cancer and ovarian cancer, HLA-DR was found to be downregulated, and correlated with immune response evasion and tumor aggressiveness [29,30,31]. These results suggested that the regulatory mechanisms of HLA-DR differed among different types of cancers.

It has been indicated by neuro-immunologists that HLA-DR antigens are expressed on glioma cells [32]. One study has reported that MHC class II antigens, mainly HLA-DR and DQ, are important for lymphocyte response, and transfection of these genes could increase the immunogenicity of glioma cells [33]. Consistently, we also observed higher mRNA expression of five HLA-DR genes that were linked with increased tumor grade. In addition, the HLA-DR score, which defined the expression of all HLA-DR genes, showed that a higher value was correlated with an increased tumor grade. Interestingly, we also observed that astrocytic tumors, which are the more aggressive histological subtype of LGGs, had a significantly higher HLA-DR score in comparison to oligodendroglial tumors.

The relationship of HLA-DR expression with patient survival in cancer was first described in large bowel carcinoma in 1993, where strong HLA-DR expression was correlated with good prognosis [34]. Since then, multiple studies in many cancers (colorectal, gastric and breast) showed a similar pattern, where tumors that were HLA-DR positive had a better prognosis [26, 28, 35]. However, our data in glioma patients was conflicting with these observations, as the glioma patients with high HLA-DR scores, i.e. high mRNA expression of HLA-DR genes, had a poor clinical outcome. Furthermore, the HLA-DR score was demonstrated to be an independent predictor for survival of LGG patients. A previous study that analyzed the protein levels of HLA-DR in glioma patients was consistent with our result [20].

Here seemed to be a contradiction. Theoretically, tumors with higher HLA-DR expression should show stronger immunogenicity, thus appear to be less aggressive and have a better prognosis. However, according to the current study, gliomas with high HLA-DR scores were precisely those with poor immunogenicity, few characterized cancer antigens, and more invasive characteristics [36], and in LGGs, patients with high HLA-DR scores had a poor survival rates. This contradictory observation regarding the correlation of a high HLA-DR score with poor clinical outcome can be attributed to the complex mechanisms of glioma immunogenicity, which have only been partially elucidated. Despite HLA-DR playing a role in shaping tumor immunogenicity, we speculate that the positive effects have been offset by much stronger negative effects. For instance, at least two principle HLA-related mechanisms of immune escape have been proposed. One is associated with overexpression of HLA-G, which we have previously demonstrated and has been found to contribute to GBM immune escape [18, 37]. The second is the complete loss or downregulation of classical HLA class I molecules, which we need to investigate further, but has been reported in other cancers [38]. Moreover, we cannot rule out the possibility of other specific yet unknown mechanisms involved in enhancing/decreasing the immunogenicity of gliomas.

In recent years, immunotherapy has become one of the major options for anti-cancer therapy. Based on the results of our survival analysis, we may initially extrapolate to determine the type of glioma patient suitable to receive immunotherapy. Our survival analysis results suggested no significant difference in the prognosis of subgroups with high and low HLA-DR scores in AGs or GBMs patients. This can be explained by the fact that the interaction between the immune system and tumor is a relatively slow process, and the rapid progression of these high-grade gliomas did not allow the immune system enough time to exert its influence. Therefore, the efficacy of immunotherapy may be limited in patients with high-grade gliomas. In addition, it has been reported that several months of Bacille Calmette-Guerin (BCG) immunotherapy, could induce a strong level of HLA-DR expression in bladder cancer [39]. Another study showed that BCG-treated patients with absent or low expression of HLA-DR had a prolonged disease-free survival in melanoma [40]. Accordingly, immunotherapy may show better efficacy on LGG patients with low HLA-DR expression.

Conclusion

In the current study, data from 305 patients having histologically confirmed gliomas were retrospectively analyzed, and we observed that the mRNA expression levels of five HLA-DR genes were higher in different grades of gliomas compared to normal brain tissue. In addition, we defined a new variable, termed the HLA-DR score, to represent the expression of the entire HLA-DR gene group. We observed that a higher HLA-DR score was correlated with both the higher tumor grade in gliomas and the more aggressive histological type in LGG. Furthermore, the HLA-DR score was also an independent predictive factor of clinical outcomes in LGG patients. Considering the complex interaction between the HLA system and gliomas, these results are just the tip of the iceberg. Further investigation into the roles of other HLA molecules in gliomas are still required, and these studies may further reveal the interaction between gliomas and the immune system.