Identification of critically carcinogenesis-related genes in basal cell carcinoma

Background Basal cell carcinoma (BCC) is a frequent malignant tumor of skin cancers with high morbidity. The objective of this study was to identify critical genes and pathways related to the carcinogenesis of BCC and gain more insights into the underlying molecular mechanisms of BCC. Materials and methods The gene expression profiles of GSE7553 and GSE103439 were downloaded from the Gene Expression Omnibus database with 19 tumors and 6 normal skin tissues. Differentially expressed genes (DEGs) were screened between BCC samples and normal tissues, followed by gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Subsequently, protein–protein interaction (PPI) network was constructed for these DEGs, and module analysis was performed. Results A total of 313 DEGs were obtained. Among them, 222 genes were upregulated and 91 genes were downregulated. Enrichment analysis indicated that the upregulated genes were significantly enriched in cell cycle and mitosis, while the downregulated genes were mainly associated with unsaturated fatty acid metabolic process and cell differentiation. In addition, TOP2A, CDK1, and CCNB1 were identified as the top three hub genes ranked by degrees in the PPI network. Meanwhile, three subnetworks were derived, which indicated that these DEGs were significantly enriched in pathways, including “cell cycle”, “extracellular matrix–receptor interaction”, “basal cell carcinoma”, and “hedgehog signaling pathway”. Conclusions The novel critical DEGs and pathways identified in this study may serve pivotal roles in the carcinogenesis of BCC and indicate more molecular targets for the treatment of BCC.

[1]  B. N. Akay,et al.  A Novel PTCH1 Frameshift Mutation Leading to Nevoid Basal Cell Carcinoma Syndrome , 2018, Cytogenetic and Genome Research.

[2]  V. Ilankovan,et al.  Epidemiology of basal cell carcinoma: a 10-year comparative study. , 2017, The British journal of oral & maxillofacial surgery.

[3]  I. Majsterek,et al.  Relation between sonic hedgehog pathway gene polymorphisms and basal cell carcinoma development in the Polish population , 2015, Archives of Dermatological Research.

[4]  Pei-Ru Wu,et al.  High nuclear/cytoplasmic ratio of Cdk1 expression predicts poor prognosis in colorectal cancer patients , 2014, BMC Cancer.

[5]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[6]  M. Barbacid,et al.  Cell cycle, CDKs and cancer: a changing paradigm , 2009, Nature Reviews Cancer.

[7]  Gary D. Bader,et al.  An automated method for finding molecular complexes in large protein interaction networks , 2003, BMC Bioinformatics.

[8]  B. Wang,et al.  Strategies to target the Hedgehog signaling pathway for cancer therapy , 2018, Medicinal research reviews.

[9]  R. Dummer,et al.  Treatment with two different doses of sonidegib in patients with locally advanced or metastatic basal cell carcinoma (BOLT): a multicentre, randomised, double-blind phase 2 trial. , 2015, The Lancet. Oncology.

[10]  Deevya L. Narayanan,et al.  Review: Ultraviolet radiation and skin cancer , 2010, International journal of dermatology.

[11]  Damian Szklarczyk,et al.  STRING v9.1: protein-protein interaction networks, with increased coverage and integration , 2012, Nucleic Acids Res..

[12]  D. Donoghue,et al.  Patched1 interacts with cyclin B1 to regulate cell cycle progression , 2001, The EMBO journal.

[13]  Juancarlos Chan,et al.  Gene Ontology Consortium: going forward , 2014, Nucleic Acids Res..

[14]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[15]  Pierre Dubus,et al.  Cdk1 is sufficient to drive the mammalian cell cycle , 2007, Nature.

[16]  Martin A Weinstock,et al.  Non-melanoma skin cancer incidence and impact of skin cancer screening on incidence. , 2014, The Journal of investigative dermatology.

[17]  S. Vranić,et al.  The role of the Hedgehog signaling pathway in cancer: A comprehensive review. , 2017, Bosnian journal of basic medical sciences.

[18]  H. Hiasa,et al.  Topoisomerases as anticancer targets. , 2018, The Biochemical journal.

[19]  J. Molnár,et al.  Selective in situ protein expression profiles correlate with distinct phenotypes of basal cell carcinoma and squamous cell carcinoma of the skin. , 2013, Histology and histopathology.

[20]  S. Antonarakis,et al.  Genomic analysis identifies new drivers and progression pathways in skin basal cell carcinoma , 2016, Nature Genetics.

[21]  Benjamin M. Bolstad,et al.  affy - analysis of Affymetrix GeneChip data at the probe level , 2004, Bioinform..

[22]  P. Harms,et al.  Tumor Architecture and Notch Signaling Modulate Drug Response in Basal Cell Carcinoma. , 2018, Cancer cell.

[23]  B. Wainwright,et al.  Patched1 functions as a gatekeeper by promoting cell cycle progression. , 2006, Cancer research.

[24]  Zhi-hua Chen,et al.  Kyoto Encyclopedia of Genes and Genomes were used for functional enrichment analysis of differentially expressed genes (DEGs). A protein‐protein interaction network was constructed, and the hub genes were subjected to module analysis and identification using Search Tool for the Retrieval , 2019 .

[25]  J. Edwards,et al.  Androgen receptor phosphorylation at serine 515 by Cdk1 predicts biochemical relapse in prostate cancer patients , 2012, British Journal of Cancer.

[26]  C. Harris,et al.  Loss of Cytoplasmic CDK1 Predicts Poor Survival in Human Lung Cancer and Confers Chemotherapeutic Resistance , 2011, PloS one.

[27]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[28]  S. Eschrich,et al.  The gene expression profiles of primary and metastatic melanoma yields a transition point of tumor progression and metastasis , 2008, BMC Medical Genomics.

[29]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[30]  Enhao Fang,et al.  Identification of breast cancer hub genes and analysis of prognostic values using integrated bioinformatics analysis. , 2017, Cancer biomarkers : section A of Disease markers.

[31]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[32]  R. Bristow,et al.  AZD5438, an inhibitor of Cdk1, 2, and 9, enhances the radiosensitivity of non-small cell lung carcinoma cells. , 2012, International journal of radiation oncology, biology, physics.

[33]  Juan Wang,et al.  Exploration of the mechanism of colorectal cancer metastasis using microarray analysis , 2017, Oncology letters.

[34]  M. Nihal,et al.  Polo-like kinase 1 (Plk1) in non-melanoma skin cancers , 2009, Cell cycle.