Cutaneous squamous cell carcinoma (cSCC) is the leading cause of death in patients with nonmelanoma skin cancers (NMSC). However, the unclear pathogenesis of cSCC limits the application of molecular targeted therapy.Three microarray datasets (GSE2503, GSE45164, and GSE66359) were downloaded from the Gene Expression Omnibus (GEO). After identifying the differentially expressed genes (DEGs) in tumor and nontumor tissues, five kinds of analyses, namely, functional annotation, protein-protein interaction (PPI) network, hub gene selection, TF-miRNA-mRNA regulatory network analysis, and ferroptosis mechanism, were performed.A total of 146 DEGs were identified with significant differences, including 113 upregulated genes and 33 downregulated genes. The enriched functions and pathways of the DEGs included microtubule-based movement, ATP binding, cell cycle, P53 signaling pathway, oocyte meiosis, and PLK1 signaling events. Nine hub genes were identified (CDK1, AURKA, RRM2, CENPE, CCNB1, KIAA0101, ZWINT, TOP2A, and ASPM). Finally, RRM2, AURKA, and SAT1 were identified as significant ferroptosis-related genes in cSCC. The differential expression of these genes has been verified in two other independent datasets.By integrated bioinformatic analysis, the hub genes identified in this study elucidated the molecular mechanism of the pathogenesis and progression of cSCC and are expected to become future biomarkers or therapeutic targets.Copyright © 2022 Wenxing Su et al.