Haoyuan Wang

A three-gene signature for prognosis in patients with MGMT promoter-methylated glioblastoma

Glioblastoma is the most malignant tumor and has high mortality rate. The methylated prompter of MGMT results in chemotherapy sensitivity for these patients. However, there are still other factors that affected the prognosis for the glioblastoma patients with similar MGMT methylation status. We developed a signature with three genes screened from the whole genome mRNA expression profile from Chinese Glioma Genome Atlas (CGGA) and RNAseq data from The Cancer Genome Atlas (TCGA). Patients with MGMT methylation in low risk group had longer survival than those in high risk group (median overall survival 1074 vs. 372 days; P = 0.0033). Moreover, the prognostic value of the signature was significant difference in cohorts stratified by MGMT methylation and chemotherapy (P=0.0473), while there is no significant difference between low and high risk group or unmethylated MGMT patients without chemotherapy. Multivariate analysis indicated that the risk score was an independent prognosis factor (P = 0.004). In conclusion, our results showed that the signature has prognostic value for patients with MGMT promoter-methylated glioblastomas based on bioinformatics analysis.

Upregulation of long noncoding RNA HOXA-AS3 promotes tumor progression and predicts poor prognosis in glioma

Long noncoding RNAs (lncRNAs) have recently emerged as new potentially promising therapeutic targets in many cancers. However, their prognostic value and biological functions associated with glioma remain to be elucidated. Here, High-throughput RNAseq was performed to detect the expression profiles of lncRNAs in 325 human glioma tissues. It was shown that a novel lncRNA HOXA-AS3 was one of the most significantly upregulated lncRNAs in glioma tissues. Quantitative PCR further verified the increased expression of HOXA-AS3 in patient samples and glioma cell lines. Uni and Multivariate Cox regression analysis revealed that HOXA-AS3 was an independent prognostic factor in glioma patients. Gene set enrichment analysis indicated that the gene sets correlated with HOXA-AS3 expression were involved in cell cycle progression and E2F targets. Functionally, HOXA-AS3 silencing resulted in proliferation arrest by altering cell cycle progression and promoting cell apoptosis, and impaired cell migration in glioma cells. Furthermore, the growth-inhibiting effect of HOXA-AS3 knockdown was also demonstrated in Xenograft mouse model. Our results highlight the important role of HOXA-AS3 in glioma progression, and indicate that HOXA-AS3 may be served as a valuable prognostic biomarker for glioma.Long noncoding RNAs (lncRNAs) have recently emerged as new potentially promising therapeutic targets in many cancers. However, their prognostic value and biological functions associated with glioma remain to be elucidated. Here, High-throughput RNAseq was performed to detect the expression profiles of lncRNAs in 325 human glioma tissues. It was shown that a novel lncRNA HOXA-AS3 was one of the most significantly upregulated lncRNAs in glioma tissues. Quantitative PCR further verified the increased expression of HOXA-AS3 in patient samples and glioma cell lines. Uni and Multivariate Cox regression analysis revealed that HOXA-AS3 was an independent prognostic factor in glioma patients. Gene set enrichment analysis indicated that the gene sets correlated with HOXA-AS3 expression were involved in cell cycle progression and E2F targets. Functionally, HOXA-AS3 silencing resulted in proliferation arrest by altering cell cycle progression and promoting cell apoptosis, and impaired cell migration in glioma cells. Furthermore, the growth-inhibiting effect of HOXA-AS3 knockdown was also demonstrated in Xenograft mouse model. Our results highlight the important role of HOXA-AS3 in glioma progression, and indicate that HOXA-AS3 may be served as a valuable prognostic biomarker for glioma.

LncRNA profile study reveals four-lncRNA signature associated with the prognosis of patients with anaplastic gliomas

Anaplastic glioma is Grade III and the median overall survival is about 37.6 months. However, there are still other factors that affect the prognosis for anaplastic glioma patients due to variable overall survival. So we screened four-lncRNA signature (AGAP2-AS1, TPT1-AS1, LINC01198 and MIR155HG) from the lncRNA expression profile from the GSE16011, CGGA and REMBRANDT datasets. The patients in low risk group had longer overall survival than high risk group (median OS 2208.25 vs. 591.30 days; P < 0.0001). Moreover, patients in the low risk group showed similar overall survival to Grade II patients (P = 0.1669), while the high risk group showed significant different to Grade IV (P = 0.0005) with similar trend. So based on the four-lncRNA, the anaplastic gliomas could be divided into grade II-like and grade IV-like groups. On the multivariate analysis, it showed the signature was an independent prognostic factor (P = 0.000). The expression of four lncRNAs in different grades showed that AGAP2-AS1, LINC01198 and MIR155HG were increased with tumor grade, while TPT1-AS1 was decreased. Knockdown of AGAP2-AS1 can inhibit the cell proliferation, migration and invasion, while increase the apoptosis cell rates in vitro. In conclusion, our results showed that the four-lncRNA signature has prognostic value for anaplastic glioma. Moreover, clinicians should conduct corresponding therapies to achieve best treatment with less side effects for two groups patients.

