Non-melanoma skin cancers (NMSC) including basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are more common kinds of skin cancer. Although these tumors share common pathological and clinical features, their similarity and heterogeneity at molecular levels are not fully elaborated yet. Here, by performing comparative analysis of gene expression profiling of BCC, SCC, and normal skin tissues, we could classify the BCC into three subtypes of classical, SCC-like, and normal-like BCCs. Functional enrichment and pathway analyses revealed the molecular characteristics of each subtype. The classical BCC showed the enriched expression and transcription signature with the activation of Wnt and Hedgehog signaling pathways, which were well known key features of BCC. By contrast, the SCC-like BCC was enriched with immune-response genes and oxidative stress-related genes. Network analysis revealed the PLAU/PLAUR as a key regulator of SCC-like BCC. The normal-like BCC showed prominent activation of metabolic processes particularly the fatty acid metabolism. The existence of these molecular subtypes could be validated in an independent dataset, which demonstrated the three subgroups of BCC with distinct functional enrichment. In conclusion, we suggest a novel molecular classification of BCC providing insights on the heterogeneous progression of BCC.
