ORIGINAL RESEARCH article

Front. Genet., 16 June 2021

Sec. Computational Genomics

Volume 12 - 2021 | https://doi.org/10.3389/fgene.2021.679097

Prognostic Value of MEG3 and Its Correlation With Immune Infiltrates in Gliomas

  • 1. Department of Neurosurgery, The Second People’s Hospital of Yibin, Yibin, China

  • 2. Department of Health Examination, The Second People’s Hospital of Yibin, Yibin, China

Abstract

Accumulating evidence has revealed that dysregulated lncRNA expression contributes to the onset and progression of cancer. However, the mechanistic role of lncRNA in glioma progression and tumor immunology remains largely unknown. This study aimed to evaluate the significance of maternally expressed gene 3 (MEG3) in the prognosis of and its immune-related roles in gliomas. The expression levels of MEG3 were analyzed using Oncomine and TIMER database. As an important imprinted gene, the copy number variation (CNV) of MEG3 in both glioblastoma multiforme (GBM) and low-grade glioma (LGG) were analyzed using GSCALite database, whereas its prognostic significance was assessed using PrognoScan and GEPIA databases. The relationship between MEG3 and tumor-infiltrated immune cells was analyzed using TIMER. Results showed that MEG3 expression was lower in most of the human cancer tissues than in the normal tissues. We also found that heterozygous deletion of MEG3 occurred more frequent than heterozygous amplification in gliomas, and mRNA expression of MEG3 was significantly positively correlated with its CNV in both the GBM and LGG group. Survival analysis showed that the CNV level of MEG3 had significant correlation with overall survival (OS) and progression-free survival (PFS) compared with wild type in LGG. Lower MEG3 expression was related with poor prognosis. Further analysis showed that in GBM, MEG3 expression level was significantly positively correlated with that of infiltrating CD8+ T cells and significantly negatively correlated with that of infiltrating dendritic cells. In LGG, MEG3 expression level was significantly negatively correlated with levels of infiltrating B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells. Univariate Cox survival analysis demonstrated that only the level of infiltrating dendritic cells significantly affected the survival time of patients with GBM, while all six types of immune cells had a significant effect on the survival time of patients with LGG. Furthermore, MEG3 expression showed strong correlations with multiple immune markers in gliomas, especially in LGG. The current findings suggest that MEG3 expression might serve as a possible prognostic marker and potential immunotherapeutic target for gliomas.

Introduction

Glioma is the most common primary malignant tumor of the brain. Despite many breakthrough analyses deciphering the underlying molecular mechanisms of gliomas, comprehensive treatment options are still lacking and the long-term survival rate of glioma patients remains poor. The molecular complexity and unique microenvironment in the brain that lead to therapeutic resistance, disease progression, and tumor recurrence add to the challenges of glioma treatment. Efforts are under way to uncover key molecular mechanisms, which may lead to the development of new therapeutics for glioma patients.

Long non-coding RNA (lncRNA) is an RNA with a length of more than 200 nt. It is known to play an important role in carcinogenesis and cancer progression. Accumulating evidence has revealed that dysregulated lncRNA expression contributes to the development and metastasis of tumors (). However, the mechanistic role of lncRNA in glioma progression and tumor immunology remains largely unknown. Maternally expressed gene 3 (MEG3), also known as gene-trap locus 2 (GTL2), is an imprinted gene. It is highly expressed in the brain, pituitary gland, placenta, and adrenal gland. MEG3 has been reported to participate in the regulation of a variety of diseases. The first evidences of the contribution of MEG3 in human cancers were obtained from pituitary non-functioning adenomas (Zhou et al., 2012). Several studies have shown downregulation of MEG3 in human cancers, such as breast cancer (Zhang et al., 2016), liver cancer (Zheng et al., 2018), gliomas (Zhang et al., 2017), lung cancer (Wu et al., 2018), squamous cell carcinoma (), and gastrointestinal cancer (). It was found that MEG3 expression was significantly reduced or completely lost in 25% of neuroblastomas (), 81% of hepatocellular cancers (), and 82% of gliomas (). A recent meta-analysis has also shown a correlation between MEG3 downregulation and poor patient outcomes ().

