Meta-Analysis of EGF-Stimulated Normal and Cancer Cell Lines to Discover EGF-Associated Oncogenic Signaling Pathways and Prognostic Biomarkers

Although epidermal growth factor (EGF) controls many crucial processes in the human body, it can increase the risk of developing cancer when overexpresses. This study focused on detecting cancer-associated genes that are dysregulated by EGF overexpression. To identify differentially expressed genes (DEGs), two independent meta-analyses with normal and cancer RNA-Seq samples treated by EGF were conducted. The new DEGs detected only via two meta-analyses were used in all downstream analyses. To reach count data, the tools of FastQC, Trimmomatic, HISAT2, SAMtools, and HTSeq-count were employed. DEGs in each individual RNA-Seq study and the meta-analysis of RNA-Seq studies were identified using DESeq2 and metaSeq R package, respectively. MCODE detected densely interconnected top clusters in the protein-protein interaction (PPI) network of DEGs obtained from normal and cancer datasets. The DEGs were then introduced to Enrichr and ClueGO/CluePedia, and terms, pathways, and hub genes enriched in Gene Ontology (GO) and KEGG and Reactome were detected. The meta-analysis of normal and cancer datasets revealed 990 and 541 new DEGs, all upregulated. A number of DEGs were enriched in protein K48-linked deubiquitination, ncRNA processing, ribosomal large subunit binding, and protein processing in endoplasmic reticulum. Hub genes overexpression (DHX33, INTS8, NMD3, OTUD4, P4HB, RPS3A, SEC13, SKP1, USP34, USP9X, and YOD1) in tumor samples were validated by TCGA and GTEx databases. Overall survival and disease-free survival analysis also confirmed worse survival in patients with hub genes overexpression. The detected hub genes could be used as cancer biomarkers when EGF overexpresses.