Identification of hub genes and regulatory factors of glioblastoma multiforme subgroups by RNA-seq data analysis

Bioinformatics Strategies to Identify Shared Molecular Biomarkers That Link Ischemic Stroke and Moyamoya Disease with Glioblastoma

Expanding data suggest that glioblastoma is accountable for the growing prevalence of various forms of stroke formation, such as ischemic stroke and moyamoya disease. However, the underlying deterministic details are still unspecified. Bioinformatics approaches are designed to investigate the relationships between two pathogens as well as fill this study void. Glioblastoma is a form of cancer that typically occurs in the brain or spinal cord and is highly destructive. A stroke occurs when a brain region starts to lose blood circulation and prevents functioning. Moyamoya disorder is a recurrent and recurring arterial disorder of the brain. To begin, adequate gene expression datasets on glioblastoma, ischemic stroke, and moyamoya disease were gathered from various repositories. Then, the association between glioblastoma, ischemic stroke, and moyamoya was established using the existing pipelines. The framework was developed as a generalized workflow to allow for the aggregation of transcriptomic gene expression across specific tissue; Gene Ontology (GO) and biological pathway, as well as the validation of such data, are carried out using enrichment studies such as protein–protein interaction and gold benchmark databases. The results contribute to a more profound knowledge of the disease mechanisms and unveil the projected correlations among the diseases.

Allopregnanolone Promotes Migration and Invasion of Human Glioblastoma Cells through the Protein Tyrosine Kinase c-Src Activation

Glioblastomas (GBs) are the most aggressive and common primary malignant brain tumors. Steroid hormone progesterone (P4) and its neuroactive metabolites, such as allopregnanolone (3α-THP) are synthesized by neural, glial, and malignant GB cells. P4 promotes cellular proliferation, migration, and invasion of human GB cells at physiological concentrations. It has been reported that 3α-THP promotes GB cell proliferation. Here we investigated the effects of 3α-THP on GB cell migration and invasion, the participation of the enzymes involved in its metabolism (AKR1C1-4), and the role of the c-Src kinase in 3α-THP effects in GBs. 3α-THP 100 nM promoted migration and invasion of U251, U87, and LN229 human-derived GB cell lines. We observed that U251, LN229, and T98G cell lines exhibited a higher protein content of AKR1C1-4 than normal human astrocytes. AKR1C1-4 silencing did not modify 3α-THP effects on migration and invasion. 3α-THP activated c-Src protein at 10 min (U251 cells) and 15 min (U87 and LN229 cells). Interestingly, the pharmacological inhibition of c-Src decreases the promoting effects of 3α-THP on cell migration and invasion. Together, these data indicate that 3α-THP promotes GB migration and invasion through c-Src activation.

Identification of key genes as predictive biomarkers for osteosarcoma metastasis using translational bioinformatics

Background

Osteosarcoma (OS) metastasis is the most common cause of cancer-related mortality, however, no sufficient clinical biomarkers have been identified. In this study, we identified five genes to help predict metastasis at diagnosis.

Methods

We performed weighted gene co-expression network analysis (WGCNA) to identify the most relevant gene modules associated with OS metastasis. An important machine learning algorithm, the support vector machine (SVM), was employed to predict key genes for classifying the OS metastasis phenotype. Finally, we investigated the clinical significance of key genes and their enriched pathways.

Results

Eighteen modules were identified in WGCNA, among which the pink, red, brown, blue, and turquoise modules demonstrated good preservation. In the five modules, the brown and red modules were highly correlated with OS metastasis. Genes in the two modules closely interacted in protein–protein interaction networks and were therefore chosen for further analysis. Genes in the two modules were primarily enriched in the biological processes associated with tumorigenesis and development. Furthermore, 65 differentially expressed genes were identified as common hub genes in both WGCNA and protein–protein interaction networks. SVM classifiers with the maximum area under the curve were based on 30 and 15 genes in the brown and red modules, respectively. The clinical significance of the 45 hub genes was analyzed. Of the 45 genes, 17 were found to be significantly correlated with survival time. Finally, 5/17 genes, including ADAP2 (P = 0.0094), LCP2 (P = 0.013), ARHGAP25 (P = 0.0049), CD53 (P = 0.016), and TLR7 (P = 0.04) were significantly correlated with the metastatic phenotype. In vitro verification, western blotting, wound healing analyses, transwell invasion assays, proliferation assays, and colony formation assays indicated that ARHGAP25 promoted OS cell migration, invasion, proliferation, and epithelial–mesenchymal transition.

