scholarly journals The Prognostic Value of Bone Morphogenetic Proteins and Their Receptors in Lung Adenocarcinoma

2021 ◽  
Vol 11 ◽  
Author(s):  
Wangyang Meng ◽  
Han Xiao ◽  
Rong Zhao ◽  
Dong Li ◽  
Kuo Li ◽  
...  

BackgroundBone morphogenetic proteins (BMPs) regulate tumor progression via binding to their receptors (BMPRs). However, the expression and clinical significance of BMPs/BMPRs in lung adenocarcinoma remain unclear due to a lack of systematic studies.MethodsThis study screened differentially expressed BMPs/BMPRs (deBMPs/BMPRs) in a training dataset combining TCGA-LUAD and GTEx-LUNG and verified them in four GEO datasets. Their prognostic value was evaluated via univariate and multivariate Cox regression analyses. LASSO was performed to construct an initial risk model. Subsequently, after weighted gene co-expression network analysis (WGCNA), differential expression analysis, and univariate Cox regression analysis, hub genes co-expressed with differentially expressed BMPs/BMPRs were filtered out to improve the risk model and explore potential mechanisms. The improved risk model was re-established via LASSO combining hub genes with differentially expressed BMPs/BMPRs as the core. In the testing cohort including 93 lung adenocarcinoma patients, immunohistochemistry (IHC) was performed to verify BMP5 protein expression and its association with prognosis.ResultsBMP2, BMP5, BMP6, GDF10, and ACVRL1 were verified as downregulated in lung adenocarcinoma. Survival analysis identified BMP5 as an independent protective prognostic factor. We also found that BMP5 was significantly correlated with EGFR expression and mutations, suggesting that BMP5 may play a role in targeted therapy. The initial risk model containing only BMP5 showed a significant correlation (HR: 1.71, 95% CI: 1.28−2.28, p: 3e-04) but low prognostic accuracy (AUC of 1-year survival: 0.6, 3-year survival: 0.6, 5-year survival: 0.63). Seventy-nine hub genes co-expressed with BMP5 were identified, and their functions were enriched in cell migration and tumor metastasis. The re-established risk model showed greater prognostic correlation (HR: 2.58, 95% CI: 1.92–3.46, p: 0) and value (AUC of 1-year survival: 0.72, 3-year survival: 0.69, and 5-year survival: 0.68). IHC results revealed that BMP5 protein was also downregulated in lung adenocarcinoma and higher expression was markedly associated with better prognosis (HR: 0.44, 95% CI: 0.23–0.85, p: 0.0145).ConclusionBMP5 is a potential crucial target for lung adenocarcinoma treatment based on significant differential expression and superior prognostic value.

2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Zhihao Wang ◽  
Kidane Siele Embaye ◽  
Qing Yang ◽  
Lingzhi Qin ◽  
Chao Zhang ◽  
...  

Abstract Background Given that dysregulated metabolism has been recently identified as a hallmark of cancer biology, this study aims to establish and validate a prognostic signature of lung adenocarcinoma (LUAD) based on metabolism-related genes (MRGs). Methods The gene sequencing data of LUAD samples with clinical information and the metabolism-related gene set were obtained from The Cancer Genome Atlas (TCGA) and Molecular Signatures Database (MSigDB), respectively. The differentially expressed MRGs were identified by Wilcoxon rank sum test. Then, univariate cox regression analysis was performed to identify MRGs that related to overall survival (OS). A prognostic signature was developed by multivariate Cox regression analysis. Furthermore, the signature was validated in the GSE31210 dataset. In addition, a nomogram that combined the prognostic signature was created for predicting the 1-, 3- and 5-year OS of LUAD. The accuracy of the nomogram prediction was evaluated using a calibration plot. Finally, cox regression analysis was applied to identify the prognostic value and clinical relationship of the signature in LUAD. Results A total of 116 differentially expressed MRGs were detected in the TCGA dataset. We found that 12 MRGs were most significantly associated with OS by using the univariate regression analysis in LUAD. Then, multivariate Cox regression analyses were applied to construct the prognostic signature, which consisted of six MRGs-aldolase A (ALDOA), catalase (CAT), ectonucleoside triphosphate diphosphohydrolase-2 (ENTPD2), glucosamine-phosphate N-acetyltransferase 1 (GNPNAT1), lactate dehydrogenase A (LDHA), and thymidylate synthetase (TYMS). The prognostic value of this signature was further successfully validated in the GSE31210 dataset. Furthermore, the calibration curve of the prognostic nomogram demonstrated good agreement between the predicted and observed survival rates for each of OS. Further analysis indicated that this signature could be an independent prognostic indicator after adjusting to other clinical factors. The high-risk group patients have higher levels of immune checkpoint molecules and are therefore more sensitive to immunotherapy. Finally, we confirmed six MRGs protein and mRNA expression in six lung cancer cell lines and firstly found that ENTPD2 might played an important role on LUAD cells colon formation and migration. Conclusions We established a prognostic signature based on MRGs for LUAD and validated the performance of the model, which may provide a promising tool for the diagnosis, individualized immuno-/chemotherapeutic strategies and prognosis in patients with LUAD.


