Impact Factor 6.429

The 5th most cited journal in Immunology

Original Research ARTICLE

Front. Immunol., 07 November 2017 | https://doi.org/10.3389/fimmu.2017.01451

CD4+ and Perivascular Foxp3+ T Cells in Glioma Correlate with Angiogenesis and Tumor Progression

imageLuyan Mu1,2, imageChanglin Yang3, imageQiang Gao4, imageYu Long2,3, imageHaitao Ge2, imageGabriel DeLeon3†, imageLinchun Jin2,3, imageYifan (Emily) Chang3†, imageElias J. Sayour3, imageJingjing Ji5, imageJie Jiang6, imagePaul S. Kubilis3, imageJiping Qi6, imageYunhe Gu6, imageJiabin Wang2, imageYuwen Song1, imageDuane A. Mitchell3, imageZhiguo Lin2* and imageJianping Huang3*
  • 1Department of Neurosurgery, The Fourth Affiliated Hospital, Harbin Medical University, Harbin, China
  • 2The Fourth Section of Department of Neurosurgery, The First Affiliated Hospital, Harbin Medical University, Harbin, China
  • 3Lillian S. Wells Department of Neurosurgery, University of Florida, Gainesville, FL, United States
  • 4Department of Geriatrics, The Second Affiliated Hospital of Harbin Medical University, Harbin, China
  • 5Department of Pathology, The Second Affiliated Hospital of Harbin Medical University, Harbin, China
  • 6Department of Pathology, The First Affiliated Hospital, Harbin Medical University, Harbin, China

Background: Angiogenesis and immune cell infiltration are key features of gliomas and their manipulation of the microenvironment, but their prognostic significance remains indeterminate. We evaluate the interconnection between tumor-infiltrating lymphocyte (TIL) and tumor blood-vasculatures in the context of glioma progression.

Methods: Paired tumor tissues of 44 patients from three tumor-recurrent groups: diffuse astrocytomas (DA) recurred as DA, DA recurred as glioblastomas (GBM), and GBM recurred as GBM were evaluated by genetic analysis, immunohistochemistry for tumor blood vessel density, TIL subsets, and clinical outcomes. These cells were geographically divided into perivascular and intratumoral TILs. Associations were examined between these TILs, CD34+ tumor blood vessels, and clinical outcomes. To determine key changes in TIL subsets, microarray data of 15-paired tumors from patients who failed antiangiogenic therapy- bevacizumab, and 16-paired tumors from chemo-naïve recurrent GBM were also evaluated and compared.

Results: Upon recurrence in primary gliomas, similar kinetic changes were found between tumor blood vessels and each TIL subset in all groups, but only CD4+ including Foxp3+ TILs, positively correlated with the density of tumor blood vessels. CD4 was the predominant T cell population based on the expression of gene-transcripts in primary GBMs, and increased activated CD4+ T cells were revealed in Bevacizumab-resistant recurrent tumors (not in chemo-naïve recurrent tumors). Among these TILs, 2/3 of them were found in the perivascular niche; Foxp3+ T cells in these niches not only correlated with the tumor vessels but were also an independent predictor of shortened recurrence-free survival (RFS) (HR = 4.199, 95% CI 1.522–11.584, p = 0.006).

Conclusion: The minimal intratumoral T cell infiltration and low detection of CD8 transcripts expression in primary GBMs can potentially limit antitumor response. CD4+ and perivascular Foxp3+ TILs associate with tumor angiogenesis and tumor progression in glioma patients. Our results suggest that combining antiangiogenic agents with immunotherapeutic approaches may help improve the antitumor efficacy for patients with malignant gliomas.

Introduction

Tumors recur in the resection margin in nearly 90% of CNS malignancies after primary surgery and adjuvant therapies (1, 2). Understanding which significant alterations occur in the tumor immune microenvironment during progression is critical for designing effective immunotherapeutic strategies. Tumor angiogenesis and immune cell infiltrations are crucial tumor-driven processes in tumors (3, 4); however, mechanistic insights regarding their interplay, as well as their prognostic significance in glioma recurrence remain indeterminate.

Astrocytic gliomas arise from astrocytes and are the most common type of glioma, representing 64% of human CNS malignant tumors (5). Nearly all the patients with primary gliomas will face tumor recurrence thus necessitating a better understanding of glioma progression. Recently, the angiogenesis inhibitor Bevacizumab has been shown to prolong recurrence-free survival (RFS), although no increase in overall survival (OS) was seen, in patients with newly diagnosed or recurrent GBM (68). This drug inhibits the formation of tumor blood vessels induced by vascular endothelial growth factor A (VEGF-A); increases T cell infiltrations (9, 10); reverses expression of inhibitory molecules associated with T cell exhaustion (11); and may directly mediate antitumor effect (12). These results imply that angiogenesis plays an important role in GBM progression, but the failure of these agents to improve patients’ OS also suggest that negative feedback pathways (such as hypoxia) may be spontaneously activated resulting in increased tumor cell invasion (13, 14).

Tumors can orchestrate complex biological networks via angiogenesis and recruitment of regulatory immune cell subsets. The VEGF/receptor axis has been shown to have strong immune regulatory properties (3, 911). Under hypoxia condition, tumor-infiltrating lymphocytes (TILs) were found to express VEGF-A (15). Evidence indicates that the tumor blood vessel/tumor endothelium can be a substantial barrier for TIL to extravasate into the intratumoral space stymieing their ability to mount a strong antitumor response (16). Distinctly separating TILs into perivascular and intratumoral populations will provide better insight for us to determine the interaction between tumor angiogenesis and T cell infiltration.

