ORIGINAL RESEARCH article
Development and Validation of Novel Biomarkers Related to M2 Macrophages Infiltration by Weighted Gene Co-Expression Network Analysis in Prostate Cancer
- 1Department of Urology, Urology Research Institute, The First Affiliated Hospital of Fujian Medical University, Fuzhou, China
- 2Fujian Key Laboratory of Precision Medicine for Cancer, The First Affiliated Hospital, Fujian Medical University, Fuzhou, China
- 3Central Lab, Key Laboratory of Radiation Biology of Fujian Higher Education Institutions, The First Affiliated Hospital, Fujian Medical University, Fuzhou, China
M2-tumor-associated macrophages (TAMs) work as a promoter in the processes of bone metastases, chemotherapy resistance, and castration resistance in prostate cancer (PCa), but how M2-TAMs affect PCa has not been fully understood. In this study, we analyzed the proportion of tumor-infiltrating immune cells using the CIBERSORT algorithm, based on samples from the Cancer Genome Atlas database. Then we performed weighted gene co-expression network analysis to examine the modules concerning infiltrated M2-TAMs. Gene Ontology analysis and pathway enrichment analysis were performed for functional annotation and a protein–protein interaction network was constructed. The International Cancer Genomics Consortium cohort was used as a validation cohort. The red module showed the most correlation with M2-TAMs in PCa. Biological processes and pathways were mainly associated with the immune-related processes, as revealed by functional annotation. Four hub genes were screened: ACSL1, DLGAP5, KIF23 and NCAPG. Further validation showed that the four hub genes had a higher expression level in tumor tissues than that in normal tissues, and they were good prognosis biomarkers for PCa. In conclusion, these findings contribute to understanding the underlying molecular mechanisms of how M2-TAMs affect PCa, and looking for the potential biomarkers and therapeutic targets for PCa patients.
Prostate cancer (PCa) is a heterogeneous disease, ranging from an asymptomatic stage to an distantly metastatic or castration-resistant stage (1, 2). Guideline recommends radical prostatectomy or radical radiotherapy as the primary treatment for most patients with clinical localized PCa, and these patients usually have a good prognosis (3). As for patients with metastatic or castration-resistant disease, however, local treatments have limited efficacy for them. They eventually die when systemic treatment fails to control the progression of the disease (4). Although the number of patients with advanced disease has been substantially reduced and various novel therapeutic agents have presented notable success for patients with metastatic and castration-resistant PCa in recent years, some of them still experience disease progression, which may be attributed to individual differences at the genetic level (5–7). Therefore, it is urgent to elucidate the molecular mechanism of the disease and look for appropriate biomarkers, which will be beneficial to the diagnosis, treatment and prognosis prediction for PCa.
In recent years, it has been recognized that the tumor microenvironment (TME) is intimately involved in tumor development and progression (8–10). Tumor-associated macrophages (TAMs) are the main immune cells in TME. Generally, the macrophages are divided into two subsets depending on different polarized status: the classically activated (M1) and the alternatively activated (M2) macrophages (11). In inflammatory condition, classically activated M1-TAMs phenotype prevails in sites of inflammation; however, tumor development gradually promotes a phenotypic switch in which TAMs represent a immunosuppressive M2-TAMs phenotype (12). The presence of M2-TAMs has been related to poor clinical prognosis in several malignant diseases and is thought to affect disease outcome by stimulating angiogenesis, promoting distant metastasis, suppressing antitumor immunity and possibly by reducing the effectiveness of certain treatments (9). It was reported that M2-TAMs worked as a promoter in the processes of bone metastases, chemotherapy resistance, and castration resistance in PCa (13–15). However, the specific interactive mechanism between M2-TAMs and PCa has not been fully clarified. It is crucial to explored the key genes and molecules in this interaction process.
Weighted gene co-expression network analysis (WGCNA) is an algorithm that construct free-scale gene co-expression networks to explore the relationships between gene sets and clinical features (16). It has been widely used to screen the hub genes associated with clinical feature in different cancer types (17). In present study, we performed WGCNA to explore the role of TMAs and screen the potential biomarkers based on PCa gene expression data. We calculated the M2-TAMs proportion in these samples by CIBERSORT algorithm and then identified important modules and hub genes associated with the proportion of infiltrated M2-TAMs. This is the first utilization of WGCNA to identify M2-TAMs-related hub genes and biomarkers of PCa.
Materials and Methods
Data Collection and Preprocessing
A total of 499 cases of PCa samples from the Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/), all of their expression profiles were downloaded. Cases were excluded due to inadequate clinical data. CIBERSORT algorithm was used to evaluate the cellular composition of immune cells. Meanwhile, clinical significance of each type of immune cells was evaluated using survival analyses and correlation analyses with clinical features. Thereafter, we screened differently expressed genes (DEGs) of samples with high and low ratio of M2-TAMs.WGCNA was performed to determine the module associated with M2-TAMs. Gene Ontology (GO) analysis and pathway enrichment analysis were performed for functional annotation of selected modules. A protein–protein interaction (PPI) network was built and hub genes were screened according to the degree of connectivity. Then, online databases were used for further validation. Besides, the additional independent cohort including 494 PCa samples from International Cancer Genomics Consortium (ICGC) database (https://dcc.icgc.org/) and cohort from TCGA database were used for survival analyses. Figure 1 shows the flowchart detailing the study design and samples.
Evaluation of Tumor-Infiltrating Immune Cells
CIBERSORT is a deconvolution algorithm that uses a set of reference gene expression values (a “signature matrix” of 547 genes) considered a minimal representation for each cell type and, based on those values, infers cell type proportions in data from bulk tumor samples of mixed cell types using support vector regression (18). Normalized gene expression data were used to infer the relative proportions of 22 types of infiltrating immune cells using the CIBERSORT algorithm. Briefly, gene expression datasets were prepared using standard annotation files and data uploaded to the CIBERSORT web portal (http://cibersort.stanford.edu/), with the algorithm run using the default signature matrix at 1,000 permutations. CIBERSORT derives a P-value for the deconvolution of each sample using Monte Carlo sampling, providing a measure of confidence in the results. Finally, we selected 292 samples with a CIBERSORT P-value of < 0.05, from the analyzed samples of 499 PCa cases.
Differentially Expressed Genes Screening
The DEGs between samples with high and low proportions of M2-TAMs in TCGA cohort were analyzed using the “limma” R package. As described previously, adjusted P-value <0.05 and |logFC| ≥2 was set as the cut-off criterion for better accuracy and significance (19, 20).
Construction of Co-Expression Network
Data quality assessment and pre-processing were performed at first. Then, WGCNA algorithm was used to construct a scale-free co-expression network for the DEGs. After that, a weighted adjacency matrix was established, as calculated by amn = |cmn|β (Where cmn is Pearson’s correlation between gene m and gene n, and where amn is adjacency between gene m and gene n). Considered as a soft-thresholding parameter, β can emphasize strong relations between genes while penalize weak correlations. We chose a proper power of β according to the mean connectivity. We transformed the adjacency into a topological overlap matrix (TOM), which is used to describe corresponding dissimilarity (1-TOM) (21). Subsequently, according to the TOM-based dissimilarity with a minimum size (gene group) of 30 for the genes dendrogram, we constructed average linkage hierarchical clustering aimed to classify genes that have similar expression profile into modules.
Identification of Modules Related to the Proportion of M2 Macrophage Cells
Two methods were used for identification of modules associated with clinical features of PCa. Gene significance (GS) defined as log10 transformation of the P-value (GS = lgP) in the linear regression between gene expression and the clinical characteristic and module significance (MS) considered as the average GS of all genes in the module were calculated, respectively. The module with the MS ranked first in all of the selected modules was regarded as the one with clinical significance. Module eigengene (ME) was defined as the first principle component in each gene modules to summarize the expression patterns of all genes into a single characteristic expression profile in a given module. The relation between each ME and clinical characteristic for identification of the relevant module was computed. Ultimately, we chose the module highly correlated to specific clinical features for further analysis.
Functional Enrichment Analysis
Metascape (http://metascape.org/) and Webgestalt (http://www.webgestalt.org/) are online databases providing a comprehensive set of functional annotation tools for researchers to better understand biological meaning behind large list of genes (22). We uploaded genes to perform GO analysis and pathway enrichment analysis. P-value <0.05 was thought statistically significant.
PPI Network and Hub Genes Selection
Search Tool for the Retrieval of Interacting Genes (STRING) is a biological database for constructing PPI networks, providing a system-wide view of interactions between each member (23). Genes of selected module were mapped to STRING to explore their relationships with each other. A combined score of >0.4 was set as the cut-off criterion as described previously (24). Then, we established PPI network using Cytoscape software (24). The genes with high degree of connectivity were selected as the hub genes.
Validation of the Hub Genes by Online Database
UALCAN (http://ualcan.path.uab.edu/) is a portal for facilitating tumor subgroup gene expression and survival analyses (25). Expression levels in mRNA and promoter methylation levels of hub genes were revealed using UALCAN. Afterwards, we used the Human Protein Atlas (http://www.proteinatlas.org) for validation in immunohistochemistry aspect.
Evaluation the Value of Hub Genes as Prognosis Biomarkers
Cohort containing 494 cases of PCa samples from ICGC database and Cohort containing 292 cases from TCGA database were used for survival analyses by Kaplan–Meier method. Patients were divided into high expression group and low expression group according to a cut-off value of mean expression of the hub genes. Afterwards, survival analyses for hub genes were performed. The hazard ratio (HR) with 95% confidence intervals and log-rank P-value was calculated and displayed.
For statistical analysis, we used GraphPad Prism 5.0 (GraphPad Software, San Diego, CA, USA) and SPSS version 22.0 software (SPSS, Chicago, IL, USA). The relationship between the proportion of M2-TAMs and clinicopathological parameters were analyzed using chi-square test. Survival curves plotted by the Kaplan–Meier method were compared to the log-rank test. P < 0.05 was considered statistically significant.
Tumor-Infiltrating M2 Macrophage in PCa
The CIBERSORT algorithm was used to investigate the 22 infiltrating immune cells subsets based on 292 PCa samples, whose clinical characteristics were shown in Figure 2A. Age, T stage, N stage, and survival status of the cases were collected. The results revealed that a higher proportion of infiltrating M2-TAMs was related to a more advanced stage (Figure 2C). Compared to patients with a low proportion of M2-TAMs, those with a high proportion of M2-TAMs had a worse prognosis (Figure 2B).
Figure 2 CIBERSORT analysis and clinical significance of M2-TAMs in PCa. (A) Relative percent of each type of immune cell in 292 PCa samples from TCGA cohort; (B) Proportion of M2-TAMs in different stages; (C) Overall survival between patients with high and low proportion of infiltrated M2-TAMs. *P < 0.05; **P < 0.01.
Using the CIBERSORT algorithm, we detected the proportion of infiltrating M2-TAMs and the mRNA expression data from 292 PCa samples whose clinicopathology characteristics were summarized in Table 1. Using the average value of the proportion of M2-TAMs as a cut-off value, cases were divided into two groups, which one was 117 cases with a high proportion of M2-TAMs, another one was 175 cases with a low proportion of M2-TAMs. A total of 2,950 DEGs (1,745 upregulated and 1,205 downregulated) were chosen for subsequent analysis under the threshold of adjusted P-value <0.05 and |logFC| ≥2.
Construction of Weighted Co-Expression Network and Identification of Key Module
The power of β = 22 (scale free R2 = 0.85) was selected as the best soft-thresholding parameter (Figures 3A, B). As can be seen from Figures 3C, D, the rationality test had a positive result. Then, the dendrogram of sample constructed shows the similarity among the samples, in which the clinical characteristics of each sample are displayed (Figure 3E). After that, nine modules were identified (Figure 3F). Two methods were used to test the relationship between each module and the M2-TAMs. Modules with a higher MS value were considered to have more connection with high infiltration of M2-TAMs, and we found that the MS of the red module was higher than those of any other modules (Figure 4A). Afterwards, the ME of the red module revealed a higher correlation with M2-TAMs than other modules (Figure 4C). Eventually, we identified the red module as the module most associated with high infiltration of M2-TAMs in PCa (Figure 4B).
Figure 3 WGCNA was performed based on PCa samples from TCGA database. (A) Analysis of the scale-free fit index for various soft-thresholding powers; (B) Analysis of the mean connectivity for various soft-thresholding powers; (C) Histogram of connectivity distribution when β = 22; (D) Checking the scale free topology when β = 22; (E) The sample dendrogram and corresponding clinical characteristics; red: high, blue: low. (F) Cluster dendrogram of 292 samples with eligible data.
Figure 4 Identification of modules associated with clinical characteristics. (A) Distribution of average gene significance and errors in the modules associated with the proportion of M2-TAMs in PCa; (B) Scatter plot of module eigengenes in red module; (C) Heatmap of the correlation between module eigengenes and different clinical characteristics of PCa.
Functional Enrichment Analysis
Functional enrichment analysis was conducted to look for the biological processes and pathways associated with red module. We used the tool of Webgestalt to perform GO analysis. GO analysis of biological process showed that genes in red module were mainly involved in response to stimulus, biological regulation, cell communication, localization and metabolic process (Figure 5A). GO analysis of cellular component showed that these genes were mainly enriched in membrane, vesicle, endomembrance system, extracellular space and protein-containing complex (Figure 5B). GO analysis of molecular function demonstrated that these genes were mainly involved in protein binding, molecular transducer activity, ion binding, carbohydrate binding and transferase activity (Figure 5C). Then we used Metascape to further investigate the relevant biological processes and pathways. The result revealed that the biological processes and pathways were mainly associated with the immune processes, such as myeloid leukocyte activation, osteoclast differentiation, negative regulation of immune system process, immune response-regulating signaling pathway and regulation of cytokine production (Figures 5D, E and Table 2).
Figure 5 Functional enrichment analysis and construction of PPI network. (A) GO analysis of biological process for genes in red module; (B) GO analysis of cellular component for genes in red module; (C) GO analysis of molecular function for genes in red module; (D) Functional enrichment analysis for genes in red module. (E) P-value of each gene in the network; (F) PPI network constructed using STRING.
Construction of PPI Network and Identification of Hub Genes
PPI network was constructed with STRING (Figure 5F), and finally four hub genes were screened for further investigation. They were ACSL1 (Acyl-CoA Synthetase Long Chain Family Member 1), DLGAP5 (DLG Associated Protein 5), KIF23 (Kinesin Family Member 23) and NCAPG (Non-SMC Condensin I Complex Subunit G). The four genes were all of significantly clinical significance. The module membership values of them were 0.91, 0.82, 0.74 and 0.87, respectively, which were significant statistically. Figure 4B showed the scatter plot of module eigengenes in red module.
Validation and Efficacy Evaluation of Hub Genes
Data set from the online database UALCAN was used for validation. All four hub genes revealed higher expression levels in PCa samples, compared with normal tissues (Figures 6A, D, G, J). Besides, these hub genes revealed lower levels of promoter methylation in PCa tissues, when compared with the normal tissues (Figures 6B, E, H, K). Immunohistochemistry staining from the Human Protein Atlas database demonstrated that protein levels of the hub genes were significantly higher in PCa tissues than those in normal tissues (Figures 6C, F, I, L).
Figure 6 Validation of the hub genes using online databases. (A, D, G, J) Expression levels of the hub genes in PCa and normal prostate tissues; (B, E, H, K) Promoter methylation levels of the hub genes in PCa and normal prostate tissues; (C, F, I, L) Immunohistochemical staining of the hub genes in PCa and normal prostate tissues. *P < 0.05.
Data sets from the TCGA database and ICGC database (Table 3) were used for survival analyses to evaluate value of hub genes as prognosis biomarkers. It revealed that high expression levels of ACSL1 (HR 0.245 [0.092–0.427], P = 0.045) were related to poor overall survival in PCa patients, as well as DLGAP5 (HR 0.279 [0.108–0.494], P = 0.026), KIF23 (HR 0.305 [0.103–0.621], P = 0.011), and NCAPG (HR 0.271 [0.102–0.485], P = 0.031), based on TCGA cohort (Figures 7A, C, E, G). In addition, survival analyses based on ICGC cohort show consistent results: high expression levels of ACSL1 (HR 0.248 [0.101–0.418], P = 0.080), DLGAP5 (HR 0.416 [0.175–0.768], P = 0.039), KIF23 (HR 0.428 [0.117–0.696], P <0.001), and NCAPG (HR 0.231 [0.082–0.494], P = 0.005) were related to poor overall survival of PCa (Figures 7B, D, F, H).
Figure 7 Survival analyses based on TCGA and ICGC cohorts to evaluate value of hub genes as prognosis biomarkers. (A, C, E, G) Overall survival between patients with high and low expression of the four hub genes based on TCGA cohorts; (B, D, F, H) Overall survival between patients with high and low expression of the four hub genes based on ICGC cohorts.
Tumor cells and the microenvironment of the local host tissue interact with each other and form the TME. The TME plays a crucial role in tumor differentiation, tumor epigenetics, infiltration metastasis, and immune escape (26). Most therapies kill the tumor cells through cytocidal or cytostatic activity. In recent years, besides targeting tumor cells, therapies targeting the TME have been explored. For example, it has been demonstrated that targeting the vasculature can reduce tumor-related immunosuppression and thus improve therapeutic effects (27). Sipuleucel-T, an antigen-specific active immunotherapy agent, can sensitize the patient’s immune system against specific antigens or molecules, to train the system to attack the tumor, and has been proved to be beneficial for patients with metastatic castration-resistant prostate cancer (mCRPC) (28). As a major immune component in TME, TAMs have a tumor-suppressing polarized status called M1 subtype and a tumor-promoting polarized status called M2 subtype. The presence of M2-TAMs was related to poor clinical outcome in many cancers, such as squamous cell carcinoma, breast cancer and renal cell carcinoma (29–31). They cannot only promote the progression of the tumors, but also affect the efficacy of anti-tumor therapy (32, 33). In PCa, patients with higher proportions of M2-TAMs in TME had a higher death rate and a shorter time to biochemical recurrence (34). Recent study suggested that M2-TAMs cocultured with PCa cells caused epithelial–mesenchymal transition and proliferation of tumor cells (35). Besides, the presence of M2-TAMs may be associated with the therapy resistance to immunotherapy for patients with mCRPC (15). M2-TAMs can secrete cytokines and pro-metastatic factors which prevent immune responses initiated by cytotoxic T cells, causing immunological silence (36). In addition, Guan reported that increased infiltration of M2-TAMs induced a drug resistance to docetaxel treatment for CRPC cells, and PLX3397 could restore the chemosensitivity by preventing the recruitment and M2-polarization of TAMs (14). However, the specific interactive mechanism between M2-TAMs and PCa is still incompletely understood.
Recent technological advances have allowed the development of various molecular prognostic indicators (37, 38). However, these tools are not necessarily universal, given that they are restricted to subsets of patients based on criteria such as hormone receptor, pathological types and nodal status (39, 40). It is urgent to explore better biomarkers or predictive tools for PCa. Based on WCGNA analysis, the modules and genes associated with the proportion of infiltrated M2-TAMs in PCa were explored in the present study. The results revealed that red module was the most relevant one. Functional annotation showed that biological processes and pathways were mainly associated with the immune-related processes, which was plausible and supported the result of WGCNA. Consequently, we explored the hub genes in red module. In the end, we identified ACSL1, DLGAP5, KIF23, and NCAPG as the hub genes with the degree of connectivity. Subsequent validation demonstrated that these four hub genes had a higher expression level in tumor tissues than that in normal tissues, and they were good prognosis biomarkers related to the proportion of infiltrated M2-TAMs for PCa.
High expression of ACSL1 was significantly related to poor outcome a variety of tumors. In hepatocellular carcinoma, ACSL1 promotes the fat metabolism of the tumor and the progression of the disease by catalyzing the ATP-dependent acylation of fatty acids into long-chain acyl CoAs (LCA-CoAs) (41). It was reported that decreased ACSL1 could suppress the protumorigenic inflammatory responses induced by granulocyte–macrophage colony-stimulating factor (GM-CSF) in breast cancer (42). In PCa, lipid metabolism regulated by ACSL1 rewires the PCa metabolome to support growth and resistance to endocrine therapies. Therapeutic strategies targeting lipid metabolism and androgen receptor are starting to emerge, providing new chances to re-sensitize tumors to endocrine therapies with lipid metabolic approaches (43). However, research regarding the effect of M2-TAMs on lipid metabolism of PCa is relatively limited. ACSL1 may be a potential target, when more studies confirm its values.
DLGAP5, which is mapped to chromosome 14q22.3, is a cell cycle regulator involved in carcinogenesis. DLGAP5 facilitates microtubules formation by working as a kinetochore protein (44). The expression of DLGAP5 is positively associated with the recurrence rate of PCa. Besides, it was observed that DLGAP5 was highly expressed in CRPC samples in many studies (45). Furthermore, Hewit’s study on PCa suggested that DLGAP5 played an important role in surviving microtubule assault from docetaxel and stabilizing spindle formation, in an androgen-regulated cell cycle system (46). Therefore, further studies on DLGAP5 were needed and were of great significance in clarifying the mechanism of how M2-TAMs cause chemoresistance and hormone resistance of PCa.
As a member of the kinesin superfamily, KIF23 was highly expressed in many kinds of tumors, such as gliomas, ovarian cancer, and gastric cancer (47, 48). In PCa, it has been proved that several members of KIF family play important roles in the process of metastasis, invasion and drug resistance (49, 50). In present study, KIF23 showed a high correlation with the proportion of infiltrated M2-TAMs, which meant that KIF23 may be a key gene in the interaction between M2-TAMs and PCa.
NCAPG is a subunit of the condensin complex, which is responsible for the stabilization and condensation of chromosomes during mitosis and meiosis (51). Studies suggested that NCAPG might be an essential oncogene of hepatocellular carcinoma and gliomas by inducing apoptosis (52). Similarly, Arai reported that abnormal expression of NCAPG promoted PCa cell aggressiveness (53). High expression level of NCAPG was observed in clinical samples of CRPC, and its expression was found to be important for PCa pathogenesis, as revealed by analysis of TCGA database. Our study suggested NCAPG might be a potential target, through which M2-TAMs affect the biological behavior of PCa cells.
There were limitations in our study that need to be recognized. First, the clinical information from the online database was not comprehensive, and the detailed information of some clinical treatment regimens also could not be obtained. Therefore, we could only evaluate the effect of these hub genes on overall survival to a certain extent, rather than that on other clinical endpoints. Second, we screened these hub genes by bioinformatics methods, where further experiments were required to verify and to determine the mechanisms underlying the process of malignant progression of PCa. Furthermore, the present study still did not identify a better mathematical model to combine all eligible hub genes together for predicting the survival. Therefore, future studies may aim for a better re-evaluation of the prognostic performance of the model for PCa.
In conclusion, we identified four hub genes which were closely related to the proportion of infiltrated M2-TAMs in PCa. These findings contribute to understanding the underlying molecular mechanisms of how M2-TAMs affect PCa, and looking for the potential biomarkers and therapeutic targets for PCa patients. However, the molecular mechanism and function of these genes need to be confirmed in further experiments.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
This study was supported by Startup Fund for scientific research, Fujian Medical University (Grant number: 2018QH1067), the Young and Middle-aged Talents Training Project of Fujian Provincial Health Commission (Grant number: 2019-ZQN-54), and Educational Research Programs for Young and Middle-aged Teachers of Education Department of Fujian Province (Grant number: JAT190204).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
1. Lin TT, Chen YH, Wu YP, Chen SZ, Li XD, Lin YZ, et al. Risk Factors for Progression to Castration-Resistant Prostate Cancer in Metastatic Prostate Cancer Patients. J Cancer (2019) 10(22):5608–13. doi: 10.7150/jca.30731
2. Xu N, Wu YP, Chen DN, Ke ZB, Cai H, Wei Y, et al. Can Prostate Imaging Reporting and Data System Version 2 Reduce Unnecessary Prostate Biopsies in Men With PSA Levels of 4-10 Ng/Ml? J Cancer Res Clin Oncol (2018) 144(5):987–95. doi: 10.1007/s00432-018-2616-6
3. Mottet N, van den Bergh RCN, Briers E, Van den Broeck T, Cumberbatch MG, De Santis M, et al. EAU-EANM-ESTRO-ESUR-SIOG Guidelines on Prostate Cancer-2020 Update. Part 1: Screening, Diagnosis, and Local Treatment With Curative Intent. Eur Urol (2021) 79(2):243–62. doi: 10.1016/j.eururo.2020.09.042
4. Cornford P, van den Bergh RCN, Briers E, Van den Broeck T, Cumberbatch MG, De Santis M, et al. EAU-EANM-ESTRO-ESUR-SIOG Guidelines on Prostate Cancer. Part II-2020 Update: Treatment of Relapsing and Metastatic Prostate Cancer. Eur Urol (2021) 79(2):263–82. doi: 10.1016/j.eururo.2020.09.046
5. Antonarakis ES, Shaukat F, Isaacsson Velho P, Kaur H, Shenderov E, Pardoll DM, et al. Clinical Features and Therapeutic Outcomes in Men With Advanced Prostate Cancer and DNA Mismatch Repair Gene Mutations. Eur Urol (2019) 75(3):378–82. doi: 10.1016/j.eururo.2018.10.009
6. Hager S, Ackermann CJ, Joerger M, Gillessen S, Omlin A. Anti-Tumour Activity of Platinum Compounds in Advanced Prostate Cancer-A Systematic Literature Review. Ann Oncol (2016) 27(6):975–84. doi: 10.1093/annonc/mdw156
7. Nava Rodrigues D, Rescigno P, Liu D, Yuan W, Carreira S, Lambros MB, et al. Immunogenomic Analyses Associate Immunological Alterations With Mismatch Repair Defects in Prostate Cancer. J Clin Invest (2018) 128(10):4441–53. doi: 10.1172/JCI121924
9. Allavena P, Sica A, Solinas G, Porta C, Mantovani A. The Inflammatory Micro-Environment in Tumor Progression: The Role of Tumor-Associated Macrophages. Crit Rev Oncol Hematol (2008) 66(1):1–9. doi: 10.1016/j.critrevonc.2007.07.004
10. Crocetto F, Boccellino M, Barone B, Di Zazzo E, Sciarra A, Galasso G, et al. The Crosstalk Between Prostate Cancer and Microbiota Inflammation: Nutraceutical Products Are Useful to Balance This Interplay? Nutrients (2020) 12(9):2648. doi: 10.3390/nu12092648
13. Kim SW, Kim JS, Papadopoulos J, Choi HJ, He J, Maya M, et al. Consistent Interactions Between Tumor Cell IL-6 and Macrophage TNF-α Enhance the Growth of Human Prostate Cancer Cells in the Bone of Nude Mouse. Int Immunopharmacol (2011) 11(7):862–72. doi: 10.1016/j.intimp.2011.01.004
14. Guan W, Hu J, Yang L, Tan P, Tang Z, West BL, et al. Inhibition of TAMs Improves the Response to Docetaxel in Castration-Resistant Prostate Cancer. Endocr Relat Cancer (2019) 26(1):131–40. doi: 10.1530/ERC-18-0284
15. Zarif JC, Baena-Del Valle JA, Hicks JL, Heaphy CM, Vidal I, Luo J, et al. Mannose Receptor-Positive Macrophage Infiltration Correlates With Prostate Cancer Onset and Metastatic Castration-Resistant Disease. Eur Urol Oncol (2019) 2(4):429–36. doi: 10.1016/j.euo.2018.09.014
16. Xu N, Chen SH, Lin TT, Cai H, Ke ZB, Dong RN, et al. Development and Validation of Hub Genes for Lymph Node Metastasis in Patients With Prostate Cancer. J Cell Mol Med (2020) 24(8):4402–14. doi: 10.1111/jcmm.15098
17. Chen YH, Chen SH, Hou J, Ke ZB, Wu YP, Lin TT, et al. Identifying Hub Genes of Clear Cell Renal Cell Carcinoma Associated With the Proportion of Regulatory T Cells by Weighted Gene Co-Expression Network Analysis. Aging (Albany NY) (2019) 11(21):9478–91. doi: 10.18632/aging.102397
20. Chen YH, Lin TT, Wu YP, Li XD, Chen SH, Xue XY, et al. Identification of Key Genes and Pathways in Seminoma by Bioinformatics Analysis. Onco Targets Ther (2019) 12:3683–93. doi: 10.2147/OTT.S199115
22. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat Commun (2019) 10(1):1523. doi: 10.1038/s41467-019-09234-6
23. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: Protein-Protein Interaction Networks, Integrated Over the Tree of Life. Nucleic Acids Res (2015) 43(Database issue):D447–52. doi: 10.1093/nar/gku1003
24. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res (2003) 13(11):2498–504. doi: 10.1101/gr.1239303
25. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, et al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia (2017) 19(8):649–58. doi: 10.1016/j.neo.2017.05.002
27. Huang Y, Yuan J, Righi E, Kamoun WS, Ancukiewicz M, Nezivar J, et al. Vascular Normalizing Doses of Antiangiogenic Treatment Reprogram the Immunosuppressive Tumor Microenvironment and Enhance Immunotherapy. Proc Natl Acad Sci USA (2012) 109(43):17561–6. doi: 10.1073/pnas.1215397109
29. Miyasato Y, Shiota T, Ohnishi K, Pan C, Yano H, Horlad H, et al. High Density of CD204-Positive Macrophages Predicts Worse Clinical Prognosis in Patients With Breast Cancer. Cancer Sci (2017) 108(8):1693–700. doi: 10.1111/cas.13287
30. Kubota K, Moriyama M, Furukawa S, Rafiul HASM, Maruse Y, Jinno T, et al. CD163+CD204+ Tumor-Associated Macrophages Contribute to T Cell Regulation Via Interleukin-10 and PD-L1 Production in Oral Squamous Cell Carcinoma. Sci Rep (2017) 7(1):1755. doi: 10.1038/s41598-017-01661-z
31. Komohara Y, Hasita H, Ohnishi K, Fujiwara Y, Suzu S, Eto M, et al. Macrophage Infiltration and Its Prognostic Relevance in Clear Cell Renal Cell Carcinoma. Cancer Sci (2011) 102(7):1424–31. doi: 10.1111/j.1349-7006.2011.01945.x
34. Zhang Q, Xia J, Wang Y, Zhang J, Ji C, Cong R, et al. Tumor Infiltrating M2 Macrophages Could Predict Biochemical Recurrence of Localized Prostate Cancer After Radical Prostatectomy. Exp Cell Res (2019) 384(1):111588. doi: 10.1016/j.yexcr.2019.111588
35. Roca H, Hernandez J, Weidner S, McEachin RC, Fuller D, Sud S, et al. Transcription Factors OVOL1 and OVOL2 Induce the Mesenchymal to Epithelial Transition in Human Cancer. PloS One (2013) 8(10):e76773. doi: 10.1371/journal.pone.0076773
37. Nicolazzo C, Busetto GM, Del Giudice F, Sperduti I, Giannarelli D, Gradilone A, et al. The Long-Term Prognostic Value of Survivin Expressing Circulating Tumor Cells in Patients With High-Risk Non-Muscle Invasive Bladder Cancer (NMIBC). J Cancer Res Clin Oncol (2017) 143(10):1971–6. doi: 10.1007/s00432-017-2449-8
38. Nicolazzo C, Busetto GM, Gradilone A, Sperduti I, Del Giudice F, Loreni F, et al. Circulating Tumor Cells Identify Patients With Super-High-Risk Non-Muscle-Invasive Bladder Cancer: Updated Outcome Analysis of a Prospective Single-Center Trial. Oncologist (2019) 24(5):612–6. doi: 10.1634/theoncologist.2018-0784
39. Ferro M, Di Lorenzo G, Vartolomei MD, Bruzzese D, Cantiello F, Lucarelli G, et al. Absolute Basophil Count Is Associated With Time to Recurrence in Patients With High-Grade T1 Bladder Cancer Receiving Bacillus Calmette-Guérin After Transurethral Resection of the Bladder Tumor. World J Urol (2020) 38(1):143–50. doi: 10.1007/s00345-019-02754-2
40. Sciarra A, Gentilucci A, Silvestri I, Salciccia S, Cattarino S, Scarpa S, et al. Androgen Receptor Variant 7 (AR-V7) in Sequencing Therapeutic Agents for Castratrion Resistant Prostate Cancer: A Critical Review. Medicine (Baltimore) (2019) 98(19):e15608. doi: 10.1097/MD.0000000000015608
41. Cui M, Wang Y, Sun B, Xiao Z, Ye L, Zhang X, et al. MiR-205 Modulates Abnormal Lipid Metabolism of Hepatoma Cells Via Targeting acyl-CoA Synthetase Long-Chain Family Member 1 (ACSL1) mRNA. Biochem Biophys Res Commun (2014) 444(2):270–5. doi: 10.1016/j.bbrc.2014.01.051
44. Shi YX, Yin JY, Shen Y, Zhang W, Zhou HH, Liu ZQ, et al. Genome-Scale Analysis Identifies NEK2, DLGAP5 and ECT2 as Promising Diagnostic and Prognostic Biomarkers in Human Lung Cancer. Sci Rep (2017) 7(1):8072. doi: 10.1038/s41598-017-08615-5
45. Yamamoto S, Takayama KI, Obinata D, Fujiwara K, Ashikari D, Takahashi S, et al. Identification of New Octamer Transcription Factor 1-Target Genes Upregulated in Castration-Resistant Prostate Cancer. Cancer Sci (2019) 110(11):3476–85. doi: 10.1111/cas.14183
46. Hewit K, Sandilands E, Martinez RS, James D, Leung HY, Bryant DM, et al. A Functional Genomics Screen Reveals a Strong Synergistic Effect Between Docetaxel and the Mitotic Gene DLGAP5 That Is Mediated by the Androgen Receptor. Cell Death Dis (2018) 9(11):1069. doi: 10.1038/s41419-018-1115-7
47. Liu Y, Chen H, Dong P, Xie G, Zhou Y, Ma Y, et al. KIF23 Activated Wnt/β-Catenin Signaling Pathway Through Direct Interaction With Amer1 in Gastric Cancer. Aging (Albany NY) (2020) 12(9):8372–96. doi: 10.18632/aging.103146
48. Li T, Li Y, Gan Y, Tian R, Wu Q, Shu G, et al. Methylation-Mediated Repression of MiR-424/503 Cluster Promotes Proliferation and Migration of Ovarian Cancer Cells Through Targeting the Hub Gene KIF23. Cell Cycle (2019) 18(14):1601–18. doi: 10.1080/15384101.2019.1624112
49. Zhang Z, Xie H, Zhu S, Chen X, Yu J, Shen T, et al. High Expression of KIF22/Kinesin-Like DNA Binding Protein (Kid) as a Poor Prognostic Factor in Prostate Cancer Patients. Med Sci Monit (2018) 24:8190–7. doi: 10.12659/MSM.912643
50. Sircar K, Huang H, Hu L, Liu Y, Dhillon J, Cogdell D, et al. Mitosis Phase Enrichment With Identification of Mitotic Centromere-Associated Kinesin as a Therapeutic Target in Castration-Resistant Prostate Cancer. PloS One (2012) 7(2):e31259. doi: 10.1371/journal.pone.0031259
52. Gong C, Ai J, Fan Y, Gao J, Liu W, Feng Q, et al. NCAPG Promotes the Proliferation of Hepatocellular Carcinoma Through PI3K/AKT Signaling. Onco Targets Ther (2019) 12:8537–52. doi: 10.2147/OTT.S217916
Keywords: M2-TAMs, prostate cancer, CIBERSORT algorithm, weighted gene co-expression network analysis, biomarkers
Citation: Xu N, Dong R-N, Lin T-T, Lin T, Lin Y-Z, Chen S-H, Zhu J-M, Ke Z-B, Huang F, Chen Y-H and Xue X-Y (2021) Development and Validation of Novel Biomarkers Related to M2 Macrophages Infiltration by Weighted Gene Co-Expression Network Analysis in Prostate Cancer. Front. Oncol. 11:634075. doi: 10.3389/fonc.2021.634075
Received: 26 November 2020; Accepted: 27 May 2021;
Published: 29 June 2021.
Edited by:Matteo Ferro, European Institute of Oncology (IEO), Italy
Reviewed by:Felice Crocetto, Federico II University Hospital, Italy
Francesco Del Giudice, Sapienza University of Rome, Italy
Copyright © 2021 Xu, Dong, Lin, Lin, Lin, Chen, Zhu, Ke, Huang, Chen and Xue. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work