The tumor microenvironment is now widely accepted as an important regulator of cancer progression and therapeutic response (). Previous studies have highlighted the association between immune cells and lncRNAs, which are responsible for differentiation, development, and activation of immune cells (). The role of lncRNAs as regulators in the immune modulation of cancer is emerging (; ). Recent studies indicate that MEG3 is involved in the regulation of CD4+ T cell activation in aplastic anemia () and the expression level of PD-L1 in aggressive endometrial cancer (Xu et al., 2020). However, to the best of our knowledge, the immune-specific functions of MEG3 expression in gliomas have not yet been evaluated. The success of immunotherapies (; ) calls for the evaluation of the role of MEG3 in gliomas. Hence, we conducted the present bioinformatics study to detect the expression levels of MEG3 in gliomas using Oncomine and “Tumor Immune Estimation Resource (TIMER)” databases, analyze CNV using GSCALite database, and analyze its prognostic significance using PrognoScan and “Gene Expression Profiling Interactive Analysis (GEPIA)” databases. We further investigated the correlation between MEG3 expression and immune cell infiltration in gliomas by using the TIMER database.

Materials and Methods

Gene Expression Analysis

The gene expression level of MEG3 in various cancers was analyzed using Oncomine ()1 and TIMER ()2 databases. The parameters were set as follows in the Oncomine database: P-value of 0.001, fold change of 2, and gene ranking of top 10%, six studies meeting the criteria were included for further analysis, including TCGA Brain, Sun Brain (), Murat Brain (), Bredel Brain2 (), Liang Brain (), and Shai Brain (). These studies reported a total of 955 cases, of which 932 are primary tumor samples, and 904 samples have MEG3 profile data. We further evaluated the expression level of MEG3 using TIMER (P < 0.05).

Copy Number Variation (CNV) Analysis

Monoallelic expression dependent on the parent of origin is a hallmark of imprinted genes; it is often lost in certain tumors (). We further analyzed the CNV of MEG3 in both GBM and LGG using GSCALite database ()3, which consists of analytic modules for data from three major sources including multi-omics data from TCGA, Genomics of Drug Sensitivity in Cancer (GDSC) (Yang et al., 2013), and Cancer Therapeutics Response Portal (CTRP) (), and normal tissue expression data from GTEx.

Prognostic Value Analysis

We further analyzed the correlation between MEG3 expression and clinical prognosis in various types of cancers using the PrognoScan database4, which is a powerful tool to investigate the prognostic value of genes (). All survival values, such as overall survival (OS), relapse-free survival (RFS), disease-specific survival (DSS), disease-free survival (DFS), and distant metastasis-free survival (DMFS), were included in the present analysis. Only those studies with corrected P < 0.05 were included in further analyses. The GEPIA database ()5 was used to generate survival curves, including OS, and RFS, based on gene expression using the log-rank test and Mantel–Cox test in both LGG and GBM.

Immune Infiltration Analysis

Tumor immune estimation resource (TIMER) is a comprehensive database designed for the analysis of immune cell infiltrates across different types of cancers (see text footnote 2) (). We initially analyzed the correlation of MEG3 expression with the abundance of six immune infiltrates, including B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells, using the TIMER “gene” module. Furthermore, we used the Kaplan–Meier method to analyze the influence of MEG3 expression and immune cell infiltration on the prognosis of GBM and LGG patients; we also constructed a multivariate Cox proportional risk model. Lastly, the correlation of MEG3 expression with different immune infiltrating cell markers was also assessed through the correlation module. These immune infiltrating cells, including T cells, B cells, monocytes, tumor-associated macrophages (TAMs), macrophages, neutrophils, natural killer (NK) cells, dendritic cells, different T-helper cells, Tregs, and exhausted T cells, and the immune cell markers were referenced as described in the previous studies (; Zhang et al., 2019). The GEPIA database was used to confirm the gene correlations identified by TIMER analysis. The tumor tissue datasets were used for the analysis, and all default setting values were accepted.

