- Open Access
New biomarker: the gene HLA-DRA associated with low-grade glioma prognosis
Chinese Neurosurgical Journal volume 8, Article number: 12 (2022)
Low-grade gliomas (LGG) are WHO grade II tumors presenting as the most common primary malignant brain tumors in adults. Currently, LGG treatment involves either or a combination of surgery, radiation therapy, and chemotherapy. Despite the knowledge of constitutive genetic risk factors contributing to gliomas, the role of single genes as diagnostic and prognostic biomarkers is limited. The aim of the current study is to discover the predictive and prognostic genetic markers for LGG.
Transcriptome data and clinical data were obtained from The Cancer Genome Atlas (TCGA) database. We first performed the tumor microenvironment (TME) survival analysis using the Kaplan-Meier method. An analysis was undertaken to screen for differentially expressed genes. The function of these genes was studied by Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Following which a protein-protein interaction network (PPI) was constructed and visualized. Univariate and multivariate COX analyses were performed to obtain the probable prognostic genes. The key genes were selected by an intersection of core and prognostic genes. A clinical correlation analysis of single-gene expression was undertaken. GSEA enrichment analysis was performed to identify the function of key genes. Finally, a single gene-related correlation analysis was performed to identify the core immune cells involved in the development of LGG.
A total of 529 transcriptome data and 515 clinical samples were obtained from the TCGA. Immune cells and stromal cells were found to be significantly increased in the LGG microenvironment. The top five core genes intersected with the top 38 prognostically relevant genes and two key genes were identified. Our analysis revealed that a high expression of HLA-DRA was associated with a poor prognosis of LGG. Correlation analysis of immune cells showed that HLA-DRA expression level was related to immune infiltration, positively related to macrophage M1 phenotype, and negatively related to activation of NK cells.
HLA-DRA may be an independent prognostic indicator and an important biomarker for diagnosing and predicting survival in LGG patients. It may also be associated with the immune infiltration phenotype in LGG.
Glioma is a common primary intracranial neuroepithelial tumor occurring in the brain and arising in glial tissue. Low-grade gliomas are generally referred to as WHO grade 3. Routine histopathological classification of gliomas . However, the identification of key molecular alterations has led to substantial changes in the updated 2016 WHO Classification of Tumors of the CNS, such as IDH 1 and 2 mutations and 1p/19q codeletion [2, 3]. As such, the term “low-grade gliomas” is now often used to refer to both grade 2 and 3 gliomas, consistent with The Cancer Genome Atlas Project categorization. Molecularly altered gliomas are a heterogeneous group of primary brain tumors that vary not only in malignancy but also in histology and genomic alterations. They may arise from neural stem cells (NSCs), NSC-derived astrocytes, or oligodendrocyte precursor cells .
Low-grade gliomas are treated through surgery, chemotherapy, and/or radiation therapy. Despite the availability of therapeutic options, the possibility of recurrence exists . Some of the factors that may contribute to ineffective treatment are the inability in getting around the blood-brain barrier for efficient drug delivery, overcoming immune-suppressive tumor microenvironment (TME), and development of drug resistance . Significant advances have been made towards assessing the tumor microenvironment and have been used for treating LGG . Recently, several immunotherapy methods have been shown to be successful in treating malignant tumors. Some of these include immune checkpoint blockade, cytokine therapy, cell therapy, and therapeutic vaccines . But response rates vary between tumor types and within tumors. The variability in response limits the non-personalized use of immunotherapy for LGG . The success of immunotherapy is thus dependent on a better understanding of the glioma-specific immune microenvironment . While the genetic basis of gliomas is well-established, the utility of gene therapy in the field is less explored which can help overcome some of the therapeutic challenges and influence relapse rates in LGG. These gaps in the existing literature necessitate more research towards identifying new biomarkers and their applicability as therapeutic agents for LGG .
The major histocompatibility complex, class II, DR alpha (HLA-DRA) is a protein-coding gene. HLA-DRA is often associated with the occurrence of Graham-Little-Piccardi-Lassueur syndrome and penicillin allergy. The relevance of Parkinson’s disease with the HLA-DRA was also verified in a cohort of the Iranian population . In addition, Lee et al. verified an intergenic variant rs9268877 between HLA-DRA and HLA-DRB contributing to the clinical course and long-term outcome of ulcerative colitis . Chu et al. verified that HLA-DRA was a potential prognostic biomarker for renal clear cell carcinoma . However, HLA-DRA has not been reported in LGG.
In this study, first, we explored the differential expression of HLA-DRA in LGG tissues. We correlated the expression of HLA-DRA and immune correlation in the LGG dataset of The Cancer Genome Atlas (TCGA) to ascertain the significance of HLA-DRA expression within the tumor samples. In addition, we performed GO, KEGG, COX, HPA, GSEA, and meta-analysis to assess the overall prognostic significance of HLA-DRA using data from TCGA databases and CGGA databases. Further, as substantial attention has been focused on the crucial role of the immune microenvironment in the progression of LGG [9, 14, 15], we also evaluated the potential correlation between HLA-DRA and immune infiltration levels in LGG by applying tools from the R software to data from the TCGA database.
Access to TCGA datasets
The gene transcriptome data tissues involving 529 samples and clinical data involving 515 samples were downloaded from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/) . Clinical data collected included gender, age, grade, survival status, and duration of survival. Our study was in line with the publication guidelines provided by TCGA .
Survival analysis of microenvironment
We compiled the transcriptome data using Perl (https://www.perl.org/). Similarly, the clinical data was summarized. ESTIMATE was used to assess the immune cell infiltration levels (measured as immune score) and stromal content (measured as stromal score) per sample. The aforementioned analysis was undertaken using the Limma and the Estimate package of the R software. Furthermore, microenvironmental survival analysis was performed using the survival package of the R software. The samples were grouped into high and low grade based on their median score. Meanwhile, the corresponding p-value and the number of patients at each time point were determined. The survival curve was plotted using the R software with the following screening criterion: adj p < 0.001. To verify the correlation between tumor microenvironment and LGG, we conducted a clinical correlation analysis of the microenvironment; the boxplot was plotted using the R package ggpubr.
Based on the ESTIMATE analysis, we sub-grouped patients as high and low, based on their immune and stromal score, respectively . Differential expression analysis was performed between the high and low immune score groups and high and low stromal score groups by using the R package, with filter conditions of log2 |fold change| = 2 and false detection rate = 0.001. Heatmap and clustering were generated by using the “pheatmap” R package . The intersection of differential gene results was obtained, and the Venn diagram of the differential gene was plotted. Subsequently, the clinical correlation was analyzed between stromal and immune cell scores.
Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis
We used the “ClusterProfiler” R package for functional enrichment analysis of differential genes to identify the potential functions and pathways . Enrichment of genes was conducted by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). p-value < 0.05 and q-value < 0.05 were considered as statistically significant.
Protein–protein interaction (PPI) network
The PPI network was built using a search tool that retrieves a database of interacting genes (STRING) , and the minimum required interaction score was set to 0.95. We analyzed and visualized the network via Cytoscape software(https://cytoscape3.7.2.org/download.html) . Network characteristics such as the number of core proteins the and interplay between different proteins were studied.
COX regression model
Univariate and multivariate COX regression analysis and forest plots were plotted using the “forest plot” R package to display the p-values and hazard ratios (HRs) for each variable. Based on the results of COX analysis, prognostic-related genes were identified.
Key gene acquisition
Through the intersection of PPI core genes and prognostic genes, we obtained the most probable key prognostic genes and plotted the Venn diagram of differential genes using the “VennDiagram” R package. To verify the prognostic correlation, we performed a survival analysis of the key gene by the survival package of the R software. LGG patients were divided into two groups of high and low expression according to the median value of gene expression, and the survival correlation was verified by comparing the postoperative differences between the high and low gene expression groups. The survival plot was drawn using “ggplot.” To verify the clinical relevance, the samples were divided into two groups according to age, gender, and grade, and the differences in the key gene expression levels were compared. A p < 0.05 was considered as statistically significant.
Human Protein Atlas
The Human Protein Atlas (HPA, https://www.proteinatlas.org/) is a protein database that provides tissue and cellular distribution information for 26,000 human proteins . In this study, we used the HPA to verify the HLA-DRA gene expression in glioma cells and glial brain cells.
Gene set enrichment analysis
To explore the effects of the function or pathway of the gene expression in the tumor, samples were divided into two groups of high and low gene expression. Enrichment of the KEGG pathways in the high- and low-expression groups was further analyzed using the GSEA (ttp://software.broadinstitute.org/gsea/downloads.jsp) tool. The filter condition was set to 50, and the first 5 pathways were selected from the 50 enrichment pathways for multi-pathway enrichment analysis.
Meta-analysis of the included studies
The RNA expression matrix and clinical data of low-grade glioma came from the TCGA database (ID: TCGALGG) and CGGA database (http://www.cgga.org.cn/, Dataset ID: mRNAseq_693 and mRNAseq_325). All data were corrected using the “limma” package in R. Clinical correlation analysis and univariate COX analysis were performed by the R. We got three COX files. Then, a meta-analysis was conducted using the R software to verify the correlation between the HLA-DRA gene and the prognosis of LGG. Heterogeneity was tested by Q-statistic and I2 statistic, I2 > 50% was considered as a significant heterogeneity, the random-effects model or the fixed-effects model was adopted, and p < 0.05 indicated that there were statistical significances.
Immunological correlation analysis
The type and the relative number of immune cells in each tumor sample were visualized and plotted as a correlation curve by CIBERSORT ; the filter condition was set to a p-value = 0.05. Sample data with p < 0.05 were retained. Differential and correlation analyses were used to verify the correlation of immune cells and gene expression. The samples were divided into two groups of high and low gene expression. The presence of immune cells was comparatively analyzed between the two groups of high and low expression. The results obtained by the two analytical methods were intersected, to finally attain expressed immune cells in response to specific genes.
All analyses were performed using the R software and Perl. A p-value < 0.05 was considered statistically significant.
Tumor microenvironment with the prognosis of LGG
Transcriptome data (n = 529) and 515 clinical data extracted from the TCGA database were included in this study. An ESTIMATE analysis was performed. We found that the immune score of the selected samples was distributed between − 1681.2 and 2466.40, and the stromal score was distributed between − 1774.32 and 1701.59. By Kaplan-Meier analysis, we verified that the tumor microenvironment significantly correlated with survival in LGG patients (Fig. 1A). The clinical relevance of the microenvironment was also verified (Fig. 1B). As shown in the figure, the presence of stromal cells and immune cells in the tumor microenvironment is associated with the survival of the patient.
Screening for differential genes
Differential analyses on stromal cells were performed to screen out upregulated genes and downregulated genes. Based on the comparative scores, we found 200 genes that were upregulated and 6 genes that were downregulated. Similarly, we performed a differential analysis on immune cells and found 145 genes that were upregulated and 65 genes that were downregulated. Following which we plotted a heatmap using the fold change and corrected p-value (Fig. 2A, B). The intersection of upregulated genomes was considered to ascertain differentially expressed genes. We found 117 genes that were upregulated (Fig. 2C).
Functional enrichment analysis of differential genes
Upon undertaking a GO and KEGG enrichment analysis of the differentially expressed genes, we found that these genes were mainly responsible for T cell activation, leukocyte-mediated immunity, positive regulation of cytokine production, and cytokine−cytokine receptor interaction.
Screening of hub gene
We obtained a protein interactome for the differential genes using the String database. The network was visualized and analyzed using the Cytoscape software. The analysis yielded 30 core genes. The COX analysis, performed on each sample, reveals 38 prognostically relevant genes. We took the top five genes with the most nodes in the PPI and the 38 genes associated with prognosis to build the intersection. The Venn diagram yielded two key genes HLA-DRA and CD74 (Fig. 2D). Since high expression of CD74 in gliomas has been validated for the association with poor prognosis and high immune infiltrates . We proceeded with our analysis with HLA-DRA. Single-gene survival analysis showed p < 0.05, indicating that HLA-DRA gene expression was related to patient survival, and the high expression group predicted poor prognosis (Fig. 3A). Single-gene clinical correlation analysis showed that HLA-DRA expression was unrelated to age and gender. We found a positive correlation between the expression of HLA-DRA and LGG (Fig. 3B–D).
Human Protein Atlas
To verify the differences in the HLA-DRA gene expression in tissues, we used the HPA to analyze the HLA-DRA gene expression in the glial cells and brain. We found that the HLA-DRA gene was not expressed in normal glial cells (patient ID: 1371; 3731; 3739) but was highly expressed in glioma cells (patient ID: 3137; 3120; 3174). The results of immunohistochemical staining were shown in Fig. 4A, B.
Gene set enrichment analysis
The biological properties of HLA-DRA in LGG were analyzed by the GSEA enrichment method. The GSEA results showed that a high expression of HLA-DRA is mainly associated with cancer-associated cellular immunity whereas low expression is mainly associated with substance metabolism.
Meta-analysis of the included studies
Clinical correlation analysis and univariate COX analysis were performed by the R. We got three COX files. Then, three COX documents were meta-analyzed by the meta package of R. The meta-analysis results were shown in Fig. 5A.
Gene-related immune cells
CIBERSORT  (http://cibersort.stanford.edu/) is a gene expression-based deconvolution algorithm that evaluates the expression changes of a group of genes relative to all other genes in the sample. Thus, the number of immune cells in each sample can be accurately enumerated through this process. The continued performance of CIBERSORT has prompted reliance on its utility to cellular heterogeneity research [25, 26]. Immune cell visualization and immune cell correlation analysis were performed by calculating the immune cell content by CIBERSORT. We found a differential expression of immune cells between the groups showing high and low gene expression. The correlation analysis verified the association between immune cell content and gene expression. Based on the differential analysis method, we found 12 types of immune cells. Based on the correlation analysis method, we found 9 types of immune cells. Of which we found an overlap of seven different immune cells, namely, plasma cells, activated natural killer (NK) cells, monocytes, macrophages M1 and M2 polarized, eosinophils, and activated mast cells.
Significant progress has been made in the treatment of glioma with postoperative chemotherapy, radiotherapy, and/or immune intervention, but the prognosis for patients with glioma continues to remain dismal. Over the past 30 years, epidemiological studies have shown that the overall survival rate of patients with glioma has been poor , and this seemingly has changed very little with advances in recent clinical practices involving immunotherapeutic approaches [28,29,30]. A better understanding of the tumor immune microenvironment is critical to improving the efficacy of current immunotherapies. Gene therapies have also been evaluated in preclinical and clinical settings [31, 32]. While off-target effects are a major concern for gene therapy, localized management of drugs and recently more advanced gene-editing tools have significantly increased the convenience and specificity of gene editing . Therefore, there is still a need to accelerate research in the treatment of gliomas.
We used data from TCGA database and CGGA database and analyzed it using tools from the R language and Perl, to provide some verification for tumor microenvironment (TME) in low-grade glioma (LGG) and to identify key genes associated with prognostic immunity. First, we calculated the immune and stromal scores for each sample using the R software and found that the stromal cell and immune scores correlated significantly with LGG prognosis. The results showed that the stromal cells and immune cell scores of grade III gliomas were higher in comparison with grade I and II gliomas. This suggests that TME is associated with LGG tumor progression and disease prognosis.
Then, we performed a differential expression analysis on the high and low immune scoring groups and the stromal scoring group and identified 416 differentially expressed genes. Of which, 117 genes were found to be upregulated and six genes were found to be downregulated. The GO and KEGG analysis revealed that most of these genes were involved in the immune processes. The results show that immune cells and stromal cell in the LGG microenvironments are closely related which are concordant with previous findings [34, 35]. Through an R analysis, we obtained 123 related but differentially expressed genes. The protein–protein interactome (PPI) network explored the network relationship between survival genes and used Cytoscape for visualization. Significantly high network nodes, IL10, CCL2, CD74, and HLA-DRA, were found to be responsible for tumorigenesis, macrophage enrichment, and immune response among others [36,37,38]. Thirty-eight prognostic genes were obtained by the COX analysis, and then, the top five genes with the most nodes were selected from the PPI network. The intersection of the two genes was used to obtain the core gene HLA-RDA. Clinical correlation analysis of HLA-DRA showed that there was higher expression of HLA-DRA in grade III, which indicated that HLA-DRA could promote tumor development and serve as a reference for the utility of gene therapy in glioma. Then, we used the HPA to verify the HLA-DRA gene expression in glioma cells and glial brain cells. The results showed that HLA-DRA was not expressed in glioma cells but significantly expressed in glioma cells. This demonstrated the possibility of HLA-DRA as a new biomarker. The GSEA results showed that a high expression of HLA-DRA is mainly enriched in cancer-related cellular immunity whereas low expression is mainly enriched in substance metabolism. To better verify the reliability of the study, we adopted multiple databases for meta-analysis. The results show that HLA-DRA is highly uniform in the two databases. The results of immune cell correlation analysis showed that HLA-DRA was associated with a variety of immune cell infiltrates. This may indicate that HLA-DRA expression is correlated with multiple immune pathways or immune-related genes. At the same time, it provides a reference direction for future research.
From what has been discussed above, our results support the association of HLA-DRA with prognosis and immune infiltration of LGG. The mechanism of HLA-DRA in gliomas is not currently known, but our study validates the relevance of HLA-DRA to the tumor immune microenvironment. Further investigation into the association of HLA-DRA in the prognosis of LGG may allow to evaluate its utility as a potential therapeutic tool.
In this study, we identified the gene HLA-DRA as an independent prognostic marker for LGG based on bioinformatics analysis. The results showed that with the occurrence of LGG, the content of stromal cells and immune cells in the microenvironment increased. The high expression of HLA-DRA is associated with the poor prognosis of LGG, and the high expression of HLA-DRA is closely related to the immune cells in the microenvironment.
Availability of data and materials
The original contributions presented in this study are contained in the article/supplementary material. Further inquiries may be made to the appropriate authors.
Louis DN, Ohgaki H, Wiestler OD, Cavenee WK, Burger PC, Jouvet A, et al. The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol. 2007;114(5):547. https://doi.org/10.1007/s00401-007-0278-6.
van den Bent MJ, Weller M, Wen PY, Kros JM, Aldape K, Chang S. A clinical perspective on the 2016 WHO brain tumor classification and routine molecular diagnostics. Neuro Oncol. 2017;19(5):614–24. https://doi.org/10.1093/neuonc/now277.
Eckel-Passow JE, Lachance DH, Molinaro AM, Walsh KM, Decker PA, Sicotte H, et al. Glioma groups based on 1p/19q, IDH, and TERT promoter mutations in tumors. New Engl J Medicine. 2015;372(26):2499–508. https://doi.org/10.1056/NEJMoa1407279.
Zong H, Verhaak RG, Canoll P. The cellular origin for malignant glioma and prospects for clinical advancements. Expert Rev Mol Diagn. 2014;12(4):383–94. https://doi.org/10.1586/erm.12.30.
Aldape K, Brindle KM, Chesler L, Chopra R, Gajjar A, Gilbert MR, et al. Challenges to curing primary brain tumours. Nat Rev Clin Oncol. 2019;16(8):509–20. https://doi.org/10.1038/s41571-019-0177-5.
Banerjee K, Núñez FJ, Haase S, McClellan BL, Faisal SM, Carney SV, et al. Current approaches for glioma gene therapy and virotherapy. Front Mol Neurosci. 2021;14:621831. https://doi.org/10.3389/fnmol.2021.621831.
Xu Y, Zhang H, Sun Q, Geng R, Yuan F, Liu B, et al. Immunomodulatory effects of tryptophan metabolism in the glioma tumor microenvironment. Front Immunol. 2021;12:730289. https://doi.org/10.3389/fimmu.2021.730289.
Christofi T, Baritaki S, Falzone L, Libra M, Zaravinos A. Current perspectives in cancer immunotherapy. Cancers. 2019;11(10):1472. https://doi.org/10.3390/cancers11101472.
Wu F, Wang ZL, Wang KY, Li GZ, Chai RC, Liu YQ, et al. Classification of diffuse lower-grade glioma based on immunological profiling. Mol Oncol. 2020;14(9):2081–95. https://doi.org/10.1002/1878-0261.12707.
Tan Y, Zhang S, Xiao Q, Wang J, Zhao K, Liu W, et al. Prognostic significance of ARL9 and its methylation in low-grade glioma. Genomics. 2020;112(6):4808–16. https://doi.org/10.1016/j.ygeno.2020.08.035.
Jamshidi J, Movafagh A, Emamalizadeh B, Zare Bidoki A, Manafi A, Ghasemi Firouzabadi S, et al. HLA-DRA is associated with Parkinson’s disease in Iranian population. Int J Immunogenet. 2014;41(6):508–11. https://doi.org/10.1111/iji.12151.
Lee H-S, Yang SK, Hong M, Jung S, Kim BM, Moon JW, et al. An intergenic variant rs9268877 between HLA-DRA and HLA-DRB contributes to the clinical course and long-term outcome of ulcerative colitis. J Crohn’s Colitis. 2018;12(9):1113–21. https://doi.org/10.1093/ecco-jcc/jjy080.
Chu G, Jiao W, Yang X, Liang Y, Li Z, Niu H. C3, C3AR1, HLA-DRA , and HLA-E as potential prognostic biomarkers for renal clear cell carcinoma. Transl Androl Urology. 2020;9(6):2640656–2656. https://doi.org/10.21037/tau-20-699.
Wu F, Li GZ, Liu HJ, Zhao Z, Chai RC, Liu YQ, et al. Molecular subtyping reveals immune alterations in IDH wild-type lower-grade diffuse glioma. J Pathology. 2020;251(3):272–83. https://doi.org/10.1002/path.5468.
Deng X, Lin D, Zhang X, Shen X, Yang Z, Yang L, et al. Profiles of immune-related genes and immune cell infiltration in the tumor microenvironment of diffuse lower-grade gliomas. J Cell Physiol. 2020;235(10):7321–31. https://doi.org/10.1002/jcp.29633.
Blum A, Wang P, Zenklusen JC. SnapShot: TCGA-analyzed tumors. Cell. 2018;173(2):530. https://doi.org/10.1016/j.cell.2018.03.059.
Shi H, Song Y, Song Z, Huang C. CKMT1B is a potential prognostic biomarker and associated with immune infiltration in lower-grade glioma. Plos One. 2021;16(1):e0245524. https://doi.org/10.1371/journal.pone.0245524.
Yu G, Wang L-G, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics J Integr Biology. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(D1):D447–52. https://doi.org/10.1093/nar/gku1003.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Thul PJ, Åkesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, et al. A subcellular map of the human proteome. Science. 2017;356(6340). https://doi.org/10.1126/science.aal3321.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7. https://doi.org/10.1038/nmeth.3337.
Xu S, Li X, Tang L, Liu Z, Yang K, Cheng Q. CD74 correlated with malignancies and immune microenvironment in gliomas. Front Mol Biosci. 2021;8:706949. https://doi.org/10.3389/fmolb.2021.706949.
Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938–45. https://doi.org/10.1038/nm.3909.
Ali HR, Chlon L, Pharoah PDP, Markowetz F, Caldas C. Patterns of immune infiltration in breast cancer and their clinical implications: a gene-expression-based retrospective study. Plos Med. 2016;13(12):e1002194. https://doi.org/10.1371/journal.pmed.1002194.
Bense RD, et al. Relevance of tumor-infiltrating immune cell composition and functionality for disease outcome in breast cancer. J Natl Cancer I. 2016;109:djw192.
Claus EB, Walsh KM, Wiencke JK, Molinaro AM, Wiemels JL, Schildkraut JM, et al. Survival and low-grade glioma: the emergence of genetic information. Neurosurg Focus. 2015;38(1):E6. https://doi.org/10.3171/2014.10.FOCUS12367.
Weller M, Butowski N, Tran DD, Recht LD, Lim M, Hirte H, et al. Rindopepimut with temozolomide for patients with newly diagnosed, EGFRvIII-expressing glioblastoma (ACT IV): a randomised, double-blind, international phase 3 trial. Lancet Oncol. 2017;18(10):1373–85. https://doi.org/10.1016/S1470-2045(17)30517-X.
Schalper KA, Rodriguez-Ruiz ME, Diez-Valle R, López-Janeiro A, Porciuncula A, Idoate MA, et al. Neoadjuvant nivolumab modifies the tumor immune microenvironment in resectable glioblastoma. Nat Med. 2019;25(3):470–6. https://doi.org/10.1038/s41591-018-0339-5.
Ahmed N, Brawley V, Hegde M, Bielamowicz K, Kalra M, Landi D, et al. HER2-specific chimeric antigen receptor–modified virus-specific T cells for progressive glioblastoma: a phase 1 dose-escalation trial. Jama Oncol. 2017;3(8):1094–101. https://doi.org/10.1001/jamaoncol.2017.0184.
Chai L-P, Wang ZF, Liang WY, Chen L, Chen D, Wang AX, et al. In vitro and in vivo effect of 5-FC combined gene therapy with TNF-α and CD suicide gene on human laryngeal carcinoma cell line Hep-2. Plos One. 2013;8(4):e61136. https://doi.org/10.1371/journal.pone.0061136.
Herman JM, Wild AT, Wang H, Tran PT, Chang KJ, Taylor GE, et al. Randomized phase III multi-institutional study of TNFerade biologic with fluorouracil and radiotherapy for locally advanced pancreatic cancer: final results. J Clin Oncol. 2013;31(7):886–94. https://doi.org/10.1200/JCO.2012.44.7516.
Lan T, Chen L, Wei X. Inflammatory cytokines in cancer: comprehensive understanding and clinical progress in gene therapy. Cells. 2021;10(1):100. https://doi.org/10.3390/cells10010100.
Yang S, Liu T, Nan H, Wang Y, Chen H, Zhang X, et al. Comprehensive analysis of prognostic immune-related genes in the tumor microenvironment of cutaneous melanoma. J Cell Physiol. 2020;235(2):1025–35. https://doi.org/10.1002/jcp.29018.
Jia D, Li S, Li D, Xue H, Yang D, Liu Y. Mining TCGA database for genes of prognostic value in glioblastoma microenvironment. Aging Albany Ny. 2018;10(4):592–605. https://doi.org/10.18632/aging.101415.
Tang J, Pan R, Xu L, Ma Q, Ying X, Zhao J, et al. IL10 hypomethylation is associated with the risk of gastric cancer. Oncol Lett. 2021;21(4):241. https://doi.org/10.3892/ol.2021.12502.
Kong X, Xu M, Cui X, Ma L, Cheng H, Hou J, et al. Potential role of macrophage phenotypes and CCL2 in the pathogenesis of Takayasu arteritis. Front Immunol. 2021;12:646516. https://doi.org/10.3389/fimmu.2021.646516.
Loreth D, Schuette M, Zinke J, Mohme M, Piffko A, Schneegans S, et al. CD74 and CD44 expression on CTCs in cancer patients with brain metastasis. Int J Mol Sci. 2021;22(13):6993. https://doi.org/10.3390/ijms22136993.
This study was supported by the National Natural Science Foundation of China (No. 82173384 and 81773161).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
About this article
Cite this article
Chen, D., Yao, J., Hu, B. et al. New biomarker: the gene HLA-DRA associated with low-grade glioma prognosis. Chin Neurosurg Jl 8, 12 (2022). https://doi.org/10.1186/s41016-022-00278-0
- Low-grade gliomas