Epigenetic reprogramming at estrogen-receptor binding sites alters 3D chromatin landscape in endocrine-resistant breast cancer

17α-Ethynylestradiol alters testicular epigenetic profiles and histone-to-protamine exchange in mice

Spermatogenesis starts with the onset of puberty within the seminiferous epithelium of the testes. It is a complex process under intricate control of the endocrine system. Physiological regulations by steroid hormones in general and by estrogens in particular are due to their chemical nature prone to be disrupted by exogenous factors acting as endocrine disruptors (EDs). 17α-Ethynylestradiol (EE2) is an environmental pollutant with a confirmed ED activity and a well-known effect on spermatogenesis and chromatin remodeling in haploid germ cells. The aim of our study was to assess possible effects of two doses (2.5ng/ml; 2.5 μg/ml) of EE2 on both histone-to-protamine exchange and epigenetic profiles during spermatogenesis performing a multi/transgenerational study in mice. Our results demonstrated an impaired histone-to-protamine exchange with a significantly higher histone retention in sperm nuclei of exposed animals, when this process was accompanied by the changes of histone post-translational modifications (PTMs) abundancies with a prominent effect on H3K9Ac and partial changes in protamine 1 promoter methylation status. Furthermore, individual changes in molecular phenotypes were partially transmitted to subsequent generations, when no direct trans-generational effect was observed. Finally, the uncovered specific localization of the histone retention in sperm nuclei and their specific PTMs profile after EE2 exposure may indicate an estrogenic effect on sperm motility and early embryonic development via epigenetic mechanisms.

Epigenetic modulation via the C-terminal tail of H2A.Z

H2A.Z-nucleosomes are present in both euchromatin and heterochromatin and it has proven difficult to interpret their disparate roles in the context of their stability features. Using an in situ assay of nucleosome stability and DT40 cells expressing engineered forms of the histone variant we show that native H2A.Z, but not C-terminally truncated H2A.Z (H2A.Z∆C), is released from nucleosomes of peripheral heterochromatin at unusually high salt concentrations. H2A.Z and H3K9me3 landscapes are reorganized in H2A.Z∆C-nuclei and overall sensitivity of chromatin to nucleases is increased. These tail-dependent differences are recapitulated upon treatment of HeLa nuclei with the H2A.Z-tail-peptide (C9), with MNase sensitivity being increased genome-wide. Fluorescence correlation spectroscopy revealed C9 binding to reconstituted nucleosomes. When introduced into live cells, C9 elicited chromatin reorganization, overall nucleosome destabilization and changes in gene expression. Thus, H2A.Z-nucleosomes influence global chromatin architecture in a tail-dependent manner, what can be modulated by introducing the tail-peptide into live cells.

Unexpected nuclear hormone receptor and chromatin dynamics regulate estrous cycle dependent gene expression

Chromatin changes in response to estrogen and progesterone are well established in cultured cells, but how they control gene expression under physiological conditions is largely unknown. To address this question, we examined in vivo estrous cycle dynamics of mouse uterus hormone receptor occupancy, chromatin accessibility and chromatin structure by combining RNA-seq, ATAC-seq, HiC-seq and ChIP-seq. Two estrous cycle stages were chosen for these analyses, diestrus (highest estrogen) and estrus (highest progesterone). Unexpectedly, rather than alternating with each other, estrogen receptor alpha (ERα) and progesterone receptor (PGR) were co-bound during diestrus and lost during estrus. Motif analysis of open chromatin followed by hypoxia inducible factor 2A (HIF2A) ChIP-seq and conditional uterine deletion of this transcription factor revealed a novel role for HIF2A in regulating diestrus gene expression patterns that were independent of either ERα or PGR binding. Proteins in complex with ERα included PGR and cohesin, only during diestrus. Combined with HiC-seq analyses, we demonstrate that complex chromatin architecture changes including enhancer switching are coordinated with ERα and PGR co-binding during diestrus and non-hormone receptor transcription factors such as HIF2A during estrus to regulate most differential gene expression across the estrous cycle.

Induction of the TEAD Co-activator VGLL1 by Estrogen Receptor-Targeted Therapy Drives Resistance in Breast Cancer

Resistance to endocrine therapies (ET) is common in estrogen receptor (ER) positive breast cancer, and most relapsed patients die with ET-resistant disease. While genetic mutations provide explanations for some relapses, mechanisms of resistance remain undefined in many cases. Drug-induced epigenetic reprogramming has been shown to provide possible routes to resistance. By analyzing histone H3 lysine 27 acetylation (H3K27ac) profiles and transcriptional reprogramming in models of ET resistance, we discovered that selective ER degraders (SERDs), such as fulvestrant, promote expression of VGLL1, a co-activator for TEAD transcription factors. VGLL1, acting via TEADs, promoted expression of genes that drive growth of fulvestrant-resistant breast cancer cells. Pharmacological disruption of VGLL1/TEAD4 interaction inhibited VGLL1/TEAD-induced transcriptional programs to prevent growth of resistant cells. EGFR was among the VGLL1/TEAD-regulated genes, and VGLL1-directed EGFR upregulation sensitized fulvestrant-resistant breast cancer cells to EGFR inhibitors. Taken together, these findings identify VGLL1 as a transcriptional driver in ET resistance and advance therapeutic possibilities for relapsed ER+ breast cancer patients.

Deciphering Pre-existing and Induced 3D Genome Architecture Changes involved in Constricted Melanoma Migration

Metastatic cancer cells traverse constricted spaces that exert forces on their nucleus and the genomic contents within. Cancerous tumors are highly heterogeneous and not all cells within them can achieve such a feat. Here, we investigated what initial genome architecture characteristics favor the constricted migratory ability of cancer cells and which arise only after passage through multiple constrictions. We identified a cell surface protein (ITGB4) whose expression correlates with increased initial constricted migration ability in human melanoma A375 cells. Sorting out this subpopulation allowed us to identify cellular and nuclear features that pre-exist and favor migration, as well as alterations that only appear after cells have passed through constrictions. We identified specific genomic regions that experienced altered genome spatial compartment profiles only after constricted migration. Our study reveals 3D genome structure contributions to both selection and induction mechanisms of cell fate change during cancer metastasis.

Machine learning enables pan-cancer identification of mutational hotspots at persistent CTCF binding sites

CCCTC-binding factor (CTCF) is an insulator protein that binds to a highly conserved DNA motif and facilitates regulation of three-dimensional (3D) nuclear architecture and transcription. CTCF binding sites (CTCF-BSs) reside in non-coding DNA and are frequently mutated in cancer. Our previous study identified a small subclass of CTCF-BSs that are resistant to CTCF knock down, termed persistent CTCF binding sites (P-CTCF-BSs). P-CTCF-BSs show high binding conservation and potentially regulate cell-type constitutive 3D chromatin architecture. Here, using ICGC sequencing data we made the striking observation that P-CTCF-BSs display a highly elevated mutation rate in breast and prostate cancer when compared to all CTCF-BSs. To address whether P-CTCF-BS mutations are also enriched in other cell-types, we developed CTCF-INSITE-a tool utilising machine learning to predict persistence based on genetic and epigenetic features of experimentally-determined P-CTCF-BSs. Notably, predicted P-CTCF-BSs also show a significantly elevated mutational burden in all 12 cancer-types tested. Enrichment was even stronger for P-CTCF-BS mutations with predicted functional impact to CTCF binding and chromatin looping. Using in vitro binding assays we validated that P-CTCF-BS cancer mutations, predicted to be disruptive, indeed reduced CTCF binding. Together this study reveals a new subclass of cancer specific CTCF-BS DNA mutations and provides insights into their importance in genome organization in a pan-cancer setting.