Author(s):  
Haihui Zhong ◽  
Jie Wang ◽  
Yaru Zhu ◽  
Yefeng Shen

Lung adenocarcinoma (LUAD) is the most common malignancy, leading to more than 1 million related deaths each year. Due to low long-term survival rates, the exploration of molecular mechanisms underlying LUAD progression and novel prognostic predictors is urgently needed to improve LUAD treatment. In our study, cancer-specific differentially expressed genes (DEGs) were identified using the robust rank aggregation (RRA) method between tumor and normal tissues from six Gene Expression Omnibus databases (GSE43458, GSE62949, GSE68465, GSE115002, GSE116959, and GSE118370), followed by a selection of prognostic modules using weighted gene co-expression network analysis. Univariate Cox regression, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses were applied to identify nine hub genes (CBFA2T3, CR2, SEL1L3, TM6SF1, TSPAN32, ITGA6, MAPK11, RASA3, and TLR6) that constructed a prognostic risk model. The RNA expressions of nine hub genes were validated in tumor and normal tissues by RNA-sequencing and single-cell RNA-sequencing, while immunohistochemistry staining from the Human Protein Atlas database showed consistent results in the protein levels. The risk model revealed that high-risk patients were associated with poor prognoses, including advanced stages and low survival rates. Furthermore, a multivariate Cox regression analysis suggested that the prognostic risk model could be an independent prognostic factor for LUAD patients. A nomogram that incorporated the signature and clinical features was additionally built for prognostic prediction. Moreover, the levels of hub genes were related to immune cell infiltration in LUAD microenvironments. A CMap analysis identified 13 small molecule drugs as potential agents based on the risk model for LUAD treatment. Thus, we identified a prognostic risk model including CBFA2T3, CR2, SEL1L3, TM6SF1, TSPAN32, ITGA6, MAPK11, RASA3, and TLR6 as novel biomarkers and validated their prognostic and predicted values for LUAD.


2020 ◽  
Author(s):  
Zhihao Wang ◽  
Kidane Siele Embaye ◽  
Qing Yang ◽  
Lingzhi Qin ◽  
Chao Zhang ◽  
...  

Abstract Background: Given that dysregulated metabolism has been recently identified as a hallmark of cancer biology, this study aims to establish and validate a prognostic signature of lung adenocarcinoma (LUAD) based on metabolism-related genes(MRGs). Methods: The gene sequencing data of LUAD samples with clinical information and the metabolism-related gene set were obtained from The Cancer Genome Atlas (TCGA) and Molecular Signatures Database (MSigDB), respectively. The differentially expressed MRGs were identified by Wilcoxon rank sum test. Then, univariate Cox regression analysis were performed to identify MRGs that related to overall survival(OS). A prognostic signature was developed by multivariate Cox regression analysis. Furthermore, the signature was validated in the GSE31210 dataset. In addition, a nomogram that combined the prognostic signature was created for predicting the 1-, 3- and 5-year OS of LUAD.The accuracy of the nomogram prediction was evaluated using a calibration plot. Finally, cox regression analysis was applied to identify the prognostic value and clinical relationship of the signature in LUAD. Results: A total of 116 differentially expressed MRGs were detected in the TCGA dataset. We found that 12 MRGs were most significantly associated with OS by using the univariate regression analysis in LUAD. Then, multivariate Cox regression analyses were applied to construct the prognostic signature, which consisted of six MRGs(ALDOA, CAT, ENTPD2, GNPNAT1, LDHA, and TYMS). The prognostic value of this signature was further successfully validated in the GSE31210 dataset. Furthermore, the calibration curve of the prognostic nomogram demonstrated good agreement between the predicted and observed survival rates for each of OS. Further analysis indicated that this signature could be an independent prognostic indicator after adjusting to other clinical factors. Finally, the signature was found to be associated with various clinicopathological features. Conclusions: We established a prognostic signature based on MRGs for LUAD and validated the performance of the model, which may provide a promising tool for the diagnosis and prognosis in patients with LUAD.


