ABSTRACT
LGG tumors are characterized by a low infiltration of immune cells, requiring therapeutic interventions to boost the immune response. We conducted a study analyzing mRNA expression datasets from the UCSC Xena web platform. To screen for upregulated genes, we sought to compare normal brain tissue with LGG tumor samples. We also used cBioportal to determine the relationship between mRNA expression levels of 513 LGG patients and their overall survival (OS) outcomes. Three tumor-associated macrophage (TAM) markers, MSR1/CD204, CD86, and CD68, exhibited a 6-fold (p < 0.0001), 8.9-fold (p < 0.0001), and 15.6-fold increase in mRNA expression levels, respectively, in LGG tumors. In addition, both TGFB1 (4.1-fold increase, p < 0.0001) and TGFB2 (2.2-fold increase, p < 0.0001) ligands were also upregulated in these tumors compared to normal brain tissue, suggesting that TGFB ligands are pivotal in establishing an immunosuppressive, angiogenic, and pro-tumorigenic TME in gliomas mediated through TAMs. In addition, mRNA upregulation of interferon-gamma receptors, IFNGR1 and IFNGR2, and the downstream signaling molecules STAT1, IRF1, and IRF5, pointed to an essential role for IFN-γ mediated remodeling of the TME. Interestingly, the mRNA expression of a tumor-associated antigen, CD276/B7-H3, showed a significant (p < 0.0001) 4.03-fold increase in tumor tissue, giving further insights into the roles of macrophages and tumor cells in supporting the immunosuppressive TME. Multivariate Cox proportional hazards models investigating the interaction of TGFB2 and activation of IFNGR2, STAT1, IRF1, or IRF5 showed that the prognostic impact of high mRNA levels (25th percentile cut-off) of TGFB2 was independent of IFNGR2, STAT1, IRF1, or IRF5 mRNA levels (TGFB2high HR (95% CI) = 4.07 (2.35-7.06), 6 (3.62-10.11), 4.38 (2.67-7.17), and 4.48 (2.82-7.12) for models with IFNGR2, STAT1, IRF1, or IRF5, respectively) and age at diagnosis. Patients with high levels of TGFB2 and IFNGR2 were over-represented by LGG patients with isocitrate dehydrogenase wild-type (IDHwt) mutation status. The prognostic impact of high levels of TGFB2 and IDH wild-type observed by the increases in hazard ratios for TGFB2 (HR (95% CI range) = 2.02 (1.05-3.89)) and IDH wild-type (HR (95% CI range) = 4.44 (1.9-10.4)) were independent predictors of survival, suggesting that risk stratification of patients identifies LGG patients with IDH wild-type and high levels of TGFB2 in the design of clinical trials. Furthermore, we have additional IRF5 and CD276/B7-H3 as prognostic markers that can also be targeted for combination therapies with TGFB2 inhibitors. In support of these findings, we demonstrated that low levels of gene methylation in TGFB2, IFNGR2, IRF1, IRF5, STAT1, and CD276 were associated with significantly worse overall survival (OS) outcomes. This suggests that potential mechanisms to increase the expression of these prognostic markers occur via the action of demethylation enzymes.