MET fusions and splicing variants is a strong adverse prognostic factor in astrocytoma, isocitrate dehydrogenase mutant.

Liu et al. describe the adverse prognostic role of MET fusions and splicing variants in astrocytoma, isocitrate dehydrogenase mutant. On this basis, MET fusions and splicing variants was suggested to be a biomarker for the diagnosis of high-grade astrocytoma, isocitrate dehydrogenase mutant.

circRNADisease v2.0: an updated resource for high-quality experimentally supported circRNA-disease associations.

circRNADisease v2.0 is an enhanced and reliable database that offers experimentally verified relationships between circular RNAs (circRNAs) and various diseases. It is accessible at http://cgga.org.cn/circRNADisease/ or http://cgga.org.cn:9091/circRNADisease/. The database currently includes 6998 circRNA-disease entries across multiple species, representing a remarkable 19.77-fold increase compared to the previous version. This expansion consists of a substantial rise in the number of circRNAs (from 330 to 4246), types of diseases (from 48 to 330) and covered species (from human only to 12 species). Furthermore, a new section has been introduced in the database, which collects information on circRNA-associated factors (genes, proteins and microRNAs), molecular mechanisms (molecular pathways), biological functions (proliferation, migration, invasion, etc.), tumor and/or cell line and/or patient-derived xenograft (PDX) details, and prognostic evidence in diseases. In addition, we identified 7 159 865 relationships between mutations and circRNAs among 30 TCGA cancer types. Due to notable enhancements and extensive data expansions, the circRNADisease 2.0 database has become an invaluable asset for both clinical practice and fundamental research. It enables researchers to develop a more comprehensive understanding of how circRNAs impact complex diseases.

Molecular and clinical characterization of PTRF in glioma via 1,022 samples.

Polymerase I and transcript release factor (PTRF) plays a role in the regulation of gene expression and the release of RNA transcripts during transcription, which have been associated with various human diseases. However, the role of PTRF in glioma remains unclear. In this study, RNA sequencing (RNA-seq) data (n = 1022 cases) and whole-exome sequencing (WES) data (n = 286 cases) were used to characterize the PTRF expression features. Gene ontology (GO) functional enrichment analysis was used to assess the biological implication of changes in PTRF expression. As a result, the expression of PTRF was associated with malignant progression in gliomas. Meanwhile, somatic mutational profiles and copy number variations (CNV) revealed the glioma subtypes classified by PTRF expression showed distinct genomic alteration. Furthermore, GO functional enrichment analysis suggested that PTRF expression was associated with cell migration and angiogenesis, particularly during an immune response. Survival analysis confirmed that a high expression of PTRF is associated with a poor prognosis. In summary, PTRF may be a valuable factor for the diagnosis and treatment target of glioma.

Genomic profiling and prognostic factors of H3 K27M-mutant spinal cord diffuse glioma.

H3 K27-altered diffuse midline glioma is a highly lethal pediatric-type tumor without efficacious treatments. Recent findings have highlighted the heterogeneity among diffuse midline gliomas with different locations and ages. Compared to tumors located in the brain stem and thalamus, the molecular and clinicopathological features of H3 K27-altered spinal cord glioma are still largely elusive, thus hindering the accurate management of patients. Here we aimed to characterize the clinicopathological and molecular features of H3 K27M-mutant spinal cord glioma in 77 consecutive cases. We found that the H3 K27M-mutant spinal cord glioma, with a median age of 35 years old, had a significantly better prognosis than H3 K27M-mutant brain tumors. We noticed a molecular heterogeneity of H3 K27M-mutant spinal cord astrocytoma via targeted sequencing with 34 cases. TP53 mutation which occurred in 58.8% of cases is mutually exclusive with PPM1D (26%) and NF1 (44%) mutations. The TP53-mutant cases had a significantly higher number of copy number variants (CNV) and a marginally higher proportion of pediatric patients (age at diagnosis <18 years old, p = 0.056). Cox regression and Kaplan-Meier curve analysis showed that the higher number of CNV events (≥3), chromosome (Chr) 9p deletion, Chr 10p deletion, ATRX mutation, CDK6 amplification, and retinoblastoma protein (RB) pathway alteration are associated with worse survival. Cox regression analysis with clinicopathological features showed that glioblastoma histological type and a high Ki-67 index (>10%) are associated with a worse prognosis. Interestingly, the histological type, an independent prognostic factor in multivariate Cox regression, can also stratify molecular features of H3 K27M-mutant spinal cord glioma, including the RB pathway, KRAS/PI3K pathway, and chromosome arms CNV. In conclusion, although all H3 K27M-mutant spinal cord diffuse glioma were diagnosed as WHO Grade 4, the histological type, molecular features representing chromatin instability, and molecular alterations associated with accelerated cell proliferative activity should not be ignored in clinical management.