The role of HOTAIR in the modulation of resistance to anticancer therapy

Leading anti-tumour therapeutic strategies typically involve surgery and radiotherapy for locally advanced (non-metastatic) cancers, while hormone therapy, chemotherapy, and molecular targeted therapy are the current treatment options for metastatic cancer. Despite the initially high sensitivity rate to anticancer therapies, a large number of patients develop resistance, leading to a poor prognosis. The mechanisms related to drug resistance are highly complex, and long non-coding RNAs appear to play a crucial role in these processes. Among these, the lncRNA homeobox transcript antisense intergenic RNA (HOTAIR), widely implicated in cancer initiation and progression, likewise plays a significant role in anticancer drug resistance. It can modulate cell activities such as proliferation, apoptosis, hypoxia, autophagy, as well as epithelial-mesenchymal transition, thereby contributing to the development of resistant tumour cells. In this manuscript, we describe different mechanisms of antitumor drug resistance in which HOTAIR is involved and suggest its potential as a therapeutic predictive biomarker for the management of cancer patients.

Oestrogen receptor positive breast cancer and its embedded mechanism: breast cancer resistance to conventional drugs and related therapies, a review

Traditional medication and alternative therapies have long been used to treat breast cancer. One of the main problems with current treatments is that there is an increase in drug resistance in the cancer cells owing to genetic differences such as mutational changes, epigenetic changes and miRNA (microRNA) alterations such as miR-1246, miR-298, miR-27b and miR-33a, along with epigenetic modifications, such as Histone3 acetylation and CCCTC-Binding Factor (CTCF) hypermethylation for drug resistance in breast cancer cell lines. Certain forms of conventional drug resistance have been linked to genetic changes in genes such as ABCB1, AKT, S100A8/A9, TAGLN2 and NPM. This review aims to explore the current approaches to counter breast cancer, the action mechanism, along with novel therapeutic methods endowing potential drug resistance. The investigation of novel therapeutic approaches sheds light on the phenomenon of drug resistance including genetic variations that impact distinct forms of oestrogen receptor (ER) cancer, genetic changes, epigenetics-reported resistance and their identification in patients. Long-term effective therapy for breast cancer includes selective oestrogen receptor modulators, selective oestrogen receptor degraders and genetic variations, such as mutations in nuclear genes, epigenetic modifications and miRNA alterations in target proteins. Novel research addressing combinational therapies including maytansine, photodynamic therapy, guajadiol, talazoparib, COX2 inhibitors and miRNA 1246 inhibitors have been developed to improve patient survival rates.

Epigenetic deregulation in breast cancer microenvironment: Implications for tumor progression and therapeutic strategies

Breast cancer comprises a substantial proportion of cancer diagnoses in women and is a primary cause of cancer-related mortality. While hormone-responsive cases generally have a favorable prognosis, the aggressive nature of triple-negative breast cancer presents challenges, with intrinsic resistance to established treatments being a persistent issue. The complexity intensifies with the emergence of acquired resistance, further complicating the management of breast cancer. Epigenetic changes, encompassing DNA methylation, histone and RNA modifications, and non-coding RNAs, are acknowledged as crucial contributors to the heterogeneity of breast cancer. The unique epigenetic landscape harbored by each cellular component within the tumor microenvironment (TME) adds great diversity to the intricate regulations which influence therapeutic responses. The TME, a sophisticated ecosystem of cellular and non-cellular elements interacting with tumor cells, establishes an immunosuppressive microenvironment and fuels processes such as tumor growth, angiogenesis, and extracellular matrix remodeling. These factors contribute to challenging conditions in cancer treatment by fostering a hypoxic environment, inducing metabolic stress, and creating physical barriers to drug delivery. This article delves into the complex connections between breast cancer treatment response, underlying epigenetic changes, and vital interactions within the TME. To restore sensitivity to treatment, it emphasizes the need for combination therapies considering epigenetic changes specific to individual members of the TME. Recognizing the pivotal role of epigenetics in drug resistance and comprehending the specificities of breast TME is essential for devising more effective therapeutic strategies. The development of reliable biomarkers for patient stratification will facilitate tailored and precise treatment approaches.

Epigenetic mechanisms of cancer progression and therapy resistance in estrogen-receptor (ER+) breast cancer

Estrogen receptor-positive (ER+) breast cancer is the most frequent breast cancer subtype. Agents targeting the ER signaling pathway have been successful in reducing mortality from breast cancer for decades. However, mechanisms of resistance to these treatments arise, especially in the metastatic setting. Recently, it has been recognized that epigenetic dysregulation is a common feature that facilitates the acquisition of cancer hallmarks across cancer types, including ER+ breast cancer. Alterations in epigenetic regulators and transcription factors (TF) coupled with changes to the chromatin landscape have been found to orchestrate breast oncogenesis, metastasis, and the development of a resistant phenotype. Here, we review recent advances in our understanding of how the epigenome dictates breast cancer tumorigenesis and resistance to targeted therapies and discuss novel therapeutic interventions for overcoming resistance.

Near-Infrared Light-Driven Nanoparticles for Cancer Photoimmunotherapy by Synergizing Immune Cell Death and Epigenetic Regulation

Histone deacetylases (HDACs) are a class of epigenetic enzymes that are closely related to tumorigenesis and suppress the expression of tumor suppressor genes. Whereas the HDACs inhibitors can release DNA into the cytoplasm and trigger innate immunity. However, the high density of chromatin limits DNA damage and release. In this study, suitable nanosized CycNHOH NPs (150 nm) and CypNHOH NPs (85 nm) efficiently accumulate at the tumor site due to the enhanced permeability and retention (EPR) effect. In addition, robust single-linear oxygen generation and good photothermal conversion efficiency under NIR laser irradiation accelerated the DNA damage process. By effectively initiating immune cell death, CypNHOH NPs activated both innate and adaptive immunity by maturing dendritic cells, infiltrating tumors with natural killer cells, and activating cytotoxic T lymphocytes, which offer a fresh perspective for the development of photo-immunotherapy.

Proteasome Inhibition Reprograms Chromatin Landscape in Breast Cancer

The 26S proteasome is the major protein degradation machinery in cells. Cancer cells use the proteasome to modulate gene expression networks that promote tumor growth. Proteasome inhibitors have emerged as effective cancer therapeutics, but how they work mechanistically remains unclear. Here, using integrative genomic analysis, we discovered unexpected reprogramming of the chromatin landscape and RNA polymerase II (RNAPII) transcription initiation in breast cancer cells treated with the proteasome inhibitor MG132. The cells acquired dynamic changes in chromatin accessibility at specific genomic loci termed differentially open chromatin regions (DOCR). DOCRs with decreased accessibility were promoter proximal and exhibited unique chromatin architecture associated with divergent RNAPII transcription. Conversely, DOCRs with increased accessibility were primarily distal to transcription start sites and enriched in oncogenic superenhancers predominantly accessible in non-basal breast tumor subtypes. These findings describe the mechanisms by which the proteasome modulates the expression of gene networks intrinsic to breast cancer biology.

SIGNIFICANCE

