Abstract
The binding of platelet-derived growth factor D (PDGF-DD) to the NKp44 receptor activates a distinct transcriptional program in primary IL-2 expanded human natural killer (NK) cells. We were interested in knowing if the PDGF-DD-NKp44 pathway of NK cell activation might play a clinically relevant role in anti-tumor immunity. In order to address this question, we determined transcriptional signatures unique to resting, IL-2 expanded, and PDGF-DD activated, NK cells, in addition to different T cell subsets, and established the abundance of these immune cell phenotypes in The Cancer Genome Atlas (TCGA) low-grade glioma (LGG) dataset using CIBERSORT. Our results show that LGG patient tumors enriched for either the PDGF-DD activated NK cell or memory CD8+ T cell phenotypes are associated with a more favorable prognosis. Combined cell phenotype analyses revealed that patients with LGG tumors enriched for the PDGF-DD activated NK cell phenotype and the CD4+ T helper cell phenotype had a more favorable prognosis. High expression of transcripts encoding members of the killer cell lectin-like receptor (KLR) family, such as KLRK1 and KLRC2, KLRC3 and KLRC4 in LGG tumors were associated with more favorable prognosis, suggesting that these NK cell family receptors may play a prominent role in LGG anti-tumor immunity. Finally, many of the TCGA findings were reciprocated in LGG patients from the Chinese Glioma Genome Atlas (CGGA) dataset. Our results provide transcriptomic evidence that PDGF-DD activated NK cells and KLR family receptors may play an important clinical role in immune surveillance of LGG.
Introduction
Diffuse and infiltrative low-grade gliomas (LGGs) are derived from the malignant transformation of astrocytes or oligodendrocytes (). Whilst grade I LGGs are readily resectable benign tumors, grade II LGG display pathologic traits and inexorably progress to high grade gliomas, such as glioblastoma (GBM), with terminal neurological decline (, ). Early surgical excision and temozolomide treatment followed by radiotherapy underpin current standards of care but are not curative (–). Despite growing understanding of LGG pathogenesis, clinical outcomes have failed to improve particularly for young adults (). Furthermore, variable rates of progression to lethal disease impede timely clinical intervention and make accurate prognoses difficult (). Thus, there is an urgent need to understand effective anti-tumor immunity in LGG. Whilst the phenotype and function of tumor-infiltrating lymphocytes (TILs) have been explored for high grade gliomas (), the prognostic value of TIL subsets and the molecular pathways of tumor recognition for LGG remain unclear.
NK cells preferentially eliminate nascent tumors that have downregulated MHC class-I (MHC-I) (–). In homeostasis, NK cells are retained in the tissues surrounding the brain parenchyma by the blood-brain-barrier (BBB) (, ). Whilst NK cells have been identified in brain tumors and the surrounding tissue microenvironment (–), the mechanisms facilitating NK cell transmigration across the BBB and activation within the brain are poorly defined, although the BBB is more permeable under inflammatory conditions (–). Glioma cell lines are readily susceptible to NK cell lysis in vitro (, ). However, in vivo studies reveal a highly vascularized tumor microenvironment that actively subverts immune control (, ) and so the significance of NK cell surveillance for gliomas remains to be fully understood (). Defining immune surveillance mechanisms in those LGG patients with enhanced survival will therefore be critical for the development of novel cancer immunotherapies.
Among others, germline-encoded activating receptors, such as KLRK1 (also known as NKG2D) and the Natural cytotoxicity receptors (NCRs), such as NKp46 (NCR1), NKp44 (NCR2), and NKp30 (NCR3), can synergize to overcome inhibitory thresholds and evoke NK cell anti-tumor functions (–). As such, NK cell anti-tumor activity is sensitive to activating receptor surface phenotype and the expression of tumor ligands (). KLRK1 recognizes a range of ligands upregulated by transformed cells, such as MHC class I chain-related sequence (MIC) A and MICB, which are major determinants of NK cell tumor cytolysis in humans (–33). Of the NCRs, NKp44 also recognizes a range of cellular and tumor-associated surface ligands, such as Nidogen-1 (34), the heparan sulfate proteoglycan, Syndecan-4 (35), a subset of HLA-DP molecules (36), a splice variant of the mixed lineage leukemia 5 (MLL5) gene (37), and proliferating cell nuclear antigen (PCNA) (38), that have all been reported to positively or negatively regulate NK cell function (39, 40). Recently, PDGF-DD was shown to induce signaling from the activating NKp44 immunoreceptor (41).
The platelet-derived growth factor (PDGF) family are comprised of four polypeptides that assemble into five dimeric isoforms, PDGF-AA, PDGF-BB, PDGF-AB, PDGF-CC, and PDGF-DD. PDGFs play essential roles in embryonic development, cell proliferation, cell migration, survival and chemotaxis by engaging PDGF receptors (PDGFRs) that are mostly expressed by mesenchymal cells (42). PDGF-DD is a potent mitogen that plays an important role in wound healing and blood vessel maturation during angiogenesis by inducing PDGFR-β signaling on mesenchymal cells (43, 44). In brain cancer, PDGF-DD binding to PDGFR-β can induce pro-tumorigenic signaling that drives glioma progression (45–49).
PDGF-DD stimulation of NKp44 induced NK cell secretion of IFN-γ and TNF that arrest tumor cell proliferation in vitro and may confer a survival benefit in GBM (41). In support of this, PDGF-DD is abundantly expressed in GBM, suggesting a novel mode of NK cell tumor surveillance (50), but the clinical significance of the NKp44-PDGF-DD pathway for anti-tumor immunity in many other human cancers including LGG remains unclear.
Here, we employed a computational approach to investigate the clinical impact of the relative enrichment of resting, IL-2 expanded, and PDGF-DD activated NK cell phenotypes in the LGG tumor microenvironment. To achieve this, we used transcriptional signatures from the three NK cell activation states and estimated their relative abundance in LGG tumor specimens from TCGA database and tested the association with curated progression-free survival (51).
Methods
Material Availability
The R codes for the analyses presented in this study are available at RAGG3D/LGG_SPANK (github.com). An overview of the methods used in this study are shown in Supplementary Figure 1.
Data Collection
Gene transcript-abundance and patient clinical information were collected from TCGA through the GDC Data Portal (52) and the CGGA (53–56). Progression-free survival information was used as a measure of clinical outcome (51). The cell-type specific transcriptional signatures were derived from a large collection of RNA sequencing samples spanning a wide range of cell types. For NK cells, experimentally derived dataset for IL-2 expanded (27 biological replicates), PDGF-DD activated via NKp44 signaling (4 biological replicates), and resting (25 biological replicates from 6 datasets) (41) were included. For other cell types, the data collected was from the following datasets: BLUEPRINT (57), Monaco et al. (58), ENCODE (59), Squires et al. (60), GSE77808 (61), Tong et al. (62), PRJNA339309 (63), GSE122325 (64), FANTOM5 (65), GSE125887 (66), GSE130379 (67), GSE130286 (68).
Transcriptional Signatures
In order to derive transcriptional signatures of 21 cell types (memory B cell, naive B cell, immature dendritic myeloid cell, immature dendritic myeloid cell, endothelial, eosinophil, epithelial, fibroblast, macrophage M1 and M2, mast cell, monocyte, neutrophil, ReNK, IL2NK, SPANK, memory CD4 T cell, memory CD8 T cell, naive CD8 T cell, gamma-delta T cell and helper T cell), a total of 592 highly curated (i.e. for which identity was confirmed in the literature), non-redundant biological replicates (including 25 ReNK samples, 27 IL2NK samples and 4 SPANK samples), have been used. Due to the sparse nature of heterogeneous data sets, the expected value and variability of gene transcription abundance was inferred for each cell type using a publicly available Bayesian statistical model (github: stemangiola/cellsig), based on a negative binomial data distribution (69). This model allows to fit incomplete data (e.g. transcript abundance of one gene for which data is available in a subset of reference biological replicates) and calculate theoretical data distributions of cell-type/gene pairs. The cell-type transcriptional marker selection was based on the pairwise comparison of each cell type within cell-type categories along a cell differentiation hierarchy (Supplementary Figure 2) (70). For example, all cell-type permutations from the root node of level one (including epithelial, endothelial, fibroblasts and immune cells) were interrogated in order to select the genes for which the theoretical transcript abundance distribution (data generated from the posterior distribution) was higher for one cell type compared to another. This was executed calculating the distance of the upper and lower 95% credible intervals, respectively (obtained from cellsig). From each comparison, the top 5, 10 and 20 genes per cell-type pair were selected from levels 1, 2, and 3 (Supplementary Figure 2), and the union of all genes was taken as overall marker gene list. This hierarchical approach favors the identification of marker genes that distinguish broad cell-type categories as well as specific activation phenotypes.
Estimation of the Association of Cell-Type Abundance With Relapse-Free Patient Survival
In order to estimate the cell type relative abundance for each biological replicate, we used the algorithm CIBERSORT (71) with our RNA sequencing-derived gene marker signature. In order to estimate the clinical relevance of NK activation phenotypes (72) (73),for each cancer-type/cell-type pair, Kaplan-Meier (KM) survival curves (74) were calculated from the median split CIBERSORT-inferred proportions through the R framework tidybulk. Percent survival vs time-to-event statistics were calculated by the Log-rank (Mantel-Cox) Test (75). Statistics of KM curves were performed by log-rank test then adjusted by the Benjamini-Hochberg (BH) procedure. A table of all p-values prior to adjustment is provided in Supplementary Table 1.
Data analysis and visualization were performed using the R environment in RStudio (76). Packages include tidyverse (77), tidybulk (73), tidyHeatmap (78), survminer (79), survival (80, 81), foreach (82), org.Hs.eg.db (83), cowplot (77), ggsci (84), GGally (85), gridExtra (86), grid (76), reshape (87), Hmisc (88), and viridis (89).
Benchmark of the Transcriptional Signatures
In order to visually evaluate the ability of the marker gene selection in segregating cell types, we first performed principal component analyses (PCA) (90) for three levels of cell differentiation: (i) the NK activated states, (ii) all NK cells, and (iii) all cell types. Briefly, the raw read counts were normalized by trimmed mean of M values (TMM) using tidybulk function scale_abundance, whilst CGGA transcripts were already normalized by transcripts per kilobase of exon model per million mapped reads (TPM) via RNA-seq by expectation Maximization (RSEM). PCA analyses were performed by reduce_dimensions(method=“PCA”) (73). To directly test whether the selected signature for PDGF-DD activated NK cells was suitable to accurately detect changes in cell abundance across samples with a censored time-to-event, we implemented a test on simulated data. First, for a selected number of patients N, we sampled the progression-free survival time from the clinical annotation of the LGG patient cohort. Cell type proportions were simulated using a Dirichlet distribution, according to a linear model with a slope value S and a progression-free survival as factor of interest. For each simulated dataset, the slope S was assigned to only one cell type, and the slope of 0 to all the others (i.e. only one cell type changing for each simulated tissue mixture). The intercept (baseline proportion) was defined to be the same for all cell types. The simulated proportions were used to compose the in-silico mixtures. For each simulated dataset, the transcriptional profile of each cell type was sampled at random from the reference dataset. In order to test the accuracy of our method against the presence of foreign cell types (of which transcriptomic signature was not included in our reference set), a proportion P of neural cells was added to the mixture. The framework tidybulk was used to infer the cell type proportions through CIBERSORT and perform a multiple cox regression on the predicted proportions (logit-transformed) (81), with progression-free survival censored time as a covariate. That is, for half of the N samples, the survival time was censored to half of its value. The significance calls were compared with the ground truth to generate a receiver operating characteristic (ROC) curve. For each simulation condition (values of N, S, and P), 63 test runs were performed with one variable cell type each. A range of simulation conditions were tested, ranging N from 250 to 1000, S from 0.2 to 1, and P from 0 to 0.8.
Results
NK Cell Phenotypes Have Unique Transcriptional Profiles
Given the innate ability of NK cells to lyse tumor cells and secrete potent anti-tumor cytokines, such as IFN-γ and TNF, prior to immunization, we hypothesized that NK cells of unique phenotype may infiltrate different cancers and confer anti-tumor immunity. Specifically, we were motivated by the recent discovery of PDGF-DD as a ligand for the activating NK cell receptor NKp44 (41) and whether this mechanism of NK cell stimulation might constitute a clinically relevant pathway of anti-tumor immunity. In order to answer this question, we gathered publicly available RNA sequencing data from 21 different immune and stromal cell types (see Methods) in order to define marker genes that may distinguish resting (91), IL-2 expanded, and PDGF-DD activated, NK cell phenotypes (41).
We next performed a principal component analysis (PCA) to determine the ability of these marker genes to segregate all three NK cell phenotypes from each other and from other cell types, such as T cells (Figure 1A). NK cell activation states are associated with the first principal components (Figure 1A, left and middle panels), and NK cells overall are defined by a unique cluster when compared with other major cell types (Figure 1A, right panel). NK cell phenotypes segregated from other cell types and from each other by PCA.
Figure 1
We next performed a benchmark for the inference of changes in the relative abundance of PDGF-DD activated NK cells in association with survival information for artificial tissue mixtures built from our reference data set (see Methods). This benchmark measured the ability of the PDGF-DD activated NK cell signature to extract clinically-relevant information from TCGA whole tissue RNA sequencing data (Figure 1B). The benchmark showed a high accuracy (area under curve) across simulation settings including magnitude of variability, sample-size, and proportion of unknown cells (please see Methods). An accuracy of 0.75 (representing the area under the ROC curve) was reached for simulation settings that match our findings on TCGA data (slope and sample size; Figures 1C, D). We refer to the different NK cell phenotypes as: resting NK cells (ReNK), IL-2 expanded NK cells (IL2NK), and the signature of PDGF-DD activated NK cells (SPANK), respectively, and the transcript abundance of each marker gene in these NK cell phenotypes are shown in Supplementary Table 2.
The SPANK Is Associated With Improved Prognosis in TCGA LGG Dataset
Previous studies have implicated NK cells in immune responses to glioma (, 92–95). We next sought the association of the ReNK, IL2NK, and SPANK phenotypes with LGG. The SPANK was more abundant than either the ReNK or IL2NK phenotypes in LGG tumors from TCGA (Figure 2A). LGG tumors enriched for the SPANK were associated with greater overall patient survival compared to the ReNK or IL2NK phenotypes (Figure 2B). These results show that tumor abundance of a distinct NK cell phenotype is associated with cancer patient survival, such as LGG. Moreover, these data also suggest that LGG tumors express PDGF-DD which may activate pro-tumorigenic pathways.
Figure 2
PDGFD Expression Is Associated With LGG Invasion and Poor Prognosis
In contrast to evoking NK cell anti-tumor functions through PDGF-DD binding to NKp44 (41), PDGF-DD (encoded by PDGFD) binding to PDGFR-β (encoded by PDGFRB) induces pro-tumorigenic signaling pathways that are detrimental for cancer patient survival (42, 47, 96, 97). Expression of a three-gene signature, comprised of TGFBI, IGFBP3, and CHI3L1, has previously been associated with glioma tumor cell invasion and migration and poor patient survival (98). Tumor expression of TGFBI, IGFBP3, and CHI3L1 were positively correlated with PDGFD and PDGFRB expression in TCGA LGG dataset, respectively (Figure 3A) (98).
Figure 3
Since PDGFD and PDGFRB were associated with genes involved in glioma tumor cell invasion and migration, we next examined the relationship between tumor expression of PDGFD or PDGFRB and LGG patient survival. LGG patients with low tumor expression of PDGFD had more favorable prognosis compared to LGG patients with high tumor expression of PDGFD (Figure 3B). Higher LGG tumor expression levels of PDGFRB alone displayed a trend towards poor survival (Figure 3B) when PDGFD expression was low, but this was not statistically significant (Figure 3C). These data show that high tumor expression of PDGFD is primarily associated with poor prognosis compared to PDGFRB expression in TCGA LGG dataset.
SPANK Abundance Mitigates the Pro-Tumorigenic Effects of PDGFD in TCGA LGG
Anti-tumor immunity would be expected to curtail pro-tumorigenic factors and benefit patient survival. In addition to pro-tumor functions, we hypothesized that tumors enriched for the SPANK would contribute to anti-tumor immunity resulting in a more favorable TCGA LGG prognosis (41). In order to assess whether the abundance of these NK cell phenotypes counteracted the pro-tumorigenic expression of PDGFD and thus improve TCGA LGG prognosis, we next determined the progression-free survival of LGG patients with tumors stratified for PDGFD expression and abundance of either the ReNK, IL2NK, or SPANK phenotypes. When LGG tumors were stratified for PDGFD expression, patients with tumors enriched for the SPANK had a more favorable prognosis compared to LGG patients with a lower tumor abundance of SPANK (e.g. compare HPDGFD/HSPANK, red KM curve, to HPDGFD/LSPANK, grey KM curve and compare LPDGFD/HSPANK, yellow curve, to LPDGFD/LSPANK, blue curve) (Figure 4). In contrast, this was not observed for either ReNK or IL2NK (Figure 4). These results show that LGG tumors enriched for the SPANK may mitigate the detrimental effect of PDGFD expression on the prognosis of TCGA LGG patient cohort compared to the ReNK or IL2NK phenotypes.
Figure 4
Memory CD8+ T Cell Abundance Mitigates the Pro-Tumorigenic Effects of PDGFD in TCGA LGG
For a given cancer, it is likely that immune subsets other than NK cells infiltrate the tumor microenvironment to elicit anti-tumor immunity, particularly T cells. We were interested in knowing whether the abundance of a given T cell subset in the LGG tumor microenvironment is associated with anti-tumor immunity. TCGA LGG tumors enriched for the memory CD8+ T cell phenotype were associated with improved prognosis, but not the naïve, γδ, CD4+ memory, or Helper, T cell phenotypes (Figure 5A). Our analyses show that TCGA LGG patients enriched for the memory CD8+ T cell phenotype have improved prognosis.
Figure 5
In order to determine whether each T cell phenotype might counteract the detrimental expression of PDGFD on LGG prognosis, similarly to the SPANK, we next determined progression-free survival for LGG patients with tumors stratified for expression of PDGFD and the abundance of each T cell phenotype, respectively (Figure 5B). Using this approach, LGG tumors enriched for the memory CD8+ T cell phenotype were associated with improved prognosis (Figure 5B). These results show that LGG tumors enriched for the memory CD8+ T cell phenotype mitigate the pro-tumorigenic effects of PDGFD because they are associated with improved prognosis.
The SPANK and Helper T cell Phenotype Are Associated With Improved LGG Prognosis
Given that the abundance of NK and T cell phenotypes differ markedly in LGG tumors (Supplementary Figure 3), we were interested in understanding the relative contribution of SPANK and T cell phenotypes for LGG prognosis, respectively. We therefore determined patient survival for TCGA LGG tumors stratified for the abundance of SPANK and each respective T cell subset (Figure 5C). Interestingly, LGG tumors enriched for the SPANK and CD4+ T helper phenotypes (TH) had improved prognosis compared to other strata e.g. compare HSPANK/HTH to either HSPANK/LTH or LSPANK/HTH or LSPANK/LTH (Figure 5C, column 1). In contrast, LGG tumors enriched for the SPANK and memory CD8+ T cell phenotypes (CD8mem) did not further improve LGG patient survival compared to other strata e.g. compare HSPANK/HCD8mem to either LSPANK/HCD8mem or HSPANK/LCD8mem (Figure 5C, column 5). We conclude that LGG tumors enriched for the SPANK and CD4+ T helper cell phenotypes are associated with improved LGG prognosis in TCGA. These results provide new insights into the possible cooperation between different NK and T cell subsets for LGG anti-tumor immunity which may inform adoptive cell therapies.
Critical Role for Killer Cell Lectin-Like Receptor Family Members in LGG Anti-Tumor Immunity
NK cells express a family of germline-encoded activating and inhibitory surface receptors that engage in cancer immune surveillance, which can also be expressed by memory CD8+ T cells. However, the NK cell family receptors most critical for anti-tumor immunity in LGG remain unclear. Given that LGG patients in TCGA with tumors enriched for the SPANK or memory CD8+ T cells were associated with improved prognosis, we were interested in analyzing whether tumor expression of transcripts encoding NK cell family receptors was also associated with improved prognosis for TCGA LGG patients. LGG tumors with high expression of the KLRK1, KLRC1, KLRC2, KLRC3, or KLRC4 transcripts encoding the NKG2D, NKG2A, NKG2C, NKG2E, and NKG2F NK cell receptors, respectively, were associated with improved prognosis (Figure 6A). In contrast, high LGG tumor expression of CD226, CD244, CRTAM, KIR2DL4, NCR1, or NCR3 encoding the DNAM-1, 2B4, CRTAM, NKp46 and NKp30 NK cell receptors, respectively, were not associated with prognosis (Supplementary Figure 4). Moreover, expression of the KLRK1, KLRC1, KLRC2, KLRC3, and KLRC4 receptor genes were overwhelmingly positively correlated with the SPANK and memory CD8+ T cell phenotypes in TCGA LGG tumors (Figure 6B). These results show that high expression of the Killer cell lectin-like receptor (KLR) family in LGG tumors is associated with improved prognosis, suggesting that expression of KLR receptors may be critical for regulating anti-tumor immunity in LGG.
Figure 6
The SPANK and Memory CD8+ T Cell Phenotypes Are Associated With KLRC2 Expression and More Favorable Prognosis in CGGA LGG Patients
Since LGG tumors enriched for expression of the SPANK and KLR family receptors were associated with improved prognosis in the TCGA LGG patient cohort, we next sought to validate these findings using another glioma patient dataset, such as the CGGA (Figure 7). Similar to TCGA LGG patient cohort, high tumor expression of PDGFD in CGGA LGG patients was associated with a poor prognosis (Figure 7A). Moreover, LGG tumors enriched for the SPANK and memory CD8+ T cell phenotypes were associated with improved prognosis when CGGA tumors were also stratified for PDGFD expression compared to the IL2NK phenotype (Figure 7B). In contrast to TCGA LGG dataset, high tumor expression of KLRC1 and KLRC2 in CGGA LGG patients were associated with improved prognosis, but not KLRC3, KLRC4 or KLRK1 (Figure 7C). Like TCGA, high expression of KLRC2, which encodes the activating NKG2C receptor, was also associated with the SPANK and memory CD8+ T cell phenotypes in LGG tumors (Figure 7D). However, in contrast to TCGA, expression of KLRC1, which encodes the inhibitory NKG2A receptor, was associated with the ReNK phenotype and not the SPANK or memory CD8+ T cell phenotypes in CGGA LGG tumors (Figure 7D). Similar to TCGA, these results show that high tumor expression of PDGFD is associated with poor CGGA LGG prognosis, and tumors enriched for the SPANK and memory CD8+ T cell phenotypes have improved prognosis. Moreover, like TCGA, high LGG tumor expression of KLRC2 is also associated with the SPANK and memory CD8+ T cell phenotypes and improved prognosis of CGGA LGG patients.
Figure 7
Discussion
The clinical relevance of NK cells in cancer immune surveillance, particularly for solid tumors, remains unclear. We hypothesized that differential tumor enrichment of NK cells in different activation states may contribute to anti-tumor immunity. To investigate this important question, we determined TS from experimental RNA-seq datasets derived from NK cells in three different functional states; ReNK, IL2NK, and SPANK. Using this unbiased approach, we found that enrichment of the SPANK phenotype in LGG tumors was associated with improved prognosis in TCGA and CGGA datasets. The SPANK was derived from RNA-seq data from NK cells that had been stimulated with PDGF-DD, suggesting that PDGF-DD is expressed in the LGG tumor microenvironment.
In addition to activating NK cells, PDGF-DD binding to PDGFR-β can induce pro-tumorigenic signaling pathways. We reasoned that high expression of the genes for PDGF-D and PDGFR-β in LGG might predict poor LGG cancer prognosis (42, 47, 96, 97). In support of this, PDGFD and PDGFRB were positively correlated with a three-gene signature (TGFBI, IGFBP3, and CHI3L1) associated with glioma tumor cell invasion and migration and poor patient survival (98). Indeed, our analysis shows that LGG patients with high tumor expression of PDGFD had a poor prognosis in both TCGA and CGGA cohorts (42, 47, 96, 97). This model awaits confirmation in mouse models of glioma to determine whether PDGF-DD/PDGFR-β signaling can be targeted to restrict glioma tumor cell migration and invasion or even progression to higher glioma grades.
Given that the expression of PDGFD was primarily associated with pro-tumor pathways and poor prognosis, we further hypothesized that enrichment of the SPANK phenotype in LGG tumors may be associated with effective anti-tumor immunity and improved prognosis. High tumor abundance of the SPANK was associated with improved prognosis when LGG tumors were stratified for the expression of PDGFD in both TCGA and CGGA patient cohorts. Our data show that the relative abundance of the SPANK may counteract the pro-tumorigenic properties of PDGFD expression and improve LGG prognosis. Interestingly, Nidogen-1 (34), the heparan sulfate proteoglycan, Syndecan-4 (35), a subset of HLA-DP molecules (36), a splice variant of the MLL5 gene (37), and PCNA (38), have all been reported to bind and regulate NKp44 signaling, and it will be interesting to determine whether expression of these latter genes in the LGG tumor microenvironment can influence the association between the NK cell phenotypes that we describe and LGG prognosis.
It is very likely that NK cells are not the sole mediators of anti-tumor immunity in vivo and our analyses revealed that enrichment of the CD8+ memory T cell phenotype in LGG tumors was also associated with improved prognosis when LGG tumors were stratified for high or low tumor expression of PDGFD in both TCGA and CGGA patient cohorts. Interestingly, stratifying LGG tumors for T cell subsets and NK cells phenotypes revealed new insights into potential cooperation between these innate and adaptive immune cell subsets that was not revealed from the analysis of these immune cell phenotypes alone. For example, LGG patients with tumors enriched for the SPANK and CD4+ T helper phenotypes had improved survival suggesting that adoptive transfer of NK cells with CD4+ T helper cells may represent a novel therapeutic approach for LGG. Again, these computational results await confirmation in pre-clinical mouse models of glioma to determine whether the adoptive transfer of NK cells with CD4+ T helper cells can restrict glioma tumor cell migration and invasion or possibly even progression to higher glioma grades.
Finally, high tumor expression of the NK cell receptor genes; KLRK1, KLRC1, KLRC2, KLRC3, and KLRC4 that encode the NKG2D, NKG2A, NKG2C, NKG2E and NKG2F, respectively, were all associated with improved prognosis and positively correlated with the SPANK and memory CD8+ T cell phenotypes, suggesting expression of these KLR NK cell receptor family gene products are important for LGG anti-tumor immunity in TCGA patient cohort.
Interestingly, the KLRC1 receptor, also known as NKG2A, is expressed as a heterodimer with CD94 on the surface of NK cells and T cells (99). CD94/NKG2A can bind to HLA-E as ligand to negatively regulate signaling from other activating KLR family members including KLRC2, known as NKG2C, which also heterodimerizes with CD94 to bind HLA-E (100, 101). Interestingly, KLRC1 has recently been shown to function as a checkpoint inhibitor that when blocked can promote NK cell and CD8+ T cell-mediated anti-tumor immunity (102, 103). Since NK cell effector function is regulated by the balance of signaling from an array of germline-encoded activating and inhibitory receptors (, 104), it is possible that signaling from the activating KLRK1 (NKG2D), NCR2 (NKp44) and KLRC2 (NKG2C) receptors, which are all associated with improved survival (Figure 6A), may cooperate to overcome any inhibitory threshold set by KLRC1 (NKG2A) in TCGA LGG patient cohort (103).
High tumor expression of KLRC1 and KLRC2 were also associated with improved survival in the CGGA LGG patient cohort, but not KLRK1, KLRC2 or KLRC4. However, the SPANK and CD8+ T cell phenotypes were associated with expression of KLRC2 and not KLRC1 in the CGGA LGG patients, suggesting the balance of signaling may favor KLRC2 activation in the CGGA LGG patient cohort compared to TCGA LGG cohort. Given our results and those from other laboratories, it will be interesting to determine the expression of HLA-E in LGG and KLR family receptors on glioma-infiltrating NK cells and CD8+ T cells in different ethnic groups and to test whether blocking the inhibitory function of KLRC1 can enhance the anti-tumor activity of NK cells and CD8+ T cells in LGG and other brain cancers (102, 103). Finally, using CIBERSORT, we have uncovered an intriguing association between tumor expression of PDGFD and tumor enrichment of the SPANK and T helper and memory CD8+ T cell signatures, that may be important for LGG patient survival. However, it is premature to conclude that the SPANK or the T cell signatures play a definitive role in LGG survival and future studies will aim to determine the biological significance of the SPANK and different T cell phenotypes and NK cell receptors, such as the KLR family, in LGG patient survival.
Funding
This work was funded by a MRFF research acceleration grant APP1162217 awarded to AB, and a University of Melbourne PhD scholarship awarded to YS. SM is funded by a Galli Early Career Research Fellowship.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
AB, SM, YS, and YP contributed to the conception and design of the study. YS and SM performed computational and statistical analysis. AB and YS wrote the first draft of the manuscript. SM and AS wrote sections of the manuscript. AB and SM contributed equally and are joint senior authors. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.668391/full#supplementary-material
Supplementary Figure 1Computational pipeline. In addition to the three NK cell functional phenotypes, bulk RNA-seq data was collected from open online resources for 21 immune and stromal cell types including five T cell subsets, two B cell subsets, two dendritic cells subsets, two macrophage subsets, eosinophils, mast cells and neutrophils, fibroblasts, endothelial and epithelial cells. Secondly, we selected marker genes by pairwise comparison and CIBERSORT to create our transcriptional signature (TS). Thirdly, we input our TS and RNA-seq matrix of TCGA cancer patients into CIBERSORT to estimate the proportion of each cell type in each patient. Finally, we performed a series of statistical analysis using the immune cell type profiles and clinical results of all TCGA cancer patients.
Supplementary Figure 2The Cell differentiation hierarchy used in marker gene selection. The higher levels are differentiated from cells in the lower levels (L1 differentiates into L2, L2 differentiates into L3, and L3 differentiates into L4).
Supplementary Figure 3Profiles of NK and T cell subset phenotypes in LGG patients. (A) Profiles showing original fractions of NK and T cell phenotypes in LGG tumors (B) Proportions of total NK and T cell abundance in LGG tumors.
Supplementary Figure 4Kaplan-Meier curves of NK cell receptor family genes in LGG patients. Survival analysis of receptor genes CD226 (encoding DNAM-1), CD244 (encoding 2B4/SLAMF4), CRTAM, KIR2DL4, NCR1 (encoding NKp46) and NCR3 (encoding NKp30) in TCGA-LGG. High expression of CD226 gene is significantly related to poor survival. Low (L) and high (H) expression groups were split by median TMM normalized counts. Statistics was performed by log-rank test.
Supplementary Table 1Log rank p values and Benjamini-Hochberg corrected p values. Log rank p-values of Kaplan-Meier curves were all corrected by BH procedure (/ = unchanged).
Supplementary Table 2NK phenotype TS transcripts. The abundance of each marker gene in the TS of NK cell phenotypes.
References
1
LouisDNPerryAReifenbergerGVon DeimlingAFigarella-BrangerDCaveneeWKet al. The 2016 World Health Organization Classification of Tumors of the Central Nervous System: A Summary. Acta Neuropathologica (2016) 131(6):803–20. doi: 10.1007/s00401-016-1545-1
2
WesselsPHWeberWEJRavenGRamaekersFCSHopmanAHNTwijnstraA. Supratentorial Grade II Astrocytoma: Biological Features and Clinical Course. Lancet Neurol (2003) 2(7):395–403. doi: 10.1016/S1474-4422(03)00434-4
3
ClausEBHorlacherAHsuLSchwartzRBDello-IaconoDTalosFet al. Survival Rates in Patients With Low-Grade Glioma After Intraoperative Magnetic Resonance Image Guidance. Cancer: Interdiscip Int J Am Cancer Soc (2005) 103(6):1227–33. doi: 10.1002/cncr.20867
4
JakolaASMyrmelKSKlosterR. Comparison of a Strategy Favoring Early Surgical Resection vs a Strategy Favoring Watchful Waiting in Low-Grade Gliomas. JAMA (2012) 25(October):1–8. doi: 10.1001/jama.2012.12807
5
SanaiNBergerMS. Glioma Extent of Resection and its Impact on Patient Outcome. Neurosurgery (2008) 62(4):753–66. doi: 10.1227/01.neu.0000318159.21731.cf
6
DuffauHTaillandierL. New Concepts in the Management of Diffuse Low-Grade Glioma: Proposal of a Multistage and Individualized Therapeutic Approach. Neuro-oncology (2015) 17(3):332–42. doi: 10.1093/neuonc/nou153
7
ClausEBWalshKMWienckeJKMolinaroAMWiemelsJLSchildkrautJMet al. Survival and Low-Grade Glioma: The Emergence of Genetic Information. Neurosurg Focus (2015) 38(1):E6. doi: 10.3171/2014.10.FOCUS12367
8
KimY-HJungT-YJungSJangW-YMoonK-SKimI-Yet al. Tumour-Infiltrating T-Cell Subpopulations in Glioblastomas. Br J Neurosurg (2012) 26(1):21–7. doi: 10.3109/02688697.2011.584986
9
KarlhoferFMRibaudoRKYokoyamaWM. MHC Class I Alloantigen Specificity of Ly-49+ IL-2-Activated Natural Killer Cells. Nature (1992) 358(6381):66–70. doi: 10.1038/358066a0
10
ColonnaMSamaridisJ. Cloning of Immunoglobulin-Superfamily Members Associated With HLA-C and HLA-B Recognition by Human Natural Killer Cells. Science (1995) 268(5209):405–8. doi: 10.1126/science.7716543
11
MalnatiMSPeruzziMParkerKCBiddisonWECicconeEMorettaAet al. Peptide Specificity in the Recognition of MHC Class I by Natural Killer Cell Clones. Science (1995) 267(5200):1016–8. doi: 10.1126/science.7863326
12
KorinBBen-ShaananTLSchillerMDubovikTAzulay-DebbyHBoshnakNTet al. High-Dimensional, Single-Cell Characterization of the Brain’s Immune Compartment. Nat Neurosci (2017) 20(9):1300. doi: 10.1038/nn.4610
13
EngelhardtBVajkoczyPWellerRO. The Movers and Shapers in Immune Privilege of the CNS. Nat Immunol (2017) 18(2):123. doi: 10.1038/ni.3666
14
DominguesPGonzález-TablasMOteroÁPascualDMirandaDRuizLet al. Tumor Infiltrating Immune Cells in Gliomas and Meningiomas. Brain Behav Immun (2016) 53:1–15. doi: 10.1016/j.bbi.2015.07.019
15
RenFZhaoQHuangLZhengYLiLHeQet al. The R132H Mutation in IDH 1 Promotes the Recruitment of NK Cells Through CX 3 CL 1/CX 3 CR 1 Chemotaxis and is Correlated With a Better Prognosis in Gliomas. Immunol Cell Biol (2019) 97(5):457–69. doi: 10.1111/imcb.12225
16
HollEKFrazierVLandaKBeasleyGHwangESNairSK. Examining Peripheral and Tumor Cellular Immunome in Patients With Cancer. Front Immunol (2019) 10:1767. doi: 10.3389/fimmu.2019.01767
17
LepennetierGHracskoZUngerMVan GriensvenMGrummelVKrumbholzMet al. Cytokine and Immune Cell Profiling in the Cerebrospinal Fluid of Patients With Neuro-Inflammatory Diseases. J Neuroinflamm (2019) 16(1):219. doi: 10.1186/s12974-019-1601-6
18
KastrukoffLFLauASTakeiFCarboneFRScalzoAA. A NK Complex-Linked Locus Restricts the Spread of Herpes Simplex Virus Type 1 in the Brains of C57BL/6 Mice. Immunol Cell Biol (2015) 93(10):877–84. doi: 10.1038/icb.2015.54
19
OwensGCGarciaAJMochizukiAYChangJWReyesSDSalamonNet al. Evidence for Innate and Adaptive Immune Responses in a Cohort of Intractable Pediatric Epilepsy Surgery Patients. Front Immunol (2019) 10:121. doi: 10.3389/fimmu.2019.00121
20
ZhangYGaoZWangDZhangTSunBMuLet al. Accumulation of Natural Killer Cells in Ischemic Brain Tissues and the Chemotactic Effect of IP-10. J Neuroinflamm (2014) 11(1):79. doi: 10.1186/1742-2094-11-79
21
SivoriSParoliniSMarcenaroECastriconiRPendeDMilloRet al. Involvement of Natural Cytotoxicity Receptors in Human Natural Killer Cell-Mediated Lysis of Neuroblastoma and Glioblastoma Cell Lines. J Neuroimmunol (2000) 107(2):220–5. doi: 10.1016/S0165-5728(00)00221-6
22
CastriconiRDagaADonderoAZonaGPolianiPLMelottiAet al. NK Cells Recognize and Kill Human Glioblastoma Cells With Stem Cell-Like Properties. J Immunol (2009) 182(6):3530–9. doi: 10.4049/jimmunol.0802845
23
CraneCAHanSJBarryJJAhnBJLanierLLet al. TGF-β Downregulates the Activating Receptor NKG2D on NK Cells and CD8+ T Cells in Glioma Patients. Neuro-oncology (2010) 12(1):7–13. doi: 10.1093/neuonc/nop009
24
AppolloniIAlessandriniFCeresaDMarubbiDGambiniEReverberiDet al. Progression From Low- to High-Grade in a Glioblastoma Model Reveals the Pivotal Role of Immunoediting. Cancer Lett (2019) 442:213–21. doi: 10.1016/j.canlet.2018.10.006
25
SedgwickAJGhazanfariNConstantinescuPMantamadiotisTBarrowAD. The Role of NK Cells and Innate Lymphoid Cells in Brain Cancer. Front Immunol (2020) 11. doi: 10.3389/fimmu.2020.01549
26
PendeDCantoniCRiveraPVitaleMCastriconiRMarcenaroSet al. Role of NKG2D in Tumor Cell Lysis Mediated by Human NK Cells: Cooperation With Natural Cytotoxicity Receptors and Capability of Recognizing Tumors of Nonepithelial Origin. Eur J Immunol (2001) 31(4):1076–86. doi: 10.1002/1521-4141(200104)31:4<1076::AID-IMMU1076>3.0.CO;2-Y
27
EagleRAJafferjiIBarrowAD. Beyond Stressed Self: Evidence for NKG2D Ligand Expression on Healthy Cells. Curr Immunol Rev (2009) 5(1):22–34. doi: 10.2174/157339509787314369
28
BrycesonYTMarchMELjunggrenH-GLongEO. Synergy Among Receptors on Resting NK Cells for the Activation of Natural Cytotoxicity and Cytokine Secretion. Blood (2006) 107(1):159–66. doi: 10.1182/blood-2005-04-1351
29
LanierLL. Up on the Tightrope: Natural Killer Cell Activation and Inhibition. Nat Immunol (2008) 9(5):495–502. doi: 10.1038/ni1581
30
BauerSGrohVWuJSteinleAPhillipsJHLanierLLet al. Activation of NK Cells and T Cells by NKG2D, a Receptor for Stress-Inducible MICA. Science (1999) 285(5428):727–9. doi: 10.1126/science.285.5428.727
31
CosmanDMüllbergJSutherlandCLChinWArmitageRFanslowWet al. ULBPs, Novel MHC Class I–related Molecules, Bind to CMV Glycoprotein UL16 and Stimulate NK Cytotoxicity Through the NKG2D Receptor. Immunity (2001) 14(2):123–33. doi: 10.1016/S1074-7613(01)00095-4
32
EagleRAFlackGWarfordAMartínez-BorraJJafferjiITraherneJAet al. Cellular Expression, Trafficking, and Function of Two Isoforms of Human ULBP5/RAET1G. PloS One (2009) 4(2):e4503. doi: 10.1371/journal.pone.0004503
33
CerwenkaABakkerABHMcClanahanTWagnerJWuJPhillipsJHet al. Retinoic Acid Early Inducible Genes Define a Ligand Family for the Activating NKG2D Receptor in Mice. Immunity (2000) 12(6):721–7. doi: 10.1016/S1074-7613(00)80222-8
34
GaggeroSBruschiMPetrettoAParodiMZottoGDLavarelloCet al. Nidogen-1 is a Novel Extracellular Ligand for the NKp44 Activating Receptor. OncoImmunology (2018) 7(9):e1470730. doi: 10.1080/2162402X.2018.1470730
35
BrusilovskyMRadinskyOCohenLYossefRShemeshABraimanAet al. Regulation of Natural Cytotoxicity Receptors by Heparan Sulfate Proteoglycans in -Cis: A Lesson From Nkp44. Eur J Immunol (2015) 45(4):1180–91. doi: 10.1002/eji.201445177
36
NiehrsAGarcia-BeltranWFNormanPJWatsonGMHölzemerAChapelAet al. A Subset of HLA-DP Molecules Serve as Ligands for the Natural Cytotoxicity Receptor Nkp44. Nat Immunol (2019) 20(9):1129–37. doi: 10.1038/s41590-019-0448-4
37
VieillardVBaychelierFDebréP. NKp44L: A New Tool for Fighting Cancer. Oncoimmunology (2014) 3(1):e27988. doi: 10.4161/onci.27988
38
KunduKGhoshSSarkarREdriABrusilovskyMGershoni-YahalomOet al. Inhibition of the NKp44-PCNA Immune Checkpoint Using a mAb to PCNA. Cancer Immunol Res (2019) 7(7):1120–34. doi: 10.1158/2326-6066.CIR-19-0023
39
BarrowADMartinCJColonnaM. The Natural Cytotoxicity Receptors in Health and Disease. Front Immunol (2019) 10:909. doi: 10.3389/fimmu.2019.00909
40
ParodiMFavoreelHCandianoGGaggeroSSivoriSMingariMCet al. NKp44-NKp44 Ligand Interactions in the Regulation of Natural Killer Cells and Other Innate Lymphoid Cells in Humans. Front Immunol (2019) 10:719. doi: 10.3389/fimmu.2019.00719
41
BarrowADEdelingMA. Natural Killer Cells Control Tumor Growth by Sensing a Growth Factor. Cell (2018) 172(3):534–548. e19. doi: 10.1016/j.cell.2017.11.037
42
AndraeJGalliniRBetsholtzC. Role of Platelet-Derived Growth Factors in Physiology and Medicine. Genes Dev (2008) 22(10):1276–312. doi: 10.1101/gad.1653708
43
MaJZhaoYSunLSunXZhaoXSunXet al. Exosomes Derived From Akt-Modified Human Umbilical Cord Mesenchymal Stem Cells Improve Cardiac Regeneration and Promote Angiogenesis via Activating Platelet-Derived Growth Factor D. Stem Cells Transl Med (2017) 6(1):51–9. doi: 10.5966/sctm.2016-0038
44
UutelaMWirzeniusMPaavonenKRajantieIHeYKarpanenTet al. PDGF-D Induces Macrophage Recruitment, Increased Interstitial Pressure, and Blood Vessel Maturation During Angiogenesis. Blood (2004) 104(10):3198–204. doi: 10.1182/blood-2004-04-1485
45
HermansonMFunaKHartmanMClaesson-WelshLHeldinC-HWestermarkBet al. Platelet-Derived Growth Factor and its Receptors in Human Glioma Tissue: Expression of Messenger RNA and Protein Suggests the Presence of Autocrine and Paracrine Loops. Cancer Res (1992) 52(11):3213–9.
46
PlateKHScholzADumontDJ. Tumor Angiogenesis and Anti-Angiogenic Therapy in Malignant Gliomas Revisited. Acta Neuropathologica (2012) 124(6):763–75. doi: 10.1007/s00401-012-1066-5
47
LokkerNASullivanCMHollenbachSJIsraelMAGieseNA. Platelet-Derived Growth Factor (PDGF) Autocrine Signaling Regulates Survival and Mitogenic Pathways in Glioblastoma Cells: Evidence That the Novel PDGF-C and PDGF-D Ligands may Play a Role in the Development of Brain Tumors. Cancer Res (2002) 62(13):3729–35.
48
KrenzlinHBeheraPLorenzVPassaroCZdiorukMNowickiMOet al. Cytomegalovirus Promotes Murine Glioblastoma Growth via Pericyte Recruitment and Angiogenesis. J Clin Invest (2019) 129(4):1671–83. doi: 10.1172/JCI123375
49
LaRochelleWJJeffersMMcDonaldWFChillakuruRAGieseNALokkerNAet al. PDGF-D, a New Protease-Activated Growth Factor. Nat Cell Biol (2001) 3:517–21. doi: 10.1038/35074593
50
TrinchieriG. Natural Killer Cells Detect a Tumor-Produced Growth Factor: A Vestige of Antiviral Resistance? Trends Immunol (2018) 39(5):357–8. doi: 10.1016/j.it.2018.02.002
51
LiuJLichtenbergTHoadleyKAPoissonLMLazarAJCherniackADet al. An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell (2018) 173(2):400–416.e11. doi: 10.1016/j.cell.2018.02.052
52
TomczakKCzerwinskaPWiznerowiczM. The Cancer Genome Atlas (TCGA): An Immeasurable Source of Knowledge. Contemp Oncol (Pozn) (2015) 19(1A):A68–77. doi: 10.5114/wo.2014.47136
53
LiuXLiYQianZSunZXuKWangKet al. A Radiomic Signature as a non-Invasive Predictor of Progression-Free Survival in Patients With Lower-Grade Gliomas. NeuroImage: Clin (2018) 20:1070–7. doi: 10.1016/j.nicl.2018.10.014
54
ZhaoZMengFWangWWangZZhangCJiangT. Comprehensive RNA-Seq Transcriptomic Profiling in the Malignant Progression of Gliomas. Sci Data (2017) 4(1):170024. doi: 10.1038/sdata.2017.24
55
BaoZSChenHMYangMYZhangCBYuKYeWLet al. RNA-Seq of 272 Gliomas Revealed a Novel, Recurrent PTPRZ1-MET Fusion Transcript in Secondary Glioblastomas. Genome Res (2014) 24(11):1765–73. doi: 10.1101/gr.165126.113
56
WangYQianTYouGPengXChenCYouYet al. Localizing Seizure-Susceptible Brain Regions Associated With Low-Grade Gliomas Using Voxel-Based Lesion-Symptom Mapping. Neuro-Oncology (2014) 17(2):282–8. doi: 10.1093/neuonc/nou130
57
FernandezJMde la TorreVRichardsonDRoyoRPuiggrosMMoncunillVet al. The BLUEPRINT Data Analysis Portal. Cell Syst (2016) 3(5):491–495 e5. doi: 10.1016/j.cels.2016.10.021
58
MonacoGLeeBXuWMustafahSHwangYYCarreCet al. RNA-Seq Signatures Normalized by mRNA Abundance Allow Absolute Deconvolution of Human Immune Cell Types. Cell Rep (2019) 26(6):1627–1640 e7. doi: 10.1016/j.celrep.2019.01.041
59
SloanCAChanETDavidsonJMMalladiVSStrattanJSHitzBCet al. ENCODE Data at the ENCODE Portal. Nucleic Acids Res (2016) 44(D1):D726–32. doi: 10.1093/nar/gkv1160
60
SquiresRBNoronhaJHuntVGarcia-SastreAMackenCBaumgarthNet al. Influenza Research Database: An Integrated Bioinformatics Resource for Influenza Research and Surveillance. Influenza Other Respir Viruses (2012) 6(6):404–16. doi: 10.1111/j.1750-2659.2011.00331.x
61
MaoYvan HoefVZhangXWennerbergELorentJWittKet al. IL-15 Activates mTOR and Primes Stress-Activated Gene Expression Leading to Prolonged Antitumor Capacity of NK Cells. Blood (2016) 128(11):1475–89. doi: 10.1182/blood-2016-02-698027
62
TongPLRoedigerBKolesnikoffNBiroMTaySSJainRet al. The Skin Immune Atlas: Three-Dimensional Analysis of Cutaneous Leukocyte Subsets by Multiphoton Microscopy. J Invest Dermatol (2015) 135(1):84–93. doi: 10.1038/jid.2014.289
63
FerraroNMDampierWWeingartenMSSpillerKL. Deconvolution of Heterogeneous Wound Tissue Samples Into Relative Macrophage Phenotype Composition via Models Based on Gene Expression. Integr Biol (Camb) (2017) 9(4):328–38. doi: 10.1039/c7ib00018a
64
WangYGellatlyKMcCauleySVangalaPKimKDerrAet al. HIV-1-Induced Cytokines Deplete Homeostatic ILCs and Expand TCF7-Dependent Memory NK Cells. bioRxiv (2019) p:221010. doi: 10.1038/s41590-020-0593-9
65
AnderssonRGebhardCMiguel-EscaladaIHoofIBornholdtJBoydMet al. An Atlas of Active Enhancers Across Human Cell Types and Tissues. Nature (2014) 507(7493):455–61. doi: 10.1038/nature12787
66
CildirGToubiaJYipKHZhouMPantHHissariaPet al. Genome-Wide Analyses of Chromatin State in Human Mast Cells Reveal Molecular Drivers and Mediators of Allergic and Inflammatory Diseases. Immunity (2019) 51(5):949–965 e6. doi: 10.1016/j.immuni.2019.09.021
67
MarquardtNKekalainenEChenPLourdaMWilsonJNScharenbergMet al. Unique Transcriptional and Protein-Expression Signature in Human Lung Tissue-Resident NK Cells. Nat Commun (2019) 10(1):3841. doi: 10.1038/s41467-019-11632-9
68
SabryMZubiakAHoodSPSimmondsPArellano-BallesteroHCournoyerEet al. Tumor- and Cytokine-Primed Human Natural Killer Cells Exhibit Distinct Phenotypic and Transcriptional Signatures. PloS One (2019) 14(6):e0218674. doi: 10.1371/journal.pone.0218674
69
PillowJWScottJG. Fully Bayesian Inference for Neural Models With Negative-Binomial Spiking. Adv Neural Inf Process Syst (2012) 3:1898–906. doi: 10.5555/2999325.2999347
70
MangiolaSStuchberyRMcCoyPChowKKurganovsNKergerMet al. Androgen Deprivation Therapy Promotes an Obesity-Like Microenvironment in Periprostatic Fat. Endocr Connect (2019) 8(5):547–58. doi: 10.1530/EC-19-0029
71
ChenBKhodadoustMSLiuCLNewmanAMAlizadehAA. Profiling Tumor Infiltrating Immune Cells With CIBERSORT. Methods Mol Biol (2018) 1711:243–59. doi: 10.1007/978-1-4939-7493-1_12
72
LundeATimmermannABlakeD. The Hazards of Mutual Fund Underperformance: A Cox Regression Analysis. J Empirical Finance (1999) 6(2):121–52. doi: 10.1016/S0927-5398(98)00013-9
73
MangiolaSMolaniaRDongRDoyleMAPapenfussAT. Tidybulk: An R Tidy Framework for Modular Transcriptomic Data Analysis. Genome Biol (2021) 22(1):42. doi: 10.1186/s13059-020-02233-7
74
BlandJMAltmanDG. Survival Probabilities (the Kaplan-Meier Method). BMJ (1998) 317(7172):1572–80. doi: 10.1136/bmj.317.7172.1572
75
MantelN. Evaluation of Survival Data and Two New Rank Order Statistics Arising in its Consideration. Cancer Chemother Rep (1966) 50(3):163–70.
76
TeamRC. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing (2020).
77
WickhamHAverickMBryanJChangWMcGowanLDAFrançoisRet al. Welcome to the Tidyverse. J Open Source Softw (2019) 4(43):1686. doi: 10.21105/joss.01686
78
MangiolaSPapenfussAT. Tidyheatmap: An R Package for Modular Heatmap Production Based on Tidy Principles. J Open Source Softw (2020) 5(52):2472. doi: 10.21105/joss.02472
79
KassambaraAKosinskiMBiecekP. Survminer: Drawing Survival Curves Using ‘Ggplot2’. (2020). Available at: https://CRAN.R-project.org/package=survminer.
80
TherneauTM. A Package for Survival Analysis in R. (2020). Available at: https://CRAN.R-project.org/package=survival.
81
TherneauTMGrambschPM. Modeling Survival Data: Extending the Cox Model. New York: Springer (2000).
82
MicrosoftWestonS. Foreach: Provides Foreach Looping Construct. (2020). Available at: https://CRAN.R-project.org/package=foreach.
83
CarlsonM. Org.Hs.eg.db: Genome Wide Annotation for Human. (2019). Available at: https://bioconductor.org/packages/release/data/annotation/html/org.Hs.eg.db.html.
84
XiaoN. Ggsci: Scientific Journal and Sci-Fi Themed Color Palettes for ‘Ggplot2’. (2018). Available at: https://CRAN.R-project.org/package=ggsci.
85
SchloerkeBCookDLarmarangeJBriatteFMarbachMThoenEet al. GGally: Extension to ‘Ggplot2’. (2020). Available at: https://CRAN.R-project.org/package=GGally.
86
AuguieB. Gridextra: Miscellaneous Functions for “Grid” Graphics. (2017). Available at: https://CRAN.R-project.org/package=gridExtra.
87
WickhamH. Reshaping Data With the Reshape Package. J Stat Softw (2007) 21(12). doi: 10.18637/jss.v021.i12
88
HarrellFEJrwith contributions from Charles Dupont and many others. Hmisc: Harrell Miscellaneous. (2020). Available at: https://CRAN.R-project.org/package=Hmisc.
89
GarnierS. Viridis: Default Color Maps From ‘Matplotlib’. (2018). Available at: https://CRAN.R-project.org/package=viridis.
90
WoldSEsbensenKGeladiP. Principal Component Analysis. Chemometrics Intelligent Lab Syst (1987) 2(1):37–52. doi: 10.1016/0169-7439(87)80084-9
91
CrinierAMilpiedPEscalièreBPiperoglouCGallusoJBalsamoAet al. High-Dimensional Single-Cell Analysis Identifies Organ-Specific Signatures and Conserved NK Cell Subsets in Humans and Mice. Immunity (2018) 49(5):971–86. doi: 10.1016/j.immuni.2018.09.009
92
Dominguez-ValentinMGras NavarroARahmanAMKumarSRetièreCUlvestadEet al. Identification of a Natural Killer Cell Receptor Allele That Prolongs Survival of Cytomegalovirus-Positive Glioblastoma Patients. Cancer Res (2016) 76(18):5326–36. doi: 10.1158/0008-5472.CAN-16-1162
93
VauléonETonyAHamlatAEtcheverryAChiforeanuDCMeneiPet al. Immune Genes are Associated With Human Glioblastoma Pathology and Patient Survival. BMC Med Genomics (2012) 5:41. doi: 10.1186/1755-8794-5-41
94
QiuHLiYChengSLiJHeCLiJ. A Prognostic Microenvironment-Related Immune Signature via ESTIMATE (PROMISE Model) Predicts Overall Survival of Patients With Glioma. Front Oncol (2020) 10:580263. doi: 10.3389/fonc.2020.580263
95
CloseHJSteadLFNsengimanaJReillyKADroopAWurdakHet al. Expression Profiling of Single Cells and Patient Cohorts Identifies Multiple Immunosuppressive Pathways and an Altered NK Cell Phenotype in Glioblastoma. Clin Exp Immunol (2020) 200(1):33–44. doi: 10.1111/cei.13403
96
BergstenEUutelaMLiXPietrasKÖstmanAHeldinC-Het al. PDGF-D is a Specific, Protease-Activated Ligand for the PDGF β-Receptor. Nat Cell Biol (2001) 3(5):512–6. doi: 10.1038/35074588
97
GialeliCNikitovicDKletsasDD TheocharisAN TzanakakisGKaramanosN. PDGF/PDGFR Signaling and Targeting in Cancer Growth and Progression: Focus on Tumor Microenvironment and Cancer-Associated Fibroblasts. Curr Pharm Design (2014) 20(17):2843–8. doi: 10.2174/13816128113199990592
98
AtkinsRJStylliSSKurganovsNMangiolaSNowellCJWareTMet al. Cell Quiescence Correlates With Enhanced Glioblastoma Cell Invasion and Cytotoxic Resistance. Exp Cell Res (2019) 374(2):353–64. doi: 10.1016/j.yexcr.2018.12.010
99
SullivanLCNguyenTHOHarpurCMStankovicSKanagarajahARKoutsakosMet al. Natural Killer Cell Receptors Regulate Responses of HLA-E-Restricted T Cells. Sci Immunol (2021) 6(58). doi: 10.1126/sciimmunol.abe9057
100
TomasecPBraudVMRickardsCPowellMBMcSharryBPGadolaSet al. Surface Expression of HLA-E, an Inhibitor of Natural Killer Cells, Enhanced by Human Cytomegalovirus Gpul40. Science (2000) 287(5455):1031. doi: 10.1126/science.287.5455.1031
101
HammerQRückertTBorstEMDunstJHaubnerADurekPet al. Peptide-Specific Recognition of Human Cytomegalovirus Strains Controls Adaptive Natural Killer Cells. Nat Immunol (2018) 19(5):453–63. doi: 10.1038/s41590-018-0082-6
102
AndréPDenisCSoulasCBourbon-CailletCLopezJArnouxTet al. Anti-NKG2A mAb is a Checkpoint Inhibitor That Promotes Anti-Tumor Immunity by Unleashing Both T and NK Cells. Cell (2018) 175(7):1731–43. doi: 10.1016/j.cell.2018.10.014
103
van MontfoortNBorstLKorrerMJSluijterMMarijtKASantegoetsSJet al. NKG2A Blockade Potentiates CD8 T Cell Immunity Induced by Cancer Vaccines. Cell (2018) 175(7):1744–1755.e15. doi: 10.1016/j.cell.2018.10.028
104
TrowsdaleJJonesDCBarrowADTraherneJA. Surveillance of Cell and Tissue Perturbation by Receptors in the LRC. Immunol Rev (2015) 267(1):117–36. doi: 10.1111/imr.12314
Summary
Keywords
NK cell, low grade glioma, NK cell receptor, anti-tumor immunity, The Cancer Genome Atlas
Citation
Sun Y, Sedgwick AJ, Palarasah Y, Mangiola S and Barrow AD (2021) A Transcriptional Signature of PDGF-DD Activated Natural Killer Cells Predicts More Favorable Prognosis in Low-Grade Glioma. Front. Immunol. 12:668391. doi: 10.3389/fimmu.2021.668391
Received
18 February 2021
Accepted
10 August 2021
Published
02 September 2021
Volume
12 - 2021
Edited by
Paul J. Norman, University of Colorado Denver, United States
Reviewed by
Maria Graciela Castro, University of Michigan, United States; Francesco Colucci, University of Cambridge, United Kingdom; Pedro Lowenstein, University of Michigan, United States
Updates
Copyright
© 2021 Sun, Sedgwick, Palarasah, Mangiola and Barrow.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Stefano Mangiola, mangiola.s@wehi.edu.au; Alexander David Barrow, alexanderdav@unimelb.adu.au
This article was submitted to NK and Innate Lymphoid Cell Biology, a section of the journal Frontiers in Immunology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.