WES data from 286 diffuse gliomas under the 2021 WHO Classification of Tumors of the Central Nervous System.

Diffuse gliomas (DGs) are the most common and lethal primary neoplasms in the central nervous system. The latest 2021 World Health Organization (WHO) Classification of Tumors of the Central Nervous System (CNS) was published in 2021, immensely changing the approach to diagnosis and decision making. As a part of the Chinese Glioma Genome Atlas (CGGA) project, our aim was to provide genomic profiling of gliomas in a Chinese cohort. Two hundred eighty six gliomas with different grades were collected over the last decade. Using the Illumina HiSeq platform, over 75.8 million high-quality 150 bp paired-end reads were generated per sample, yielding a total of 43.4 billion reads. We also collected each patient's clinical and pathological information and used it to annotate their genetic data. All patients were diagnosed and classified by neuro-pathologist under the 2021 WHO classification. This dataset provides an important reference for researchers and will significantly advance our understanding of gliomas.

YTHDF2 facilitates UBXN1 mRNA decay by recognizing METTL3-mediated m6A modification to activate NF-κB and promote the malignant progression of glioma.

BACKGROUND: The prognosis for diffuse gliomas is very poor and the mechanism underlying their malignant progression remains unclear. Here, we aimed to elucidate the role and mechanism of the RNA N6,2'-O-dimethyladenosine (m6A) reader, YTH N6-methyladenosine RNA binding protein 2 (YTHDF2), in regulating the malignant progression of gliomas. METHODS: YTHDF2 mRNA levels and functions were assessed using several independent datasets. Western blotting, quantitative polymerase chain reaction, and immunohistochemistry were used to evaluate the expression levels of YTHDF2 and other molecules in human and mouse tumor tissues and cells. Knockdown and overexpression were used to evaluate the effects of YTHDF2, methyltransferase-like 3 (METTL3), and UBX domain protein 1 (UBXN1) on glioma malignancy in cell and orthotopic xenograft models. RNA immunoprecipitation (RIP), methylated RIP, and RNA stability experiments were performed to study the mechanisms underlying the oncogenic role of YTHDF2. RESULTS: YTHDF2 expression was positively associated with a higher malignant grade and molecular subtype of glioma and poorer prognosis. YTHDF2 promoted the malignant progression of gliomas in both in vitro and in vivo models. Mechanistically, YTHDF2 accelerated UBXN1 mRNA degradation via METTL3-mediated m6A, which, in turn, promoted NF-κB activation. We further revealed that UBXN1 overexpression attenuated the oncogenic effect of YTHDF2 overexpression and was associated with better survival in patients with elevated YTHDF2 expression. CONCLUSIONS: Our findings confirmed that YTHDF2 promotes the malignant progression of gliomas and revealed important insight into the upstream regulatory mechanism of NF-κB activation via UBXN1 with a primary focus on m6A modification.

METTL3 enhances the stability of MALAT1 with the assistance of HuR via m6A modification and activates NF-κB to promote the malignant progression of IDH-wildtype glioma.