2020 ◽  
Author(s):  
Zhihao Wang ◽  
Kidane Siele Embaye ◽  
Qing Yang ◽  
Lingzhi Qin ◽  
Chao Zhang ◽  
...  

Abstract Background: Given that dysregulated metabolism has been recently identified as a hallmark of cancer biology, this study aims to establish and validate a prognostic signature of lung adenocarcinoma (LUAD) based on metabolism-related genes(MRGs).Methods: The gene sequencing data of LUAD samples with clinical information and the metabolism-related gene set were obtained from The Cancer Genome Atlas (TCGA) and Molecular Signatures Database (MSigDB), respectively. The differentially expressed MRGs were identified by Wilcoxon rank sum test. Then, univariate Cox regression analysis were performed to identify MRGs that related to overall survival(OS). A prognostic signature was developed by multivariate Cox regression analysis. Furthermore, the signature was validated in the GSE31210 dataset. In addition, a nomogram that combined the prognostic signature was created for predicting the 1-, 3- and 5-year OS of LUAD.The accuracy of the nomogram prediction was evaluated using a calibration plot. Finally, cox regression analysis was applied to identify the prognostic value and clinical relationship of the signature in LUAD.Results: A total of 116 differentially expressed MRGs were detected in the TCGA dataset. We found that 12 MRGs were most significantly associated with OS by using the univariate regression analysis in LUAD. Then, multivariate Cox regression analyses were applied to construct the prognostic signature, which consisted of six MRGs (ALDOA, CAT, ENTPD2, GNPNAT1, LDHA, and TYMS). The prognostic value of this signature was further successfully validated in the GSE31210 dataset. Furthermore, the calibration curve of the prognostic nomogram demonstrated good agreement between the predicted and observed survival rates for each of OS. Further analysis indicated that this signature could be an independent prognostic indicator after adjusting to other clinical factors. The signature was found to be associated with various clinicopathological features. Finally, we confirmed six MRGs protein and mRNA expression in six lung cancer cells and firstly found that ENTPD2 might played an important role on lung adenocarcinoma cell lines clone formation and migration.Conclusions: We established a prognostic signature based on MRGs for LUAD and validated the performance of the model, which may provide a promising tool for the diagnosis and prognosis in patients with LUAD.


Author(s):  
Yongmei Wang ◽  
Guimin Zhang ◽  
Ruixian Wang

Background: This study aims to explore the prognostic values of CT83 and CT83-related genes in lung adenocarcinoma (LUAD). Methods: We downloaded the mRNA profiles of 513 LUAD patients (RNA sequencing data) and 246 NSCLC patients (Affymetrix Human Genome U133 Plus 2.0 Array) from TCGA and GEO databases. According to the median expression of CT83, the TCGA samples were divided into high and low expression groups, and differential expression analysis between them was performed. Functional enrichment analysis of differential expression genes (DEGs) was conducted. Univariate Cox regression analysis and LASSO Cox regression analysis were performed to screen the optimal prognostic DEGs. Then we established the prognostic model. A Nomogram model was constructed to predict the overall survival (OS) probability of LUAD patients. Results: CT83 expression was significantly correlated to the prognosis of LUAD patients. A total of 59 DEGs were identified, and a predictive model was constructed based on six optimal CT83-related DEGs, including CPS1, RHOV, TNNT1, FAM83A, IGF2BP1, and GRIN2A, could effectively predict the prognosis of LUAD patients. The nomogram could reliably predict the OS of LUAD patients. Moreover, the six important immune checkpoints (CTLA4, PD1, IDO1, TDO2, LAG3, and TIGIT) were closely correlated with the Risk Score, which was also differentially expressed between the LUAD samples with high and low-Risk Scores, suggesting that the poor prognosis of LUAD patients with high-Risk Score might be due to the immunosuppressive microenvironments. Conclusion: A prognostic model based on six optimal CT83 related genes could effectively predict the prognosis of LUAD patients.


2021 ◽  
Author(s):  
Rui Feng ◽  
Jian Li ◽  
Weiling Xuan ◽  
Hanbo Liu ◽  
Dexin Cheng ◽  
...  

