Exploring and comparing of the gene expression and methylation differences between lung adenocarcinoma and squamous cell carcinoma

Lung cancer is one of the most frequently diagnosed malignant tumors and the main reason for cancer-related death around the world, whereas nonsmall cell lung cancer that consists two subtypes: lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC) is responsible for an estimated 85% of all lung cancers. The current study aimed to explore gene expression and methylation differences between LUAD and LUSC. EdgeR was used to identify differentially regulated genes between normal and cancer in the LUAD and LUSC extracted from The Cancer Genome Atlas (TCGA), respectively, whereas SAM was used to find genes with differential methylation between normal and cancer in the LUAD and LUSC, respectively. Finally, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to analyze the function which these genes enriched in. A total of 391 genes with opposite methylation patterns in LUAD and LUSC and four functional pathways were obtained (false discovery rate (FDR) < 0.1). These pathways mainly included fat digestion and absorption, phenylalanine metabolism, bile secretion, and so on, which were related to the airframe nutrition metabolic pathway. Moreover, two genes CTSE (cathepsin E) and solute carrier family 5 member 7 (SLC5A7) were also found, among which CTSE was overexpressed and hypomethylated in LUAD corresponding to normal lung tissues, whereas SLC5A7 showed the opposite in LUAD. In conclusion, this study investigated the differences between the gene expression and methylation patterns in LUAD and LUSC, and explored their different biological characteristics. Further understanding of these differences may promote the discovery and development of new, accurate strategies for the prevention, diagnosis, and treatment of lung cancer.     

The increased expression and aberrant methylation of SHC1 in non–small cell lung cancer: Integrative analysis of clinical and bioinformatics databases

Despite the previous evidence showing that SHC adaptor protein 1 (SHC1) could encode three distinct isoforms (p46SHC, p52SHC and p66SHC) that function in different activities such as regulating life span and Ras activation, the precise underlying role of SHC1 in lung cancer also remains obscure. In this study, we firstly found that SHC1 expression was up-regulated both in lung adenocarcinoma (LUAD) and in lung squamous cell carcinoma (LUSC) tissues. Furthermore, compared to patients with lower SHC1 expression, LUAD patients with higher expression of SHC1 had poorer overall survival (OS). Moreover, higher expression of SHC1 was also associated with worse OS in patients with stages 1 and 2 but not stage 3 lung cancer. Significantly, the analysis showed that SHC1 methylation level was associated with OS in lung cancer patients. It seemed that the methylation level at specific probes within SHC1 showed negative correlations with SHC1 expression both in LUAD and in LUSC tissues. The LUAD and LUSC patients with hypermethylated SHC1 at cg12473916 and cg19356022 probes had a longer OS. Therefore, it is reasonable to conclude that SHC1 has a potential clinical significance in LUAD and LUSC patients.     

In silico assessment of EpCAM transcriptional expression and determination of the prognostic biomarker for human lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC)

Epithelial cell adhesion molecule (EpCAM) is a transmembrane glycoprotein which is involved in cell signaling, proliferation, maturation, and movement, all of which are crucial for the proper development of cells and tissues. Cleavage of the EpCAM protein leads to the up-regulation of c-myc, e-fabp, and cyclins A and E which promote tumorigenesis. EpCAM can act as potential diagnostic and prognostic biomarker for different types of cancers as it is also found to be expressed in epithelia and epithelial-derived neoplasms. Hence, we aimed to analyze the EpCAM gene expression and any associated feedback in the patients of two major types of lung cancer (LC) i.e., lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC), based on the publicly available online databases. In this study, server-based gene expression analysis represents the up-regulation of EpCAM in both LUAD and LUSC subtypes as compared to the corresponding normal tissues. Besides, the histological sections revealed the over-expression of EpCAM protein in cancerous tissues by depicting strong staining signals. Furthermore, mutation analysis suggested missense as the predominant type of mutation both in LUAD and LUSC in the EpCAM gene. A significant correlation (P-value < 0.05) between the higher EpCAM expression and lower patient survival was also found in this study. Finally, the co-expressed genes were identified with their ontological features and signaling pathways associated in LC development. The overall study suggests EpCAM to be a significant biomarker for human LC prognosis.     

Identification of a novel glycolysis-related gene signature that can predict the survival of patients with lung adenocarcinoma

