Abstract
Neuroendocrine differentiation (NED) is associated with WNT signaling activation and can be significantly observed after failure of androgen-deprivation therapy (ADT) for prostatic adenocarcinomas. Cytokine signaling is stimulated in NED prostate cancer; however, how ADT-upregulated WNT signaling promotes activation of cytokine signaling and contributes to NED of prostate cancer is poorly understood. In this study, we identified ADT-mediated upregulation of transcription factor 7 like 1 (TCF7L1), which increases the cytokine response and enhances NED of prostate cancer through interleukin (IL)-8/C-X-C motif chemokine receptor type 2 (CXCR2) signaling activation. ADT induced the secretion of WNT4 which upon engagement of TCF7L1 in prostate cancer cells, enhanced IL-8 and CXCR2 expressions. TCF7L1 directly binds to the regulatory sequence region of IL-8 and CXCR2 through WNT4 activation, thus upregulating IL-8/CXCR2 signaling-driven NED and cell motility. Analysis of prostate tissue samples collected from small-cell neuroendocrine prostate cancer (SCPC) and castration-resistant prostate cancer (CRPC) tumors showed an increased intensity of nuclear TCF7L1 associated with CXCR2. Our results suggest that induction of WNT4/TCF7L1 results in increased NED and malignancy in prostate cancer that is linked to dysregulation of androgen receptor signaling and activation of the IL-8/CXCR2 pathway.
Similar content being viewed by others
Introduction
The androgen receptor (AR) signaling pathway was shown to play an essential role in prostate cancer (PCa) [1]. Although targeting the AR is initially clinically effective, most tumors become resistant to androgen-deprivation therapy (ADT) within a few years, due to progression to hormone-refractory or castration-resistant PCa (CRPC) [2]. Reduced AR expression and signaling support increased neuroendocrine (NE) differentiation (NED) of PCa and is associated with the emergence of an aggressive phenotype [3]. Loss of AR signaling in NED PCa (NEPC) leads to decreased levels of serum prostate-specific antigen (PSA), which increases difficulty in early diagnoses [3]. Therefore, identifying consensus molecular mechanisms of NED that are modulated by inhibition of AR signaling is an urgent issue in PCa therapy. In NEPC induced by ADT, several models were proposed [4,5,6], including WNT [7], SRC [8], and leukemia inhibitory factor (LIF)/interleukin (IL)-6/signal transduction and activator of transcription 3 (STAT3) [9, 10] signaling activation, tumor suppressor (e.g., retinoblastoma 1 (RB1), p53, or phosphatase and tensin homolog (PTEN)) loss of function [11], epigenetic disorders [12], upregulation of transcription factors (e.g., MYCN [13], FOXA2 [14], or ZBTB46 [15]), epithelial-mesenchymal transition (EMT)-stemness crosstalk [16], and hypoxia [17].
It was reported that androgen deprivation promotes an immune response and causes an increase in interleukin (IL)-8 in human PCa cells [18]. IL-8 is thought to be one of the factors contributing to the progression of the castration-resistant state of PCa [19]. IL-8 is specifically expressed by NE tumor cells of human PCa tissues [20]. Upregulation of IL-8 and its receptor, C-X-C chemokine receptor type 2 (CXCR2), and their autocrine actions are critical factors in promoting the progression and metastasis of colon cancer cells [21]. Overexpression of CXCR2 was reported to promote invasion and metastasis of lung adenocarcinomas [22]. Despite IL-8 and CXCR2 having been found to be highly expressed in NE cells of PCa [20], the mechanisms through which upregulation of IL-8/CXCR2 signaling promotes NED and malignant progression of prostate adenocarcinomas post-ADT are not clear.
WNT is a family of signaling proteins that has pivotal roles in neurodevelopmental processes and tumor progression [23]. Activation of the canonical WNT pathway leads to transactivation of target genes, which are involved in the NED of PCa [24]. Inhibition of WNT signaling has the potential to reduce the NED of PCa that drives cancer progression and thus could improve therapeutic outcomes [25]. Nevertheless, activation of AR signaling can repress WNT-mediated transcription induced by androgen in PCa cells [26]. Those findings suggest that inhibition of AR expression might activate the highly upregulated WNT signaling pathway. Therefore, the relationship between the AR and WNT signaling pathway needs to be further explored in order to characterize possible NED mechanisms in PCa after ADT.
Transcription factor (TCF) 7 like 1 (TCF7L1) is a member of the TCF/lymphoid enhancer factor (LEF) family of transcription factors that participate in maintaining stem cell pluripotency [27], homeostasis of skin epithelial tissues [28], determination of cell lineages during stomach formation [29], and vertebrate brain development [30]. It was shown that the growth and invasive capability of highly metastatic breast cancer cells were impaired by ectopic TCF7L1 expression [31]. Overexpression of TCF7L1 was also found to induce the growth of colorectal tumors [32]. TCF7L1 is a candidate oncogene and prognostic marker for skin squamous cell carcinoma (SCC) [33]. Herein, we aimed to understand if TCF7L1 is highly upregulated in prostate tumors after ADT, compared to hormone-naïve prostate patients. We investigated whether there is an association between TCF7L1 overexpression and the NED potential post-ADT in PCa patients. Because activation of IL-8/CXCR2 signaling may increase the induction of CRPC after ADT [34], we showed that TCF7L1 could be the activator of IL-8 and CXCR2 in NED PCa cells through WNT4 activation. We studied the crosstalk mechanism of WNT4/TCF7L1 and IL-8/CXCR2 signals in the NED of PCa to understand the pathogenic mechanism of TCF7L1-dependent activation that predisposes PCa resistance to ADT.
Results
TCF7L1 is upregulated in CXCR2+ NE cells and correlated with CXCR2 in small-cell NE PCa (SCPC) and CRPC samples
As to the association between WNT signaling and CXCR2, it is known that CXCR2 may participate in regulating the lineage commitment of human pluripotent stem cells by activating the β-catenin pathway [35]; however, the mechanism by which CXCR2 regulates NED of PCa via WNT signaling is not clear. Although CXCR2 is highly expressed in NE cells of PCa [20], it is still unclear how WNT signaling upregulates CXCR2 and promotes NED of PCa after ADT. We analyzed RNA-sequencing (RNA-Seq) data (GSE114326) expression in a set of prostatic tumor samples collected from CXCR2+ NE cells compared to CXCR2− luminal cells [36]. We focused on the TCF family because we previously showed that TCF transcription factors may promote the malignant progression of PCa after ADT and induce treatment resistance [37], but their NED function is unknown. Unexpectedly, we observed that TCF7L1, a member of the TCF family, was significantly upregulated in CXCR2+ NE cells compared to TCF7 and TCF7L2 (Fig. 1A, Supplementary Fig. S1A). In PCa cells, the full-length TCF7, TCF7L1, and TCF7L2 complementary (c)DNAs were transiently transfected into LNCaP cells and analyzed for CXCR2, IL-8, and NE markers (chromogranin A (CHGA), chromogranin B (CHGB), synaptophysin (SYP), and enolase 2 (ENO2)). We found that overexpression of TCF7L1 significantly increased expressions of CXCR2, IL-8, and NE markers, but TCF7 and TCF7L2 only partially increased the expression of IL-8 or CXCR2, but not the expressions of NE markers (Supplementary Fig. S1B–D). A geneset enrichment analysis (GSEA) confirmed that TCF7L1 was highly expressed in PCa tissues with upregulated CXCR2+ NE-responsive gene signatures [36] and NE-responsive gene signatures [38] in The Cancer Genome Atlas (TCGA) PCa dataset (Fig. 1B). Furthermore, TCF7L1 upregulation was positively associated with the gene signature of neural development (gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG)) and negatively associated with the gene signature of the androgen response (GO, pathway interaction database (PID), Wang [39], and Nelson [40]) as confirmed by a GSEA from TCGA PCa database (Fig. 1C), suggesting that TCF7L1 upregulation after ADT is involved in NED progression. NE cells were shown to accumulate during PCa progression after ADT through interfering with NED [41, 42], and they undergo progression into advanced PCa as a result of SCPC [43]. We sought to validate the positive correlation between CXCR2 and TCF7L1 in advanced PCa tissue samples. We measured TCF7L1 and CXCR2 protein levels in SCPC samples collected from Taipei Medical University-Wan Fang Hospital (Taipei, Taiwan) by immunohistochemical (IHC) staining using consecutive tissue sections, and found that there was a positive relationship between TCF7L1 and CXCR2 (Fig. 1D, E). We next analyzed a CRPC tissue microarray (TMA) collected by the Department of Pathology, Duke University School of Medicine (Durham, NC, USA), and we found a positive association between TCF7L1 and CXCR2 in CRPC samples (Fig. 1F, G). We further validated TCF7L1 expression with gene signatures reflecting IL-8/CXCR2 signaling components in TCGA PCa dataset and found that tissues expressing high levels of TCF7L1 were robustly associated with a gene signature of upregulated IL-8/CXCR2 responsiveness (PID, Fig. 1H). This finding suggests that the NED function is regulated by IL-8/CXCR2 and is connected to upregulation of TCF7L1 in PCa after ADT.
Activation of IL-8/CXCR2 signaling is upregulated by TCF7L1
Despite a WNT-stimulated cytokine response in the microenvironment being required to maintain aggressiveness and metastatic progression in several types of cancer [44], there is little knowledge about the upregulation of TCF7L1 associated with cytokine responses in PCa. We validated our hypotheses by a human cytokine and chemokine protein array analysis in PC3 cells with TCF7L1-knockdown (KD). We identified that a significant reduction in IL-8 was associated with TCF7L1-KD based on both cellular and secreted levels (Fig. 2A–C). We confirmed that TCF7L1-KD significantly decreased messenger (m)RNA expressions of IL-8 and CXCR2 in AR-negative PC3 cells (Fig. 2D). We stably overexpressed TCF7L1 in AR-positive LNCaP and C4-2 cells, and found that IL-8 and CXCR2 mRNAs were upregulated according to a real-time reverse-transcription quantitative polymerase chain reaction (RT-qPCR) analysis (Fig. 2E). In order to understand the relationship between TCF7L1 and CXCR2, we downloaded chromatin immunoprecipitation (ChIP)-sequence data from the Gene Expression Omnibus (GEO) (GSE80331) based on an Illumina HiSeq analysis and analyzed a putative TCF7L1-binding site located in the CXCR2 regulatory sequence. We found that there was a possible consensus TCF7L1-binding element in the CXCR2 regulatory sequence located around chr 2:218,989,000 to 218,990,000 (Supplementary Fig. S2A). Next, we looked for sequences resembling the TCF7L1 response element in the putative IL-8 and CXCR2 regulatory sequence regions and respectively found five and three candidate binding elements for TCF7L1 in the IL-8 and CXCR2 regulatory sequences (Fig. 2F). To test whether TCF7L1 directly binds to the regulatory sequence of IL-8 and CXCR2, we performed ChIP assays and used a qPCR to determine how much chromatin of IL-8 and CXCR2 was perceptible by an anti-TCF7L1 antibody vs. a positive control anti-acetyl-H3 antibody, from nuclear extracts of LNCaP cells stably expressing an empty vector (EV) or TCF7L1 cDNA vector. We found significantly increased TCF7L1 binding at the putative TCF7L1-binding site on the responsive element (RE) 2, RE3, and RE4 sites of the IL-8 gene and on the RE1 and RE3 sites of the CXCR2 gene in cells with TCF7L1 expression (Fig. 2G, H). We performed reporter assays with a DNA construct containing individual wild-type (WT) and a mutant of the putative TCF7L1 response elements from the IL-8 and CXCR2 regulatory sequences (Supplementary Fig. S2B). We found that LNCaP cells overexpressing TCF7L1 had significantly increased reporter gene activities at the RE2, RE3, and RE4 sites on the IL-8 gene and the RE1 and RE3 sites on the CXCR2 gene compared to EV-transfected cells (Fig. 2I, J). Nevertheless, decreased reporter activity was seen in cells transfected with the reporter construct containing TCF7L1-RE mutants regardless of TCF7L1 overexpression (Fig. 2I, J). Moreover, decreased IL-8 and CXCR2 reporter construct activities from the RE2, RE3, and RE4 sites on the IL-8 gene and from the RE1 and RE3 sites on the CXCR2 gene were detected in PC3 cells with TCF7L1-KD (Fig. 2K, L). These data are consistent with a mechanism whereby TCF7L1 promotes IL-8 and CXCR2 transcription by direct physical interactions between TCF7L1 and the regulatory sequences of IL-8 and CXCR2.
Suppression of AR signaling upregulates TCF7L1, which is associated with IL-8/CXCR2-driven NED PCa
Since ADT is one of the main causes of NED [3], we next determined whether upregulation of TCF7L1 was mediated by ADT. We collected PCa tissues from the same PCa patients before and after ADT in Taipei Medical University-Wan Fang Hospital (Taipei, Taiwan). IHC results showed that TCF7L1 levels were higher in PCa patients after ADT compared to the same patients before ADT (Fig. 3A, B). We validated mRNA expression levels of TCF7L1 in an RNA-Seq dataset (GSE48403) in paired PCa samples pre- and post-ADT, and found that post-ADT patients had increased TCF7L1 expression in the GSE48403 dataset (Fig. 3C), consistent with our data. Moreover, overexpression of TCF7L1 was observed in an AR antagonist-resistant cell line (C4-2-MDVR, which are C4-2 cells exposed to long-term treatment with MDV3100), and that had increased NE marker and decreased androgen-responsive gene expressions compared to parental C4-2 cells (Fig. 3D). C4-2-MDVR cells also showed increased IL-8 and CXCR2 mRNA expressions (Fig. 3D), suggesting that inhibition of AR signaling might upregulate IL-8 and CXCR2 expressions. To mimic an androgen-deprivation condition, AR-positive LNCaP and C4-2 cells were cultured with charcoal-stripped serum (CSS)-containing medium and further responded to activate AR signaling by AR ligand dihydrotestosterone (DHT) treatment. Notably, cells treated with CSS-containing medium exhibited increased mRNA and protein levels of TCF7L1, IL-8, CXCR2, and NE markers, and decreased levels of androgen-responsive markers compared to cells treated with fetal bovine serum (FBS)-containing medium (Fig. 3E, F). Importantly, TCF7L1, IL-8, CXCR2, and NE markers decreased, and androgen-responsive markers increased in DHT-treated cells (Fig. 3E, F), which supports TCF7L1/IL-8/CXCR2 upregulating NED function in PCa cells as observed post-ADT. To determine the tumorigenic role of TCF7L1 in PCa, we used a TMA obtained from the Department of Pathology at Duke University School of Medicine (Durham, NC, USA), which was comprised of normal and PCa-progressed tumor tissues. Interestingly, high-grade tumors had moderately high TCF7L1 expression compared to low-grade and normal tissues (Fig. 3G, H). We also analyzed eight SCPC samples, and found that most of those cases had higher TCF7L1 expression (Fig. 3G, H), consistent with samples collected from Taipei Medical University-Wan Fang Hospital (Taipei, Taiwan). These results supported upregulation of TCF7L1 promoting progression of NED in PCa.
TCF7L1 promotes the NED of PCa cells
Next, we investigated the functional role of TCF7L1 in the progression of NED. We checked protein expression levels of TCF7L1, NE markers (CHGA and ENO2), and androgen-responsive genes (NKX3-1 and KLK3) in a panel of PCa cell lines. We found that AR-suppressed PC3 cells, the NE-like NE-1-8 cell line [45] (LNCaP cells with long-term androgen withdrawal), and an SCPC-like LASCPC01 cell line [13] exhibited higher expressions of TCF7L1 and NE markers and lower expressions of androgen-responsive genes than did the AR-positive LNCaP and C4-2 prostate adenocarcinoma cell lines (Fig. 4A). To determine if TCF7L1 overexpression can induce the NED of AR-positive cells, we generated LNCaP and C4-2 cells that stably expressed TCF7L1 cDNA and found that the forced expression of TCF7L1 significantly increased mRNA and protein levels of NE markers, and reduced androgen-responsive gene expressions (Fig. 4B, C). In addition, decreased mRNA and protein levels of TCF7L1 and NE markers were found in PC3 and LASCPC01 cells exposed to the TCF7L1 short hairpin (sh)RNA vector compared to a nontarget control (NC) shRNA (Fig. 4D, E). We further knocked down TCF7L1 in AR-positive LNCaP and C4-2 cells and treated those cells with CSS-containing medium to mimic ADT conditions. Results demonstrated that AR-positive cells treated with CSS-containing medium had increased NE markers; however, TCF7L1-KD cells that underwent ADT had decreased expressions of TCF7L1 and NE markers, and constitutively expressed androgen-responsive genes (Fig. 4F). Moreover, MDV3100-resistant cells had increased NE marker and decreased androgen-responsive gene expressions, whereas TCF7L1-KD in C4-2-MDVR cells significantly decreased NE markers and increased levels of androgen-responsive genes (Fig. 4G). These results support TCF7L1 expression being upregulated after AR signaling inhibition, and its overexpression is involved in the NED of PCa cells.
ADT-resistance increases TCF7L1 and upregulates NED through increased WNT4
To determine the specific WNT ligand involved in ADT resistance of PCa cells, we analyzed expression profiles of WNT signaling pathway-related genes in MDV3100-resistant C4-2 cells and parental C4-2 cells via an RT [18] profiler PCR array analysis. We found that upregulation of WNT3A, WNT4, and WNT9B mRNAs was enriched in C4-2-MDVR cells compared to parental C4-2 cells (Fig. 5A). The RT-qPCR analysis confirmed that significant increases in WNT3A, WNT4, and WNT9B mRNAs in C4-2-MDVR cells were associated with TCF7L1 compared to parental C4-2 cells (Fig. 5B). In order to test the association of the upregulation of WNT3A, WNT4, or WNT9B with TCF7L1 expression and AR signal inactivation, we examined expressions of WNT3A, WNT4, and WNT9B in AR-positive LNCaP cells in response to ADT. We found that WNT4, but not WNT3A or WNT9B, was upregulated after androgen withdrawal (CSS treatment) and was downregulated in the presence of DHT (Fig. 5C), suggesting that WNT4 overexpression was linked to AR signaling inactivation in response to ADT. Treatment of LNCaP and C4-2 cells with the WNT4 recombinant protein led to a dose-dependent increase in TCF7L1, which was related to increased expressions of IL-8, CXCR2, and NE markers and decreased expressions of androgen-responsive genes (Fig. 5D). Moreover, we examined expressions of IL-8, CXCR2, and NE markers in cells treated with WNT4 protein after pretreatment with the porcupine inhibitor, IWP-O1. It is known that IWP-O1 can prevent secretion of the WNT protein [46]. Our results indicated that IWP-O1 downregulated expressions of IL-8, CXCR2, and NE markers. We also noticed that although cells are treated with IWP-O1, WNT4 can still increase the expression of IL-8, CXCR2, and NE markers compared to treatment with IWP-O1 alone (Supplementary Fig. S1E). These results support an association between activated TCF7L1 and the IL-8/CXCR2 signaling pathway through WNT4 upregulation. Although MDV3100-resistant cells showed induction of TCF7L1 and NE markers and reductions of androgen-responsive genes, we found that when WNT4 was knocked down, MDV3100-resistant cells had reduced expressions of TCF7L1 and NE markers and increased expressions of androgen-responsive genes (Fig. 5E). Consistently, CSS-treated cells had increased TCF7L1 and NE markers and decreased androgen-responsive gene expressions, while in cells with WNT4-KD, levels of TCF7L1 and NE markers were reduced, and those of androgen-responsive genes were increased. (Fig. 5F). These data suggest that WNT4 stimulation is related to activation of TCF7L1-driven NED in PCa after ADT. We next knocked down TCF7L1 and further treated cells with the WNT4 protein, and found that TCF7L1-KD cells exhibited abolished expressions of WNT4-mediated IL-8 and CXCR2 (Fig. 5G, H), supporting that WNT4 upregulation associated with IL8/CXCR2 expression is TCF7L1-dependent in PCa cells. It was shown that targeting WNT3A signaling reduces β-catenin expression and may be considered to have therapeutic potential in NEPC [47]. We further investigated whether WNT3A can induce IL-8/CXCR2 expression in a TCF7L1-dependent manner by treating TCF7L1-KD cells with the WNT3A protein. We found that WNT3A could still upregulate IL-8 and CXCR2 expressions, regardless of whether TCF7L1 was knocked down in cells (Supplementary Fig. S1F). This result indicates that WNT3A may also mediate IL-8/CXCR2 expression, but through a TCF7L1-independent pathway. In addition, LNCaP and C4-2 cells treated with the WNT4 protein had increased NE marker protein levels and decreased androgen-responsive genes; however, cells with TCF7L1-KD did not show these effects (Fig. 5I). To test the effect of WNT4 on transcriptional activities of the IL-8 and TCF7L1 genes, results from the ChIP assay showed that TCF7L1-binding signals of the IL-8 and CXCR2 genes increased in cells in response to WNT4 protein treatment, whereas TCF7L1-KD decreased this effect (Supplementary Fig. S2C, D). Moreover, increased IL-8 and CXCR2 reporter activities were found in cells with WNT4 protein treatment; however, TCF7L1-KD attenuated the ability of WNT4 to upregulate IL-8 and CXCR2 reporter activities according to reporter assays (Supplementary Fig. S2E, F). These data are consistent with our hypothesis that WNT4 upregulation is associated with TCF7L1-mediated transcription of the IL-8 and CXCR2 genes.
WNT4/TCF7L1 enhances cell migration and invasion of PCa cells
In order to evaluate the contribution of WNT4/TCF7L1 to the malignant progression of human prostate adenocarcinomas, AR-positive C4-2 and LNCaP cells stably expressing TCF7L1 were treated with the WNT4 protein. The roles of WNT4/TCF7L1 in driving cell migration, invasion, and proliferation were further tested. Results showed that cells with TCF7L1 overexpression or WNT4 treatment had significantly increased cell migration and invasion compared to phosphate-buffered saline (PBS)-treated cells expressing an EV, while TCF7L1-expressing cells treated with WNT4 had the strongest cell migration and invasion abilities (Fig. 6A–D). These results confirmed the synergistic effect of WNT4 and TCF7L1 of increasing cell migration and invasion. Interestingly, we found that cells with TCF7L1 overexpression did not exhibit an increased proliferation rate even in cells treated with the WNT4 protein (Supplementary Fig. S3A, B), indicating that TCF7L1 mainly affects cell migration rather than cell proliferation. In addition, AR-negative PC3 cells with WNT4 protein treatment had increased cell migration and invasion capabilities according to cell migration and invasion assays; however, cells expressing TCF7L1-KD did not exhibit these effects (Fig. 6E–H). Notably, there were no significant changes in cell proliferation of PC3 cells expressing TCF7L1-KD, even in cells with WNT4 treatment (Supplementary Fig. S3C). These results were supported by further in vivo experiments. Mice that were administered subcutaneous injections of PC3 cells harboring TCF7L1-KD showed no changes in tumor growth rates or tumor weights compared to mice injected with cells carrying the NC shRNA (Supplementary Fig. S3D–F). To verify the ability of TCF7L1 to upregulate IL-8/CXCR2 and mediate cellular NED, tumors were harvested and examined for TCF7L1, IL-8, CXCR2, and NE and proliferation markers by IHC. Results showed that tumors from PC3 cells with TCF7L1-KD had decreased TCF7L1, IL-8, CXCR2, and NE marker expressions compared to NC-expressing cells, but there was no significant difference in cell proliferation markers (Fig. 6I, J), confirming that TCF7L1 is an important factor that induces IL-8/CXCR2-associated NED in PCa cells. These data suggest that WNT4/TCF7L1 promotes a variety of cell motility properties and the NED of PCa cells. Taken together, we demonstrated that inhibiting AR activities increased TCF7L1 expression, and induction of TCF7L1 through upregulation of WNT4 relieved activation of the oncogenic roles of IL-8 and CXCR2, subsequently leading to the development of malignant progression and the NED of PCa (Fig. 6K).
Discussion
Evidence indicates that AR inhibition in PCa after ADT can activate the WNT signaling pathway [26, 48, 49]. Exploring the relationship between AR pathway inhibition and WNT transcriptional activity is important for patients with PCa who may require different management strategies for their disease [50]. We studied WNT ligands that are upregulated after ADT, which may be related to the progression of NED in PCa cells undergoing MDV3100 resistance. Indeed, ß-catenin and various WNT ligands, including WNT4, WNT5A, WNT7B, and WNT11, were shown to play roles in the development of ADT resistance in PCa cells [51,52,53]. In addition, it was reported that WNT7B and WNT11 can increase expressions of NE markers in PCa cells [24, 54]. Although secretion of WNT ligands is common in PCa after ADT treatment [7], each ligand may induce ADT resistance and NED progression through different signaling pathways. Our results demonstrated the role of WNT4 in upregulating TCF7L1-dependent NED and linked TCF7L1 to IL-8/CXCR4-driven NED of PCa.
Activation of WNT signaling is highly associated with NED due to their impacts on neurodevelopment [44]. Our data support the drug-resistance role of WNT4 in PCa. We indentified the regulatory mechanism underlying the impact of WNT4 in TCF7L1-driven IL-8/CXCR2 signaling on the NED of PCa after ADT. Our model was supported by ectopic WNT4 expression and the addition of exogenous WNT4. Both of these had indirect effects on WNT4 binding to secreted frizzled (FZD)-related protein (sFRP) family proteins or by replacing other endogenous WNT proteins with the opposite effect of binding to FZD receptors [55, 56]. We focused on WNT4 because it is more highly expressed in ADT-resistant cells. Our data confirmed that WNT4 gene silencing in MDV3100-resistant C4-2 and CSS-treated LNCaP cells reduced expressions of TCF7L1 and NE markers, which provides support for the link between endogenous WNT4 and TCF7L1 expressions in PCa cells under AR deprivation. Our study suggests that inhibition of AR-increased WNT4 upregulation activates TCF7L1-driven IL-8/CXCR2 signaling which is associated with PCa NED progression. We explored the regulatory mechanisms of TCF7L1 which link suppression of the AR pathway with upregulated WNT4 signaling, and defined the roles of TCF7L1 in this crosstalk leading to progression of NED PCa by upregulating IL-8 and CXCR2.
The relationship between cytokine signaling and PCa progression has been investigated for several years [57]. However, the connection between WNT signaling and cytokine responses in NED PCa cells is still unclear, and effective cytokine-targeted therapies for NED PCa are also still unavailable. IL-8 is known as a downstream effector of WNT signaling, which is essential for cancer progression, metastasis, and development of chemoresistance [58]. As ovarian tumor cells undergo the EMT, those cells activate WNT signaling to promote IL-8 expression and further reinforce an IL-8/CXCR2/WNT feedback loop [59]. The IL-8/CXCR2/WNT feedback loop was also observed in the development of chemoresistance and metastasis in estrogen receptor (ER)-positive breast cancer, in which WNT signaling is specifically mediated by WNT3A coupled with SNAI1/nuclear factor (NF)-κB signaling activation [60]. In PCa, we demonstrated that IL-8/CXCR2 signaling leads to interactions with WNT4/TCF7L1 in the NED progressive signaling pathway, thereby mediating the ability of aggressive phenotypes of PCa post-ADT. Our study provides a functional link between low AR activities and increased TCF7L1, which may be associated with WNT4 stimulation that contributes to TCF7L1-driven PCa NED progression. This supports low AR activity resulting in elevation of the WNT pathway at low androgen levels and in turn promoting castration-refractory and NED PCa progression [25, 26, 48, 49]. Given that IL-8/CXCR2 overexpression in NE cells enhances a high metastatic potential [36], we demonstrated that the activated WNT4/TCF7L1 pathway promotes NED and cytokine responses of PCa through upregulating IL-8/CXCR2 signaling.
In summary, our report supports the increased aggressiveness of PCa cells after ADT possibly being accompanied by cellular NED. We demonstrated that ADT may induce the upregulation of WNT4/TCF7L1 signaling, leading to the activation of IL-8/CXCR2 to drive the progression of NED. Inhibition of AR signaling is sufficient to increase WNT4 and TCF7L1 expressions, which implies a dominant biological function of WNT4/TCF7L1 for NED progression of PCa in the WNT signaling regulatory network. Our data provide a mechanism whereby high levels of TCF7L1 post-ADT result in increased IL-8 and CXCR2 levels through upregulation of WNT4, in turn leading to NED progression in PCa. A novel regulatory role for the candidate prostate tumor promoter, TCF7L1, was identified. We demonstrated that the interaction of TCF7L1 with the IL-8/CXCR2 signaling pathway requires the presence of components of WNT4 signaling. This could explain why crosstalk between the cytokine response and WNT signaling is frequently associated with NED PCa progression.
Materials and methods
Cell lines and cell culture
The human PCa LNCaP, C4-2, and PC3 cell lines were acquired from American Type Culture Collection (ATCC; Manassas, VA, USA). All cells were cultured in RPMI 1640 medium (Thermo-Fisher, Waltham, MA, USA) supplemented with 10% fetal bovine serum (FBS). The NE-like NE-1-8 cell line was purchased from ATCC and cultured in RPMI 1640 medium supplemented with 5% charcoal-stripped serum (CSS) and Dulbecco’s modified Eagle medium (DMEM) supplemented with 10% FBS, 0.005 mg/ml bovine insulin (Sigma–Aldrich, Darmstadt, Germany), and 10 nM dihydrotestosterone (DHT; Sigma–Aldrich). The LASCPC01 cell line was obtained from ATCC and was cultured in HITES medium which was composed by RPMI 1640 medium, 1× insulin-transferrin-selenite (ITS-G, Thermo-Fisher), 2 mM GlutaMAX (Thermo-Fisher), 10 μM hydrocortisone (Sigma–Aldrich), 10 μM β-estradiol (Sigma–Aldrich), 1× penicillin/streptomycin (Thermo-Fisher), and 5% FBS (Thermo-Fisher). The MDV3100-resistant C4-2-MDVR cell line is a viable cell line generated by growing C4-2 cells under selective pressure of 20 μM enzalutamide (MDV3100) for 6 months. For DHT treatment, cells were treated with 10 nM DHT (Sigma–Aldrich) for 24 h in 10% CSS-containing medium. All cell lines tested negative for mycoplasma contamination. For WNT recombinant protein treatment, cells were treated with 100 ng/ml of human WNT4 or WNT3A (R&D Systems, Minneapolis, MN, USA) for 48 h, in 10% CSS-containing medium. For the porcupine inhibitor IWP-O1 treatment, cells were treated with 1 μM IWP-O1 (Sigma–Aldrich) for 48 h in 10% FBS-containing medium. All cellular experiments were carried out within 20 passages to ensure consistency and uniformity.
Real-time reverse-transcription quantitative polymerase chain reaction (RT-qPCR)
The experimental procedures are detailed in “Supplementary Information”.
Immunohistochemical (IHC) staining
Clinical samples consisted of 11 independent SCPC tumors collected from Taipei Medical University-Wan Fang Hospital (Taipei, Taiwan). PCa TMA sections, including 12 normal prostatic epithelial samples, 77 primary prostate adenocarcinomas, and eight SCPCs, and CRPC TMA sections, including 25 CRPC patients, were provided by Duke University School of Medicine (Durham, NC, USA). Eighteen PCa samples from the same patients before and after ADT were collected from Taipei Medical University-Wan Fang Hospital (Taipei, Taiwan). Tissue samples were used in accordance with the U.S. Common Rule and the Declaration of Helsinki, and their use was approved by the Duke University School of Medicine Institutional Review Board (protocol ID: Pro00070193) and the Taipei Medical University Joint Institutional Review Board (approval no.: N202001017). IHC was performed using TCF7L1 (14519-1-AP; Proteintech, Rosemont, IL, USA) and CXCR2 (555932; BD Biosciences, San Jose, CA, USA) antibodies at respective 1:150 and 1:200 dilutions as described in “Supplementary Materials and Methods”.
Cytokine array assay
A cytokine assay was performed using the human cytokine array kit (ARY005B, R&D Systems) following the protocol in the user’s manual as described in “Supplementary Materials and Methods”.
Chromatin immunoprecipitation (ChIP) assay
The experimental procedures are detailed in “Supplementary Information”.
Promoter reporter assay
The experimental procedures are detailed in “Supplementary Information”.
Western blot analysis
The experimental procedures are detailed in “Supplementary Information”.
Polymerase chain reaction (PCR)-array analysis
The experimental procedures are detailed in “Supplementary Information”.
Migration and invasion assays
The experimental procedures are detailed in “Supplementary Information”.
Cell proliferation assay
The experimental procedures are detailed in “Supplementary Information”.
Tumorigenicity assays in mice
Animal work was performed in accordance with a protocol approved by the Taipei Medical University Animal Care and Use Committee (approval no.: LAC-2019-0465, Taipei, Taiwan). Five-week-old male nude mice (NLAC, Taipei, Taiwan) were randomly subcutaneous injected with TCF7L1 short hairpin (sh)RNA or nontarget control (NC) vector-transfected PC3 cells in 50% MatrigelTM (BD Biosciences). The tumor volume was measured weekly with calipers, and the tumor volume was calculated using the following formula: tumor volume = (4/3) × (L/2) (W/2)2, where L is the length and W the width. The results are presented as the mean ± standard error (SE) for each experimental group. At the end of the 40th day, tumors of mice injected with PC3 cells expressing TCF7L1-knockdown (KD) or the control vector were collected, and the weights were measured. Tumors were fixed in 4% paraformaldehyde (PFA) in PBS, and IHC staining was performed.
Statistical analysis
Experimental results were collected and plotted by GraphPad Prism v8.0 (GraphPad Software, San Diego, CA, USA). Data from each group were compared using an analysis of variance (ANOVA) and Dunnett’s test for a post-hoc analysis compared to the control groups. All data are presented as the mean ± standard error of the mean (SEM). Data groups significantly differing from the control were appropriately labeled for comparisons among three or more groups. A log-rank test was used for the survival curve analysis. p values of <0.05 were considered statistically significant.
References
Yuan X, Cai C, Chen S, Chen S, Yu Z, Balk SP. Androgen receptor functions in castration-resistant prostate cancer and mechanisms of resistance to new agents targeting the androgen axis. Oncogene. 2014;33:2815–25.
Katzenwadel A, Wolf P. Androgen deprivation of prostate cancer: leading to a therapeutic dead end. Cancer Lett. 2015;367:12–17.
Beltran H, Hruszkewycz A, Scher HI, Hildesheim J, Isaacs J, Yu EY, et al. The role of lineage plasticity in prostate cancer therapy resistance. Clin Cancer Res. 2019;25:6916–24.
Berchuck JE, Viscuse PV, Beltran H, Aparicio A. Clinical considerations for the management of androgen indifferent prostate cancer. Prostate Cancer Prostatic Dis. 2021;24:623–37.
Tong D. Unravelling the molecular mechanisms of prostate cancer evolution from genotype to phenotype. Crit Rev Oncol Hematol. 2021;163:103370.
Wang Y, Wang Y, Ci X, Choi SYC, Crea F, Lin D, et al. Molecular events in neuroendocrine prostate cancer development. Nat Rev Urol. 2021.18:581–96. https://pubmed.ncbi.nlm.nih.gov/34290447/.
Bland T, Wang J, Yin L, Pu T, Li J, Gao J, et al. WLS-Wnt signaling promotes neuroendocrine prostate cancer. iScience. 2021;24:101970.
Gururajan M, Cavassani KA, Sievert M, Duan P, Lichterman J, Huang JM, et al. SRC family kinase FYN promotes the neuroendocrine phenotype and visceral metastasis in advanced prostate cancer. Oncotarget. 2015;6:44072–83.
Liu YN, Niu S, Chen WY, Zhang Q, Tao Y, Chen WH, et al. Leukemia inhibitory factor promotes castration-resistant prostate cancer and neuroendocrine differentiation by activated ZBTB46. Clin Cancer Res. 2019;25:4128–40.
Sun F, Zhang ZW, Tan EM, Lim ZLR, Li Y, Wang XC, et al. Icaritin suppresses development of neuroendocrine differentiation of prostate cancer through inhibition of IL-6/STAT3 and Aurora kinase A pathways in TRAMP mice. Carcinogenesis. 2016;37:701–11.
Mu P, Zhang Z, Benelli M, Karthaus WR, Hoover E, Chen CC, et al. SOX2 promotes lineage plasticity and antiandrogen resistance in TP53- and RB1-deficient prostate cancer. Science. 2017;355:84–88.
Zhang Y, Zheng D, Zhou T, Song H, Hulsurkar M, Su N, et al. Androgen deprivation promotes neuroendocrine differentiation and angiogenesis through CREB-EZH2-TSP1 pathway in prostate cancers. Nat Commun. 2018;9:4080.
Lee JK, Phillips JW, Smith BA, Park JW, Stoyanova T, McCaffrey EF, et al. N-Myc drives neuroendocrine prostate cancer initiated from human prostate epithelial cells. Cancer Cell. 2016;29:536–47.
Park JW, Lee JK, Witte ON, Huang J. FOXA2 is a sensitive and specific marker for small cell neuroendocrine carcinoma of the prostate. Mod Pathol. 2017;30:1262–72.
Chen WY, Wen YC, Lin SR, Yeh HL, Jiang KC, Chen WH, et al. Nerve growth factor interacts with CHRM4 and promotes neuroendocrine differentiation of prostate cancer and castration resistance. Commun Biol. 2021;4:22.
Soundararajan R, Paranjape AN, Maity S, Aparicio A, Mani SA. EMT, stemness and tumor plasticity in aggressive variant neuroendocrine prostate cancers. Biochim Biophys Acta Rev Cancer. 2018;1870:229–38.
Guo H, Ci X, Ahmed M, Hua JT, Soares F, Lin D, et al. ONECUT2 is a driver of neuroendocrine prostate cancer. Nat Commun. 2019;10:278.
Araki S, Omori Y, Lyn D, Singh RK, Meinbach DM, Sandman Y, et al. Interleukin-8 is a molecular determinant of androgen independence and progression in prostate cancer. Cancer Res. 2007;67:6854–62.
Gueron G, De Siervi A, Vazquez E. Advanced prostate cancer: reinforcing the strings between inflammation and the metastatic behavior. Prostate Cancer Prostatic Dis. 2012;15:213–21.
Huang J, Yao JL, Zhang L, Bourne PA, Quinn AM, di Sant’Agnese PA, et al. Differential expression of interleukin-8 and its receptors in the neuroendocrine and non-neuroendocrine compartments of prostate cancer. Am J Pathol. 2005;166:1807–15.
Lee YS, Choi I, Ning Y, Kim NY, Khatchadourian V, Yang D, et al. Interleukin-8 and its receptor CXCR2 in the tumour microenvironment promote colon cancer growth, progression and metastasis. Br J Cancer. 2012;106:1833–41.
Saintigny P, Massarelli E, Lin S, Ahn YH, Chen Y, Goswami S, et al. CXCR2 expression in tumor cells is a poor prognostic factor and promotes invasion and metastasis in lung adenocarcinoma. Cancer Res. 2013;73:571–82.
Niehrs C. The complex world of WNT receptor signalling. Nat Rev Mol Cell Biol. 2012;13:767–79.
Moparthi L, Pizzolato G, Koch S. Wnt activator FOXB2 drives the neuroendocrine differentiation of prostate cancer. Proc Natl Acad Sci USA. 2019;116:22189–95.
Shan J, Al-Muftah MA, Al-Kowari MK, Abuaqel SWJ, Al-Rumaihi K, Al-Bozom I, et al. Targeting Wnt/EZH2/microRNA-708 signaling pathway inhibits neuroendocrine differentiation in prostate cancer. Cell Death Discov. 2019;5:139.
Truica CI, Byers S, Gelmann EP. Beta-catenin affects androgen receptor transcriptional activity and ligand specificity. Cancer Res. 2000;60:4709–13.
Yi F, Pereira L, Hoffman JA, Shy BR, Yuen CM, Liu DR, et al. Opposing effects of Tcf3 and Tcf1 control Wnt stimulation of embryonic stem cell self-renewal. Nat Cell Biol. 2011;13:762–70.
Nguyen H, Merrill BJ, Polak L, Nikolova M, Rendl M, Shaver TM, et al. Tcf3 and Tcf4 are essential for long-term homeostasis of skin epithelia. Nat Genet. 2009;41:1068–75.
Hoffman JA, Wu CI, Merrill BJ. Tcf7l1 prepares epiblast cells in the gastrulating mouse embryo for lineage specification. Development. 2013;140:1665–75.
Andoniadou CL, Signore M, Young RM, Gaston-Massuet C, Wilson SW, Fuchs E, et al. HESX1- and TCF3-mediated repression of Wnt/beta-catenin targets is required for normal development of the anterior forebrain. Development. 2011;138:4931–42.
Slyper M, Shahar A, Bar-Ziv A, Granit RZ, Hamburger T, Maly B, et al. Control of breast cancer growth and initiation by the stem cell-associated transcription factor TCF3. Cancer Res. 2012;72:5613–24.
Murphy M, Chatterjee SS, Jain S, Katari M, DasGupta R. TCF7L1 modulates colorectal cancer growth by inhibiting expression of the tumor-suppressor gene EPHB3. Sci Rep. 2016;6:28299.
He Y, Huang H, Farischon C, Li D, Du Z, Zhang K, et al. Combined effects of atorvastatin and aspirin on growth and apoptosis in human prostate cancer cells. Oncol Rep. 2017;37:953–60.
Maynard JP, Ertunc O, Kulac I, Baena-Del Valle JA, De Marzo AM, Sfanos KS. IL8 expression is associated with prostate cancer aggressiveness and androgen receptor loss in primary and metastatic prostate cancer. Mol Cancer Res. 2020;18:153–65.
Jung JH, Kang KW, Kim J, Hong SC, Park Y, Kim BS. CXCR2 inhibition in human pluripotent stem cells induces predominant differentiation to mesoderm and endoderm through repression of mTOR, β-catenin, and hTERT activities. Stem Cells Dev. 2016;25:1006–19.
Li Y, He Y, Butler W, Xu L, Chang Y, Lei K, et al. Targeting cellular heterogeneity with CXCR2 blockade for the treatment of therapy-resistant prostate cancer. Sci Transl Med. 2019;11:eaax0428. https://pubmed.ncbi.nlm.nih.gov/31801883/.
Siu MK, Chen WY, Tsai HY, Chen HY, Yin JJ, Chen CL, et al. TCF7 is suppressed by the androgen receptor via microRNA-1-mediated downregulation and is involved in the development of resistance to androgen deprivation in prostate cancer. Prostate Cancer Prostatic Dis. 2017;20:172–8.
Beltran H, Prandi D, Mosquera JM, Benelli M, Puca L, Cyrta J, et al. Divergent clonal evolution of castration-resistant neuroendocrine prostate cancer. Nat Med. 2016;22:298–305.
Wang G, Jones SJ, Marra MA, Sadar MD. Identification of genes targeted by the androgen and PKA signaling pathways in prostate cancer cells. Oncogene. 2006;25:7311–23.
Nelson PS, Clegg N, Arnold H, Ferguson C, Bonham M, White J, et al. The program of androgen-responsive genes in neoplastic prostate epithelium. Proc Natl Acad Sci USA. 2002;99:11890–5.
Ge R, Wang Z, Montironi R, Jiang Z, Cheng M, Santoni M, et al. Epigenetic modulations and lineage plasticity in advanced prostate cancer. Ann Oncol. 2020;31:470–9.
Li Z, Chen CJ, Wang JK, Hsia E, Li W, Squires J, et al. Neuroendocrine differentiation of prostate cancer. Asian J Androl. 2013;15:328–32.
Conteduca V, Oromendia C, Eng KW, Bareja R, Sigouros M, Molina A, et al. Clinical features of neuroendocrine prostate cancer. Eur J Cancer. 2019;121:7–18.
Sfanos KS, Yegnasubramanian S, Nelson WG, De Marzo AM. The inflammatory microenvironment and microbiome in prostate cancer development. Nat Rev Urol. 2018;15:11–24.
Zhang XQ, Kondrikov D, Yuan TC, Lin FF, Hansen J, Lin MF. Receptor protein tyrosine phosphatase alpha signaling is involved in androgen depletion-induced neuroendocrine differentiation of androgen-sensitive LNCaP human prostate cancer cells. Oncogene. 2003;22:6704–16.
You L, Zhang C, Yarravarapu N, Morlock L, Wang X, Zhang L, et al. Development of a triazole class of highly potent Porcn inhibitors. Bioorg Med Chem Lett. 2016;26:5891–5.
Bonci D, Coppola V, Musumeci M, Addario A, Giuffrida R, Memeo L, et al. The miR-15a-miR-16-1 cluster controls prostate cancer by targeting multiple oncogenic activities. Nat Med. 2008;14:1271–7.
Chesire DR, Isaacs WB. Ligand-dependent inhibition of beta-catenin/TCF signaling by androgen receptor. Oncogene. 2002;21:8453–69.
Wan X, Liu J, Lu JF, Tzelepi V, Yang J, Starbuck MW, et al. Activation of beta-catenin signaling in androgen receptor-negative prostate cancer cells. Clin Cancer Res. 2012;18:726–36.
Murillo-Garzon V, Kypta R. WNT signalling in prostate cancer. Nat Rev Urol. 2017;14:683–96.
Miyamoto DT, Zheng Y, Wittner BS, Lee RJ, Zhu H, Broderick KT, et al. RNA-Seq of single prostate CTCs implicates noncanonical Wnt signaling in antiandrogen resistance. Science. 2015;349:1351–6.
Pal SK, Swami U, Agarwal N. Characterizing the Wnt pathway in advanced prostate cancer: when, why, and how. Eur Urol. 2020;77:22–23.
Zhang Z, Cheng L, Li J, Farah E, Atallah NM, Pascuzzi PE, et al. Inhibition of the Wnt/β-catenin pathway overcomes resistance to enzalutamide in castration-resistant prostate cancer. Cancer Res. 2018;78:3147–62.
Uysal-Onganer P, Kawano Y, Caro M, Walker MM, Diez S, Darrington RS, et al. Wnt-11 promotes neuroendocrine-like differentiation, survival and migration of prostate cancer cells. Mol Cancer. 2010;9:55.
Lescher B, Haenig B, Kispert A. sFRP-2 is a target of the Wnt-4 signaling pathway in the developing metanephric kidney. Dev Dyn. 1998;213:440–51.
MacDonald BT, Tamai K, He X. Wnt/beta-catenin signaling: components, mechanisms, and diseases. Dev Cell. 2009;17:9–26.
Stark T, Livas L, Kyprianou N. Inflammation in prostate cancer progression and therapeutic targeting. Transl Androl Urol. 2015;4:455–63.
Ha H, Debnath B, Neamati N. Role of the CXCL8-CXCR1/2 axis in cancer and inflammatory diseases. Theranostics. 2017;7:1543–88.
Long X, Ye Y, Zhang L, Liu P, Yu W, Wei F, et al. IL-8, a novel messenger to cross-link inflammation and tumor EMT via autocrine and paracrine pathways (Review). Int J Oncol. 2016;48:5–12.
Lim S, Becker A, Zimmer A, Lu J, Buettner R, Kirfel J. SNAI1-mediated epithelial-mesenchymal transition confers chemoresistance and cellular plasticity by regulating genes involved in cell death and stem cell maintenance. PLoS ONE. 2013;8:e66558.
Funding
This work was jointly supported by grants from the Ministry of Science and Technology of Taiwan (MOST109-2326-B-038-001-MY3 and MOST110-2622-B-038-002 to YNL, and MOST109-2314-B-038-105 to YCW), Taipei Medical University-Wan Fang Hospital (109TMU-WFH-01 to WYC), and the National Health Research Institute, Taiwan (NHRI-EX109-10702BI to YNL).
Author information
Authors and Affiliations
Contributions
YCW, YNL, and WYC designed the experiments and supervised the project. YCW, HLY, WHC, KCJ, and SRL performed the experiments. YCW, WYC, and JH provided the human prostate cancer samples. WYC performed histomorphometric analysis. HLY constructed databases and performed statistical and computational analysis. MH assisted to the animal experiments. All the authors analyzed and interpreted the data. YCW, YNL, and WYC wrote, reviewed, and/or revised the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wen, YC., Liu, YN., Yeh, HL. et al. TCF7L1 regulates cytokine response and neuroendocrine differentiation of prostate cancer. Oncogenesis 10, 81 (2021). https://doi.org/10.1038/s41389-021-00371-6
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41389-021-00371-6
This article is cited by
-
Wnt Signaling and Therapeutic Resistance in Castration-Resistant Prostate Cancer
Current Pharmacology Reports (2023)
-
Competing endogenous RNA networks related to prognosis in chronic lymphocytic leukemia: comprehensive analyses and construction of a novel risk score model
Biomarker Research (2022)