Although several studies have previously evaluated lymphocyte infiltrations in glioma patients (1722), those studies have focused on non-paired patient samples, which could potentially introduce significant variations. Since gliomas are heterogeneous, district genetic or phenotypic characteristics have been observed in cells from the same tumor (23, 24), which may be conserved in recurrence. Distinctions between primary and recurrent tumors in the same patient may enable a better understanding of transition/transformation during tumor progression. In this study, we studied 44 paired tumor samples (before and after of progression) in three categories of recurrence from astrocytic glioma patients and probed the interconnection between differential geographic T cell subsets and tumor blood vessels. Additionally, we also analyzed 15-paired GBMs samples from patients who had failed bevacizumab therapy and determined distinct TIL changes within these tumors. Our data suggest that CD4+ and perivascular Foxp3+ TILs impact angiogenesis and tumor recurrence in patients with gliomas.

Materials and Methods

Patient Population

This study included 44 patients who underwent surgery for primary and recurrent glioma by MRI imaging at the same hospital from 2005 to 2014 (Table S2 in Supplementary Material). The average period between diagnosis and surgeries was 5.4 days (range 2–10 days). The inclusion criteria included the following: greater than 18 years of age at the time of primary diagnosis, a histopathological diagnosis of diffuse astrocytomas (DA) or GBM, and the availability of formalin-fixed, paraffin-embedded (FFPE) tumor tissue blocks. The patients were separated into three groups: (1) the DA–DA (DA, diffuse astrocytoma) group (DAs that recur as DAs, 15 pairs); (2) the GBM–GBM group (GBMs that recur as GBMs, 15 pairs); and (3) the DA–GBM group (DAs that recur as GBMs, 14 pairs). None of the 44 patients received treatments before their primary surgery. Subsequently, patients received standard radiotherapy (60 Gy), while some received concomitant–adjuvant chemotherapy (temozolomide/other). The recurrence-free survival (RFS) was defined as the time from the date of the primary surgery to the first MRI-confirmed tumor recurrence. All recurrent patients underwent recurrent surgery. OS was defined as the interval between primary surgery and death or last follow-up. Survival after recurrence (SR) was defined from the recurrent surgery to the date of death or the last follow-up. The mean follow-up period was 1,355.1 days (range 289–3,520 days), during which 36 patients died, and 8 patients remained alive. No patients were lost to follow-up. Tumor size was calculated based on enhanced MR or CT images (on the layer with the maximum amount of tumor) as follows: long diameter (centimeter) × wide diameter (centimeter) × thickness (centimeter) × 0.5 (25). All specimens were selected after histopathological review to ensure that we would have sufficient tissue for analysis. Treatment information was obtained from the time of primary surgery, including any adjuvant treatment. The latest updated OS information was obtained on February 13, 2016. For patients who were still alive, the date of the last follow-up was used in the data analysis. All subjects gave written informed consent in accordance with the Declaration of Helsinki and research protocols were reviewed and approved by the Committee on Human Research at the Harbin Medical University, China.

Microarray Analysis of Gene Expression

To generate high-quality RNA from paired FFPE samples, a Sensation Plus™ FFPE Amplification and WT Labeling Kit was used to extract and amplify total RNA that was derived from whole tumor FFPE samples. The total RNA obtained from each sample was quantified using a NanoDrop ND-1000, and we selected three of the highest-quality paired samples from each group. The samples were used for labeling and submitted for hybridization array scanning using an Affymetrix Gene Chip® Human Transcriptome Array 2.0. The microarray was processed, and the data were analyzed by Beijing Biolancet Technology, Co., Ltd. Nine pairs of samples were processed by Affymetrix HTA2.0 kit. After sequencing and QC, raw microarray data (CEL) were analyzed by Affymetrix Expression Console 1.4 software. HTA-2_0 library and annotation files were downloaded from Expression Console’s build-in database. Raw data were further normalized by the RMA-algorithm (Table S3 in Supplementary Material). Differentially expressed genes between samples were identified using fold-change filtering to identify distinct gene expression profiles between samples. For gene enrichment analysis, paired recurrent and primary gliomas samples’ microarray data were normalized by the log2 algorithm of genes mean value. Then, Gene Set Enrichment Analysis (GSEA) (26, 27) was applied to analyze the normalized data by groups. Enrichment Score, Normalized Enrichment Score, False Discovery Rate, and Nominal p-Value were reported by GSEA. Interpreting GSEA Results of the above four key statistics can be found in http://software.broadinstitute.org/gsea/doc/GSEAUserGuideFrame.html.

Histological Analysis, Immunohistochemistry, and Immunofluorescence

A systematic neuro-pathological review was performed based on the 2007 World Health Organization classification guidelines for CNS tumors (28). Paired samples of slides were reviewed to determine their histological classifications. Immunohistochemical staining was performed for CD3 (clone LN10, dilution 1:200, Quanhui, China), CD4 (clone UMAB64, dilution 1:50; ZSGB-BIO, China), CD8 (clone SP16, dilution 1:100; ZSGB-BIO, China), Foxp3 (clone mAbcam 450, dilution 1:50; Abcam), and CD34 (clone EP88, dilution 1:150; ZSGB-BIO, China). The secondary antibodies came as an immunohistochemical kit (KIT-5930, Maxim, China) and were incubated for 40 min at room temperature. Each staining batch included a negative control that was processed without the primary antibody, a biological negative control that consisted of normal brain tissues (the donor had died of myocardial infarction), and a positive control (tonsil tissues).

The IHC analyses were performed using a quantitative approach under a light microscope (Leica SP2, Leica Optical Co. Ltd., Germany). Subpopulations of TILs expressing CD3, CD4, Foxp3, or CD8 were divided into intratumoral fractions, which contained no vessels in each high powered field (HPF) (40× objective and 10× eyepiece), and perivascular fraction, which contained more than one vessel in each HPF (29, 30). The mean counts of CD34+ vascular circles were obtained from 10 consecutive HPF by three pathologists blinded to outcome data. Similarly, for TIL frequencies, the results were determined from a mean of 10 consecutive HPFs of typical tumor regions. The average scores counted by three experienced pathologists blinded to the clinical background were recorded as the final result. When there was a large difference in scores between observers, the score was re-evaluated to reach an agreement. Results from CD3+, CD4+, FoxP3+, CD8+ cell, and CD34+ vascular circle counts in each patient were used in the statistical analysis.

