AUTHOR=Yu Wanli , Wu Pengfei , Wang Fang , Miao Li , Han Bo , Jiang Zhiqun 

TITLE=Construction of Novel Methylation-Driven Gene Model and Investigation of PARVB Function in Glioblastoma

JOURNAL=Frontiers in Oncology

VOLUME=11

YEAR=2021

URL=https://www.frontiersin.org/journals/oncology/articles/10.3389/fonc.2021.705547

DOI=10.3389/fonc.2021.705547

ISSN=2234-943X

ABSTRACT=<sec><title>Background</title><p>Glioblastoma multiforme (GBM) is characterized by widespread genetic and transcriptional heterogeneity. Aberrant DNA methylation plays a vital role in GBM progression by regulating gene expression. However, little is known about the role of methylation and its association with prognosis in GBM. Our aim was to explore DNA methylation-driven genes (DMDGs) and provide evidence for survival prediction and individualized treatment of GBM patients.</p></sec><sec><title>Methods</title><p>Use of the MethylMix R package identified DMDGs in GBM. The prognostic signature of DMDGs based on the risk score was constructed by multivariate Cox regression analysis. Receiver operating characteristics (ROC) curve and C-index were applied to assess the predictive performance of the DMDG prognostic signature. The predictive ability of the multigene signature model was validated in TCGA and CGGA cohorts. Finally, the role of DMDG β-Parvin (PARVB) was explored <italic>in vitro</italic>.</p></sec><sec><title>Results</title><p>The prognostic signature of DMDGs was constructed based on six genes (<italic>MDK</italic>, <italic>NMNAT3</italic>, <italic>PDPN</italic>, <italic>PARVB</italic>, <italic>SERPINB1</italic>, and <italic>UPP1</italic>). The low-risk cohort had significantly better survival than the high-risk cohort (p &lt; 0.001). The area under the curve of the ROC of the six-gene signature was 0.832, 0.927, and 0.980 within 1, 2, and 3 years, respectively. The C-index of 0.704 indicated superior specificity and sensitivity. The six-gene model has been demonstrated to be an independent prognostic factor for GBM. In addition, joint survival analysis indicated that the <italic>MDK</italic>, <italic>NMNAT3</italic>, <italic>PARVB</italic>, <italic>SERPINB1</italic>, and <italic>UPP1</italic> genes were significantly associated with prognosis and therapeutic targets for GBM. Importantly, our DMDG prognostic model was more suitable and accurate for low-grade gliomas. Finally, we verified that PARVB induced epithelial-mesenchymal transition partially through the JAK2/STAT3 pathway, which in turn promoted GBM cell proliferation, migration, and invasion.</p></sec><sec><title>Conclusion</title><p>This study demonstrated the potential value of the prognostic signature of DMDGs and provided important bioinformatic and potential therapeutic target data to facilitate individualized treatment for GBM, and to elucidate the specific mechanism by which PARVB promotes GBM progression.</p></sec>