- Open Access
HLA class II molecule HLA-DRA identifies immuno-hot tumors and predicts the therapeutic response to anti-PD-1 immunotherapy in NSCLC
BMC Cancer volume 22, Article number: 738 (2022)
Immune checkpoint blockade (ICB) only works well for a certain subset of patients with non-small cell lung cancer (NSCLC). Therefore, biomarkers for patient stratification are desired, which can suggest the most beneficial treatment.
In this study, three datasets (GSE126044, GSE135222, and GSE136961) of immunotherapy from the Gene Expression Omnibus (GEO) database were analyzed, and seven intersected candidates were extracted as potential biomarkers for ICB followed by validation with The Cancer Genome Atlas (TCGA) dataset and the in-house cohort data.
Among these candidates, we found that human leukocyte antigen-DR alpha (HLA-DRA) was downregulated in NSCLC tissues and both tumor and immune cells expressed HLA-DRA. In addition, HLA-DRA was associated with an inflamed tumor microenvironment (TME) and could predict the response to ICB in NSCLC. Moreover, we validated the predictive value of HLA-DRA in immunotherapy using an in-house cohort. Furthermore, HLA-DRA was related to the features of inflamed TME in not only NSCLC but also in most cancer types.
Overall, HLA-DRA could be a promising biomarker for guiding ICB in NSCLC.
Lung cancer is the most common malignancy, and non-small cell lung cancer (NSCLC) accounts for approximately 85% of lung cancers . The prognosis of patients with lung cancer is significantly poor . The treatment landscape for advanced NSCLC has been changing rapidly. Recent advancements in the molecular profiling of NSCLC enable targeted molecular therapy with EGFR tyrosine kinase inhibitors (TKI) and immune checkpoint inhibitors [3, 4]. Accordingly, patients with EGFR-mutated NSCLC have achieved better progression-free survival (PFS) and overall survival (OS). On the other hand, patients harboring wild-type EGFR and overexpressed PD-L1 may benefit from immune checkpoint blockade (ICB) therapy .
As one of the significant features of malignant tumors, immune escape plays an important role in the oncogenesis, progression, and therapeutic resistance of malignant tumors . In the past decade, immunotherapy targeting immune checkpoints has opened up a new path for treating malignant tumors . Thus far, immunotherapy drugs represented by PD-1/PD-L1 inhibitors have been widely used for treating a variety of malignant tumors, such as NSCLC, triple-negative breast cancer, melanoma, and Hodgkin lymphoma [8,9,10,11]. Among these cancer types, immunotherapy is the most widely used in NSCLC. However, not all patients could benefit from immunotherapy. Therefore, appropriate biomarkers to identify the dominant population are very critical in clinical application.
Currently, the expression of PD-L1 in tumor cells is considered the most reasonable biomarker for selecting the dominant population in anti-PD-1 immunotherapy. The results of multiple clinical trials in advanced NSCLC support the effectiveness of positive expression of PD-L1 in tumor cells in predicting the efficacy of immunotherapy [5, 12, 13]. Nevertheless, CheckMate017 and OAK studies reported that patients with negative PD-L1 expression can also benefit from immunotherapy [14, 15]. In summary, it can be found that the PD-L1 expression level of tumor cells is unreliable as a biomarker for predicting the efficacy of immunotherapy [16, 17]. Our previous studies have been devoted to the study of immunotherapy biomarkers, and proposed that deglycosylated PD-L1 and IFITM3 could be potential biomarkers for immunotherapy [18, 19].
This study aimed to identify novel therapeutic biomarkers for immunotherapy in NSCLC. First, we collected three public datasets of immunotherapy from the Gene Expression Omnibus (GEO) database and screened potential biomarkers for immunotherapy. Given the expression patterns and prognostic values of intersected candidates, human leukocyte antigen-DR alpha (HLA-DRA) was selected for further analysis and validation with the Cancer Genome Atlas (TCGA) dataset and the in-house cohort. Overall, our findings demonstrated that HLA-DRA could be a novel biomarker in NSCLC, which may make up for the deficiency of PD-L1 in the selection of the dominant population of immunotherapy.
Acquisition of public data
The preprocessed RNA-sequencing (RNA-seq) data and immunotherapeutic responses of GSE126044 , GSE135222 , and GSE136961  were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). The pan-cancer (TOIL RSEM tpm) and NSCLC (IlluminaHiSeq) gene expression profiles as well as clinical annotations of the TCGA dataset were obtained from the UCSC Xena (https://xenabrowser.net/datapages/). The abbreviations for various cancer types are exhibited in Supplementary Table S1.
DGEs screening and enrichment analysis
To identify differentially expressed genes (DEGs) between responders and non-responders, student t test was used and genes with P value ≤0.05 were deemed to be candidates. Next, common candidates between GSE126044, GSE135222, and GSE136961 datasets were extracted using Venn analysis. To analyze the biological processes (BP) of intersected genes in these three GEO datasets, the DAVID tool was applied. DAVID (https://david.ncifcrf.gov/) is a widely used gene functional annotation tool . The human genome (homo sapiens) was selected as the background variable. BP terms were deemed to be statistically significant when the P value was ≤0.05 and the top 10 terms were retained.
TISCH (http://tisch.comp-genomics.org/gallery/) is a single-cell RNA-sequencing (scRNA-seq) database focusing on the tumor microenvironment (TME) and provides detailed cell-type annotation at the single-cell level . In this research, we used TISCH to analyze the cell subpopulation patterns of seven intersected genes from three GEO datasets in NSCLC. Default options were used for all parameters.
Immunological correlation of HLA-DRA in NSCLC
The immunological characteristics of TME in NSCLC contained immunomodulators, tumor purity, infiltration levels of tumor-infiltrating immune cells (TIICs), and the expression of inhibitory immune checkpoints. Firstly, we studied the expression of 150 immunomodulators, including MHC, receptors, chemokines, immunoinhibitors, and immunostimulators. In addition, tumor purity was obtained from previous research and its correlation with HLA-DRA was assessed . Moreover, the correlations between HLA-DRA and immune checkpoints levels were evaluated. To avoid the calculation error caused by various algorithms, we used two independent algorithms, including TIMER  and EPIC  to estimate TIICs abundance. To investigate the associations between HLA-DRA and anti-tumor immunity in NSCLC, we divided the samples into the high and low HLA-DRA groups at its median expression, and then analyzed the differences in immunological characteristics of TME in these aspects.
Correlation between HLA-DRA and immunotherapeutic response
According to a previous report, immunophenoscore (IPS) was calculated to predict therapeutic response to immunotherapy . The IPS values of NSCLC patients were obtained from the Cancer Immunome Atlas (TCIA) website (http://tcia.at/home/). In addition, we also computed the T cell inflamed score based on the expression and weighting coefficient of 18 genes . Referring to previous research , we obtained several gene-sets which were correlated with the immunotherapeutic response as well as specific gene markers of biological processes related to anti-tumor immunity (Supplementary Table S2). The enrichment scores of these gene signatures were calculated utilizing the R package “GSVA” . To evaluate the role of HLA-DRA in predicting immunotherapeutic responses, the associations between HLA-DRA expression and these aspects were evaluated.
Pan-cancer analysis of immunological correlation of HLA-DRA
To evaluate the pan-cancer immunological correlation of HLA-DRA, we collected pan-cancer expression of 150 immunomodulators, including MHC, receptors, chemokines, immunoinhibitors, and immunostimulators. Then, the correlations between HLA-DRA and tumor purity as well as TIICs abundance were also assessed.
The NSCLC tissue microarray (TMA, Cat. HLugA150CS03) and the pan-cancer TMA (Cat. HOrgC120PG04) were purchased from Outdo BioTech (Shanghai, China). A total of 75 NSCLC and para-tumor samples were contained in this research. The HOrgC120PG04 microarray contained 11 kinds of cancers with 2–6 tumor samples and para-tumor or normal samples per type. Ethical approval for the use of TMAs was granted by the Clinical Research Ethics Committee in Outdo Biotech (Shanghai, China). In addition, 16 cancer patients who received anti-PD-1 immunotherapy monotherapy or a combination of chemotherapy were recruited by The Affiliated Wuxi People’s Hospital of Nanjing Medical University. The therapeutic response was evaluated according to the RECIST 1.1 criterion. Among them, 8 patients reached partial response (PR), and 8 patients were stable disease (SD). Ethical approval for the collection of tissue sections was granted by the Clinical Research Ethics Committee, The Affiliated Wuxi People’s Hospital of Nanjing Medical University.
Immunohistochemistry (IHC) staining and semi-quantitative scoring
IHC staining was conducted on the above sections according to the standardized procedures. Sections were retrieved by EDTA. The primary antibodies used were as follows: anti-HLA-DRA (1:2500 dilution, Cat. 17,221–1-AP, ProteinTech) and anti-PD-L1 (Ready-to-use, Cat. GT2280, GeneTech). Staining was visualized with DAB and hematoxylin counterstain, and stained sections were captured using Aperio Digital Pathology Slide Scanners. For semi-quantitative analysis, the stained sections were independently evaluated by two pathologists according to the evaluation standard on a 12-point scale by calculating the immunoreactivity score (IRS) .
All statistical analyses were conducted using SPSS 26.0 software or R language. All data are presented as means ± SDs. The difference between the two groups was analyzed by parametric Student’s t-test or non-parametric Mann Whitney test. Survival analysis was performed by log-rank test. Correlation analysis between two variables was analyzed by Pearson test. All statistical tests were two-sided, and P value < 0.05 was considered statistically significant and labeled with *P < 0.05; **P < 0.01; ***P < 0.001.
Identification of potential biomarkers for immunotherapy in NSCLC
To discover novel biomarkers to identify the responsive population to the particular immunotherapy, we downloaded GSE126044, GSE135222, and GSE136961 datasets from the GEO database. DEGs were extracted with P value ≤0.05. A total of 1605 DEGs in the GSE126044 dataset (Fig. 1A), 1166 DEGs in the GSE135222 dataset (Fig. 1B), and 60 DEGs in the GSE136961 dataset (Fig. 1C) were identified. Subsequently, common candidates in these three datasets were extracted using Venn analysis, and a total of seven genes (CD2, HAVCR2, HLA-C, HLA-DRA, HLA-E, KLRK1, and TYROBP) were selected (Fig. 1D). Next, we analyzed the potential functions of intersected genes using the DAVID tool. The results showed that these genes mainly participated in immune-related processes, such as the innate immune response, regulation of immune response, antigen processing and presentation, etc. (Fig. 1E, Supplementary Table S3). Overall, these findings revealed seven potential biomarkers for immunotherapy in NSCLC.
Expression, atlas of cell subtypes, and prognostic values of potential biomarkers in NSCLC
We analyzed the expression of potential biomarkers in NSCLC. In the TCGA dataset, the expression levels of CD2, HAVCR2, HLA-C, HLA-DRA, HLA-E, KLRK1, and TYROBP were all significantly downregulated in tumor tissues compared with para-tumor tissues (Fig. 2A). We also investigated the atlas of cell subtypes of these candidates. The results showed that CD2 was mainly expressed in CD4+ and CD8+ T cells, HAVCR2 was mainly expressed in exhausted T cells, three HLA molecules were expressed in almost all cell types, TYROBP was mainly expressed in monocytes and macrophages, and KLRK1 was lowly expressed in all cell types (Fig. 2B). As to their prognostic values, seven genes could not predict the clinical outcome in NSCLC and lung squamous cell carcinoma (LUSC), while CD2 and HLA-DRA were associated with the improved prognosis in lung adenocarcinoma (LUAD) (Figs. 2C-E). Moreover, the expression levels of these genes were significantly correlated with each other in GSE126044, GSE135222, and GSE136961 datasets (Supplementary Figs. S1A-C). In view of the above findings, we selected HLA-DRA for further investigation and validation.
HLA-DRA is associated with an inflamed TME in NSCLC
Given the correlation between HLA-DRA and immunotherapeutic response, we explored the immunological role of HLA-DRA in NSCLC using the TCGA cohort. Many chemokines, paired receptors, MHC molecules, immunoinhibitors, and immunostimulators were notably upregulated in the high-HLA-DRA group (Fig. 3A). Considering that these chemokines and receptors recruit effector TIICs, we speculated that HLA-DRA was correlated with increased TIICs in TME. Interestingly, HLA-DRA was negatively correlated with tumor purity (Fig. 3B). Furthermore, HLA-DRA was positively correlated with TIIC abundance estimated by TIMER and EPIC algorithms (Fig. 3C). Moreover, HLA-DRA was positively correlated with expression levels of immune checkpoints in NSCLC (Fig. 3D). To validate the above findings, an NSCLC TMA cohort was used (Fig. 4A). The results showed that the expression of HLA-DRA protein was significantly downregulated in tumor tissues (Figs. 4B-C). Furthermore, the current NSLCL cohort was classified into low- and high-HLA-DRA expression groups based on its median level (IRS ≤ 2 vs. IRS > 3). We found that PD-L1 expression was higher in the high-HLA-DRA group (Figs. 4D-E). Furthermore, HLA-DRA was positively correlated with PD-L1 expression in the NSCLC cohort (Fig. 4F). Overall, HLA-DRA is highly correlated with the inflamed TME and identified immuno-hot tumors in NSCLC.
HLA-DRA predicts the response to immunotherapy in NSCLC
We further validated whether NSCLC patients with high HLA-DRA expression exhibited high responses to immunotherapy. IPS was a surrogate of the response to immunotherapy, and we found that patients with high HLA-DRA expression had notably high IPS (Fig. 5A). T cell-inflamed score was developed using IFN-γ-related mRNA profiles to predict response to anti-PD-1 therapy. The results showed that HLA-DRA was positively correlated with T cell-inflamed scores in NSCLC (Fig. 5B). Moreover, the immunotherapy-related enrichment scores, such as IFN-γ signature, APM signal, and hypoxia, were positively correlated with HLA-DRA expression, while FGFR3-coexpressed genes, PPAR-γ network, WNT/β-catenin network, and VEGFA pathway were negatively correlated HLA-DRA expression (Fig. 5C). To further validate the association between HLA-DRA expression and immunotherapeutic response, a small-scale immunotherapy cohort was used (Fig. 5D). The results showed that HLA-DRA and PD-L1 were significantly upregulated in the patients with the response of PR than SD (Figs. 5E-F). Collectively, public and in-house cohorts support that HLA-DRA could be a novel biomarker for predicting the response to anti-PD-1 immunotherapy in NSCLC.
Extension of the immunological role of HLA-DRA in pan-cancer
Our data illustrated that HLA-DRA was associated with an inflamed TME and could identify immuno-hot tumors in NSCLC. However, the immunological role of HLA-DRA in other cancer types was unclear. Next, we analyzed the correlations between HLA-DRA and chemokine, MHC, receptor, immunoinhibitors, and immunostimulators. Except for several cancer types, HLA-DRA was positively correlated with the expression levels of these immunomodulators (Fig. 6A). Moreover, HLA-DRA was negatively correlated with tumor purity but positively correlated with TIIC levels in most cancer types (Figs. 6B-C). HLA-DRA was also positively correlated with the expression levels of immune checkpoints in pan-cancer (Supplementary Figs. S2). We next collected a pan-cancer TMA cohort to validate the above findings (Fig. 6D). The results showed that HLA-DRA was highly correlated with PD-L1 expression in the current cohort (Fig. 6E-F). Taken together, the data suggested that HLA-DRA is a pan-cancer classifier for immuno-hot tumors except for a few tumor types.
Anti-PD-1 immunotherapy exhibited encouraging efficacy for advanced NSCLC, but not all patients could benefit from immunotherapy. Predictive signatures to identify the dominant population are urgent for the field. TME, consisting of multiple cell types and the interplay between these cells via cytokines, chemokines, and growth factors , lays the basis for determining whether immunotherapy is effective . According to the features of TME, tumors can be classified into immuno-cold or immuno-hot. Immuno-cold tumors are featured with immunosuppressive TME and insensitive immunotherapy, but immuno-hot tumors represent the high response to immunotherapy, which is accompanied by active T cell infiltration . Therefore, distinguishing immuno-hot or immuno-cold tumors using practical biomarkers is a significant strategy to demarcate the response to immunotherapy.
In the current research, we obtained three public RNA-seq datasets from immunotherapy cohorts to identify novel biomarkers for immunotherapy, and reported that CD2, HAVCR2, HLA-C, HLA-DRA, HLA-E, KLRK1, and TYROBP might be novel biomarkers for immunotherapy in NSCLC. In view of expression, atlas of cell subtypes, and prognostic values, we selected HLA-DRA for further investigation and validation. We found that HLA-DRA was associated with an inflamed TME and identified immuno-hot tumors in NSCLC. Additionally, HLA-DRA could predict the response to immunotherapy from the point of multiple signatures and in our recruited cohort. Moreover, pan-cancer investigation and validation revealed that HLA-DRA was a pan-cancer classifier for immuno-hot tumors.
Previous research has proved that HLA-DR participates in the inhibition of tumor growth. HLA-DR presents tumor-associated antigens (TAA) that are recognized by CD4+ T cells, which then produce cytokines, including interleukins and interferon-γ, to suppress tumor growth [36, 37]. HLA-DRA is a subunit of HLA-DR, which plays a critical role in human cancers. Overexpression of HLA-DRA was reported in hepatocellular cancer , colorectal cancer , and cervical cancer , while it was decreased in breast cancer . In this report, we found that HLA-DRA was significantly decreased in NSCLC. In addition, HLA-DRA has also been shown to function as a prognostic biomarker for clinical outcomes [42, 43]. Nevertheless, no research has uncovered the predictive role of HLA-DRA for immunotherapy in NSCLC.
In addition to the observational results from the public and in-house cohorts, we found that HLA-DRA was associated with multiple immunotherapy-related pathways. HLA-DRA was positively correlated with IFN-γ signature, APM signal, and hypoxia, while was negatively correlated with FGFR3-coexpressed genes, PPAR-γ network, WNT/β-catenin network, and VEGFA pathway. IFN-γ, APM, and hypoxia pathways were correlated with high response to immunotherapy [44, 45], while FGFR3, β-catenin, PPAR-γ, and VEGFA pathways were activated in the non-inflamed tumors [46, 47]. The interplays between HLA-DRA and IFN-γ signature as well as the APM signal were clear [48, 49]. For example, IFN-γ could increase the expression of the HLA-DRA . However, whether HLA-DRA could regulate other immunotherapy-related pathways needs to be further studied.
This study identified potential biomarkers to predict the response to immunotherapy. After systematic analysis, HLA-DRA could a promising biomarker that is associated with an inflamed TME and may predict the response to immunotherapy in NSCLC. Furthermore, the pan-cancer analysis revealed that HLA-DRA identifies high immunogenicity. Overall, HLA-DRA might be a promising biomarker for guiding immunotherapy in NSCLC.
Availability of data and materials
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors. In addition, original data for bioinformatics analysis could be downloaded from corresponding platforms.
Non-small cell lung cancer
Tyrosine kinase inhibitor
Immune checkpoint blockade
Gene Expression Omnibus
Human leukocyte antigen-DR alpha
The Cancer Genome Atlas
Differentially expressed gene
Tumor-infiltrating immune cell
The Cancer Immunome Atlas
Lung squamous cell carcinoma
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2022. CA Cancer J Clin. 2022;72(1):7–33.
Zhang S, He Y, Xuan Q, Ling X, Men K, Zhao X, et al. TMEM139 prevents NSCLC metastasis by inhibiting lysosomal degradation of E-cadherin. Cancer Sci. 2022.
Passaro A, Mok T, Peters S, Popat S, Ahn MJ, de Marinis F. Recent advances on the role of EGFR tyrosine kinase inhibitors in the management of NSCLC with uncommon, non exon 20 insertions. EGFR Mutations J Thorac Oncol. 2021;16(5):764–73.
Anagnostou VK, Brahmer JR. Cancer immunotherapy: a future paradigm shift in the treatment of non-small cell lung cancer. Clin Cancer Res. 2015;21(5):976–84.
Mok TSK, Wu YL, Kudaba I, Kowalski DM, Cho BC, Turna HZ, et al. Pembrolizumab versus chemotherapy for previously untreated, PD-L1-expressing, locally advanced or metastatic non-small-cell lung cancer (KEYNOTE-042): a randomised, open-label, controlled, phase 3 trial. Lancet. 2019;393(10183):1819–30.
Anichini A, Perotti VE, Sgambelluri F, Mortarini R. Immune escape mechanisms in non small cell lung Cancer. Cancers (Basel). 2020;12(12).
Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12(4):252–64.
Robert C, Schachter J, Long GV, Arance A, Grob JJ, Mortier L, et al. Pembrolizumab versus Ipilimumab in Advanced Melanoma. N Engl J Med. 2015;372(26):2521–32.
Herbst RS, Baas P, Kim DW, Felip E, Perez-Gracia JL, Han JY, et al. Pembrolizumab versus docetaxel for previously treated, PD-L1-positive, advanced non-small-cell lung cancer (KEYNOTE-010): a randomised controlled trial. Lancet. 2016;387(10027):1540–50.
Nanda R, Chow LQ, Dees EC, Berger R, Gupta S, Geva R, et al. Pembrolizumab in patients with advanced triple-negative breast Cancer: phase Ib KEYNOTE-012 study. J Clin Oncol. 2016;34(21):2460–7.
Armand P, Shipp MA, Ribrag V, Michot JM, Zinzani PL, Kuruvilla J, et al. Programmed Death-1 blockade with Pembrolizumab in patients with classical Hodgkin lymphoma after Brentuximab Vedotin failure. J Clin Oncol. 2016;34(31):3733–9.
Brahmer JR, Rodriguez-Abreu D, Robinson AG, Hui R, Csoszi T, Fulop A, et al. Health-related quality-of-life results for pembrolizumab versus chemotherapy in advanced, PD-L1-positive NSCLC (KEYNOTE-024): a multicentre, international, randomised, open-label phase 3 trial. Lancet Oncol. 2017;18(12):1600–9.
Langer CJ, Gadgeel SM, Borghaei H, Papadimitrakopoulou VA, Patnaik A, Powell SF, et al. Carboplatin and pemetrexed with or without pembrolizumab for advanced, non-squamous non-small-cell lung cancer: a randomised, phase 2 cohort of the open-label KEYNOTE-021 study. Lancet Oncol. 2016;17(11):1497–508.
Brahmer J, Reckamp KL, Baas P, Crino L, Eberhardt WE, Poddubskaya E, et al. Nivolumab versus docetaxel in advanced squamous-cell non-small-cell lung Cancer. N Engl J Med. 2015;373(2):123–35.
Rittmeyer A, Barlesi F, Waterkamp D, Park K, Ciardiello F, von Pawel J, et al. Atezolizumab versus docetaxel in patients with previously treated non-small-cell lung cancer (OAK): a phase 3, open-label, multicentre randomised controlled trial. Lancet. 2017;389(10066):255–65.
Cha JH, Chan LC, Li CW, Hsu JL, Hung MC. Mechanisms controlling PD-L1 expression in Cancer. Mol Cell. 2019;76(3):359–70.
Chen N, Higashiyama N, Hoyos V. Predictive biomarkers of immune checkpoint inhibitor response in breast Cancer: looking beyond Tumoral PD-L1. Biomed. 2021;9(12).
Mei J, Xu J, Yang X, Gu D, Zhou W, Wang H, et al. A comparability study of natural and deglycosylated PD-L1 levels in lung cancer: evidence from immunohistochemical analysis. Mol Cancer. 2021;20(1):11.
Cai Y, Ji W, Sun C, Xu R, Chen X, Deng Y, et al. Interferon-induced transmembrane protein 3 shapes an inflamed tumor microenvironment and identifies Immuno-hot tumors. Front Immunol. 2021;12(3162).
Cho JW, Hong MH, Ha SJ, Kim YJ, Cho BC, Lee I, et al. Genome-wide identification of differentially methylated promoters and enhancers associated with response to anti-PD-1 therapy in non-small cell lung cancer. Exp Mol Med. 2020;52(9):1550–63.
Jung H, Kim HS, Kim JY, Sun JM, Ahn JS, Ahn MJ, et al. DNA methylation loss promotes immune evasion of tumours with high mutation and copy number load. Nat Commun. 2019;10(1):4278.
Hwang S, Kwon AY, Jeong JY, Kim S, Kang H, Park J, et al. Immune gene signatures for predicting durable clinical benefit of anti-PD-1 immunotherapy in patients with non-small cell lung cancer. Sci Rep. 2020;10(1):643.
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.
Sun D, Wang J, Han Y, Dong X, Ge J, Zheng R, et al. TISCH: a comprehensive web resource enabling interactive single-cell transcriptome visualization of tumor microenvironment. Nucleic Acids Res. 2021;49(D1):D1420–D30.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, et al. The immune landscape of Cancer. Immunity. 2018;48(4):812–30 e14.
Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–W14.
Racle J, de Jonge K, Baumgaertner P, Speiser DE, Gfeller D. Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data. Elife. 2017;6.
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer Immunogenomic analyses reveal genotype-Immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62.
Ayers M, Lunceford J, Nebozhyn M, Murphy E, Loboda A, Kaufman DR, et al. IFN-gamma-related mRNA profile predicts clinical response to PD-1 blockade. J Clin Invest. 2017;127(8):2930–40.
Hu J, Yu A, Othmane B, Qiu D, Li H, Li C, et al. Siglec15 shapes a non-inflamed tumor microenvironment and predicts the molecular subtype in bladder cancer. Theranostics. 2021;11(7):3089–108.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Mei J, Liu Y, Yu X, Hao L, Ma T, Zhan Q, et al. YWHAZ interacts with DAAM1 to promote cell migration in breast cancer. Cell Death Dis. 2021;7(1):221.
Duan Q, Zhang H, Zheng J, Zhang L. Turning cold into hot: firing up the tumor microenvironment. Trends Cancer. 2020;6(7):605–18.
Zemek RM, De Jong E, Chin WL, Schuster IS, Fear VS, Casey TH, et al. Sensitization to immune checkpoint blockade through activation of a STAT1/NK axis in the tumor microenvironment. Sci Transl Med. 2019;11(501).
Gajewski TF. The next hurdle in Cancer immunotherapy: overcoming the non-T-cell-inflamed tumor microenvironment. Semin Oncol. 2015;42(4):663–71.
Sabbatino F, Liguori L, Polcaro G, Salvato I, Caramori G, Salzano FA, et al. Role of human leukocyte antigen system as a predictive biomarker for checkpoint-based immunotherapy in Cancer patients. Int J Mol Sci. 2020;21(19).
Seliger B, Kloor M, Ferrone S. HLA class II antigen-processing pathway in tumors: molecular defects and clinical relevance. Oncoimmunology. 2017;6(2):e1171447.
Matoba K, Iizuka N, Gondo T, Ishihara T, Yamada-Okabe H, Tamesa T, et al. Tumor HLA-DR expression linked to early intrahepatic recurrence of hepatocellular carcinoma. Int J Cancer. 2005;115(2):231–40.
Matsushita K, Takenouchi T, Shimada H, Tomonaga T, Hayashi H, Shioya A, et al. Strong HLA-DR antigen expression on cancer cells relates to better prognosis of colorectal cancer patients: possible involvement of c-myc suppression by interferon-gamma in situ. Cancer Sci. 2006;97(1):57–63.
Balakrishnan CK, Tye GJ, Balasubramaniam SD, Kaur G. CD74 and HLA-DRA in cervical carcinogenesis: potential targets for antitumour therapy. Medicina. 2022;58(2).
da Silva GB, Silva TG, Duarte RA, Neto NL, Carrara HH, Donadi EA, et al. Expression of the classical and nonclassical HLA molecules in breast Cancer. Int J Breast Cancer. 2013;2013:250435.
Piao XM, Kang HW, Jeong P, Byun YJ, Lee HY, Kim K, et al. A prognostic immune predictor, HLA-DRA, plays diverse roles in non-muscle invasive and muscle invasive bladder cancer. Urol Oncol 2021;39(4):237 e21- e29.
Axelrod ML, Cook RS, Johnson DB, Balko JM. Biological consequences of MHC-II expression by tumor cells in Cancer. Clin Cancer Res. 2019;25(8):2392–402.
Bretz AC, Parnitzke U, Kronthaler K, Dreker T, Bartz R, Hermann F, et al. Domatinostat favors the immunotherapy response by modulating the tumor immune microenvironment (TIME). J Immunother Cancer. 2019;7(1):294.
Hong S, Zhang Y, Cao M, Lin A, Yang Q, Zhang J, et al. Hypoxic characteristic genes predict response to immunotherapy for urothelial carcinoma. Front Cell Dev Biol. 2021;9:762478.
Sweis RF, Spranger S, Bao R, Paner GP, Stadler WM, Steinberg G, et al. Molecular drivers of the non-T-cell-inflamed tumor microenvironment in urothelial bladder Cancer. Cancer Immun Res. 2016;4(7):563–8.
Wang Q, Gao J, Di W, Wu X. Anti-angiogenesis therapy overcomes the innate resistance to PD-1/PD-L1 blockade in VEGFA-overexpressed mouse tumor models. Cancer Immunol Immunother. 2020;69(9):1781–99.
Antohe I, Tanasa MP, Dascalescu A, Danaila C, Titieanu A, Zlei M, et al. The MHC-II antigen presentation machinery and B7 checkpoint ligands display distinctive patterns correlated with acute myeloid leukaemias blast cells HLA-DR expression. Immunobiology. 2021;226(1):152049.
Basham TY, Merigan TC. Recombinant interferon-gamma increases HLA-DR synthesis and expression. J Immunol. 1983;130(4):1492–4.
This work was partially supported by the Natural Science Foundation of Jiangsu Province (BK20210068), the Top Talent Support Program for Young and Middle-aged People of Wuxi Municipal Health Commission (HB2020003), the High-end Medical Expert Team of the 2019 Taihu Talent Plan (2019-THRCTD-1), and the Beijing Great Physician Commonweal Foundation (DY-Tumot2022-J001).
Ethics approval and consent to participate
Ethical approval for the study of tissue microarray slides was granted by the Clinical Research Ethics Committee, Outdo Biotech (Shanghai, China). Ethical approval for the collection of tissue sections was granted by the Clinical Research Ethics Committee, The Affiliated Wuxi People’s Hospital of Nanjing Medical University. All methods in this study were carried out in accordance with relevant guidelines and regulations. Informed consent was obtained from all subjects.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Mei, J., Jiang, G., Chen, Y. et al. HLA class II molecule HLA-DRA identifies immuno-hot tumors and predicts the therapeutic response to anti-PD-1 immunotherapy in NSCLC. BMC Cancer 22, 738 (2022). https://doi.org/10.1186/s12885-022-09840-6