Association of tamoxifen resistance and lipid reprogramming in breast cancer.

BACKGROUND: Tamoxifen treatment of estrogen receptor (ER)-positive breast cancer reduces mortality by 31%. However, over half of advanced ER-positive breast cancers are intrinsically resistant to tamoxifen and about 40% will acquire the resistance during the treatment. METHODS: In order to explore mechanisms underlying endocrine therapy resistance in breast cancer and to identify new therapeutic opportunities, we created tamoxifen-resistant breast cancer cell lines that represent the luminal A or the luminal B. Gene expression patterns revealed by RNA-sequencing in seven tamoxifen-resistant variants were compared with their isogenic parental cells. We further examined those transcriptomic alterations in a publicly available patient cohort. RESULTS: We show that tamoxifen resistance cannot simply be explained by altered expression of individual genes, common mechanism across all resistant variants, or the appearance of new fusion genes. Instead, the resistant cell lines shared altered gene expression patterns associated with cell cycle, protein modification and metabolism, especially with the cholesterol pathway. In the tamoxifen-resistant T-47D cell variants we observed a striking increase of neutral lipids in lipid droplets as well as an accumulation of free cholesterol in the lysosomes. Tamoxifen-resistant cells were also less prone to lysosomal membrane permeabilization (LMP) and not vulnerable to compounds targeting the lipid metabolism. However, the cells were sensitive to disulfiram, LCS-1, and dasatinib. CONCLUSION: Altogether, our findings highlight a major role of LMP prevention in tamoxifen resistance, and suggest novel drug vulnerabilities associated with this phenotype.

Systematic drug screening reveals specific vulnerabilities and co-resistance patterns in endocrine-resistant breast cancer.

BACKGROUND: The estrogen receptor (ER) inhibitor tamoxifen reduces breast cancer mortality by 31 % and has served as the standard treatment for ER-positive breast cancers for decades. However, 50 % of advanced ER-positive cancers display de novo resistance to tamoxifen, and acquired resistance evolves in 40 % of patients who initially respond. Mechanisms underlying resistance development remain poorly understood and new therapeutic opportunities are urgently needed. Here, we report the generation and characterization of seven tamoxifen-resistant breast cancer cell lines from four parental strains. METHODS: Using high throughput drug sensitivity and resistance testing (DSRT) with 279 approved and investigational oncology drugs, exome-sequencing and network analysis, we for the first time, systematically determine the drug response profiles specific to tamoxifen resistance. RESULTS: We discovered emerging vulnerabilities towards specific drugs, such as ERK1/2-, proteasome- and BCL-family inhibitors as the cells became tamoxifen-resistant. Co-resistance to other drugs such as the survivin inhibitor YM155 and the chemotherapeutic agent paclitaxel also occurred. CONCLUSION: This study indicates that multiple molecular mechanisms dictate endocrine resistance, resulting in unexpected vulnerabilities to initially ineffective drugs, as well as in emerging co-resistances. Thus, combatting drug-resistant tumors will require patient-tailored strategies in order to identify new drug vulnerabilities, and to understand the associated co-resistance patterns.

Transient cyclical methylation of promoter DNA.

Methylation of CpG dinucleotides is generally associated with epigenetic silencing of transcription and is maintained through cellular division. Multiple CpG sequences are rare in mammalian genomes, but frequently occur at the transcriptional start site of active genes, with most clusters of CpGs being hypomethylated. We reported previously that the proximal region of the trefoil factor 1 (TFF1, also known as pS2) and oestrogen receptor alpha (ERalpha) promoters could be partially methylated by treatment with deacetylase inhibitors, suggesting the possibility of dynamic changes in DNA methylation. Here we show that cyclical methylation and demethylation of CpG dinucleotides, with a periodicity of around 100 min, is characteristic for five selected promoters, including the oestrogen (E2)-responsive pS2 gene, in human cells. When the pS2 gene is actively transcribed, DNA methylation occurs after the cyclical occupancy of ERalpha and RNA polymerase II (polII). Moreover, we report conditions that provoke methylation cycling of the pS2 promoter in cell lines in which pS2 expression is quiescent and the proximal promoter is methylated. This coincides with a low-level re-expression of ERalpha and of pS2 transcripts.

Inhibition of histone methyltransferase DOT1L silences ERα gene and blocks proliferation of antiestrogen-resistant breast cancer cells.

Breast cancer (BC) resistance to endocrine therapy results from constitutively active or aberrant estrogen receptor α (ERα) signaling, and ways to block ERα pathway in these tumors are sought after. We identified the H3K79 methyltransferase DOT1L as a novel cofactor of ERα in BC cell chromatin, where the two proteins colocalize to regulate estrogen target gene transcription. DOT1L blockade reduces proliferation of hormone-responsive BC cells in vivo and in vitro, consequent to cell cycle arrest and apoptotic cell death, with widespread effects on ER-dependent gene transcription, including ERα and FOXA1 gene silencing. Antiestrogen-resistant BC cells respond to DOT1L inhibition also in mouse xenografts, with reduction in ERα levels, H3K79 methylation, and tumor growth. These results indicate that DOT1L is an exploitable epigenetic target for treatment of endocrine therapy-resistant ERα-positive BCs.

Epigenetic reprogramming of HOXC10 in endocrine-resistant breast cancer.

Resistance to aromatase inhibitors (AIs) is a major clinical problem in the treatment of estrogen receptor (ER)-positive breast cancer. In two breast cancer cell line models of AI resistance, we identified widespread DNA hyper- and hypomethylation, with enrichment for promoter hypermethylation of developmental genes. For the homeobox gene HOXC10, methylation occurred in a CpG shore, which overlapped with a functional ER binding site, causing repression of HOXC10 expression. Although short-term blockade of ER signaling caused relief of HOXC10 repression in both cell lines and breast tumors, it also resulted in concurrent recruitment of EZH2 and increased H3K27me3, ultimately transitioning to increased DNA methylation and silencing of HOXC10. Reduced HOXC10 in vitro and in xenografts resulted in decreased apoptosis and caused antiestrogen resistance. Supporting this, we used paired primary and metastatic breast cancer specimens to show that HOXC10 was reduced in tumors that recurred during AI treatment. We propose a model in which estrogen represses apoptotic and growth-inhibitory genes such as HOXC10, contributing to tumor survival, whereas AIs induce these genes to cause apoptosis and therapeutic benefit, but long-term AI treatment results in permanent repression of these genes via methylation and confers resistance. Therapies aimed at inhibiting AI-induced histone and DNA methylation may be beneficial in blocking or delaying AI resistance.

Reanalysis of RNA-sequencing data reveals several additional fusion genes with multiple isoforms.

RNA-sequencing and tailored bioinformatic methodologies have paved the way for identification of expressed fusion genes from the chaotic genomes of solid tumors. We have recently successfully exploited RNA-sequencing for the discovery of 24 novel fusion genes in breast cancer. Here, we demonstrate the importance of continuous optimization of the bioinformatic methodology for this purpose, and report the discovery and experimental validation of 13 additional fusion genes from the same samples. Integration of copy number profiling with the RNA-sequencing results revealed that the majority of the gene fusions were promoter-donating events that occurred at copy number transition points or involved high-level DNA-amplifications. Sequencing of genomic fusion break points confirmed that DNA-level rearrangements underlie selected fusion transcripts. Furthermore, a significant portion (>60%) of the fusion genes were alternatively spliced. This illustrates the importance of reanalyzing sequencing data as gene definitions change and bioinformatic methods improve, and highlights the previously unforeseen isoform diversity among fusion transcripts.

KRAS oncogene rearrangements and gene fusions: unexpected rare encounters in late-stage prostate cancers.

Wang and colleagues identify a fusion between UBE2L3 and KRAS in a subset of metastatic prostate cancers.

Identification of fusion genes in breast cancer by paired-end RNA-sequencing.

BACKGROUND: Until recently, chromosomal translocations and fusion genes have been an underappreciated class of mutations in solid tumors. Next-generation sequencing technologies provide an opportunity for systematic characterization of cancer cell transcriptomes, including the discovery of expressed fusion genes resulting from underlying genomic rearrangements. RESULTS: We applied paired-end RNA-seq to identify 24 novel and 3 previously known fusion genes in breast cancer cells. Supported by an improved bioinformatic approach, we had a 95% success rate of validating gene fusions initially detected by RNA-seq. Fusion partner genes were found to contribute promoters (5' UTR), coding sequences and 3' UTRs. Most fusion genes were associated with copy number transitions and were particularly common in high-level DNA amplifications. This suggests that fusion events may contribute to the selective advantage provided by DNA amplifications and deletions. Some of the fusion partner genes, such as GSDMB in the TATDN1-GSDMB fusion and IKZF3 in the VAPB-IKZF3 fusion, were only detected as a fusion transcript, indicating activation of a dormant gene by the fusion event. A number of fusion gene partners have either been previously observed in oncogenic gene fusions, mostly in leukemias, or otherwise reported to be oncogenic. RNA interference-mediated knock-down of the VAPB-IKZF3 fusion gene indicated that it may be necessary for cancer cell growth and survival. CONCLUSIONS: In summary, using RNA-sequencing and improved bioinformatic stratification, we have discovered a number of novel fusion genes in breast cancer, and identified VAPB-IKZF3 as a potential fusion gene with importance for the growth and survival of breast cancer cells.

Systematic analysis of microRNAs targeting the androgen receptor in prostate cancer cells.

Androgen receptor (AR) is expressed in all stages of prostate cancer progression, including in castration-resistant tumors. Eliminating AR function continues to represent a focus of therapeutic investigation, but AR regulatory mechanisms remain poorly understood. To systematically characterize mechanisms involving microRNAs (miRNAs), we conducted a gain-of function screen of 1129 miRNA molecules in a panel of human prostate cancer cell lines and quantified changes in AR protein content using protein lysate microarrays. In this way, we defined 71 unique miRNAs that influenced the level of AR in human prostate cancer cells. RNA sequencing data revealed that the 3'UTR of AR (and other genes) is much longer than currently used in miRNA target prediction programs. Our own analyses predicted that most of the miRNA regulation of AR would target an extended 6 kb 3'UTR. 3'UTR-binding assays validated 13 miRNAs that are able to regulate this long AR 3'UTR (miR-135b, miR-185, miR-297, miR-299-3p, miR-34a, miR-34c, miR-371-3p, miR-421, miR-449a, miR-449b, miR-634, miR-654-5p, and miR-9). Fifteen AR downregulating miRNAs decreased androgen-induced proliferation of prostate cancer cells. In particular, analysis of clinical prostate cancers confirmed a negative correlation of miR-34a and miR-34c expression with AR levels. Our findings establish that miRNAs interacting with the long 3'UTR of the AR gene are important regulators of AR protein levels, with implications for developing new therapeutic strategies to inhibit AR function and androgen-dependent cell growth.

E2-mediated cathepsin D (CTSD) activation involves looping of distal enhancer elements.

Estrogen receptor alpha (ERalpha) is a ligand dependent transcription factor that regulates the expression of target genes through interacting with cis-acting estrogen response elements (EREs). However, only a minority of ERalpha binding sites are located within the proximal promoter regions of responsive genes. Here we report the characterization of an ERE located 9kbp upstream of the TSS of the cathepsin D gene (CTSD) that up-regulates CTSD expression upon estrogen stimulation in MCF-7 cells. Using ChIP, we show recruitment of ERalpha and phosphorylated PolII at the CTSD distal enhancer region. Moreover, we determine the kinetics of transient CpG methylation on the promoter region of CTSD and for the first time, at a distal enhancer element. We show that ERalpha is crucial for long-distance regulation of CTSD expression involving a looping mechanism.