Lung cancer is one of the most malignant cancers worldwide, and lung adenocarcinoma (LUAD) is the most common histologic subtype. Thousands of biomarkers related to the survival and prognosis of patients with this cancer type have been investigated through database mining; however, the prediction effect of a single gene biomarker is not satisfactorily specific or sensitive. Thus, the present study aimed to develop a novel gene signature of prognostic values for patients with LUAD. Using a data-mining method, we performed expression profiling of 1145 mRNAs in large cohorts with LUAD (n = 511) from The Cancer Genome Atlas database. Using the Gene Set Enrichment Analysis, we selected 198 genes related to GLYCOLYSIS, which is the most important enrichment gene set. Moreover, these genes were identified using Cox proportional regression modeling. We established a risk score staging system to predict the outcome of patients with LUAD and subsequently identified four genes (AGRN, AKR1A1, DDIT4, and HMMR) that were closely related to the prognosis of patients with LUAD. The identified genes allowed us to classify patients into the high-risk group (with poor outcome) and low-risk group (with better outcome). Compared with other clinical factors, the risk score has a better performance in predicting the outcome of patients with LUAD, particularly in the early stage of LUAD. In conclusion, we developed a four-gene signature related to glycolysis by utilizing the Cox regression model and a risk staging model for LUAD, which might prove valuable for the clinical management of patients with LUAD.     

Epigenetic loss of putative tumor suppressor SFRP3 correlates with poor prognosis of lung adenocarcinoma patients

Secreted frizzled related protein 3 (SFRP3) contains a cysteine-rich domain (CRD) that shares homology with Frizzled CRD and regulates WNT signaling. Independent studies showed epigenetic silencing of SFRP3 in melanoma and hepatocellular carcinoma. Moreover, a tumor suppressive function of SFRP3 was shown in androgen-independent prostate and gastric cancer cells. The current study is the first to investigate SFRP3 expression and its potential clinical impact on non-small cell lung carcinoma (NSCLC). WNT signaling components present on NSCLC subtypes were preliminary elucidated by expression data of The Cancer Genome Atlas (TCGA). We identified a distinct expression signature of relevant WNT signaling components that differ between adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC). Of interest, canonical WNT signaling is predominant in LUAD samples and non-canonical WNT signaling is predominant in LUSC. In line, high SFRP3 expression resulted in beneficial clinical outcome for LUAD but not for LUSC patients. Furthermore, SFRP3 mRNA expression was significantly decreased in NSCLC tissue compared to normal lung samples. TCGA data verified the reduction of SFRP3 in LUAD and LUSC patients. Moreover, DNA hypermethylation of SFRP3 was evaluated in the TCGA methylation dataset resulting in epigenetic inactivation of SFRP3 expression in LUAD, but not in LUSC, and was validated by pyrosequencing of our NSCLC tissue cohort and in vitro demethylation experiments. Immunohistochemistry confirmed SFRP3 protein downregulation in primary NSCLC and indicated abundant expression in normal lung tissue. Two adenocarcinoma gain-of-function models were used to analyze the functional impact of SFRP3 on cell proliferation and regulation of CyclinD1 expression in vitro. Our results indicate that SFRP3 acts as a novel putative tumor suppressor gene in adenocarcinoma of the lung possibly regulating canonical WNT signaling.     

Exploration of the prognostic prediction value of the PANoptosis-based risk score and its correlation with tumor immunity in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a common cancer with high mortality worldwide. PANoptosis is a novel inflammatory programmed cell death modality with the characteristics of pyroptosis, apoptosis and necroptosis. It is necessary to explore PANoptosis-related genes in LUAD patients and offer evidence for prognosis prediction and therapeutic strategies. Single-cell RNA sequencing data and RNA expression profiles of LUAD patients from The Cancer Genome Atlas and Gene Expression Omnibus databases are used to screen PANoptosis-related differential genes for the construction of a risk model. Fifteen PANoptosis-related markers with prognostic value were identified by Least Absolute Shrinkage and Selection Operator (LASSO)–Cox regression analysis. Kaplan–Meier analysis and receiver operating characteristic curve analysis further demonstrated the significant predictive capability. Immune infiltration, Single Nucleotide Variants (SNV) mutations, and clinical drug susceptibility were analyzed. In conclusion, a risk model of 15 PANoptosis-related genes has significant value in prognostic prediction for LUAD and has potential to direct clinical therapeutic strategies during the treatment.     