Statistical Analysis

The results generated using the Oncomine and TIMER databases are shown with the P-value. Survival curves were generated using PrognoScan and Kaplan–Meier plots. The results of Kaplan–Meier plots and PrognoScan are accompanied with the hazard ratio (HR) and P or Cox P-values from a log-rank test. The correlation of gene expression was evaluated by Spearman’s correlation and a P < 0.05 was considered statistically significant.

Results

Analysis of MEG3 Expression Levels

MEG3 expression level in different human cancers and normal tissues were analyzed using Oncomine database. It was observed that MEG3 expression was lower in most tumor tissues than in normal tissues. In brain and CNS cancers, six studies reported the downregulation of MEG3, whereas no study reported its upregulation. A comprehensive analysis of 15 studies yielded a median rank = 793 and P = 6.29E-7, indicating that MEG3 was downregulated in glioma tissues (Figure 1A). Furthermore, we evaluated the MEG3 expression level through the TIMER database, finding that MEG3 expression was significantly lower in BLCA (bladder urothelial carcinoma), BRCA (breast invasive carcinoma), ESCA (esophageal carcinoma), GBM (glioblastoma multiforme), KICH (kidney chromophobe), KIRP (kidney renal papillary cell carcinoma), PRAD (prostate adenocarcinoma), UCEC (uterine corpus endometrial carcinoma), than in normal controls. By contrast, MEG3 expression was significantly higher in LUAD (lung adenocarcinoma) than in normal tissue (Figure 1B). We also further analyzed the expression of MEG3 in different grades of gliomas, finding that the expression level of MEG3 in grade 4 was significantly lower than that in other grades of gliomas and normal control, it decreased with the increasing of malignant degree of gliomas (Figure 1C). However, no significant difference was observed with respect to sex (Figure 1D). Finally, the expression of MEG3 showed a significant difference in LGG and GBM compared with that in normal control through the GEPIA database (Supplementary Figure 1). Taken together, these results were consistent with previous studies that MEG3 might serve as a tumor suppressor in most human cancers, including gliomas ().

FIGURE 1

Analysis of CNV of MEG3 Expression

Imprinted gene plays an important role in different biological processes, and loss of imprinting has been found in various cancers (). By analyzing GSCALite database, we found that the heterozygous deletion of MEG3 occurred in 27.73 and 20.08% of GBM and LGG, respectively (Figures 2A,C), and Pearson correlation analysis showed that mRNA expression of MEG3 was significantly positively correlated with its CNV in both GBM (r = 0.32, P < 0.05) and LGG group (r = 0.16, P < 0.05) (Figure 2B), which indicated that mRNA expression of MEG3 was significantly affected by CNV. Survival analysis showed that CNV level of MEG3 had significant correlation with OS and progression-free survival (PFS) compared with wild type in LGG, but not in GBM (Figure 2D).

FIGURE 2

Analysis of the Clinical Prognostic Impact of MEG3 Expression

PrognoScan was used to investigate the potential prognostic impact of MEG3 expression level in patients suffering from different types of cancer. The results are summarized in Figure 3. Notably, MEG3 expression was significantly correlated with the prognosis of a total of nine types of cancers, including bladder cancer, brain cancer, breast cancer, head and neck cancer, blood cancer, colorectal cancer, eye cancer, lung cancer, and ovarian cancer. We found that lower MEG3 expression was often associated with a poorer prognosis in these cancer patients.

FIGURE 3