Understanding the role of N6-methyladenosine (m6A) in tumorigenesis and stem cell maintenance is an emerging field in glioma research. However, it is necessary to study the function of m6A in IDH-mutation and IDH-wildtype gliomas separately. Here, we aimed to elucidate the role and mechanism of the m6A writer METTL3 in regulating the malignant progression of IDH-wildtype gliomas. We demonstrated that METTL3 expression is positively associated with a higher malignant grade and poorer prognosis of IDH-wildtype gliomas but not IDH-mutant gliomas. METTL3 could also promote the malignant progression of gliomas in both in vitro and in vivo models. Mechanistically, METTL3 upregulated MALAT1 expression by enhancing its stability via m6A modification. We further revealed that HuR was essential for METTL3-mediated MALAT1 stabilization, and upregulated MALAT1 subsequently activated NF-κB. Taken together, our findings confirmed that METTL3 promoted the malignant progression of IDH-wildtype gliomas and revealed important insight into the upstream regulatory mechanism of MALAT1 and NF-κB with a primary focus on m6A modification.

Chinese Glioma Genome Atlas (CGGA): A Comprehensive Resource with Functional Genomic Data from Chinese Glioma Patients.

Gliomas are the most common and malignant intracranial tumors in adults. Recent studies have revealed the significance of functional genomics for glioma pathophysiological studies and treatments. However, access to comprehensive genomic data and analytical platforms is often limited. Here, we developed the Chinese Glioma Genome Atlas (CGGA), a user-friendly data portal for the storage and interactive exploration of cross-omics data, including nearly 2000 primary and recurrent glioma samples from Chinese cohort. Currently, open access is provided to whole-exome sequencing data (286 samples), mRNA sequencing (1018 samples) and microarray data (301 samples), DNA methylation microarray data (159 samples), and microRNA microarray data (198 samples), and to detailed clinical information (age, gender, chemoradiotherapy status, WHO grade, histological type, critical molecular pathological information, and survival data). In addition, we have developed several tools for users to analyze the mutation profiles, mRNA/microRNA expression, and DNA methylation profiles, and to perform survival and gene correlation analyses of specific glioma subtypes. This database removes the barriers for researchers, providing rapid and convenient access to high-quality functional genomic data resources for biological studies and clinical applications. CGGA is available at http://www.cgga.org.cn.

Transcriptional Characteristics of IDH-Wild Type Glioma Subgroups Highlight the Biological Processes Underlying Heterogeneity of IDH-Wild Type WHO Grade IV Gliomas.

Isocitric dehydrogenase (IDH)-wild type diffuse gliomas, which have a poorer prognosis than their IDH-mutant counterparts, are also accompanied with high heterogeneity. Here, we aimed to identify the key biological processes associated with the three groups of IDH-wild type diffuse gliomas in 323 patients. By The Consortium to Inform Molecular and Practical Approaches to CNS Tumor Taxonomy (cIMPACT-NOW) update 3 recommendation, the three groups are Group A, diffuse astrocytic glioma, World Health Organization (WHO) grade II/III; Group B, diffuse astrocytic glioma, with one (or more) of the three genetic alterations: TERT promoter mutation, EGFR gene amplification, gain of chromosome 7 combined with loss of chromosome 10, WHO grade IV; and Group C, glioblastoma, WHO grade IV. Consistent with their histologic and genetic molecular features, we successfully identified that biological activities associated with "cell cycle" and "cell mitosis" are significantly elevated in Group B compared with Group A; microenvironment-related hallmarks "angiogenesis" and "hypoxia," and biological processes of "extracellular matrix," "immune response," and "positive regulation of transcriptional activities" were more enriched in Group C than Group B. We also constructed a nine-gene signature from differentially expressed genes among the three groups to further stratify the WHO grade IV gliomas (Groups B and C) whose survival cannot be clearly stratified by current classification systems. This signature was an independent prognosis factor for WHO grade IV gliomas and had better prognostic value than other known factors in both training and validation dataset. In addition, the signature risk score was positively correlated with the amount of infiltrated immune cells, expression of immune checkpoints, and the genes enriched in biological processes of "immune response," "cell cycle," and "extracellular matrix." The bioinformatic analysis results were also validated by immunohistochemistry and patient-derived cell proliferation assay. Overall, our findings revealed the key biological processes underlying the new classifications of IDH-wild type diffuse glioma. Meanwhile, we constructed a signature, which could properly stratify the prognosis, cell proliferation activates, extracellular matrix-mediated biological activities, and immune-microenvironment of IDH-wild type WHO grade IV gliomas.