Identification of a panel of mitotic spindle-related genes as a signature predicting survival in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the most malignant tumor types worldwide. Our objective was to identify a genetic signature that could predict the prognosis of patients with LUAD. We extracted gene data sets from The Cancer Genome Atlas and obtained differentially expressed genes that were highly expressed at every stage. These genes were analyzed using gene set enrichment analysis to obtain four biological processes associated with LUAD. Subsequently, Cox univariate and multivariate analyses were performed to generate four optimized models (G2M checkpoint, E2F targets, mitotic spindle, and glycolysis). We identified a mitotic spindle-related signature (KIF15, BUB1, CCNB2, CDK1, KIF4A, DLGAP5, ECT2, and ANLN), which could be an independent prognostic indicator, to predict the prognosis of patients with LUAD. This new discovery should offer opportunities to explore the pathogenesis of LUAD and prove clinically useful in predicting LUAD patient prognosis.     

Integrated analysis of gene expression and DNA methylation datasets identified key genes and a 6-gene prognostic signature for primary lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the main subtype of non-small cell lung cancer with a poor survival prognosis. In our study, gene expression, DNA methylation, and clinicopathological data of primary LUAD were utilized to identify potential prognostic markers for LUAD, which were recruited from The Cancer Genome Atlas (TCGA) database. Univariate regression analysis showed that there were 21 methylation-associated DEGs related to overall survival (OS), including 9 down- and 12 up-regulated genes. The 12 up-regulated genes with hypomethylation may be risky genes, whereas the other 9 down-regulated genes with hypermethylation might be protective genes. By using the Step-wise multivariate Cox analysis, a methylation-associated 6-gene (consisting of CCL20, F2, GNPNAT1, NT5E, B3GALT2, and VSIG2) prognostic signature was constructed and the risk score based on this gene signature classified patients into high- or low-risk groups. Patients of the high-risk group had shorter OS than those of the low-risk group in both the training and validation cohort. Multivariate Cox analysis and the stratified analysis revealed that the risk score was an independent prognostic factor for LUAD patients. The methylation-associated gene signature may serve as a prognostic factor for LUAD patients and the represent hypermethylated or hypomethylated genes might be potential targets for LUAD therapy.     

DNA methylation profiling to predict recurrence risk in stage Ι lung adenocarcinoma: Development and validation of a nomogram to clinical management

Increasing evidence suggested DNA methylation may serve as potential prognostic biomarkers; however, few related DNA methylation signatures have been established for prediction of lung cancer prognosis. We aimed at developing DNA methylation signature to improve prognosis prediction of stage I lung adenocarcinoma (LUAD). A total of 268 stage I LUAD patients from the Cancer Genome Atlas (TCGA) database were included. These patients were separated into training and internal validation datasets. GSE39279 was used as an external validation set. A 13‐DNA methylation signature was identified to be crucially relevant to the relapse‐free survival (RFS) of patients with stage I LUAD by the univariate Cox proportional hazard analysis and the least absolute shrinkage and selection operator (LASSO) Cox regression analysis and multivariate Cox proportional hazard analysis in the training dataset. The Kaplan‐Meier analysis indicated that the 13‐DNA methylation signature could significantly distinguish the high‐ and low‐risk patients in entire TCGA dataset, internal validation and external validation datasets. The receiver operating characteristic (ROC) analysis further verified that the 13‐DNA methylation signature had a better value to predict the RFS of stage I LUAD patients in internal validation, external validation and entire TCGA datasets. In addition, a nomogram combining methylomic risk scores with other clinicopathological factors was performed and the result suggested the good predictive value of the nomogram. In conclusion, we successfully built a DNA methylation‐associated nomogram, enabling prediction of the RFS of patients with stage I LUAD.     

Systematic construction and validation of an immune prognostic model for lung adenocarcinoma

