scholarly journals Identification of Prognostic Genes in Neuroblastoma in Children by Weighted Gene Coexpression Network Analysis

2021 ◽  
Vol 2021 ◽  
pp. 1-13
Author(s):  
Jun Yang ◽  
Ying Zhang ◽  
Jiaying Zhou ◽  
Shaohua Wang

Background. Neuroblastoma is a malignant neuroendocrine tumor from the sympathetic nervous system, the most common extracranial tumor in children. Identifying potential prognostic markers of neuroblastoma can provide clues for early diagnosis, recurrence, and treatment. Methods. RNA sequence data and clinical features of 147 neuroblastomas were obtained from the TARGET (Therapeutically Applicable Research to Generate Effective Treatments project) database. Application weighted gene coexpression network analysis (WGCNA) was used to construct a free-scale gene coexpression network, to study the interrelationship between its potential modules and clinical features, and to identify hub genes in the module. We performed Lasso regression and Cox regression analyses to identify the three most important genes and develop a new prognostic model. Data from the GSE85047 cohort verified the predictive accuracy of the prognostic model. Results. 14 coexpression modules were constructed using WGCNA. Brown coexpression modules were found to be significantly associated with disease survival status. Multivariate Cox analysis was performed on genes from univariate Cox regression and Lasso regression analyses using the Cox proportional hazards regression model. Finally, we constructed a three-gene prognostic model: risk score = (0.003812659 ∗ CKB) + (−0.152376975 ∗ expDST) + (0.032032815 ∗ expDUT). The prognosis of samples in the high-risk group was significantly poorer than that of samples in the low-risk group ( P = 1.225 e − 06 ). The risk model was also regarded as an independent predictor of prognosis (HR = 1.632; 95% CI = 1.391–1.934; P < 0.001 ). Conclusion. Our study constructed a neuroblastoma coexpressing gene module and identified a prognostic potential risk model for prognosis in neuroblastoma.

2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Xiaoqiang Gu ◽  
Qiqi Zhang ◽  
Xueying Wu ◽  
Yue Fan ◽  
Jianxin Qian

Abstract Background Pancreatic adenocarcinoma (PAAD) is a nonimmunogenic tumor, and very little is known about the relationship between the host immune response and patient survival. We aimed to develop an immune prognostic model (IPM) and analyze its relevance to the tumor immune profiles of patients with PAAD. Methods We investigated differentially expressed genes between tumor and normal tissues in the TCGA PAAD cohort. Immune-related genes were screened from highly variably expressed genes with weighted gene correlation network analysis (WGCNA) to construct an IPM. Then, the influence of IPM on the PAAD immune profile was comprehensively analyzed. Results A total of 4902 genes highly variably expressed among primary tumors were used to construct a weighted gene coexpression network. One hundred seventy-five hub genes in the immune-related module were used for machine learning. Then, we established an IPM with four core genes (FCGR2B, IL10RA, and HLA-DRA) to evaluate the prognosis. The risk score predicted by IPM was an independent prognostic factor and had a high predictive value for the prognosis of patients with PAAD. Moreover, we found that the patients in the low-risk group had higher cytolytic activity and lower innate anti-PD-1 resistance (IPRES) signatures than patients in the high-risk group. Conclusions Unlike the traditional methods that use immune-related genes listed in public databases to screen prognostic genes, we constructed an IPM through WGCNA to predict the prognosis of PAAD patients. In addition, an IPM prediction of low risk indicated enhanced immune activity and a decreased anti-PD-1 therapeutic response.


2021 ◽  
Vol 11 ◽  
Author(s):  
Qingchuan Chen ◽  
Yuen Tan ◽  
Chao Zhang ◽  
Zhe Zhang ◽  
Siwei Pan ◽  
...  

BackgroundGastric cancer (GC) is a major public health problem worldwide. In recent decades, the treatment of gastric cancer has improved greatly, but basic research and clinical application of gastric cancer remain challenges due to the high heterogeneity. Here, we provide new insights for identifying prognostic models of GC.MethodsWe obtained the gene expression profiles of GSE62254 containing 300 samples for training. GSE15459 and TCGA-STAD for validation, which contain 200 and 375 samples, respectively. Weighted gene co-expression network analysis (WGCNA) was used to identify gene modules. We performed Lasso regression and Cox regression analyses to identify the most significant five genes to develop a novel prognostic model. And we selected two representative genes within the model for immunohistochemistry staining with 105 GC specimens from our hospital to verify the prediction efficiency. Moreover, we estimated the correlation coefficient between our model and immune infiltration using the CIBERSORT algorithm. The data from GSE15459 and TCGA cohort validated the robustness and predictive accuracy of this prognostic model.ResultsOf the 12 gene modules identified, 1,198 green-yellow module genes were selected for further analysis. Multivariate Cox analysis was performed on genes from univariate Cox regression and Lasso regression analysis using the Cox proportional hazards regression model. Finally, we constructed a five gene prognostic model: Risk Score = [(-0.7547) * Expression (ARHGAP32)] + [(-0.8272) * Expression (KLF5)] + [1.09 * Expression (MAMLD1)] + [0.5174 * Expression (MATN3)] + [1.66 * Expression (NES)]. The prognosis of samples in the high-risk group was significantly poorer than that of samples in the low-risk group (p = 6.503e-11). The risk model was also regarded as an independent predictor of prognosis (HR, 1.678, p &lt; 0.001). The observed correlation with immune cells suggested that this risk model could potentially predict immune infiltration.ConclusionThis study identified a potential risk model for prognosis and immune infiltration prediction in GC using WGCNA and Cox regression analysis.