Effect of Zhenxin Xingshui Yizhi Fang on Aß25-35 induced expression of related transporters in HBMEC cell model.

ETHNOPHARMACOLOGICAL RELEVANCE: Aß (ß-amyloid) deposition and abnormal transport were suggested to be risk factors for Alzheimer's disease (AD). Zhenxin Xingshui Yizhi Fang (XSF), an ancient prescription in traditional Chinese medicine, was first recorded in Qianjin Yifang for treating palpitation, hypnosia, amnesia. It is reported that XSF could improve mice learning memory ability, reduce the deposition of senile plaques in hippocampus of rat brain. In this study, the neuroprotective effect of XSF against Aß25-35-induced apoptosis in cultured human brain microvascular endothelial cells (HBMEC) and its potential mechanism were investigated. MATERIALS AND METHODS: HBMEC cells were treated with Aß25-35 to established neurotoxic cell model. After that, the cells were treated with 125, 250, 500 µg/mL XSF to observe the protective effect. The viability of HBMEC cells were evaluated by MTT assay, the Aß25-35-induced apoptosis was characterized by Hoechst-33258 and the activity of cysteinyl aspartate specific proteinase-3. The expression level of Aß1-42 in cells induced by Aß25-35 was measured by human Aß1-42 kit. Protein and mRNA expression levels of advanced glycation end products (RAGE), low density lipoprotein receptor-related protein 1 (LRP1), glucose transporter 1 and 3 (GLUT1 and GLUT3) were assayed by capillary electrophoresis immunoassay and quantitative real-time polymerase chain reaction analyses. RESULTS: In Aß25-35 induced neurotoxic cells, the percentage of apoptotic cells, the concentration of Aß1-42 and CASPASE-3 activity, protein and mRNA expression levels of RAGE increased significantly, but that of LRP1, GLUT1 and GLUT3 significantly decreased. XSF could inhibit the apoptotic of cells, reduced the concentration of Aß1-42 and CASPASE-3 expression, downregulate RAGE and upregulate LRP1, GLUT1 and GLUT3 expression. CONCLUSION: The results suggest that XSF can reduce the cytotoxicity of HBMEC induced by Aß25-35, inhibit apoptosis, and regulate the transport of Aß on BBB and energy metabolism disorder in HBMEC.

Siglecs, Novel Immunotherapy Targets, Potentially Enhance The Effectiveness of Existing Immune Checkpoint Inhibitors in Glioma Immunotherapy.

BACKGROUND: Inhibitors of immune checkpoints have shown little effect in clinical trials involving glioma patients. Here, we explored novel targets for use in future treatments. Previous studies showed the sialic acid-binding Ig-like lectin (Siglec) family to have a specific role in immunosuppression. We aimed to study the characteristics and immune function of Siglec family members. METHODS: Transcriptome data from 1024 glioma samples and 1551 glioma single cells were used in our study. Clinical and molecular pathology information was also included. Statistical, bioinformatical methods, and single-cell sequencing analysis were applied to investigate the role of Siglec family members. RESULTS: Siglecs-5, -7, -9, and -16 showed a significant correlation with immunosuppression in glioma. They are typically expressed in higher grade, IDH-wildtype, and mesenchymal subtype gliomas. Siglec-5, -7, and -9 had a similar immune function to TIM-3, while Siglec-16 was similar to PD-L1, suppressing tumor immunity via different mechanisms. Joint use of Siglec-inhibitors and immune checkpoint inhibitors could prolong the survival of glioma patients. CONCLUSION: Siglec-5, -7, -9, and -16 suppressed tumor immunity in different ways. Joint usage of inhibitors may be an effective means to improve the efficacy of glioma immunotherapy.