Role of KCNB1 in the prognosis of gliomas and autophagy modulation

Increasing evidence suggests that ion channel genes play an important role in the progression of gliomas. However, the mechanisms by which ion channel genes influence the progression of glioma are not fully understood. We identified KCNB1 as a novel ion gene, associated with malignant progression and favorable overall survival (OS) and progression-free survival (PFS) in glioma patients from three datasets (CGGA, GSE16011 and REMBRANDT). Moreover, we characterized a novel function of autophagy induction accompanied by increased apoptosis and reduced proliferation and invasion of glioma cells for KCNB1. KEGG pathway analysis and in vitro studies suggested that the ERK pathway is involved in KCNB1-mediated regulation of autophagy, which was confirmed by inhibition of KCNB1-induced autophagy by using a selective ERK1/2 inhibitor (U0126) or siERK1/2. In vivo studies showed that KCNB1 induced autophagy while inhibiting tumor growth and increasing survival. Overall, our studies define KCNB1 as a novel prognostic factor for gliomas that exerts its tumor suppressive function through autophagy induction.

Integrated Analysis of LncRNA-mRNA Co-Expression Profiles in Patients with Moyamoya Disease

Moyamoya disease (MMD) is an idiopathic disease associated with recurrent stroke. However, the pathogenesis of MMD remains unknown. Therefore, we performed long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles in blood samples from MMD patients (N = 15) and healthy controls (N = 10). A total of 880 differentially expressed lncRNAs (3649 probes) and 2624 differentially expressed mRNAs (2880 probes) were obtained from the microarrays of MMD patients and healthy controls (P < 0.05; Fold Change >2.0). Gene ontology (GO) and pathway analyses showed that upregulated mRNAs were enriched for inflammatory response, Toll-like receptor signaling pathway, chemokine signaling pathway and mitogen-activated protein kinase (MAPK) signaling pathway among others, while the downregulated mRNAs were enriched for neurological system process, digestion, drug metabolism, retinol metabolism and others. Our results showed that the integrated analysis of lncRNA-mRNA co-expression networks were linked to inflammatory response, Toll-like signaling pathway, cytokine-cytokine receptor interaction and MAPK signaling pathway. These findings may elucidate the pathogenesis of MMD, and the differentially expressed genes could provide clues to find key components in the MMD pathway.

Classifying lower grade glioma cases according to whole genome gene expression

OBJECTIVE To identify a gene-based signature as a novel prognostic model in lower grade gliomas. RESULTS A gene signature developed from HOXA7, SLC2A4RG and MN1 could segregate patients into low and high risk score groups with different overall survival (OS), and was validated in TCGA RNA-seq and GSE16011 mRNA array datasets. Receiver operating characteristic (ROC) was performed to show that the three-gene signature was more sensitive and specific than histology, grade, age, IDH1 mutation and 1p/19q co-deletion. Gene Set Enrichment Analysis (GSEA) and GO analysis showed high-risk samples were associated with tumor associated macrophages (TAMs) and highly invasive phenotypes. Moreover, HOXA7-siRNA inhibited migration and invasion in vitro, and downregulated MMP9 at the protein level in U251 glioma cells. METHODS A cohort of 164 glioma specimens from the Chinese Glioma Genome Atlas (CGGA) array database were assessed as the training group. TCGA RNA-seq and GSE16011 mRNA array datasets were used for validation. Regression analyses and linear risk score assessment were performed for the identification of the three-gene signature comprising HOXA7, SLC2A4RG and MN1. CONCLUSIONS We established a three-gene signature for lower grade gliomas, which could independently predict overall survival (OS) of lower grade glioma patients with higher sensitivity and specificity compared with other clinical characteristics. These findings indicate that the three-gene signature is a new prognostic model that could provide improved OS prediction and accurate therapies for lower grade glioma patients.

Aberrant expression of lncRNAs and mRNAs in patients with intracranial aneurysm

Intracranial aneurysm (IA) is pathological dilatations of the cerebral artery and rupture of IAs can cause subarachnoid hemorrhage, which has a high ratio of fatality and morbidity. However, the pathogenesis of IAs remains unknown. We performed long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles in IA tissues and superficial temporal arteries (STAs). A total of 4129 differentially expressed lncRNAs and 2926 differentially expressed mRNAs were obtained from the microarrays (P < 0.05). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that up-regulated mRNAs were enriched in immune response, inflammatory response, regulation of immune response and lysosome, et al; while the down-regulated mRNAs were enriched in muscle contraction, smooth muscle contraction, cGMP-PKG signaling pathway and vascular smooth muscle contraction, et al. The lncRNA-mRNA co-expression networks were represented in immune response, inflammatory response, muscle contraction and vascular smooth muscle contraction. These findings may gain insight in the pathogenesis of IAs and provide clues to find key roles for IA patients.

