Pathogenic genes related to the progression of actinic keratoses to cutaneous squamous cell carcinoma

Abstract

Background: Actinic keratosis (AK) is an incipient form of cutaneous squamous cell carcinoma (cSCC). Understanding the differentially expressed genes between AK and cSCC states would be helpful for the early prevention and treatment of cSCC. Consequently, this study aimed to screen the key genes associated with the progression of AK to cSCC.

Methods: The microarray dataset GSE45216 was downloaded from the Gene Expression Omnibus, which included 10 AK and 30 primary cSCC skin tissue samples. Differentially expressed genes (DEGs) in cSCC samples, compared to those in AK, were identified. Gene co-expression relationships were investigated, followed by miRNA prediction. The potential functions of the co-expressed genes were predicted by gene ontology (GO) and pathway enrichment analyses. In addition, the transcription factors and drug molecules, significantly related to the co-expressed genes, were obtained.

Results: A total of 320 DEGs were identified in the cSCC group, relative to the AK group. Moreover, 96 DEGs and 2,390 connecting edges were identified in the gene co-expression network. An miRNA regulatory network was constructed, including 96 DEGs and 16 miRNAs. In addition, three co-expression network modules were obtained; EIF4EBP1, SNX17, PRPF4, NXT1, and UBA5 were significant nodes in the modules.

Conclusions: EIF4EBP1, SNX17, PRPF4, NXT1, and UBA5 may be the pathogenic genes contributing to the development of cSCC from AK.