A Novel DNA Methylation-Based Signature Can Predict the Responses of MGMT Promoter Unmethylated Glioblastomas to Temozolomide.

Glioblastoma (GBM) is the most malignant glioma, with a median overall survival (OS) of 14-16 months. Temozolomide (TMZ) is the first-line chemotherapy drug for glioma, but whether TMZ should be withheld from patients with GBMs that lack O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation is still under debate. DNA methylation profiling holds great promise for further stratifying the responses of MGMT promoter unmethylated GBMs to TMZ. In this study, we studied 147 TMZ-treated MGMT promoter unmethylated GBM, whose methylation information was obtained from the HumanMethylation27 (HM-27K) BeadChips (n = 107) and the HumanMethylation450 (HM-450K) BeadChips (n = 40) for training and validation, respectively. In the training set, we performed univariate Cox regression and identified that 3,565 CpGs were significantly associated with the OS of the TMZ-treated MGMT promoter unmethylated GBMs. Functional analysis indicated that the genes corresponding to these CpGs were enriched in the biological processes or pathways of mitochondrial translation, cell cycle, and DNA repair. Based on these CpGs, we developed a 31-CpGs methylation signature utilizing the least absolute shrinkage and selection operator (LASSO) Cox regression algorithm. In both training and validation datasets, the signature identified the TMZ-sensitive GBMs in the MGMT promoter unmethylated GBMs, and only the patients in the low-risk group appear to benefit from the TMZ treatment. Furthermore, these identified TMZ-sensitive MGMT promoter unmethylated GBMs have a similar OS when compared with the MGMT promoter methylated GBMs after TMZ treatment in both two datasets. Multivariate Cox regression demonstrated the independent prognostic value of the signature in TMZ-treated MGMT promoter unmethylated GBMs. Moreover, we also noticed that the hallmark of epithelial-mesenchymal transition, ECM related biological processes and pathways were highly enriched in the MGMT unmethylated GBMs with the high-risk score, indicating that enhanced ECM activities could be involved in the TMZ-resistance of GBM. In conclusion, our findings promote our understanding of the roles of DNA methylation in MGMT umethylated GBMs and offer a very promising TMZ-sensitivity predictive signature for these GBMs that could be tested prospectively.

RNA processing genes characterize RNA splicing and further stratify lower-grade glioma.

BACKGROUND: Aberrant expression of RNA processing genes may drive the alterative RNA profile in lower-grade gliomas (LGGs). Thus, we aimed to further stratify LGGs based on the expression of RNA processing genes. METHODS: This study included 446 LGGs from The Cancer Genome Atlas (TCGA, training set) and 171 LGGs from the Chinese Glioma Genome Atlas (CGGA, validation set). The least absolute shrinkage and selection operator (LASSO) Cox regression algorithm was conducted to develop a risk-signature. The receiver operating characteristic (ROC) curves and Kaplan-Meier curves were used to study the prognosis value of the risk-signature. RESULTS: Among the tested 784 RNA processing genes, 276 were significantly correlated with the OS of LGGs. Further LASSO Cox regression identified a 19-gene risk-signature, whose risk score was also an independently prognosis factor (P<0.0001, multiplex Cox regression) in the validation dataset. The signature had better prognostic value than the traditional factors "age", "grade" and "WHO 2016 classification" for 3- and 5-year survival both two datasets (AUCs > 85%). Importantly, the risk-signature could further stratify the survival of LGGs in specific subgroups of WHO 2016 classification. Furthermore, alternative splicing events for genes such as EGFR and FGFR were found to be associated with the risk score. mRNA expression levels for genes, which participated in cell proliferation and other processes, were significantly correlated to the risk score. CONCLUSIONS: Our results highlight the role of RNA processing genes for further stratifying the survival of patients with LGGs and provide insight into the alternative splicing events underlying this role.