Our study provides a strong basis for understanding the mechanisms by which proteasome inhibitors exert anticancer effects. We find open chromatin regions that change during proteasome inhibition, are typically accessible in non-basal breast cancers.

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.

Hi-C profiling in tissues reveals 3D chromatin-regulated breast tumor heterogeneity and tumor-specific looping-mediated biological pathways

Current knowledge in three-dimensional (3D) chromatin regulation in normal and disease states was mostly accumulated through Hi-C profiling in in vitro cell culture system. The limitations include failing to recapitulate disease-specific physiological properties and often lacking clinically relevant disease microenvironment. In this study, we conduct tissue-specific Hi-C profiling in a pilot cohort of 12 breast tissues comprising of two normal tissues (NTs) and ten ER+ breast tumor tissues (TTs) including five primary tumors (PTs), and five tamoxifen-treated recurrent tumors (RTs). We find largely preserved compartments, highly heterogeneous topological associated domains (TADs) and intensively variable chromatin loops among breast tumors, demonstrating 3D chromatin-regulated breast tumor heterogeneity. Further cross-examination identifies RT-specific looping-mediated biological pathways and suggests CA2, an enhancer-promoter looping (EPL)-mediated target gene within the bicarbonate transport metabolism pathway, might play a role in driving the tamoxifen resistance. Remarkably, the inhibition of CA2 not only impedes tumor growth both in vitro and in vivo , but also reverses chromatin looping. Our study thus yields significant mechanistic insights into the role and clinical relevance of 3D chromatin architecture in breast cancer endocrine resistance.

Characterisation and reproducibility of the HumanMethylationEPIC v2.0 BeadChip for DNA methylation profiling

BACKGROUND

The Illumina family of Infinium Methylation BeadChip microarrays has been widely used over the last 15 years for genome-wide DNA methylation profiling, including large-scale and population-based studies, due to their ease of use and cost effectiveness. Succeeding the popular HumanMethylationEPIC BeadChip (EPICv1), the recently released Infinium MethylationEPIC v2.0 BeadChip (EPICv2) claims to extend genomic coverage to more than 935,000 CpG sites. Here, we comprehensively characterise the reproducibility, reliability and annotation of the EPICv2 array, based on bioinformatic analysis of both manifest data and new EPICv2 data from diverse biological samples.

RESULTS

We find a high degree of reproducibility with EPICv1, evidenced by comparable sensitivity and precision from empirical cross-platform comparison incorporating whole genome bisulphite sequencing (WGBS), and high correlation between technical sample replicates, including between samples with DNA input levels below the manufacturer's recommendation. We provide a full assessment of probe content, evaluating genomic distribution and changes from previous array versions. We characterise EPICv2's new feature of replicated probes and provide recommendations as to the superior probes. In silico analysis of probe sequences demonstrates that probe cross-hybridisation remains a significant problem in EPICv2. By mapping the off-target sites at single nucleotide resolution and comparing with WGBS we show empirical evidence for preferential off-target binding.

CONCLUSIONS

Overall, we find EPICv2 a worthy successor to the previous Infinium methylation microarrays, however some technical issues remain. To support optimal EPICv2 data analysis we provide an expanded version of the EPICv2 manifest to aid researchers in understanding probe design, data processing, choosing appropriate probes for analysis and for integration with methylation datasets from previous versions of the Infinium Methylation BeadChip.

Sex-specific genetic architecture of blood pressure

The genetic and genomic basis of sex differences in blood pressure (BP) traits remain unstudied at scale. Here, we conducted sex-stratified and combined-sex genome-wide association studies of BP traits using the UK Biobank resource, identifying 1,346 previously reported and 29 new BP trait-associated loci. Among associated loci, 412 were female-specific (Pfemale ≤ 5 × 10-8; Pmale > 5 × 10-8) and 142 were male-specific (Pmale ≤ 5 × 10-8; Pfemale > 5 × 10-8); these sex-specific loci were enriched for hormone-related transcription factors, in particular, estrogen receptor 1. Analyses of gene-by-sex interactions and sexually dimorphic effects identified four genomic regions, showing female-specific associations with diastolic BP or pulse pressure, including the chromosome 13q34-COL4A1/COL4A2 locus. Notably, female-specific pulse pressure-associated loci exhibited enriched acetylated histone H3 Lys27 modifications in arterial tissues and a female-specific association with fibromuscular dysplasia, a female-biased vascular disease; colocalization signals included Chr13q34: COL4A1/COL4A2, Chr9p21: CDKN2B-AS1 and Chr4q32.1: MAP9 regions. Sex-specific and sex-biased polygenic associations of BP traits were associated with multiple cardiovascular traits. These findings suggest potentially clinically significant and BP sex-specific pleiotropic effects on cardiovascular diseases.

The potential of epigenetic therapy to target the 3D epigenome in endocrine-resistant breast cancer

Three-dimensional (3D) epigenome remodeling is an important mechanism of gene deregulation in cancer. However, its potential as a target to counteract therapy resistance remains largely unaddressed. Here, we show that epigenetic therapy with decitabine (5-Aza-mC) suppresses tumor growth in xenograft models of pre-clinical metastatic estrogen receptor positive (ER+) breast tumor. Decitabine-induced genome-wide DNA hypomethylation results in large-scale 3D epigenome deregulation, including de-compaction of higher-order chromatin structure and loss of boundary insulation of topologically associated domains. Significant DNA hypomethylation associates with ectopic activation of ER-enhancers, gain in ER binding, creation of new 3D enhancer-promoter interactions and concordant up-regulation of ER-mediated transcription pathways. Importantly, long-term withdrawal of epigenetic therapy partially restores methylation at ER-enhancer elements, resulting in a loss of ectopic 3D enhancer-promoter interactions and associated gene repression. Our study illustrates the potential of epigenetic therapy to target ER+ endocrine-resistant breast cancer by DNA methylation-dependent rewiring of 3D chromatin interactions, which are associated with the suppression of tumor growth.

Topological regulation of the estrogen transcriptional response by ZATT-mediated inhibition of TOP2B activity

Human type-II topoisomerases, TOP2A and TOP2B, remove transcription associated DNA supercoiling, thereby affecting gene-expression programs, and have recently been associated with 3D genome architecture. Here, we study the regulatory roles of TOP2 paralogs in response to estrogen, which triggers an acute transcriptional induction that involves rewiring of genome organization. We find that, whereas TOP2A facilitates transcription, as expected for a topoisomerase, TOP2B limits the estrogen response. Consistent with this, TOP2B activity is locally downregulated upon estrogen treatment to favor the establishment and stabilization of regulatory chromatin contacts, likely through an accumulation of DNA supercoiling. We show that estrogen-mediated inhibition of TOP2B requires estrogen receptor α (ERα), a non-catalytic function of TOP2A, and the action of the atypical SUMO-ligase ZATT. This mechanism of topological transcriptional-control, which may be shared by additional gene-expression circuits, highlights the relevance of DNA topoisomerases as central actors of genome dynamics.

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

Spatial -omics technologies: the new enterprise in 3D breast cancer models

The fields of tissue bioengineering, -omics, and spatial biology are advancing rapidly, each offering the opportunity for a paradigm shift in breast cancer research. However, to date, collaboration between these fields has not reached its full potential. In this review, we describe the most recently generated 3D breast cancer models regarding the biomaterials and technological platforms employed. Additionally, their biological evaluation is reported, highlighting their advantages and limitations. Specifically, we focus on the most up-to-date -omics and spatial biology techniques, which can generate a deeper understanding of the biological relevance of bioengineered 3D breast cancer in vitro models, thus paving the way towards truly clinically relevant microphysiological systems, improved drug development success rates, and personalised medicine approaches.