The following antibodies were used for two-color immunofluorescence: rabbit anti-CD34 and mouse anti-CD3 monoclonal antibodies (as described above). The primary antibodies were incubated overnight at 4°C. The sections were then incubated with goat anti-mouse fluorescein-isothiocyanate-conjugated secondary antibodies and goat anti-rabbit rhodamine-conjugated secondary antibodies for 2 h at room temperature. The nuclei were stained with DAPI. Finally, the sections were visualized using a confocal microscope (as described above).

Statistical Analysis

The raw data were tabulated using Microsoft Excel (Microsoft Corp.) and analyzed using nonparametric tests and matched-paired-tests (Wilcoxon tests or Mann–Whitney U tests). Chi-square tests were used to identify differences in chemotherapy and sex. Spearman’s rank order correlation analyses were performed to detect significant associations with positive marker expression. A Kaplan–Meier survival analysis was used to evaluate differences in RFS, SR, and OS. To adjust for potential confounders, Cox proportional hazards models were used to evaluate hazard ratios (HRs) for recurrence or death according to the number of identified TIL subpopulations and clinical features. All statistical analyses were performed using SPSS 13.0 (SPSS Inc. Chicago, IL, USA) and GraphPad Prism 6.0 (GraphPad Software Inc., La Jolla, CA, USA). All tests used to determine the level of significance were two-sided. A (two-tailed) p-value threshold of 0.05 was considered to indicate statistical significance (*p < 0.05, **p < 0.01, ***p < 0.001).

Results

Similar Kinetic Changes between Tumor-Associated Blood Vessels and TIL Subsets upon Recurrence

To learn more about overall genetic changes before and after recurrence across the three recurrent groups (DA–DA, DA–GBM, and GBM-GBM), three pairs of FFPE glioma samples from each group were analyzed by Microarray. A distinct pattern of gene expression was observed before and after each recurrence (Figure S1A in Supplementary Material). Among these altered genes upon recurrence, VEGF-A was found increased by 10.2-fold after tumor recurrence in the DA–GBM recurrent group (Table S1 in Supplementary Material). Distinct clinical outcomes were also revealed in the three groups; the longest RFS was the DA to GBM recurrence (Table S2 and Figure S1B in Supplementary Material). Based on these results, blood vessels were measured using a vascular endothelial cell marker, CD34. Tumor-infiltrating T cells, such as CD3+, CD4+, and CD8+ T cells were determined by immunohistochemistry (Figure 1A; Figure S2A–C in Supplementary Material). Compared to other two groups, the primary tumors in the GBM-GBM group showed relatively higher on the average density of CD34+ circles and T cells, while the DA–GBM group displayed a wider-range of T cell infiltrations (Figures 1B–E). Similar kinetics of CD34+ circles and infiltrating T cells were observed when these parameters were compared before and after tumor recurrence/progression; compared to recurrent GBM, significantly increased CD34+ circles and T cell subsets were found in secondary GBM (Figures 1F–I).

FIGURE 1
www.frontiersin.org

Figure 1. A similar kinetic movement of tumor angiogenesis and T cell infiltrations upon tumor recurrence [e.g., diffuse astrocytomas (DA)–DA, 15 pairs, DA–GBM, 15 pairs and GBM-GBM, 14 pairs]. (A) A representative fluorescent immunohistochemistry staining of tumor blood vessels measured by CD34 staining (green, red arrow) and tumor-infiltrating CD3+ T cells (red, white arrow) in a surgical resected primary GBM sample. (B–E) The baseline counts of CD34+ circles, CD3+, CD4+, and CD8+ cells in tumors. Counts were carried out for average numbers [10 consecutive high powered fields (HPFs)] of these marker expressing cells in primary tumors among three groups. (F) Enhancement of CD34+ circles after the recurrence. A calculation was carried out by comparing the increase in CD34+ circles between recurrent and primary tumors. (G–I) A similar measurement was also carried out for the average numbers of infiltrating CD3+, CD4+, and CD8+ T cells in these patients’ tumors, respectively. The significance between two groups was measured using Mann–Whitney U test, *p < 0.05.

CD4+ T Cells Are Associated with Tumor Blood Vessels

To test which infiltrating T cell subset associated most closely with tumor angiogenesis, we next did a Spearman’s rank correlation analysis between the density of tumor vessels and T cell infiltrating subsets. CD4+ T cells correlated strongly with CD34+ circles in both primary and recurrent tumors across all three groups; notwithstanding primary tumors of DA–GBM, there was no association with CD8+ T cells (Figures 2A–C).

FIGURE 2
www.frontiersin.org

Figure 2. A correlation was found between the density of CD34 circles and CD4+, but not for CD8+ T cells, in primary and recurrent tumors. The correlation CD4+ or CD8+ T cell counts with CD34 circles [average cell or vascular circle counts were obtained from 10 consecutive high powered fields (HPFs)] in primary or recurrent tumors in diffuse astrocytomas (DA)–DA group (A), GBM–GBM group (B), and DA–GBM group (C). Spearman’s rank order correlation analyses were performed and graphs with a correlation (Spearman’s rho ≥ 0.5, p < 0.05) are indicated.

CD4 Transcript Expression Is Abundant in Primary GBMs, and Activated CD4+ T Cells Are Enriched in Bevacizumab Resistant Tumors