Co-expression of mitosis-regulating genes contributes to malignant progression and prognosis in oligodendrogliomas.

The clinical prognosis of patients with glioma is determined by tumor grades, but tumors of different subtypes with equal malignancy grade usually have different prognosis that is largely determined by genetic abnormalities. Oligodendrogliomas (ODs) are the second most common type of gliomas. In this study, integrative analyses found that distribution of TCGA transcriptomic subtypes was associated with grade progression in ODs. To identify critical gene(s) associated with tumor grades and TCGA subtypes, we analyzed 34 normal brain tissue (NBT), 146 WHO grade II and 130 grade III ODs by microarray and RNA sequencing, and identified a co-expression network of six genes (AURKA, NDC80,CENPK, KIAA0101, TIMELESS and MELK) that was associated with tumor grades and TCGA subtypes as well as Ki-67 expression. Validation of the six genes was performed by qPCR in additional 28 ODs. Importantly, these genes also were validated in four high-grade recurrent gliomas and the initial lower-grade gliomas resected from the same patients. Finally, the RNA data on two genes with the highest discrimination potential (AURKA and NDC80) and Ki-67 were validated on an independent cohort (5 NBTs and 86 ODs) by immunohistochemistry. Knockdown of AURKA and NDC80 by siRNAs suppressed Ki-67 expression and proliferation of gliomas cells. Survival analysis showed that high expression of the six genes corporately indicated a poor survival outcome. Correlation and protein interaction analysis provided further evidence for this co-expression network. These data suggest that the co-expression of the six mitosis-regulating genes was associated with malignant progression and prognosis in ODs.

Genome-wide analysis of differentially expressed lncRNAs and mRNAs in primary gonadotrophin adenomas by RNA-seq

Recently, long non-coding RNAs (lncRNAs) have received increased research interest owing to their participation via distinct mechanisms in the biological processes of nonfunctional pituitary adenomas. However, changes in the expression of lncRNAs in gonadotrophin adenoma, which is the most common nonfunctional pituitary adenomas, have not yet been reported. In this study, we performed a genome-wide analysis of lncRNAs and mRNAs obtained from gonadotrophin adenoma patients’ samples and normal pituitary tissues using RNA-seq. The differentially expressed lncRNAs and mRNAs were identified using fold-change filtering. We identified 839 lncRNAs and 1015 mRNAs as differentially expressed. Gene Ontology analysis indicated that the biological functions of differentially expressed mRNAs were related to transcription regulator activity and basic metabolic processes. Ingenuity Pathway Analysis was performed to identify 64 canonical pathways that were significantly enriched in the tumor samples. Furthermore, to investigate the potential regulatory roles of the differentially expressed lncRNAs on the mRNAs, we constructed general co-expression networks for 100 coding and 577 non-coding genes that showed significantly correlated expression patterns in tumor cohort. In particular, we built a special sub-network of co-expression involving 186 lncRNAs interacting with 15 key coding genes of the mTOR pathway, which might promote the pathogenesis of gonadotrophin tumor. This is the first study to explore the patterns of genome-wide lncRNAs expression and co-expression with mRNAs, which might contribute to the molecular pathogenesis of gonadotrophin adenoma.

A three ion channel genes-based signature predicts prognosis of primary glioblastoma patients and reveals a chemotherapy sensitive subtype.

Increasing evidence suggests that ion channels not only regulate electric signaling in excitable cells but also play important roles in the development of brain tumor. However, the roles of ion channels in glioma remain controversial. In the present study, we systematically analyzed the expression patterns of ion channel genes in a cohort of Chinese patients with glioma using RNAseq expression profiling. First, a molecular signature comprising three ion channel genes (KCNN4, KCNB1 and KCNJ10) was identified using Univariate Cox regression and two-tailed students t test conducted in overall survival (OS) and gene expression. We assigned a risk score based on three ion channel genes to each primary Glioblastoma multiforme (pGBM) patient. We demonstrated that pGBM patients who had a high risk of unfavorable outcome were sensitive to chemotherapy. Next, we screened the three ion genes-based signature in different molecular glioma subtypes. The signature showed a Mesenchymal subtype and wild-type IDH1 preference. Gene ontology (GO) analysis for the functional annotation of the signature showed that patients with high-risk scores tended to exhibit the increased expression of proteins associated with apoptosis, immune response, cell adhesion and motion and vasculature development. Gene Set Enrichment Analysis (GSEA) results showed that pathways associated with negative regulation of programmed cell death, cell proliferation and locomotory behavior were highly expressed in the high-risk group. These results suggest that ion channel gene expression could improve the subtype classification in gliomas at the molecular level. The findings in the present study have been validated in two independent cohorts.