Proteasome inhibition reprograms chromatin landscape in breast cancer

The 26S proteasome is the major protein degradation machinery in cells. Cancer cells use the proteasome to modulate gene expression networks that promote tumor growth. Proteasome inhibitors have emerged as effective cancer therapeutics, but how they work mechanistically remains unclear. Here, using integrative genomic analysis, we discovered unexpected reprogramming of the chromatin landscape and RNAPII transcription initiation in breast cancer cells treated with the proteasome inhibitor MG132. The cells acquired dynamic changes in chromatin accessibility at specific genomic loci termed Differentially Open Chromatin Regions (DOCRs). DOCRs with decreased accessibility were promoter proximal and exhibited unique chromatin architecture associated with divergent RNAPII transcription. Conversely, DOCRs with increased accessibility were primarily distal to transcription start sites and enriched in oncogenic super enhancers predominantly accessible in non-basal breast tumor subtypes. These findings describe the mechanisms by which the proteasome modulates the expression of gene networks intrinsic to breast cancer biology.

Highlights

Proteasome inhibition uncovers de novo Differential Open Chromatin Regions (DOCRs) in breast cancer cells. Proteasome inhibitor sensitive promoters exhibit a distinctive chromatin architecture with discrete transcription initiation patterns.Proteasome inhibition reprograms accessibility of super enhancers.Proteasome inhibitor sensitive super enhancers distinguish basal from non-basal breast cancer subtypes.

Osteosarcoma Multi-Omics Landscape and Subtypes

Osteosarcoma (OS) is the most common primary bone malignancy that exhibits remarkable histologic diversity and genetic heterogeneity. The complex nature of osteosarcoma has confounded precise molecular categorization, prognosis, and prediction for this disease. In this study, we performed a comprehensive multiplatform analysis on 86 osteosarcoma tumors, including somatic copy-number alteration, gene expression and methylation, and identified three molecularly distinct and clinically relevant subtypes of osteosarcoma. The subgrouping criteria was validated on another cohort of osteosarcoma tumors. Previously unappreciated osteosarcoma-type-specific changes in specific genes' copy number, expression and methylation were revealed based on the subgrouping. The subgrouping and novel gene signatures provide insights into refining osteosarcoma therapy and relationships to other types of cancer.

Combination of multiple omics techniques for a personalized therapy or treatment selection

Despite targeted therapies and immunotherapies have revolutionized the treatment of cancer patients, only a limited number of patients have long-term responses. Moreover, due to differences within cancer patients in the tumor mutational burden, composition of the tumor microenvironment as well as of the peripheral immune system and microbiome, and in the development of immune escape mechanisms, there is no "one fit all" therapy. Thus, the treatment of patients must be personalized based on the specific molecular, immunologic and/or metabolic landscape of their tumor. In order to identify for each patient the best possible therapy, different approaches should be employed and combined. These include (i) the use of predictive biomarkers identified on large cohorts of patients with the same tumor type and (ii) the evaluation of the individual tumor with "omics"-based analyses as well as its ex vivo characterization for susceptibility to different therapies.

Recapitulation of patient-specific 3D chromatin conformation using machine learning

Regulatory networks containing enhancer-gene edges define cellular states. Multiple efforts have revealed these networks for reference tissues and cell lines by integrating multi-omics data. However, the methods developed cannot be applied for large patient cohorts due to the infeasibility of chromatin immunoprecipitation sequencing (ChIP-seq) for limited biopsy material. We trained machine-learning models using chromatin interaction analysis with paired-end tag sequencing (ChIA-PET) and high-throughput chromosome conformation capture combined with chromatin immunoprecipitation (HiChIP) data that can predict connections using only assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA-seq data as input, which can be generated from biopsies. Our method overcomes limitations of correlation-based approaches that cannot distinguish between distinct target genes of given enhancers or between active vs. poised states in different samples, a hallmark of network rewiring in cancer. Application of our model on 371 samples across 22 cancer types revealed 1,780 enhancer-gene connections for 602 cancer genes. Using CRISPR interference (CRISPRi), we validated enhancers predicted to regulate ESR1 in estrogen receptor (ER)+ breast cancer and A1CF in liver hepatocellular carcinoma.

The EstroGene Database Reveals Diverse Temporal, Context-Dependent, and Bidirectional Estrogen Receptor Regulomes in Breast Cancer