Abstract Background Hepatocellular carcinoma (HCC) is a prevalent primary liver cancer and the main cause of cancer mortality. Its high complexity and dismal prognosis bring dramatic difficulty to treatment. Due to the disclosed dual functions of autophagy in cancer development, understanding autophagy-related genes devotes into seeking novel biomarkers for HCC. Methods Differential expression of genes in normal and tumor groups was analyzed to acquire autophagy-related genes in HCC. GO and KEGG pathway analyses were conducted on these genes. Genes were then screened by univariate regression analysis. The screened genes were subjected to multivariate Cox regression analysis to build a prognostic model. The model was validated by ICGC validation set. Results Altogether, 42 autophagy-related differential genes were screened by differential expression analysis. Enrichment analysis showed that they were mainly enriched in pathways including regulation of autophagy and cell apoptosis. Genes were screened by univariate analysis and multivariate Cox regression analysis to build a prognostic model. The model was constituted by 6 feature genes: EIF2S1, BIRC5, SQSTM1, ATG7, HDAC1, FKBP1A. Validation confirmed the accuracy and independence of this model in predicting HCC patient’s prognosis. Conclusion A total of 6 feature genes were identified to build a prognostic risk model. This model is conducive to investigating interplay between autophagy-related genes and HCC prognosis.


2021 ◽  
Vol 1 (3) ◽  
pp. 77-87
Author(s):  
Gong Xiao ◽  
Qiongjing Yuan ◽  
Wei Wang

Background: Multiple myeloma (MM) is one of the most common cancers of the blood system. N6-methyladenosine (m6A) plays an important role in cancer progression. We aimed to investigate the prognostic relevance of the m6A score in multiple myeloma through a series of bioinformatics analyses. Methods: The microarray dataset GSE4581 and GSE57317 used in this study were downloaded from the Gene Expression Omnibus (GEO) database. The m6A score was calculated using the GSVA package. The Random forests, univariate Cox regression analysis and Lasso analyses were performed for the differentially expressed genes (DEGs). Kaplan–Meier analysis and an ROC curve were used to diagnose the effectiveness of the model. Results: The GSVA R software package was used to predict the function. A total of 21 m6A genes were obtained, and 286 DEGs were identified between high and low m6A score groups. The risk model was constructed and composed of PRX, LBR, RB1, FBXL19-AS1, ARSK, MFAP3L, SLC44A3, UNC119 and SHCBP1. Functional analysis of risk score showed that with the increase in the risk score, Activated CD4 T cells, Memory B cells and Type 2 T helper cells were highly infiltrated. Conclusions: Immune checkpoints such as HMGB1, TGFB1, CXCL9 and HAVCR2 were significantly positively correlated with the risk score. We believe that the m6A score has a certain prognostic value in multiple myeloma.


2020 ◽  
Author(s):  
Xiang Zhou ◽  
Keying Zhang ◽  
Fa Yang ◽  
Chao Xu ◽  
Jianhua Jiao ◽  
...  

Abstract Background: Hepatocellular carcinoma (HCC) is a disease with higher morbidity, mortality, and poor prognosis in the whole world. Understanding the crosslink between HCC and the immune system is essential for people to uncover a few potential and valuable therapeutic strategies. This study aimed to reveal the correlation between HCC and immune-related genes and establish a clinical evaluation model. Methods: We had analyzed the clinical information consisted of 373 HCC and 49 normal samples from the cancer genome atlas (TCGA). The differentially expressed genes (DEGs) were selected by the Wilcoxon test and the immune-related differentially expressed genes (IRDEGs) in DEGs were identified by matching DEGs with immune-related genes downloaded from the ImmPort database. Furthermore, the univariate Cox regression analysis and multivariate Cox regression analysis were performed to construct a prognostic risk model. Then, twenty-two types of tumor immune-infiltrating cells (TIICs) were downloaded from Tumor Immune Estimation Resource (TIMER) and were used to construct the correlational graphs between the TIICs and risk score by the CIBERSORT. Subsequently, the transcription factors (TFs) were gained in the Cistrome website and the differentially expressed TFs (DETFs) were achieved. Finally, the KEGG pathway analysis and GO analysis were performed to further understand the molecular mechanisms between DETFs and PDIRGs.Results: In our study, 5839 DEGs, 326 IRDEGs, and 31 prognosis-related IRDEGs (PIRDEGs) were identified. And 8 optimal PIRDEGs were employed to construct a prognostic risk model by multivariate Cox regression analysis. The correlation between risk genes and clinical characterizations and TIICs has verified that the prognostic model was effective in predicting the prognosis of HCC patients. Finally, several important immune-related pathways and molecular functions of the eight PIRDEGs were significantly enriched and there was a distinct association between the risk IRDEGs and TFs. Conclusion: The prognostic risk model showed a more valuable predicting role for HCC patients, and produced many novel therapeutic targets and strategies for HCC.