Conclusion

We identified five genes, namely ADAP2, LCP2, ARHGAP25, CD53, and TLR7, as candidate biomarkers for the prediction of OS metastasis; ARHGAP25 inhibits MG63 OS cell growth, migration, and invasion in vitro, indicating that ARHGAP25 can serve as a promising specific and prognostic biomarker for OS metastasis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02308-w.

Description of a CSF-Enriched miRNA Panel for the Study of Neurological Diseases

Background: The study of circulating miRNAs in CSF has gained tremendous attention during the last years, as these molecules might be promising candidates to be used as biomarkers and provide new insights into the disease pathology of neurological disorders. Objective: The main aim of this study was to describe an OpenArray panel of CSF-enriched miRNAs to offer a suitable tool to identify and characterize new molecular signatures in different neurological diseases. Methods: Two hundred and fifteen human miRNAs were selected to be included in the panel, and their expression and abundance in CSF samples were analyzed. In addition, their stability was studied in order to propose suitable endogenous controls for CSF miRNA studies. Results: miR-143-3p and miR-23a-3p were detected in all CSF samples, while another 80 miRNAs were detected in at least 70% of samples. miR-770-5p was the most abundant miRNA in CSF, presenting the lowest mean Cq value. In addition, miR-26b-5p, miR-335-5p and miR-92b-3p were the most stable miRNAs and could be suitable endogenous normalizers for CSF miRNA studies. Conclusions: These OpenArray plates might be a suitable and efficient tool to identify and characterize new molecular signatures in different neurological diseases and would improve the yield of miRNA detection in CSF.

Identification of novel drug resistance mechanisms by genomic and transcriptomic profiling of glioblastoma cells with mutation-activated EGFR

AIMS

Epidermal growth factor receptor (EGFR) is not only involved in carcinogenesis, but also in chemoresistance. We characterized U87.MGΔEGFR glioblastoma cells with constitutively active EGFR due to deletion at the ligand binding domain in terms of gene expression profiling and chromosomal aberrations. Wild-type U87.MG cells served as control.

MATERIALS AND METHODS

RNA sequencing and network analyses (Ingenuity Pathway Analysis) were performed to identify novel drug resistance mechanisms related to expression of mutation activated EGFR. Chromosomal aberrations were characterized by multicolor fluorescence in situ hybridization (mFISH) and array comparative genomic hybridization (aCGH).

KEY FINDINGS

U87.MGΔEGFR cells presented much more chromosomal aberrations, amplifications and deletions than wild-type U87.MG cells. Still, both cell lines were near-triploid. Numerous genes were overexpressed in U87.MGΔEGFR cells, some of which have been already linked to drug resistance. PXDN, which is associated with epithelial mesenchymal transition, was the most upregulated gene (901.8-fold). TENM1 was 331.6-fold upregulated, and it was previously reported to modulate neural development. EGFR-AS1 (161.2-fold upregulated) has been reported to increase the EGFR mRNA stability and its expression - in accordance with that of EGFR - was upregulated (85.5-fold). In addition to well-known resistance genes, numerous novel genes and genomic aberrations were identified. ANGPT2 upregulation and CPM downregulation were validated by Western blotting.

SIGNIFICANCE

Transcriptomics and genomics analyses in U87.MGΔEGFR cells unraveled a range of novel drug resistance mechanisms including apoptosis, DNA repair, ferroptosis, glutathione related gene activities, heat shock, oxidative stress, transcription factor activities, which may have important implications for future treatment strategies.

Effect of Grilled Nux Vomica on Differential RNA Expression Profile of Gastrocnemius Muscle and Toll‑Like Receptor 4 (TLR-4)/Nuclear Factor kappa B (NF-κB) Signaling in Experimental Autoimmune Myasthenia Gravis Rats