We further examined the potential impact of the MEG3 expression on survival rate in LGG and GBM patients using the GEPIA database. Interestingly, the results illustrated that MEG3 expression level had no significant correlation with OS and RFS in LGG and GBM (Supplementary Figure 2). These results confirmed the differential prognostic value of MEG3 in specific types of cancer.

Analysis of Immune Infiltration

“Tumor-infiltrating lymphocytes” are a key predictor of survival as well as response to immunotherapy in the cancer patients (; ). TIMER database was used to evaluate whether the expression of MEG3 in glioma is correlated with immune infiltration. Interestingly, we found that in GBM, the MEG3 expression level was significantly positively correlated with level of infiltrating CD8+ T cells (r = 0.206, P = 2.28e-05) and significantly negatively correlated with the level of infiltrating dendritic cells (r = –0.18, P = 2.17e-04). On the other hand, in LGG, the MEG3 expression level was significantly negatively correlated with levels of infiltrating B cells (r = –0.309, P = 5.22-12), CD8+ T cells (r = –0.299, P = 2.39e-11), CD4+ T cells (r = –0.243, P = 7.80e-08), macrophages (r = –0.297, P = 4.51e-11), neutrophils (r = –0.22, P = 1.35e-06), and dendritic cells (r = –0.326, P = 2.87e-13) (Figure 4 and Table 1). Next, to evaluate the correlation between immune cell infiltration and the prognosis of glioma patients, “Survival” module was used to generate Kaplan–Meier plots. Of note, we found that only dendritic cell infiltration was significantly correlated with GBM prognosis, while infiltration of all six types of immune cells showed significant correlation with LGG prognosis (Figure 5). Lastly, to explore the clinical relevance of immune cell subsets in gliomas, Cox proportional hazard model was constructed. It was observed that only dendritic cell infiltration was significantly associated with OS of patients with GBM as revealed by the univariate Cox survival analysis. Six types of immune cells significantly affected the survival time of patients with LGG; however, MEG3 expression did not significantly affect the survival time in these patients (Table 2). In addition, multivariate Cox survival analysis revealed age, macrophages, MEG3 expression, and neutrophils to be independent prognostic biomarkers for LGG (Table 3), and age and dendritic cells to be the independent prognostic biomarkers for GBM (Table 4). These findings suggest that MEG3 plays an important immune-related role in gliomas, particularly in LGG.

FIGURE 4

TABLE 1

CancerVariablePartial.corP
LGGDendritic Cell–0.326272.87E-13
LGGB Cell–0.308655.22E-12
LGGCD8+ T Cell–0.299292.39E-11
LGGMacrophage–0.296774.51E-11
LGGCD4+ T Cell–0.243147.80E-08
LGGNeutrophil–0.219571.35E-06
GBMCD8+ T Cell0.2055792.28E-05
GBMDendritic Cell–0.179950.000217
GBMPurity0.088220.071239
GBMNeutrophil0.0469960.337813
GBMB Cell–0.04530.355574
LGGPurity0.0388660.396039
GBMCD4+ T Cell0.0202120.680316
GBMMacrophage–0.002120.96552

MEG3 expression associated with immune infiltration level in GBM and LGG.

FIGURE 5

TABLE 2

CancerVariableP
LGGNeutrophil5.83E-06
LGGMacrophage9.20E-06
LGGB Cell4.25E-05
LGGCD4+ T Cell0.00046
LGGDendritic Cell0.000634
GBMDendritic Cell0.001617
LGGCD8+ T Cell0.009535
GBMB Cell0.105243
GBMMEG30.262632
LGGMEG30.541297
GBMCD4+ T Cell0.629122
GBMNeutrophil0.741002
GBMMacrophage0.755369
GBMCD8+ T Cell0.805153

Univariate analysis of the correlation of MEG3 expression and immune infiltrates with OS in glioma.

TABLE 3

