Dataset Information

Seven key hub genes identified by gene co-expression network in cutaneous squamous cell carcinoma.

ABSTRACT:

Background

Cutaneous squamous cell carcinoma (cSCC) often follows actinic keratosis (AK) and is the second most common skin cancer worldwide. To reduce metastasis risk, it is important to diagnose and treat cSCC early. This study aimed to identify hub genes associated with cSCC and AK.

Methods

This study used three datasets GSE45216, GSE98774, and GSE108008. We combined samples from the GSE45216 and GSE98774 datasets into the new dataset GSE45216-98774. We applied a weighted gene co-expression network analysis (WGCNA) to investigate key modules and hub genes associated with cSCC and AK. We considered the hub genes found in both the GSE45216-98774 and GSE108008 datasets as validated hub genes. We tested whether the expression of hub genes could predict patient survival outcomes in other cancers using TCGA pan-cancer data.

Results

We identified modules most relevant to cSCC and AK. Additionally, we identified and validated seven hub genes of cSCC: GATM, ARHGEF26, PTHLH, MMP1, POU2F3, MMP10 and GATA3. We did not find validated hub genes for AK. Each hub gene was significantly associated with the survival of various cancer types. Only GATA3 was significantly associated with melanoma survival.

Conclusions

We applied WGCNA to find seven hub genes that play important roles in cSCC tumorigenesis. These results provide new insights that help explain the pathogenesis of cSCC. These hub genes may become biomarkers or therapeutic targets for accurate diagnosis and treatment of cSCC in the future.

SUBMITTER: Chen H

PROVIDER: S-EPMC8306372 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundLaryngeal squamous cell carcinoma (LSCC) is the prominent cancer in head and neck, which greatly affects life quality of patients. The pathogenesis of LSCC is not clear. Presently, the LSCC treatments include chemotherapy, surgery and radiotherapy; however, these methods have poor efficacy in patients with recurrent and persistent cancer. Therefore, the study identified the hub genes accompanied with LSCC, which may be a potential therapeutic target in the future.MethodsWe extracted whole transcriptome high-throughput sequencing (HTS) LSCC data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases and calculate differentially expressed genes (DEGs) between LSCC and normal samples using statistical software RStudio. Through weighted gene co-expression network analysis (WGCNA), enrichment examination of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) functions, and examination of protein-protein interaction (PPI) network, we obtained network hub genes and validated the hub genes prognostic value and expression levels of protein.ResultsThrough analysis of differential gene expression, from the GEO and TCGA databases 2,139 and 2,774 DEGs were obtained, respectively, 13 and 15 modules were screened from TCGA-LSCC and GSE127165 datasets by WGCNA, respectively. The most significant positive and negative correlation modules in the WGCNA and DEG lists were overlapped, and overall 36 co-expressed overlapping genes were retrieved. Through enrichment analysis of GO and KEGG, it was found that the gene functions were highly concentrated in cell junction assembly, basement membrane, extracellular matrix (ECM) structural constituent etc., and the pathways were mainly concentrated in ECM receptor interaction, focal adhesion, small cell lung cancer, and toxoplasmosis. Through analysis of PPI network analysis, 10 network hub genes (SNAI2, ITGA6, LAMB3, LAMC2, CAV1, COL7A1, GJA1, EHF, OAT, and GPT) were obtained. Finally, survival analysis and protein expression validation of these genes confirmed that low OAT expression and high CAV1 expression remarkably influenced the survival of patient's prognosis with LSCC.ConclusionsWe recognized the hub genes and key modules nearly associated to LSCC and these genes were validated by survival analysis and the database of Human Protein Atlas (HPA), which is of high importance for unveiling the pathogenesis of LSCC and probing for new precise biological marker and potential therapeutic targets.

Project description:Laryngeal squamous cell carcinoma (LSCC) is the most common malignant tumor in the head and neck, and can seriously affect the daily life of patients. To study the mechanisms of LSCC, the microarray of GSE51958 was analyzed in the present study. GSE51958 was downloaded from Gene Expression Omnibus, and included a collection of LSCC tissue samples and matched adjacent non-cancerous tissue samples from 10 patients. Differentially-expressed genes (DEGs) were identified using limma package. Next, a weighted co-expression network was constructed for the DEGs by WGCNA package in R. Modules of the weighted co-expression network were obtained through constructing a hierarchical clustering tree using the hybrid dynamic shear tree method. Using the clusterProfiler package, the potential functions of DEGs in the modules correlated with LSCC were predicted by pathway enrichment analysis. In total, 959 DEGs were screened from the LSCC samples compared with the adjacent non-cancerous samples, including 553 upregulated and 406 downregulated genes. The appointed black, brown, gray, pink and yellow modules were screened for the DEGs in the weighted co-expression network. For the DEGs in the brown and yellow modules, the enriched pathways were cytokine-cytokine receptor interaction and metabolic pathways, respectively. The DEGs in the pink module were involved in the majority of pathways. With high connectivity degrees in the pink module, TPX2, microtubule-associated (TPX2; degree, 25), minichromosome maintenance complex component 2 (MCM2; degree, 25), ubiquitin-like with PHD and ring finger domains 1 (UHRF1; degree, 22), cyclin-dependent kinase 2 (CDK2; degree, 20) and protein regulator of cytokinesis 1 (PRC1; degree, 20) may be involved in LSCC. Overall, In conclusion, from the integrated bioinformatics analysis of genes that may be associated with LSCC, 959 DEGs were obtained from LSCC samples compared with adjacent non-cancerous samples, and TPX2, MCM2, UHRF1, CDK2 and PRC1 were found to hold a possible association with the disease.