BACKGROUND Myasthenia gravis (MG) is a progressive autoimmune disorder caused by the production of antibodies directed against acetylcholine receptors (AChRs), resulting in muscle weakness and fatigue. This study aimed to explore the effect and mechanism of grilled nux vomica (GNV) in experimental autoimmune myasthenia gravis (EAMG) rats. MATERIAL AND METHODS Rat 97-116 peptides were used to mediate disease in the EAMG model in SPF female Lewis rats. The treatment groups received grilled nux vomica (75 mg/kg, 150 mg/kg, and 225 mg/kg). The autoantibody and inflammatory cytokines levels were measured by enzyme-linked immunosorbent assay (ELISA). RNA profiling was performed on high-dose and model group rats. Profiling results and TLR-4/NF-kappaB signaling were validated by q-PCR and Western blot analysis. RESULTS The results showed that GNV could attenuate the symptoms of EAMG rats. There was a decreased level of AChR-ab, IFN-γ, TNF-alpha, IL-2, IL-4, and IL-17 levels, and an increased level of TGF-ß1. In total, 235 differentially expressed genes (DEGs), consisting of 175 upregulated DEGs and 60 downregulated DEGs, were identified. Functional annotation demonstrated that DEGs were largely associated with leukocyte cell-cell adhesion, NF-kappa B signaling pathway, muscle contraction, and cardiac muscle contraction pathway. Rac2, Itgb2, Lcp2, Myl3, and Tnni1 were considered as hub genes with a higher degree value in the protein-protein interaction (PPI) network. The q-PCR and Western blot results of hub genes were consistent with RNA profiles. GNV treatment also significantly reduced the TLR-4 and NF-kappaB p65 protein expression in EAMG rats. CONCLUSIONS These results indicate that grilled nux vomica ameliorates EAMG by depressing the TLR-4/NF-kappaB signaling pathway, and hub genes may serve as potential targets for MG treatment.

Grade II/III Glioma Microenvironment Mining and Its Prognostic Merit

OBJECTIVE

The tumor microenvironment greatly influences tumor formation, invasion, and progression. The ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues) algorithm quantifies stromal and immune components in a tumor, reflecting the tumor microenvironment. This study aimed to explore key prognostic genes in a grade II/III glioma microenvironment.

METHODS

We obtained stromal/immune scores for the Cancer Genome Atlas (TCGA) grade II/III glioma cohort from the online ESTIMATE portal. The associations of stromal/immune scores with clinicopathologic characteristics and overall survival of patients with grade II/III glioma were assessed by the Mann-Whitney U test and the Kaplan-Meier method, respectively. Functional enrichment analysis and protein-protein interaction network assessments were employed to analyze differentially expressed genes (DEGs). The top 7 genes with 5 or more edges in the protein-protein interaction network were selected. For validation, CGGA grade II/III glioma data were analyzed.

RESULTS

The results showed that elevated stromal/immune/ESTIMATE score was significantly associated with poor survival of patients with TCGA grade II/III glioma. Functional enrichment analysis showed that DEGs were associated with immune cell regulation, extracellular matrix, cytokine activation, and receptor binding. The selected DEGs (interleukin-10, beta-2 microglobulin, C-C motif chemokine ligand 5, cluster of differentiation 74, human leukocyte antigen-DRA, lymphocyte cytosolic protein 2, and myxovirus resistance protein 1) showed prognostic values in patients with grade II/III glioma of the TCGA and CGGA database.

CONCLUSIONS

Stromal/immune/ESTIMATE scores have prognostic values in patients with grade II/III glioma. The selected DEGs, including interleukin-10, beta-2 microglobulin, C-C motif chemokine ligand 5, cluster of differentiation 74, human leukocyte antigen-DRA, lymphocyte cytosolic protein 2, and myxovirus resistance protein 1, associated with tumor immunity and microenvironment, have prognostic values in grade II/III glioma. Further investigation of these genes could provide novel insights into the tumor microenvironment of glioma.

Identification of key candidate genes and pathways in glioblastoma by integrated bioinformatical analysis

