Abstract
Immune cell infiltration varies widely between different glioblastomas (GBMs). The underlying mechanism, however, remains unknown. Here we show that TGF-beta regulates proliferation, migration, and tumorigenicity of mesenchymal GBM cancer stem cells (CSCs) in vivo and in vitro. In contrast, proneural GBM CSCs resisted TGF-beta due to TGFR2 deficiency. In vivo, a substantially increased infiltration of immune cells was observed in mesenchymal GBMs, while immune infiltrates were rare in proneural GBMs. On a functional level, proneural CSC lines caused a significantly stronger TGF-beta-dependent suppression of NKG2D expression on CD8+ T and NK cells in vitro providing a mechanistic explanation for the reduced immune infiltration of proneural GBMs. Thus, the molecular subtype of CSCs TGF-beta-dependently contributes to the degree of immune infiltration.
Introduction
G
The translational relevance of the molecular subtypes remains vague. A study investigating the response to immunotherapy with adjuvated autologous tumor lysate-pulsed dendritic cells found that patients with mesenchymal GBMs had more infiltrating immune cells and showed longer survival than patients with proneural GBMs. In addition, only mesenchymal GBMs responded to immunotherapy with increased immune infiltration. This finding, however, still needs to be confirmed by a larger number of samples [9]. Thus, the lack of infiltrating immune cells and productive immune responses in patients with proneural GBMs might be due to increased immune-paralyzing effects of the respective tumors [10]. To understand the different degrees of immune paralysis in both subtypes, we used the cell culture system of proneural-like and mesenchymal-like CSC lines and identified the functional role of TGF-beta signaling as a crucial difference.
Materials and Methods
Culture of primary GBM cells
Generation and culture of the CSC lines used has been described previously [7, 6]. Metabolic activity of CSC lines was determined using AlamarBlue (Biosource) according to the manufacturer's instructions. For clonogenicity assays, cells were seeded in 48-well plates and treated as indicated. Medium supplemented with TGF-beta was replaced twice a week.
In vivo tumor model
Cells were intracranially injected into T-lymphocyte-deficient NMRI(nu/nu)-mice as described previously [7]. Procedures were conducted in accordance with German laws governing animal care. Animals were sacrificed 10 weeks after tumor inoculation and 10-μm sections of snap-frozen tumor samples were analyzed. For the CSC line R28, the tumor size was determined by bioluminescence 3 weeks after implantation as described previously [11].
Quantitative real-time PCR
RNA was isolated and cDNA was synthesized using standard protocols. Real-time PCR was performed using the following primer pairs: TGFR2: forward 5′-AGCAACTGCAGCATCACCTC-3′; reverse 5′-GCACTTTGGAGAAGCAGCATC-3′; 18s-RNA: forward 5′-CGGCTACCACATCCAAGGAA-3′; reverse 5′-GCTGGAATTACCGCGGCT-3′. TGF-beta 1: forward 5′- GGCGATACCTCAGCAACCGG-3′; reverse 5′-TCGGCGGCCGGTAGTGAACC-3′; TGF-beta 2: forward 5′-CACCATAAAGACAGGAACCTG-3′; reverse 5′-GGAGGTGCCATCAATACCTGC-3′.
Immunohistochemistry/immunocytochemistry/western blot
The tumor samples were graded by a senior neuropathologist (Joachim Weis) according to WHO classification criteria. The use of the samples had been approved by the local ethics committee. Four-micrometer sections of paraffin-embedded GBM samples were stained using the following antibodies: TGF-beta 1 (Acris:DM1047; 1:100), TGF-beta 2 (Acris:AP15815PU; 1:25), Olig 2 (Sigma-Aldrich; 1:1,250), CD8 (C8/144B, Dako; 1:100), CD68 (KP1, Dako; 1:1,000), DLL3 (St. Crzu; 1:200), NeuN (Millipore; 1:1,000), YKL-40 (Teco Medical group 1:500), CD44 (R&D; 1:200), and VEGF (St. Cruz 1:200). Antibodies used for immunocytochemistry are given in Ref. [7]. All GBM samples were scored blinded by either Markus Riemenschneider or Christoph P. Beier. The index to assign GBM as mesenchymal or proneural tumor was calculated as follows: For each marker (proneural: Olig2, DLL3, NeuN; mesenchymal: YKL-40, CD44, VEGF), 5 positions per slide were scored from 0 to 4 quantifying the number of tumor cells expressing the respective antigen (0: no expression; 1: <10% expression; 2: 10%–49% expression, 3: 50%–90%, 4: >90%). The sum of the means of the mesenchymal markers was then subtracted from the sum of the means of the proneural markers. Tumors with an index above 0 were assigned as proneural; if the index was below 0, the tumor was considered as mesenchymal. For scoring TGF-beta 1/2 and CD68, an identical semiquantitative score was employed. A different score was used for CD8+ cells: 0: no CD8+ cells, 1: 1–10 cells per visual field, 2: 20–29 cells per visual field, 3: 30–50 cells per visual field, 4: more than 50 CD8+ cells per visual field. For pSmad2 western blot, cell lysates were analyzed using pSmad2 (1:1,000; Calbiochem) and beta-actin (1:5,000; Abcam) antibodies.
Immunological assays and flow cytometry
Cells were stained using CD133/2-PE (clone293C3, Miltenyi Biotech; control: mIgG2b-PE; Invitrogen). To assess NKG2D expression, freshly isolated peripheral blood mononuclear cells were cocultured at a 10:1 ratio with CSCs or supernatant in a CSC medium [12]. About 48 h later, NKG2D levels were analyzed (BAT221; anti-NKG2D-PE; Miltenyi Biotech); control: IgG-PE isotype (BioLegend)). CD8+ T cells and NK cells were identified by costaining with anti-CD3-FITC, anti-CD8-PE, and anti-CD56-APC (Immunotools). SD-208 was obtained from Tocris.
Statistics and bioinformatics
Acquisition of microarray data is described in Ref. [7]. Analysis was performed using R and Bioconductor. Expression values were corrected and normalized using the robust multiarray average function (Affymetrix package). Differentially expressed genes were computed using linear models (Bioconductor package limma). We selected 100 genes (TGF target genes) showing the highest log FC (absolute logFC>1.47) and clustered gene samples [13] using complete linkage. Euclidean distances were computed between all samples only based on the TGF-beta target genes. To evaluate the stability of the clustering, consensus clustering was performed. A similar approach was used for the 24-gene signature [8]. The expression indices were calculated per sample by fitting a standard additive model with an independent gene, and sample effects using Tukey's median polish procedure. Parts of the analysis were performed using BRB ArrayTools developed by Dr. Richard Simon and BRB-ArrayTools Development Team.
Results
Differential responsiveness of proneural-like and mesenchymal-like GBM CSC lines
To characterize the role of TGF-beta signaling in both types of GBM CSC lines, we used the cell culture model of CSC lines that closely resembles the initial tumor and recapitulate the intertumoral heterogeneity in vitro [14,15]. TGF-beta-associated transcripts were differentially regulated between CD133+ (proneural-like, type I) CSC lines and CD133− (mesenchymal-like, type II) CSC lines [8]. To further understand these differences, we used 6 mesenchymal-like and 7 proneural-like CSC lines [7,8] that were assigned based on a recently published 24-gene signature (Supplementary Fig. S1C; Supplementary Data are available online at

Responsiveness to TGF-beta in proneural- and mesenchymal-like cancer stem cell (CSC) lines.
TGF-beta-induced effects on mesenchymal- and proneural-like CSC lines
With regard to biological effects, TGF-beta 1, 2, or the combination of both behaved similarly in all experiments (not shown). TGF-beta incubation did affect neither migration, proliferation, or tumorigenicity, nor stem cell properties of those proneural-like CSC lines, which lacked SMAD2-phosphorylation and TGFR2 expression (Fig. 2A–D and Supplementary Fig. S3A–C). Also, the proneural-like CSC line R11, which expresses TGFR2 and showed SMAD2-phosphorylation, was found to resist TGF-beta-induced effects. In this cell line, however, TGF-beta might still modulate the self-renewal of GBM via induction of LIF, as recently described by Penuelas et al. [16]. Nevertheless, as the growth medium for GBM CSCs is routinely supplemented with LIF, LIF-dependent growth-promoting effects triggered by TGF-beta have been obscured for R11 (Supplementary Fig. S2B). In the mesenchymal-like CSC lines, in contrast, TGF-beta either suppressed or promoted migration, proliferation, and tumorigenicity (Fig. 2A–D). Moreover, as opposed to tumorigenicity, other specific stem cell properties like clonogenicity, proportion of CD133+ cells, or the differentiation profile remained unaltered after in vitro treatment with TGF-beta (Supplementary Fig. S3A–C). An analysis of factors like DAB2 that may decide whether TGF-beta acts as tumor suppressor or promoter [17] would, however, require a much larger number of cell lines and therefore exceeds the scope of the present study.

Comparison of TGF-beta effects in proneural-like and mesenchymal-like CSC lines.
Mesenchymal, but not proneural GBMs express TGF-beta-induced transcripts
To identify TGF-beta responsive GBMs, we compared gene expression profiles of TGF-beta responsive CSC lines (R8/R53) with and without TGF-beta treatment and determined a set of 100 TGF-beta target genes (Supplementary Table S1). We then analyzed 80 previously published microarrays from well-characterized GBM samples [13]. Clustering based on the TGF-beta target genes unveiled 2 groups: a TGF-betaresponsive cluster characterized by expression patterns typical for activated TGF-beta and a TGF-betaunresponsive cluster, suggesting that autocrine TGF-beta signaling is only active in a subgroup (Fig. 3A). We then clustered the same set of 80 GBMs based on the 24-gene signature differentiating between proneural-like and mesenchymal-like CSCs in vivo and in vitro [8] and again obtained 2 stable clusters. To compare the expression of TGF-beta target genes and 24 signature genes, we calculated 2 expression indices. The TGF-beta target genes and the 24-gene signature index anticorrelate (r=−0.885) with most TGF-betaresponsive GBMs resembling mesenchymal-like CSC lines. Conversely, TGF-betaunresponsive GBMs mainly displayed a proneural-like CSC background (Fig. 3B).

Only GBMs resembling mesenchymal-like CSCs express TGF-beta-induced transcripts.
TGF-beta-dependent immunosuppression by proneural-like CSCs
The combination of TGF-beta resistance and TGF-beta expression in proneural-like CSCs suggests an immunosuppressive function of TGF-beta in this subtype. As the respective GBM patients are deceased, we could not compare cytotoxic and proliferative T-cell assays in a syngeneic context. Instead, we investigated the TGF-beta-induced downregulation of NKG2D receptor expression on CD8+ T-cells and NK-cells in vitro [18]. The TGFR1 kinase inhibitor SD-208 specifically blocks TGF-beta signaling and thereby unmasks TGF-beta-dependent effects [19]. In the absence of immune cells, CSC lines did not secrete active TGF-beta, because coculture of PBMC with supernatants of CSC lines remained without effect (Fig. 4C). Cocultures of PBMC with CSC lines revealed an increased TGF-beta-dependent downregulation of NKG2D when cocultured with proneural-like CSC lines (Fig. 4A, B, D). Thus, TGF-beta-resistant proneural-like CSC lines exert enhanced, TGF-beta-dependent immunosuppression in 2 highly sensitive functional tests in vitro. We then investigated xenotransplants (Fig. 2D) from mesenchymal-like and proneural-like CSC lines for infiltration with CD68+ microglia. Tumors formed by both mesenchymal-like CSC lines, but not the tumors formed by the proneural-like CSC line were infiltrated by microglia further supporting the concept of differential immunosuppressive properties of proneural- and mesenchymal-like CSCs (Fig. 4E). One possible explanation for the different bioactivity of TGF-beta between the 2 subtypes could be the differential expression of Latent TGF-binding protein 1 (LTBP1), which retains TGF-beta in the extracellular matrix in an inactivate state (Fig. 4F).

Proneural-like CSC lines downregulate NKG2D receptor expression on CD56+ NK and CD8+ T cells in a TGF-beta-dependent manner. 1.5×106 PBMC were cocultured with 1.5×105 proneural-like (R18, R54) or mesenchymal-like CSC lines (R8, R53)
Immune infiltration of GBM correlates to the molecular subtype
To further substantiate this hypothesis of an increased immunosuppression by proneural CSCs, we stained a series of GBMs against markers for cytotoxic T-lymphocytes (CD8+) and CD68+ microglia (Fig. 5A). The degree of immune infiltration varied. To differentiate between proneural and mesenchymal GBMs, we used a set of marker genes (proneural: Olig2, DLL3, NeuN; mesenchymal: YKL-40, CD44, VEGF, Supplementary Fig. S1A). We calculated an index based on the average expression of the marker genes (Supplementary Fig. S1B), indicating the assignment to proneural and mesenchymal GBMs. The assignment to the mesenchymal phenotype significantly correlated with the amount of CD8+ (CD8: r=0.68) and CD68+ cells (r=0.78, Fig. 5B–D), suggesting that proneural and mesenchymal GBMs differ with respect to infiltration by immune cells.

Proneural and mesenchymal GBMs show different levels of infiltration by immune cells.
Discussion
The relationship of GBM to proneural-like or mesenchymal-like CSCs indicates a profound biological difference in the function of TGF-beta in these tumors. While tumors expressing transcripts of proneural-like CSCs show TGF-beta expression in vitro and in vivo, they are not susceptible toward pro- or antiproliferative effects of TGF-beta. As TGF-beta is the most potent immunosuppressive cytokine known, it may suppress the function and invasion of immune cells in proneural GBMs [9,20]. We could substantiate this hypothesis using functional assays showing a significantly stronger TGF-beta-dependent suppression of NKG2D receptor expression on NK- and T-cells by proneural-like CSC lines as compared to mesenchymal controls (Fig. 4A–D). In line with these results, proneural GBMs suppressed invasion of immune cells more efficiently as compared to mesenchymal GBMs (Fig. 5B–D). These findings bear substantial implication, because they suggest that therapies aiming at the attenuation of TGF-cells beta-mediated immunosuppression are likely to be essential for every immunotherapy applied to tumors maintained by proneural-like CSCs. Conversely, the strong immunosuppression by TGF-beta provides a possible explanation for the completely refractory phenotype of proneural GBMs in an otherwise highly promising vaccination study on GBM patients [9].
Conversely, the functions of TGF-beta in GBMs maintained by mesenchymal-like CSCs are manifold and TGF-beta may either inhibit or promote tumor growth [17]. While the suppression of invading immune cells may be reduced due to the increased expression of Latent TGF-beta binding protein 1, it may also function as tumor suppressor in at least a subgroup, which calls for an extremely cautious application of TGF-beta targeting therapies to mesenchymal GBMs.
In summary, the CSC subtype determines the degree of immune infiltration and local immune paralysis via a different degree of TGF-beta-mediated suppression of infiltrating immune cells. Thus, the identification of the CSC subtype in GBMs based on the 24-gene signature will help to guide both TGF-beta targeting and immune therapies by identifying GBMs with strong TGF-beta-dependent immune paralysis. The CSC-guided subclassification of GBMs therefore provides a new approach to guide both immune and TGF-beta targeting therapies.
Footnotes
Acknowledgments
This work was supported by a grant of the Bavarian Research Foundation to D.B., U.B., and P.K. and a career award to young female scientists of the RWTH Aachen, Medical School to D.B. Parts of the project were funded by BayGene and the NGFNplus Brain Tumor Network to C.P.B. and P.H. (Subproject 7A and 7B No. 01GS0887 and 01GS1105) and the in-house research support of the RWTH Aachen, Medical School (START-Program) to D.B.
Author Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
