Correlation of alteration of HLA-F expression and clinical characterization in 593 brain glioma samples
Journal of Neuroinflammation volume 16, Article number: 33 (2019)
Human gliomas are highly fatal tumors with a significant feature of immune suppression. The association of the immune system in gliomas is gradually revealed, and immunotherapy is expected to improve the survival of glioma patients. In-depth understanding of the immune microenvironment of gliomas and their associated immunotherapy was increased exponentially in recent years. Gliomas provide clinical targets for immunotherapy during the search of key regulators of immune response. Our study focused on the human leukocyte antigen (HLA) system that is responsible for regulating the immune system, and discovered the relationship between HLA-F expression and clinical prognosis in gliomas.
A total of 593 patients with gliomas were included in our research. Of these, 325 patients were from the Chinese Glioma Genome Atlas (CGGA) and 268 were from the GSE 16011 set. Kaplan-Meier (KM) analysis was performed to explore the prognostic value of HLA-F. t test analysis was used to find the distribution difference in various groups. R language packages are used for other statistical computations and figure drawing.
HLA-F was negatively correlated with overall survival (OS) in all grades of glioma and glioblastoma (GBM). Moreover, HLA-F was enriched in GBM and isocitrate dehydrogenase 1 wild-type (IDH1 wt) group and considered HLA-F as a mesenchymal subtype marker. Pearson correlation test showed that HLA-F was correlated with other HLA-I molecules.
HLA-F expression was positively correlated with malignant phenotype and negatively correlated with OS, indicating that HLA-F could predict the immune state of gliomas and might be a clinical target of glioma immunotherapy.
Glioma is the most common and malignant brain tumor, accounting for 46% of all intracranial tumors and 2% of all adult cancers [1, 2]. Based on the 2007 World Health Organization (WHO) classification of gliomas, these are classified from WHO grade I to WHO grade IV . The invasion of tumor cells is gradually increased with the progression of tumor grade . Despite the advanced therapeutic methods, the survival time of gliomas show no significant improvement from the past four decades; especially for glioblastoma (GBM), the median survival time is nearly 14 months [5, 6]. The unfavorable prognosis is largely attributed to the escape of immune surveillance and suppression of the immune system . Hence, it is necessary to uncover the regulators of the immune system in gliomas to improve the patients’ survival time.
Although the central nervous system (CNS) is originally considered to be an immune-privileged site, more recent studies discover that it is an immunologically specialized site [8, 9]. Immunotherapy by reversing the suppressed immune status in tumors becomes a promising anticancer therapy . Regulators of immune escape are prevented effectively by current immunotherapies to eliminate gliomas . These defects in surveillance of gliomas could be due to the range of alterations of human leukocyte antigen (HLA) class I alleles [12, 13], dysfunction of T cell receptor (TCR) [14, 15], and abnormal expression of co-stimulatory ligands on glioma cells.
HLA-F is located at the terminal end of chromosome 6, is one of the non-classical HLA class Ib molecules, and is considered as a tolerogenic molecule . HLA-F is a ligand for KIR receptors in natural killer cells (NK cells) and CD8+ T cells. HLA-F interacts physically and functionally with three KIR receptors: KIR3DL2, KIR2DS4, and KIR3DS1. HLA-F regulates immunity during infection , pregnancy , and autoimmunity  through NK cell receptors.
HLA-F was reported to be involved in the immunomodulatory effect of visceral system tumors. For example, Lin et al. reported that HLA-F expression was negatively associated with the overall survival (OS) in non-small cell lung cancer patients . Ishigami et al. showed that HLA-F expression was upregulated in gastric cancer . But the function of HLA-F in gliomas remains unknown. In our research, we used the Chinese Glioma Genome Atlas (CGGA) RNA-seq set as a training set and validated the results in the GSE 16011 array set. We hypothesiz that HLA-F plays as an oncogene in glioma and is suggested as a new potential biotarget for therapy.
Patients and data collection
HLA-F expression is downloaded from CGGA which includes 109 grade II, 72 grade III, and 144 grade IV and from the GSE 16011 array set which includes 24 grade II, 85 grade III, and 159 grade IV. To further validate the protein expression of HLA-F in different grades, we randomly selected 15 different grade glioma patients (5 patients for each grade) from Beijing Ditan Hospital for immunohistochemistry (IHC). All the patients are histologically diagnosed by two neuropathologists according to the 2007 WHO classification guidelines, and clinical information is downloaded from each website. This research is approved by the Ethics Committee of Beijing Ditan Hospital.
Detection of biomarkers of gliomas
Isocitrate dehydrogenase 1 (IDH1) mutations are the most common mutations in gliomas, especially in low-grade gliomas and secondary GBM . In the CGGA RNA-seq set, the genomic DNA was isolated from frozen tissues using a QIAamp DNA Mini Kit (Qiagen) according to the manufacturer’s protocol . IDH1 status in the CGGA RNA-seq set was tested by pyrosequencing. TCGA subtype was based on previous research .
Immunohistochemistry for HLA-F expression
IHC was performed to detect HLA-F protein expression. The tumor tissues excised during the operation were immediately placed in 10% formalin for fixation, followed by dehydration, paraffin embedding, and sectioning. Anti-HLA-F antibody (Abcam) was used at a dilution of 1:100. Each slide was individually reviewed and scored by two independent neuropathologists, HLA-F protein expression (semi-quantitative scoring) = expression intensity × expression area. Expression intensity was scored using a 4-point scale from 0 to 3. Expression area was scored using a 5-point scale from 0 to 4.
Related signature identification
Significantly related genes with HLA-F expression were retrieved by using the Pearson correlation analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyzed the related gene sets and relevant biological function on the DAVID website (http://david.ncifcrf.gov/). In addition, the heatmap package of R language was used to list the genes positively related with HLA-F expression.
We performed Kaplan-Meier analysis to explore the prognostic value of HLA-F. Cox proportional hazards model analysis was performed to verify HLA-F as an independent prognostic factor. OS was estimated from the date of diagnosis to the date of either death or last follow-up. Student’s t test was performed to explore the expression distribution in different groups. Pearson correlation analysis was used to find genes related to HLA-F expression. R language packages (ggplot2, pheatmap, pROC, and corrgram) are used for other statistical computations and drawing figures. All differences were considered statistically significant at the level of p < 0.05.
HLA-F mRNA expression is upregulated in high-grade gliomas and downregulated in IDH1 mutation gliomas
Firstly, 325 samples with HLA-F expression profile and whole follow-up information are obtained from the CGGA RNA-seq set. The GSE 16011 array set was used for validation. In Fig. 1a, the results showed that HLA-F mRNA expression was positively associated with grade (p < 0.05). Moreover, HLA-F expression was higher in the isocitrate dehydrogenase 1 wild-type (IDH1 wt) group than in the IDH1 mutation (IDH1 mut) group (p < 0.05, Fig. 1b). Fortunately, all the above results were validated in 268 samples from the GSE 16011 array set (Fig. 1c, d). To further explore the relation between tumor grade and HLA-F protein expression, immunohistochemistry (IHC) was performed (Fig. 2a, b), and the results showed that the protein expression of HLA-F (semi-quantitative scoring) was higher in high-grade gliomas (p < 0.05).
HLA-F is a potential marker for mesenchymal molecular subtype gliomas
To explore the expression distribution of HLA-F, we evaluated the expression pattern of HLA-F in different molecular subtypes as defined by the Cancer Genome Atlas (TCGA) network . The results (Fig. 3a, b) showed that HLA-F was significantly upregulated in the mesenchymal molecular subtype than in other subtypes in both the CGGA RNA-seq set and GSE 16011 array set. To further validate this finding, receiver operating characteristic curves (ROC) for HLA-F expression and mesenchymal subtype of all grade gliomas are performed. Interestingly, the areas under the curves (AUC) are 80.5% and 76.5% in the CGGA RNA-seq set and GSE 16011 array set, respectively (Fig. 3c, d).
HLA-F is an independent predictive marker of OS in patients with malignant gliomas
All 593 patients from the CGGA RNA-seq set and GSE 16011 array set have HLA-F mRNA expression information, and the median value was used as the cutoff point (0.2 for CGGA, 10.8 for GSE 16011). Then, we evaluated the survival time of each group and found that the patients with higher HLA-F expression have remarkably shorter OS in all grade gliomas and GBM in the CGGA RNA-seq set (Fig. 4a, b; p < 0.0001 for all grade gliomas, p = 0.0587 for GBM). Similar results were obtained from the GSE 16011 array set (Fig. 4c, d; p < 0.0001 for all grade gliomas, p = 0.0238 for GBM). Then, to verify HLA-F expression as an independent prognostic factor (Table 1), univariate and multivariate Cox regression analyses were performed in CGGA RNA-seq. The results of univariate regression showed that HLA-F (p < 0.001) along with grade (p < 0.0001), age (p < 0.0001), and IDH1 status (p < 0.0001) predicted the OS in all grade gliomas. In multivariate regression, HLA-F showed significant results (p = 0.027) after adjusting for grade (p < 0.0001), age (p = 0.687), and IDH1 status (p = 0.023). The above results suggested that HLA-F plays an important role in the progression of gliomas. Next, biologic function analysis should be performed to further validate our findings.
HLA-F is associated with immune functions in gliomas
To investigate the relationship of HLA-F and another HLA-I molecular family, we performed Pearson correlation analysis in the CGGA RNA-seq set and GSE 16011 array set. The results (Fig. 5a, b) showed that the HLA-F expression was significantly associated with HLA-A, HLA-B, HLA-C, and HLA-E. To explore the biological processes associated with HLA-F expression in gliomas, Pearson correlation analysis between HLA-F expression and other genes in whole genome gene profiling of 325 patients in the CGGA RNA-seq set was performed. The results revealed that 183 genes (R > 0.5) are positively related to HLA-F expression (Fig. 6a). Biological process (BP) and KEGG analysis are also performed through the DAVID website (https://david.ncifcrf.gov). The results are listed in Fig. 6b, wherein 5 BPs for 183 genes were observed including antigen processing and presentation, defense response to virus, innate immune response, negative regulation of viral genome replication, and positive regulation of catalytic activity. As for KEGG analysis, the proteasome, Epstein-Barr virus infection, antigen processing and presentation, insulin signaling pathway, and phagosome were included. In summary, all the above results indicated that HLA-F could affect glioma-related immune activities.
Gliomas are the most common type of primary malignant tumors and possess high cancer-associated mortality rate . According to the Stupp protocol, surgery and postoperative radiotherapy plus chemotherapy have improved the survival of patients with gliomas; however, the outcome of most of the glioma patients remained poor . Slow progression in the extension of patients’ lives highlighted the urgent need for new and effective treatments for gliomas .
Immunotherapies for gliomas were considered to have a bright future. Anti-glioma immunotherapy is a general term that encompassed various strategies that are intended to stimulate the patient’s innate and adaptive immune system against gliomas and promote immune-mediated anti-glioma responses . The TCGA network describes a robust gene expression-based molecular classification of gliomas into proneural, neural, classical, and mesenchymal subtypes. Gliomas in the mesenchymal group are more malignant than others . In this study, we retrospectively analyzed 593 glioma patients from the CGGA RNA-seq set and GSE 16011 array set. Our findings indicated that HLA-F significantly predicted the OS in all grade gliomas and GBM. Moreover, the mRNA expression level remarkably varied in different grades and IDH1 status groups. Next, we found that HLA-F was upregulated in the mesenchymal group, enlightening HLA-F as a biomarker for the mesenchymal subtype, a more malignant phenotype and worse prognosis.
HLA complex, also known as the major histocompatibility complex (MHC), consists of three groups: HLA-I (e.g., HLA-A, HLA-E, and HLA-F), HLA-II (e.g., HLA-DR, HLA-DP, and HLA-DQ), and HLA-III . Former studies indicated that HLA-F expression could protect the trophoblasts from NK- and CTL-mediated lysis in the maternal uterine environment . Moreover, HLA-F, as an important immunosuppressive molecule, is exhibited differently in various types of tumors. In hepatocellular carcinoma, HLA-F expression was significantly correlated with the degree of lymphatic or venous invasion . Furthermore, in esophageal squamous cell carcinoma, higher HLA-F expression indicated worse survival .
HLA-A, HLA-B, and HLA-C mainly present antigens to CD8 T cells and regulate various immunologic functions. HLA-E and HLA-F mainly regulate NK cell function through NK cell receptors. HLA-G is highly expressed on the surface of cytotrophoblasts and placenta and presented self-polypeptides. We found that HLA-F was correlated with most of the HLA complexes, such as HLA-A, HLA-B, HLA-C, and HLA-E. A variety of HLA molecules were enriched in gliomas to participate in the immunosuppressive state of CD8 and NK cells. Furthermore, HLA-F-positive gliomas were associated with immune response, indicating that the function of HLA-F in gliomas involved the regulation of immune status of tumor tissues.
Lack of HLA-F could render gliomas to be susceptible to elimination by NK cells. In our opinion, during the development of gliomas, glioma-specific and microenvironment cue-induced HLA-F expression could be modulated on glioma cells. Higher HLA-F expression in glioma cell surface provided immunoinhibitory effects and protected malignant glioma cells from elimination by NK cells, inducing the growth of glioma cells. HLA-F may act as a novel target of glioma immunotherapy.
In conclusion, to the best of our knowledge, this was the first study to investigate the HLA-F expression levels in gliomas. The results demonstrated that HLA-F could be considered as a novel immunotherapy target in gliomas.
Area under the curve
Chinese Glioma Genome Atlas
Central nervous system
Human leukocyte antigen
Isocitrate dehydrogenase 1
Kyoto Encyclopedia of Genes and Genomes
- NK cells:
Natural killer cells
World Health Organization
Hu H, Wang Z, Li M, et al. Gene expression and methylation analyses suggest DCTD as a prognostic factor in malignant glioma. Sci Rep. 2017;7:11568.
Sun X, Ma X, Wang J, et al. Glioma stem cells-derived exosomes promote the angiogenic ability of endothelial cells through miR-21/VEGF signal. Oncotarget. 2017;8:36137–48.
Fan X, Wang Y, Zhang C, et al. ADAM9 expression is associate with glioma tumor grade and histological type, and acts as a prognostic factor in lower-grade gliomas. Int J Mol Sci. 2016;17:1276.
Liu Y, Hu H, Wang K, et al. Multidimensional analysis of gene expression reveals TGFB1I1-induced EMT contributes to malignant progression of astrocytomas. Oncotarget. 2014;5:12593–606.
Chen B, Liang T, Yang P, et al. Classifying lower grade glioma cases according to whole genome gene expression. Oncotarget. 2016;7:74031–42.
Huang K, Fang C, Yi K, et al. The role of PTRF/Cavin1 as a biomarker in both glioma and serum exosomes. Theranostics. 2018;8:1540–57.
Fan X, Wang Y, Zhang C, Liu X, Qian Z, Jiang T. Human leukocyte antigen-G overexpression predicts poor clinical outcomes in low-grade gliomas. J Neuroimmunol. 2016;294:27–31.
Li X, Wang B, Gu L, et al. Tim-3 expression predicts the abnormal innate immune status and poor prognosis of glioma patients. Clin Chim Acta. 2018;476:178–84.
Louveau A, Harris TH, Kipnis J. Revisiting the mechanisms of CNS immune privilege. Trends Immunol. 2015;36:569–77.
Huang B, Zhang H, Gu L, et al. Advances in immunotherapy for glioblastoma multiforme. J Immunol Res. 2017;2017:3597613.
Wei J, Nduom EK, Kong LY, et al. MiR-138 exerts anti-glioma efficacy by targeting immune checkpoints. Neuro-Oncology. 2016;18:639–48.
Wiendl H, Mitsdoerffer M, Hofmeister V, et al. A functional role of HLA-G expression in human gliomas: an alternative strategy of immune escape. J Immunol. 2002;168:4772–80.
Wischhusen J, Friese MA, Mittelbronn M, Meyermann R, Weller M. HLA-E protects glioma cells from NKG2D-mediated immune responses in vitro: implications for immune escape in vivo. J Neuropathol Exp Neurol. 2005;64:523–8.
Li G, Wang Z, Zhang C, et al. Molecular and clinical characterization of TIM-3 in glioma through 1,024 samples. Oncoimmunology. 2017;6:e1328339.
Liu S, Wang Z, Wang Y, et al. PD-1 related transcriptome profile and clinical outcome in diffuse gliomas. Oncoimmunology. 2018;7:e1382792.
Laaribi AB, Hannachi N, Ben Yahia H, et al. Human leukocyte antigen (HLA-F) polymorphism is associated with chronic HBV infection. 3 Biotech. 2018;8:49.
Garcia-Beltran WF, Holzemer A, Martrus G, et al. Open conformers of HLA-F are high-affinity ligands of the activating NK-cell receptor KIR3DS1. Nat Immunol. 2016;17:1067–74.
Burrows CK, Kosova G, Herman C, et al. Expression quantitative trait locus mapping studies in mid-secretory phase endometrial cells identifies HLA-F and TAP2 as fecundability-associated genes. PLoS Genet. 2016;12:e1005858.
Song S, Miranda CJ, Braun L, et al. Major histocompatibility complex class I molecules protect motor neurons from astrocyte-induced toxicity in amyotrophic lateral sclerosis. Nat Med. 2016;22:397–403.
Lin A, Zhang X, Ruan YY, Wang Q, Zhou WJ, Yan WH. HLA-F expression is a prognostic factor in patients with non-small-cell lung cancer. Lung Cancer. 2011;74:504–9.
Ishigami S, Arigami T, Okumura H, et al. Human leukocyte antigen (HLA)-E and HLA-F expression in gastric cancer. Anticancer Res. 2015;35:2279–85.
Wang HY, Tang K, Liang TY, et al. The comparison of clinical and biological characteristics between IDH1 and IDH2 mutations in gliomas. J Exp Clin Cancer Res. 2016;35:86.
Verhaak RG, Hoadley KA, Purdom E, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. 2010;17:98–110.
Liu Z, Han H, He X, et al. Expression of the galectin-9-Tim-3 pathway in glioma tissues is associated with the clinical manifestations of glioma. Oncol Lett. 2016;11:1829–34.
Stupp R, Mason WP, van den Bent MJ, et al. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl J Med. 2005;352:987–96.
Lowenstein PR, Baker GJ, Castro MG. Cracking the glioma-NK inhibitory code: toward successful innate immunotherapy. Oncoimmunology. 2014;3:e965573.
Boussiotis VA, Charest A. Immunotherapies for malignant glioma. Oncogene. 2018;37:1121–41.
Nakajima M, Takahashi A, Kou I, et al. New sequence variants in HLA class II/III region associated with susceptibility to knee osteoarthritis identified by genome-wide association study. PLoS One. 2010;5:e9723.
Szekeres-Bartho J. Immunological relationship between the mother and the fetus. Int Rev Immunol. 2002;21:471–95.
Xu Y, Han H, Zhang F, Lv S, Li Z, Fang Z. Lesion human leukocyte antigen-F expression is associated with a poor prognosis in patients with hepatocellular carcinoma. Oncol Lett. 2015;9:300–4.
Zhang X, Lin A, Zhang JG, et al. Alteration of HLA-F and HLA I antigen expression in the tumor is associated with survival in patients with esophageal squamous cell carcinoma. Int J Cancer. 2013;132:82–9.
We thank the CGGA and GSE 16011 network for sharing large amounts of data. We thank Beijing Ditan Hospital Specimen Bank Project for specimen storage.
This work is supported by grant from Beijing Medical Management Bureau Cultivation Plan (PX2018079) and the National Natural Science Foundation of China (NO. 81871013) and Beijing Tiantan Hospital Young Scientist Program (YSP 201705).
Availability of data and materials
Please contact author for data requests
Ethics approval and consent to participate
All the procedures in this study were approved by the ethics committees of all hospitals, and written informed consent was obtained from all patients.
Consent for publication
Written informed consent for publication of their clinical details and/or clinical images was obtained from the patient/parent/guardian/relative of the patient. A copy of the consent form is available for review by the Editor of this journal.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Feng, E., Liang, T., Wang, X. et al. Correlation of alteration of HLA-F expression and clinical characterization in 593 brain glioma samples. J Neuroinflammation 16, 33 (2019). https://doi.org/10.1186/s12974-019-1418-3