The CIC-ERF co-deletion underlies fusion-independent activation of ETS family member, ETV1, to drive prostate cancer progression
Abstract
Human prostate cancer can result from chromosomal rearrangements that lead to aberrant ETS gene expression. The mechanisms that lead to fusion-independent ETS factor upregulation and prostate oncogenesis remain relatively unknown. Here, we show that two neighboring transcription factors, Capicua (CIC) and ETS2 repressor factor (ERF), which are co-deleted in human prostate tumors can drive prostate oncogenesis. Concurrent CIC and ERF loss commonly occur through focal genomic deletions at chromosome 19q13.2. Mechanistically, CIC and ERF co-bind the proximal regulatory element and mutually repress the ETS transcription factor, ETV1. Targeting ETV1 in CIC and ERF-deficient prostate cancer limits tumor growth. Thus, we have uncovered a fusion-independent mode of ETS transcriptional activation defined by concurrent loss of CIC and ERF.
Editor's evaluation
This study provides insight into a potentially new genetically defined subset of prostate tumors driven by concurrent loss of the ERF and CIC tumor suppressor genes, in the absence of the canonical fusion event involving TMPRSS2 (around 10% of all cases). The work both validates previous findings and provides new data that support a compelling overall conclusion that combined ERF and CIC loss promotes prostate tumorigenesis by increasing expression of the oncogenic driver ETV1. This is an important study based on convincing evidence, that will be of interest to researchers in the field of prostate cancer.
https://doi.org/10.7554/eLife.77072.sa0Introduction
Prostate cancer (PCa) is the most common solid tumor malignancy in men. Activation of ETS transcription factors (TFs), ERG, ETV1, ETV4, and ETV5, are present in approximately 60% of PCa, underscoring their importance in prostate oncogenesis (Nelson et al., 2003). In human PCa, ETS TFs are most commonly activated through gene rearrangements that fuse the androgen responsive gene, TMPRSS2, to either ERG, ETV1, ETV4, or ETV5 (Sizemore et al., 2017). Beyond ETS TF fusions, little is known about the underlying molecular mechanisms that lead to increased expression of wildtype (WT) ETS factors, which confer aggressive malignant phenotypes and associate with poor clinical outcomes in fusion negative PCa patients (Baena et al., 2013).
Capicua (CIC) is a High-mobility group (HMG) box TF that silences ETV1, ETV4, and ETV5 through direct target gene repression (Kim et al., 2021; Simón-Carrasco et al., 2018). CIC is frequently altered in human cancer, where it functionally suppresses tumor growth and metastasis (Ahmad et al., 2019; Bettegowda et al., 2011; Choi et al., 2015; Dissanayake et al., 2011; Kim et al., 2018; Okimoto et al., 2017; Seim et al., 2017; Yoshiya et al., 2021). Notably, in PCa, CIC is commonly altered through genomic loss (homozygous and heterozygous deletion) in ~10% of PCa patients (Abida et al., 2019; Grasso et al., 2012; Hieronymus et al., 2014; Huang et al., 2017; Robinson et al., 2015; Cancer Genome Atlas Research Network, 2015) and inactivation of CIC de-represses ETV1, ETV4, and ETV5 transcription to promote tumor progression (Bettegowda et al., 2011; Choi et al., 2015; Kim et al., 2018; Okimoto et al., 2017). Leveraging mutational profiling data from multiple PCa cohorts, we previously observed concurrent loss of the ETS2 repressor factor (ERF) in CIC-deficient prostate tumors (Huang et al., 2017). Combinatorial loss is most commonly the result of focal deletions (homozygous and heterozygous) at the 19q13.2 locus, where CIC and ERF are physically adjacent (long and short isoforms of CIC are separated from ERF by approximately 15 and 30 kb, respectively) to one another in the genome. Since ERF is a transcriptional repressor that binds ETS DNA motifs (Bose et al., 2017), we hypothesized that in a fusion independent manner, CIC and ERF cooperate to mutually suppress ETS target genes in PCa.
Through an integrative genomic and functional analysis, we mechanistically show that CIC and ERF directly bind and co-repress a proximal ETV1 regulatory element limiting PCa progression. Concomitant loss of CIC and ERF de-represses ETV1-mediated transcriptional programs and confer ETV1 dependence in multiple PCa model systems. Thus, we reveal a fusion-independent mechanism to de-repress ETS-mediated PCa progression and potentially uncover a therapeutic approach to target CIC-ERF co-deleted PCa.
Results
CIC is a TF that directly suppresses ETV1, ETV4, and ETV5 TF family members (Futran et al., 2015; Jiménez et al., 2012; Kim et al., 2021; Okimoto et al., 2017). CIC silences target genes through direct binding of a highly conserved DNA-binding motif (T[G/C]AAT[G/A]AA; Figure 1A; Ajuria et al., 2011; Futran et al., 2015; Jiménez et al., 2012). CIC is commonly altered in multiple human cancer subtypes where it suppresses tumor growth and metastasis (Kim et al., 2021). CIC is located on chromosome 19q13.2, directly adjacent to another transcriptional repressor, namely the ERF (Figure 1B). ERF binds and competes for ETS TF-binding sites (GGAA-motifs) and is frequently altered in human PCa, predominantly through focal deletions (Bose et al., 2017; Hou et al., 2020). We thus hypothesized that concurrent loss of CIC and ERF may de-repress an ETS-driven transcriptional program that drives PCa progression in a fusion-independent manner. To explore this, we first queried 15 PCa datasets curated on cBioPortal (Cerami et al., 2012; Gao et al., 2013) and identified a high cooccurrence rate (p<0.001, two-sided Fisher’s exact test [FET]) for CIC (10%) and ERF (12%) homozygous and heterozygous deletions (Figure 1C), suggesting that concurrent loss occurs through focal copy number change at the 19q13.2 locus. Through analysis of these clinically annotated specimens, we observed that the CIC-ERF co-deletion was present at an increased frequency in PCa with higher Gleason scores and later tumor stages when compared to CIC-ERF replete tumors (Figure 1D). In order to understand the association between CIC and/or ERF alterations in specific PCa cohorts, we stratified published datasets to identify patients that represent primary PCa (Fraser et al., 2017; Hieronymus et al., 2014; ‘ Cancer Genome Atlas Research Network, 2015) (PNAS 2014, n=272; Cell 2015, n=333; Nature 2017, n=477 primary PCas) and metastatic castrate resistant prostate cancer (mCRPC) (Abida et al., 2019; Grasso et al., 2012; Robinson et al., 2015) (Nature 2012, n=50; Cell 2015, n=150; PNAS 2019, n=429 mCRPCs). This analysis revealed enrichment of CIC and ERF alterations including the CIC-ERF co-deletion in mCRPC samples (Figure 1E, Figure 1—source data 1, Figure 1—figure supplement 1). Importantly, CIC-ERF co-deleted tumors clustered as a subgroup when compared to the more well-characterized molecular subsets including ERG, ETV1, ETV4, SPOP, and FOXA1 altered PCas, suggesting a distinct molecular subtype of PCa (Figure 1F). In order to explore clinical outcomes of patients harboring CIC-ERF co-deleted tumors, we performed a survival analysis using the aforementioned PCa datasets and observed significantly worse outcomes in patients who harbored the CIC-ERF co-deletion (p=0.001, disease-free survival [DFS] [ERF-CIC co-deletion 25 events/90 total; no ERF-CIC co-deletion 153 events/910 total] and p=0.01, progression-free survival [PFS] [ERF-CIC co-deletion 16 events/52 total; no ERF-CIC co-deletion 77 events/442 total]; Figure 1G). One limitation of this analysis is that individual studies included overlapping tumors from the same patient (i.e. TCGA PanCancer and Firehose Legacy cohorts). Despite this overlap, we utilized this insight to formulate a hypothesis that CIC and ERF are co-deleted with increasing frequency in mCRPC and that the CIC-ERF co-deletion is associated with worse clinical outcomes in PCa patients.
Our initial findings provided rationale to explore the genetic and functional relationship between CIC and ERF. Independently, CIC and ERF have previously been reported to promote malignant phenotypes, including tumor growth and metastasis in multiple human cancer subsets (Bettegowda et al., 2011; Bose et al., 2017; Choi et al., 2015; Huang et al., 2017; Kawamura-Saito et al., 2006; Kim et al., 2018; Bunda et al., 2019; Okimoto et al., 2019; Okimoto et al., 2017; Simón-Carrasco et al., 2017; Wong et al., 2019; Yang et al., 2017). Since our clinical data indicated that combinatorial loss of CIC and ERF was associated with worse patient outcomes, we hypothesized that CIC and ERF loss may cooperate to enhance PCa progression. To investigate this, we engineered ERF, CIC, or both CIC and ERF-deficient immortalized prostate epithelial cells (PNT2) and performed a series of in vitro and in vivo experiments to test the combinatorial effect of the CIC-ERF co-deletion. Compared to single gene loss of CIC or ERF, genetic inactivation of both CIC and ERF increased colony formation (Figure 2A–B, Figure 2—figure supplement 1A) and spheroid formation (Figure 2C–D) in PNT2 cells. CIC-ERF co-deletion also enhanced malignant phenotypes including cellular viability, invasiveness, and migratory capacity in PNT2 (Figure 2E–G). Additionally, genetic silencing of both CIC and ERF increased the frequency of subcutaneous tumor xenograft formation in severe-combined immunodeficient (SCID) mice compared to control (Figure 2H, Figure 2—figure supplement 1B). Thus, our findings demonstrate that the combination of CIC and ERF loss augments the transformation of PNT2 and promotes malignant phenotypes.
To assess the functional role of CIC and ERF in the context of human PCa progression, we leveraged two genetically annotated, androgen-insensitive PCa cell lines, DU-145 (moderate metastatic potential) and PC-3 (high-metastatic potential). DU-145 cells harbor a loss-of-function ERF mutation (ERFA132S) (Huang et al., 2017) and express functional WT CIC, (Figure 3—figure supplement 1A). By comparison, PC-3 cells are deficient in CIC (homozygous deletion) and retain functional ERF (Figure 3—figure supplement 1B; Cerami et al., 2012; Gao et al., 2013; Seim et al., 2017). Thus, these cell line models provided isogenic systems to functionally interrogate the role of CIC and ERF in human PCa. Specifically, genetic reconstitution of ERF into ERF-deficient DU-145 cells decreased colony formation in both CIC proficient (parental cells) and CIC knockout (KO) conditions (Figure 3A–B, Figure 3—figure supplement 1C, D). While CIC loss did not enhance colony formation in ERF-deficient DU-145 cells, it significantly increased tumor cell viability, invasion, and migratory capacity compared to control (Figure 3C–D, Figure 3—figure supplement 1E). Importantly, we observed that ERF expression in CIC KO DU-145 cells rescued the CIC-mediated effects on viability and migration/invasion (Figure 3C–D, Figure 3—figure supplement 1E). These findings indicate that rescuing ERF can partially restore the functional effects of CIC loss in DU-145 cells. To further understand if ERF could suppress tumor growth in vivo, we reconstituted WT ERF into CIC proficient and deficient (CIC KO) DU-145 cells and generated subcutaneous xenografts in immunodeficient mice (NU/J). Consistent with our in vitro data, CIC KO increased the tumor growth rate in vivo and genetic reconstitution of WT ERF partially suppressed tumor growth compared to DU-145 CIC KO cells (Figure 3E). Moreover, genetic reconstitution of ERF into CIC proficient DU-145 cells suppressed the tumor growth rate in vivo (Figure 3E). We next used PC-3 cells, to further test how ERF and CIC functionally interact in the context of human PCa. We first noted that ERF overexpression (OE) or reconstitution of CIC alone decreased PC-3 colony formation, with combinatorial ERF OE and CIC rescue having the most significant reduction compared to parental PC-3 cells (Figure 3F–G, Figure 3—figure supplement 1F, G). Moreover, ERF and CIC expression had a similar impact on decreasing PC-3 viability, invasion, and migratory capacity (Figure 3H–I, Figure 3—figure supplement 1H). Similarly, genetic suppression of ERF resulted in a significant increase of colony formation, viability, and invasion compared to parental PC3 cells (Figure 3F–I, Figure 3—figure supplement 1I). Interestingly, OE of ERF in mice bearing CIC deficient PC-3 tumor xenografts decreased tumor growth compared to PC-3 parental and PC-3 cells expressing WT CIC (genetic rescue of CIC; Figure 3J). Since we consistently observed that ERF expression could partially rescue the effects of CIC loss in PCa, we hypothesized that WT CIC and ERF potentially cooperate to limit PCa progression.
In order to mechanistically define how CIC and ERF (two TFs with known repressor function) were interacting to functionally regulate PCa, we performed chromatin immunoprecipitation followed by sequencing (ChIP-Seq) using a validated CIC antibody (Lin et al., 2020; Okimoto et al., 2019; Okimoto et al., 2017) in PNT2, and compared this to a publicly available ERF ChIP-Seq dataset in VCaP PCa cells (Bose et al., 2017), we were unsuccessful at pulling down ERF in PNT2 cells. This analysis identified 178 high-confidence (False-discovery rate (FDR) ≤0.05) CIC peaks that mapped to 130 annotated genes including known targets, ETV1, ETV4, and ETV5. Globally, CIC peaks were localized to distinct genomic regions including promoters (38.6%), untranslated regions (UTRs) (1.14%), introns (19.9%), and distal intergenic regions (38.64%). Interestingly, the distribution of ERF peaks was similar to CIC, with 32.8% in promoters, 1.2% in UTRs, 29.2% intronic, and 33.8% in distal intergenic regions (Figure 4A). Next, through a comparative ChIP-Seq analysis (Figure 4B), we identified 91 shared CIC and ERF target genes. Importantly, we focused on genes with shared CIC- and ERF-binding sites to potentially explain the functional cooperativity that we observed in our prior studies. In order to narrow down potential candidates, we performed Functional Clustering Analyses through the DAVID bioinformatics database (https://david.ncifcrf.gov) (Huang et al., 2009; Sherman et al., 2022) using the 91 shared CIC and ERF target genes (Supplementary file 1). Among these putative CIC and ERF targets, the PEA3 (ETV1, ETV4, and ETV5) TFs (known oncogenic drivers in PCa; Baena et al., 2013; Kim et al., 2016; Oh et al., 2012) were found to be the most highly enriched family. These findings suggested that CIC and ERF may co-regulate PEA3 family members through direct transcriptional control.
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/77072/elife-77072-fig4-data1-v1.zip
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/77072/elife-77072-fig4-data2-v1.zip
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/77072/elife-77072-fig4-data3-v1.zip
In order to further identify how CIC and/or ERF impact PEA3 TF expression, we performed quantitative real time polymerase chain reaction (qRT-PCR) in PNT2 cells, assessing ETV1, ETV4, or ETV5 mRNA levels in response to genetic silencing of CIC and/or ERF. As expected, CIC KO and/or combinatorial CIC and ERF loss in PNT2 cells (ERF and CIC WT) consistently increased ETV1, ETV4, and ETV5 levels compared to control (Figure 4—figure supplement 1A-F). In contrast, genetic silencing of ERF consistently increased ETV1 mRNA expression, but not ETV4 or ETV5 (Figure 4C, Figure 4—figure supplement 1G-J). These findings indicated that ETV1 may be a shared transcriptional target of CIC and ERF in prostate cells. In order to test a potential proteomic interaction between CIC and ERF, we performed co-immunoprecipitation (Co-IP) experiments but did not observe binding when either ERF or CIC was pulled down in HEK293T cells (Figure 4—figure supplement 1K-L). This led us to hypothesize that CIC and ERF may potentially bind (but not interact) to distinct regions along the ETV1 regulatory element. To explore this, we first localized CIC (TGAATGGA) and ERF (GGAA) DNA binding motifs within the proximal upstream regulatory element of ETV1 and independently confirmed CIC and ERF occupancy of the ETV1 promoter through ChIP-PCR (Figure 4D–G, Figure 4—figure supplement 1M, P). To extend these findings into the context of PCa, we reconstituted ERF in ERF deficient DU-145 cells and this consistently decreased ETV1 expression (not ETV4 or ETV5) in both CIC proficient and CIC deficient settings (Figure 4H, Figure 4—figure supplement 1Q-T). Moreover, ERF Knockdown (KD) or ERF OE in CIC deficient PC-3 cells increased and decreased ETV1 mRNA expression, respectively (Figure 4I–J). Since ETV1 is a known target of CIC (Dissanayake et al., 2011; Jiménez et al., 2012), we focused on further validating ETV1 as a molecular target of ERF. To this end, we engineered a ETV1 luciferase based promoter assays and observed a decrease in luciferase activity following ERF expression in 293T cells (Figure 4K). These genetic tools further validate that ERF can suppress ETV1 expression through direct transcriptional silencing of the ETV1 promoter and identifies ETV1 as an ERF target.
Consistent with a repressor function, ERF loss was previously shown to transcriptionally associate with ETV1-regulated gene set signatures (Huang et al., 2017). Yet it remains unclear if ERF can directly regulate ETV1 and how combinatorial loss of CIC and ERF controls ETV1-mediated (or other ETS family members) transcriptional programs. To explore this, we generated a signature gene set of upregulated genes from our dual CIC and ERF deficient PNT2 cells (ERF KD +CIC KO) and projected the Cancer Genome Atlas PCa (TCGA-PRAD) dataset onto the transcriptional space of these signature gene sets using the ssGSEA module (Version 10.0.9) on GenePattern (Reich et al., 2006). We found that CIC and ERF loss were significantly associated with the ETV1-regulated gene set (Information coefficient (IC)=0.619, p=0.0009), but was anti-correlated with the TMPRSS2-ERG fusion signature gene set (Setlur et al., 2008; Figure 4L, IC = –0.45, p=0.0009). Thus, the enrichment of the ETV1-regulated gene set signature was shared between ERF loss alone (Huang et al., 2017) and CIC-ERF dual suppression (Figure 4L). In contrast, combinatorial CIC and ERF loss negatively correlated with the TMPRSS2-ERG fusion signature gene set, which was not consistent with our prior studies using ERF KD alone (Huang et al., 2017). These findings led us to hypothesize that the dual suppression of CIC and ERF may increase ETV1-mediated transcriptional programs in PCa. This was supported by two major lines of experimental and conceptual evidence including: (1) dual suppression of ETV1 expression and ETV1 mediated transcriptional output by CIC and ERF; and (2) the majority of tumors derived from PCa patients that harbor ERF deletions also contain deletions in CIC.
In order to demonstrate enhanced survival dependence on ETV1 in PCa cells, we genetically silenced CIC in ERF-deficient DU-145 cells and assessed drug sensitivity to an ETV1 inhibitor (BRD32048) (Pop et al., 2014). We validated the pharmacologic effect of ETV1 inhibition (BRD32048) through suppression of known downstream target genes (ATAD2 and ID2) using qPCR in LNCaP cells (Figure 5—figure supplement 1A-C). Since BRD32048 was previously shown to decrease invasiveness in ETV1-fusion positive PCa cells (LNCaP), but not significantly impact tumor cell viability, we unexpectedly observed that silencing CIC in DU-145 cells mildly enhanced sensitivity to BRD32048 (Figure 5A). To further confirm these findings and to mitigate potential off-target effects, we silenced ETV1 using siRNA in PC3 and DU-145 cells expressing Crispr-based sgRNA targeting CIC (Figure 5—figure supplement 1D-F). Consistent with our pharmacologic studies, the viability of DU-145 CIC KO cells was decreased upon genetic ETV1 inhibition (Figure 5B). Similarly, pharmacologic and genetic ETV1 inhibition decreased invasiveness of CIC deficient DU-145 cells (Figure 5C–D). These findings indicate that loss of CIC in ERF deficient PCa cells can potentially modulate the sensitivity to ETV1-directed therapies. To expand these findings, we also consistently observed a decrease in viability and invasiveness upon ETV1 inhibition (BRD32048) in PNT2 cells with dual ERF and CIC suppression and in PC-3 cells with ERF KD (Figure 5—figure supplement 2A-J). These in vitro findings indicate that ETV1 inhibition in CIC and ERF deficient prostate cells can suppress invasion and potentially limit viability. Further studies targeting ETV1 in patient derived specimens that harbor endogenous CIC-ERF co-deletions is warranted. Collectively, our data indicate that CIC and ERF may cooperate to silence ETV1 transcriptional programs, limiting ETV1-mediated PCa progression.
Discussion
Molecular and functional subclassification of human PCa has revealed a dependence on ETS family TFs including ERG, ETV1, ETV4, and ETV5 (Feng et al., 2014; Oh et al., 2012). The predominant mode of ETS activation in PCa is through chromosomal rearrangements that fuse ERG, ETV1, ETV4, and ETV5 to the androgen-regulated TMPRSS2 gene, leading to fusion oncoproteins that drive oncogenesis (Clark and Cooper, 2009; Feng et al., 2014; Helgeson et al., 2008; Tomlins et al., 2006; Tomlins et al., 2005). Interestingly, recent data indicate that ETV1, ETV4, and ETV5 are upregulated in a fusion-independent manner and are associated with poor clinical outcomes in PCa patients (Baena et al., 2013; Hermans et al., 2008). Our study focused on understanding the molecular mechanisms that drive fusion-independent upregulation of ETS family members and we reveal a new molecular subclass of PCa defined by a co-deletion of two TFs, CIC and ERF.
CIC is a TF that directly silences ETV1, ETV4, and ETV5 transcription through direct repression at proximal regulatory sites (Jiménez et al., 2012). We observed that in ~10–12% of human PCa, CIC and ERF are co-deleted through focal homozygous or heterozygous deletions. It has been recently shown that ERF competes for ETS DNA binding motifs and our studies identify cooperative regulation of key target genes between CIC and ERF. Specifically, through ChIP-Seq analysis coupled with a series of in vitro and in vivo studies, we identify a coordinated binding of CIC and ERF to the proximal ETV1 regulatory element that physically and functionally regulates ETV1 expression. Therefore, we reveal ETV1 as a novel ERF target gene. Rescuing ERF in CIC deficient PCa cells decreases ETV1 expression and limits malignant phenotypes including viability, migration and invasion.
The 19q13.2 locus contains CIC and ERF, which are physically adjacent and oriented in opposing directions. The long and short isoforms of CIC are separated from ERF by ~15 and 30 kb, respectively. Thus, future studies directed at defining the genome topology and epigenetic states, both within and around this highly conserved region are warranted. In particular, studies aimed at mapping topology associated domains and key histone marks can potentially reveal shared upstream regulatory elements including enhancers or superenhancers that may co-regulate CIC and ERF in concert. These findings could reveal non-genetic mechanisms to functionally regulate CIC and ERF expression in a coordinated fashion.
Beyond PCa, we and others have observed that CIC and ERF are co-deleted in a subset of stomach adenocarcinoma (LeBlanc et al., 2017; Okimoto et al., 2017). Thus, combined CIC and ERF loss are not entirely specific for PCa. We speculate that since prostate (Baena et al., 2013) and potentially stomach adenocarcinoma (Keld et al., 2011) are highly dependent on PEA3 transcriptional dysregulation to enhance tumor progression, the dual loss of CIC and ERF may, in part, represent an alternative mode to de-repress ETS mediated transcriptional programs in these cancer subsets.
The use of ETV1 inhibitors has been limited to preclinical studies (Pop et al., 2014). These studies have largely focused on direct targeting of ETV1 fusion oncoproteins in PCa (Pop et al., 2014). Our findings indicate that patients with CIC-ERF deficient PCa may have a survival dependence on ETV1 and that these patients may potentially benefit from ETV1 directed therapies. A limitation of this targeted approach in PCa patients is that we did not find a statistically significant difference in ETV1 transcript levels in CIC-ERF co-deleted tumors compared to CIC-ERF replete tumors in the TCGA-PRAD (n=455) dataset. One potential explanation is that ETV1 is commonly upregulated in human PCas through mechanisms beyond CIC-ERF loss. Thus, while our data support the upregulation and potential induced dependence on ETV1 in our CIC-ERF deficient systems, it remains unclear if this will translate beyond our cell-line based models into PCa patients that harbor CIC-ERF co-deletions. Thus, larger studies that aim to evaluate: (1) ETV1 mRNA and protein expression levels; (2) the biological significance of ETV1 function; and (3) the clinical application of ETV1 inhibitors in patients with endogenous CIC-ERF co-deleted tumors (compared to CIC-ERF WT tumors) is warranted in this subset of PCa. Collectively, we have uncovered a molecular subset of PCa defined by a co-deletion of CIC and ERF and further demonstrate a mechanism-based strategy to potentially limit tumor progression through ETV1 inhibition in this subset of human PCa.
Materials and methods
Cell lines, drug, and reagents
Request a detailed protocolCell lines were cultured as recommended by the American Type Culture Collection (ATCC) or Sigma. DU-145, PC-3, and HEK293T cells were purchased and authenticated (STR profiling) by ATCC. PNT2 cells were purchased and authenticated (STR profiling) by Sigma. All cell lines were tested for mycoplasma using the e-myco PLUS PCR detection kit (Boca Scientific, Cat#25237). HEK293T cells are human embryonic kidney cells which are commonly used for transfection. PNT2 ERF KD (shERFA1, shERFB1), PNT2 CIC KO (sgCIC1, sgCIC2) and PNT2 ERF KD +CIC KO were derived from parental PNT2 cells. shRNAs targeting ERF to develop PNT2 shERFA1 and PNT2 shERFB1 were obtained from Sigma-Aldrich: TRCN000001391, TRCN0000013912. Puromycin (1 μg/ml) was used as a selection reagent. Two sgRNAs targeting CIC were previously validated and were gifts from William Hahn, Addgene (#74959 and #74953). These sgRNAs were used to develop PNT2 sgCIC1 and PNT2 sgCIC2 cells. Blasticidin (10 μg/ml) was used as a selection agent. PNT2shERFA1 +sgCIC1 and PNT2shERFA1 +sgCIC2 were developed from the combination of the above two shRNA and sgRNAs. All PNT2 cells were grown in DMEM media supplemented with 10% FBS, 100 IU/ml penicillin and 100 μ g/ml streptomycin.
DU-145 ERF, DU-145 CIC KO (sgCIC1, sgCIC2) and DU-145 ERF +CIC KO (ERF +sgCIC1, ERF +sgCIC2) were derived from parental DU-145 cell line. Lentiviral GFP-tagged ERF (GeneCopoeia, EX-S0501-Lv122) was used to develop DU-145 ERF cells with puromycin (1 μg/ml) as a selection marker. The above-mentioned two sgRNAs were used to develop DU-145 sgCIC1 and DU-145 sgCIC2 cells with blasticidin (15 μg/ml) as the selection agent. DU-145 ERF +sgCIC1 and DU-145 ERF +sgCIC2 were developed from the combination of the above two.
PC-3 ERF KD (shERFA1), PC-3 ERF OE, PC-3 CIC OE, PC-3 (ERF +CIC) OE cells were derived from parental PC-3 cell line. shRNAs targeting ERF (Sigma-Aldrich: TRCN000001391) and lentiviral GFP-tagged ERF (GeneCopoeia, EX-S0501-Lv122) were used to develop PC-3 shERFA1 and PC-3 ERF OE, respectively. PC-3 CIC OE cells were developed using CIC-Myc-tag plasmid purchased from Origene (CAT#: RC215209).
Geneticin (250 μg/ml) was used as a selection agent. PC-3 (ERF +CIC) OE cells were developed using a combination of ERF-GFP and CIC-Myc overexpressing plasmid.
All DU-145 and PC-3 cells were grown in RPMI 1640 media supplemented with 10% FBS, 100 I U/ml penicillin and 100 μ g/ml streptomycin, respectively. All cell lines were maintained at 37 °C in a humidified atmosphere at 5% CO2. All the above mentioned stable cell lines were validated by analyzing the expression of CIC and ERF using qPCR and western blot analysis (Figure 2—figure supplement 1A, Figure 3—figure supplement 1A-C, I).
BRD32048 is an ETV1 inhibitor that was purchased from Sigma-Aldrich (CAT#:SML1346).
Analysis of PCa datasets from cBioPortal
Request a detailed protocol15 PCa datasets (see Figure 1c for individual studies) were queried for alterations in CIC and ERF using the cBioPortal platform ‘query by gene’ function. Stratification into ‘ERF-CIC No co-deletion’ and ‘ERF-CIC co-deletion’ was performed in cBioPortal and associated with ‘Gleason Score’, ‘AJCC Primary Tumor T Stage’, and DFS and PFS using the ‘Plots’ and ‘Comparison/Survival’ functions. p-values for comparison between Gleason and Tumor Stage were calculated using FET and survival curves were calculated by Log-rank test.
For the CIC and ERF mutational analysis in primary prostate vs. metastatic castrate resistant PCa, we selected studies that purely represented primary prostate tumors (PNAS 2014, Cell 2014, Nature 2017) and advanced mCRPC (Nature 2012, Cell 2015, PNAS 2019).
Colony formation assay
Request a detailed protocolEqual number of cells from cell lines (500–600 cells/well) were seeded in a 6-well plate. Cells were allowed to form colonies for 7 days. At day 7, cells were fixed and stained with 0.5% crystal violet solution after washing with PBS and performed in triplicate.
Finally, the colonies with >50 cells were counted under an imageJ software.
Tumorsphere assay
Request a detailed protocolApproximately 25,000 cells from different groups were cultured in tumorsphere media at 37 °C and 5% CO2 for 7 days. Tumorsphere medium contains serum free DMEM /F12 supplemented with 10 ng/ml FGF (fibroblast growth factor), 20 ng/ml EGF (epidermal growth factor), 1xITS (Insulin-Transferrin-Selenium) and B27 supplement. On day 7, images of different areas of the wells were taken using confocal. The size of the sphere was calculated using Fiji (ImageJ) software in all the tested groups. Each group consisted of three replicate wells and at least 6 images (n>/=6).
Subcutaneous tumor xenograft assays
Request a detailed protocolFour week old male SCID mice were purchased from Jackson Laboratory. Mice were kept under specific pathogen-free conditions and facilities were approved by the UCSF IACUC. To prepare cell suspensions, PNT2 and its other genetic variants (PNT2 ERF KD, PNT2 CIC KO, PNT2 ERF KD +CIC KO) were briefly trypsinized, quenched with 10% FBS RPMI media and resuspended in PBS. Cells were pelleted again and mixed with PBS/Matrigel matrix (1:1) for a final concentration of 0.1×105 cells/µl. A 100 μl cell suspension containing 1×106 cells were injected (s.c.) in the right and left flanks of immunodeficient mice (n=10/group). Mice were observed for tumor formation in different groups over 7 weeks. For other subcutaneous xenografts, four week old male mice (NU/J) were purchased from Jackson Laboratory and were six-eight weeks old at time of experiment. 1.0×106 DU-145 cells and its variants (DU-145 ERF, DU-145 CICKO and DU-145 ERF +CICKO) were resuspended in PBS/Matrigel (1:1) matrix and injected s.c. into the right and left flanks of nude mice (n=10/group). Tumor volume was measured twice per week using Vernier caliper. Tumor volume was determined using caliper measurements of tumor length (L) and width (W) according to the formula V = (L X W2) X 0.52.
Similarly, 1.0×106 PC-3 tumor cells and its variants (PC-3 ERFOE, CICOE) were injected subcutaneously in flanks of male nude mice (NU/J), n=10/group and tumor volume was monitored in different groups. Mice body weight was measured in all the experiments throughout the study. At the end of the experiment, mice were sacrificed by CO2 overdose in accordance with IACUC guidelines.
Viability assays
Request a detailed protocol5000 cells were plated in a 12 well plate. Crystal violet staining was performed 5 days after cell plating with 3 replicates per group. CellTiter Glo experiments were performed according to the manufacturer’s protocol. In brief, cells were plated in a 96-well plate, and analyzed on a Spectramax microplate reader (Molecular Devices) at different days. Each assay was performed with six replicate wells.
Transwell invasion assays
Request a detailed protocolRPMI with 10% FBS was added to the bottom well of a trans-well chamber. 2.5×10^4 cells resuspended in serum free media was then added to the top 8 µm pore matrigel coated (invasion) or non-coated (migration) trans-well insert (BD Biosciences). After 20 hr, non-invading cells on the apical side of inserts were scraped off and the trans-well membrane was fixed in methanol for 15 min and stained with Crystal Violet for 30 min. The basolateral surface of the membrane was visualized with a Zeiss Axioplan II immunofluorescent microscope at 5×. Each trans-well insert was imaged in five distinct regions at 5×and performed in triplicate. % invasion was calculated by dividing the mean # of cells invading through Matrigel membrane / mean # of cells migrating through control insert.
Wound healing assays
Request a detailed protocolCells were plated at a density of 0.5×106 cells/well and incubated to form a monolayer in 6-well dishes. Once a uniform monolayer was formed, wound was created by scratching the monolayer with a 1 ml sterile tip. Floating cells were removed by washing the cells with PBS three times. Images of the wound were taken at this point using bright field microscope and considered as a 0 hr time point. Furthermore, media was added in all the wells and cells were left to migrate either for 24 hr (DU-145 cells) or 48 hr (PNT2 ond PC-3 cells). At end point, wound was imaged again using bright field microscope. The wound area at different points was calculated using ImageJ software. Each group consisted of at least three replicate wells.
DAVID functional clustering analysis
Request a detailed protocolThe 91 shared candidate target genes between CIC and ERF identified through our ChIPSeq analysis in PNT2 and VCaP cells, respectively, were used as an input list for analysis using the DAVID Bioinformatics Database (https://david.ncifcrf.gov).
Real-time quantitative PCR (RT-qPCR)
Request a detailed protocolRT-qPCR was performed in PNT2, DU145 and PC3 cells. Isolation and purification of RNA was performed using RNeasy Mini Kit (Qiagen). 500 ng of total RNA was used in a reverse transcriptase reaction with the SuperScript III firststrand synthesis system (Invitrogen). Quantitative PCR included three replicates per cDNA sample. Human CIC (Cat#. Hs00943425_g1), ERF (Cat#. Hs01100070_g1), ETV1 (Cat#. Hs00951951_m1), ETV4 (Cat#. Hs00383361_g1), ETV5 (Cat#. Hs00927557_m1), and endogenous controls GAPDH (Cat#. Hs02758991_g1) were amplified with Taqman gene expression assay (Applied Biosystems). Expression data were acquired using an ABI Prism 7900HT Sequence Detection System (Thermo Fisher Scientific). Expression of each target was calculated using the 2−ΔΔCt method and expressed as relative mRNA expression.
Chip-Seq and PCR
Request a detailed protocolChIP was performed on PNT2 and DU-145 cells with the SimpleChIP Enzymatic Chromatin IP kit, Cell Signaling Technology #9003 in accordance with the manufacturer’s protocol. The antibodies used for IP were as follows: CIC (Thermo Fisher Scientific – PA146018) and ERF (Thermo Fisher Scientific –PA530237). Paired-end 150 bp (PE150) sequencing on an Illumina HiSeq platform was then performed. ChIP-Seq peak calls were identified through Mode-based Analysis of ChIP-Seq (MACS). For ETV1 ChIP-PCR validation, primers were designed in the proximal regulatory element of ETV1. The promoter primer sequences are listed in supplementary experimental methods:
VCaP ERF ChIP-Seq was previously performed (Bose et al., 2017) and publicly available in the GEO database: GSE83653. For this analysis, the samples used are listed in supplementary experimental methods.
Promoter primer sequences for ChIP-PCR and ChIP-Seq are provided in the Supplementary Methods (Supplementary file 2).
GSM2612455_INPUT red replicate two luciferase promoter assay
Request a detailed protocol293T and DU-145 cells were split into a 96 well plate to achieve 50–70% confluence the day of transfection. LightSwitch luciferase assay system (SwitchGear Genomics) was used per the manufacturer’s protocol. Briefly, a mixture containing FuGENE 6 transfection reagent, 50 ng Luciferase GoClone ETV1 promoter (S720645) plasmid DNA, 50 ng of either control (empty) vector or fully sequenced ERF cDNA (GeneCopoeia [EX-S0501Lv122]), were added to each well. All transfections were performed in quintuplicate. The plates were assessed for luciferase activity after 48 hr of treatment.
Single-sample gene set enrichment analysis (ssGSEA)
Request a detailed protocolGene-level expression of our dual CIC and ERF deficient PNT2 cells are computed using RSEM (Setlur et al., 2008) (Version 1.3.3) and log-2 normalized. The signature gene set of dual CIC and ERF deficient PNT2 cells (ERF KD +CIC KO) is defined as the top upregulated genes compared to the CIC and ERF WT cells. The ssGSEA module on GenePattern was then used to project the TCGA-PRAD dataset onto the transcriptional space defined by both the ERF KD +CIC KO signature gene set and previously established gene sets including the ETV1-regulated gene set, ERF KD signature and TMPRSS2-ERG fusion signature. The ssGSEA enrichment scores of the ERF KD +CIC KO signature gene set for the TCGA-PRAD samples are compared with the scores of the other signature gene sets and visualized in heatmaps and used for downstream association analyses.
Gene knockdown (KD) and OE assays
Request a detailed protocolON-TARGET plus Scramble, ETV1 (L-003801-00-0005) and CIC (L-015185-01-0005) siRNAs were obtained from GE Dharmacon and transfection was performed with Dharmafect transfection reagent per manufacturer recommendations. ETV1 inhibitorBRD32048 (SML1346) was purchased from Millipore Sigma. Lentiviral GFP-tagged ERF was obtained from GeneCopoeia (EX-S0501-Lv122). pCMV-CIC with myc-tag was purchased from Origene and validated previously.
Experimental plan for experiments using ETV1 inhibitor
Request a detailed protocolCells were plated at a density of 0.2×106 in a 6-well plate. Next day, cells were transfected with siScramble and siCIC. After 48 hr, 25 μM and 50 μM of BMS32048 (ETV1 inhibitor) was added to both the transfected groups and control groups were treated with DMSO. 24 hr post ETV1 inhibitor treatment, cells from all the groups were replated for viability assay in a 24-well plate (n=3) and transwell invasion assay (n=3). Viability assay was performed using crystal violet. The relative density of crystal violet and % invasion was calculated comparing siCIC with siSCRM conditions in each respective drug concentration.
Western blot analysis
Request a detailed protocolAdherent cells were washed and lysed with RIPA buffer supplemented with proteinase and phosphatase inhibitors. Proteins were separated by SDS-PAGE, transferred to nitrocellulose membranes and blotted with antibodies recognizing: CIC (Thermo Fisher Scientific –PA146018), ERF (Thermo Fisher Scientific –PA530237), ETV1 (Thermo Fisher Scientific –MA515461), HSP90 (Cell Signaling– 4874 S), Actin (Cell Signaling – 4970 S). All immunoblots represent at least two independent experiments.
Co-immunoprecipitation assays
Request a detailed protocol293T cells were transfected with GFP-tagged ERF for 48 hr, lysed, quantified, and incubated with either IgG (Cell Signaling Technology; 2729) fused to Dynabeads Protein G (Thermo Fisher Scientific; 10004D) or anti-GFP magnetic beads overnight at 4 °C. Proteins were separated by SDS-PAGE, transferred to nitrocellulose membranes, and blotted with antibodies recognizing GFP or CIC. Myc-tagged CIC was transfected into HEK293T cells for 48 hr, lysed, quantified, and incubated with either IgG (Cell Signaling Technology; 2729) fused to Dynabeads Protein G (Thermo Fisher Scientific; 10004D) or anti-myc magnetic beads overnight at 4 °C. Proteins were separated by SDS-PAGE, transferred to nitrocellulose membranes, and blotted with antibodies recognizing Myc or ERF.
Statistical analysis
Request a detailed protocolExperimental data are presented as mean +/-Standard Deviation (SD) or Standard Error of the Mean (SEM). P-values derived for all in-vitro experiments were calculated using two-tailed student’s t-test or one-way ANOVA. The detailed statistical analysis performed for each experiment is defined in the figure legends.
Study approval
Request a detailed protocolFor tumor xenograft studies, specific pathogen-free conditions and facilities were approved by the American Association for Accreditation of Laboratory Animal Care. Surgical procedures were reviewed and approved by the UCSF Institutional Animal Care and Use Committee (IACUC), protocol #AN178670-03.
Data availability
All sequencing data including RNASeq and ChIP seq have been deposited in GEO under accession codes GSE216732 and GSE216318, respectively.
-
NCBI Gene Expression OmnibusID GSE216732. Decoding the Protein Composition of Whole Nucleosomes with Nuc-MS.
-
NCBI Gene Expression OmnibusID GSE216318. The CIC-ERF co-deletion underlies fusion independent activation of ETS family member, ETV1, to drive prostate cancer progression.
-
NCBI Gene Expression OmnibusID GSE83653. Loss of Function Mutations in ETS2 Repressor Factor (ERF) Reveal a Balance Between Positive and Negative ETS Factors Controlling Prostate Oncogenesis.
References
-
CIC protein instability contributes to tumorigenesis in glioblastomaNature Communications 10:661.https://doi.org/10.1038/s41467-018-08087-9
-
ETS gene fusions in prostate cancerNature Reviews. Urology 6:429–439.https://doi.org/10.1038/nrurol.2009.127
-
Molecular pathways: targeting ETS gene fusions in cancerClinical Cancer Research 20:4442–4448.https://doi.org/10.1158/1078-0432.CCR-13-0275
-
The capicua repressor--a general sensor of RTK signaling in development and diseaseJournal of Cell Science 125:1383–1391.https://doi.org/10.1242/jcs.092965
-
Histone demethylase JMJD2A drives prostate tumorigenesis through transcription factor ETV1The Journal of Clinical Investigation 126:706–720.https://doi.org/10.1172/JCI78132
-
Prostate cancerThe New England Journal of Medicine 349:366–381.https://doi.org/10.1056/NEJMra021562
-
ETV1, 4 and 5: an oncogenic subfamily of ETS transcription factorsBiochimica et Biophysica Acta 1826:1–12.https://doi.org/10.1016/j.bbcan.2012.02.002
-
Inactivation of capicua drives cancer metastasisNature Genetics 49:87–96.https://doi.org/10.1038/ng.3728
-
CIC-DUX4 oncoprotein drives sarcoma metastasis and tumorigenesis via distinct regulatory programsThe Journal of Clinical Investigation 129:3401–3406.https://doi.org/10.1172/JCI126366
-
A small molecule that binds and inhibits the ETV1 transcription factor oncoproteinMolecular Cancer Therapeutics 13:1492–1502.https://doi.org/10.1158/1535-7163.MCT-13-0689
-
Whole-genome sequence of the metastatic PC3 and lncap human prostate cancer cell linesG3: Genes, Genomes, Genetics 7:1731–1741.https://doi.org/10.1534/g3.117.039909
-
Estrogen-dependent signaling in a molecularly distinct subclass of aggressive prostate cancerJournal of the National Cancer Institute 100:815–825.https://doi.org/10.1093/jnci/djn150
-
Inactivation of capicua in adult mice causes T-cell lymphoblastic lymphomaGenes & Development 31:1456–1468.https://doi.org/10.1101/gad.300244.117
-
The ETS family of oncogenic transcription factors in solid tumoursNature Reviews. Cancer 17:337–351.https://doi.org/10.1038/nrc.2017.20
-
Impact of capicua on pancreatic cancer progressionAnnals of Surgical Oncology 28:3198–3207.https://doi.org/10.1245/s10434-020-09339-z
Article and author information
Author details
Funding
National Cancer Institute (K08CA222625)
- Ross A Okimoto
National Cancer Institute (5R37CA255453)
- Ross A Okimoto
National Cancer Institute (P30CA082103)
- Ross A Okimoto
Benioff Initiative for Prostate Cancer Research Awards
- Felix Y Feng
National Research Foundation of Korea (NRF- 2020R1A6A3A03039483)
- Ji Won Kim
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank members of the Okimoto and Huang labs and funding support from the following sources: A National Cancer Institute (NCI) Cancer Center Support Grant (P30CA082103) and Benioff Initiative for Prostate Cancer Research Awards (N.G and R.A.O), the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF- 2020R1A6A3A03039483)(J.W.K) and NCI K08-CA222625 and R37CA255453 awards to R.A.O.
Ethics
For tumor xenograft studies, specific pathogen-free conditions and facilities were approved by the American Association for Accreditation of Laboratory Animal Care. Surgical procedures were reviewed and approved by the UCSF Institutional Animal Care and Use Committee (IACUC), protocol #AN178670-03.
Copyright
© 2022, Gupta et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 755
- views
-
- 94
- downloads
-
- 4
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
Metastasis is the leading cause of cancer-related mortality. Paneth cells provide stem cell niche factors in homeostatic conditions, but the underlying mechanisms of cancer stem cell niche development are unclear. Here, we report that Dickkopf-2 (DKK2) is essential for the generation of cancer cells with Paneth cell properties during colon cancer metastasis. Splenic injection of Dkk2 knockout (KO) cancer organoids into C57BL/6 mice resulted in a significant reduction of liver metastases. Transcriptome analysis showed reduction of Paneth cell markers such as lysozymes in KO organoids. Single-cell RNA sequencing analyses of murine metastasized colon cancer cells and patient samples identified the presence of lysozyme positive cells with Paneth cell properties including enhanced glycolysis. Further analyses of transcriptome and chromatin accessibility suggested hepatocyte nuclear factor 4 alpha (HNF4A) as a downstream target of DKK2. Chromatin immunoprecipitation followed by sequencing analysis revealed that HNF4A binds to the promoter region of Sox9, a well-known transcription factor for Paneth cell differentiation. In the liver metastatic foci, DKK2 knockout rescued HNF4A protein levels followed by reduction of lysozyme positive cancer cells. Taken together, DKK2-mediated reduction of HNF4A protein promotes the generation of lysozyme positive cancer cells with Paneth cell properties in the metastasized colon cancers.
-
- Cancer Biology
- Computational and Systems Biology
This study investigates the variability among patients with non-small cell lung cancer (NSCLC) in their responses to immune checkpoint inhibitors (ICIs). Recognizing that patients with advanced-stage NSCLC rarely qualify for surgical interventions, it becomes crucial to identify biomarkers that influence responses to ICI therapy. We conducted an analysis of single-cell transcriptomes from 33 lung cancer biopsy samples, with a particular focus on 14 core samples taken before the initiation of palliative ICI treatment. Our objective was to link tumor and immune cell profiles with patient responses to ICI. We discovered that ICI non-responders exhibited a higher presence of CD4+ regulatory T cells, resident memory T cells, and TH17 cells. This contrasts with the diverse activated CD8+ T cells found in responders. Furthermore, tumor cells in non-responders frequently showed heightened transcriptional activity in the NF-kB and STAT3 pathways, suggesting a potential inherent resistance to ICI therapy. Through the integration of immune cell profiles and tumor molecular signatures, we achieved an discriminative power (area under the curve [AUC]) exceeding 95% in identifying patient responses to ICI treatment. These results underscore the crucial importance of the interplay between tumor and immune microenvironment, including within metastatic sites, in affecting the effectiveness of ICIs in NSCLC.