Screening and Identification of Four Prognostic Genes Related to Immune Infiltration and G-Protein Coupled Receptors Pathway in Lung Adenocarcinoma

Methylation and transcriptome analyses construct a prognostic model and reveal the suppressor role of VMO1 in lung adenocarcinoma

BACKGROUND

DNA methylation is an important epigenetic mechanism of gene regulation. The aberrant DNA methylation has been found to play an important role in the initiation and progression of tumors.

RESULTS

Transcriptome and DNA methylation data of lung adenocarcinoma (LUAD) patients were co-analyzed and 95 methylation-driven genes (MDGs) was found in relation to LUAD. A prognostic model based on 3 MDGs (GMNN, SPINK2 and VMO1) was constructed by Univariate and Multivariate cox regression analyses. The risk score generated from the prognostic model could be used to classify LUAD patients into high and low risk groups. Furthermore, it was found that the risk score was associated with tumor microenvironment (TME) and clinical characteristics (survival status and T stage) of patients. Interestingly, we identified and validated that the patients in the low-risk group responded better to immunotherapy treatment. Then, a nomogram model based on the risk score and clinical characteristics was established which showed significant prediction value. The down-regulation and hypermethylation levels of vitelline membrane outer layer protein 1 homolog (VMO1) were verified in paired LUAD tumor and non-tumor tissues by pyrosequencing assay and RT-qPCR. Furthermore, MTT, migration and wound healing assays were performed with lentivector-mediated ectopic over-expression and 5-Aza-dC demethylation followed by siRNA rescue experiments to investigate the role of VMO1 in LUAD cells. Our results indicated that VMO1 could inhibit proliferation and migration of A549 and NCI-H1299 cells.

CONCLUSIONS

In summary, our experiments constructed a prognostic model with high capacity for risk prediction in LUAD patients. VMO1 had a malignant suppressor role in LUAD cells. The correlation between risk score and TME might elucidate a potential mechanism of oncogenesis and provide an avenue for further therapeutic targets.

Multi-omics identification of GPCR gene features in lung adenocarcinoma based on multiple machine learning combinations

Background: Lung adenocarcinoma is a common malignant tumor that ranks second in the world and has a high mortality rate. G protein-coupled receptors (GPCRs) have been reported to play an important role in cancer; however, G protein-coupled receptor-associated features have not been adequately investigated. Methods: In this study, GPCR-related genes were screened at single-cell and bulk transcriptome levels based on AUcell, single-sample gene set enrichment analysis (ssGSEA) and weighted gene co-expression network (WGCNA) analysis. And a new machine learning framework containing 10 machine learning algorithms and their multiple combinations was used to construct a consensus G protein-coupled receptor-related signature (GPCRRS). GPCRRS was validated in the training set and external validation set. We constructed GPCRRS-integrated nomogram clinical prognosis prediction tools. Multi-omics analyses included genomics, single-cell transcriptomics, and bulk transcriptomics to gain a more comprehensive understanding of prognostic features. We assessed the response of risk subgroups to immunotherapy and screened for personalized drugs targeting specific risk subgroups. Finally, the expression of key GPCRRS genes was verified by RT-qPCR. Results: In this study, we identified 10 GPCR-associated genes that were significantly associated with the prognosis of lung adenocarcinoma by single-cell transcriptome and bulk transcriptome. Univariate and multivariate showed that the survival rate was higher in low risk than in high risk, which also suggested that the model was an independent prognostic factor for LUAD. In addition, we observed significant differences in biological function, mutational landscape, and immune cell infiltration in the tumor microenvironment between high and low risk groups. Notably, immunotherapy was also relevant in the high and low risk groups. In addition, potential drugs targeting specific risk subgroups were identified. Conclusion: In this study, we constructed and validated a lung adenocarcinoma G protein-coupled receptor-related signature, which has an important role in predicting the prognosis of lung adenocarcinoma and the effect of immunotherapy. It is hypothesized that LDHA, GPX3 and DOCK4 are new potential targets for lung adenocarcinoma, which can achieve breakthroughs in prognosis prediction, targeted prevention and treatment of lung adenocarcinoma and provide important guidance for anti-tumor.

A more novel and powerful prognostic gene signature of lung adenocarcinoma determined from the immune cell infiltration landscape

Background

Lung adenocarcinoma (LUAD) is the leading histological subtype of lung cancer worldwide, causing high mortality each year. The tumor immune cell infiltration (ICI) is closely associated with clinical outcome with LUAD patients. The present study was designed to construct a gene signature based on the ICI of LUAD to predict prognosis.

Methods