2021 ◽  
Vol 20 ◽  
pp. 153303382110414
Author(s):  
Xiaoyong Li ◽  
Jiaqong Lin ◽  
Yuguo pan ◽  
Peng Cui ◽  
Jintang Xia

Background: Liver progenitor cells (LPCs) play significant roles in the development and progression of hepatocellular carcinoma (HCC). However, no studies on the value of LPC-related genes for evaluating HCC prognosis exist. We developed a gene signature of LPC-related genes for prognostication in HCC. Methods: To identify LPC-related genes, we analyzed mRNA expression arrays from a dataset (GSE57812 & GSE 37071) containing LPCs, mature hepatocytes, and embryonic stem cell samples. HCC RNA-Seq data from The Cancer Genome Atlas (TCGA) were used to explore the differentially expressed genes (DEGs) related to prognosis through DEG analysis and univariate Cox regression analysis. Lasso and multivariate Cox regression analyses were performed to construct the LPC-related gene prognostic model in the TCGA training dataset. This model was validated in the TCGA testing set and an external dataset (International Cancer Genome Consortium [ICGC] dataset). Finally, we investigated the relationship between this prognostic model with tumor-node-metastasis stage, tumor grade, and vascular invasion of HCC. Results: Overall, 1770 genes were identified as LPC-related genes, of which 92 genes were identified as DEGs in HCC tissues compared with normal tissues. Furthermore, we randomly assigned patients from the TCGA dataset to the training and testing cohorts. Twenty-six DEGs correlated with overall survival (OS) in the univariate Cox regression analysis. Lasso and multivariate Cox regression analyses were performed in the TCGA training set, and a 3-gene signature was constructed to stratify patients into 2 risk groups: high-risk and low-risk. Patients in the high-risk group had significantly lower OS than those in the low-risk group. Receiver operating characteristic curve analysis confirmed the signature's predictive capacity. Moreover, the risk score was confirmed to be an independent predictor for patients with HCC. Conclusion: We demonstrated that the LPC-related gene signature can be used for prognostication in HCC. Thus, targeting LPCs may serve as a therapeutic alternative for HCC.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e8128 ◽  
Author(s):  
Cheng Yue ◽  
Hongtao Ma ◽  
Yubai Zhou

Background Lung cancer has the highest morbidity and mortality worldwide, and lung adenocarcinoma (LADC) is the most common pathological subtype. Accumulating evidence suggests the tumor microenvironment (TME) is correlated with the tumor progress and the patient’s outcome. As the major components of TME, the tumor-infiltrated immune cells and stromal cells have attracted more and more attention. In this study, differentially expressed immune and stromal signature genes were used to construct a TME-related prognostic model for predicting the outcomes of LADC patients. Methods The expression profiles of LADC samples with clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). The differentially expressed genes (DEGs) related to the TME of LADC were identified using TCGA dataset by Wilcoxon rank sum test. The prognostic effects of TME-related DEGs were analyzed using univariate Cox regression. Then, the least absolute shrinkage and selection operator (LASSO) regression was performed to reduce the overfit and the number of genes for further analysis. Next, the prognostic model was constructed by step multivariate Cox regression and risk score of each sample was calculated. Then, survival and Receiver Operating Characteristic (ROC) analyses were conducted to validate the model using TCGA and GEO datasets, respectively. The Kyoto Encyclopedia of Genes and Genomes analysis of gene signature was performed using Gene Set Enrichment Analysis (GSEA). Finally, the overall immune status, tumor purity and the expression profiles of HLA genes of high- and low-risk samples was further analyzed to reveal the potential mechanisms of prognostic effects of the model. Results A total of 93 TME-related DEGs were identified, of which 23 DEGs were up-regulated and 70 DEGs were down-regulated. The univariate cox analysis indicated that 23 DEGs has the prognostic effects, the hazard ratio ranged from 0.65 to 1.25 (p < 0.05). Then, seven genes were screened out from the 23 DEGs by LASSO regression method and were further analyzed by step multivariate Cox regression. Finally, a three-gene (ADAM12, Bruton Tyrosine Kinase (BTK), ERG) signature was constructed, and ADAM12, BTK can be used as independent prognostic factors. The three-gene signature well stratified the LADC patients in both training (TCGA) and testing (GEO) datasets as high-risk and low-risk groups, the 3-year area under curve (AUC) of ROC curves of three GEO sets were 0.718 (GSE3141), 0.646 (GSE30219) and 0.643 (GSE50081). The GSEA analysis indicated that highly expressed ADAM12, BTK, ERG mainly correlated with the activation of pathways involving in focal adhesion, immune regulation. The immune analysis indicated that the low-risk group has more immune activities and higher expression of HLA genes than that of the high-risk group. In sum, we identified and constructed a three TME-related DEGs signature, which could be used to predict the prognosis of LADC patients.