CoefHR (95%CI_l–95%CI_u)P-valueSig
Age0.061.062 (1.045–1.079)0***
Macrophage6.607740.156 (12.152–45,083.396)0.002**
MEG30.1691.184 (1.022–1.373)0.024*
Neutrophil-8.3670 (0–0.956)0.049*
CD8_Tcell7.0571160.386 (0.832–1,617,980.188)0.056
Dendritic2.43211.383 (0.177–731.417)0.252
Gendermale0.1681.183 (0.787–1.776)0.419
B_cell1.6975.455 (0.008–3945.616)0.614
Purity0.1681.183 (0.44–3.178)0.739
CD4_Tcell-1.0080.365 (0–1071.192)0.805
raceBlack15.9418,377,586.947 (0–Inf)0.994
raceWhite16.15610,383,164.77 (0–Inf)0.994

Multivariate analysis of the correlation of MEG3 expression and immune infiltrates with OS in LGG.

*P < 0.01; **P < 0.001; ***P < 0.0001.

TABLE 4

CoefHR (95%CI_l–95%CI_u)P-valueSig
Age0.0291.030 (1.021–1.038)0***
Dendritic0.5611.753 (1.275–2.410)0.001**
CD8_Tcell0.3541.424 (0.907–2.235)0.124
B_cell-0.480.619 (0.323–1.186)0.148
RaceWhite0.5511.736 (0.705–4.272)0.23
Macrophage0.4381.550 (0.721–3.331)0.262
MEG30.0991.104 (0.891–1.368)0.366
Gendermale0.0821.085 (0.863–1.365)0.484
RaceBlack0.2891.335 (0.482–3.701)0.578
CD4_Tcell0.171.186 (0.540–2.603)0.671
Neutrophil0.1161.123 (0.446–2.825)0.806
Purity0.051.051 (0.539–2.050)0.884

Multivariate analysis of the correlation of MEG3 expression and immune infiltrates with OS in GBM.

**P < 0.001; ***P < 0.0001.

Assessment of Correlation Between MEG3 and Immune Marker Expression

To further investigate the relationship between MEG3 and level of immune cell infiltration, we assessed the correlation between MEG3 expression and immune marker genes of various immune cells in glioma using TIMER and GEPIA databases. Different subsets of immune cells, including CD8+ T cells, T cells, B cells, monocytes, TAMs, M1 and M2 macrophages, neutrophils, NK cells, and dendritic cells, were analyzed in LGG and GBM, using GBM as the control group. Functional T cells, such as Th1 cells, Th2 cells, Tfh cells, Th17 cells, and Tregs, and exhausted T cells were also analyzed. The results showed that, after adjustment by purity, MEG3 expression level was significantly negatively correlated with most immune marker sets of various immune cells in LGG. However, interestingly, we found that only dendritic cell markers (HLA-DRA, CD11c) were significantly negatively correlated with MEG3 expression level in the GBM control group (Table 5). Similar results were also obtained while analyzing the correlation between MEG3 expression and the above-stated markers by using the GEPIA database (Table 6). Therefore, these results further confirm the findings that MEG3 is specifically correlated with immune infiltrating cells in gliomas, especially in LGG; this suggests that MEG3 plays a vital role in immune escape in the glioma microenvironment.

TABLE 5