2021 ◽  
Vol 11 ◽  
Author(s):  
Fen Liu ◽  
Zongcheng Yang ◽  
Lixin Zheng ◽  
Wei Shao ◽  
Xiujie Cui ◽  
...  

BackgroundGastric cancer is a common gastrointestinal malignancy. Since it is often diagnosed in the advanced stage, its mortality rate is high. Traditional therapies (such as continuous chemotherapy) are not satisfactory for advanced gastric cancer, but immunotherapy has shown great therapeutic potential. Gastric cancer has high molecular and phenotypic heterogeneity. New strategies for accurate prognostic evaluation and patient selection for immunotherapy are urgently needed.MethodsWeighted gene coexpression network analysis (WGCNA) was used to identify hub genes related to gastric cancer progression. Based on the hub genes, the samples were divided into two subtypes by consensus clustering analysis. After obtaining the differentially expressed genes between the subtypes, a gastric cancer risk model was constructed through univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox regression analysis. The differences in prognosis, clinical features, tumor microenvironment (TME) components and immune characteristics were compared between subtypes and risk groups, and the connectivity map (CMap) database was applied to identify potential treatments for high-risk patients.ResultsWGCNA and screening revealed nine hub genes closely related to gastric cancer progression. Unsupervised clustering according to hub gene expression grouped gastric cancer patients into two subtypes related to disease progression, and these patients showed significant differences in prognoses, TME immune and stromal scores, and suppressive immune checkpoint expression. Based on the different expression patterns between the subtypes, we constructed a gastric cancer risk model and divided patients into a high-risk group and a low-risk group based on the risk score. High-risk patients had a poorer prognosis, higher TME immune/stromal scores, higher inhibitory immune checkpoint expression, and more immune characteristics suitable for immunotherapy. Multivariate Cox regression analysis including the age, stage and risk score indicated that the risk score can be used as an independent prognostic factor for gastric cancer. On the basis of the risk score, we constructed a nomogram that relatively accurately predicts gastric cancer patient prognoses and screened potential drugs for high-risk patients.ConclusionsOur results suggest that the 7-gene signature related to tumor progression could predict the clinical prognosis and tumor immune characteristics of gastric cancer.


2020 ◽  
Vol 11 ◽  
Author(s):  
Xin Qiu ◽  
Qin-Han Hou ◽  
Qiu-Yue Shi ◽  
Hai-Xing Jiang ◽  
Shan-Yu Qin

BackgroundIntratumoral oxidative stress (OS) has been associated with the progression of various tumors. However, OS has not been considered a candidate therapeutic target for pancreatic cancer (PC) owing to the lack of validated biomarkers.MethodsWe compared gene expression profiles of PC samples and the transcriptome data of normal pancreas tissues from The Cancer Genome Atlas (TCGA) and Genome Tissue Expression (GTEx) databases to identify differentially expressed OS genes in PC. PC patients’ gene profile from the Gene Expression Omnibus (GEO) database was used as a validation cohort.ResultsA total of 148 differentially expressed OS-related genes in PC were used to construct a protein-protein interaction network. Univariate Cox regression analysis, least absolute shrinkage, selection operator analysis revealed seven hub prognosis-associated OS genes that served to construct a prognostic risk model. Based on integrated bioinformatics analyses, our prognostic model, whose diagnostic accuracy was validated in both cohorts, reliably predicted the overall survival of patients with PC and cancer progression. Further analysis revealed significant associations between seven hub gene expression levels and patient outcomes, which were validated at the protein level using the Human Protein Atlas database. A nomogram based on the expression of these seven hub genes exhibited prognostic value in PC.ConclusionOur study provides novel insights into PC pathogenesis and provides new genetic markers for prognosis prediction and clinical treatment personalization for PC patients.


Sign in / Sign up

Export Citation Format

Share Document