Lung adenocarcinoma (LUAD), the most common non‐small‐cell lung cancer, is characterized by a dense lymphocytic infiltrate, which indicates that the immune system plays an active role in the development and growth of this cancer. However, no investigations to date have proposed robust models for predicting survival outcome for patients with LUAD in terms of tumour immunology. A total of 761 LUAD patients were included in this study, in which the database of The Cancer Genome Atlas (TCGA) was utilized for discovery, and the Gene Expression Omnibus (GEO) database was utilized for validation. Bioinformatics analysis and R language tools were utilized to construct an immune prognostic model and annotate biological functions. Lung adenocarcinoma showed a weakened immune phenotype compared with adjacent normal tissues. Immune‐related gene sets were profiled, an immune prognostic model based on 2 immune genes (ANLN and F2) was developed with the TCGA database to distinguish cases as having a low or high risk of unfavourable prognosis, and the model was verified with the GEO database. The model was prognostically significant in stratified cohorts, including stage I‐II, stage III‐IV and epidermal growth factor receptor (EGFR) mutant subsets, and was considered to be an independent prognostic factor for LUAD. Furthermore, the low‐ and high‐risk groups showed marked differences in tumour‐infiltrating leucocytes, tumour mutation burden, aneuploidy and PD‐L1 expression. In conclusion, an immune prognostic model was proposed for LUAD that is capable of independently identifying patients at high risk for poor survival, suggesting a relationship between local immune status and prognosis.     

Molecular characterization of lung cancer: A two-miRNA prognostic signature based on cancer stem-like cells related genes

Lung cancer is one of the deadliest cancers worldwide. To increase the survival rate of lung cancer, it is necessary to explore specific prognosis markers. More and more evidence finds that noncoding RNA is closely associated with the survival of lung cancer, and cancer stem cells (CSCs) also play a significant role in the progress of lung cancer. The objective of this study is to find CSLCs genes that affect the prognosis of lung cancer. The differential expression of long noncoding RNAs (lncRNAs), microRNAs (miRNAs), messenger RNAs (mRNAs) in the Cancer Genome Atlas (TCGA) database and differential expression data from microarray of CD326⁺ and CD326⁻ A549 cell are intersected to identify stable and consistent expression genes (2 lncRNAs, 15 miRNAs, and 134 mRNAs). The intersection of lncRNAs and miRNAs is analyzed by univariate and multivariate Cox regression to obtained prognostic genes. Two miRNAs (miR-30b-5p and miR-29c-3p) are significantly correlated with the overall survival rate. Then using these two miRNAs to construct a risk score model as a prognosis signature of lung cancer. Subsequently, we analyzed the association between two miRNAs and clinical information of lung cancer patients, of which T stage, Neoplasm cancer and risk score (P < .05) can be used as independent prognostic indicators of lung cancer. Finally, target genes of 2 miRNAs and 134 mRNAs were annotated with Gene Ontology and analyzed with Kyoto Encyclopedia of Genes and Genomes pathway, and verified with the GEO database. In summary, this study illustrates the role of miRNAs in the promotion of lung cancer by CSCs, which is important to find molecular biomarkers of lung cancer.     

Systematic profiling of invasion-related gene signature predicts prognostic features of lung adenocarcinoma

Due to the high heterogeneity of lung adenocarcinoma (LUAD), molecular subtype based on gene expression profiles is of great significance for diagnosis and prognosis prediction in patients with LUAD. Invasion-related genes were obtained from the CancerSEA database, and LUAD expression profiles were downloaded from The Cancer Genome Atlas. The ConsensusClusterPlus was used to obtain molecular subtypes based on invasion-related genes. The limma software package was used to identify differentially expressed genes (DEGs). A multi-gene risk model was constructed by Lasso-Cox analysis. A nomogram was also constructed based on risk scores and meaningful clinical features. 3 subtypes (C1, C2 and C3) based on the expression of 97 invasion-related genes were obtained. C3 had the worst prognosis. A total of 669 DEGs were identified among the subtypes. Pathway enrichment analysis results showed that the DEGs were mainly enriched in the cell cycle, DNA replication, the p53 signalling pathway and other tumour-related pathways. A 5-gene signature (KRT6A, MELTF, IRX5, MS4A1 and CRTAC1) was identified by using Lasso-Cox analysis. The training, validation and external independent cohorts proved that the model was robust and had better prediction ability than other lung cancer models. The gene expression results showed that the expression levels of MS4A1 and KRT6A in tumour tissues were higher than in normal tissues, while CRTAC1 expression in tumour tissues was lower than in normal tissues. The 5-gene signature prognostic stratification system based on invasion-related genes could be used to assess prognostic risk in patients with LUAD.