DescriptionGene markersLGG
GBM
None
Purity
None
Purity
CorPCorPCorPCorP
CD8+ T cellCD8A0.0850.0540.1060.0210.0210.7930.0190.822
CD8B–0.180***–0.156***–0.0830.308–0.0940.277
T cell (general)CD3D–0.204***–0.178**–0.213*–0.2650.018
CD3E–0.223***–0.213***–0.150.064–0.1820.034
CD2–0.25***–0.238***–0.1860.021–0.229*
B cellCD19–0.15***–0.136**–0.0470.562–0.0750.385
CD79A–0.0670.127–0.0580.205–0.0170.831–0.0620.474
MonocyteCD86–0.328***–0.359***–0.1820.024–0.2120.013
CD115 (CSF1R)–0.25***–0.27***–0.1130.165–0.1520.076
TAMCCL2–0.221***–0.23***–0.110.174–0.1120.193
CD68–0.366***–0.379***–0.1070.187–0.1010.240
IL10–0.257***–0.268***–0.1550.056–0.1530.075
M1 macrophageINOS (NOS2)0.155**0.183***0.0110.8970.0160.851
IRF5–0.2***–0.211***0.060.4630.0650.449
COX2(PTGS2)0.0850.0550.0690.1320.0150.8530.0340.689
M2 macrophageCD163–0.254***–0.258***–0.1360.094–0.150.081
VSIG4–0.275***–0.292***–0.150.065–0.1420.099
MS4A4A–0.341***–0.349***–0.227*–0.2260.008
NeutrophilsCD66b (CEACAM8)0.0240.580.0340.4530.1430.0780.1620.058
CD11b (ITGAM)–0.228***–0.25***0.0180.830–0.0010.990
CCR7–0.0920.036–0.0780.090–0.0840.3–0.1130.190
Natural killer cellKIR2DL1–0.0730.097–0.1020.026–0.0730.373–0.0490.568
KIR2DL3–0.0610.168–0.0480.2920.040.6260.0170.845
KIR2DL4–0.15**–0.153**–0.2060.011–0.254*
KIR3DL1–0.0260.551–0.0170.716–0.060.462–0.0660.443
KIR3DL2–0.126*–0.149*–0.1390.087–0.130.129
KIR3DL30.0340.4460.0210.6540.0890.2720.0560.515
KIR2DS4–0.0860.052–0.0920.044–0.0590.469–0.0630.467
Dendritic cellHLA-DPB1–0.322***–0.326***–0.2070.010–0.234*
HLA-DQB1–0.252***–0.255***–0.1640.043–0.1710.045
HLA-DRA–0.335***–0.348***–0.242*–0.285**
HLA-DPA1–0.338***–0.342***–0.220.006–0.23*
BDCA-1(CD1C)–0.235***–0.23***–0.0680.404–0.0950.268
BDCA-4(NRP1)–0.0820.062–0.0880.055–0.010.9–0.010.909
CD11c (ITGAX)–0.050.26–0.0380.4020.296**0.318**
Th1T-bet (TBX21)–0.0050.901–0.0010.9880.1640.0430.1360.114
STAT40.249***0.273***0.1480.0670.1740.042
STAT1–0.162**–0.174***–0.0370.651–0.0530.540
IFN-g (IFNG)–0.1010.022–0.10.0290.0610.4510.040.640
TNF-a (TNF)–0.1030.020–0.1160.0110.0670.4090.0490.573
Th2GATA3–0.174***–0.149*–0.0860.293–0.1330.120
STAT60.0010.9890.0050.9170.1460.0720.1340.117
STAT5A–0.2***–0.22***–0.0690.397–0.0780.368
IL130.413***0.414***0.212*0.23*
TfhBCL60.0950.0300.0810.0750.204*0.1960.022
IL210.0050.9170.010.822–0.0060.942–0.0180.834
Th17STAT3–0.258***–0.281***0.0320.6910.0140.868
IL17A–0.0140.743–0.0320.4880.0560.4950.0550.524
TregFOXP30.148**0.168**–0.0430.596–0.0490.566
CCR8–0.0340.441–0.0240.602–0.1190.141–0.1330.122
STAT5B0.0440.3160.0250.5890.1620.0450.1510.078
TGFb (TGFB1)–0.177***–0.191***0.0620.4480.0260.762
T cell exhaustionPD-1 (PDCD1)–0.124*–0.1150.0120.0210.7980.0180.836
CTLA40.0380.4560.0580.203–0.0450.585–0.0590.492
LAG30.119*0.125*0.1150.1580.110.200
TIM-3 (HAVCR2)–0.3***–0.325***–0.1070.186–0.1240.148
GZMB–0.1480.275–0.0390.394–0.1090.18–0.1290.133