Survivin is a prognostic indicator in glioblastoma and may be a target of microRNA-218.

Accumulating evidence has revealed that survivin expression is associated with a malignant phenotype and poor prognosis in glioma. Survivin is also a potential target of microRNA (miRNA/miR)-218. The aim of the present study was to investigate the expression and function of survivin in glioblastoma, and to examine the association between survivin and miR-218. For that purpose, survivin mRNA levels were analyzed in 144 frozen samples of glioblastoma using whole-genome RNA sequencing. In vitro cell proliferation, migration, invasion and apoptosis assays were performed, and survivin expression was detected by western blotting. The results revealed that the mRNA expression levels of survivin were negatively and significantly associated with overall survival in glioblastoma. Further in vitro analyses suggested that miR-218 may inhibit the expression of survivin. Expression of miR-218 in the LN229 cell line was significantly lower than that in the immortalized human gliocyte HEB cell line. miR-218 markedly inhibited tumor cell proliferation, migration and invasion capacities, and decreased apoptosis. miR-218 also inhibited the expression of survivin. These results indicated that survivin may be a target of miR-218 and could serve as a predictive biomarker.

Systematically characterize the clinical and biological significances of 1p19q genes in 1p/19q non-codeletion glioma.

1p/19q codeletion, which leads to the abnormal expression of 1p19q genes in oligodendroglioma, is associated with chemosensitivity and favorable prognosis. Here, we aimed to explore the clinical implications of 1p19q gene expression in 1p/19q non-codel gliomas. We analyzed expression of 1p19q genes in 668 1p/19q non-codel gliomas obtained from The Cancer Genome Atlas (n = 447) and the Chinese Glioma Genome Atlas (n = 221) for training and validation, respectively. The expression of 1p19q genes was significantly correlated with the clinicopathological features and overall survival of 1p/19q non-codel gliomas. Then, we derived a risk signature of 25 selected 1p19q genes that not only had prognosis value in total 1p/19q non-codel gliomas but also had prognosis value in stratified gliomas. The prognosis value of the risk signature was superior than known clinicopathological features in 1p/19q non-codel gliomas and was also highly associated with the following features: loss of CDKN2A/B copy number in mutant-IDH-astrocytoma; telomerase reverse transcriptase (TERT) promoter mutation, combined chromosome 7 gain/chromosome 10 loss and epidermal growth factor receptor amplification in wild-type-IDH-astrocytoma; classical and mesenchymal subtypes in glioblastoma. Furthermore, genes enriched in the biological processes of cell division, extracellular matrix, angiogenesis significantly correlated to the signature risk score, and this is also supported by the immunohistochemistry and cell biology experiments. In conclusion, the expression profile of 1p19q genes is highly associated with the malignancy and prognosis of 1p/19q non-codel gliomas. A 25-1p19q-gene signature has powerfully predictive value for both malignant molecular pathological features and prognosis across distinct subgroups of 1p/19q non-codel gliomas.

Systematically profiling the expression of eIF3 subunits in glioma reveals the expression of eIF3i has prognostic value in IDH-mutant lower grade glioma.

