Lethal genes have not been systematically analyzed in breast cancer which may have significant prognostic value. The current study aims to investigate vital genes related to cell viability by analyzing the CRISPR-cas9 screening data, which may provide novel therapeutic target for patients.
Genes differentially expressed between tumor and normal tissue from the Cancer Genome Atlas (TCGA) and genes related to cell viability by CRISPR-cas9 screening from Depmap (Cancer Dependency Map) were overlapped. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis was conducted to identify which pathways of overlapped genes were enriched. GSE21653 set was randomized into training and internal validation dataset at a ratio of 3:1, and external validation was performed in GSE20685 set. The least absolute shrinkage and selection operator (LASSO) regression was used to construct a signature to predict recurrence-free survival (RFS) of breast cancer patients. Univariate and multivariate Cox regression were used to evaluate the prognostic value of this signature. Differentially expressed genes (DEGs) between high-risk and low-risk patients were then analyzed to identify the main pathways regulated by this signature. Weighted correlation network analysis (WGCNA) was conducted to recognize modules correlated with high risk. Enrichment analysis was then used to identify pathways regulated by genes shared in the overlapped genes, DEGs, and WGCNA.
A total of 86 oncogenes were upregulated in TCGA database and overlapped with lethal genes in Depmap database, which were enriched in cell cycle pathway. A total of 51 genes were included in the gene signature based on LASSO regression, and the median risk score of 2.36 was used as cut-off to separate low-risk patients from high-risk patients. High-risk patients showed worse RFS compared with low-risk patients in internal training, internal validation, and external validation dataset. Time-dependent receiver operating characteristic curves of 3 and 5 years indicated that risk score was superior to tumor stage, age, and PAM50 in both entire and external validation datasets. Cell cycle was the main different pathway between the high-risk and low-risk groups. Meanwhile, cell cycle was also the main pathway enriched in the 25 genes which were shared among 86 genes, DEGs, and WGCNA.
Cell cycle pathway, identified by CRISPR-cas9 screening, was a key pathway regulating cell viability, which has significant prognostic values and can serve as a new target for breast cancer patient treatment.