Downloaded the raw data of three cohorts of the TCGA-LUAD, GSE72094, and GSE68465 and treat them as training cohort, validation cohort one, and validation cohort two for this research. Unsupervised clustering detailed grouped LUAD cases of the training cohort based on the ICI profile. The univariate Cox regression and Kaplan-Meier was adopted to identify potential prognostic genes from the differentially expressed genes recognized from the ICI clusters. A risk score-based prognostic signature was subsequently developed using LASSO-penalized Cox regression analysis. The Kaplan-Meier analysis, Cox analysis, ROC, IAUC, and IBS were constructed to assess the ability to predict the prognosis and effects of clinical variables in another two independent validation cohorts. More innovatively, we searched similar papers in the most recent year and made comprehensive comparisons with ours. GSEA was used to discover the related signaling pathway. The immune relevant signature correlation identification and immune infiltrating analysis were used to evaluate the potential role of the signature for immunotherapy and recognize the critical immune cell that can influence the signature's prognosis capability.

Results

A signature composed of thirteen gene including ABCC2, CCR2, CERS4, CMAHP, DENND1C, ECT2, FKBP4, GJB3, GNG7, KRT6A, PCDH7, PLK1, and VEGFC, was identified as significantly associated with the prognosis in LUAD patients. The thirteen-gene signature exhibited independence in evaluating the prognosis of LUAD patients in our training and validation cohorts. Compared to our predecessors, our model has an advantage in predictive power. Nine well know immunotherapy targets, including TBX2, TNF, CTLA4, HAVCR2, GZMB, CD8A, PRF1, GZMA, and PDCD1 were recognized correlating with our signature. The mast cells were found to play vital parts in backing on the thirteen-gene signature's outcome predictive capacity.

Conclusions

Collectively, the current study indicated a robust thirteen-gene signature that can accurately predict LUAD prognosis, which is superior to our predecessors in predictive ability. The immune relevant signatures, TBX2, TNF, CTLA4, HAVCR2, GZMB, CD8A, PRF1, GZMA, PDCD1, and mast cells infiltrating were found closely correlate with the thirteen-gene signature's power.

Integrative analysis and experiments to explore angiogenesis regulators correlated with poor prognosis, immune infiltration and cancer progression in lung adenocarcinoma

Angiogenesis is the process of capillary sprouting from pre-existing vessels and it plays a critical role in the carcinogenic process of lung adenocarcinoma (LUAD). However, the association of angiogenesis regulators with the prognosis and progression of LUAD needs to be further elucidated. In this study, we adopted differential expression analysis, Cox proportional hazards (PH) regression analysis and experimental validation to identify angiogenesis regulators correlated with a poor prognosis, immune infiltration and cancer progression in LUAD. These results showed that the diagnostic and prognostic models based on COL5A2 and EPHB2 served as independent biomarkers with superior predictive ability. The patients in the high-risk group exhibited a worse prognosis in the TCGA cohort (P < 0.001, HR = 1.72, 95% CI 1.28-2.30), GSE310210 cohort (P = 0.005, HR = 2.87, 95% CI 1.46-5.61), and GSE31019 cohort (P = 0.01, HR = 2.14, 95% CI 1.19-3.86) than patients in the low-risk group. The high prognostic risk patients had a higher TMB (P < 0.001); higher fractions of M0 macrophages, neutrophils, NK cells resting, and T cells CD4 memory activated (P < 0.05); and higher expression of immune checkpoints PD-1, PDL-1, PDL-2, and B7H3 (P < 0.001). Patients in the high-risk group were more sensitive to chemotherapeutic drugs and molecular targeted drugs such as cisplatin, doxorubicin, gefitinib, and bosutinib (P < 0.0001). In addition, inhibition of COL5A2 and EPHB2 effectively suppressed the proliferation and migration of LUAD cells. The current study identified angiogenesis regulators as potential biomarkers and therapeutic targets for LUAD and may help to further optimize cancer therapy.

An Aging-Related Gene Signature-Based Model for Risk Stratification and Prognosis Prediction in Lung Adenocarcinoma

Aging is an inevitable time-dependent process associated with a gradual decline in many physiological functions. Importantly, some studies have supported that aging may be involved in the development of lung adenocarcinoma (LUAD). However, no studies have described an aging-related gene (ARG)-based prognosis signature for LUAD. Accordingly, in this study, we analyzed ARG expression data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). After LASSO and Cox regression analyses, a six ARG-based signature (APOC3, EPOR, H2AFX, MXD1, PLCG2, and YWHAZ) was constructed using TCGA dataset that significantly stratified cases into high- and low-risk groups in terms of overall survival (OS). Cox regression analysis indicated that the ARG signature was an independent prognostic factor in LUAD. A nomogram based on the ARG signature and clinicopathological factors was developed in TCGA cohort and validated in the GEO dataset. Moreover, to visualize the prediction results, we established a web-based calculator yurong.shinyapps.io/ARGs_LUAD/. Calibration plots showed good consistency between the prediction of the nomogram and actual observations. Receiver operating characteristic curve and decision curve analyses indicated that the ARG nomogram had better OS prediction and clinical net benefit than the staging system. Taken together, these results established a genetic signature for LUAD based on ARGs, which may promote individualized treatment and provide promising novel molecular markers for immunotherapy.