Single-cell sequencing reveals MYC targeting gene MAD2L1 is associated with prostate cancer bone metastasis tumor dormancy.

Background: Among malignant tumors, bone metastasis is frequently associated with prostate cancer which is seen in about 80% of patients. During cancer treatments, some tumor cells switch to a "dormant mode" to help tumor cells avoid attack from the immune system and anti-tumor therapies. In this dormant mode, tumor cells can be resuscitated, causing cancer to reoccur. The generally accepted explanation for this phenomenon is that the tumor cells have spread to the bone marrow before treatment and are dormant in the bone marrow. However, the key mechanism for inducing and maintaining the dormancy of these prostate cancer disseminated tumor cells in the bone marrow is still unclear. Therefore, studying the dormancy mechanism of tumor cells in bone metastasis is of great significance for the treatment and the prevention of recurrence of prostate cancer.

Methods: We obtained single-cell RNA-seq data of tumors from mouse models of prostate cancer bone metastasis mouse model numbered (GSE147150) from the GEO database, and obtained RNA-seq expression data and clinical information from The Cancer Genome Atlas Program (TCGA) of prostate cancer patients from the USCS Xena database. Screening of differential genes and annotation of GO functions were performed separately. Subsequently, the screened differential genes were compared and analyzed with 50 classic Hallmark signaling pathways, and the prognosis analysis of prostate cancer patients in TCGA data was performed to discover the key genes of the dormant mechanism of tumor cells in bone metastasis, and obtain new biomarkers that can be used to predict the prognosis of patients.

Results: A total of 378 differentially expressed genes were screened, of which 293 were significantly up-regulated and 85 were significantly down-regulated. Among them, the up-regulated genes were mainly related to the immune response, and the down-regulated genes were mainly related to the cell cycle. Through GSVA (Gene set variation analysis), it is found that there are differences in a total of 3 signal pathways: COMPLEMENT, MYC_TARGETS_V1 and MYC_TARGETS_V2. By comparing and analyzing the significantly down-regulated genes in dormant tumor cells with MYC_TARGETS_V1, MYC_TARGETS_V2, three significantly down-regulated genes were obtained: Ccna2, Mad2L1 and Plk1.

Conclusion: In summary, our findings indicate that the MYC targeting gene Mad2L1 is potentially related to the dormancy mechanism of prostate cancer. At the same time, Mad2L1, a gene associated with dormant prostate cancer cells, may be used as a biomarker for prognostic survival.