We found that CD4 transcript was predominantly expressed in primary GBM (>29-fold higher than CD8 alpha and beta chains), based on the RNA-seq analysis solely culled from TCGA database (TCGA Research Network: http://cancergenome.nih.gov/, Figure 3A). No similar results were observed for some few other cancer types, such as kidney renal clear cell, testicular germ cell tumor using TCGA Pan-Cancer datasets (30 different cancer types) (Figure 3B). To confirm our observation that infiltrating CD4+ T cells may be associated with tumor angiogenesis, microarray data of paired bevacizumab resistant tumors from previously published reports were analyzed and paired chemo-naïve tumors were used as a comparison control (31, 32). Only signals of activated CD4+ T cells, not other cell populations, were significantly increased in tumors after a recurrence of bevacizumab-treated patients (Figure 3C); no such trend found in the chemo-naïve recurrent patient tumors (Figure 3D).

FIGURE 3
www.frontiersin.org

Figure 3. CD4 full-length transcript expression is abundant in primary GBM and enriched in bevacizumab-resistant tumors. (A) CD4 was a predominant gene transcript product (measured by RNA-seq by Expectation-Maximization, RSEM) of T cells in primary GBM compared with transcripts of CD8 alpha (CD8A) and beta (CD8B). (B) Gene profile of CD4, CD8A, and CD8B in other tumors, including GBM. RNAseq data of TCGA pan-cancer gene expression were used, and CD4, CD8A, and CD8B gene expressions from primary tumors were clustered as Heatmap using Subio Platform after logarithmic transformation and normalization. (C) Activated CD4+ T cells (highlighted) were significantly enriched in recurrent tumors after Bevacizumab treatment had failed in GBM patients. Paired surgical tumor samples before and after Bevacizumab treatment derived from 15 patients were analyzed by microarray (32). (D) No CD4 enrichment was in paired GBMs without any treatment. Microarray data from 16 pairs of tumor samples from recurrent GBM without any therapy were analyzed (33). The raw data were enriched between baseline and recurrent tumors using Gene Set Enrichment Analysis (GSEA) software (Broad Institute). The classification of the cell populations was based on the T cell-related metagenes (26).

Increasing Levels of Perivascular TILs after Glioma Recurrence

T cell tumor trafficking and tumor angiogenesis have been shown to be linked biological phenomena (3). We found both factors had an upwards trend upon recurrence; we next separated tumor-infiltrating T cells into perivascular and intratumoral CD3+, CD4+, CD8+, and Foxp3+ T cell subsets (Figure S2A–C in Supplementary Material). We demonstrated that there were more perivascular than intratumoral TILs in primary and recurrent tumors. In primary tumors, only one-third of the TILs were located in the intratumoral space, and two-thirds of the TILs were in the perivascular niche; these findings reproduce previously published reports (30). When comparing the cell numbers before and after recurrence, the perivascular T cells were enriched after tumor progression across all three groups; the intratumoral T cells, however, were only increased in recurrent tumors of the DA–GBM group (Figures 4A–F).

FIGURE 4
www.frontiersin.org

Figure 4. A significant increase of perivascular tumor-infiltrating lymphocytes after tumor recurrence. Perivascular and intratumoral infiltrating CD3+, CD4+, CD8+, and Foxp3+ T cells before and after recurrence were evaluated from randomly selected consecutive 10 high powered fields (HPFs). (A,B) diffuse astrocytomas (DA)–DA, (C,D) DA–GBM, and (E,F) GBM–GBM. Same scales are plotted for all the recurrent groups. The significances were determined respectively using paired t-test.

Perivascular Foxp3+ Tumor-Infiltrating T Cells Associate with Angiogenesis and Is an Independent Factor That Predicts Glioma Progression/Recurrence

The association between tumor-infiltrating CD4+ T cells and tumor angiogenesis led us to further test the association between angiogenesis and a subset of immune inhibitory CD4+ T cells or regulatory T cells (Tregs). We found that perivascular Foxp3+ (Figures 5A–B upper panel), and CD4+ T cells (Figures 5A–B, lower panel) in the primary tumors correlated with tumor blood vessels. No similar trend was found for the CD8+ T cells in both geographic locations (data not shown). We sought to confirm the observation by using gene enrichment for Treg signature genes using paired samples. The results indicate that tumor recurrence significantly increases CD4 (similar to the results from bevacizumab resistant tumors) and Tregs signals in these tumors, whereas other cell subsets, i.e., NK, Th17, show contrasting results (Figure 5C and data not shown). To evaluate the correlation between differential geographic T cell subsets and clinical outcomes (RFS, OS), Kaplan–Meier survival curves were plotted using these tested parameters against RFS and OS. We found that only perivascular CD4+ T cells and Tregs were associated with shorter RFS (p = 0.007 for CD34, p = 0.01 for CD4+ T cells and p = 0.001 for Foxp3+ T cells). No association was observed for these factors with respect to OS (Figures 5D–E). Additionally, there was no association found for perivascular or intratumoral CD8+ T cells with respect to RFS and OS (Figure S3 in Supplementary Material). Next, to assess if key T cell subsets independently predict glioma progression, Cox multivariate regression analyses were performed (by including factors, such as age, sex, chemo/radiotherapy, CD34+ circles and perivascular/intratumoral infiltrating T cell subsets). The results concluded that only perivascular Foxp3+ T cells were found to be an independent predictor of shortened RFS (HR = 4.199, 95% CI 1.522–11.584, p = 0.006) when all the patients were included.

FIGURE 5
www.frontiersin.org

Figure 5. Perivascular Foxp3+ T cells associated with angiogenesis and tumor progression. (A,B) Perivascular Foxp3+ (upper panel) or CD4+ T cells (lower panel) were positively associated with CD34+ circles, respectively. (C) Regulatory T cells (Tregs) signature genes enriched by the tumor recurrence. The microarray data described in Figure S1 in Supplementary Material were input into Gene Set Enrichment Analysis (GSEA) for an analysis of activated CD4 T cells and Tregs. (D) Association of perivascular or intratumoral CD4+ T cells with RFS (left) and overall survival (OS) (right). (E) Association of intratumoral or perivascular Foxp3+ T cells with RFS and OS. The association between the markers and RFS, OS is shown using Kaplan–Meier survival curves, the median value was used for stratifying the low and high in each recurrent group. The log-rank test was used to compare the differences.

Discussion

Several clinical and experimental studies have evaluated the correlation between TILs and survival in patients with gliomas. Some reports indicate a positive correlation between the abundance of CD3+ or CD8+ TILs and patients survival (19, 20, 34, 35), whereas others report conflicting observations (36). For Tregs, no consensus has been reached for patients with gliomas (22, 35, 37). Nearly all the reports were based on analysis from un-paired patient populations. In this study, we collected 44 paired patient tumor samples and compared changes in gene expression before and after tumor recurrence. Compared to the other two groups, the most significant change in gene expression after recurrence was seen in the DA–GBM group (secondary GBM). Intriguingly, this group had the longest RFS, wider ranges of T cell infiltration in primary tumors, the greatest change in tumor blood vessel and TIL subsets post-recurrence, and the strongest correlation with the tumor vascularization, suggesting that disease progression in this group is unique. The kinetics of CD34+ vascular circles and the abundance of T cell subsets after recurrence showed an identical trend which led us to hypothesize that an interconnection exists between neovascularization and TILs. Indeed, we found a positive correlation between CD4+ TILs and the CD34+ circles, but not for CD8+ or CD3+ T cells (CD3+ T cell data not shown). To figure out why only CD4+ T cells correlated, we analyzed RNA-seq data of 155 primary GBMs culled from TCGA. In these tumors, CD4 expression was approximately 29-fold higher than CD8 transcripts (CD8A and B). This difference is not expected to be due to an error from sample preparation and analysis because the same analysis of other cancer types demonstrated higher expression of CD8 transcripts then CD4. The observation led us to speculate that the CD8+ T cells seeing in tumor tissues by IHC may be inadequate for the RNA-seq detection. Additionally, results from our recent study indicate that selective apoptosis of CD8+ T cells occurs in GBMs (38). Next, we utilized public data from Bevacizumab-treated paired patient samples; 15 patients had microarray data for their tumors before and after the bevacizumab treatment (31, 32). The results demonstrated that genes of activated CD4+ T cells were significantly increased in tumors after the recurrence, whereas other cell populations presented less or no change. These results confirmed our observation that CD4+ T cells are preferentially linked with tumor angiogenesis. Moreover, the gene profiles obtained by RNA-seq analysis of isolated CD4+ and CD8+ T cells also suggests that the activated CD4+ T cells express more pro-angiogenic molecules such as a CD40 ligand and Aquaporin 3 (39, 40). The evidence from previous reports supports our hypothesis. In an ischemia animal model, CD4-deficient mice have impaired capability to undergo angiogenesis (41). Although we were unable to be sure that FoxP3+ cells were CD4+ T cells by IHC staining (due to technical difficulties), we predict these cells to be Tregs based on: (1) morphology; (2) geographic location (i.e., perivascular association); (3) tumor immunosuppressive microenvironment limiting passage of activated T cells; and (4) gene enrichment studies demonstrating enhanced activated CD4+ cells and Tregs upon tumor recurrence. Ongoing studies in animal models have linked glioma angiogenesis and Treg involvement; combining VEGF and CD25 blockades restores the IFN-γ production by T cells previously suppressed by gliomas and significantly prolongs the OS in mice compared with single drug treatment (Long et al., manuscript in preparation). These data suggest that CD4+ T cells, including Tregs, may be involved in glioma angiogenesis.

T cell infiltration into the Virchow–Robin space is distinct from T cell infiltration into intratumoral spaces (42). Moreover, intratumoral TILs appear to associate with more favorable clinical outcomes (43). In gliomas, however, more investigations are necessary to draw definitive conclusions (30, 44). Our results showed that only one-third of all TILs are located in the intratumoral space while two-thirds of them (presumably CD4+ T cells, including Tregs) surrounding the tumor vessels. We found that Tregs highly correlated with the density of blood vessels in primary tumors, but only in the perivascular zone. Importantly, these perivascular Tregs cells were identified as an independent risk factor for tumor recurrence.

Recent genetic classification schemas have advanced our understanding of molecular characteristics and factors that impact tumor progression and survival (4547). The status of IDH, TERT promoter mutations, and the deletion of chromosome 1p/19q was found to significantly impact patient clinical outcome (45, 48, 49). Beyond the scope of this report, we also examined our patients’ IDH mutation status; among the 44 patients, 22 tumors harbored the mutation, while the other 22 were IDH wild type (mostly from the GBM-GBM group). We uncovered that IDH mutations impact on the tumor immune landscape, and affect survival outcomes (Mu et al., manuscript under review). We included all patient data in Figure 5, and when patients in the GBM-GBM groups were removed, the main conclusion of Foxp3+ T cells as the independent risk factor for tumor recurrence remained true. Presumably, these Foxp3+ T cells are CD4 positive since it was the predominant T cell expressed transcript observed in GBM, and we also have found that CD8+ T cells are apoptotic in GBM (50). Thus, these CD4+ Foxp3+ T cells not only play the key role in pro-immunosuppression but also possess the pro-angiogenic function of the CD4+ T cells. The dual effects of these cells in primary tumors make them a strong player in the promotion of tumor progression in juxtaposition with the extremely low expression of CD8 transcripts in primary GBMs, which can be a major obstacle in tumor treatment.

In summary, only one-third of TILs were found in the intratumoral space with minimal expression of CD8 transcripts in primary tumors, thus limiting the overall strength of the antitumor response. The predominant population of CD4+ T cells may promote tumor angiogenesis, and in conjunction with perivascular CD4+ Tregs predispose tumor recurrence/progression in patients with gliomas.

Author Contributions

Conception and design: LM, ZL, and JH. Development of methodology: LM, CY, QG, YL, HG, YC, LJ, JQ, JJi, JJiang, YG, JW, and YS. Data analysis: LM, CY, PK, and JH. Analysis and interpretation of data: LM, CY, ES, ZL, and JH. Writing, review, and revision of the manuscript: LM, YC, GL, ES, DM, ZL, and JH. Study supervision: JH.

Conflict of Interest Statement

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.

Funding

This research was supported by Ministry of Science and Technology of the People’s Republic of China (2014DFA31630), the Specialized Research Fund for the Doctoral Program of Higher Education (SRFDP) (No. 20132307110029), The Applied Technology Research and Development Program of Heilongjiang Province (No. GA15C108), the National Natural Science Foundation of China (NSFC) (No. 81571646), the Heilongjiang Natural Science Foundation (H2016027, H2016034), the Health and Family Planning Commission of Heilongjiang Province (No. 2014-313), Harbin medical university scientific research innovation fund(2016LCZX06), and Wells Endowment (00107592).

Supplementary Material

The Supplementary Material for this article can be found online at http://www.frontiersin.org/article/10.3389/fimmu.2017.01451/full#supplementary-material.

Figure S1. Global gene profile changes and clinical outcomes for patients with DA to DA (DA-DA); DA to GBM (DA-GBM) and GBM to GBM (GBM-GBM) recurrence.

Figure S2. Representative images of CD34+ circles and intratumoral and perivascular CD3+, CD4+, CD8+, and Foxp3+ TIL cell subsets in 3 recurrent groups.

Figure S3. No association between CD8+ T cells with RFP in primary tumors.

Table S1. Upregulated genes after DA-GBM.

Table S2. Patients’ demographic and clinical characteristics.

Table S3. Microarray analysis of paired surgical resected tumors (3 pairs DA to DA, 3 pairs GBM to GBM and 3 pairs DA to GBM).

References

1. Lemée JM, Clavreul A, Menei P. Intratumoral heterogeneity in glioblastoma: don’t forget the peritumoral brain zone. Neuro Oncol (2015) 17(10):1322–32. doi:10.1093/neuonc/nov119

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Hamard L, Ratel D, Selek L, Berger F, van der Sanden B, Wion D. The brain tissue response to surgical injury and its possible contribution to glioma recurrence. J Neurooncol (2016) 128:1–8. doi:10.1007/s11060-016-2096-y

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Motz GT, Coukos G. The parallel lives of angiogenesis and immunosuppression: cancer and other tales. Nat Rev Immunol (2011) 11:702–11. doi:10.1038/nri3064

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Das S, Marsden PA. Angiogenesis in glioblastoma. N Engl J Med (2013) 369:1561–3. doi:10.1056/NEJMcibr1309402

CrossRef Full Text | Google Scholar

5. Ostrom QT, Gittleman H, Fulop J, Liu M, Blanda R, Kromer C, et al. CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2008-2012. Neuro Oncol (2015) 17:iv1–62. doi:10.1093/neuonc/nov189

CrossRef Full Text | Google Scholar

6. Gilbert MR, Dignam JJ, Armstrong TS, Wefel JS, Blumenthal DT, Vogelbaum MA, et al. A randomized trial of bevacizumab for newly diagnosed glioblastoma. N Engl J Med (2014) 370:699–708. doi:10.1056/NEJMoa1308573

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Chinot OL, Wick W, Mason W, Henriksson R, Saran F, Nishikawa R, et al. Bevacizumab plus radiotherapy–temozolomide for newly diagnosed glioblastoma. N Engl J Med (2014) 370:709–22. doi:10.1056/NEJMoa1308345

CrossRef Full Text | Google Scholar

8. Friedman HS, Prados MD, Wen PY, Mikkelsen T, Schiff D, Abrey LE, et al. Bevacizumab alone and in combination with irinotecan in recurrent glioblastoma. J Clin Oncol (2009) 27:4733–40. doi:10.1200/JCO.2008.19.8721

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Shrimali RK, Yu Z, Theoret MR, Chinnasamy D, Restifo NP, Rosenberg SA. Antiangiogenic agents can increase lymphocyte infiltration into tumor and enhance the effectiveness of adoptive immunotherapy of cancer. Cancer Res (2010) 70:6171–80. doi:10.1158/0008-5472.CAN-10-0153

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Hamzah J, Jugold M, Kiessling F, Rigby P, Manzur M, Marti HH, et al. Vascular normalization in Rgs5-deficient tumours promotes immune destruction. Nature (2008) 453:410–4. doi:10.1038/nature06868

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Voron T, Colussi O, Marcheteau E, Pernot S, Nizard M, Pointet A-L, et al. VEGF-A modulates expression of inhibitory checkpoints on CD8+ T cells in tumors. J Exp Med (2015) 212:139–48. doi:10.1084/jem.20140559

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Kim KJ, Li B, Winer J, Armanini M, Gillett N, Phillips HS, et al. Inhibition of vascular endothelial growth factor-induced angiogenesis suppresses tumour growth in vivo. Nature (1993) 362:841–4. doi:10.1038/362841a0

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Casanovas O, Hicklin DJ, Bergers G, Hanahan D. Drug resistance by evasion of antiangiogenic targeting of VEGF signaling in late-stage pancreatic islet tumors. Cancer Cell (2005) 8:299–309. doi:10.1016/j.ccr.2005.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Keunen O, Johansson M, Oudin A, Sanzey M, Rahim SAA, Fack F, et al. Anti-VEGF treatment reduces blood supply and increases tumor cell invasion in glioblastoma. Proc Natl Acad Sci U S A (2011) 108:3749–54. doi:10.1073/pnas.1014480108

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Freeman MR, Schneck FX, Gagnon ML, Corless C, Soker S, Niknejad K, et al. Peripheral blood T lymphocytes and lymphocytes infiltrating human cancers express vascular endothelial growth factor: a potential role for T cells in angiogenesis. Cancer Res (1995) 55:4140–5.

PubMed Abstract | Google Scholar

16. Motz GT, Santoro SP, Wang L-P, Garrabrant T, Lastra RR, Hagemann IS, et al. Tumor endothelium FasL establishes a selective immune barrier promoting tolerance in tumors. Nat Med (2014) 20:607–15. doi:10.1038/nm.3541

CrossRef Full Text | Google Scholar

17. Yu JS, Lee PK, Ehtesham M, Samoto K, Black KL, Wheeler CJ. Intratumoral T cell subset ratios and Fas ligand expression on brain tumor endothelium. J Neurooncol (2003) 64:55–61. doi:10.1007/BF02700020

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Heimberger AB, Abou-Ghazal M, Reina-Ortiz C, Yang DS, Sun W, Qiao W, et al. Incidence and prognostic impact of FoxP3+ regulatory T cells in human gliomas. Clin Cancer Res (2008) 14:5166–72. doi:10.1158/1078-0432.CCR-08-0320

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yang I, Tihan T, Han SJ, Wrensch MR, Wiencke J, Sughrue ME, et al. CD8+ T-cell infiltrate in newly diagnosed glioblastoma is associated with long-term survival. J Clin Neurosci (2010) 17:1381–5. doi:10.1016/j.jocn.2010.03.031

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Kmiecik J, Poli A, Brons NH, Waha A, Eide GE, Enger PØ, et al. Elevated CD3+ and CD8+ tumor-infiltrating immune cells correlate with prolonged survival in glioblastoma patients despite integrated immunosuppressive mechanisms in the tumor microenvironment and at the systemic level. J Neuroimmunol (2013) 264:71–83. doi:10.1016/j.jneuroim.2013.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Rutledge WC, Kong J, Gao J, Gutman DA, Cooper LA, Appin C, et al. Tumor-infiltrating lymphocytes in glioblastoma are associated with specific genomic alterations and related to transcriptional class. Clin Cancer Res (2013) 19:4951–60. doi:10.1158/1078-0432.CCR-13-0551

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Sayour EJ, McLendon P, McLendon R, De Leon G, Reynolds R, Kresak J, et al. Increased proportion of FoxP3+ regulatory T cells in tumor infiltrating lymphocytes is associated with tumor recurrence and reduced survival in patients with glioblastoma. Cancer Immunol Immunother (2015) 64:419–27. doi:10.1007/s00262-014-1651-7

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Patel AP, Tirosh I, Trombetta JJ, Shalek AK, Gillespie SM, Wakimoto H, et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science (2014) 344:1396–401. doi:10.1126/science.1254257

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Meyer M, Reimand J, Lan X, Head R, Zhu X, Kushida M, et al. Single cell-derived clonal analysis of human glioblastoma links functional and genomic heterogeneity. Proc Natl Acad Sci U S A (2015) 112:851–6. doi:10.1073/pnas.1320611111

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Han S, Zhang C, Li Q, Dong J, Liu Y, Huang Y, et al. Tumour-infiltrating CD4(+) and CD8(+) lymphocytes as predictors of clinical outcome in glioma. Br J Cancer (2014) 110:2560–8. doi:10.1038/bjc.2014.162

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Angelova M, Charoentong P, Hackl H, Fischer ML, Snajder R, Krogsdam AM, et al. Characterization of the immunophenotypes and antigenomes of colorectal cancers reveals distinct tumor escape mechanisms and novel targets for immunotherapy. Genome Biol (2015) 16:64. doi:10.1186/s13059-015-0620-6

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A (2005) 102:15545–50. doi:10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Louis DN, Ohgaki H, Wiestler OD, Cavenee WK, Burger PC, Jouvet A, et al. The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol (2007) 114:97–109. doi:10.1007/s00401-007-0278-6

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Yeung JT, Hamilton RL, Ohnishi K, Ikeura M, Potter DM, Nikiforova MN, et al. LOH in the HLA class I region at 6p21 is associated with shorter survival in newly diagnosed adult glioblastoma. Clin Cancer Res (2013) 19:1816–26. doi:10.1158/1078-0432.CCR-12-2861

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Yang I, Han SJ, Sughrue ME, Tihan T, Parsa AT. Immune cell infiltrate differences in pilocytic astrocytoma and glioblastoma: evidence of distinct immunological microenvironments that reflect tumor biology: laboratory investigation. J Neurosurg (2011) 115:505–11. doi:10.3171/2011.4.JNS101172

CrossRef Full Text | Google Scholar

31. Hu Y-L, DeLay M, Jahangiri A, Molinaro AM, Rose SD, Carbonell WS, et al. Hypoxia-induced autophagy promotes tumor cell survival and adaptation to antiangiogenic treatment in glioblastoma. Cancer Res (2012) 72:1773–83. doi:10.1158/0008-5472.CAN-11-3831

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Jahangiri A, De Lay M, Miller LM, Carbonell WS, Hu Y-L, Lu K, et al. Gene expression profile identifies tyrosine kinase c-Met as a targetable mediator of antiangiogenic therapy resistance. Clin Cancer Res (2013) 19:1773–83. doi:10.1158/1078-0432.CCR-12-1281

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Phillips HS, Kharbanda S, Chen R, Forrest WF, Soriano RH, Wu TD, et al. Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell (2006) 9:157–73. doi:10.1016/j.ccr.2006.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Kim Y-H, Jung T-Y, Jung S, Jang W-Y, Moon K-S, Kim I-Y, et al. Tumour-infiltrating T-cell subpopulations in glioblastomas. Br J Neurosurg (2012) 26:21–7. doi:10.3109/02688697.2011.584986

CrossRef Full Text | Google Scholar

35. Lohr J, Ratliff T, Huppertz A, Ge Y, Dictus C, Ahmadi R, et al. Effector T-cell infiltration positively impacts survival of glioblastoma patients and is impaired by tumor-derived TGF-β. Clin Cancer Res (2011) 17:4296–308. doi:10.1158/1078-0432.CCR-10-2557

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Safdari H, Hochberg FH, Richardson EP. Prognostic value of round cell (lymphocyte) infiltration in malignant gliomas. Surg Neurol (1985) 23:221–6. doi:10.1016/0090-3019(85)90086-2

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Thomas AA, Fisher JL, Rahme GJ, Hampton TH, Baron U, Olek S, et al. Regulatory T cells are not a strong predictor of survival for patients with glioblastoma. Neuro Oncol (2015) 17:801–9. doi:10.1093/neuonc/nou363

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Jin L, Ge H, Long Y, Yang C, Chang YE, Mu L, et al. CD70, a novel target of CAR-T-cell therapy for gliomas. Neuro Oncol (2017). doi:10.1093/neuonc/nox116

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Melter M, Reinders MEJ, Sho M, Pal S, Geehan C, Denton MD, et al. Ligation of CD40 induces the expression of vascular endothelial growth factor by endothelial cells and monocytes and promotes angiogenesis in vivo. Blood (2000) 96:3801–8.

PubMed Abstract | Google Scholar

40. Saadoun S, Papadopoulos MC, Hara-Chikuma M, Verkman AS. Impairment of angiogenesis and cell migration by targeted aquaporin-1 gene disruption. Nature (2005) 434:786–92. doi:10.1038/nature03460

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Stabile E, Burnett MS, Watkins C, Kinnaird T, Bachis A, la Sala A, et al. Impaired arteriogenic response to acute hindlimb ischemia in CD4-knockout mice. Circulation (2003) 108:205–10. doi:10.1161/01.CIR.0000079225.50817.71

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Graf MR, Sauer JT, Merchant RE. Tumor infiltration by myeloid suppressor cells in response to T cell activation in rat gliomas. J Neurooncol (2005) 73:29. doi:10.1007/s11060-007-9442-z

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Fridman WH, Pagès F, Sautès-Fridman C, Galon J. The immune contexture in human tumours: impact on clinical outcome. Nat Rev Cancer (2012) 12:298–306. doi:10.1038/nrc3245

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Hewedi IH, Radwan NA, Shash LS, Elserry TH. Perspectives on the immunologic microenvironment of astrocytomas. Cancer Manag Res (2013) 5:293. doi:10.2147/CMAR.S48942

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Cancer Genome Atlas Research Network, Brat DJ, Verhaak RG, Aldape KD, Yung WK, Salama SR, et al. Comprehensive, integrative genomic analysis of diffuse lower-grade gliomas. N Engl J Med (2015) 372:2481–98. doi:10.1056/NEJMoa1402121

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Foote MB, Papadopoulos N, Diaz LA Jr. Genetic classification of gliomas: refining histopathology. Cancer Cell (2015) 28:9–11. doi:10.1016/j.ccell.2015.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Verhaak RGW, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell (2010) 17:98–110. doi:10.1016/j.ccr.2009.12.020

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Yan H, Parsons DW, Jin G, McLendon R, Rasheed BA, Yuan W, et al. IDH1 and IDH2 mutations in gliomas. N Engl J Med (2009) 360:765–73. doi:10.1056/NEJMoa0808710

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Eckel-Passow JE, Lachance DH, Molinaro AM, Walsh KM, Decker PA, Sicotte H, et al. Glioma groups based on 1p/19q, IDH, and TERT promoter mutations in tumors. N Engl J Med (2015) 372:2499–508. doi:10.1056/NEJMoa1407279

CrossRef Full Text | Google Scholar

50. Ge H, Mu L, Jin L, Yang C, Chang YE, Long Y, et al. Tumor associated CD70 expression is involved in promoting tumor migration and macrophage infiltration in GBM. Int J Cancer (2017) 141:1434–44. doi:10.1002/ijc.30830

CrossRef Full Text | Google Scholar

Keywords: gliomas, angiogenesis, tumor-infiltrating lymphocytes, progression, recurrence

Citation: Mu L, Yang C, Gao Q, Long Y, Ge H, DeLeon G, Jin L, Chang Y, Sayour EJ, Ji J, Jiang J, Kubilis PS, Qi J, Gu Y, Wang J, Song Y, Mitchell DA, Lin Z and Huang J (2017) CD4+ and Perivascular Foxp3+ T Cells in Glioma Correlate with Angiogenesis and Tumor Progression. Front. Immunol. 8:1451. doi: 10.3389/fimmu.2017.01451

Received: 13 August 2017; Accepted: 17 October 2017;
Published: 07 November 2017

Edited by:

Salem Chouaib, Institut Gustave Roussy, France

Reviewed by:

Cara Haymaker, University of Texas MD Anderson Cancer Center, United States
Amod Ashok Sarnaik, Moffitt Cancer Center, United States

Copyright: © 2017 Mu, Yang, Gao, Long, Ge, DeLeon, Jin, Chang, Sayour, Ji, Jiang, Kubilis, Qi, Gu, Wang, Song, Mitchell, Lin and Huang. 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) or licensor 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: Zhiguo Lin, linzhiguo@hotmail.com;
Jianping Huang, jianping.huang@neurosurgery.ufl.edu

Present address: Gabriel De Leon, Department of Surgery, Duke University Medical Center, Division of Applied Therapeutics, Durham, North Carolina, United States;
Yifan (Emily) Chang, Laboratory of Molecular and Cellular Neuroscience, The Rockefeller University, New York, NY, United States