Research Article - Journal of Cellular Signaling (2021) Volume 2, Issue 4
Chronic IL-1 Exposed AR+ PCa Cell Lines Show Conserved Loss of IL-1 Sensitivity and Evolve Both Conserved and Unique Differential Gene Expression Profiles
Shayna E. Thomas-Jardin1, Mohammed S. Kanchwala2, Haley Dahl1, Vivian Liu1, Rohan Ahuja1, Reshma Soundharrajan1, Nicole Roos1, Sydney Diep1, Amrit Sandhu1, Chao Xing2,3,4, Nikki A. Delk1*
1Biological Sciences Department, The University of Texas at Dallas, Richardson, TX 75080, USA
2McDermott Center of Human Growth and Development, The University of Texas Southwestern Medical Center, Dallas, TX 75390,USA
3Department of Bioinformatics, The University of Texas Southwestern Medical Center, Dallas, TX 75390, USA
4Department of Population and Data Sciences, The University of Texas Southwestern Medical Center, Dallas, TX 75390, USA
- *Corresponding Author:
- Nikki A. Delk
Received date: August 26, 2021,; Accepted date: October 29, 2021
Citation: Thomas-Jardin SE, Kanchwala MS, Dahl H, Liu V, Ahuja R, Soundharrajan R, et al. Chronic IL-1 Exposed AR+
PCa Cell Lines Show Conserved Loss of IL-1 Sensitivity and Evolve Both Conserved and Unique Differential Gene Expression Profiles. J Cell Signal 2021; 2(4): 248-260.
Copyright: © 2021 Thomas-Jardin SE, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Introduction: Inflammation drives prostate cancer (PCa) progression. While inflammation is a cancer hallmark, the underlying mechanisms mediating inflammation-induced PCa are still under investigation. Interleukin-1 (IL-1) is an inflammatory cytokine that promotes cancer progression, including PCa metastasis and castration resistance. We previously found that acute IL-1 exposure represses PCa androgen receptor (AR) expression concomitant with the upregulation of pro-survival proteins, causing de novo accumulation of castration-resistant PCa cells. However, acute inflammation is primarily anti-tumorigenic, while chronic inflammation is pro-tumorigenic. Thus, using the LNCaP PCa cell line as model, we found that PCa cells can evolve insensitivity to chronic IL-1 exposure, restoring AR and AR activity and acquiring castration resistance. In this paper we expanded our chronic IL-1 model to include the MDA-PCa-2b PCa cell line to investigate the response to acute versus chronic IL-1 exposure and to compare the gene expression patterns that evolve in the LNCaP and MDA-PCa-2b cells chronically exposed to IL-1.
Methods: We chronically exposed MDA-PCa-2b cells to IL-1a or IL-1ß for several months to establish sublines. Once established, we determined subline sensitivity to exogenous IL-1 using cell viability assay, RT-qPCR and western blot. RNA sequencing was performed for parental and subline cells and over representation analysis (ORA) for geneset enrichment of biological process/pathway was performed.
Results: MDA-PCa-2b cells repress AR and AR activity in response to acute IL-1 exposure and evolve insensitivity to chronic IL-1 exposure. While cell biological and molecular response to acute IL-1 signaling is primarily conserved in LNCaP and MDA-PCa-2b cells, including upregulation of NF-?B signaling and downregulation of cell proliferation, the LNCaP and MDA-PCa-2b cells evolve conserved and unique molecular responses to chronic IL-1 signaling that may promote or support tumor progression.
Conclusions: Our chronic IL-1 subline models can be used to identify underlying molecular mechanisms that mediate IL-1-induced PCa progression.
Chronic inflammation, Interleukin-1, IL-1, Prostate cancer, Gene expression
Prostate cancer (PCa) is the second leading cause of cancerrelated death in American men behind lung cancer . Metastatic PCa is responsible for this high burden of mortality, as the five-year survival rate drops from almost 100% for PCa confined to the prostate to 30% for metastatic disease .
Metastatic PCa is currently incurable, and with the annual
burden of metastatic PCa expected to rise 42% by 2025 ,
the need to understand the underlying mechanisms of PCa
disease initiation and progression cannot be overstated.
As part of immunosurveillance, immunes cells infiltrate the tumor microenvironment and secrete cytotoxic and cytostatic inflammatory cytokines [3,4]. Thus, acute inflammation is anti-tumorigenic, but if left unresolved, chronic inflammation ensues [3,4]. Chronic inflammation promotes cancer [3,4], including PCa initiation and progression . While chronic inflammation in the tumor microenvironment is an established hallmark of cancer , the underlying mechanisms of
inflammation-induced disease are not fully elucidated.
Interleukin-1 (IL-1) is an inflammatory cytokine secreted by tumor cells, myeloid precursor cells, macrophages, and neutrophils in the tumor microenvironment [7,8]. IL-1 induces epithelial to mesenchymal transition, metastasis and bone colonization, angiogenesis, and autocrine and paracrine secretion of pro-tumorigenic molecules such as VEGF, IL-6, and IL-8 [9-14]. As such, the FDA approved IL-1 receptor antagonist (IL-1Ra), anakinra, is in clinical trials for breast, pancreatic, and colorectal cancer to block IL-1-induced angiogenesis, metastasis, and chemotherapy resistance (clinicaltrials.gov). In PCa, IL-1 promotes castration resistance, anti-androgen resistance, neuroendocrine differentiation, metastasis, and bone colonization [8-13,15,16]. IL-1 levels are elevated in cancer patients , including in PCa patient serum  and IL-1 levels correlate with advanced Gleason score in primary PCa tumors . We have shown that IL-1Ra can attenuate bone marrow stromal cell paracrine-induced IL-1 signaling in PCa cell lines in vitro . However, we have recently shown that LNCaP PCa cells chronically exposed to IL-1 for several months evolve resistance to exogenous IL-1 and acquire protumorigenic phenotypes  that, as a consequence, can no longer be mitigated by IL-1Ra; thus, emphasizing the need for alternative target molecules and therapies against chronic IL- 1-induced tumorigenicity.
In this study, we compared our RNA sequencing and bioinformatics analysis of LNCaP and MDA-PCa-2b sublines isolated from cultures chronically exposed for several months to IL-1 to identify the conserved and divergent gene expression changes and predicted pathways that are known to promote tumorigenicity, and thus, could underlie chronic IL-1-induced PCa progression and serve as novel therapeutic targets. Both LNCaP and MDA-PCa-2b sublines evolved insensitivity to exogenous IL-1-induced cytotoxicity, cytostaticity, and intracellular signaling and show conserved basal gene expression in a relatively small subset of genes that encode for proteins implicated in cancer, including PCa. LNCaP and MDAPCa- 2b chronic IL-1 sublines share 55 basally high genes and 17 basally low genes. Many of the shared genes reportedly have functional roles in castration resistance and/or correlate with PCa disease progression. Notably, while constitutive NF- κB signaling, a canonical IL-1-induced intracellular pathway, promotes PCa castration resistance [20,21], the LNCaP and
MDA-PCa-2b sublines do not show constitutive upregulation
of canonical NF-κB target genes, suggesting that chronic IL-1
can drive NF-κB-independent PCa disease. Taken together,
chronic IL-1 exposure selects for cytotoxic- and cytostaticresistant
cells that evolve a gene expression profile favoring
Materials and Methods
LNCaP (ATCC, Manassas, VA; CRL-1740), MDA-PCa-2b (ATCC, Manassas, VA; CRL-2422) and subline prostate cancer (PCa) cell lines were maintained in a 37°C, 5.0% (v/v) CO2 growth chamber. LNCaP and LNCaP sublines, LNas1 and LNbs1, were cultured in DMEM (Gibco/Thermo Scientific; 1185-092) supplemented with 10% (v/v) fetal bovine serum essence (FBE; Seradigm; 3100-500), 0.4 mM L-glutamine (L-glut; Gibco/ Invitrogen; 25030-081), and 10 U/ml penicillin G sodium and 10 mg/ml streptomycin sulfate (pen-strep; Gibco/ Invitrogen; 15140-122). MDA-PCa-2b and MDA-PCa-2b sublines, MDA-αs1 and MDA-βs1, were cultured in BRFF-HPC1, (#0403, AthenaES, Baltimore, MD) supplemented with 20% (v/v) FB essence, 0.4mM L-glut, and 10 U/mL Gibco penicillin G sodium and 10 mg/mL streptomycin sulfate (pen-strep; #15140-122, Thermo Fisher Scientific).
Chronic IL-1 subline generation and maintenance
MDA-PCa-2b cells were maintained in HPC1/20% FB Essence (FBE) containing 0.5 ng/ml IL-1α (Gold Bio, St. Louis, MO; 1110-01A-10) or IL-1β (Gold Bio, St. Louis, MO; 1110-01B-10) for approximately 4 months; and the proliferative colonies that emerged were expanded and termed MDA-PCa-2b IL-1α subline (MDA-αs) and MDA-PCa-2b IL-1β subline (MDA-βs). We generated two sets of sublines (MDA- αs1, MDA- αs2, MDA- βs1, MDA- βs2). During subline expansion, the sublines were maintained in 0.5 ng/ml IL-1 up to an additional 4 months. During subline generation and expansion, MDA-PCa-2b parental cells were cultured in vehicle control (0.1% BSA) in 1 x PBS) alongside the sublines. Similar passage number for parental and subline cells were characterized. For phenotypic characterization, established sublines are cultured in the absence of IL-1 for at least one month to ensure the sublines are stable. Cell line authentication was previously performed for LNCaP parental, LNas1, LNbs1 . Cell line authentication was performed for MDA-PCa-2b parental 1, MDA-αs1 and MDA- βs1 was performed by STR profiling by the DNA Genotyping Core, UTSW Medical Center. All cell lines used in the study are commercially available and purchased from American Tissue Culture Collection (ATCC). We hereby confirm that none of the used cell lines require any ethics approval for their use.
Human recombinant IL?1α (GoldBio, St. Louis, MO; 1110- 01A-100) or IL?1β (GoldBio, St. Louis, MO; 1110-01B-100) were resuspended in 0.1% bovine serum albumin (BSA) (Thermo Fisher Scientific, Fair Lawn, NJ; BP 1600?1) in 1X phosphate buffered saline (PBS; Corning, Manassas, VA; 21-040-CM). Cells were treated with vehicle control (0.1% BSA in 1X PBS) or IL?1 added to growth medium.
RNA isolation and reverse transcription quantitative PCR
Total RNA was extracted, reverse transcribed, and analyzed by RT-QPCR as previously described . Primer sequences for genes of interest are listed below. Gene of interest cycle times (CT) were normalized to the β-actin. Relative mRNA levels were calculated using the 2-ΔΔCT method. Primer sequences 5’-3’: Beta actin (β-actin), forward GATGAGATTGGCATGGCT TT, reverse CACCTTCACCGGTCCAGTTT; Mitochondrial Superoxide Dismutase 2 (SOD2), forward GGCCTACGTGAACAACCTGA, reverse GTTCTCCACCACCGTTAGGG; Lipocalin-2 (LCN2), forward TCACCCTCTACGGGAGAACC, reverse GGGACAGGGA AGACGATGTG; NADPH Oxidase 1 (NOX1), forward TGGTTGT TTGGTTAGGGCTGAAT, reverse AGTGTGCGGCTGCAAAACC; Interleukin-1 Receptor Type 1(IL-1R1), forward TGGGGAAGG GTCTACCTCTG, reverse TCCCCAACGTAGTCATCCCT; Follistatin Like 1 (FSTL1), forward CAACCACTGTGAACTGCATCG, reverse GATAGCAAACAACTGGGCTGG; T-cell receptor gamma alternate reading frame protein (TARP), forward GCTGAA ACAAAGCTCCAGAAGG, reverse TCCCAGAATCGTGTTGCTCT; Cancer/Testis Antigen 2 (CTAG2), forward ACATCACGATGCCTT TC TCGTC, reverse AGCAGTCAGTCGGATAAACAGT; Prostaglandin Dehydrogenase 1 (HPGD), forward GCTTAAGGGCGCCAAGGTA,
RNA-sequencing (RNA-seq) analysis
LNCaP cells treated with vehicle control or 25 ng/ml IL-1β  (GSE105088) and untreated LNCaP parental, LNas1 and LNbs1 (GSE142706) were previously sequenced. RNA-seq was performed for MDA-PCa-2b parental 1 cells treated with vehicle control (0.1% BSA in 1xPBS) or 25 ng/ml IL-1α or IL-1β for 3 days and MDA-αs1 and MDA-βs1 grown in vehicle control for 3 days (GSE181229). RNA-seq was performed by the Genome Center at the University of Texas at Dallas (Richardson, TX). Total RNA library was prepared using Illumina Truseq Stranded Total RNA prep Gold kit (Illumina). The prepared libraries were sequenced on an Illumina NextSeq 500 platform (San Diego, CA) with 75bp single-end reads. Fastq files were checked for quality using fastqc (v0.11.2)  and fastq_screen (v0.4.4)  and were then mapped to hg19 (UCSC version from igenomes) using STAR . Read counts were generated using feature Counts  and differential expression analysis was performed using edgeR [26,27]. Differential gene expression
lists were generated using the following cut-offs: log2 counts per million (CPM) ≥ 0, log2 fold change (FC) ≥ 0.6 or ≤ -0.6, false
discovery rate (FDR) ≤ 0.05, adjusted p-value ≤ 0.05. Pathway
analysis was conducted using QIAGEN’s Ingenuity Pathway
Analysis (IPA) tool (http://www.qiagen.com/ingenuity) RNAseq
datasets generated for this study are available at GEO
NCBI, accession GSE181229.
Over representation analysis
Over representation analysis (ORA) was performed using the WEB-based Gene SeT AnaLysis Toolkit (WebGestalt) site (http:/ www.webgestalt.org/) . Gene lists were analyzed using the geneontology biological process functional database. The parameters for the enrichment analysis are minimum/ maximum genes in a category = 5/2000, FDR method = BH, number of categories to report = top 40. Induced and repressed gene lists were uploaded separately to identify, respectively, upregulated and downregulated biological processes.
Protein was isolated from cells using NP40 lysis buffer (0.5% NP40 [US Biological, Salem, MA; N3500], 50 mM of Tris [pH 7.5], 150 mM of NaCl, 3 mM of MgCl2, 1X protease inhibitors [Roche, Mannheim, Germany; 05892953001]). Protein concentration was measured using the Pierce BCA Protein Assay Kit (Thermo Fisher Scientific, Waltham, MA; 23227). For Western blot analysis, equal protein concentrations were loaded onto and separated in 12% (wt/vol) sodium dodecyl sulfate polyacrylamide gel (40% acrylamide/bisacrylamide solution; Bio?Rad, Hercules, CA; 161?0148). Proteins were transferred from the gel to 0.45 μm pore size nitrocellulose membrane (Maine Manufacturing, Sanford, ME; 1215471) and total protein visualized using Ponceau S (Amresco, Radnor, PA; K793). The membrane was blocked with 2.5% (wt/vol) BSA (Thermo Fisher Scientific, Waltham, MA; BP 1600?1) in 1X trisbuffered saline with Tween 20 (TBST; 20 mM of Tris, pH 7.6, 150 mM of NaCl, 0.05% Tween?20). Primary and secondary antibodies were diluted in 2.5% BSA in 1X TBST. Protein blot bands were visualized using Clarity Western ECL Substrate (Bio?Rad, Hercules, CA; 1705061) or SuperSignal West Femto Maximum Sensitivity Substrate (Thermo Scientific, Rockford, IL; 34095) and imaged using Amersham Imager 600 (GE, Marlborough, MA). Primary antibodies: AR (Cell Signaling, Danvers, MA; D6F11), β-actin (Santa Cruz, Santa Cruz, CA; sc- 69879), LCN2 (Cell Signaling, Danvers, MA; D4M8L), NKX3.1(Cell Signaling, D2Y1A), PARP (Cell Signaling, Danvers, MA; 9532S), PSA (Cell Signaling, Danvers, MA; 5365S), SOD2 (Abgent, San Diego, CA; AM7579a). Secondary antibodies: Sheep anti-mouse (Jackson ImmunoResearch Laboratories; 515-035-062), goat anti-rabbit (Abnova; PAB10822).
bromide) viability assay
MTT assay (Trevigen, Gaithersburg, MD; 4890-25-K) was performed according to manufacturer’s instructions. Cell viability was quantified as the optical density (OD) read at wavelengths 540 nm and 650 nm. The final OD was calculated as follows: OD 540 nm – OD 650 nm. OD was measured using the Cytation3 Imaging Reader (BioTek, Winooski, VT).
Statistical significance for MTT and cell counts was determined using unpaired student t-test calculated using Microsoft Excel. P-values of ≤ 0.05 were considered statistically significant and denoted by asterisks (*p ≤ 0.05; **p ≤ 0.05; ***p ≤ 0.005). Graphs are shown as the average of a minimum of n = 3 biological replicates ± standard deviation (STDEV).
MDA-PCa-2b cells evolve insensitivity to chronic IL-1
To investigate the effect of chronic IL-1 exposure on MDAPCa- 2b cells, cells were cultured in 0.5 ng/ml IL-1α or IL-1β for 3-4 months, as previously described for LNCaP cells . IL-1α and IL-1β (hereinafter, collectively referred to as IL-1) are the two main IL-1 cytokine family members, both bind the IL-1 receptor and have similar biological activity . As with LNCaP cells [16,29], IL-1 is cytotoxic (Figure 1A) and cytostatic (Figure 1B) for MDA-PCa-2b cells; but after 3-4 months cultured in IL-1, proliferative colonies emerged and were expanded to generate sublines. For rigor, we generated two sets of MDAPCa- 2b chronic IL-1 sublines that we termed MDA-PCa-2b IL-1α subline 1 (MDA-αs1), MDA-PCa-2b IL-1α subline 2 (MDA- αs2), MDA-PCa-2b IL-1β subline 1 (MDA-βs1), and MDA-PCa-
2b IL-1β subline 2 (MDA-βs2).
To demonstrate that the MDA-PCa-2b chronic IL-1 sublines acquired insensitivity to IL-1, the parental and subline cells were given an acute IL-1 exposure of 25 ng/ml IL-1 for 3 or 6 days. Acute IL-1 exposure reduced cell viability (Figure 1A) and proliferation (Figure 1B), reduced full-length PARP (indicative of cell death activation), and upregulated the protein accumulation of the canonical IL-1-induced genes, SuperOxide Dismutase (SOD2) and Lipocalin 2 (LCN2) (Figure 1C).
Acute IL-1 also represses AR levels and activity in AR+ PCa cell lines [19,30]. As such, acute IL-1 treatment repressed AR and
NKX3.1 (canonical AR-induced gene) protein accumulation
in the parental cells, but not in the subline cells (Figure 1C).
Finally, while acute IL-1 exposure represses AR accumulation
and activity in MDA-PCa-2b parental cells, as observed for
LNCaP , the chronically exposed subline cells have restored
AR accumulation and activity (Figure 1C). Taken together,
these results demonstrate that chronic IL-1 exposure selects
for MDA-PCa-2b cells that evolve insensitivity to exogenous
IL-1. Of note, mRNA levels of the IL-1 receptor, IL-1R1, are not
induced by IL-1 in parental cells (Figure 2A), indicating that basal IL-1R1 levels are sufficient to mediate IL-1 intracellular signaling in parental cells. Furthermore, IL-1 does not show a differential effect on IL-1R1 mRNA levels in parental versus subline cells (Figure 2A) and IL-1R1 mRNA levels are comparable in parental versus subline cells (Figure 2B), suggesting subline insensitivity to exogenous IL-1 is independent of IL-1R1 levels.
Over representation analysis (ORA) of differentially expressed genes predicts acute and chronic IL-1 response in the LNCaP and MDA-PCa-2b cells
RNA sequencing was previously performed for LNCaP cells treated acutely for 3 days with 25 ng/ml IL-1β  (GSE105088) and, here, we performed RNA sequencing for MDA-PCa- 2b cells treated acutely for 3 days with 25 ng/ml IL-1α or IL- 1β (GSE181229). LNCaP and MDA-PCa-2b share 927 genes upregulated by acute IL-1 exposure and share 747 repressed genes (Supplementary Table 1; Supplementary Figure 1A). We performed gene set ORA against the Gene Ontology biological process database using WebGestalt functional enrichment analysis web tool . IL-1 is a cytotoxic and cytostatic inflammatory cytokine that signals through NF-κB [16,29,31,32]; thus, as expected, ORA of the shared upregulated genes predicts that cytokine and immune response pathways, including NF-κB signaling, are activated by acute IL-1 exposure in LNCaP and MDA-PCa-2b cells and ORA of the shared repressed genes predicts down regulation of the cell cycle and DNA, RNA, and protein biosynthesis (Supplementary Table 1; Supplemental Figure 1B). The ORA results are in line with our experimental data showing the LNCaP and MDA-PCa-2b cells upregulate IL-1-NF-κB induced genes [33-37] LCN2, NOX1, and SOD2 (Figure 3) and show reduced viability and proliferation
in response to IL-1 (Figure 1).
We also previously performed RNA sequencing for basal gene expression in the LNCaP chronic IL-1 sublines  (GSE142706) and, here, we performed RNA sequencing for basal gene expression in the MDA-PCa-2b chronic IL-1 sublines (GSE181229). We performed ORA of the suite of basally high genes in LNCaP chronic IL-1 sublines and in MDAPCa- 2b chronic IL-1 sublines to predict upregulated biological processes (Supplementary Table 1, Tab 7 and Tab 10). In line with our experimental data showing that the LNCaP and MDA-PCa-2b chronic IL-1 sublines do not have constitutively elevated mRNA levels of NFkB-induced genes (Figure 2B), NF- κB signaling is not represented in the biological pathways predicted by ORA of the high basal gene suites (Supplementary Table 1, Tab 7 and Tab 10). ORA of basally low genes in LNCaP chronic IL-1 sublines or in MDA-PCa-2b chronic IL-1 sublines predicts that inflammatory responses are repressed, as well as the negative regulation of cell proliferation and growth (Supplementary Table 1, Tab 8 and Tab 11). Taken together, our ORA and experimental results indicate the LNCaP and MDA-PCa-2b parental cells show conserved response to acute IL-1, activating IL-1-NF-κB intracellular signaling and showing cytostatic and cytotoxic cell responses. Furthermore, following chronic IL-1 exposure, both LNCaP and MDA-PCa-2b cells lose canonical IL-1-NF-κB signaling and negative regulation of cell proliferation and growth.
Conserved and unique gene expression patterns in
LNCaP and MDA-PCa-2b chronic IL-1 sublines
Given that IL-1 promotes cancer progression , we compared the basal differential gene expression of LNCaP  and MDA-PCa-2b chronic IL-1 subline cells (Supplementary Table 1) to identify predicted tumorigenic pathways. Compared to LNCaP parental cells, there are 2954 basally high genes in LNas1 and LNbs1 chronic IL-1 sublines and 629 basally low genes (Supplementary Table 1; Supplementary Figure 2A). Compared to MDA-PCa-2b parental cells, there are 204 basally high genes in MDA-αs1 and MDA-βs1 chronic IL-1 sublines and 184 basally low genes (Supplementary Table 1; Supplementary Figure 2A). Of the differentially expressed genes in each subline set, the LNCaP and MDAPCa- 2b chronic IL-1 sublines share expression patterns in a select subset of genes. LNas1, LNbs1, MDA-αs1, and MDA- βs1 sublines share 55 basally high genes and 17 basally low genes (Supplementary Table 1; Supplementary Figure 2A). We confirmed expression of several basally high genes in the LNCaP and MDA-PCa-2b chronic IL-1 sublines that are also associated with PCa progression, CTAG2, FSTL1, HPGD and TARP (Figure 4) [38-45]. Of note, among the independently
generated MDA-PCa-2b chronic IL-1 sublines, the basal mRNA
levels of the select confirmed genes is variable. We believe this is due to clonal heterogeneity and genetic drift during
IL-1 chronic subculturing. Nevertheless, we see conserved
basally high mRNA levels for the selected genes which we will
functionally investigate in future studies.
We performed ORA for the 55 basally high and 17 basally low genes which predicted several constitutively upregulated or repressed pathways in both the LNCaP and MDA-PCa-2b chronic IL-1 sublines (Supplementary Table 1; Supplementary Figure 2B). For example, ORA of the basally high gene set predicts that the chronic IL-1 sublines evolved constitutive expression of genes that regulate insulin signaling and control glucose homeostasis and ORA of the basally low gene set predicts that the chronic IL-1 sublines evolved constitutive repression of genes that regulate abiotic stress response (Supplementary Table 1; Supplementary Figure 2). Notably, the genes involved in regulating insulin signaling and glucose homeostasis, DPP4, FKBP1B, ICAM1, IGF1, IGF1R, INHBB, PTPRN2, STXBP5L, and VWA2, and the genes that regulate abiotic stress response, HPN, DNAJB1, FOS, IER5, HSPA6, and ARNT2 (Supplementary Table 1, Tab 13 and Tab 14), are also correlated with and/or promote anti- or pro-tumorigenic phenotypes, including in PCa.
In addition to the shared gene expression pattern, the LNCaP and MDA-PCa-2b chronic IL-1 sublines show cell line-unique changes in gene expression that were acquired in response to chronic IL-1 exposure (Supplementary Table 1). We performed ORA for basal gene expression in LNas1 and LNbs1 (Supplementary Table 1, Tab 7 and Tab 8) and performed ORA for basal gene expression in MDA-αs1 and MDA-βs1 (Supplementary Table 1, Tab 10 and Tab 11). As expected, ORA predicts that LNCaP and MDA-PCa-2b chronic IL-1 sublines evolved distinct biological processes. For example, ORA predicts that cell projection and neurogenesis biological processes are upregulated and cell adhesion is downregulated in the LNas1 and LNbs1 sublines (Supplementary Table 1); while in the MDA-as1 and MDA-bs1 sublines, interferon gamma signaling is predicted to be activated and VEGF signaling is predicted to be downregulated (Supplementary Table 1). Thus, while LNCaP and MDA-PCa-2b chronic IL-1 sublines share conserved insensitivity to IL-1-induced intracellular signaling  (Figures 1 and 3), the broader response to chronic IL-1 exposure is unique among the cell lines. The divergent response to chronic IL-1 likely reflects genetic differences in LNCaP and MDA-PCa-2b cell lines, as could be expected for tumor cell response to chronic inflammation among different patients. Taken together, our models reveal potential therapeutic targets for individuals whose PCa disease progression evolves from chronic inflammation.
We have developed novel models to explore the effect of chronic inflammation on PCa progression. Specifically, we generated LNCaP  and MDA-PCa-2b sublines from cultures chronically exposed to IL-1 for several months.
Given that IL-1 is cytotoxic and cytostatic [16,29] (Figure 1),
chronic IL-1 exposure selected for viable, proliferative PCa
cells that evolved IL-1 insensitivity. Interestingly, we find that
castration-resistant cell lines have reduced or no sensitivity
to IL-1-induced cytotoxicity and cytostaticity [16,29] (data
not shown) and we and other labs have found that chronic
IL-1 exposure can select for castration-resistant PCa cells
[15,16]. Chronic IL-1 sublines can also develop resistance to
other cytotoxic and cytostatic inflammatory cytokines .
Therefore, we speculate that chronic IL-1 exposure in the
PCa tumor microenvironment imposes selective pressure on
tumor cells to acquire constitutive pro-survival mechanisms
that provide a growth advantage against cellular stresses, such
as castration or acute cytostatic and cytotoxic inflammation.
Thus, we can use our subline models to elucidate strategies
that PCa cells elicit to escape inflammatory cytokine toxicity
and growth inhibition.
Canonical IL-1 signaling occurs through the NF-κB transcription factor  and constitutive NF-κB activity induces PCa castration resistance [20,21]. Therefore, interventions targeting IL-1/IL-1R1 binding  or inhibiting NF-κB intracellular pathway activation  could mitigate NF-κB-dependent PCa castration resistance. But for cells that evolve resistance to IL-1 due to chronic IL-1 exposure and, likewise, do not acquire constitutive NF-κB signaling, IL-1-NF- κB-targeting strategies would be ineffective. Therefore, our chronic IL-1 cell line models can be used to identify alternative therapeutic targets.
While both LNCaP and MDA-PCa-2b cell lines evolve insensitivity to chronic IL-1 exposure, the underlying mechanism of insensitivity remains unknown. LNCaP chronic IL-1 sublines also develop insensitivity to TNFα  which also signals through NF-κB, suggesting that the resistance mechanism is downstream of the IL-1/IL-1R1 receptor/ligand interaction in LNCaP cells. IL-1R1 mRNA levels are comparable in LNCaP parental and sublines cells  and comparable in MDA-PCa-2b parental and subline cells (Figure 2), implying that the IL-1 resistance mechanism is independent of IL-1R1 levels. However, it will be important to determine efficacy of the IL-1/IL-1R1 protein-protein interaction in MDA-PCa-2b cells and given the different genetic backgrounds of LNCaP and MDA-PCa-2b cell lines, it is possible that the mechanism(s) governing IL-1 insensitivity in the LNCaP versus MDA-PCa- 2b chronic IL-1 sublines are not conserved. We will use our subline models to determine how cells can acquire resistance to chronic IL-1 exposure. Uncovering the mechanisms may help mitigate chronic IL-1-induced castration resistance or PCa progression.
The LNCaP and MDA-PCa-2b chronic IL-1 sublines share a select basal gene expression profile of 55 basally high and 17 basally low genes (Supplementary Table 1). Significantly, many of these genes have reported roles in cancer and may be critical for inflammation-induced PCa progression. For example, the shared basally high genes that we confirmed by RT-qPCR, CTAG2, FSTL1, HPGD and TARP are pro-tumorigenic. CTAG2 (cancer-testis antigen-2) encodes for the autoimmunogenic tumor antigen, LAGE-1, expressed by cancer cells, including PCa . Functionally, CTAG2 was shown to promote breast
cancer cell invasion . Phase 1 clinical trial results are
pending for a combined NY-ESO-1/LAGE-1 vaccine for patients
with androgen-independent metastatic PCa (NCT00616291).
Metastatic castration resistant PCa (mCRPC) patients give an
NY-ESO-1 (CTAG1) peptide vaccine in a Phase 1 clinical trial
showed slowed PSA doubling time , thus, the combined NY-ESO-1/LAGE-1 vaccine is expected to be efficacious. FSTL1 (Follistatin-like protein 1) is a proinflammatory glycoprotein  that has been shown to promote invasion and metastasis
in colorectal cancer  and chemoresistance and stemness
in breast cancer . FSTL1 is an AR-induced gene and is elevated in CRPC . HPGD (15-Hydroxyprostaglandin dehydrogenase) is an anti-inflammatory enzyme that catabolizes prostaglandins . HPGD is an AR-induced gene  and is elevated in primary and metastatic PCa tissue . HPGD siRNA silencing reduced LNCaP cell viability in vitro . TARP (T-cell receptor gamma alternate reading frame protein) is an autoimmunogenic tumor antigen expressed by PCa cells . TARP is an AR-induced gene  and TARP overexpression in the PC3 PCa cell line reduced IL-1b expression and enhanced cell proliferation . A Phase 1 clinical trial showed
that TARP peptide vaccination slows PSA doubling time and
tumor growth rate in Stage D0 PCa patients with biochemical
PSA recurrence (NCT00972309) . Taken together, the
conserved molecular response to chronic IL-1 exposure may
have functional significance in PCa tumorigenicity which we
can elucidate using our chronic IL-1 cell line models.
Indeed, the individual genes modulated by chronic IL-1 exposure that reportedly function in PCa disease can be explored as novel therapeutic targets, but it will also be important to interrogate the signaling networks that integrate the global changes in chronic IL-1-induced gene expression, including those that are unique to the different cell backgrounds. Our ORA bioinformatic analysis provides putative tumor promoting pathways to investigate that would be acquired during chronic IL-1 exposure, such as insulin and interferon gamma signaling, glucose homeostasis and metabolism and cell adhesion.
The basal changes in gene expression in parental versus chronic IL-1 subline cells suggests that chronic IL-1 exposure drives constitutive genetic or epigenetic changes over time. LNCaP and MDA-PCa-2b intracellular signaling, gene expression patterns, ORA-predicted cell biological processes, and cell behavior in response to acute IL-1 exposure are conserved, but the gene expression patterns and ORApredicted cell biological processes diverge as the different cell lines evolve insensitivity to chronic IL-1 signaling. This suggests that genetic and epigenetic reprogramming in response to chronic IL-1 exposure would evolve differently in different patients. Therefore, we can begin to uncover the various ways in which chronic IL-1 exposure remodels chromatin using our chronic IL-1 subline models in the different cell line backgrounds. Studies are underway to elucidate chronic IL-1- induced chromatin remodeling, the subsequent constitutive changes in signaling networks, and the contribution of these altered signaling networks to PCa progression.
As stated, while LNCaP and MDA-PCa-2b cell lines both develop insensitivity to chronic IL-1 exposure, the two cell line backgrounds evolved primarily disparate changes in gene expression, likely due to underlying genetic differences in the cell line backgrounds. Notably, however, chronic exposure to IL-1α or IL-1β both cause similar gene expression changes within a cell line background, indicating that a conserved mechanism leads to chronic IL-1 insensitivity and changes in intracellular signaling within the cell line. Therefore, we can identify master regulators of the chronic IL-1 molecular response downstream of the IL-1 ligand/receptor interaction using our subline models.
Although the LNCaP and MDA-PCa-2b chronic IL-1 sublines have lost sensitivity to IL-1-induced intracellular signaling, the sublines are still sensitive to other stimuli such as IL-6 or serum starvation (data not shown). Thus, the context of the tumor microenvironment and conventional therapy will be very important in determining the cell biological and physiological consequence of IL-1-insensitive PCa cell subpopulations in PCa progression. Those studies are underway.
Taken together, we propose the following model (Figure 5) wherein acute IL-1 exposure is anti-tumorigenic, reducing cell viability and proliferation and repressing AR and AR activity in androgen dependent PCa cells. If the acute inflammation is left unresolved, the PCa cell subpopulation that survives the initial acute IL-1 assault evolves insensitivity to chronic IL-1 exposure and emerges more fit to contribute to PCa tumor growth and progression. It remains unknow if the conserved loss of sensitivity to chronic IL-1 exposure among PCa cell lines is due to the same or different intracellular mechanism; but the culmination of epigenetic and gene expression changes that result from chronic IL-1 exposure appear to diverge among PCa cell lines and likely reflects the inherent genomic difference among cell lines. To this end, while AR and AR activity are restored in both LNCaP  and MDA-PCa-2b chronic IL-1 sublines (Figure 1), subsequent AR dependency and therapeutic response may differ based on the differing epigenetic and gene expression landscapes that evolve. Thus, we can take advantage of our current cell line models and expand them to patient derived tumor models to understand how chronic IL-1 inflammation can affect PCa disease and treatment in individuals.
We would like to thank the Delk and Xing lab members and
funding sources for support.
Conflict of Interests
The authors declare no potential conflicts of interest
Shayna E. Thomas-Jardin, led execution of experiments and analysis, prepared and edited manuscript; Mohammed S. Kanchwala, led bioinformatics analysis, prepared and edited manuscript; Haley Dahl assisted with cell line generation, cell biology analysis, protein and RNA analysis; Vivian Liu, Rohan Ahuja, Reshma Soundharrajan, Nicole Roos assisted with protein and RNA analysis; Sydney Diep, Amrit Sandhu assisted with bioinformatics analysis; Chao Xing oversaw bioinformatics analysis; Nikki A. Delk, led conceptualization and analysis of all data, prepared and edited manuscript.
Cell lines were authenticated. The content is solely the responsibility of the authors and does not necessarily represent the official views of the University of Texas at Dallas or the National Institutes of Health.
Data Availability Statement
The data and generated cell lines that support the findings of this study are available on request from the corresponding author. RNA sequencing data is publicly available at GEO NCBI, accession GSE181229.
The project was supported by funding from The University of Texas at Dallas (Delk) and by the Geographic Management of Cancer Health Disparities Program (GMaP) Region 3, a Supplement of the National Institutes of Health grant P30 CA118100 (Willman, CL). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
- Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J
- Kelly SP, Anderson WF, Rosenberg PS, Cook MB. Past, Current, and
Future Incidence Rates and Burden of Metastatic Prostate Cancer in
the United States. Eur Urol Focus. 2018;4:121-27.
- Zhao H, Wu L, Yan G, Chen Y, Zhou M, Wu Y, et al. Inflammation and
tumor progression: signaling pathways and targeted intervention.
Signal Transduct Target Ther. 2021;6:263.
- Korniluk A, Koper O, Kemona H, Dymicka-Piekarska V. From
inflammation to cancer. Irish J Med Sci (1971 -). 2017;186:57-62.
- Cai T, Santi R, Tamanini I, Galli IC, Perletti G, Bjerklund Johansen TE,
et al. Current Knowledge of the Potential Links between Inflammation
and Prostate Cancer. Int J Mol Sci. 2019;20:3833.
- Hanahan D, Weinberg R a. Hallmarks of cancer: the next
generation. Cell. 2011;144:646-74.
- Voronov E, Carmi Y, Apte RN. The role IL-1 in tumor-mediated
angiogenesis. Front Physiol. 2014;5:114.
- Yu Y, Zhang Q, Ma C, Yang X, Lin R, Zhang H, et al. Mesenchymal
stem cells recruited by castration-induced inflammation activation
accelerate prostate cancer hormone resistance via chemokine ligand
5 secretion. Stem Cell Res Ther. 2018;9:1-12.
- Chiao J, Hsieh T, Xu W, Sklarew R, Kancherla R. Development
of human prostate cancer cells to neuroendocrine-like cells by
interleukin-1. Int J Oncol. 1999;15:1033-40.
- Diaz M, Abdul M, Hoosein N. Modulation of neuroendocrine
differentiation in prostate cancer by interleukin-1 and -2. Prostate.
- Li Y, Wang L, Pappan L, Galliher-Beckley A, Shi J. IL-1ß promotes
stemness and invasiveness of colon cancer cells through Zeb1
activation. Mol Cancer. 2012;11:87.
- Liu Q, Russell MR, Shahriari K, Jernigan DL, Lioni MI, Garcia FU,
et al. Interleukin-1ß promotes skeletal colonization and progression
of metastatic prostate cancer cells with neuroendocrine features.
Cancer Res. 2013;73:3297-305.
- Shahriari K, Shen F, Worrede-Mahdi A, Liu Q, Gong Y, Garcia FU,
et al. Cooperation among heterogeneous prostate cancer cells in the
bone metastatic niche. Oncogene. 2016/12/19. 2017;36:2846-56.
- Perrier S, Caldefie-Chézet F, Vasson MP. IL-1 family in breast
cancer: Potential interplay with leptin and other adipocytokines.
FEBS Letters. 2009;83(2):259-65.
- Staverosky JA, Zhu X, Ha S, Logan SK. Anti-androgen resistance
in prostate cancer cells chronically induced by interleukin-1ß. Am J
Clin Exp Urol. 2013;1:53-65.
- Dahl HC, Kanchwala M, Thomas-Jardin SE, Sandhu A, Kanumuri
P, Nawas AF, et al. Chronic IL-1 exposure drives LNCaP cells to
evolve androgen and AR independence. Batra SK, editor. PLoS One.
- Litmanovich A, Khazim K, Cohen I. The Role of Interleukin-1 in
the Pathogenesis of Cancer and its Potential as a Therapeutic Target
in Clinical Practice. Oncol Ther. 2018;6:109-27.
- Tazaki E, Shimizu N, Tanaka R, Yoshizumi M, Kamma H, Imoto S, et
al. Serum cytokine profiles in patients with prostate carcinoma. Exp
Ther Med. 2011;2:887-91.
- Thomas-Jardin SE, Kanchwala MS, Jacob J, Merchant S, Meade RK,
Gahnim NM, et al. Identification of an IL-1-induced gene expression pattern in AR+ PCa cells that mimics the molecular phenotype of ARPCa
cells. Prostate. 2018;78:595-606.
- Jin RJ, Lho Y, Connelly L, Wang Y, Yu X, Jean L Saint, et al. The
nuclear factor-kappaB pathway controls the progression of prostate
cancer to androgen-independent growth. Cancer Res. 2008;68:6762-9.
- Jeong JH, Park SJ, Dickinson SI, Luo JL. A Constitutive Intrinsic
Inflammatory Signaling Circuit Composed of miR-196b, Meis2,
PPP3CC, and p65 Drives Prostate Cancer Castration Resistance. Mol
- Andrews S. FastQC: a quality control tool for high throughput
sequence data. 2010. Available: http://www.bioinformatics.
- Wingett S. FastQ Screen: quality control tool to screen a library of
sequences in FastQ format against a set of sequence databases. 2011.
- Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al.
STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15-
- Liao Y, Smyth GK, Shi W. FeatureCounts: An efficient general
purpose program for assigning sequence reads to genomic features.
- R Core Development Team. R: A language and environment for
statistical computing. Vienna, Austria. 2019.
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: A Bioconductor
package for differential expression analysis of digital gene expression
data. Bioinformatics. 2009;26:139-40.
- Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene
set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res.
- Thomas-Jardin SE, Dahl H, Kanchwala MS, Ha F, Jacob J,
Soundharrajan R, et al. RELA is sufficient to mediate interleukin-1
repression of androgen receptor expression and activity in an LNCaP
disease progression model. Prostate. 2020;80:133-45.
- Chang MA, Patel V, Gwede M, Morgado M, Tomasevich K, Fong
ELL, et al. IL-1ß Induces p62/SQSTM1 and Represses Androgen
Receptor Expression in Prostate Cancer Cells. J Cell Biochem.
- Oeckinghaus A, Ghosh S. The NF-kappaB family of transcription
factors and its regulation. Cold Spring Harb Perspect Biol.
- Hiscott J, Marois J, Garoufalis J, D’Addario M, Roulston A, Kwan I,
et al. Characterization of a functional NF-kappa B site in the human
interleukin 1 beta promoter: evidence for a positive autoregulatory
loop. Mol Cell Biol. 1993;13:6231-40.
- Tesoriere L, Attanzio A, Allegra M, Gentile C, Livrea MA.
Indicaxanthin inhibits NADPH oxidase (NOX)-1 activation and NF kB-dependent release of inflammatory mediators and prevents the
increase of epithelial permeability in IL-1b-exposed Caco-2 cells. Br J
Nutr . 2014;111(3):415-23.
- Manea A, Tanase LI, Raicu M, Simionescu M. Transcriptional
regulation of NADPH oxidase isoforms, Nox1 and Nox4, by nuclear
factor-kappaB in human aortic smooth muscle cells. Biochem Biophys
Res Commun. 2010;396:901-7.
- Zhao P, Stephens JM. STAT1, NF-?B and ERKs play a role in
the induction of lipocalin-2 expression in adipocytes. Mol Metab.
- Sommer G, Weise S, Kralisch S, Lossner U, Bluher M, Stumvoll M,
et al. Lipocalin-2 is induced by interleukin-1beta in murine adipocytes
in vitro. J Cell Biochem. 2009;106:103-8.
- Xu Y, Kiningham KK, Devalaraja MN, Yeh C-C, Majima H, Kasarskis
EJ, et al. An Intronic NF-kappaB Element Is Essential for Induction of the
Human Manganese Superoxide Dismutase Gene by Tumor Necrosis
Factor-alpha and Interleukin-1beta. DNA Cell Biol. 1999;18:709-22.
- Sonpavde G, Wang M, Peterson LE, Wang HY, Joe T, Mims MP, et
al. HLA-restricted NY-ESO-1 peptide immunotherapy for metastatic
castration resistant prostate cancer. Invest New Drugs. 2014;32:235–
- Su S, Parris AB, Grossman G, Mohler JL, Wang Z, Wilson EM.
Up-Regulation of Follistatin-Like 1 By the Androgen Receptor and
Melanoma Antigen-A11 in Prostate Cancer. Prostate. 2017;77:505-16.
- Kim T-H, Park J-M, Kim M-Y, Ahn Y-H. The role of CREB3L4 in the
proliferation of prostate cancer cells. Sci Rep . 2017;7:45300.
- Vainio P, Gupta S, Ketola K, Mirtti T, Mpindi JP, Kohonen P, et
al. Arachidonic acid pathway members PLA2G7, HPGD, EPHX2, and
CYP4F8 identified as putative novel therapeutic targets in prostate
cancer. Am J Pathol. 2011;178:525-36.
- Epel M, Carmi I, Soueid-Baumgarten S, Oh SK, Bera T, Pastan I,
et al. Targeting TARP, a novel breast and prostate tumor-associated
antigen, with T cell receptor-like human recombinant antibodies. Eur
J Immunol. 2008;38(6):1706-20.
- Cheng W-S, Giandomenico V, Pastan I, Essand M. Characterization
of the Androgen-Regulated Prostate-Specific T Cell Receptor-Chain
Alternate Reading Frame Protein (TARP) Promoter. 2003;144(8):3433-
- Wolfgang CD, Essand M, Lee B, Pastan I. T-cell receptor gamma
chain alternate reading frame protein (TARP) expression in prostate
cancer cells leads to an increased growth rate and induction of
caveolins and amphiregulin. Cancer Res. 2001;61:8122-6.
- Wood L V., Fojo A, Roberson BD, Hughes MSBB, Dahut W, Gulley
JL, et al. TARP vaccination is associated with slowing in PSA velocity
and decreasing tumor growth rates in patients with Stage D0 prostate
cancer. Oncoimmunology. 2016;5:1-13.
- Thomas-Jardin S, Dahl H, Nawas AF, Bautista M, Delk NA. NF-kB
Signaling Promotes Castration-Resistant Prostate Cancer Initiation
and Progression. Pharmacol Ther. 2020;211:107538.
- Lethé B, Lucas S, Michaux L, De Smet C, Godelaine D, Serrano
A, et al. LAGE-1, a new gene with tumor specificity. Int J cancer.
- Maine EA, Westcott JM, Prechtl AM, Dang TT, Whitehurst AW,
Pearson GW. The cancer-testis antigens SPANX-A/C/D and CTAG2
promote breast cancer invasion. Oncotarget. 2016;7:14708-26.
- Chaly Y, Hostager B, Smith S, Hirsch R. Follistatin-like protein 1
and its role in inflammation and inflammatory diseases. Immunol
- Gu C, Wang X, Long T, Wang X, Zhong Y, Ma Y, et al. FSTL1 interacts
with VIM and promotes colorectal cancer metastasis via activating
the focal adhesion signalling pathway. Cell Death Dis . 2018;9(6):654.
- Cheng S, Huang Y, Lou C, He Y, Zhang Y, Zhang Q. FSTL1 enhances
chemoresistance and maintains stemness in breast cancer cells via
integrin ß3/Wnt signaling under miR-137 regulation. 2019;20(3):328-
- Schmidleithner L, Thabet Y, Schönfeld E, Köhne M, Sommer D,
Abdullah Z, et al. Enzymatic Activity of HPGD in Treg Cells Suppresses
Tconv Cells to Maintain Adipose Tissue Homeostasis and Prevent
Metabolic Dysfunction. Immunity. 2019;50:1232-48.e14.