Analysis and identification of novel biomarkers involved in neuroblastoma via integrated bioinformatics

Metabolomics-transcriptomics joint analysis: unveiling the dysregulated cell death network and developing a diagnostic model for high-grade neuroblastoma

High-grade neuroblastoma (HG-NB) exhibits a significantly diminished survival rate in comparison to low-grade neuroblastoma (LG-NB), primarily attributed to the mechanism of HG-NB is unclear and the lacking effective therapeutic targets and diagnostic model. Therefore, the current investigation aims to study the dysregulated network between HG-NB and LG-NB based on transcriptomics and metabolomics joint analysis. Meanwhile, a risk diagnostic model to distinguish HG-NB and LG-NB was also developed. Metabolomics analysis was conducted using plasma samples obtained from 48 HG-NB patients and 36 LG-NB patients. A total of 39 metabolites exhibited alterations, with 20 showing an increase and 19 displaying a decrease in HG-NB. Additionally, transcriptomics analysis was performed on NB tissue samples collected from 31 HG-NB patients and 20 LG-NB patients. Results showed that a significant alteration was observed in a total of 1,199 mRNAs in HG-NB, among which 893 were upregulated while the remaining 306 were downregulated. In particular, the joint analysis of both omics data revealed three aberrant pathways, namely the cAMP signaling pathway, PI3K-Akt signaling pathway, and TNF signaling pathway, which were found to be associated with cell death. Notably, a diagnostic model for HG-NB risk classification was developed based on the genes MGST1, SERPINE1, and ERBB3 with an area under the receiver operating characteristic curve of 0.915. In the validation set, the sensitivity and specificity were determined to be 75.0% and 80.0%, respectively.

Melatonin suppresses tumor proliferation and metastasis by targeting GATA2 in endometrial cancer

Endometrial cancer (EC) is a reproductive system disease that occurs in perimenopausal and postmenopausal women. However, its etiology is unclear. Melatonin (MT) has been identified as a therapeutic agent for EC; however, its exact mechanism remains unclear. In the present study, we determined that GATA-binding protein 2 (GATA2) is expressed at low levels in EC and regulated by MT. MT upregulates the expression of GATA2 through MT receptor 1A (MTNR1A), whereas GATA2 can promote the expression of MTNR1A by binding to its promoter region. In addition, in vivo and in vitro experiments showed that MT inhibited the proliferation and metastasis of EC cells by upregulating GATA2 expression. The protein kinase B (AKT) pathway was also affected. In conclusion, these findings suggest that MT and GATA2 play significant roles in EC development.

P3H4 and PLOD1 expression associates with poor prognosis in bladder cancer

PURPOSE

The prolyl 3-hydroxylase family member 4 gene (P3H4) is involved in the development of human cancers. The association of P3H4 with bladder cancer (BC) prognosis is unclear. This study aimed to analyze the association of P3H4 with BC prognosis.

METHODS

RNA-Seq data were downloaded from The Cancer Genome Atlas project and BC microarray datasets (GSE13507, GSE31684, and GSE32548) were downloaded from the Gene Expression Omnibus database. We analyzed the differences in P3H4 expression levels between BC tumors and non-tumor tissues and between samples with different clinical information. The association of P3H4 and P3H4-related genes with BC prognosis and the possibility of using P3H4 expression as a prognostic biomarker in BC patients were also analyzed. RevMan was used to perform the meta-analysis.

RESULTS

P3H4 was upregulated in BC tissues compared with the adjacent non-tumor tissues (p = 4.06e-08). Univariate Cox regression analysis and meta-analysis showed that high P3H4 expression level contributed to a poor BC prognosis (Hazard ratio, HR = 1.348, 95% CI 1.140-1.594, p = 4.89e-04; meta-analysis: HR = 1.45, 95% CI 1.10-1.91; p = 9.00e-03). Among the genes related to P3H4, the PLOD1 gene was closely associated with P3H4 expression (r = 0.620, p = 2.49e-44). Also, a meta-analysis showed that PLOD1 expression was associated with a poor prognosis in BC patients (HR = 1.77, 95% CI 1.31-2.38; p = 2.00e-04).

CONCLUSIONS

The P3H4 and PLOD1 genes might be used as reliable prognostic biomarkers for BC.

Deep analysis of neuroblastoma core regulatory circuitries using online databases and integrated bioinformatics shows their pan-cancer roles as prognostic predictors

AIM

Neuroblastoma is a heterogeneous childhood cancer derived from the neural crest. The dual cell identities of neuroblastoma include Mesenchymal (MES) and Adrenergic (ADRN). These identities are conferred by a small set of tightly-regulated transcription factors (TFs) binding super enhancers, collectively forming core regulatory circuitries (CRCs). The purpose of this study was to gain a deep understanding of the role of MES and ADRN TFs in neuroblastoma and other cancers as potential indicators of disease prognosis, progression, and relapse.

METHODS

To that end, we first investigated the expression and mutational profile of MES and ADRN TFs in neuroblastoma. Moreover, we established their correlation with neuroblastoma risk groups and overall survival while establishing their extended networks with long non-coding RNAs (lncRNAs). Furthermore, we analysed the pan-cancer expression and mutational profile of these TFs and their correlation with patient survival and finally their network connectivity, using a panel of bioinformatic tools including GEPIA2, human pathology atlas, TIMER2, Omicsnet, and Cytoscape.

RESULTS

We show the association of multiple MES and ADRN TFs with neuroblastoma risk groups and overall survival and find significantly higher expression of various MES and ADRN TFs compared to normal tissues and their association with overall survival and disease-free survival in multiple cancers. Moreover, we report the strong correlation of the expression of these TFs with the infiltration of stromal and immune cells in the tumour microenvironment and with stemness and metastasis-related genes. Furthermore, we reveal extended pan-cancer networks comprising these TFs that influence the tumour microenvironment and metastasis and may be useful indicators of cancer prognosis and patient survival.