BACKGROUND: Abnormal expression of the eukaryotic initiation factor 3 (eIF3) subunits plays critical roles in tumorigenesis and progression, and also has potential prognostic value in cancers. However, the expression and clinical implications of eIF3 subunits in glioma remain unknown. METHODS: Expression data of eIF3 for patients with gliomas were obtained from the Chinese Glioma Genome Atlas (CGGA) (n = 272) and The Cancer Genome Atlas (TCGA) (n = 595). Cox regression, the receiver operating characteristic (ROC) curves and Kaplan-Meier analysis were used to study the prognostic value. Gene oncology (GO) and gene set enrichment analysis (GSEA) were utilized for functional prediction. RESULTS: In both the CGGA and TCGA datasets, the expression levels of eIF3d, eIF3e, eIF3f, eIF3h and eIF3l highly were associated with the IDH mutant status of gliomas. The expression of eIF3b, eIF3i, eIF3k and eIF3m was increased with the tumor grade, and was associated with poorer overall survival [All Hazard ratio (HR) > 1 and P < 0.05]. By contrast, the expression of eIF3a and eIF3l was decreased in higher grade gliomas and was associated with better overall survival (Both HR < 1 and P < 0.05). Importantly, the expression of eIF3i (located on chromosome 1p) and eIF3k (Located on chromosome 19q) were the two highest risk factors in both the CGGA [eIF3i HR = 2.068 (1.425-3.000); eIF3k HR = 1.737 (1.166-2.588)] and TCGA [eIF3i HR = 1.841 (1.642-2.064); eIF3k HR = 1.521 (1.340-1.726)] databases. Among eIF3i, eIF3k alone or in combination, the expression of eIF3i was the more robust in stratifying the survival of glioma in various pathological subgroups. The expression of eIF3i was an independent prognostic factor in IDH-mutant lower grade glioma (LGG) and could also predict the 1p/19q codeletion status of IDH-mutant LGG. Finally, GO and GSEA analysis showed that the elevated expression of eIF3i was significantly correlated with the biological processes of cell proliferation, mRNA processing, translation, T cell receptor signaling, NF-κB signaling and others. CONCLUSIONS: Our study reveals the expression alterations during glioma progression, and highlights the prognostic value of eIF3i in IDH-mutant LGG.

Molecular and clinical characterization of TMEM71 expression at the transcriptional level in glioma.

BACKGROUND: Glioma is the most common and aggressive type of primary brain tumor in adults. Although radiotherapy and chemotherapy are used in the treatment of glioma, survival remains unsatisfactory. Chemoresistance is one of the primary reasons for the poor prognosis of glioma. Several studies have demonstrated that glioma stem cells (GSC) may be one of the reasons for chemoresistance. In this article, we attempt to search for a new biomarker related to GSC and chemoresistance in glioma. METHODS: We used three datasets (GSE23806, COSMIC, and CGGA) to search for the genes related to GSC, temozolomide (TMZ) resistance, and overall survival. The selected gene was investigated with respect to the relationship between mRNA levels and clinical characteristics in the CGGA and TCGA dataset. Gene ontology (GO) analysis was used for bioinformatics analysis. Kaplan-Meier survival analysis and Cox regression analysis were used for survival analysis. RESULTS: The transmembrane protein 71 (TMEM71) gene was selected for further research. TMEM71 was highly expressed in GSCs and TMZ-resistant cells. The TMEM71 mRNA levels increased with increasing grades of glioma. In IDH-wild-type and MGMT-unmethylated samples, TMEM71 was overexpressed. The TMEM71 transcript levels were also increased significantly in mesenchymal subtype gliomas. GO analysis demonstrated that TMEM71 was related to the immune and inflammatory responses, cell proliferation, cell migration, chemotaxis, and the response to drugs. Specifically, PD-1, PD-L1, TIM-3, and B7-H3 were tightly associated with TMEM71 expression. This result indicates that TMEM71 may play an important role in the immune response. More importantly, high expression of TMEM71 was correlated with short survival time in both glioma and glioblastoma patients. CONCLUSION: In summary, TMEM71 expression was increased in GBM and associated with immune response. Our study suggests that TMEM71 may function as an oncogene and serve as a new effective therapeutic target for the treatment of glioma.