2019 ◽  
Vol 49 (10) ◽  
pp. 1195-1206 ◽  
Author(s):  
Aiping Tian ◽  
Ke Pu ◽  
Boxuan Li ◽  
Min Li ◽  
Xiaoguang Liu ◽  
...  

2020 ◽  
Vol 16 ◽  
pp. 117693432092056
Author(s):  
Shuping Qu ◽  
Qiuyuan Shi ◽  
Jing Xu ◽  
Wanwan Yi ◽  
Hengwei Fan

This study was aimed at revealing the dynamic regulation of mRNAs, long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) in hepatocellular carcinoma (HCC) and to identify HCC biomarkers capable of predicting prognosis. Differentially expressed mRNAs (DEmRNAs), lncRNAs, and miRNAs were acquired by comparing expression profiles of HCC with normal samples, using an expression data set from The Cancer Genome Atlas. Altered biological functions and pathways in HCC were analyzed by subjecting DEmRNAs to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. Gene modules significantly associated with disease status were identified by weighted gene coexpression network analysis. An lncRNA-mRNA and an miRNA-mRNA coexpression network were constructed for genes in disease-related modules, followed by the identification of prognostic biomarkers using Kaplan-Meier survival analysis. Differential expression and association with the prognosis of 4 miRNAs were verified in independent data sets. A total of 1220 differentially expressed genes were identified between HCC and normal samples. Differentially expressed mRNAs were significantly enriched in functions and pathways related to “plasma membrane structure,” “sensory perception,” “metabolism,” and “cell proliferation.” Two disease-associated gene modules were identified. Among genes in lncRNA-mRNA and miRNA-mRNA coexpression networks, 9 DEmRNAs and 7 DEmiRNAs were identified to be potential prognostic biomarkers. MIMAT0000102, MIMAT0003882, and MIMAT0004677 were successfully validated in independent data sets. Our results may advance our understanding of molecular mechanisms underlying HCC. The biomarkers may contribute to diagnosis in future clinical practice.


2020 ◽  
Vol 20 (1) ◽  
Author(s):  
Baiyang Yu ◽  
Jianbin Liu ◽  
Di Wu ◽  
Ying Liu ◽  
Weijian Cen ◽  
...  

An amendment to this paper has been published and can be accessed via the original article.


2019 ◽  
Vol 2019 ◽  
pp. 1-11 ◽  
Author(s):  
Yan Li ◽  
Xiao_nan He ◽  
Chao Li ◽  
Ling Gong ◽  
Min Liu

Background. Identification of potential molecular targets of acute myocardial infarction is crucial to our comprehensive understanding of the disease mechanism. However, studies of gene coexpression analysis via jointing multiple microarray data of acute myocardial infarction still remain restricted. Methods. Microarray data of acute myocardial infarction (GSE48060, GSE66360, GSE97320, and GSE19339) were downloaded from Gene Expression Omnibus database. Three data sets without heterogeneity (GSE48060, GSE66360, and GSE97320) were subjected to differential expression analysis using MetaDE package. Differentially expressed genes having upper 25% variation across samples were imported in weighted gene coexpression network analysis. Functional and pathway enrichment analyses were conducted for genes in the most significant module using DAVID. The predicted microRNAs to regulate target genes in the most significant module were identified using TargetScan. Moreover, subpathway analyses using iSubpathwayMiner package and GenCLiP 2.0 were performed on hub genes with high connective weight in the most significant module. Results. A total of 1027 differentially expressed genes and 33 specific modules were screened out between acute myocardial infarction patients and control samples. Ficolin (collagen/fibrinogen domain containing) 1 (FCN1), CD14 molecule (CD14), S100 calcium binding protein A9 (S100A9), and mitochondrial aldehyde dehydrogenase 2 (ALDH2) were identified as critical target molecules; hsa-let-7d, hsa-let-7b, hsa-miR-124-3, and hsa-miR-9-1 were identified as potential regulators of the expression of the key genes in the two biggest modules. Conclusions. FCN1, CD14, S100A9, ALDH2, hsa-let-7d, hsa-let-7b, hsa-miR-124-3, and hsa-miR-9-1 were identified as potential candidate regulators in acute myocardial infarction. These findings might provide new comprehension into the underlying molecular mechanism of disease.


Sign in / Sign up

Export Citation Format

Share Document