As one of the most successful cancer therapeutic targets, estrogen receptor-α (ER/ESR1) has been extensively studied over the past few decades. Sequencing technological advances have enabled genome-wide analysis of ER action. However, comparison of individual studies is limited by different experimental designs, and few meta-analyses are available. Here, we established the EstroGene database through unified processing of data from 246 experiments including 136 transcriptomic, cistromic, and epigenetic datasets focusing on estradiol (E2)-triggered ER activation across 19 breast cancer cell lines. A user-friendly browser (https://estrogene.org/) was generated for multiomic data visualization involving gene inquiry under user-defined experimental conditions and statistical thresholds. Notably, annotation of metadata associated with public datasets revealed a considerable lack of experimental details. Comparison of independent RNA-seq or ER ChIP-seq data with the same design showed large variability and only strong effects could be consistently detected. Temporal estrogen response metasignatures were defined, and the association of E2 response rate with temporal transcriptional factors, chromatin accessibility, and heterogeneity of ER expression was evaluated. Unexpectedly, harmonizing 146 E2-induced transcriptomic datasets uncovered a subset of genes harboring bidirectional E2 regulation, which was linked to unique transcriptional factors and highly associated with immune surveillance in the clinical setting. Furthermore, the context dependent E2 response programs were characterized in MCF7 and T47D cell lines, the two most frequently used models in the EstroGene database. Collectively, the EstroGene database provides an informative and practical resource to the cancer research community to uniformly evaluate key reproducible features of ER regulomes and unravels modes of ER signaling.

SIGNIFICANCE

A resource database integrating 246 publicly available ER profiling datasets facilitates meta-analyses and identifies estrogen response temporal signatures, a bidirectional program, and model-specific biases.

Identification of metabolic pathways contributing to ER+ breast cancer disparities using a machine-learning pipeline

African American (AA) women in the United States have a 40% higher breast cancer mortality rate than Non-Hispanic White (NHW) women. The survival disparity is particularly striking among (estrogen receptor positive) ER+ breast cancer cases. The purpose of this study is to examine whether there are racial differences in metabolic pathways typically activated in patients with ER+ breast cancer. We collected pretreatment plasma from AA and NHW ER+ breast cancer cases (AA n = 48, NHW n = 54) and cancer-free controls (AA n = 100, NHW n = 48) to conduct an untargeted metabolomics analysis using gas chromatography mass spectrometry (GC-MS) to identify metabolites that may be altered in the different racial groups. Unpaired t-test combined with multiple feature selection and prediction models were employed to identify race-specific altered metabolic signatures. This was followed by the identification of altered metabolic pathways with a focus in AA patients with breast cancer. The clinical relevance of the identified pathways was further examined in PanCancer Atlas breast cancer data set from The Cancer Genome Atlas Program (TCGA). We identified differential metabolic signatures between NHW and AA patients. In AA patients, we observed decreased circulating levels of amino acids compared to healthy controls, while fatty acids were significantly higher in NHW patients. By mapping these metabolites to potential epigenetic regulatory mechanisms, this study identified significant associations with regulators of metabolism such as methionine adenosyltransferase 1A (MAT1A), DNA Methyltransferases and Histone methyltransferases for AA individuals, and Fatty acid Synthase (FASN) and Monoacylglycerol lipase (MGL) for NHW individuals. Specific gene Negative Elongation Factor Complex E (NELFE) with histone methyltransferase activity, was associated with poor survival exclusively for AA individuals. We employed a comprehensive and novel approach that integrates multiple machine learning and statistical methods, coupled with human functional pathway analyses. The metabolic profile of plasma samples identified may help elucidate underlying molecular drivers of disproportionately aggressive ER+ tumor biology in AA women. It may ultimately lead to the identification of novel therapeutic targets. To our knowledge, this is a novel finding that describes a link between metabolic alterations and epigenetic regulation in AA breast cancer and underscores the need for detailed investigations into the biological underpinnings of breast cancer health disparities.

SERPINA3-ANKRD11-HDAC3 pathway induced aromatase inhibitor resistance in breast cancer can be reversed by HDAC3 inhibition

Endocrine resistance is a major challenge for breast cancer therapy. To identify the genes pivotal for endocrine-resistance progression, we screened five datasets and found 7 commonly dysregulated genes in endocrine-resistant breast cancer cells. Here we show that downregulation of serine protease inhibitor clade A member 3 (SERPINA3) which is a direct target gene of estrogen receptor α contributes to aromatase inhibitor resistance. Ankyrin repeat domain containing 11 (ANKRD11) works as a downstream effector of SERPINA3 in mediating endocrine-resistance. It induces aromatase inhibitor insensitivity by interacting with histone deacetylase 3 (HDAC3) and upregulating its activity. Our study suggests that aromatase inhibitor therapy downregulates SERPINA3 and leads to the ensuing upregulation of ANKRD11, which in turn promotes aromatase inhibitor resistance via binding to and activating HDAC3. HDAC3 inhibition may reverse the aromatase inhibitor resistance in ER-positive breast cancer with decreased SERPINA3 and increased ANKRD11 expression.

Comprehensive multiomics analyses reveal pervasive involvement of aberrant cohesin binding in transcriptional and chromosomal disorder of cancer cells

Chromatin organization, whose malfunction causes various diseases including cancer, is fundamentally controlled by cohesin. While cancer cells have been found with mutated or misexpressed cohesin genes, there is no comprehensive survey about the presence and role of abnormal cohesin binding in cancer cells. Here, we systematically identified ∼1% of cohesin-binding sites (701-2,633) as cancer-aberrant binding sites of cohesin (CASs). We integrated CASs with large-scale transcriptomics, epigenomics, 3D genomics, and clinical information. CASs represent tissue-specific epigenomic signatures enriched for cancer-dysregulated genes with functional and clinical significance. CASs exhibited alterations in chromatin compartments, loops within topologically associated domains, and cis-regulatory elements, indicating that CASs induce dysregulated genes through misguided chromatin structure. Cohesin depletion data suggested that cohesin binding at CASs actively regulates cancer-dysregulated genes. Overall, our comprehensive investigation suggests that aberrant cohesin binding is an essential epigenomic signature responsible for dysregulated chromatin structure and transcription in cancer cells.

Stratification of Tamoxifen Synergistic Combinations for the Treatment of ER+ Breast Cancer

Breast cancer alone accounts for the majority of cancer deaths among women, with the most commonly diagnosed subtype being estrogen receptor positive (ER+). Survival has greatly improved for patients with ER+ breast cancer, due in part to the development of antiestrogen compounds, such as tamoxifen. While treatment of the primary disease is often successful, as many as 30% of patients will experience recurrence and metastasis, mainly due to developed endocrine therapy resistance. In this study, we discovered two tamoxifen combination therapies, with simeprevir and VX-680, that reduce the tumor burden in animal models of ER+ breast cancer more than either compound or tamoxifen alone. Additionally, these tamoxifen combinations reduced the expression of HER2, a hallmark of tamoxifen treatment, which can facilitate acquisition of a treatment-resistant phenotype. These combinations could provide clinical benefit by potentiating tamoxifen treatment in ER+ breast cancer.

LPAD: using network construction and label propagation to detect topologically associating domains from Hi-C data

With the development of chromosome conformation capture technique, the study of spatial conformation of a genome based on Hi-C technique has made a quantum leap. Previous studies reveal that genomes are folded into hierarchy of three-dimensional (3D) structures associated with topologically associating domains (TADs), and detecting TAD boundaries is of great significance in the chromosome-level analysis of 3D genome architecture. In this paper, we propose a novel TAD identification method, LPAD, which first extracts node correlations from global interactions of chromosomes based on the random walk with restart and then builds an undirected graph from Hi-C contact matrix. Next, LPAD designs a label propagation-based approach to discover communities and generates TADs. Experimental results verify the effectiveness and quality of TAD detections compared with existing methods. Furthermore, experimental evaluation of chromatin immunoprecipitation sequencing data shows that LPAD performs high enrichment of histone modifications remarkably nearby the TAD boundaries, and these results demonstrate LPAD's advantages on TAD identification accuracy.

BAP18 facilitates CTCF-mediated chromatin accessible to regulate enhancer activity in breast cancer

The estrogen receptor alpha (ERα) signaling pathway is a crucial target for ERα-positive breast cancer therapeutic strategies. Co-regulators and other transcription factors cooperate for effective ERα-related enhancer activation. Recent studies demonstrate that the transcription factor CTCF is essential to participate in ERα/E2-induced enhancer transactivation. However, the mechanism of how CTCF is achieved remains unknown. Here, we provided evidence that BAP18 is required for CTCF recruitment on ERα-enriched enhancers, facilitating CTCF-mediated chromatin accessibility to promote enhancer RNAs transcription. Consistently, GRO-seq demonstrates that the enhancer activity is positively correlated with BAP18 enrichment. Furthermore, BAP18 interacts with SMARCA1/BPTF to accelerate the recruitment of CTCF to ERα-related enhancers. Interestingly, BAP18 is involved in chromatin accessibility within enhancer regions, thereby increasing enhancer transactivation and enhancer-promoter looping. BAP18 depletion increases the sensitivity of anti-estrogen and anti-enhancer treatment in MCF7 cells. Collectively, our study indicates that BAP18 coordinates with CTCF to enlarge the transactivation of ERα-related enhancers, providing a better understanding of BAP18/CTCF coupling chromatin remodeling and E-P looping in the regulation of enhancer transcription.

Rewiring of the 3D genome during acquisition of carboplatin resistance in a triple-negative breast cancer patient-derived xenograft

Changes in the three-dimensional (3D) structure of the genome are an emerging hallmark of cancer. Cancer-associated copy number variants and single nucleotide polymorphisms promote rewiring of chromatin loops, disruption of topologically associating domains (TADs), active/inactive chromatin state switching, leading to oncogene expression and silencing of tumor suppressors. However, little is known about 3D changes during cancer progression to a chemotherapy-resistant state. We integrated chromatin conformation capture (Hi-C), RNA-seq, and whole-genome sequencing obtained from triple-negative breast cancer patient-derived xenograft primary tumors (UCD52) and carboplatin-resistant samples and found increased short-range (< 2 Mb) interactions, chromatin looping, formation of TAD, chromatin state switching into a more active state, and amplification of ATP-binding cassette transporters. Transcriptome changes suggested the role of long-noncoding RNAs in carboplatin resistance. Rewiring of the 3D genome was associated with TP53, TP63, BATF, FOS-JUN family of transcription factors and led to activation of aggressiveness-, metastasis- and other cancer-related pathways. Integrative analysis highlighted increased ribosome biogenesis and oxidative phosphorylation, suggesting the role of mitochondrial energy metabolism. Our results suggest that 3D genome remodeling may be a key mechanism underlying carboplatin resistance.

PHF6 recruits BPTF to promote HIF-dependent pathway and progression in YAP-high breast cancer

Background

Aberrant epigenetic remodeling events contribute to progression and metastasis of breast cancer (Bca). The specific mechanims that epigenetic factors rely on to mediate tumor aggressiveness remain unclear. We aimed to elucidate the roles of epigenetic protein PHF6 in breast tumorigenesis.

Methods

Published datasets and tissue samples with PHF6 staining were used to investigate the clinical relevance of PHF6 in Bca. CCK-8, clony formation assays were used to assess cell growth capacity. Cell migration and invasion abilities were measured by Transwell assay. The gene mRNA and protein levels were measured by quantitative real-time PCR and western blot. Chromatin immunoprecipitation (ChIP)-qPCR assays were used to investigate transcriptional relationships among genes. The Co-immunoprecipitation (Co-IP) assay was used to validate interactions between proteins. The CRISPR/Cas9 editing technology was used to construct double HIF knockout (HIF-DKO) cells. The subcutaneous xenograft model and orthotopic implantation tumor model were used to asess in vivo tumor growth.

Results

In this study, we utilized MTT assay to screen that PHF6 is required for Bca growth. PHF6 promotes Bca proliferation and migration. By analyzing The Cancer Genome Atlas breast cancer (TCGA-Bca) cohort, we found that PHF6 was significantly higher in tumor versus normal tissues. Mechanistically, PHF6 physically interacts with HIF-1α and HIF-2α to potentiate HIF-driven transcriptional events to initiate breast tumorigenesis. HIF-DKO abolished PHF6-mediated breast tumor growth, and PHF6 deficiency in turn impaired HIF transcriptional effects. Besides, hypoxia could also rely on YAP activation, but not HIF, to sustain PHF6 expressions in Bca cells. In addition, PHF6 recuits BPTF to mediate epigenetic remodeling to augment HIF transcriptional activity. Targeting PHF6 or BPTF inhibitor (AU1) is effective in mice models. Lastly, PHF6 correlated with HIF target gene expression in human breast tumors, which is an independent prognostic regulator.

Conclusions

Collectively, this study identified PHF6 as a prognostic epigenetic regulator for Bca, functioning as a HIF coactivator. The fundamental mechanisms underlying YAP/PHF6/HIF axis in breast tumors endowed novel epigenegtic targets for Bca treatment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-023-04031-8.

Bypassing anti-PD-(L)1 therapy: Mechanisms and management strategies

Resistance to immune checkpoint inhibitors (ICIs) is a major issue of the current era in cancer immunotherapy. Immune evasion is a multi-factorial event, which occurs generally at a base of cold immunity. Despite advances in the field, there are still unsolved challenges about how to combat checkpoint hijacked by tumor cells and what are complementary treatment strategies to render durable anti-tumor outcomes. A point is that anti-programed death-1 receptor (PD-1)/anti-programmed death-ligand 1 (PD-L1) is not the solo path of immune escape, and responses in many types of solid tumors to the PD-1/PD-L1 inhibitors are not satisfactory. Thus, seeking mechanisms inter-connecting tumor with its immune ecosystem nearby unravel more about resistance mechanisms so as to develop methods for sustained reinvigoration of immune activity against cancer. In this review, we aimed to discuss about common and specific paths taken by tumor cells to evade immune surveillance, describing novel detection strategies, as well as suggesting some approaches to recover tumor sensitivity to the anti-PD-(L)1 therapy based on the current knowledge.

Epigenetically programmed resistance to chemo- and immuno-therapies

Resistance to cancer treatments remains a major barrier in developing cancer cures. While promising combination chemotherapy treatments and novel immunotherapies have improved patient outcomes, resistance to these treatments remains poorly understood. New insights into the dysregulation of the epigenome show how it promotes tumor growth and resistance to therapy. By altering control of gene expression, tumor cells can evade immune cell recognition, ignore apoptotic cues, and reverse DNA damage induced by chemotherapies. In this chapter, we summarize the data on epigenetic remodeling during cancer progression and treatment that enable cancer cell survival and describe how these epigenetic changes are being targeted clinically to overcome resistance.

ESR1 gene amplification and MAP3K mutations are selected during adjuvant endocrine therapies in relapsing Hormone Receptor-positive, HER2-negative breast cancer (HR+ HER2- BC)

BACKGROUND

Previous studies have provided a comprehensive picture of genomic alterations in primary and metastatic Hormone Receptor (HR)-positive, Human Epidermal growth factor Receptor 2 (HER2)-negative breast cancer (HR+ HER2- BC). However, the evolution of the genomic landscape of HR+ HER2- BC during adjuvant endocrine therapies (ETs) remains poorly investigated.

METHODS AND FINDINGS

We performed a genomic characterization of surgically resected HR+ HER2- BC patients relapsing during or at the completion of adjuvant ET. Using a customized panel, we comprehensively evaluated gene mutations and copy number variation (CNV) in paired primary and metastatic specimens. After retrieval and quality/quantity check of tumor specimens from an original cohort of 204 cases, 74 matched tumor samples were successfully evaluated for DNA mutations and CNV analysis. Along with previously reported genomic alterations, including PIK3CA, TP53, CDH1, GATA3 and ESR1 mutations/deletions, we found that ESR1 gene amplification (confirmed by FISH) and MAP3K mutations were enriched in metastatic lesions as compared to matched primary tumors. These alterations were exclusively found in patients treated with adjuvant aromatase inhibitors or LHRH analogs plus tamoxifen, but not in patients treated with tamoxifen alone. Patients with tumors bearing MAP3K mutations in metastatic lesions had significantly worse distant relapse-free survival (hazard ratio [HR] 3.4, 95% CI 1.52-7.70, p value 0.003) and worse overall survival (HR 5.2, 95% CI 2.10-12.8, p-value < 0.001) independently of other clinically relevant patient- and tumor-related variables.

CONCLUSIONS

ESR1 amplification and activating MAP3K mutations are potential drivers of acquired resistance to adjuvant ETs employing estrogen deprivation in HR+ HER2- BC. MAP3K mutations are associated with worse prognosis in patients with metastatic disease.

Multiomics in primary and metastatic breast tumors from the AURORA US network finds microenvironment and epigenetic drivers of metastasis

The AURORA US Metastasis Project was established with the goal to identify molecular features associated with metastasis. We assayed 55 females with metastatic breast cancer (51 primary cancers and 102 metastases) by RNA sequencing, tumor/germline DNA exome and low-pass whole-genome sequencing and global DNA methylation microarrays. Expression subtype changes were observed in ~30% of samples and were coincident with DNA clonality shifts, especially involving HER2. Downregulation of estrogen receptor (ER)-mediated cell-cell adhesion genes through DNA methylation mechanisms was observed in metastases. Microenvironment differences varied according to tumor subtype; the ER+/luminal subtype had lower fibroblast and endothelial content, while triple-negative breast cancer/basal metastases showed a decrease in B and T cells. In 17% of metastases, DNA hypermethylation and/or focal deletions were identified near HLA-A and were associated with reduced expression and lower immune cell infiltrates, especially in brain and liver metastases. These findings could have implications for treating individuals with metastatic breast cancer with immune- and HER2-targeting therapies.

Bioinformatic Evaluation of Features on Cis-regulatory Elements at 6q25.1

Eukaryotic non-coding regulatory features contribute significantly to cellular plasticity which on aberration leads to cellular malignancy. Enhancers are cis-regulatory elements that contribute to the development of resistance to endocrine therapy in estrogen receptor (ER)-positive breast cancer leading to poor clinical outcome. ER is vital for therapeutic targets in ER-positive breast cancer. Here, we review and report the different regulatory features present on ER with the objective to delineate potential mechanisms which may contribute to development of resistance. The UCSC Genome Browser, data mining, and bioinformatics tools were used to review enhancers, transcription factors (TFs), histone marks, long non-coding RNAs (lncRNAs), and variants residing in the non-coding region of the ER gene. We report 7 enhancers, 3 of which were rich in TF-binding sites and histone marks in a cell line-specific manner. Furthermore, some enhancers contain estrogen resistance variants and sites for lncRNA. Our review speculates putative models suggesting potential aberrations in gene regulation and expression if these regulatory landscapes and assemblies are altered. This review gives an interesting perspective in designing integrated in vitro studies including non-coding elements to study development of endocrine resistance in ER-positive breast cancer.

AKTIP loss is enriched in ERα-positive breast cancer for tumorigenesis and confers endocrine resistance

Recurrent deletion of 16q12.2 is observed in luminal breast cancer, yet the causal genomic alterations in this region are largely unknown. In this study, we identify that loss of AKTIP, which is located on 16q12.2, drives tumorigenesis of estrogen receptor alpha (ERα)-positive, but not ERα-negative, breast cancer cells and is associated with poor prognosis of patients with ERα-positive breast cancer. Intriguingly, AKTIP-depleted tumors have increased ERα protein level and activity. Cullin-associated and neddylation-dissociated protein 1 (CAND1), which regulates the cullin-RING E3 ubiquitin ligases, protects ERα from cullin 2-dependent proteasomal degradation. Apart from ERα signaling, AKTIP loss triggers JAK2-STAT3 activation, which provides an alternative survival signal when ERα is inhibited. AKTIP-depleted MCF7 cells and ERα-positive patient-derived organoids are more resistant to ERα antagonists. Importantly, the resistance can be overcome by co-inhibition of JAK2/STAT3. Together, our results highlight the subtype-specific functional consequences of AKTIP loss and provide a mechanistic explanation for the enriched AKTIP copy-number loss in ERα-positive breast cancer.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

Subtype-specific 3D genome alteration in acute myeloid leukaemia

Acute myeloid leukaemia (AML) represents a set of heterogeneous myeloid malignancies, and hallmarks include mutations in epigenetic modifiers, transcription factors and kinases1-5. The extent to which mutations in AML drive alterations in chromatin 3D structure and contribute to myeloid transformation is unclear. Here we use Hi-C and whole-genome sequencing to analyse 25 samples from patients with AML and 7 samples from healthy donors. Recurrent and subtype-specific alterations in A/B compartments, topologically associating domains and chromatin loops were identified. RNA sequencing, ATAC with sequencing and CUT&Tag for CTCF, H3K27ac and H3K27me3 in the same AML samples also revealed extensive and recurrent AML-specific promoter-enhancer and promoter-silencer loops. We validated the role of repressive loops on their target genes by CRISPR deletion and interference. Structural variation-induced enhancer-hijacking and silencer-hijacking events were further identified in AML samples. Hijacked enhancers play a part in AML cell growth, as demonstrated by CRISPR screening, whereas hijacked silencers have a downregulating role, as evidenced by CRISPR-interference-mediated de-repression. Finally, whole-genome bisulfite sequencing of 20 AML and normal samples revealed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Treatment of AML cells with a DNA hypomethylating agent and triple knockdown of DNMT1, DNMT3A and DNMT3B enabled the manipulation of DNA methylation to revert 3D genome organization and gene expression. Overall, this study provides a resource for leukaemia studies and highlights the role of repressive loops and hijacked cis elements in human diseases.

Constricted migration is associated with stable 3D genome structure differences in cancer cells

To spread from a localized tumor, metastatic cancer cells must squeeze through constrictions that cause major nuclear deformations. Since chromosome structure affects nucleus stiffness, gene regulation, and DNA repair, here, we investigate the relationship between 3D genome structure and constricted migration in cancer cells. Using melanoma (A375) cells, we identify phenotypic differences in cells that have undergone multiple rounds of constricted migration. These cells display a stably higher migration efficiency, elongated morphology, and differences in the distribution of Lamin A/C and heterochromatin. Hi-C experiments reveal differences in chromosome spatial compartmentalization specific to cells that have passed through constrictions and related alterations in expression of genes associated with migration and metastasis. Certain features of the 3D genome structure changes, such as a loss of B compartment interaction strength, are consistently observed after constricted migration in clonal populations of A375 cells and in MDA-MB-231 breast cancer cells. Our observations suggest that consistent types of chromosome structure changes are induced or selected by passage through constrictions and that these may epigenetically encode stable differences in gene expression and cellular migration phenotype.

EWS/FLI mediated reprogramming of 3D chromatin promotes an altered transcriptional state in Ewing sarcoma

Ewing sarcoma is a prototypical fusion transcription factor-associated pediatric cancer that expresses EWS/FLI or a highly related FET/ETS chimera. EWS/FLI dysregulates transcription to induce and maintain sarcomagenesis, but the mechanisms utilized are not fully understood. We therefore sought to define the global effects of EWS/FLI on chromatin conformation and transcription in Ewing sarcoma cells using a well-validated ‘knock-down/rescue’ model of EWS/FLI function in combination with next generation sequencing assays to evaluate how the chromatin landscape changes with loss, and recovery, of EWS/FLI expression. We found that EWS/FLI (and EWS/ERG) genomic localization is largely conserved across multiple patient-derived Ewing sarcoma cell lines. This EWS/FLI binding signature is associated with establishment of topologically-associated domain (TAD) boundaries, compartment activation, enhancer-promoter looping that involve both intra- and inter-TAD interactions, and gene activation. In addition, EWS/FLI co-localizes with the loop-extrusion factor cohesin to promote chromatin loops and TAD boundaries. Importantly, local chromatin features provide the basis for transcriptional heterogeneity in regulation of direct EWS/FLI target genes across different Ewing sarcoma cell lines. These data demonstrate a key role of EWS/FLI in mediating genome-wide changes in chromatin configuration and support the notion that fusion transcription factors serve as master regulators of three-dimensional reprogramming of chromatin.

Super-enhancer-controlled positive feedback loop BRD4/ERα-RET-ERα promotes ERα-positive breast cancer

Estrogen and estrogen receptor alpha (ERα)-induced gene transcription is tightly associated with ERα-positive breast carcinogenesis. ERα-occupied enhancers, particularly super-enhancers, have been suggested to play a vital role in regulating such transcriptional events. However, the landscape of ERα-occupied super-enhancers (ERSEs) as well as key ERα-induced target genes associated with ERSEs remain to be fully characterized. Here, we defined the landscape of ERSEs in ERα-positive breast cancer cell lines, and demonstrated that bromodomain protein BRD4 is a master regulator of the transcriptional activation of ERSEs and cognate ERα target genes. RET, a member of the tyrosine protein kinase family of proteins, was identified to be a key ERα target gene of BRD4-regulated ERSEs, which, in turn, is vital for ERα-induced gene transcriptional activation and malignant phenotypes through activating the RAS/RAF/MEK2/ERK/p90RSK/ERα phosphorylation cascade. Combination therapy with BRD4 and RET inhibitors exhibited additive effects on suppressing ERα-positive breast cancer both in vitro and in vivo, comparable with that of standard endocrine therapy tamoxifen. Furthermore, combination therapy re-sensitized a tamoxifen-resistant ERα-positive breast cancer cell line to tamoxifen treatment. Taken together, our data uncovered the critical role of a super-enhancer-associated positive feedback loop constituting BRD4/ERα–RET–ERα in ERα-positive breast cancer, and suggested that targeting components in this loop would provide a new therapeutic avenue for treating ERα-positive breast cancer in the clinic.

Reorganization of 3D chromatin architecture in doxorubicin-resistant breast cancer cells

Background: Doxorubicin resistance remains a major therapeutic challenge leading to poor survival prognosis and treatment failure in breast cancer. Although doxorubicin induces massive changes in the transcriptional landscape are well known, potential diagnostic or therapeutic targets associated with the reorganization of three-dimensional (3D) chromatin architecture have not yet been systematically investigated.

Methods: Here we performed in situ high-throughput chromosome conformation capture (Hi-C) on parental and doxorubicin-resistant MCF7 (MCF7-DR) human breast cancer cells, followed by integrative analysis of HiC, ATAC-seq, RNA-seq and TCGA data.

Results: It revealed that A/B compartment switching was positively correlated to genome-wide differential gene expression. The genome of MCF7-DR cells was spatially reorganized into smaller topologically associating domains (TADs) and chromatin loops. We also revealed the contribution of increased chromatin accessibility and potential transcription factor families, including CTCF, AP-1 and bHLH, to gained TADs or loops. Intriguingly, we observed two condensed genomic regions (∼20 kb) with decreased chromatin accessibility flanking TAD boundaries, which might play a critical role in the formation or maintenance of TADs. Finally, combining data from TCGA, we identified a number of gained and lost enhancer-promoter interactions and their corresponding differentially expressed genes involved in chromatin organization and breast cancer signaling pathways, including FA2H, FOXA1 and JRKL, which might serve as potential treatment targets for breast cancer.

Conclusion: These data uncovered a close connection between 3D genome reorganization, chromatin accessibility as well as gene transcription and provide novel insights into the epigenomic mechanisms involving doxorubicin resistance in breast cancer.

ACBD3 Bioinformatic Analysis and Protein Expression in Breast Cancer Cells

ACBD3 overexpression has previously been found to correlate with worse prognosis for breast cancer patients and, as an incredibly diverse protein in both function and cellular localisation, ACBD3 may have a larger role in breast cancer than previously thought. This study further investigated ACBD3′s role in breast cancer. Bioinformatic databases were queried to characterise ACBD3 expression and mutation in breast cancer and to investigate how overexpression affects breast cancer patient outcomes. Immunohistochemistry was carried out to examine ACBD3 location within cells and tissue structures. ACBD3 was more highly expressed in breast cancer than in any other cancer or matched normal tissue, and expression over the median level resulted in reduced relapse-free, overall, and distant metastasis-free survival for breast cancer patients as a whole, with some differences observed between subtypes. IHC analysis found that ACBD3 levels varied based on hormone receptor status, indicating that ACBD3 could be a candidate biomarker for poor patient prognosis in breast cancer and may possibly be a biomarker for ER signal reprogramming of precancerous breast tissue.

CTCF and Its Partners: Shaper of 3D Genome during Development

The 3D genome organization and its dynamic modulate genome function, playing a pivotal role in cell differentiation and development. CTCF and cohesin, acting as the core architectural components involved in chromatin looping and genome folding, can also recruit other protein or RNA partners to fine-tune genome structure during development. Moreover, systematic screening for partners of CTCF has been performed through high-throughput approaches. In particular, several novel protein and RNA partners, such as BHLHE40, WIZ, MAZ, Aire, MyoD, YY1, ZNF143, and Jpx, have been identified, and these partners are mostly implicated in transcriptional regulation and chromatin remodeling, offering a unique opportunity for dissecting their roles in higher-order chromatin organization by collaborating with CTCF and cohesin. Here, we review the latest advancements with an emphasis on features of CTCF partners and also discuss the specific functions of CTCF-associated complexes in chromatin structure modulation, which may extend our understanding of the functions of higher-order chromatin architecture in developmental processes.

Genetic variation as a long-distance modulator of RAD21 expression in humans

Somatic mutations and changes in expression of RAD21 are common in many types of cancer. Moreover, sub-optimal levels of RAD21 expression in early development can result in cohesinopathies. Altered RAD21 levels can result directly from mutations in the RAD21 gene. However, whether DNA variants outside of the RAD21 gene could control its expression and thereby contribute to cancer and developmental disease is unknown. In this study, we searched for genomic variants that modify RAD21expression to determine their potential to contribute to development or cancer by RAD21 dysregulation. We searched 42,953,834 genomic variants for a spatial-eQTL association with the transcription of RAD21. We identified 123 significant associations (FDR < 0.05), which are local (cis) or long-distance (trans) regulators of RAD21 expression. The 123 variants co-regulate a further seven genes (AARD, AKAP11, GRID1, KCNIP4, RCN1, TRIOBP, and USP32), enriched for having Sp2 transcription factor binding sites in their promoter regions. The Sp2 transcription factor and six of the seven genes had previously been associated with cancer onset, progression, and metastasis. Our results suggest that genome-wide variation in non-coding regions impacts on RAD21 transcript levels in addition to other genes, which then could impact on oncogenesis and the process of ubiquitination. This identification of distant co-regulation of oncogenes represents a strategy for discovery of novel genetic regions influencing cancer onset and a potential for diagnostics.

Assessing comparative importance of DNA sequence and epigenetic modifications on gene expression using a deep convolutional neural network

Gene expression is regulated at both transcriptional and post-transcriptional levels. DNA sequence and epigenetic modifications are key factors which regulate gene transcription. Understanding their complex interactions and their respective contributions to gene expression regulation remains a challenge in biological studies. We have developed iSEGnet, a framework of deep convolutional neural network to predict mRNA abundance using the information on DNA sequences as well as epigenetic modifications within genes and their cis-regulatory regions. We demonstrate that our framework outperforms other machine learning models in terms of predicting mRNA abundance using transcriptional and epigenetic profiles from six distinct cell lines/types chosen from the ENCODE. The analysis from the learned models also reveals that specific regions around promotors and transcription termination sites are most important for gene expression regulation. Using the method of Integrated Gradients, we identify narrow segments in these regions which are most likely to impact gene expression for a specific epigenetic modification. We further show that these identified segments are enriched in known active regulatory regions by comparing the transcription factor binding sites obtained via ChIP-seq. Moreover, we demonstrate how iSEGnet can uncover potential transcription factors that have regulatory functions in cancer using two cancer multi-omics data.