key: cord-0003608-6m4q219k authors: Gao, Mengqi; Lin, Yi; Liu, Xing; Li, Yiming; Zhang, Chuanbao; Wang, Zheng; Wang, Zhiliang; Wang, Yulin; Guo, Zongze title: ISG20 promotes local tumor immunity and contributes to poor survival in human glioma date: 2018-10-31 journal: Oncoimmunology DOI: 10.1080/2162402x.2018.1534038 sha: 57b88add1645a89614742dc32b0e11922c5f37d2 doc_id: 3608 cord_uid: 6m4q219k Recent evidence has confirmed that a mutation of the isocitrate dehydrogenase (IDH) gene occurs early in gliomagenesis and contributes to suppressed immunity. The present study aimed to determine the candidate genes associated with IDH mutation status that could serve as biomarkers of immune suppression for improved prognosis prediction. Clinical information and RNA-seq gene expression data were collected for 932 glioma samples from the CGGA and TCGA databases, and differentially expressed genes in both lower-grade glioma (LGG) and glioblastoma (GBM) samples were identified according to IDH mutation status. Only one gene, interferon-stimulated exonuclease gene 20 (ISG20), with reduced expression in IDH mutant tumors, demonstrated significant prognostic value. ISG20 expression level significantly increased with increasing tumor grade, and its high expression was associated with a poor clinical outcome. Moreover, increased ISG20 expression was associated with increased infiltration of monocyte-derived macrophages and neutrophils, and suppressed adaptive immune response. ISG20 expression was also positively correlated with PD-1, PD-L1, and CTLA4 expression, along with the levels of several chemokines. We conclude that ISG20 is a useful biomarker to identify IDH-mediated immune processes in glioma and may serve as a potential therapeutic target. Glioma is the most common and lethal type of malignancy in the primary central nervous system (CNS). 1 Although patients with low-grade gliomas (LGGs) have a more favorable prognosis than those with glioblastomas (GBMs), many tend to progress to a higher grade, leading to poor survival. 2 Nevertheless, the outcome of glioma patients is highly variable, even among those with the same tumor grade. 3 Recent analyses demonstrated that an IDH1/2 mutation, encoding isocitrate dehydrogenase (IDH) gene, occurs early in gliomagenesis, affecting a common glial precursor cell population. 4 Patients with tumors harboring an IDH1/2 mutation (IDHmut) show significantly longer survival than those expressing wild-type IDH1/2 (IDHwt). [5] [6] [7] IDH mutation leads to a CpG island methylator phenotype (CIMP) by modulating the methylation patterns on a genome-wide scale, changing transcriptional programs and altering the differentiation state. 8 CIMP is associated with microsatellite instability and longer survival in several cancers. 6, 7, [9] [10] [11] [12] [13] IDHmut and IDHwt tumors differ with regards to various biological processes, including immune cell infiltration. [14] [15] [16] [17] Human IDH1-mutant gliomas have less infiltrating immune cells than IDH1-wild type gliomas, with global depletion of immune infiltrates, including microglia, macrophages, dendritic cells, B cells, and T cells. Accordingly, early IDHmut glioma progenitor cells have suppressed immunity compared with IDHwt cells, 4, 15, 18 which may be responsible for their improved clinical outcomes. 15 Moreover, IDHmut tumors have reduced expression of cytotoxic T lymphocyte-associated genes and interferon (IFN)-γ-inducible chemokines, as well as suppressed accumulation of T cells in the tumor compared with IDHwt tumors. 16 IDHwt gliomas are also characterized by more prominent regulatory T cell infiltration and higher programmed death-ligand 1 (PD-L1) expression levels than IDHmut cases. 17 Although IDH status clearly appears to affect the immune state and progression of glioma, the underlying mechanisms remain unclear. To elucidate these mechanisms and identify the candidate prognostic and/or therapeutic markers, we investigated the differential expression of immune-related genes and their role in glioma progression. In particular, we collected clinical and transcriptome (RNA-seq) data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) databases, including 932 glioma samples. We then determined the differentially expressed immune-related genes according to IDH mutation status, analyzed separately for LGG and GBM cases, and performed pathway enrichment analysis for functional annotation. To identify the differentially expressed immune-related genes according to IDH status, we compared their expression levels between IDHwt and IDHmut tumors. Genes were analyzed in four groups: LGG in the CGGA database (CGGA-LGG), GBM in the CGGA database (CGGA-GBM), LGG in the TCGA database (TCGA-LGG), and GBM in the TCGA database (TCGA-GBM). Twelve genes with upregulated expression in IDHmut gliomas and 71 genes with upregulated expression in IDHwt tumors were found to be significant across all four cohorts ( Figure 1 , Table S1 ). The prognostic value of these genes was further evaluated (Table S2) . Only one gene, ISG20, with up-regulated expression in IDHwt tumors, was found to have a consistently significant influence on patient survival across all groups. High ISG20 expression is related to higher tumor malignancy in glioma To further investigate the clinical significance of ISG20 in glioma, we analyzed the expression of ISG20 according to the World Health Organization (WHO) classification and molecular subtypes. We found that ISG20 expression level increased with increasing tumor grade (Figure 2(a); Fig. S1A ). The mesenchymal subtype of GBM showed the highest ISG20 expression level compared with the other subtypes, whereas the proneural subtype had the lowest ISG20 level (Figure 2 (b); Fig. S1B ). Moreover, ISG20 expression was up-regulated in 1p/ 19q-non-co-deleted gliomas compared with those with 1p/19qco-deletion (Figure 2(c); Fig. S1C ). In addition, ISG20 expression was significantly higher in IDHwt tumors compared with IDHmut tumors across all tumor grades (Figure 2 (d-f), Fig. S1D-F) . Overall, these results suggest that ISG20 expression is suppressed in IDH mutant tumors and is associated with higher tumor malignancy. To further determine the prognostic value of ISG20, we carried out Cox regression analyses in the CGGA and TCGA databases, respectively. In the CGGA database, ISG20 emerged as an independent prognostic factor of patient survival (hazard ratio = 1.058, 95%CI = 1.024-1.093, p = 0.001), independent of age, tumor grade, chemo-and radiation Figure 1 . Numbers of differentially expressed genes according to IDH status in four groups of samples. Genes were analyzed in four cohorts: low-grade gliomas in CGGA cohort (CGGA LGGs), glioblastomas in CGGA cohort (CGGA GBMs), low-grade gliomas in TCGA cohort (TCGA LGGs), and glioblastomas in TCGA cohort (TCGA GBMs). Genes that were both significant in Student's t-test (p < 0.05) and SAM (FDR < 0.01) were counted. Each digit represents the number of significantly upregulated genes in the corresponding group. The overlapped genes across four groups were selected for further survival analysis. therapy, IDH mutation, 1p19q-co-deletion status, and MGMT promoter methylation (Table 1) . Analysis in TCGA database showed the same result (hazard ratio = 1.250, 95%CI = 1.070-1.460, p = 0.005, Table 1 ). The samples were further stratified based on IDH status and tumor grades. As shown in Figure 3 and Fig. S2 , high ISG20 expression (i.e. higher than the median ISG20 level) was associated with shorter patient survival across all subgroups, except for grade II glioma (Figure 3, Fig. S2 ). To examine the consistency of ISG20 protein expression with gene expression findings, we used immunohistochemistry and measured ISG20 expression in tumor specimens from 43 glioma patients. In line with the results from the bioinformatics analyses, protein expression of ISG20 showed the same results, with increased expression in higher tumor grade and decreased expression in IDH mutant tumors. Representative immunohistochemical staining of ISG20 in gliomas is also illustrated (Figure 2 (e,f). To better understand the mechanism of ISG20-promoted glioma progression, we determined the genes showing correlated expression with ISG20 among the whole transcriptome. Among the total 21,468 genes analyzed, 570 genes in the CGGA and 1925 genes in TCGA showed significant correlations with ISG20 expression (r > 0.5, Table S3 ). Subsequent pathway analyses showed that these genes were involved in a wide variety of immune responses, including innate immune response, T cell activation, inflammatory response, antigen processing and presentation, hematopoietic cell lineage, leukocyte migration, cytokine-cytokine receptor interaction, lysosome activation, and type I and II IFN signaling pathway (Figure 4; Fig. S3A ). ISG20 also appears to be involved in essential metabolic pathways, such as the regulation of protein processing and modification, proteasome synthesis, and galactose metabolism (data not shown). Collectively, these results indicated that ISG20 is highly involved with the immune response in the glioma microenvironment. To validate this hypothesis, we further investigated which pathways were associated with ISG20 expression. The Gene Set Enrichment Analysis (GSEA) highlighted that the high expression of ISG20 was significantly enriched in multiple immune processes, especially in innate immunity (Fig S3B; Table S4 ). The correlationship of ISG20 expression with the above immune processes was further validated in detail using Gene Set Variation Analysis (GSVA; Fig S3C ; Table S4 ). To further explore the connection of ISG20 and immunity in the local tumor microenvironment, we performed gene set variation analysis (GSVA) and analyzed the enrichment scores of various immune cell-characterized gene sets in ISG20-associated genes. We found that as the ISG20 expression level increases, the numbers of macrophages and neutrophils that infiltrates the tumor also rise, while the numbers of central memory T cells (Tcm), follicular helper T cells (Tfh), effector memory T cells (Tem), and gamma delta T cells (Tgd) decreased ( Figure 5(a); Fig. S4A ). However, NK cells, CD4 + T cells, and CD8 + T cells were not correlated with ISG20 expression. Mutual relationship analysis of ISG20related immune cells was further conducted to identify their clusters according to the extent of ISG20 expression ( Figure 5 (b), Fig. S4B ), demonstrating that macrophages and neutrophils were tightly correlated, both of which were positively correlated with ISG20 expression; while subsets of T cells were correlated with each other, and were negatively correlated with ISG20 levels. To further confirm these initial findings, we next studied the intratumoral immune cell infiltrates from 43 glioma patients with immunostaining and compared their immune cell profiles relative to ISG20 expression ( Figure 5(c,d) ). In line with the bioinformatic findings, the composition of the immune infiltrates differed between the low ISG20 expression group and the high ISG20 expression group, with a more substantial number of macrophages and lower number of the T cell subsets (T cells, Tfh cells, memory T cells) in high ISG20 expression tumors. Moreover, microglia and cytotoxic T cells were not associated with the ISG20 protein level. We also investigated the association of ISG20 and chemokine levels in the tumor. As expected, the ISG20 expression level increased in parallel with the levels of many chemokines, including CCL2, CCL5, CCL23, CCL26, CCR1, CCR2, CCR7, CCRL2, CXCL16, CXCR4, and CXCR6 (r > 0.4, p < 0.05; data not shown). Since PD-1/PD-L1 and CTLA-4 are important immune inhibitory molecules, we further determined the association of ISG20 expression with these levels. We found that PD-1, PD-L1, and CTLA-4 expression levels increased in parallel with increasing ISG20 levels (r > 0.45; p < 0.0001) (Fig. S5 ). Mounting evidence suggests that glioma with an IDH mutation has suppressed immunity. 16 Inflammation may promote the proliferation and survival of malignant cells and metastasis, contribute to overturning adaptive immunity, and change tumor responses to chemotherapeutic agents. 19 A large proportion of the glioma tumor microenvironment consists of an inflammatory infiltrate predominated by macrophages, which are thought to be subverted by glioblastoma cells for tumor growth. 20, 21 In gliomas, IDHwt tumors are associated with infiltration of more immune-related components, contributing to poor survival. However, the associated molecular mechanisms remain poorly understood, and thus more extensive research in this area could provide insight toward improved diagnosis and treatment. [22] [23] [24] [25] Toward this end, we have identified that ISG20 may play a key role in mediating the IDH-related immune response in the glioma microenvironment. The ISG20 gene, located at chromosome 15q26.1, encodes a 20-kDa protein. It was first introduced in 1997 by Celine Gongora and colleagues as an IFN-induced promyelocytic leukemia nuclear body (PML-NB)-associated protein. 26 PML-NB is reported to be extensively involved in oncogenesis and gene transcription, [27] [28] [29] and the function differs according to its binding protein. [30] [31] [32] [33] The endogenous ISG20 protein is present both in the nucleolus and in the Cajal bodies, and participates in the maturation of small nucleolar RNAs and ribosomal RNAs, and in ribosome biogenesis, [34] [35] [36] [37] [38] [39] leading to the control of RNA stability. 40,41 ISG20 can be induced by type I and type II IFNs, 42,43 and its expression has been shown to be elevated during infection [44] [45] [46] and in several types of cancers. 33, [47] [48] [49] [50] [51] However, the exact function and mechanism of action of this protein remain elusive. In this study, we found that ISG20 expression is associated with many chemokines, leading to tumor infiltration of a variety of immune cells. CCL2/CCR2 is known for recruiting monocytes to the sites of inflammation produced by either tissue injury or infection, and augments the accumulation of regulatory Foxp3(+)CD4(+) T cells and of nitric oxide-and YM-1-expressing macrophages and microglia. 52 CCL5 induces chemotaxis in T cells and monocytes. 53-56 CCL23 might play vital roles in inflammation through the recruitment of macrophages and dendritic cells. 57 CCRL2/CXCR2 is the main neutrophil attractor in vitro, 58 and is involved in the control of both innate and adaptive immune responses. 59 CXCR6, a chemokine receptor for CXCL16 that is expressed on a subset of CD4 + T helper 1 cells and natural killer T cells, is involved in lymphocyte homing and modulates the development and progression of atherosclerosis. 60 Collectively, elevated expression of these chemokines may contribute to glioma progression by recruiting macrophages and neutrophils to the local tumor environment. Nevertheless, we did not find a significant association of microglia with ISG20 expression, suggesting that the ISG20-associated macrophages were monocyte-derived instead of residual microglial cells. 61 Moreover, we found that the majority of adaptive immune cells subtypes, such as Tgd, Tcm, Tfh, and Tem were less abundant, with increased expression of ISG20, which suggested that high ISG20 expression is associated with inefficient antitumor immunity. However, we did not find any correlation between ISG20 and natural killer cells, CD4 + T cells, and CD8 + T cells. We also detected a positive correlation of ISG20 levels with PD-1/PD-L1 and CTLA4 expression, further inhibiting T cell function and leading to tumor evasion of immune responses. 62 The association of ISG20 expression with patient survival did not reach statistical significance for WHO grade II tumors and in the IDHmut tumors in TCGA cohort, although we found the same trend as the other groups. This may be due to the relatively longer survival of these patients. Further studies with more extended follow-up are needed to determine the prognostic value of ISG20 in these subgroups of patients. In conclusion, our observations suggest that ISG20 strongly correlates with the tumor-induced immune response and infiltration of monocyte-derived macrophages and neutrophils, which increases immunity in the tumor environment. Accordingly, ISG20 may serve as a useful prognostic marker for glioma and as a potential therapeutic target. Further studies are warranted to confirm these results and unveil the underlying mechanism. The CGGA (http://www.cgga.org.cn) and TCGA (http://can cergenome.nih.gov) datasets were downloaded online. The transcriptome expression data of 932 gliomas were collected, ranging from WHO grade II to grade IV. In the CGGA dataset, there were 325 samples, including 144 GBM and 181 LGG samples. 63 In the TCGA dataset, 150 GBM and 457 LGG samples were available. Information on age, gender, diagnosis, WHO grade, chemo-and radiation therapy regimen, molecular data, and patient prognoses was also collected. Tumor tissues were also collected from 43 glioma patients who underwent surgery at the First Hospital of China Medical University between 2017 and 2018, including 19 WHO grade II patients, 12 grade III patients, and 12 GBM patients (Table S5 ). This study was approved by the ethics committee of the First Hospital of China Medical University, the institutional review boards of Capital Medical University, and the Beijing Institute for Brain Disorders Brain Tumor Center. Written informed consent was obtained from all patients. The IDH mutation status of samples from the CGGA dataset was determined by aligning to the human reference genome (Hg19 Refseq) using IGV (version 2.3.93), which was previously shown to have good correlation with pyrosequencing results. 64, 65 The IDH status was available directly from the TCGA dataset. The 1p/19q co-deletion status was downloaded from the two datasets respectively. Detection of IDH1 and IDH2 status of 43 glioma patients in China Medical University was performed in the pathology department of the hospital by immunohistochemistry (antihuman IDH1R132H H09) and sequencing respectively. Immunohistochemistry was performed using formalin-fixed, paraffin-embedded tissues. Four μm-thick sections were cut and dewaxed in xylene, rinsed in graded ethanol, and rehydrated in distilled water. After antigen retrieval with sodium citrate buffer (10 mM sodium citrate, pH 6.0), endogenous peroxidase activity was blocked with 3% H 2 O 2 . Then the sections were added with the appropriate antibody at appropriate dilution, then DAB staining solution was added. The sections were counter-stained with hematoxylin, dehydrated, and sealed. For immunofluorescence, antigen was retrieved with EDTA buffer (1 mM Tris/EDTA, pH 9.0) and endogenous fluorescence was eliminated with AutoFluo Quencher (Servicebio, Cat# G1221). The details of the antibodies used for immunofluorescence and immunohistochemistry are given in Table S6 . Staining for each marker was scored using Image-pro plus (v.6.0). Integrated optical density (IOD) to area ratio was calculated for each marker to assess the staining intensity. Immune-related genes were collected from the canonical biological pathways in the Molecular Signatures Database v. 4.0 (MSigDB, http://www.broad.mit.edu/gsea/msigdb/), 66,67 combined with genes identified in three publications from the Journal of Allergy and Clinical Immunology, [68] [69] [70] and genes selected among those involved in canonical immune processes indicated in Janeway's Immunobiology (9th edition). A total of 829 immune genes were confirmed by removing overlapping genes. R language (v. 3.4.3), SPSS software (v. 22.0), and GraphPad Prism (v. 7.0) for Windows were used for statistical analyses and generating figures. GBM and LGG samples from the CGGA and TCGA datasets were analyzed respectively. Genes with significantly different expression between groups with distinct IDH status were estimated by a two-tailed Student's t-test. The Significance Analysis of Microarrays (SAM) package of R was performed to control the FDR. Values of p < 0.05 and FDR < 0.01 were considered statistically significant. Genes that showed consistent differential expression in both the GBM and LGG cohorts from the two datasets were extracted. We then calculated the prognostic value of these differentially expressed genes using the survival package of R. A multivariate Cox proportional hazard model was performed for evaluating the independent prognostic variables. Kaplan-Meier curves were employed to depict survival distributions. Genes and immune cells correlated with ISG20 expression were explored by Pearson's correlation coefficient (r) using R considering the effect of variant grades or IDH status. An absolute r-value of greater than 0.4 was considered to indicate a forcefully significant correlation with ISG20. Gene annotation and pathway analyses were performed by Database for Annotation, Visualization and Integrated Discovery (DAVID v6.8 http://david.abcc.n-cifcrf.gov/), STRING (version 10.5 https://string-db.org/), and Kyoto Encyclopedia of Genes and Genomes (KEGG http://www. kegg.jp/kegg/pathway.html). Gene Set Enrichment Analysis (GSEA) was used to explore biological functions associated with ISG20. Gene Set Variation Analysis (GSVA package of R http://www.bioconductor.org/) was implemented to validate the relationship of ISG20 and the candidate functions. To further define the immune cell subpopulations affected by ISG20 expression, we used the GSVA to explore the relationship between ISG20 and the predefined, highly distinctive transcriptional profile of each immune cell type. [71] [72] [73] [74] Twenty-four types of immune cells with corresponding gene signatures were analyzed (Table S7) . A threshold of an absolute correlation coefficient (r) greater than 0.3 with a p-value < 0.05 was used for selecting immune cells significantly correlated with ISG20. Recurrence and malignant degeneration after resection of adult hemispheric low-grade gliomas CGCG clinical practice guidelines for the management of adult diffuse gliomas IDH1 mutations are early events in the development of astrocytomas and oligodendrogliomas Using the molecular classification of glioblastoma to inform personalized treatment A CpG island methylator phenotype in acute myeloid leukemia independent of IDH mutations and associated with a favorable outcome CpG island methylator phenotype is associated with response to adjuvant irinotecan-based therapy for stage III colon cancer IDH1 mutation is sufficient to establish the glioma hypermethylator phenotype The CpG island methylator phenotype (CIMP) in colorectal cancer Prognostic and predictive value of CpG island methylator phenotype in patients with locally advanced nonmetastatic sporadic colorectal cancer The global DNA methylation surrogate LINE-1 methylation is correlated with MGMT promoter methylation and is a better prognostic factor for glioma Breast cancer methylomes establish an epigenomic foundation for metastasis Methylation profile of triple-negative breast carcinomas Mutant IDH1 regulates the tumor-associated immune system in gliomas IDH mutant gliomas escape natural killer cell immune surveillance by downregulation of NKG2D ligand expression Isocitrate dehydrogenase mutations suppress STAT1 and CD8+ T cell accumulation in gliomas Correlation of immune phenotype with IDH mutation in diffuse glioma An integrated genomic analysis of human glioblastoma multiforme Cancer-related inflammation Glioblastoma-associated microglia and macrophages: targets for therapies to improve prognosis Macrophages in experimental and human brain tumors. Part 2: studies of the macrophage content of human brain tumors Interferon-alpha enhances CD95L-induced apoptosis of human malignant glioma cells CCNU-dependent potentiation of TRAIL/Apo2L-induced apoptosis in human glioma cells is p53-independent but may involve enhanced cytochrome c release MiR-15a/16 deficiency enhances anti-tumor immunity of glioma-infiltrating CD8+ T cells through targeting mTOR Cancer immunotherapy with recombinant poliovirus induces IFN-dominant activation of dendritic cells and tumor antigen-specific CTLs Molecular cloning of a new interferon-induced PML nuclear body-associated protein Promyelocytic leukemia (PML) nuclear bodies are protein structures that do not accumulate RNA PML bodies: a meeting place for genomic loci? The transcriptional role of PML and the nuclear body Transcriptional regulation is affected by subnuclear targeting of reporter plasmids to PML nuclear bodies The mechanisms of PML-nuclear body formation Cellular proteins PML and Daxx mediate an innate antiviral defense antagonized by the adenovirus E4 ORF3 protein Role of SUMO-1-modified PML in nuclear body formation The exonuclease ISG20 mainly localizes in the nucleolus and the Cajal (Coiled) bodies and is associated with nuclear SMN protein-containing complexes Three conserved members of the RNase D family have unique and overlapping functions in the processing of 5S, 5.8S, U4, U5, RNase MRP and RNase P RNAs in yeast ISG20, an actor of the innate immune response U3 snoRNP associates with fibrillarin a component of the scleroderma clumpy nucleolar domain Mammalian and yeast U3 snoRNPs are matured in specific and related nuclear compartments New clues to the function of the Cajal body Increased rate of degradation of c-myc mRNA in interferontreated Daudi cells RNase-L-dependent destabilization of interferoninduced mRNAs. A role for the 2-5A system in attenuation of the interferon response A unique ISRE, in the TATA-less human Isg20 promoter, confers IRF-1-mediated responsiveness to both interferon type I and type II Type I interferon-regulated gene expression and signaling in murine mixed glial cells lacking signal transducers and activators of transcription 1 or 2 or interferon regulatory factor 9 Genomewide transcriptome induced by Porphyromonas gingivalis LPS supports the notion of host-derived periodontal destruction and its association with systemic diseases Herpes virus induced proteasomedependent degradation of the nuclear bodies-associated PML and Sp100 proteins Identification and validation of genes involved in cervical tumourigenesis Tyrosine kinase chromosomal translocations mediate distinct and overlapping gene regulation events Treatment of patients with disseminated colorectal cancer with recombinant human alpha 2-interferon. Studies on the immune system Suppression of natural killer cell cytotoxicity in the peripheral blood of patients receiving interferon therapy Interferon-alpha and cancer: mechanisms of action and new perspectives of clinical use Transgenic CCL2 expression in the central nervous system results in a dysregulated immune response and enhanced lethality after coronavirus infection Selective attraction of monocytes and T lymphocytes of the memory phenotype by cytokine RANTES Sánchez-Mateos P. Migration of human blood dendritic cells across endothelial cell monolayers: adhesion molecules and chemokines involved in subset-specific transmigration Alpha and beta chemokines induce NK cell migration and enhance NK-mediated cytolysis Cytokine RANTES released by thrombin-stimulated platelets is a potent attractant for human eosinophils Tissue proteases convert CCL23 into potent monocyte chemoattractants in patients with chronic rhinosinusitis The atypical receptor CCRL2 is required for CXCR2-dependent neutrophil recruitment and tissue damage CCRL2, a fringe member of the atypical chemoattractant receptor family CXCR6 promotes atherosclerosis by supporting T-cell homing, interferon-gamma production, and macrophage accumulation in the aortic wall Identification of a unique TGF-beta-dependent molecular and functional signature in microglia PD-1 and its ligands in tolerance and immunity Comprehensive RNA-seq transcriptomic profiling in the malignant progression of gliomas Integrative genomics viewer Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes 1. Overview of the immune response 2. Cytokines and chemokines Cytokines and chemokines Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer Identification of T cell-restricted genes, and signatures for different T cell responses, using a comprehensive collection of microarray datasets Distinct transcriptional profiles in ex vivo CD4+ and CD8+ T cells are established early in human immunodeficiency virus type 1 infection and are characterized by a chronic interferon response as well as extensive transcriptional changes in CD8 Gene and protein characteristics reflect functional diversity of CD56dim and CD56bright NK cells The authors report no conflict of interest.