Correlation analysis between MEG3 and relate genes and markers of immune cells in TIMER.

*P < 0.01; **P < 0.001; ***P < 0.0001.

TABLE 6

DescriptionGene markersLGG
RP
Dendritic cellHLA-DPB1–0.11*
HLA-DQB1–0.0580.19
HLA-DRA–0.094*
HLA-DPA1–0.080.07
BDCA-1(CD1C)–0.0770.081
BDCA-4(NRP1)0.16***
CD11c (ITGAX)–0.0320.47
MonocyteCD86–0.19***
CD115 (CSF1R)–0.17***
TAMCCL20.0550.21
CD68–0.19***
IL10–0.0580.18
M1 MacrophageIRF5–0.13***
COX2(PTGS2)0.0360.41
M2 MacrophageVSIG4–0.13***
MS4A4A–0.11***

Correlation analysis between MEG3 and relate genes and markers of dendritic cell, monocyte, TAM, and macrophages in GEPIA.

*P < 0.01; ***P < 0.0001.

Discussion

In the present study, we analyzed the expression level of the lncRNA MEG3 in glioma and other tumors using the Oncomine database. The results showed that MEG3 was expressed at a lower level in most cancer types, including gliomas, than in normal tissues. This observation is consistent with the finding in the previous study reporting that MEG3 might have a tumor suppressive role in gliomas (). We also found that heterozygous deletion of MEG3 occurred more frequent than heterozygous amplification in gliomas; mRNA expression of MEG3 was significantly positively correlated with its CNV in both GBM and LGG groups. Survival analysis showed that the CNV level of MEG3 had significant correlation with OS and PFS in LGG compared with wild type. Analysis pertaining to the prognostic impact showed that increased expression of MEG3 was correlated with improved survival in multiple types of human cancers. Through immune infiltration analysis, we found that MEG3 expression was positively correlated with the degree of CD8+ T cell and dendritic cell infiltration in GBM as well as with the degree of dendritic cell, B cell, CD8+ T cell, macrophage, CD4+ T cell, and neutrophil infiltration in LGG. We further found that infiltration of all these immune cells was significantly associated with LGG prognosis, while only dendritic cell infiltration was significantly associated with GBM prognosis. Furthermore, our analyses showed that levels of immune cell infiltration and diverse immune marker sets were correlated with the level of MEG3 expression in gliomas, especially in LGG.

The key finding of this study is that MEG3 expression is correlated with diverse immune cell infiltration levels in gliomas, particularly in LGG. Due to the restrictions imposed by the blood–brain barrier (BBB), the brain has long been considered as an immune privileged organ; nonetheless, the brain is now proposed to be an immunologically distinct organ (). Our understanding of the brain tumor microenvironment remains limited. Immune cells constitute an important component of the glioma microenvironments and they can reach up to 50% of the total cell mass in some tumors (). Gliomas are enriched in diverse immune cell types including neutrophils and myeloid-derived suppressor cells (). On the other hand, brain metastases are enriched in lymphocytes and neutrophils. Moreover, brain metastases of different origin showed distinct immune landscape (). The different immunotherapeutic effects in certain tumors may be attributed to differences in the disease- and cell-type-specific tumor microenvironments. A previous study showed that downregulation of MEG3 promotes angiogenesis after ischemic brain injury (). BBB disruption following stroke promotes inflammation by enabling leukocytes, T cells, and other immune cells to migrate across the BBB (), which may indirectly explain the potential mechanism of the role of MEG3 in the immune microenvironment. The tumor immune microenvironment influences cancer immune escape, response to immunotherapy, and the survival rate of patients (Zhang et al., 2020). The observed correlation between MEG3 and the expression of certain immunological marker genes suggests that, in gliomas, MEG3 may interact with infiltrated immune cells within the tumor microenvironment. Of note, evidence that any therapeutic intervention administered to glioma patients has a major effect on survival is lacking. Our analysis would help in an efficient implementation of immunotherapy as a part of the standard care for patients with glioma.