Glioblastoma (GBM), characterized by high morbidity and mortality, is one of the most common lethal diseases worldwide. To identify the molecular mechanisms that contribute to the development of GBM, three cohort profile datasets (GSE50161, GSE90598 and GSE104291) were integrated and thoroughly analyzed; these datasets included 57 GBM cases and 22 cases of normal brain tissue. The current study identified differentially expressed genes (DEGs), and analyzed potential candidate genes and pathways. Additionally, a DEGs-associated protein-protein interaction (PPI) network was established for further investigation. Then, the hub genes associated with prognosis were identified using a Kaplan-Meier analysis based on The Cancer Genome Atlas database. Firstly, the current study identified 378 consistent DEGs (240 upregulated and 138 downregulated). Secondly, a cluster analysis of the DEGs was performed based on functions of the DEGs and signaling pathways were analyzed using the enrichment analysis tool on DAVID. Thirdly, 245 DEGs were identified using PPI network analysis. Among them, two co-expression modules comprising of 30 and 27 genes, respectively, and 35 hub genes were identified using Cytoscape MCODE. Finally, Kaplan-Meier analysis of the hub genes revealed that the increased expression of calcium-binding protein 1 (CABP1) was negatively associated with relapse-free survival. To summarize, all enriched Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways may participate in mechanisms underlying GBM occurrence and progression, however further studies are required. CABP1 may be a key gene associated with the biological process of GBM development and may be involved in a crucial mechanism of GBM progression.

Candidate Biomarkers and Molecular Mechanism Investigation for Glioblastoma Multiforme Utilizing WGCNA

To reveal the potential molecular mechanism of glioblastoma multiforme (GBM) and provide the candidate biomarkers for GBM gene therapy. Microarray dataset GSE50161 was obtained from GEO database. The differentially expressed genes (DEGs) were identified between GBM samples and control samples, followed by the module partition analysis based on WGCNA. Then, the pathway and functional enrichment analyses of DEGs were performed. The hub genes were further investigated, followed by the survival analysis and data validation. A total of 1913 DEGs were investigated between two groups, followed by analysis of 5 modules using WGCNA. These DEGs were mainly enriched in functions like inflammatory response. The hub genes including upregulated N-Myc and STAT Interactor (NMI), Capping Actin Protein-Gelsolin Like (CAPG), and Proteasome Subunit Beta 8 (PSMB8) were revealed as potential liquid biopsy molecules for GBM diagnose. Moreover, Nucleolar and Spindle Associated Protein 1 (NUSAP1) and G Protein-Coupled Receptor 65 (GPR65) were outstanding genes in survival analysis. Our results suggested that CPNE6, HAPLN2, CMTM3, NMI, CAPG, and PSMB8 might be used as potential molecules for liquid biopsy of GBM. NUSAP1 and GPR65 might be novel prognostic targets for GBM gene therapy. Furthermore, the upregulated NMI might play an important role in GBM progression via inflammatory response.

MicroRNAs as biomarkers for human glioblastoma: progress and potential

Glioblastoma multiforme (GBM) is the most common malignant glioma. Despite innovative research efforts in tumor therapy, the outcome for most diagnosed patients remains poor; therefore, early diagnosis of GBM is the most effective method for achieving better patient outcomes. In recent years, combined research efforts including cellular, molecular, genetic, and bioinformatics methods have been used to investigate GBM, and the results show that variations in miRNA expression occur in GBM tissues and biological fluids. Some highly stable miRNAs circulate in the blood and cerebrospinal fluid (CSF) of both healthy individuals and diagnosed patients, thus raising the possibility that miRNAs may serve as novel diagnostic markers. In addition, increased understanding of the miRNA and mRNA interactions involved in GBM progression may lead to discovering predictive biomarkers, some of which are clinically relevant for targeted therapy and predicting prognosis. However, as this field is relatively new, some studies have yielded conflicting results. To progress in the field, different advanced techniques must be combined, including bioinformatics methods and molecular and cellular techniques. In addition, we must overcome the various challenges in non-invasive GBM biomarker detection. Here, we discuss the progress and potential of miRNAs as biomarkers for GBM and related signaling pathways. Studying the clinical relevance and applicability of these biomarkers may alter GBM patient diagnosis and treatment.

Glioblastoma in natalizumab-treated multiple sclerosis patients

We present two natalizumab‐treated multiple sclerosis patients who developed glioblastoma multiforme (GBM) with variable outcomes. One patient had an isocitrate dehydrogenase (IDH)‐wildtype GBM with aggressive behavior, who declined treatment and died 13 weeks after symptoms onset. The other patient underwent resection of an IDH‐mutant secondary GBM that arose from a previously diagnosed grade II astrocytoma. He is still alive 5 years after the diagnosis of GBM. JC virus was not detected in either case. Whether natalizumab played a role in the development of GBM in those patients deserves further investigation.