CONCLUSION

Our meta-analysis shows the significance of MES and ADRN TFs as indicators of patient prognosis and the putative utility of these TFs as potential novel biomarkers.

Identify potential miRNA-mRNA regulatory networks contributing to high-risk neuroblastoma

Neuroblastoma (NB) is a common tumor in children, usually in the retroperitoneum. After various treatments, low- and intermediate-risk patients have achieved good results, but the prognosis of high-risk patients is still very poor. Therefore, it is necessary to find new effective targets for the treatment of high-risk patients. In this study, comprehensive bioinformatics analysis was used to identify the differentially expressed genes (DEG and DEM) between high-risk patients and non-high-risk patients, and it was identified that ADRB2 may affect the survival status of high-risk patients due to miR -30a-5p regulation. The GSE49710, GSE73517, and GSE121513 datasets were downloaded from the Gene Expression Synthesis (GEO) database, and DEG and DEM were selected. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were applied to the selected DEGs. The STRING database and Cytoscape software were used to construct protein-protein interaction (PPI) networks and perform modular analysis of the DEGs. The TARGET data set containing information on overall survival days were used for the prognostic analysis of central genes. We identified a total of 255 DEGs from GSE49710 and GSE73517, and 193 DEMs from GSE121513. We identified the 5 most important central genes from the PPI network, performed a prognostic analysis in the target data set, and verified their expression using RT-qPCR to select the most important ADRB2 gene to predict miRNA. Integrating the differential miRNA predicted by miRDB and miRSystem and GSE121513 between the targeted miRNA and the prognosis, miR-30a-5p was finally identified as the targeted miRNA of ADRB2.

Screening and identification of potential biomarkers and therapeutic drugs in melanoma via integrated bioinformatics analysis

Melanoma is a highly aggressive malignant skin tumor with a high rate of metastasis and mortality. In this study, a comprehensive bioinformatics analysis was used to clarify the hub genes and potential drugs. Download the GSE3189, GSE22301, and GSE35388 microarray datasets from the Gene Expression Omnibus (GEO), which contains a total of 33 normal samples and 67 melanoma samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) approach analyze DEGs based on the DAVID. Use STRING to construct protein-protein interaction network, and use MCODE and cytoHubba plug-ins in Cytoscape to perform module analysis and identified hub genes. Use Gene Expression Profile Interactive Analysis (GEPIA) to assess the prognosis of genes in tumors. Finally, use the Drug-Gene Interaction Database (DGIdb) to screen targeted drugs related to hub genes. A total of 140 overlapping DEGs were identified from the three microarray datasets, including 59 up-regulated DEGs and 81 down-regulated DEGs. GO enrichment analysis showed that these DEGs are mainly involved in the biological process such as positive regulation of gene expression, positive regulation of cell proliferation, positive regulation of MAP kinase activity, cell migration, and negative regulation of the apoptotic process. The cellular components are concentrated in the membrane, dendritic spine, the perinuclear region of cytoplasm, extracellular exosome, and membrane raft. Molecular functions include protein homodimerization activity, calmodulin-binding, transcription factor binding, protein binding, and cytoskeletal protein binding. KEGG pathway analysis shows that these DEGs are mainly related to protein digestion and absorption, PPAR signaling pathway, signaling pathways regulating stem cells' pluripotency, and Retinol metabolism. The 23 most closely related DEGs were identified from the PPI network and combined with the GEPIA prognostic analysis, CDH3, ESRP1, FGF2, GBP2, KCNN4, KIT, SEMA4D, and ZEB1 were selected as hub genes, which are considered to be associated with poor prognosis of melanoma closely related. Besides, ten related drugs that may have therapeutic effects on melanoma were also screened. These newly discovered genes and drugs provide new ideas for further research on melanoma.

Seven Novel Genes Related to Cell Proliferation and Migration of VHL-Mutated Pheochromocytoma

Pheochromocytoma, as a neuroendocrine tumor with the highest genetic correlation in all types of tumors, has attracted extensive attention. Von Hipper Lindau (VHL) has the highest mutation frequency among the genes associated with pheochromocytoma. However, the effect of VHL on the proteome of pheochromocytoma remains to be explored. In this study, the VHL knockdown (VHL-KD) PC12 cell model was established by RNA interference (shRNA). We compared the proteomics of VHL-KD and VHL-WT PC12 cell lines. The results showed that the expression of 434 proteins (VHL shRNA/WT > 1.3) changed significantly in VHL-KD-PC12 cells. Among the 434 kinds of proteins, 83 were involved in cell proliferation, cell cycle and cell migration, and so on. More importantly, among these proteins, we found seven novel key genes, including Connective Tissue Growth Factor (CTGF), Syndecan Binding Protein (SDCBP), Cysteine Rich Protein 61 (CYR61/CCN1), Collagen Type III Alpha 1 Chain (COL3A1), Collagen Type I Alpha 1 Chain (COL1A1), Collagen Type V Alpha 2 Chain (COL5A2), and Serpin Family E Member 1 (SERPINE1), were overexpressed and simultaneously regulated cell proliferation and migration in VHL-KD PC12 cells. Furthermore, the abnormal accumulation of HIF2α caused by VHL-KD significantly increased the expression of these seven genes during hypoxia. Moreover, cell-counting, scratch, and transwell assays demonstrated that VHL-KD could promote cell proliferation and migration, and changed cell morphology. These findings indicated that inhibition of VHL expression could promote the development of pheochromocytoma by activating the expression of cell proliferation and migration associated genes.