Functional studies indicate that lncRNAs may act as oncogenes or tumor suppressors in various cancers (). The results obtained from PrognoScan database analysis were consistent with the proposed tumor suppressor role for MEG3. Several studies have shown a correlation between MEG3 expression and patients’ survival. Low MEG3 expression level was correlated with poor prognosis (Zhao et al., 2018). A recent meta-analysis has verified that low expression of MEG3 is significantly associated with poor prognosis for patients with different types of cancer (), although GEPIA database analysis showed that changes in MEG3 expression had no significant effect on the OS and DFS in both LGG and GBM. Our study showed that immune cell infiltration has close correlation with prognosis in glioma patients. The association between MEG3 expression and immune cell infiltration indirectly indicated that MEG3 may be a useful predictor of glioma patient survival. Furthermore, using the PrognoScan database, we observed that MEG3 expression was significantly correlated with the prognosis of a total of nine cancer types, including gliomas. These findings warrant validation in a larger cohort as previous sample size was limited. MEG3 inhibits cancer progression through several independent mechanisms, including regulation of the tumor suppressor genes p53 and Rb (; ), inhibition of angiogenesis (), acting as a competitive endogenous (ce)RNA (; Zhang et al., 2017; Zhang and Guo, 2019; ), and induction of EMT and invasion via autophagy (; Yang et al., 2020). Although the role and prognostic value of the dysregulation of MEG3 in gliomas has not yet been fully evaluated, the current findings suggest that MEG3 may regulate tumor progression via immune-based mechanisms. Further studies are needed to establish its exact molecular mechanism in the progression of gliomas.

The tumor database-based studies have advantages of a larger sample size and higher reliability over traditional studies; such studies also provide a preliminary foundation for further studies. Although our study provides insights into understanding the potential role of MEG3 in tumor immunology and its use as a cancer biomarker, it is limited in terms of lack of in vitro and in vivo experiments to validate the relationship between MEG3 expression and infiltration of immune cells in gliomas. Therefore, further research is needed to verify the role of MEG3 in gliomas.

Conclusion

Taken together, the current findings suggest that the expression of MEG3 might serve as a possible prognostic biomarker and potential immunotherapeutic target for gliomas.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

XX analyzed the data and wrote the manuscript. YY reviewed the data and manuscript. ZZ and YS analyzed the data, edited the manuscript, and supervised the study. All authors contributed to the article and approved the final submitted version of the manuscript.

Funding

This study was supported in part by the Sichuan Science and Technology Program (Grant No. 2021YJ0164).

Acknowledgments

We would like to thank Editage (www.editage.com) for English language editing.

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/fgene.2021.679097/full#supplementary-material

Supplementary Figure 1

MEG3 expression levels in GBM and LGG in the GEPIA database.

Supplementary Figure 2

Kaplan-Meier survival curves comparing the low and high expression of MEG3 in LGG and GBM using the GEPIA database.

References

Summary

Keywords

MEG3, glioma, immune infiltration, prognosis, biomarker

Citation

Xu X, Zhong Z, Shao Y and Yi Y (2021) Prognostic Value of MEG3 and Its Correlation With Immune Infiltrates in Gliomas. Front. Genet. 12:679097. doi: 10.3389/fgene.2021.679097

Received

11 March 2021

Accepted

10 May 2021

Published

16 June 2021

Volume

12 - 2021

Edited by

Guohua Huang, Shaoyang University, China

Reviewed by

Luca Zammataro, Yale University, United States; Enrique Medina-Acosta, State University of the North Fluminense Darcy Ribeiro, Brazil

Updates

Copyright

*Correspondence: Yong Yi,

This article was submitted to Computational Genomics, a section of the journal Frontiers in Genetics

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics