Hallmarks of aromatase inhibitor drug resistance revealed by epigenetic profiling in breast cancer

Nucleosome reorganisation in breast cancer tissues

BACKGROUND

Nucleosome repositioning in cancer is believed to cause many changes in genome organisation and gene expression. Understanding these changes is important to elucidate fundamental aspects of cancer. It is also important for medical diagnostics based on cell-free DNA (cfDNA), which originates from genomic DNA regions protected from digestion by nucleosomes.

RESULTS

We have generated high-resolution nucleosome maps in paired tumour and normal tissues from the same breast cancer patients using MNase-assisted histone H3 ChIP-seq and compared them with the corresponding cfDNA from blood plasma. This analysis has detected single-nucleosome repositioning at key regulatory regions in a patient-specific manner and common cancer-specific patterns across patients. The nucleosomes gained in tumour versus normal tissue were particularly informative of cancer pathways, with ~ 20-fold enrichment at CpG islands, a large fraction of which marked promoters of genes encoding DNA-binding proteins. The tumour tissues were characterised by a 5-10 bp decrease in the average distance between nucleosomes (nucleosome repeat length, NRL), which is qualitatively similar to the differences between pluripotent and differentiated cells. This effect was correlated with gene activity, differential DNA methylation and changes in local occupancy of linker histone variants H1.4 and H1X.

CONCLUSIONS

Our study offers a novel resource of high-resolution nucleosome maps in breast cancer patients and reports for the first time the effect of systematic decrease of NRL in paired tumour versus normal breast tissues from the same patient. Our findings provide a new mechanistic understanding of nucleosome repositioning in tumour tissues that can be valuable for patient diagnostics, stratification and monitoring.

Breast cancer risk SNPs converge on estrogen receptor binding sites commonly shared between breast tumors to locally alter estrogen signalling output

Estrogen Receptor alpha (ERα) is the main driver and prime drug target in luminal breast. ERα chromatin binding is extensively studied in cell lines and a limited number of human tumors, using consensi of peaks shared among samples. However, little is known about inter-tumor heterogeneity of ERα chromatin action, along with its biological implications. Here, we use a large set of ERα ChIP-seq data from 70 ERα+ breast cancers to explore inter-patient heterogeneity in ERα DNA binding, to reveal a striking inter-tumor heterogeneity of ERα action. Interestingly, commonly-shared ERα sites showed the highest estrogen-driven enhancer activity and were most-engaged in long-range chromatin interactions. In addition, the most-commonly shared ERα-occupied enhancers were enriched for breast cancer risk SNP loci. We experimentally confirm SNVs to impact chromatin binding potential for ERα and its pioneer factor FOXA1. Finally, in the TCGA breast cancer cohort, we could confirm these variations to associate with differences in expression for the target gene. Cumulatively, we reveal a natural hierarchy of ERα-chromatin interactions in breast cancers within a highly heterogeneous inter-tumor ERα landscape, with the most-common shared regions being most active and affected by germline functional risk SNPs for breast cancer development.

The crucial role of epigenetic regulation in breast cancer anti-estrogen resistance: Current findings and future perspectives

Breast cancer (BC) cell de-sensitization to Tamoxifen (TAM) or other selective estrogen receptor (ER) modulators (SERM) is a complex process associated with BC heterogeneity and the transformation of ER signalling. The most influential resistance-related mechanisms include modifications in ER expression and gene regulation patterns. During TAM/SERM treatment, epigenetic mechanisms can effectively silence ER expression and facilitate the development of endocrine resistance. ER status is efficiently regulated by specific epigenetic tools including hypermethylation of CpG islands within ER promoters, increased histone deacetylase activity in the ER promoter, and/or translational repression by miRNAs. Over-methylation of the ER α gene (ESR1) promoter by DNA methyltransferases was associated with poor prognosis and indicated the development of resistance. Moreover, BC progression and spreading were marked by transformed chromatin remodelling, post-translational histone modifications, and expression of specific miRNAs and/or long non-coding RNAs. Therefore, targeted inhibition of histone acetyltransferases (e.g. MYST3), deacetylases (e.g. HDAC1), and/or demethylases (e.g. lysine-specific demethylase LSD1) was shown to recover and increase BC sensitivity to anti-estrogens. Indicated as a powerful molecular instrument, the administration of epigenetic drugs can regain ER expression along with the activation of tumour suppressor genes, which can in turn prevent selection of resistant cells and cancer stem cell survival. This review examines recent advances in the epigenetic regulation of endocrine drug resistance and evaluates novel anti-resistance strategies. Underlying molecular mechanisms of epigenetic regulation will be discussed, emphasising the utilization of epigenetic enzymes and their inhibitors to re-program irresponsive BCs.

Tumor cell intrinsic and extrinsic features predict prognosis in estrogen receptor positive breast cancer

Although estrogen-receptor-positive (ER+) breast cancer is generally associated with favorable prognosis, clinical outcome varies substantially among patients. Genomic assays have been developed and applied to predict patient prognosis for personalized treatment. We hypothesize that the recurrence risk of ER+ breast cancer patients is determined by both genomic mutations intrinsic to tumor cells and extrinsic immunological features in the tumor microenvironment. Based on the Cancer Genome Atlas (TCGA) breast cancer data, we identified the 72 most common genomic aberrations (including gene mutations and indels) in ER+ breast cancer and defined sample-specific scores that systematically characterized the deregulated pathways intrinsic to tumor cells. To further consider tumor cell extrinsic features, we calculated immune infiltration scores for six major immune cell types. Many individual intrinsic features are predictive of patient prognosis in ER+ breast cancer, and some of them achieved comparable accuracy with the Oncotype DX assay. In addition, statistical learning models that integrated these features predicts the recurrence risk of patients with significantly better performance than the Oncotype DX assay (our optimized random forest model AUC = 0.841, Oncotype DX model AUC = 0.792, p = 0.04). As a proof-of-concept, our study indicates the great potential of genomic and immunological features in prognostic prediction for improving breast cancer precision medicine. The framework introduced in this work can be readily applied to other cancers.

A systematic review on active sites and functions of PIM-1 protein

The Proviral Integration of Molony murine leukemia virus (PIM)-1 protein contributes to the solid cancers and hematologic malignancies, cell growth, proliferation, differentiation, migration, and other life activities. Many studies have related these functions to its molecular structure, subcellular localization and expression level. However, recognition of specific active sites and their effects on the activity of this constitutively active kinase is still a challenge. Based on the close relationship between its molecular structure and functional activity, this review covers the specific residues involved in the binding of ATP and different substrates in its catalytic domain. This review then elaborates on the relevant changes in protein conformation and cell functions after PIM-1 binds to different substrates. Therefore, this intensive study can improve the understanding of PIM-1-regulated signaling pathways by facilitating the discovery of its potential phosphorylation substrates.

Anastrozole Regulates Fatty Acid Synthase in Breast Cancer

Our previous matched case-control study of postmenopausal women with resected early-stage breast cancer revealed that only anastrozole, but not exemestane or letrozole, showed a significant association between the 6-month estrogen concentrations and risk of breast cancer. Anastrozole, but not exemestane or letrozole, is a ligand for estrogen receptor α. The mechanisms of endocrine resistance are heterogenous and with the new mechanism of anastrozole, we have found that treatment of anastrozole maintains fatty acid synthase (FASN) protein level by limiting the ubiquitin-mediated FASN degradation, leading to increased breast cancer cell growth. Mechanistically, anastrozole decreases the guided entry of tail-anchored proteins factor 4 (GET4) expression, resulting in decreased BCL2-associated athanogene cochaperone 6 (BAG6) complex activity, which in turn, prevents RNF126-mediated degradation of FASN. Increased FASN protein level can induce a negative feedback loop mediated by the MAPK pathway. High levels of FASN are associated with poor outcome only in patients with anastrozole-treated breast cancer, but not in patients treated with exemestane or letrozole. Repressing FASN causes regression of breast cancer cell growth. The anastrozole-FASN signaling pathway is eminently targetable in endocrine-resistant breast cancer.

An Optimized ChIP-Seq Protocol to Determine Chromatin Binding of Estrogen Receptor Beta

Estrogen regulates transcription through two nuclear receptors, ERα and ERβ, in a tissue and cellular-dependent manner. Both the receptors bind estrogen and activate transcription through direct or indirect interactions with DNA. Revealing their interactions with the chromatin is key to understanding their transcriptional activities and their biological functions. Chromatin-immunoprecipitation followed by sequencing (ChIP-Seq) is a powerful technique to map protein-DNA interactions at precise genomic locations. The genome-wide binding of ERα has been extensively studied. Similar studies of ERβ, however, have been more difficult, in part due to a lack of endogenous expression in cell lines and lack of specific antibodies. In this chapter, we provide an optimized stepwise ChIP protocol for a well-validated ERβ antibody, which is applicable for ChIP-Seq analysis of cell lines with exogenous expression of ERβ.

A gentle introduction to understanding preclinical data for cancer pharmaco-omic modeling

A central goal of precision oncology is to administer an optimal drug treatment to each cancer patient. A common preclinical approach to tackle this problem has been to characterize the tumors of patients at the molecular and drug response levels, and employ the resulting datasets for predictive in silico modeling (mostly using machine learning). Understanding how and why the different variants of these datasets are generated is an important component of this process. This review focuses on providing such introduction aimed at scientists with little previous exposure to this research area.

BMI1 activates P-glycoprotein via transcription repression of miR-3682-3p and enhances chemoresistance of bladder cancer cell

Chemoresistance is the most significant reason for the failure of cancer treatment following radical cystectomy. The response rate to the first-line chemotherapy of cisplatin and gemcitabine does not exceed 50%. In our previous research, elevated BMI1 (B-cell specific Moloney murine leukemia virus integration region 1) expression in bladder cancer conferred poor survival and was associated with chemoresistance. Herein, via analysis of The Cancer Genome Atlas database and validation of clinical samples, BMI1 was elevated in patients with bladder cancer resistant to cisplatin and gemcitabine, which conferred tumor relapse and progression. Consistently, BMI1 was markedly increased in the established cisplatin- and gemcitabine-resistant T24 cells (T24/DDP&GEM). Functionally, BMI1 overexpression dramatically promoted drug efflux, enhanced viability and decreased apoptosis of bladder cancer cells upon treatment with cisplatin or gemcitabine, whereas BMI1 downregulation reversed this effect. Mechanically, upon interaction with p53, BMI1 was recruited on the promoter of miR-3682-3p gene concomitant with an increase in the mono-ubiquitination of histone H2A lysine 119, leading to transcription repression of miR-3682-3p gene followed by derepression of ABCB1 (ATP binding cassette subfamily B member 1) gene. Moreover, suppression of P-glycoprotein by miR-3682-3p mimics or its inhibitor XR-9576, could significantly reverse chemoresistance of T24/DDP&GEM cells. These results provided a novel insight into a portion of the mechanism underlying BMI1-mediated chemoresistance in bladder cancer.

Synthetic lethality-mediated precision oncology via the tumor transcriptome

Precision oncology has made significant advances, mainly by targeting actionable mutations in cancer driver genes. Aiming to expand treatment opportunities, recent studies have begun to explore the utility of tumor transcriptome to guide patient treatment. Here, we introduce SELECT (synthetic lethality and rescue-mediated precision oncology via the transcriptome), a precision oncology framework harnessing genetic interactions to predict patient response to cancer therapy from the tumor transcriptome. SELECT is tested on a broad collection of 35 published targeted and immunotherapy clinical trials from 10 different cancer types. It is predictive of patients' response in 80% of these clinical trials and in the recent multi-arm WINTHER trial. The predictive signatures and the code are made publicly available for academic use, laying a basis for future prospective clinical studies.

A plug and play microfluidic platform for standardized sensitive low-input chromatin immunoprecipitation

Epigenetic profiling by chromatin immunoprecipitation followed by sequencing (ChIP-seq) has become a powerful tool for genome-wide identification of regulatory elements, for defining transcriptional regulatory networks, and for screening for biomarkers. However, the ChIP-seq protocol for low-input samples is laborious and time-consuming and suffers from experimental variation, resulting in poor reproducibility and low throughput. Although prototypic microfluidic ChIP-seq platforms have been developed, these are poorly transferable as they require sophisticated custom-made equipment and in-depth microfluidic and ChIP expertise, while lacking parallelization. To enable standardized, automated ChIP-seq profiling of low-input samples, we constructed microfluidic PDMS-based plates capable of performing 24 sensitive ChIP reactions within 30 min of hands-on time and 4.5 h of machine-running time. These disposable plates can be conveniently loaded into a widely available controller for pneumatics and thermocycling. In light of the plug and play (PnP) ChIP plates and workflow, we named our procedure PnP-ChIP-seq. We show high-quality ChIP-seq on hundreds to a few thousand of cells for all six post-translational histone modifications that are included in the International Human Epigenome Consortium set of reference epigenomes. PnP-ChIP-seq robustly detects epigenetic differences on promoters and enhancers between naive and more primed mouse embryonic stem cells (mESCs). Furthermore, we used our platform to generate epigenetic profiles of rare subpopulations of mESCs that resemble the two-cell stage of embryonic development. PnP-ChIP-seq allows nonexpert laboratories worldwide to conveniently run robust, standardized ChIP-seq, whereas its high throughput, consistency, and sensitivity pave the way toward large-scale profiling of precious sample types such as rare subpopulations of cells or biopsies.

Comparative analysis of the AIB1 interactome in breast cancer reveals MTA2 as a repressive partner which silences E-Cadherin to promote EMT and associates with a pro-metastatic phenotype

Steroid regulated cancer cells use nuclear receptors and associated regulatory proteins to orchestrate transcriptional networks to drive disease progression. In primary breast cancer, the coactivator AIB1 promotes estrogen receptor (ER) transcriptional activity to enhance cell proliferation. The function of the coactivator in ER+ metastasis however is not established. Here we describe AIB1 as a survival factor, regulator of pro-metastatic transcriptional pathways and a promising actionable target. Genomic alterations and functional expression of AIB1 associated with reduced disease-free survival in patients and enhanced metastatic capacity in novel CDX and PDX ex-vivo models of ER+ metastatic disease. Comparative analysis of the AIB1 interactome with complementary RNAseq characterized AIB1 as a transcriptional repressor. Specifically, we report that AIB1 interacts with MTA2 to form a repressive complex, inhibiting CDH1 (encoding E-cadherin) to promote EMT and drive progression. We further report that pharmacological and genetic inhibition of AIB1 demonstrates significant anti-proliferative activity in patient-derived models establishing AIB1 as a viable strategy to target endocrine resistant metastasis. This work defines a novel role for AIB1 in the regulation of EMT through transcriptional repression in advanced cancer cells with a considerable implication for prognosis and therapeutic interventions.

A comprehensive enhancer screen identifies TRAM2 as a key and novel mediator of YAP oncogenesis

BACKGROUND

Frequent activation of the co-transcriptional factor YAP is observed in a large number of solid tumors. Activated YAP associates with enhancer loci via TEAD4-DNA-binding protein and stimulates cancer aggressiveness. Although thousands of YAP/TEAD4 binding-sites are annotated, their functional importance is unknown. Here, we aim at further identification of enhancer elements that are required for YAP functions.

RESULTS

We first apply genome-wide ChIP profiling of YAP to systematically identify enhancers that are bound by YAP/TEAD4. Next, we implement a genetic approach to uncover functions of YAP/TEAD4-associated enhancers, demonstrate its robustness, and use it to reveal a network of enhancers required for YAP-mediated proliferation. We focus on EnhancerTRAM2, as its target gene TRAM2 shows the strongest expression-correlation with YAP activity in nearly all tumor types. Interestingly, TRAM2 phenocopies the YAP-induced cell proliferation, migration, and invasion phenotypes and correlates with poor patient survival. Mechanistically, we identify FSTL-1 as a major direct client of TRAM2 that is involved in these phenotypes. Thus, TRAM2 is a key novel mediator of YAP-induced oncogenic proliferation and cellular invasiveness.

CONCLUSIONS

YAP is a transcription co-factor that binds to thousands of enhancer loci and stimulates tumor aggressiveness. Using unbiased functional approaches, we dissect YAP enhancer network and characterize TRAM2 as a novel mediator of cellular proliferation, migration, and invasion. Our findings elucidate how YAP induces cancer aggressiveness and may assist diagnosis of cancer metastasis.

Estrogen Receptor on the move: Cistromic plasticity and its implications in breast cancer

Estrogen Receptor (ERα) is a hormone-driven transcription factor, critically involved in driving tumor cell proliferation in the vast majority of breast cancers (BCas). ERα binds the genome at cis-regulatory elements, dictating the expression of a large spectrum of responsive genes in 3D genomic space. While initial reports described a rather static ERα chromatin binding repertoire, we now know that ERα DNA interactions are highly versatile, altered in breast tumor development and progression, and deviate between tumors from patients with differential outcome. Multiple cellular signaling cascades are known to impinge on ERα genomic function, changing its cistrome to retarget the receptor to other regions of the genome and reprogram its impact on breast cell biology. This review describes the current state-of-the-art on which factors manipulate the ERα cistrome and how this alters the response to both endogenous and exogenous hormonal stimuli, ultimately impacting BCa cell progression and response to commonly used therapeutic interventions. Novel insights in ERα cistrome dynamics may pave the way for better patient diagnostics and the development of novel therapeutic interventions, ultimately improving cancer care and patient outcome.

Androgen receptor signalling in macrophages promotes TREM-1-mediated prostate cancer cell line migration and invasion

The androgen receptor (AR) is the master regulator of prostate cancer (PCa) development, and inhibition of AR signalling is the most effective PCa treatment. AR is expressed in PCa cells and also in the PCa-associated stroma, including infiltrating macrophages. Macrophages have a decisive function in PCa initiation and progression, but the role of AR in macrophages remains largely unexplored. Here, we show that AR signalling in the macrophage-like THP-1 cell line supports PCa cell line migration and invasion in culture via increased Triggering Receptor Expressed on Myeloid cells-1 (TREM-1) signalling and expression of its downstream cytokines. Moreover, AR signalling in THP-1 and monocyte-derived macrophages upregulates IL-10 and markers of tissue residency. In conclusion, our data suggest that AR signalling in macrophages may support PCa invasiveness, and blocking this process may constitute one mechanism of anti-androgen therapy.

Anti-androgen therapy inhibits prostate cancer (PC) progression, and is thought to act directly on cancer cells. Here the authors show that androgen receptor is expressed on normal and PC-associated macrophages, and its stimulation alters macrophage secretome to promote migration of cultured PC cell lines.

Methylation in HOX Clusters and Its Applications in Cancer Therapy

HOX genes are commonly known for their role in embryonic development, defining the positional identity of most structures along the anterior–posterior axis. In postembryonic life, HOX gene aberrant expression can affect several processes involved in tumorigenesis such as proliferation, apoptosis, migration and invasion. Epigenetic modifications are implicated in gene expression deregulation, and it is accepted that methylation events affecting HOX gene expression play crucial roles in tumorigenesis. In fact, specific methylation profiles in the HOX gene sequence or in HOX-associated histones are recognized as potential biomarkers in several cancers, helping in the prediction of disease outcomes and adding information for decisions regarding the patient’s treatment. The methylation of some HOX genes can be associated with chemotherapy resistance, and its identification may suggest the use of other treatment options. The use of epigenetic drugs affecting generalized or specific DNA methylation profiles, an approach that now deserves much attention, seems likely to be a promising weapon in cancer therapy in the near future. In this review, we summarize these topics, focusing particularly on how the regulation of epigenetic processes may be used in cancer therapy.

Endonuclease FEN1 Coregulates ERα Activity and Provides a Novel Drug Interface in Tamoxifen-Resistant Breast Cancer

Estrogen receptor α (ERα) is a key transcriptional regulator in the majority of breast cancers. ERα-positive patients are frequently treated with tamoxifen, but resistance is common. In this study, we refined a previously identified 111-gene outcome prediction-classifier, revealing FEN1 as the strongest determining factor in ERα-positive patient prognostication. FEN1 levels were predictive of outcome in tamoxifen-treated patients, and FEN1 played a causal role in ERα-driven cell growth. FEN1 impacted the transcriptional activity of ERα by facilitating coactivator recruitment to the ERα transcriptional complex. FEN1 blockade induced proteasome-mediated degradation of activated ERα, resulting in loss of ERα-driven gene expression and eradicated tumor cell proliferation. Finally, a high-throughput 465,195 compound screen identified a novel FEN1 inhibitor, which effectively blocked ERα function and inhibited proliferation of tamoxifen-resistant cell lines as well as ex vivo-cultured ERα-positive breast tumors. Collectively, these results provide therapeutic proof of principle for FEN1 blockade in tamoxifen-resistant breast cancer. SIGNIFICANCE: These findings show that pharmacologic inhibition of FEN1, which is predictive of outcome in tamoxifen-treated patients, effectively blocks ERα function and inhibits proliferation of tamoxifen-resistant tumor cells.

OSbrca: A Web Server for Breast Cancer Prognostic Biomarker Investigation With Massive Data From Tens of Cohorts

Potential prognostic mRNA biomarkers are exploited to assist in the clinical management and treatment of breast cancer, which is the first life-threatening tumor in women worldwide. However, it is technically challenging for untrained researchers to process high dimensional profiling data to screen and validate the potential prognostic values of genes of interests in multiple cohorts. Our aim is to develop an easy-to-use web server to facilitate the screening, developing, and evaluating of prognostic biomarkers in breast cancers. Herein, we collected more than 7,400 cases of breast cancer with gene expression profiles and clinical follow-up information from The Cancer Genome Atlas and Gene Expression Omnibus data, and built an Online consensus Survival analysis web server for Breast Cancers, abbreviated OSbrca, to generate the Kaplan–Meier survival plot with a hazard ratio and log rank P-value for given genes in an interactive way. To examine the performance of OSbrca, the prognostic potency of 128 previously published biomarkers of breast cancer was reassessed in OSbrca. In conclusion, it is highly valuable for biologists and clinicians to perform the preliminary assessment and validation of novel or putative prognostic biomarkers for breast cancers. OSbrca could be accessed at http://bioinfo.henu.edu.cn/BRCA/BRCAList.jsp.

A CRISPR-Cas9 screen identifies essential CTCF anchor sites for estrogen receptor-driven breast cancer cell proliferation

Estrogen receptor α (ERα) is an enhancer activating transcription factor, a key driver of breast cancer and a main target for cancer therapy. ERα-mediated gene regulation requires proper chromatin-conformation to facilitate interactions between ERα-bound enhancers and their target promoters. A major determinant of chromatin structure is the CCCTC-binding factor (CTCF), that dimerizes and together with cohesin stabilizes chromatin loops and forms the boundaries of topologically associated domains. However, whether CTCF-binding elements (CBEs) are essential for ERα-driven cell proliferation is unknown. To address this question in a global manner, we implemented a CRISPR-based functional genetic screen targeting CBEs located in the vicinity of ERα-bound enhancers. We identified four functional CBEs and demonstrated the role of one of them in inducing chromatin conformation changes in favor of activation of PREX1, a key ERα target gene in breast cancer. Indeed, high PREX1 expression is a bona-fide marker of ERα-dependency in cell lines, and is associated with good outcome after anti-hormonal treatment. Altogether, our data show that distinct CTCF-mediated chromatin structures are required for ERα- driven breast cancer cell proliferation.

Exogenous ERα Expression in the Mammary Epithelium Decreases Over Time and Does Not Contribute to p53-Deficient Mammary Tumor Formation in Mice

Approximately 75% of all breast cancers express the nuclear hormone receptor estrogen receptor α (ERα). However, the majority of mammary tumors from genetically engineered mouse models (GEMMs) are ERα-negative. To model ERα-positive breast cancer in mice, we exogenously introduced expression of mouse and human ERα in an existing GEMM of p53-deficient breast cancer. After initial ERα expression during mammary gland development, expression was reduced or lost in adult glands and p53-deficient mammary tumors. Chromatin immunoprecipitation (ChIP)-sequencing analysis of primary mouse mammary epithelial cells (MMECs) derived from these models, in which expression of the ERα constructs was induced in vitro, confirmed interaction of ERα with the DNA. In human breast and endometrial cancer, and also in healthy breast tissue, DNA binding of ERα is facilitated by the pioneer factor FOXA1. Surprisingly, the ERα binding sites identified in primary MMECs, but also in mouse mammary gland and uterus, showed an high enrichment of ERE motifs, but were devoid of Forkhead motifs. Furthermore, exogenous introduction of FOXA1 and GATA3 in ERα-expressing MMECs was not sufficient to promote ERα-responsiveness of these cells. Together, this suggests that species-specific differences in pioneer factor usage between mouse and human are dictated by the DNA sequence, resulting in ERα-dependencies in mice that are not FOXA1 driven. These species-specific differences in ERα-biology may limit the utility of mice for in vivo modeling of ERα-positive breast cancer.

Potential of epigenetic events in human thyroid cancer

Thyroid cancer remains the highest prevailing endocrine malignancy, and its incidence rate has progressively increased in the previous years. Above 95% of thyroid tumor are follicular cells types of carcinoma in which are considered invasive type of tumor. The pathogenesis and molecular mechanism of thyroid tumors are yet remains elucidated, in spite of activating RET, RAS and BRAF carcinogenesis have been well introduced. Nemours molecular alterations have been defined and have revealed promise for their diagnostic, prognostic and therapeutic capacity but still need further confirmation. Among different types of mechanisms, the current article reviews the importance of epigenetic modifications in thyroid cancer. Increasing data from previous reports demonstrate that acquired epigenetic abnormalities together with genetic changes plays an important role in alteration of gene expression patterns. Aberrant DNA methylation has been well known in the CpG regions and profile of microRNAs (mi-RNAs) expression also involved in cancer development. In addition, the gene expression through epigenetic control contribution to thyroid cancer is analyzed and it is semi considered in the clinic. However the epigenetic of the thyroid cancer is yet remains in its early stages, and it carries encouraging potential thyroid cancer detections in its early stages, assessment of prognosis and targeted cancer treatment.

Decoding the Inversion Symmetry Underlying Transcription Factor DNA-Binding Specificity and Functionality in the Genome

Understanding why a transcription factor (TF) binds to a specific DNA element in the genome and whether that binding event affects transcriptional output remains a great challenge. In this study, we demonstrate that TF binding in the genome follows inversion symmetry (IS). In addition, the specific DNA elements where TFs bind in the genome are determined by internal IS within the DNA element. These DNA-binding rules quantitatively define how TFs select the appropriate regulatory targets from a large number of similar DNA elements in the genome to elicit specific transcriptional and cellular responses. Importantly, we also demonstrate that these DNA-binding rules extend to DNA elements that do not support transcriptional activity. That is, the DNA-binding rules are obeyed, but the retention time of the TF at these non-functional DNA elements is not long enough to initiate and/or maintain transcription. We further demonstrate that IS is universal within the genome. Thus, IS is the DNA code that TFs use to interact with the genome and dictates (in conjunction with known DNA sequence constraints) which of those interactions are functionally active.

High expression of Linc00959 predicts poor prognosis in breast cancer

Background

Accumulating studies have focused on the oncogenic roles of the newly identified lncRNAs in human cancers. The aim of this study was to examine the expression pattern of Linc00959 in BC and to evaluate its biological role and clinical significance in prediction of prognosis.

Methods

Expression of Linc00959 was detected in 290 BC tissues by quantitative reverse-transcription polymerase chain reaction (qRT-PCR). We analyzed the relationship between Linc00959 expression and clinic pathological features of BC patients. The correlation was calculated by SPSS software.

Results

Our results revealed that Linc00959 expression was correlated with ER status (p = 0.005), PR status (p = 0.036), Ki67 (p = 0.025) and HER2 status (p = 0.009). The Kaplan–Meier survival curves indicated that the overall survival (OS) (p = 0.022) and relapse-free survival (RFS) (p = 0.002) were significantly poor in high Linc00959 expression BC patients (p = 0.023). Furthermore, the survival analysis by Cox regression showed that Linc00959 served as an independent prognostic marker in breast cancer (p = 0.004).

Conclusion

Our studies indicate that Linc00959 is significantly associated with poor prognosis and may represent a new marker of prognosis in breast cancer.

Electronic supplementary material

The online version of this article (10.1186/s12935-019-0748-7) contains supplementary material, which is available to authorized users.

Optimized ChIP-seq method facilitates transcription factor profiling in human tumors

This study presents an optimized ChIP-seq protocol to enhance transcription factor profiling in human tumours, enabling the analysis of highly challenging samples, including core needle biopsies.

Chromatin immunoprecipitation (ChIP)-seq analyses of transcription factors in clinical specimens are challenging due to the technical limitations and low quantities of starting material, often resulting in low enrichments and poor signal-to-noise ratio. Here, we present an optimized protocol for transcription factor ChIP-seq analyses in human tissue, yielding an ∼100% success rate for all transcription factors analyzed. As proof of concept and to illustrate general applicability of the approach, human tissue from the breast, prostate, and endometrial cancers were analyzed. In addition to standard formaldehyde fixation, disuccinimidyl glutarate was included in the procedure, greatly increasing data quality. To illustrate the sensitivity of the optimized protocol, we provide high-quality ChIP-seq data for three independent factors (AR, FOXA1, and H3K27ac) from a single core needle prostate cancer biopsy specimen. In summary, double-cross-linking strongly improved transcription factor ChIP-seq quality on human tumor samples, further facilitating and enhancing translational research on limited amounts of tissue.

Up-regulation of HOXB cluster genes are epigenetically regulated in tamoxifen-resistant MCF7 breast cancer cells

Tamoxifen (TAM) is commonly used to treat estrogen receptor (ER)-positive breast cancer. Despite the remarkable benefits, resistance to TAM presents a serious therapeutic challenge. Since several HOX transcription factors have been proposed as strong candidates in the development of resistance to TAM therapy in breast cancer, we generated an in vitro model of acquired TAM resistance using ER-positive MCF7 breast cancer cells (MCF7-TAMR), and analyzed the expression pattern and epigenetic states of HOX genes. HOXB cluster genes were uniquely up-regulated in MCF7-TAMR cells. Survival analysis of in slico data showed the correlation of high expression of HOXB genes with poor response to TAM in ER-positive breast cancer patients treated with TAM. Gain- and loss-of-function experiments showed that the overexpression of multi HOXB genes in MCF7 renders cancer cells more resistant to TAM, whereas the knockdown restores TAM sensitivity. Furthermore, activation of HOXB genes in MCF7-TAMR was associated with histone modifications, particularly the gain of H3K9ac. These findings imply that the activation of HOXB genes mediate the development of TAM resistance, and represent a target for development of new strategies to prevent or reverse TAM resistance.

Integrative epigenetic taxonomy of primary prostate cancer

The Androgen Receptor (AR) is the key-driving transcription factor in prostate cancer, tightly controlled by epigenetic regulation. To date, most epigenetic profiling has been performed in cell lines or limited tissue samples. Here, to comprehensively study the epigenetic landscape, we perform RNA-seq with ChIP-seq for AR and histone modification marks (H3K27ac, H3K4me3, H3K27me3) in 100 primary prostate carcinomas. Integrative molecular subtyping of the five data streams revealed three major subtypes of which two were clearly TMPRSS2-ERG dictated. Importantly, we identify a third subtype with low chromatin binding and activity of AR, but with high activity of FGF and WNT signaling. While positive for neuroendocrine-hallmark genes, these tumors were copy number-neutral with low mutational burden, significantly depleted for genes characteristic of poor-outcome associated luminal B-subtype. We present a unique resource on transcriptional and epigenetic control in prostate cancer, revealing tight control of gene regulation differentially dictated by AR over three subtypes.

The Androgen Receptor (AR) is the main driver of prostate cancer and functions in conjunction with chromatin modifications. Here, the authors comprehensively profile 100 primary prostate carcinomas by sequencing RNA transcripts in combination with ChIP-sequencing for AR and the active histone marks H3K27ac, H3K4me3 and repressive mark H3K27me3.

The Divergent Function of Androgen Receptor in Breast Cancer; Analysis of Steroid Mediators and Tumor Intracrinology

Androgen receptor (AR) is the most widely expressed steroid receptor protein in normal breast tissue and is detectable in approximately 90% of primary breast cancers and 75% of metastatic lesions. However, the role of AR in breast cancer development and progression is mired in controversy with evidence suggesting it can either inhibit or promote breast tumorigenesis. Studies have shown it to antagonize estrogen receptor alpha (ERα) DNA binding, thereby preventing pro-proliferative gene transcription; whilst others have demonstrated AR to take on the mantle of a pseudo ERα particularly in the setting of triple negative breast cancer. Evidence for a potentiating role of AR in the development of endocrine resistant breast cancer has also been mounting with reports associating high AR expression with poor response to endocrine treatment. The resurgence of interest into the function of AR in breast cancer has resulted in various emergent clinical trials evaluating anti-AR therapy and selective androgen receptor modulators in the treatment of advanced breast cancer. Trials have reported varied response rates dependent upon subtype with overall clinical benefit rates of ~19-29% for anti-androgen monotherapy, suggesting that with enhanced patient stratification AR could prove efficacious as a breast cancer therapy. Androgens and AR have been reported to facilitate tumor stemness in some cancers; a process which may be mediated through genomic or non-genomic actions of the AR, with the latter mechanism being relatively unexplored in breast cancer. Steroidogenic ligands of the AR are produced in females by the gonads and as sex-steroid precursors secreted from the adrenal glands. These androgens provide an abundant reservoir from which all estrogens are subsequently synthesized and their levels are undiminished in the event of standard hormonal therapeutic intervention in breast cancer. Steroid levels are known to be altered by lifestyle factors such as diet and exercise; understanding their potential role in dictating the function of AR in breast cancer development could therefore have wide-ranging effects in prevention and treatment of this disease. This review will outline the endogenous biochemical drivers of both genomic and non-genomic AR activation and how these may be modulated by current hormonal therapies.

The miR 495-UBE2C-ABCG2/ERCC1 axis reverses cisplatin resistance by downregulating drug resistance genes in cisplatin-resistant non-small cell lung cancer cells

Cisplatin (DDP) resistance has become the leading cause of mortality in non-small cell lung cancer (NSCLC). miRNA dysregulation significantly contributes to tumor progression. In this study, we found that miR-495 was significantly downregulated in lung cancer tissue specimens. This study aimed to elucidate the functions, direct target genes, and molecular mechanisms of miR-495 in lung cancer. miR-495 downregulated its substrate UBE2C through direct interaction with UBE2C 3′- untranslated region. UBE2C is a proto-oncogene activated in lung cancer; however, its role in chemotherapeutic resistance is unclear. Herein, UBE2C expression levels were higher in DDP-resistant NSCLC cells; this was associated with the proliferation, invasion, and DDP resistance in induced cisplatin-resistant NSCLC cells. Furthermore, epithelial–mesenchymal transitions (EMT) contributed to DDP resistance. Moreover, UBE2C knockdown downregulated vimentin. In contrast, E-cadherin was upregulated. Importantly, miR-495 and UBE2C were associated with cisplatin resistance. We attempted to evaluate their effects on cell proliferation and cisplatin resistance. We also performed EMT, cell migration, and invasion assays in DDP-resistant NSCLC cells overexpressing miR-495 and under-expressing UBE2C. Furthermore, in silico assays coupled with western blotting and luciferase assays revealed that UBE2C directly binds to the 5′-UTR of the drug-resistance genes ABCG2 and ERCC1. Furthermore, miR-495 downregulated ABCG2 and ERCC1 via regulation of UBE2C. Together, the present results indicate that the miR495-UBE2C-ABCG2/ERCC1 axis reverses DDP resistance via downregulation of anti-drug genes and reducing EMT in DDP-resistant NSCLC cells.

Loss of androgen receptor signaling in prostate cancer-associated fibroblasts (CAFs) promotes CCL2- and CXCL8-mediated cancer cell migration

Fibroblasts are abundantly present in the prostate tumor microenvironment (TME), including cancer‐associated fibroblasts (CAFs) which play a key role in cancer development. Androgen receptor (AR) signaling is the main driver of prostate cancer (PCa) progression, and stromal cells in the TME also express AR. High‐grade tumor and poor clinical outcome are associated with low AR expression in the TME, which suggests a protective role of AR signaling in the stroma against PCa development. However, the mechanism of this relation is not clear. In this study, we isolated AR‐expressing CAF‐like cells. Testosterone (R1881) exposure did not affect CAF‐like cell morphology, proliferation, or motility. PCa cell growth was not affected by culturing in medium from R1881‐exposed CAF‐like cells; however, migration of PCa cells was inhibited. AR chromatin immune precipitation sequencing (ChIP‐seq) was performed and motif search suggested that AR in CAF‐like cells bound the chromatin through AP‐1‐elements upon R1881 exposure, inducing enhancer‐mediated AR chromatin interactions. The vast majority of chromatin binding sites in CAF‐like cells were unique and not shared with AR sites observed in PCa cell lines or tumors. AR signaling in CAF‐like cells decreased expression of multiple cytokines; most notably CCL2 and CXCL8 and both cytokines increased migration of PCa cells. These results suggest direct paracrine regulation of PCa cell migration by CAFs through AR signaling.

Reprogramming of the estrogen responsive transcriptome contributes to tamoxifen-dependent protection against tumorigenesis in the p53 null mammary epithelial cells

The tumor suppressor gene p53 is frequently mutated in human breast cancer and is a marker for poor prognosis and resistance to chemotherapy. Transplantation of p53 null mouse mammary epithelium into syngeneic wild-type mice leads to normal mammary gland development followed by spontaneous mammary tumors that recapitulate many of the phenotypic, molecular and genetic features of human breast cancer. Transient exposure of p53 null mice to the anti-estrogen, tamoxifen leads to sustained and robust protection against tumor development. However the mechanism underlying this anti-tumor activity remains poorly understood. Here we demonstrate that transient exposure to tamoxifen leads to a reduction in mammary ductal side-branching and epithelial cell proliferation after tamoxifen withdrawal. Global gene expression analysis showed that transient tamoxifen exposure leads to persistent changes in the expression of a subset of estrogen regulated gene signatures in mammary epithelial cells (MECs). Among these was the protein tyrosine phosphatase, non-receptor type 5 (Ptpn5). We show that Ptpn5 is a novel tamoxifen regulated target gene which is upregulated in MECs after transient tamoxifen exposure and displays tumor suppressor activity in human breast cancer cells. Further, PTPN5 expression is strongly associated with good clinical outcome in tamoxifen treated human breast cancer patients suggesting that PTPN5 may represent a novel biomarker of tamoxifen response in human breast cancer.

Luminal lncRNAs Regulation by ERα-Controlled Enhancers in a Ligand-Independent Manner in Breast Cancer Cells

Estrogen receptor-α (ERα) is a ligand-inducible protein which mediates estrogenic hormones signaling and defines the luminal BC phenotype. Recently, we demonstrated that even in absence of ligands ERα (apoERα) binds chromatin sites where it regulates transcription of several protein-coding and lncRNA genes. Noteworthy, apoERα-regulated lncRNAs marginally overlap estrogen-induced transcripts, thus representing a new signature of luminal BC genes. By the analysis of H3K27ac enrichment in hormone-deprived MCF-7 cells, we defined a set of Super Enhancers (SEs) occupied by apoERα, including one mapped in proximity of the DSCAM-AS1 lncRNA gene. This represents a paradigm of apoERα activity since its expression is largely unaffected by estrogenic treatment, despite the fact that E2 increases ERα binding on DSCAM-AS1 promoter. We validated the enrichment of apoERα, p300, GATA3, FoxM1 and CTCF at both DSCAM-AS1 TSS and at its associated SE by ChIP-qPCR. Furthermore, by analyzing MCF-7 ChIA-PET data and by 3C assays, we confirmed long range chromatin interaction between the SE and the DSCAM-AS1 TSS. Interestingly, CTCF and p300 binding showed an enrichment in hormone-depleted medium and in the presence of ERα, elucidating the dynamics of the estrogen-independent regulation of DSCAM-AS1 expression. The analysis of this lncRNA provides a paradigm of transcriptional regulation of a luminal specific apoERα regulated lncRNA.

Characterizing steroid hormone receptor chromatin binding landscapes in male and female breast cancer

Male breast cancer (MBC) is rare and poorly characterized. Like the female counterpart, most MBCs are hormonally driven, but relapse after hormonal treatment is also noted. The pan-hormonal action of steroid hormonal receptors, including estrogen receptor alpha (ERα), androgen receptor (AR), progesterone receptor (PR), and glucocorticoid receptor (GR) in this understudied tumor type remains wholly unexamined. This study reveals genomic cross-talk of steroid hormone receptor action and interplay in human tumors, here in the context of MBC, in relation to the female disease and patient outcome. Here we report the characterization of human breast tumors of both genders for cistromic make-up of hormonal regulation in human tumors, revealing genome-wide chromatin binding landscapes of ERα, AR, PR, GR, FOXA1, and GATA3 and enhancer-enriched histone mark H3K4me1. We integrate these data with transcriptomics to reveal gender-selective and genomic location-specific hormone receptor actions, which associate with survival in MBC patients.

Neoadjuvant tamoxifen synchronizes ERα binding and gene expression profiles related to outcome and proliferation

Estrogen receptor alpha (ERα)-positive breast cancers are frequently treated with tamoxifen, but resistance is common. It remains elusive how tamoxifen resistance occurs and predictive biomarkers for treatment outcome are needed. Because most biomarker discovery studies are performed using pre-treatment surgical resections, the effects of tamoxifen therapy directly on the tumor cell in vivo remain unexamined. In this study, we assessed DNA copy number, gene expression profiles and ERα/chromatin binding landscapes on breast tumor specimens, both before and after neoadjuvant tamoxifen treatment. We observed neoadjuvant tamoxifen treatment synchronized ERα/chromatin interactions and downstream gene expression, indicating that hormonal therapy reduces inter-tumor molecular variability. ERα-synchronized sites are associated with dynamic FOXA1 action at these sites, which is under control of growth factor signaling. Genes associated with tamoxifen-synchronized sites are capable of differentiating patients for tamoxifen benefit. Due to the direct effects of therapeutics on ERα behavior and transcriptional output, our study highlights the added value of biomarker discovery studies after neoadjuvant drug exposure.

Structure-activity relationships and docking studies of synthetic 2-arylindole derivatives determined with aromatase and quinone reductase 1

In our ongoing effort of discovering anticancer and chemopreventive agents, a series of 2-arylindole derivatives were synthesized and evaluated toward aromatase and quinone reductase 1 (QR1). Biological evaluation revealed that several compounds (e.g., 2d, IC50 = 1.61 μM; 21, IC50 = 3.05 μM; and 27, IC50 = 3.34 μM) showed aromatase inhibitory activity with half maximal inhibitory concentration (IC50) values in the low micromolar concentrations. With regard to the QR1 induction activity, 11 exhibited the highest QR1 induction ratio (IR) with a low concentration to double activity (CD) value (IR = 8.34, CD = 2.75 μM), while 7 showed the most potent CD value of 1.12 μM. A dual acting compound 24 showed aromatase inhibition (IC50 = 9.00 μM) as well as QR1 induction (CD = 5.76 μM) activities. Computational docking studies using CDOCKER (Discovery Studio 3.5) provided insight in regard to the potential binding modes of 2-arylindoles within the aromatase active site. Predominantly, the 2-arylindoles preferred binding with the 2-aryl group toward a small hydrophobic pocket within the active site. The C-5 electron withdrawing group on indole was predicted to have an important role and formed a hydrogen bond with Ser478 (OH). Alternatively, meta-pyridyl analogs may orient with the pyridyl 3'-nitrogen coordinating with the heme group.

DNA Sequence Constraints Define Functionally Active Steroid Nuclear Receptor Binding Sites in Chromatin

Gene regulatory programs are encoded in the sequence of the DNA. Since the completion of the Human Genome Project, millions of gene regulatory elements have been identified in the human genome. Understanding how each of those sites functionally contributes to gene regulation, however, remains a challenge for nearly every field of biology. Transcription factors influence cell function by interpreting information contained within cis-regulatory elements in chromatin. Whereas chromatin immunoprecipitation-sequencing has been used to identify and map transcription factor-DNA interactions, it has been difficult to assign functionality to the binding sites identified. Thus, in this study, we probed the transcriptional activity, DNA-binding competence, and functional activity of select nuclear receptor mutants in cellular and animal model systems and used this information to define the sequence constraints of functional steroid nuclear receptor cis-regulatory elements. Analysis of the architecture within sNR chromatin interacting sites revealed that only a small fraction of all sNR chromatin-interacting events is associated with transcriptional output and that this functionality is restricted to elements that vary from the consensus palindromic elements by one or two nucleotides. These findings define the transcriptional grammar necessary to predict functionality from regulatory sequences, with a multitude of future implications.

Dissecting the genomic activity of a transcriptional regulator by the integrative analysis of omics data

In the study of genomic regulation, strategies to integrate the data produced by Next Generation Sequencing (NGS)-based technologies in a meaningful ensemble are eagerly awaited and must continuously evolve. Here, we describe an integrative strategy for the analysis of data generated by chromatin immunoprecipitation followed by NGS which combines algorithms for data overlap, normalization and epigenetic state analysis. The performance of our strategy is illustrated by presenting the analysis of data relative to the transcriptional regulator Estrogen Receptor alpha (ERα) in MCF-7 breast cancer cells and of Glucocorticoid Receptor (GR) in A549 lung cancer cells. We went through the definition of reference cistromes for different experimental contexts, the integration of data relative to co-regulators and the overlay of chromatin states as defined by epigenetic marks in MCF-7 cells. With our strategy, we identified novel features of estrogen-independent ERα activity, including FoxM1 interaction, eRNAs transcription and a peculiar ontology of connected genes.

Membrane and Nuclear Estrogen Receptor Alpha Actions: From Tissue Specificity to Medical Implications

Estrogen receptor alpha (ERα) has been recognized now for several decades as playing a key role in reproduction and exerting functions in numerous nonreproductive tissues. In this review, we attempt to summarize the in vitro studies that are the basis of our current understanding of the mechanisms of action of ERα as a nuclear receptor and the key roles played by its two activation functions (AFs) in its transcriptional activities. We then depict the consequences of the selective inactivation of these AFs in mouse models, focusing on the prominent roles played by ERα in the reproductive tract and in the vascular system. Evidence has accumulated over the two last decades that ERα is also associated with the plasma membrane and activates non-nuclear signaling from this site. These rapid/nongenomic/membrane-initiated steroid signals (MISS) have been characterized in a variety of cell lines, and in particular in endothelial cells. The development of selective pharmacological tools that specifically activate MISS and the generation of mice expressing an ERα protein impeded for membrane localization have begun to unravel the physiological role of MISS in vivo. Finally, we discuss novel perspectives for the design of tissue-selective ER modulators based on the integration of the physiological and pathophysiological roles of MISS actions of estrogens.

A review of estrogen receptor/androgen receptor genomics in male breast cancer

Male breast cancer is a rare disease, of which little is known. In contrast to female breast cancer, the very vast majority of all cases are positive for estrogen receptor alpha (ERα), implicating the function of this steroid hormone receptor in tumor development and progression. Consequently, adjuvant treatment of male breast cancer revolves around inhibition of ERα. In addition, the androgen receptor (AR) gradually receives more attention as a relevant novel target in breast cancer treatment. Importantly, the rationale of treatment decision making is strongly based on parallels with female breast cancer. Yet, prognostic indicators are not necessarily the same in breast cancer between both genders, complicating translatability of knowledge developed in female breast cancer toward male patients. Even though ERα and AR are expressed both in female and male disease, are the genomic functions of both steroid hormone receptors conserved between genders? Recent studies have reported on mutational and epigenetic similarities and differences between male and female breast cancer, further suggesting that some features are strongly conserved between the two diseases, whereas others are not. This review critically discusses the recent developments in the study of male breast cancer in relation to ERα and AR action and highlights the potential future studies to further elucidate the genomic regulation of this rare disease.

Estrogen receptor α wields treatment-specific enhancers between morphologically similar endometrial tumors

The DNA-binding sites of estrogen receptor α (ERα) show great plasticity under the control of hormones and endocrine therapy. Tamoxifen is a widely applied therapy in breast cancer that affects ERα interactions with coregulators and shifts the DNA-binding signature of ERα upon prolonged exposure in breast cancer. Although tamoxifen inhibits the progression of breast cancer, it increases the risk of endometrial cancer in postmenopausal women. We therefore asked whether the DNA-binding signature of ERα differs between endometrial tumors that arise in the presence or absence of tamoxifen, indicating divergent enhancer activity for tumors that develop in different endocrine milieus. Using ChIP sequencing (ChIP-seq), we compared the ERα profiles of 10 endometrial tumors from tamoxifen users with those of six endometrial tumors from nonusers and integrated these results with the transcriptomic data of 47 endometrial tumors from tamoxifen users and 64 endometrial tumors from nonusers. The ERα-binding sites in tamoxifen-associated endometrial tumors differed from those in the tumors from nonusers and had distinct underlying DNA sequences and divergent enhancer activity as marked by histone 3 containing the acetylated lysine 27 (H3K27ac). Because tamoxifen acts as an agonist in the postmenopausal endometrium, similar to estrogen in the breast, we compared ERα sites in tamoxifen-associated endometrial cancers with publicly available ERα ChIP-seq data in breast tumors and found a striking resemblance in the ERα patterns of the two tissue types. Our study highlights the divergence between endometrial tumors that arise in different hormonal conditions and shows that ERα enhancer use in human cancer differs in the presence of nonphysiological endocrine stimuli.

Genome-wide epigenomic profiling for biomarker discovery

A myriad of diseases is caused or characterized by alteration of epigenetic patterns, including changes in DNA methylation, post-translational histone modifications, or chromatin structure. These changes of the epigenome represent a highly interesting layer of information for disease stratification and for personalized medicine. Traditionally, epigenomic profiling required large amounts of cells, which are rarely available with clinical samples. Also, the cellular heterogeneity complicates analysis when profiling clinical samples for unbiased genome-wide biomarker discovery. Recent years saw great progress in miniaturization of genome-wide epigenomic profiling, enabling large-scale epigenetic biomarker screens for disease diagnosis, prognosis, and stratification on patient-derived samples. All main genome-wide profiling technologies have now been scaled down and/or are compatible with single-cell readout, including: (i) Bisulfite sequencing to determine DNA methylation at base-pair resolution, (ii) ChIP-Seq to identify protein binding sites on the genome, (iii) DNaseI-Seq/ATAC-Seq to profile open chromatin, and (iv) 4C-Seq and HiC-Seq to determine the spatial organization of chromosomes. In this review we provide an overview of current genome-wide epigenomic profiling technologies and main technological advances that allowed miniaturization of these assays down to single-cell level. For each of these technologies we evaluate their application for future biomarker discovery. We will focus on (i) compatibility of these technologies with methods used for clinical sample preservation, including methods used by biobanks that store large numbers of patient samples, and (ii) automation of these technologies for robust sample preparation and increased throughput.

Androgen receptor profiling predicts prostate cancer outcome

Prostate cancer is the second most prevalent malignancy in men. Biomarkers for outcome prediction are urgently needed, so that high-risk patients could be monitored more closely postoperatively. To identify prognostic markers and to determine causal players in prostate cancer progression, we assessed changes in chromatin state during tumor development and progression. Based on this, we assessed genomewide androgen receptor/chromatin binding and identified a distinct androgen receptor/chromatin binding profile between primary prostate cancers and tumors with an acquired resistance to therapy. These differential androgen receptor/chromatin interactions dictated expression of a distinct gene signature with strong prognostic potential. Further refinement of the signature provided us with a concise list of nine genes that hallmark prostate cancer outcome in multiple independent validation series. In this report, we identified a novel gene expression signature for prostate cancer outcome through generation of multilevel genomic data on chromatin accessibility and transcriptional regulation and integration with publically available transcriptomic and clinical datastreams. By combining existing technologies, we propose a novel pipeline for biomarker discovery that is easily implementable in other fields of oncology.

The Estrogen Receptor α-Cistrome Beyond Breast Cancer

Although many tissues express estrogen receptor (ER)α, most studies focus on breast cancer where ERα occupies just a small fraction of its total repertoire of potential DNA-binding sites, based on sequence. This raises the question: Can ERα occupy these other potential binding sites in a different context? Ligands, splice variants, posttranslational modifications, and acquired mutations of ERα affect its conformation, which may alter chromatin interactions. To date, literature describes the DNA-binding sites of ERα (the ERα cistrome) in breast, endometrium, liver, and bone, in which the receptor mainly binds to enhancers. Chromosomal boundaries provide distinct areas for dynamic gene regulation between tissues, where the usage of enhancers deviates. Interactions of ERα with enhancers and its transcriptional complex depend on the proteome, which differs per cell type. This review discusses the biological variables that influence ERα cistromics, using reports from human specimens, cell lines, and mouse tissues, to assess whether ERα genomics in breast cancer can be translated to other tissue types.

LRG1 mRNA expression in breast cancer associates with PIK3CA genotype and with aromatase inhibitor therapy outcome

BACKGROUND

PIK3CA is the most frequent somatic mutated oncogene in estrogen receptor (ER) positive breast cancer. We previously observed an association between PIK3CA genotype and aromatase inhibitors (AI) treatment outcome. This study now evaluates whether expression of mRNAs and miRs are linked to PIK3CA genotype and are independently related to AI therapy response in order to define potential expressed biomarkers for treatment outcome.

MATERIALS AND METHODS

The miR and mRNA expression levels were evaluated for their relationship with the PIK3CA genotype in two breast tumor datasets, i.e. 286 luminal cancers from the TCGA consortium and our set of 84 ER positive primary tumors of metastatic breast cancer patients who received first line AI. BRB Array tools class comparison was performed to define miRs and mRNAs whose expression associate with PIK3CA exon 9 and 20 status. Spearman correlations established miR-mRNA pairs and mRNAs with related expression. Next, a third dataset of 25 breast cancer patients receiving neo-adjuvant letrozole was evaluated, to compare expression levels of identified miRs and mRNAs in biopsies before and after treatment. Finally, to identify potential biomarkers miR and mRNA levels were related with overall survival (OS) and progression free survival (PFS) after first-line AI therapy.

RESULTS

Expression of 3 miRs (miR-449a, miR-205-5p, miR-301a-3p) and 9 mRNAs (CCNO, FAM81B, LRG1, NEK10, PLCL1, PGR, SERPINA3, SORBS2, VTCN1) was related to the PIK3CA status in both datasets. All except miR-301a-3p had an increased expression in tumors with PIK3CA mutations. Validation in a publicly available dataset showed that LRG1, PGR, and SERPINA3 levels were decreased after neo-adjuvant AI-treatment. Six miR-mRNA pairs correlated significantly and stepdown analysis of all 12 factors revealed 3 mRNAs (PLCL1, LRG1, FAM81B) related to PFS. Further analyses showed LRG1 and PLCL1 expression to be unrelated with luminal subtype and to associate with OS and with PFS, the latter independent from traditional predictive factors.

CONCLUSION

We showed in two datasets of ER positive and luminal breast tumors that the expression of 3 miRs and 9 mRNAs associate with the PIK3CA status. Expression of LRG1 is independent of luminal (A or B) subtype, decreased after neo-adjuvant AI-treatment, and is proposed as potential biomarker for AI therapy outcome.

Therapeutic Implications of Cellular Heterogeneity and Plasticity in Breast Cancer

Cellular heterogeneity represents one of the greatest challenges in cancer therapeutics. In many malignancies, this heterogeneity is generated during tumor evolution through a combination of genetic alterations and epigenetic events that recapitulate normal developmental processes including stem cell self-renewal and differentiation. Many, if not most, tumors display similar hierarchal organization, at the apex of which are "stem-like cells" that drive tumor growth, mediate metastasis, and contribute to treatment resistance. Using breast cancer as a model, we discuss how an improved understanding of tumor cellular heterogeneity and plasticity may lead to development of more effective therapeutic strategies.

Endocrine therapy resistance in estrogen receptor (ER)-positive breast cancer

Estrogen receptor (ER)-positive breast cancer represents the majority (∼70%) of all breast malignancies. In this subgroup of breast cancers, endocrine therapies are effective both in the adjuvant and recurrent settings, although resistance remains a major issue. Several high-throughput approaches have been used to elucidate mechanisms of resistance and to derive potential predictive markers or alternative therapies. In this review, we cover the state-of-the-art of endocrine-resistance biomarker discovery with regard to the latest technological developments, and discuss current opportunities and restrictions for their implementation into a clinical setting.

Analysis of Paired Primary-Metastatic Hormone-Receptor Positive Breast Tumors (HRPBC) Uncovers Potential Novel Drivers of Hormonal Resistance

We sought to identify genetic variants associated with disease relapse and failure to hormonal treatment in hormone-receptor positive breast cancer (HRPBC). We analyzed a series of HRPBC with distant relapse, by sequencing pairs (n = 11) of tumors (primary and metastases) at >800X. Comparative genomic hybridization was performed as well. Top hits, based on the frequency of alteration and severity of the changes, were tested in the TCGA series. Genes determining the most parsimonious prognostic signature were studied for their functional role in vitro, by performing cell growth assays in hormonal-deprivation conditions, a setting that mimics treatment with aromatase inhibitors. Severe alterations were recurrently found in 18 genes in the pairs. However, only MYC, DNAH5, CSFR1, EPHA7, ARID1B, and KMT2C preserved an independent prognosis impact and/or showed a significantly different incidence of alterations between relapsed and non-relapsed cases in the TCGA series. The signature composed of MYC, KMT2C, and EPHA7 best discriminated the clinical course, (overall survival 90,7 vs. 144,5 months; p = 0.0001). Having an alteration in any of the genes of the signature implied a hazard ratio of death of 3.25 (p<0.0001), and early relapse during the adjuvant hormonal treatment. The presence of the D348N mutation in KMT2C and/or the T666I mutation in the kinase domain of EPHA7 conferred hormonal resistance in vitro. Novel inactivating mutations in KMT2C and EPHA7, which confer hormonal resistance, are linked to adverse clinical course in HRPBC.