- Open Access
High copy number variations, particular transcription factors, and low immunity contribute to the stemness of prostate cancer cells
Journal of Translational Medicine volume 19, Article number: 206 (2021)
Tumor metastasis is the main cause of death of cancer patients, and cancer stem cells (CSCs) is the basis of tumor metastasis. However, systematic analysis of the stemness of prostate cancer cells is still not abundant. In this study, we explore the effective factors related to the stemness of prostate cancer cells by comprehensively mining the multi-omics data from TCGA database.
Based on the prostate cancer transcriptome data in TCGA, gene expression modules that strongly relate to the stemness of prostate cancer cells are obtained with WGCNA and stemness scores. Copy number variation of stemness genes of prostate cancer is calculated and the difference of transcription factors between prostate cancer and normal tissues is evaluated by using CNV (copy number variation) data and ATAC-seq data. The protein interaction network of stemness genes in prostate cancer is constructed using the STRING database. Meanwhile, the correlation between stemness genes of prostate cancer and immune cells is analyzed.
Prostate cancer with higher Gleason grade possesses higher cell stemness. The gene set highly related to prostate cancer stemness has higher CNV in prostate cancer samples than that in normal samples. Although the transcription factors of stemness genes have similar expressions, they have different contributions between normal and prostate cancer tissues; and particular transcription factors enhance the stemness of prostate cancer, such as PUM1, CLOCK, SP1, TCF12, and so on. In addition, the lower tumor immune microenvironment is conducive to the stemness of prostate cancer. CD8 + T cells and M1 macrophages may play more important role in the stemness of prostate cancer than other immune cells in the tumor microenvironment. Finally, EZH2 is found to play a central role in stemness genes and is negatively correlated with resting mast cells and positively correlated with activated memory CD4 + T cells.
Based on the systematic and combined analysis of multi-omics data, we find that high copy number variation, specific transcription factors, and low immune microenvironment jointly contribute to the stemness of prostate cancer cells. These findings may provide us new clues and directions for the future research on stemness of prostate cancer.
Cancer stem cells (CSCs) are a few stemness-like cells with the ability of self-renewal and differentiation in cancers . They play an important role in the occurrence and development of cancers, especially closely related to cancer metastasis [2,3,4]. In prostate cancer, the stemness of cancer cells (including prostate cancer stem cells, PCSC) is closely related to the metastasis of prostate cancer . In the metastasis process of prostate cancer, PCSCs initiates EMT (epithelial-mesenchymal transition) to form fibroblast-like cells and then enter the blood. With the circulation system, prostate cancer cells migrate to other tissues (such as bone tissue and lymph tissue) and grow into tumor tissues, which leads to tumor metastasis (cancer cell spreading).
As we known, there are many factors which relate to the stemness of CSC cells, including both intracellular factors (such as stemness-related genes) and microenvironment of cancer tissues (such as immune cells in the tumor microenvironment) [3, 6, 7]. In prostate cancer (PCa), it has been reported that the immune cells (especially CD8 + T cells and macrophages) in the microenvironment of PCa are closely related to the metastasis of PCa cells [8, 9]. The number of immune cells around the early PCa tissue will be decreased with the growth of the cancer tissue, which results in the decrease of the immunity in the microenvironment of PCa . With the development of PCa to the later stage (Gleason score to 6–10), some immune cells (such as related T cells) in the microenvironment can reverse to promote or enhance the growth of PCa and help cancer cell metastasis .
Although increasing evidences have shown the relationship between stemness and metastasis in PCa cells, few studies on the stemness regulation of PCa cells are reported [12, 13]. Lots of factors affecting the stemness of PCa cells remain unclear and need to be investigated.
The bioinformatics method based on the TCGA database has been increasingly used to analyze the molecular basis of prostate cancer development and clinical patient prognosis [14,15,16]. Using appropriate analysis software and methods to explore a variety of large-size data of clinical specimen from the TCGA database (including transcriptome sequencing data, gene sequencing data, ATAC-seq data, etc.), the molecular basis of prostate tumorigenesis, development of PCa, and the prognosis of patients may be figured out [17, 18]. Therefore, the bioinformatics analysis of TCGA data can provide clues and directions for both the basic experimental research and the clinical cancer treatment in PCa.
In this study, we firstly obtained the most important stemness-related modules and genes in prostate cancer cells from transcriptome data and OCLR scores (a method for scoring the stemness of tumors) . Further analysis showed that stemness-related genes had higher CNV in prostate cancer than that in prostate normal samples. Transcription factors of stemness genes enhanced the stemness of prostate cancer cells. The stemness of prostate cancer cells was negatively correlated with the immune response; and low immune scores were beneficial to the prostate cancer stemness. EZH2 was found to play a central role in these stemness genes. All our multi-omics analyzing results might provide some theoretical clues for us to experimentally investigate the factors affecting the PCa cell stemness and its relationships between PCSC and PCa metastasis.
Analysis of clinical data
We obtained the clinical information of prostate cancer from TCGA database and then divided prostate cancer samples into five grades according to the Gleason score and named Gleason grade 1 (Gleason score 6), Gleason grade 2 (Gleason score (3 + 4)), Gleason grade 3 (Gleason score (4 + 3)), Gleason grade 4 (Gleason score 8), Gleason grade 5 (Gleason score 9 or 10), respectively.
Analysis of transcriptome (RNA-seq) data
Based on OCLR stemness scores and Gleason classification of TCGA prostate cancer clinical data, correlation analysis between PCa cell stemness and Gleason grade was carried out. By combing the stemness score, WGCNA  analysis was performed on the transcriptome data of prostate cancer in TCGA. Differential expression analysis of genes in the WGCNA results that mostly related to PCa cell stemness was performed and presented in a heatmap. Further, the normalized expression of stemness genes was obtained by using between-array normalization of the limma  package to analyze the transcriptome data of 33 samples (GSE104786) from GEO database and then presented in a heatmap. Also, the transcriptome data from SRA database (Normal samples SRR7651698, SRR7651699, SRR7651700), PCa cell lines (SRR7651715, SRR7651716, SRR7651717, SRR7651718), and other small cell prostate cancer (SRR7651719, SRR7651720) were aligned and quantified by using HISAT2  and HTSeq , respectively. Finally, the normalized expression data was obtained by using TPM (Transcript per million) normalization method.
Analysis of gene CNV data
The genes most related to stemness in WGCNA results were screened, and their locations in genome were obtained by local Perl script method. Combined with the CNV data of TCGA, the local Perl script was used to screen the segments containing the locations of important stemness genes. The prostate cancer samples were classified by Gleason grade and then the CNVs of essential stemness genes in each type of prostate cancer sample were calculated by using GISTIC2.0 .
Analysis of ATAC-seq data
From TCGA and SRA database, we got the ATAC-seq data, including normal prostate samples (SRR7651660, SRR7651661, SRR7651662), prostate cancer samples (SRR7651675, SRR7651676), and TCGA ATAC-seq alignment results of prostate cancer samples. Then 2 kb data of upstream and downstream of the important stemness genes were obtained by using Bowtie2 software . According to the results, the TSS signal intensity of stemness genes was drown by using deeptools . After analyzing the above alignment data using MACS  and HOMER , information about transcription regulators of the stemness genes was obtained; and the importance of each transcriptional factor was further obtained from TCGA transcriptome data by PCA analysis with R language. For example, the upstream transcriptional factor of EZH2 was drawn using Sushi  package in R language.
Analysis of tumor microenvironment and immune infiltration
The tumor microenvironment score and tumor immune infiltration score were calculated by using ESTIMATE  and CIBERSORT , respectively. Then, the correlation between stemness score and the immune score of PCa cells was calculated with R by combining with Gleason grade of prostate cancer. Based on CIBERSORT analysis results, the distribution map of immune cell components of different Gleason grades in prostate cancer was drawn with R language. The correlation network among immune cells was also drawn with R language according to the correlation and significance of different immune cells.
Protein interaction network and the correlation between stemness genes and immune infiltration
After getting and analyzing PPI (Protein–Protein Interaction) data from the STRING (functional protein association networks) database, the PPI data of essential stemness genes were obtained. The interaction network of proteins of essential stemness genes was drawn with Cytoscape . GO and KEGG enrichment analysis of important stemness genes were also carried out by using ClusterProfiler . Finally, by using R language to jointly analyze and calculate the transcription data of essential stemness genes and the consequences of immune infiltration, the correlation between the critical stemness genes and immune cells was obtained.
Stemness-related key genes are highly expressed in high Gleason grade prostate cancer samples
Based on the joint analysis of the OCLR's results and transcriptome data of PCa clinical specimen in TCGA, it’s found that the stemness of PCa cells of Gleason grade 3–5 was significantly higher than that of Gleason grade 1–2 (Fig. 1a), Besides, by using WGCNA (Weighted Gene Co-expression Network Analysis) to analyze the PCa transcriptome data for cancer cell stemness-related gene transcripts, we got 30 transcript modules related to PCa cell stemness (Fig. 1b), and samples were clustered into different modules (Additional file 1: Fig. S1a.). About 88.2% of the 2158 known cell stemness-related genes  were found in the WGCNA analysis results (Additional file 2: Table S1; Additional file 1: Fig. S1c). Among these modules, we found that the MEmagenta module had the most positive correlation with the stemness of prostate cancer (Fig. 1b: Additional file 1: Fig. S1b), suggesting that the genes in the MEmagenta module (Additional file 3: Table S2) might be important and play a vital role in the stemness of prostate cancer.
Further analysis showed that the expression of genes in the MEmagenta module in PCa cancer specimen was generally higher than that in the normal prostate specimen; and the expression of these stemness genes in PCa increased with Gleason grade, and the highest expression was found at Gleason grade 5 (Fig. 1c). In addition, we analyzed the data from 33 samples of GEO database and found that expressions of the most stemness genes in small cell prostate cancer (SCPC) were higher than those in prostate adenocarcinoma (PRAD) (Additional file 1: Fig. S1d). Meanwhile, analysis of the transcriptome data of different prostate (normal or PRAD) cell lines and SCPC cell lines was also carried out, and it's found that the expression of stemness gene in PCa cells (PRAD cell lines) was higher than that in normal cells, as well as the expression of stemness genes in SCPC cell lines was slightly higher than in PRAD cell lines (Additional file 1: Fig. S1e).
In the genes of the MEmagenta module (Additional file 3: Table S2), we found that the known genes closely related to cancer stemness in non-prostate cancer, such as BRCA1 , EZH2 , FOXM1 , CDC20 , and CDCA8 , were all clustered in this module, indicating these genes were also closely related to the stemness of prostate cancer cells. From the data of Additional file 3, it’s confirmed that EZH2 was the most significant gene in the correlation with the stemness of PCa cells.
The CNV of stemness-related genes was elevated in prostate cancer
The segments containing all MEmagenta module genes were obtained from TCGA and then analyzed by using GISTIC (Genome Identification of Significant Targets in Cancer). Analytical results showed that most of MEmagenta module genes had higher variations in prostate cancer than those in normal prostate tissue, indicating that the stemness of PCa cells was influenced by the CNVs (Copy Number Variations) of these genes (Fig. 2a and Additional file 4: Fig. S2a). In the two kinds of gene CNVs, amplification and deletion (Additional file 5: Table S3 and Additional file 6: Table S4), the changes of deletion CNVs were more than those of amplification CNVs in the MEmagenta module genes of PCa samples; and both amplification and deletion CNVs in PCa samples were much higher than those in prostate normal samples (Fig. 2b, c, and Additional file 4: Fig. S2b, c). By jointly analyzing GISTIC results with Gleason grade of PCa, it’s found that CNVs of MEmagenta module genes also increased with the Gleason grade (Fig. 2d), suggesting that CNVs of MEmagenta module genes were highly related to the malignancy of PCa.
In details of CNVs of some genes in the MEmagenta module, it was found that the CNVs (amplification) and expression of SKA3 and RUVBL1, which promoted tumor metastasis and played a role in the development of stem cells [37,38,39,40], were all increased with the Gleason grade (Additional file 7: Fig. S3a. b); and the CNVs (deletion) and expression of MCM6 and CENPH, which enhanced cancer cell proliferation, stemness, and metastasis and promoted cancer development [41,42,43,44,45], were also increased with the increase of Gleason grade (Additional file 7: Fig. S3c, d).
Transcription factors of stemness-related genes enhanced the stemness of PCa cells
We downloaded the ATAC-seq data of the MEmagenta module genes in PCa samples from TCGA and analyzed the sequence data of the 2 kb range of the transcription start site (TSS). From the results, although both normal and tumor prostate samples had the open chromatin signal, the open chromatin signal of all PCa samples (including the different Gleason grades of PCa and small cell PCa) was weaker than that of the normal prostate samples (Fig. 3a); and the open chromatin signal in small cell PCa was the weakest one (Fig. 3a). These results indicated that the open chromatin signal of the MEmagenta module genes in PCa might differ from those in normal prostate tissues.
According to the ATAC-seq data analysis results, it was found that most transcription factors of stemness genes in prostate samples were the same (Additional file 8: Fig. S4b). After filtering the 238 common transcription factors according to the genome information and transcriptome data, the expression profile of 163 transcription factors was finally obtained. According to the expression of the 163 transcription factors, 3 clusters were obtained by clustering the transcriptome data of the 163 transcription factors (Additional file 8: Fig. S4a). The first cluster had 86 transcription factors (such as FOXA1, HOXB13, ERG1, et al.), the second cluster had 44 transcription factors (such as SOX15, FOXM1, SOX2, et al.), and the third cluster had 33 transcription factors (such as NANOG, SOX3, FOXA2, et al.) (Additional file 8: Fig. S4a). In the top ten results of principal component analysis (PCA) of the 163 transcription factors, stemness-related PUM1 [46, 47], CLOCK , SP1, and TCF12 played a major positive regulatory role in PCa cell stemness (Fig. 3b and Additional file 8: Fig. S4c). In comparison, IRF3  negatively correlated with other 9 transcription factors and played a negative regulatory role in PCa cell stemness (Fig. 3b). As we know, IRF3 (interferon regulatory factor 3) signaling played an essential role in TLR3-mediated apoptosis in LNCaP cells through the activation of the intrinsic and extrinsic apoptotic pathways , suggesting that the immune system might play a role in suppressing the stemness of PCa cells. Furthermore, our analysis results showed that AR, FOXA1, NFYA, CTCF, and FOXO3 might enhance the stemness of PCa cells, while FOXF1 might be negatively correlated with these transcription factors (Fig. 3c).
In normal prostate samples, we found that the major transcription regulators of stemness genes were HOXB4, NFYB, and TFE3 (Additional file 9: Fig. S5a, c), which differed from those in prostate cancer. From the PCA analysis results, we found that the roles of FOXA1, NFYA, and FOXP1 in regulating stemness genes in normal prostate samples were changed by comparing those in prostate cancer samples (Additional file 9: Fig. S5b and Fig. 3c). These results indicated that the same transcription factors might play different roles in regulating the cell stemness in prostate normal and cancer samples.
Our results of transcriptional regulator analysis showed that the upstream of the EZH2, the most relevant gene to the stemness of PCa cells, could be significantly bound by NFY (Fig. 3d). As a transcription factor, NFY regulates the self-renewal of hematopoietic stem cells  and promotes the self-renewal and expansion of prostate cancer cells and their stemness . Hence EZH2 might play a stemness role in prostate cancer.
Immunological microenvironment negatively related to the stemness of PCa cells
By scoring the stemness and immunity of the immune microenvironment of PCa with different Gleason grades in clinical samples, there was a negative correlation between the PCa cell stemness and the immunity of the microenvironment of PCa in all clinical samples. The correlation coefficient between stemness and immunity of the immune microenvironment in PCa with Gleason grade 1–2 was almost the same as that in PCa with Gleason grade 3–5 (Fig. 4a). By analyzing and scoring the stromal and immune cells in the PCa microenvironment, we found that the scores of stromal and immune cells in the PCa microenvironment were all inversely related to the stemness of PCa cells (Additional file 10: Fig. S6a). Although the immune score generally increased with the increase of Gleason grade, the stemness of PCa samples with low-immunity was significantly higher than that of PCa samples with high-immunity in the same Gleason grade PCa samples (Additional file 10: Fig. S6b). After clustering the MEmangenta module genes based on cell stemness and immune scores of PCa, we found that PCa cell stemness was negatively related to the immunity of the PCa microenvironment in the clinical samples with high expression of the MEmagenta module genes (Additional file 10: Fig. S6d).
By analyzing and scoring the transcriptome data from TCGA for the immune infiltration of PCa clinical samples, the number of 22 types of immune cells in prostate cancer with different Gleason grades was obtained (Fig. 4b). From the results, it’s found that most types of immune cells in the microenvironment of PCa were increased with the Gleason grade increase; and both CD8 + T cells  and macrophage M1  were significantly increased in all types of immune cells, while the plasma cells  were reduced with the Gleason grade increase (Fig. 4b), indicating that plasma cells (B cells) in the microenvironment of PCa played a vital role in anti-PCa immunity.
Based on the correlation between different types of immune cells, immune cells were divided into 4 clusters by using clustering method. There were 15 immune cells in cluster-A, 3 immune cells in cluster-B, 2 immune cells in cluster-C and 2 immune cells in cluster-D; and correlation between different immune cells in cluster-A was more than that in other clusters, suggesting that cluster-A was more complicated in this immune cell network (Fig. 4c and Additional file 11: Table S5). For example, there was a negative correlation between activated NK cells and resting NK cells, which is most significant among different types of immune cells (Additional file 10: Fig. S6c), and the number of activated NK cells increased with the Gleason grade increase in PCa (Fig. 4b), indicating that activated NK cells might play a significant inhibitory function on prostate cancer stemness. The most significant positive correlation was between the activated dendritic cells and the memory B cells, and the number of both activated dendritic cells and memory B cells was all increased with the Gleason grade increase (Fig. 4c), which indicated that these two types of immune cells might play essential roles in inhibiting the stemness of PCa cells.
The protein interaction network of stemness-related genes and its relationship with prostate cancer immune cells
By screening and analyzing human protein interaction data containing MEmagenta module genes from the STRING database, we not only constructed the important protein–protein interaction network in the proteins of MEmagenta module genes, but also found that EZH2 interacted directly with 17 proteins. In the EZH2-related 17 protein–protein interactions, we also found that EZH2 could regulate the entire protein interaction network of MEmagenta module stemness genes by mainly interacting with CENPA, BUB1B, and PARP1 (Fig. 5a) [52,53,54,55]. Furthermore, function enrichment analysis of the MEmagenta module stemness genes revealed that most functions of these genes were concentrated in cell mitosis, and the most significant functional pathway was related to cell cycle (Additional file 12: Fig. S7a, b), suggesting that these stemness genes might involve in the regulation of PCa stem cell mitosis.
From the analysis results of the relationship between the expression of the MEmagenta module stemness genes and immune infiltration of PCa, we found that different immune cells had different effects on the expression of stemness genes in PCa. Different types of immune cells could affect the same stemness gene expression, and one stemness gene expression could also reversely affect different types of immune cells (Additional file 13: Table S6). In the correlations between the expression of stemness genes and immune cells, we found that expression of most stemness genes was positively correlated to memory B cells and naive B cells and negatively correlated to plasma cells (Fig. 5b). With the increase of PCa Gleason grade, the number of B cells increased and the number of plasma cells decreased. (Fig. 4b). These results indicated that B cells might play the opposite effects on PCa cell stemness in different conditions (it's consistent with reference ). Besides, the number of activated NK cells and memory CD4 + T cells were all increased with the PCa Gleason grade increase; and the expression of the MEmagenta module stemness genes was positively correlated to resting NK cells and memory CD4 + T cells, while negatively correlated to activated NK cells (Figs. 4b and 5b).
Expression of EZH2, the most relevant gene to PCa cell stemness, was most positively correlated to the activated memory CD4 + T cells and negatively correlated to the resting Mast cells (Fig. 5c, d). Most types of immune cells positively correlated to the expression of the EZH2 were T cells and B cells (Fig. 5c), suggesting that T cells and B cells were the critical immune cells in regulating PCa cell stemness by controlling the expression of the EZH2. In addition, with the increased PCa Gleason grade, the number of resting mast cells decreased and the number of activated mast cells increased (Fig. 4b). The resting mast cells were the most negatively correlated to the expression of the EZH2 (Fig. 5c, d). These results indicated that the immunity and immune cells of the microenvironment of PCa played an essential role in the tumorigenesis and development of prostate cancer.
It had been identified that cancer stem cells were closely related to cancer cell migration and played critical roles in the metastasis of PCa in clinic . The research on the stemness of PCa cells and their influencing factors was inevitably becoming a key direction of tumor research. In this study, we found the gene module (MEmagenta module) most related to stemness of prostate cancer cells by analyzing multi-omics of PCa clinical samples and cell lines, and also found that the CNVs and transcriptional regulators of the MEmagenta module genes were the critical influencing factors related to the stemness of PCa cells by conjoint analysis of transcriptome data, gene CNV data and ATAC-seq data of prostate cancer. In the analysis results of the immune microenvironment of prostate cancer, it's found that the immune microenvironment of PCa negatively correlated with the stemness of PCa cells, which played an important role in the stemness of PCa cells. The immune cells in the immune microenvironment of PCa had different correlations with the expression of stemness genes of the MEmagenta module, suggesting that the immune microenvironment of PCa and its immune cells were also the essential influencing factors related to the stemness of PCa cells.
It's reported that tumor cell stemness could affect the immune response of the tumor immune microenvironment and then result in tumor heterogeneity . Although there have been many studies on the relationship between tumor cell stemness and tumor development in other tumors, few studies on the relationship between tumor cell stemness and PCa development, especially on the factors affecting and regulating the cell stemness of prostate cancer. In this study, by fully using a variety of bioinformatics analysis software and methods to analyze multiple omics data, we not only obtained two factors related to PCa cell stemness, including the CNVs of stemness-related genes and the immunity of tumor microenvironment, but also got the gene set most relevant to stemness of PCa cells. Many critical genes included in this gene set and their transcription regulators were important in regulating and affecting PCa metastasis via influencing PCa cell stemness.
EZH2 not only promotes the formation of cancer stem cells, but also expands the aggressive cancer cell population and leads to cancer progression . EZH2 can also co-regulate prostate cancer stem cell properties with BRCA1 . RUBVL1 is a gene relevant to the stemness of prostate cancer cells, and its copy number increases with the tumor progresses. It is reported that RUVBL1 is essential for the survival of hematopoietic stem cells , and its gene copy number is increased in head and neck squamous cancers . MYBL2 can help DNA double-strand repair in hematopoietic stem cells  and the emergence of CNV leads to the occurrence and development of cancer [61, 62]. AURKB could determine the identity of embryonic stem cells , and the CNV (deletion) also contributes to the formation of aggressive tumors . PUM1 regulates the expression of hematopoietic stem cells  and promotes the migration of cancer cells . CLOCK regulates the biological clock of cancer stem cells and promotes the self-renewal of cancer cells . HMGCS1 can promote cancer development  and affect the function of NK cells . SUV39H1 attenuates the apoptosis of cancer cells  and enhances the immune escape of tumor cells . In addition to the reported genes related to cancer cell stemness, our results also showed many genes that had not been reported to be related to cancer cell stemness, especially to PCa cell stemness, such as TEDC2, TMEM132A, and VARS, etc.
Up to now, there are few studies on the relationship between cell stemness regulatory factors and tumor malignancy in PCa. From our multi-omics analysis results, the CNVs of cancer cell stemness-related genes were closely related to the malignant degree of prostate cancer (Gleason grade). The CNVs of genes positively related to PCa cell stemness were also positively correlated with the degree of malignant prostate cancer (Gleason grade); conversely, the CNVs of genes negatively associated with PCa cell stemness were also negatively correlated with Gleason grades of prostate cancers. The correlations between transcriptional regulators of PCa cell stemness genes and prostate tumor malignancy were similar to the correlation between CNVs of PCa cell stemness genes and prostate tumor malignancy.
Although we found that the CNVs and transcriptional regulations of PCa cell stemness genes and the immune infiltration of prostate cancer were all important factors in influencing the stemness of PCa cells and further regulating the development of PCa and we also analyzed and obtained the key gene set and genes in regulating the stemness of PCa cells, the experimental evidence and the detailed mechanism of the correlation and regulation among CNVs of PCa cell stemness genes, transcription regulators of PCa cell stemness genes and immunity of PCa microenvironment were still unclear and needed to be further investigated in the future.
Prostate cancer samples with high Gleason grade had high CNVs of stemness genes and high stemness score. Low immunity of immune microenvironment was beneficial to the stemness of prostate cancer cells. Stemness of PCa cells was positively correlated with the malignancy of PCa. The transcription factors of most of important stemness genes were important in promoting the stemness of PCa cells via enhancing the expression of stemness genes in prostate cancer. However, its detailed molecular mechanism still needed to be future identified with in vivo and in vitro experiments.
Availability of data and materials
All data generated or analyzed during this study are available.
Prostate cancer stem cells
Weighted gene co-expression network analysis
Functional protein association networks
Copy number variants
Cancer stem cells
Small cell prostate cancer
Friedmann-Morvinski D, Verma IM. Dedifferentiation and reprogramming: origins of cancer stem cells. EMBO Rep. 2014;15:244–53.
Malta TM, Sokolov A, Gentles AJ, Burzykowski T, Poisson L, Weinstein JN, Kaminska B, Huelsken J, Omberg L, Gevaert O, et al. Machine learning identifies stemness features associated with oncogenic dedifferentiation. Cell. 2018;173:338-354 e315.
Skvortsov S, Skvortsova I-I, Tang DG, Dubrovska A. Concise review: prostate cancer stem cells: current understanding. Stem Cells. 2018;36:1457–74.
Rycaj K, Tang DG. Cell-of-origin of cancer versus cancer stem cells: assays and interpretations. Cancer Res. 2015;75:4003.
Mei W, Lin X, Kapoor A, Gu Y, Zhao K, Tang D. The contributions of prostate cancer stem cells in prostate cancer initiation and metastasis. Cancers. 2019;11:434.
Liao C-P, Adisetiyo H, Liang M, Roy-Burman P. Cancer Stem cells and microenvironment in prostate cancer progression. Hormones Cancer. 2010;1:297–305.
Krueger TE, Thorek DLJ, Meeker AK, Isaacs JT, Brennen WN. Tumor-infiltrating mesenchymal stem cells: Drivers of the immunosuppressive tumor microenvironment in prostate cancer? Prostate. 2019;79:320–30.
Wang W, Green M, Choi JE, Gijón M, Kennedy PD, Johnson JK, Liao P, Lang X, Kryczek I, Sell A, et al. CD8+ T cells regulate tumour ferroptosis during cancer immunotherapy. Nature. 2019;569:270–4.
Mills CD, Lenz LL, Harris RA. A breakthrough: macrophage-directed cancer immunotherapy. Cancer Res. 2016;76:513.
Ihle CL, Owens P. Integrating the immune microenvironment of prostate cancer induced bone disease. Mol Carcinog. 2020;59:822–9.
Shiao SL, Chu GC, Chung LW. Regulation of prostate cancer progression by the tumor microenvironment. Cancer Lett. 2016;380:340–8.
Weng CC, Ding PY, Liu YH, Hawse JR, Subramaniam M, Wu CC, Lin YC, Chen CY, Hung WC, Cheng KH. Mutant Kras-induced upregulation of CD24 enhances prostate cancer stemness and bone metastasis. Oncogene. 2019;38:2005–19.
Pai VC, Hsu CC, Chan TS, Liao WY, Chuu CP, Chen WY, Li CR, Lin CY, Huang SP, Chen LT, Tsai KK. ASPM promotes prostate cancer stemness and progression by augmenting Wnt-Dvl-3-beta-catenin signaling. Oncogene. 2019;38:1340–53.
Xu N, Wu YP, Yin HB, Xue XY, Gou X. Molecular network-based identification of competing endogenous RNAs and mRNA signatures that predict survival in prostate cancer. J Transl Med. 2018;16:274.
Zhu X, Gou X, Zhou M. Nomograms predict survival advantages of gleason score 3+4 Over 4+3 for prostate cancer: a SEER-based study. Front Oncol. 2019;9:646.
Cancer Genome Atlas Research N. The molecular taxonomy of primary prostate cancer. Cell. 2015;163:1011–25.
He Z, Tang F, Lu Z, Huang Y, Lei H, Li Z, Zeng G. Analysis of differentially expressed genes, clinical value and biological pathways in prostate cancer. Am J Transl Res. 2018;10:1444–56.
He Z, Duan X, Zeng G. Identification of potential biomarkers and pivotal biological pathways for prostate cancer using bioinformatics analysis methods. PeerJ. 2019;7:e7872.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformat. 2008;9:559.
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37:907–15.
Anders S, Pyl PT, Huber W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.
Mermel CH, Schumacher SE, Hill B, Meyerson ML, Beroukhim R, Getz G. GISTIC20 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 2011;12:41.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Ramírez F, Dündar F, Diehl S, Grüning BA, Manke T. deepTools: a flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42:W187–91.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.
Phanstiel DH, Boyle AP, Araya CL, Snyder MP, Sushi R. flexible, quantitative and integrative genomic visualizations for publication-quality multi-panel figures. Bioinformatics. 2014;30:2808–10.
Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, Treviño V, Shen H, Laird PW, Levine DA, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12:453–7.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
Miranda A, Hamilton PT, Zhang AW, Pattnaik S, Becht E, Mezheyeuski A, Bruun J, Micke P, de Reynies A, Nelson BH. Cancer stemness, intratumoral heterogeneity, and immune response across cancers. Proc Natl Acad Sci. 2019;116:9020.
Gorodetska I, Lukiyanchuk V, Peitzsch C, Kozeretska I, Dubrovska A. BRCA1 and EZH2 cooperate in regulation of prostate cancer stem cell phenotype. Int J Cancer. 2019;145:2974–85.
Song I-S, Jeong YJ, Jeong SH, Heo HJ, Kim HK, Bae KB, Park Y-H, Kim SU, Kim J-M, Kim N, et al. FOXM1-Induced PRX3 regulates stemness and survival of colon cancer cells via maintenance of mitochondrial function. Gastroenterology. 2015;149:1006-1016.e1009.
Zhang Q, Huang H, Liu A, Li J, Liu C, Sun B, Chen L, Gao Y, Xu D, Su C. Cell division cycle 20 (CDC20) drives prostate cancer progression via stabilization of β-catenin in cancer stem-like cells. EBioMedicine. 2019;42:397–407.
Dai C, Miao CX, Xu XM, Liu LJ, Gu YF, Zhou D, Chen LS, Lin G, Lu GX. Transcriptional activation of human CDCA8 gene regulated by transcription factor NF-Y in embryonic stem cells and cancer cells. J Biol Chem. 2015;290:22423–34.
Naef V, Monticelli S, Corsinovi D, Mazzetto MT, Cellerino A, Ori M. The age-regulated zinc finger factor ZNF367 is a new modulator of neuroblast proliferation during embryonic neurogenesis. Sci Rep. 2018;8:11836.
Hu R, Wang MQ, Niu WB, Wang YJ, Liu YY, Liu LY, Wang M, Zhong J, You HY, Wu XH, et al. SKA3 promotes cell proliferation and migration in cervical cancer by activating the PI3K/Akt signaling pathway. Cancer Cell Int. 2018;18:183.
Taniuchi K, Furihata M, Iwasaki S, Tanaka K, Shimizu T, Saito M, Saibara T. RUVBL1 directly binds actin filaments and induces formation of cell protrusions to promote pancreatic cancer cell invasion. Int J Oncol. 2014;44:1945–54.
Bereshchenko O, Mancini E, Luciani L, Gambardella A, Riccardi C, Nerlov C. Pontin is essential for murine hematopoietic stem cell survival. Haematologica. 2012;97:1291–4.
Liu M, Hu Q, Tu M, Wang X, Yang Z, Yang G, Luo R. MCM6 promotes metastasis of hepatocellular carcinoma via MEK/ERK pathway and serves as a novel serum biomarker for early recurrence. J Exp Clin Cancer Res. 2018;37:10.
Barton KM, Levine EM. Expression patterns and cell cycle profiles of PCNA, MCM6, cyclin D1, cyclin A2, cyclin B1, and phosphorylated histone H3 in the developing mouse retina. Dev Dyn. 2008;237:672–82.
Tomonaga T, Matsushita K, Ishibashi M, Nezu M, Shimada H, Ochiai T, Yoda K, Nomura F. Centromere protein H is up-regulated in primary human colorectal cancer and its overexpression induces aneuploidy. Cancer Res. 2005;65:4683–9.
Zhang JP, Zhang H, Wang HB, Li YX, Liu GH, Xing S, Li MZ, Zeng MS. Down-regulation of Sp1 suppresses cell proliferation, clonogenicity and the expressions of stem cell markers in nasopharyngeal carcinoma. J Transl Med. 2014;12:222.
Zhao WF, Wang HB, Xie B, Hu LJ, Xu LH, Kuang BH, Li MZ, Zhang X. Sp1 and Sp3 are involved in the full transcriptional activity of centromere protein H in human nasopharyngeal carcinoma cells. FEBS J. 2012;279:2714–26.
Uyhazi KE, Yang Y, Liu N, Qi H, Huang XA, Mak W, Weatherbee SD, Song X, Lin H. Pumilio proteins exert distinct biological functions and multiple modes of post-transcriptional regulation in embryonic stem cell pluripotency and early embryogenesis. BioRxiv. 2019. https://doi.org/10.1101/751909.
Guan X, Chen S, Liu Y. Wang L-l, Zhao Y, Zong Z-H: PUM1 promotes ovarian cancer proliferation, migration and invasion. Biochem Biophys Res Commun. 2018;497:313–8.
Janich P, Pascual G, Merlos-Suárez A, Batlle E, Ripperger J, Albrecht U. Cheng H-YM, Obrietan K, Di Croce L, Benitah SA: The circadian molecular clock creates epidermal stem cell heterogeneity. Nature. 2011;480:209–14.
Gambara G, Desideri M, Stoppacciaro A, Padula F, De Cesaris P, Starace D, Tubaro A, del Bufalo D, Filippini A, Ziparo E, Riccioli A. TLR3 engagement induces IRF-3-dependent apoptosis in androgen-sensitive prostate cancer cells and inhibits tumour growth in vivo. J Cell Mol Med. 2015;19:327–39.
Blum R, Gupta R, Burger PE, Ontiveros CS, Salm SN, Xiong X, Kamb A, Wesche H, Marshall L, Cutler G, et al. Molecular signatures of prostate stem cells reveal novel signaling pathways and provide insights into prostate cancer. PLoS ONE. 2009;4:e5722.
Wouters MCA, Nelson BH. Prognostic significance of tumor-infiltrating b cells and plasma cells in human cancer. Clin Cancer Res. 2018;24:6125.
Arco A, Edgar BA, Erhardt S. In vivo analysis of centromeric proteins reveals a stem cell-specific asymmetry and an essential role in differentiated, non-proliferating cells. Cell Rep. 2018;22:1982–93.
Behnan J, Grieg Z, Joel M, Ramsness I, Stangeland B. Gene knockdown of CENPA reduces sphere forming ability and stemness of glioblastoma initiating cells. Neuroepigenetics. 2016;7:6–18.
Venere M, Miller TE, Rich JN. Mitotic control of cancer stem cells. Cancer Discov. 2013;3:141.
Hsieh M-H, Chen Y-T, Chen Y-T, Lee Y-H, Lu J, Chien C-L, Chen H-F, Ho H-N, Yu C-J, Wang Z-Q, Teng S-C. PARP1 controls KLF4-mediated telomerase expression in stem cells and cancer cells. Nucleic Acids Res. 2017;45:10492–503.
Sharonov GV, Serebrovskaya EO, Yuzhakova DV, Britanova OV, Chudakov DM. B cells, plasma cells and antibody repertoires in the tumour microenvironment. Nat Rev Immunol. 2020;20:294–307.
Horning AM, Wang Y, Lin CK, Louie AD, Jadhav RR, Hung CN, Wang CM, Lin CL, Kirma NB, Liss MA, et al. Single-Cell RNA-seq reveals a subpopulation of prostate cancer cells with enhanced cell-cycle-related transcription and attenuated androgen response. Cancer Res. 2018;78:853–64.
Chang CJ, Hung MC. The role of EZH2 in tumour progression. Br J Cancer. 2012;106:243–7.
Lin D, Lin B, Bhanot H, Riou R, Abt NB, Rajagopal J, Saladi SV. RUVBL1 is an amplified epigenetic factor promoting proliferation and inhibiting differentiation program in head and neck squamous cancers. Oral Oncol. 2020;111:104930.
Bayley R, Blakemore D, Cancian L, Dumon S, Volpe G, Ward C, Almaghrabi R, Gujar J, Reeve N, Raghavan M, et al. MYBL2 supports DNA double strand break repair in hematopoietic stem cells. Cancer Res. 2018;78:5767–79.
Musa J, Aynaud MM, Mirabeau O, Delattre O, Grunewald TG. MYBL2 (B-Myb): a central regulator of cell proliferation, cell survival and differentiation involved in tumorigenesis. Cell Death Dis. 2017;8:e2895.
Nakajima T, Yasui K, Zen K, Inagaki Y, Fujii H, Minami M, Tanaka S, Taniwaki M, Itoh Y, Arii S, et al. Activation of B-Myb by E2F1 in hepatocellular carcinoma. Hepatol Res. 2008;38:886–95.
Shin J, Kim TW, Kim H, Kim HJ, Suh MY, Lee S, Lee HT, Kwak S, Lee SE, Lee JH, et al. Aurkb/PP1-mediated resetting of Oct4 during the cell cycle determines the identity of embryonic stem cells. Elife. 2016;5:e10877.
Hegyi K, Egervari K, Sandor Z, Mehes G. Aurora kinase B expression in breast carcinoma: cell kinetic and genetic aspects. Pathobiology. 2012;79:314–22.
Naudin C, Hattabi A, Michelet F, Miri-Nezhad A, Benyoucef A, Pflumio F, Guillonneau F, Fichelson S, Vigon I, Dusanter-Fourt I, Lauret E. PUMILIO/FOXP1 signaling drives expansion of hematopoietic stem/progenitor and leukemia cells. Blood. 2017;129:2493–506.
Dong Z, Zhang G, Qu M, Gimple RC, Wu Q, Qiu Z, Prager BC, Wang X, Kim LJY, Morton AR, et al. Targeting glioblastoma stem cells through disruption of the circadian clock. Cancer Discov. 2019;9:1556–73.
Ashida S, Kawada C, Inoue K. Stromal regulation of prostate cancer cell growth by mevalonate pathway enzymes HMGCS1 and HMGCR. Oncol Lett. 2017;14:6533–42.
Assmann N, O’Brien KL, Donnelly RP, Dyck L, Zaiatz-Bittencourt V, Loftus RM, Heinrich P, Oefner PJ, Lynch L, Gardiner CM, et al. Srebp-controlled glucose metabolism is essential for NK cell functional responses. Nat Immunol. 2017;18:1197–206.
Watson GW, Wickramasekara S, Palomera-Sanchez Z, Black C, Maier CS, Williams DE, Dashwood RH, Ho E. SUV39H1/H3K9me3 attenuates sulforaphane-induced apoptotic signaling in PC3 prostate cancer cells. Oncogenesis. 2014;3:e131.
Lu C, Yang D, Klement JD, Oh IK, Savage NM, Waller JL, Colby AH, Grinstaff MW, Oberlies NH, Pearce CJ, et al. SUV39H1 represses the expression of cytotoxic T-Lymphocyte effector genes to promote colon tumor immune evasion. Cancer Immunol Res. 2019;7:414–27.
Thanks to Dr. Wei Guanyun for her suggestions on data analysis.
This research was funded by the National Natural Science Foundation of China (NNSFC) (Grant No. 81872104 and 81472415).
Ethics approval and consent to participate
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Figure S1. Screen and verification of stemness gene modules by WGCNA analysis. a WGCNA analysis results. Branches with different colors corresponding to different modules. b Analysis of the correlation among genes, modules and stemness. c The stemness gene obtained from jointly analyzing WGCNA analysis results and the most known stemness gene markers. d and e The expression of stemness genes in prostate cancer tissues and cell lines, respectively.
Table S1. The relationship between genes and color modules in WGCNA analysis results.
Table S2. Results of differential expression of stemness genes in MEmagenta module.
: Figure S2. CNVs of MEmagenta module stemness genes in normal prostate tissue samples. a Overview of CNV in normal prostate tissue samples. b and c The CNVs with amplification or with deletion, respectively.
Table S3. The result of copy number amplification in the genome location of the stemness genes in the MEmagenta module.
Table S4. The result of copy number deletion in the genome location of the stemness genes in the MEmagenta module.
: Figure S3 CNVs and expression of stemness genes all increased with the increase of Gleason grades. a and b The expression and CNVs (amplification) of SKA3 and RUVBL1 increased with increase of Gleason grades. c and d The expression and CNVs (deletion) of MCM6 and CENPH increased with increase of Gleason grades.
: Figure S4 Transcriptional regulators of stemness genes in prostate cancer. a Expression heatmap of transcriptional regulators of stemness genes in PCa samples with different Gleason grades. b Venn diagram of transcriptional regulators of stemness genes in in PCa samples with different Gleason grades. c The importance of transcriptional regulators of stemness genes in prostate cancer was obtained by using PCA analysis.
: Figure S5 Analysis of transcriptional regulators of stemness genes in normal samples. a The main transcriptional regulators were obtained by using PCA analysis. b Based on PCA analysis, the importance of known transcriptional regulators of stemness genes in normal samples was obtained. c Based on PCA analysis, the importance of transcriptional regulators of stemness genes in normal samples was obtained.
: Figure S6 Prostate tumor immune microenvironment is negatively correlated with stemness. a Stromal cell score and immunity score of the PCa microenvironment were negatively correlated with PCa cell stemness score. b The immune score increased with the Gleason grade increase, and the low immunity score had a higher stemness score in PCa with different Gleason grades. c The correlations between different types of immune cells in 22 types of immune cell. d The expression heatmap of stemness genes in PCa samples with high/low immune score, different stemness score, and different Gleason grade.
Table S5. Results of correlation between different immune cells.
: Figure S7 GO and KEGG enrichment analysis of stemness gene in MEmagenta module. a The results of GO enrichment analysis. b The results of KEGG enrichment analysis.
Table S6. Correlation results between stemness genes and different immune cells in the MEmagenta module.
About this article
Cite this article
Dai, Z., Liu, P. High copy number variations, particular transcription factors, and low immunity contribute to the stemness of prostate cancer cells. J Transl Med 19, 206 (2021). https://doi.org/10.1186/s12967-021-02870-x
- Stemness of prostate cancer
- Immune infiltration