A CTCF-independent role for cohesin in tissue-specific transcription

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.

Cohesin composition and dosage independently affect early development in zebrafish

Cohesin, a chromatin-associated protein complex with four core subunits (Smc1a, Smc3, Rad21 and either Stag1 or 2), has a central role in cell proliferation and gene expression in metazoans. Human developmental disorders termed 'cohesinopathies' are characterized by germline variants of cohesin or its regulators that do not entirely eliminate cohesin function. However, it is not clear whether mutations in individual cohesin subunits have independent developmental consequences. Here, we show that zebrafish rad21 or stag2b mutants independently influence embryonic tailbud development. Both mutants have altered mesoderm induction, but only homozygous or heterozygous rad21 mutation affects cell cycle gene expression. stag2b mutants have narrower notochords and reduced Wnt signaling in neuromesodermal progenitors as revealed by single-cell RNA sequencing. Stimulation of Wnt signaling rescues transcription and morphology in stag2b, but not rad21, mutants. Our results suggest that mutations altering the quantity versus composition of cohesin have independent developmental consequences, with implications for the understanding and management of cohesinopathies.

Transcription decouples estrogen-dependent changes in enhancer-promoter contact frequencies and spatial proximity

How enhancers regulate their target genes in the context of 3D chromatin organization is extensively studied and models which do not require direct enhancer-promoter contact have recently emerged. Here, we use the activation of estrogen receptor-dependent enhancers in a breast cancer cell line to study enhancer-promoter communication at two loci. This allows high temporal resolution tracking of molecular events from hormone stimulation to efficient gene activation. We examine how both enhancer-promoter spatial proximity assayed by DNA fluorescence in situ hybridization, and contact frequencies resulting from chromatin in situ fragmentation and proximity ligation, change dynamically during enhancer-driven gene activation. These orthogonal methods produce seemingly paradoxical results: upon enhancer activation enhancer-promoter contact frequencies increase while spatial proximity decreases. We explore this apparent discrepancy using different estrogen receptor ligands and transcription inhibitors. Our data demonstrate that enhancer-promoter contact frequencies are transcription independent whereas altered enhancer-promoter proximity depends on transcription. Our results emphasize that the relationship between contact frequencies and physical distance in the nucleus, especially over short genomic distances, is not always a simple one.

Churros: a Docker-based pipeline for large-scale epigenomic analysis

The epigenome, which reflects the modifications on chromatin or DNA sequences, provides crucial insight into gene expression regulation and cellular activity. With the continuous accumulation of epigenomic datasets such as chromatin immunoprecipitation followed by sequencing (ChIP-seq) data, there is a great demand for a streamlined pipeline to consistently process them, especially for large-dataset comparisons involving hundreds of samples. Here, we present Churros, an end-to-end epigenomic analysis pipeline that is environmentally independent and optimized for handling large-scale data. We successfully demonstrated the effectiveness of Churros by analyzing large-scale ChIP-seq datasets with the hg38 or Telomere-to-Telomere (T2T) human reference genome. We found that applying T2T to the typical analysis workflow has important impacts on read mapping, quality checks, and peak calling. We also introduced a useful feature to study context-specific epigenomic landscapes. Churros will contribute a comprehensive and unified resource for analyzing large-scale epigenomic data.

Peak-agnostic high-resolution cis-regulatory circuitry mapping using single cell multiome data

Single same cell RNAseq/ATACseq multiome data provide unparalleled potential to develop high resolution maps of the cell-type specific transcriptional regulatory circuitry underlying gene expression. We present CREMA, a framework that recovers the full cis-regulatory circuitry by modeling gene expression and chromatin activity in individual cells without peak-calling or cell type labeling constraints. We demonstrate that CREMA overcomes the limitations of existing methods that fail to identify about half of functional regulatory elements which are outside the called chromatin 'peaks'. These circuit sites outside called peaks are shown to be important cell type specific functional regulatory loci, sufficient to distinguish individual cell types. Analysis of mouse pituitary data identifies a Gata2-circuit for the gonadotrope-enriched disease-associated Pcsk1 gene, which is experimentally validated by reduced gonadotrope expression in a gonadotrope conditional Gata2-knockout model. We present a web accessible human immune cell regulatory circuit resource, and provide CREMA as an R package.

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.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

PAXIP1 and STAG2 converge to maintain 3D genome architecture and facilitate promoter/enhancer contacts to enable stress hormone-dependent transcription

How steroid hormone receptors (SHRs) regulate transcriptional activity remains partly understood. Upon activation, SHRs bind the genome together with a co-regulator repertoire, crucial to induce gene expression. However, it remains unknown which components of the SHR-recruited co-regulator complex are essential to drive transcription following hormonal stimuli. Through a FACS-based genome-wide CRISPR screen, we functionally dissected the Glucocorticoid Receptor (GR) complex. We describe a functional cross-talk between PAXIP1 and the cohesin subunit STAG2, critical for regulation of gene expression by GR. Without altering the GR cistrome, PAXIP1 and STAG2 depletion alter the GR transcriptome, by impairing the recruitment of 3D-genome organization proteins to the GR complex. Importantly, we demonstrate that PAXIP1 is required for stability of cohesin on chromatin, its localization to GR-occupied sites, and maintenance of enhancer-promoter interactions. In lung cancer, where GR acts as tumor suppressor, PAXIP1/STAG2 loss enhances GR-mediated tumor suppressor activity by modifying local chromatin interactions. All together, we introduce PAXIP1 and STAG2 as novel co-regulators of GR, required to maintain 3D-genome architecture and drive the GR transcriptional programme following hormonal stimuli.

CRISPR screens in sister chromatid cohesion defective cells reveal PAXIP1-PAGR1 as regulator of chromatin association of cohesin

The cohesin complex regulates higher order chromosome architecture through maintaining sister chromatid cohesion and folding chromatin by DNA loop extrusion. Impaired cohesin function underlies a heterogeneous group of genetic syndromes and is associated with cancer. Here, we mapped the genetic dependencies of human cell lines defective of cohesion regulators DDX11 and ESCO2. The obtained synthetic lethality networks are strongly enriched for genes involved in DNA replication and mitosis and support the existence of parallel sister chromatid cohesion establishment pathways. Among the hits, we identify the chromatin binding, BRCT-domain containing protein PAXIP1 as a novel cohesin regulator. Depletion of PAXIP1 severely aggravates cohesion defects in ESCO2 mutant cells, leading to mitotic cell death. PAXIP1 promotes global chromatin association of cohesin, independent of DNA replication, a function that cannot be explained by indirect effects of PAXIP1 on transcription or DNA repair. Cohesin regulation by PAXIP1 requires its binding partner PAGR1 and a conserved FDF motif in PAGR1. PAXIP1 co-localizes with cohesin on multiple genomic loci, including active gene promoters and enhancers. Possibly, this newly identified role of PAXIP1-PAGR1 in regulating cohesin occupancy on chromatin is also relevant for previously described functions of PAXIP1 in transcription, immune cell maturation and DNA repair.

Context-dependent perturbations in chromatin folding and the transcriptome by cohesin and related factors

Cohesin regulates gene expression through context-specific chromatin folding mechanisms such as enhancer-promoter looping and topologically associating domain (TAD) formation by cooperating with factors such as cohesin loaders and the insulation factor CTCF. We developed a computational workflow to explore how three-dimensional (3D) structure and gene expression are regulated collectively or individually by cohesin and related factors. The main component is CustardPy, by which multi-omics datasets are compared systematically. To validate our methodology, we generated 3D genome, transcriptome, and epigenome data before and after depletion of cohesin and related factors and compared the effects of depletion. We observed diverse effects on the 3D genome and transcriptome, and gene expression changes were correlated with the splitting of TADs caused by cohesin loss. We also observed variations in long-range interactions across TADs, which correlated with their epigenomic states. These computational tools and datasets will be valuable for 3D genome and epigenome studies.

Oestrogen treatment restores dentate gyrus development in premature newborns by IGF1 regulation

Prematurely-born infants cared for in the neonatal units suffer from memory and learning deficits. Prematurity diminishes neurogenesis and synaptogenesis in the hippocampal dentate gyrus (DG). This dysmaturation of neurons is attributed to elevated PSD95, NMDR2A, and IGF1 levels. Since oestrogen treatment plays key roles in the development and plasticity of DG, we hypothesized that 17β-estradiol (E2) treatment would ameliorate neurogenesis and synaptogenesis in the DG, reversing cognitive deficits in premature newborns. Additionally, E2-induced recovery would be mediated by IGF1 signalling. These hypotheses were tested in a rabbit model of prematurity and nonmaternal care, in which premature kits were gavage-fed and reared by laboratory personnel. We compared E2- and vehicle-treated preterm kits for morphological, molecular, and behavioural parameters. We also treated kits with oestrogen degrader, RAD1901, and assessed IGF1 signalling. We found that E2 treatment increased the number of Tbr2+ and DCX+ neuronal progenitors and increased the density of glutamatergic synapses in the DG. E2 treatment restored PSD95 and NMDAR2A levels and cognitive function in preterm kits. Transcriptomic analyses showed that E2 treatment contributed to recovery by influencing interactions between IGF1R and neurodegenerative, as well as glutamatergic genes. ERα expression was reduced on completion of E2 treatment at D7, followed by D30 elevation. E2-induced fluctuation in ERα levels was associated with a reciprocal elevation in IGF1/2 expression at D7 and reduction at D30. ERα degradation by RAD1901 treatment enhanced IGF1 levels, suggesting ERα inhibits IGF1 expression. E2 treatment alleviates the prematurity-induced maldevelopment of DG and cognitive dysfunctions by regulating ERα and IGF1 levels.

Identification of Two Subsets of Subcompartment A1 Associated with High Transcriptional Activity and Frequent Loop Extrusion

Three-dimensional genome organization has been increasingly recognized as an important determinant of the precise regulation of gene expression in mammalian cells, yet the relationship between gene transcriptional activity and spatial subcompartment positioning is still not fully comprehended. Here, we first utilized genome-wide Hi-C data to infer eight types of subcompartment (labeled A1, A2, A3, A4, B1, B2, B3, and B4) in mouse embryonic stem cells and four primary differentiated cell types, including thymocytes, macrophages, neural progenitor cells, and cortical neurons. Transitions of subcompartments may confer gene expression changes in different cell types. Intriguingly, we identified two subsets of subcompartments defined by higher gene density and characterized by strongly looped contact domains, named common A1 and variable A1, respectively. We revealed that common A1, which includes highly expressed genes and abundant housekeeping genes, shows a ~2-fold higher gene density than the variable A1, where cell type-specific genes are significantly enriched. Thus, our study supports a model in which both types of genomic loci with constitutive and regulatory high transcriptional activity can drive the subcompartment A1 formation. Special chromatin subcompartment arrangement and intradomain interactions may, in turn, contribute to maintaining proper levels of gene expression, especially for regulatory non-housekeeping genes.

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.

The role of loop extrusion in enhancer-mediated gene activation

Gene expression patterns in complex multicellular organisms are regulated by enhancers, which communicate with their target gene promoters in three-dimensional (3D) chromatin structures. Despite advances in our understanding of the mechanisms that organize mammalian genomes into compartments and topologically associating domains (TADs), it is not well understood how specific interactions between enhancers and promoters are controlled in this 3D context. In this review, we give an overview of recent evidence that shows that a process of loop extrusion plays an important role in the regulation of enhancer-promoter communication and discuss recent insights into the molecular mechanism by which loop extrusion contributes to enhancer-mediated gene activation.

Full circle: a brief history of cohesin and the regulation of gene expression

The cohesin complex has a range of crucial functions in the cell. Cohesin is essential for mediating chromatid cohesion during mitosis, for repair of double-strand DNA breaks, and for control of gene transcription. This last function has been the subject of intense research ever since the discovery of cohesin's role in the long-range regulation of the cut gene in Drosophila. Subsequent research showed that the expression of some genes is exquisitely sensitive to cohesin depletion, while others remain relatively unperturbed. Sensitivity to cohesin depletion is also remarkably cell type- and/or condition-specific. The relatively recent discovery that cohesin is integral to forming chromatin loops via loop extrusion should explain much of cohesin's gene regulatory properties, but surprisingly, loop extrusion has failed to identify a 'one size fits all' mechanism for how cohesin controls gene expression. This review will illustrate how early examples of cohesin-dependent gene expression integrate with later work on cohesin's role in genome organization to explain mechanisms by which cohesin regulates gene expression.

Role of ZNF143 and Its Association with Gene Expression Patterns, Noncoding Mutations, and the Immune System in Human Breast Cancer

The function of noncoding sequence variations at ZNF143 binding sites in breast cancer cells is currently not well understood. Distal elements and promoters, also known as cis-regulatory elements, control the expression of genes. They may be identified by functional genomic techniques and sequence conservation, and they frequently show cell- and tissue-type specificity. The creation, destruction, or modulation of TF binding and function may be influenced by genetic modifications at TF binding sites that affect the binding affinity. Therefore, noncoding mutations that affect the ZNF143 binding site may be able to alter the expression of some genes in breast cancer. In order to understand the relationship among ZNF143, gene expression patterns, and noncoding mutations, we adopted an integrative strategy in this study and paid close attention to putative immunological signaling pathways. The immune system-related pathways ErbB, HIF1a, NF-kB, FoxO, JAK-STAT, Wnt, Notch, cell cycle, PI3K-AKT, RAP1, calcium signaling, cell junctions and adhesion, actin cytoskeleton regulation, and cancer pathways are among those that may be significant, according to the overall analysis.

CohesinDB: a comprehensive database for decoding cohesin-related epigenomes, 3D genomes and transcriptomes in human cells

Cohesin is a multifunctional protein responsible for transcriptional regulation and chromatin organization. Cohesin binds to chromatin at tens of thousands of distinct sites in a conserved or tissue-specific manner, whereas the function of cohesin varies greatly depending on the epigenetic properties of specific chromatin loci. Cohesin also extensively mediates cis-regulatory modules (CRMs) and chromatin loops. Even though next-generation sequencing technologies have provided a wealth of information on different aspects of cohesin, the integration and exploration of the resultant massive cohesin datasets are not straightforward. Here, we present CohesinDB (https://cohesindb.iqb.u-tokyo.ac.jp), a comprehensive multiomics cohesin database in human cells. CohesinDB includes 2043 epigenomics, transcriptomics and 3D genomics datasets from 530 studies involving 176 cell types. By integrating these large-scale data, CohesinDB summarizes three types of 'cohesin objects': 751 590 cohesin binding sites, 957 868 cohesin-related chromatin loops and 2 229 500 cohesin-related CRMs. Each cohesin object is annotated with locus, cell type, classification, function, 3D genomics and cis-regulatory information. CohesinDB features a user-friendly interface for browsing, searching, analyzing, visualizing and downloading the desired information. CohesinDB contributes a valuable resource for all researchers studying cohesin, epigenomics, transcriptional regulation and chromatin organization.

A cohesin-associated gene score may predict immune checkpoint blockade in hepatocellular carcinoma

Stromal antigen 1 (STAG1), a component of cohesion, is overexpressed in various cancers, but it is unclear whether it has a role in the transcriptional regulation of hepatocellular carcinoma (HCC). To test this hypothesis, here, we screened global HCC datasets and performed multiscale embedded gene co-expression network analysis to identify the potential functional modules of differentially expressed STAG1 co-expressed genes. The putative transcriptional targets of STAG1 were identified using chromatin immunoprecipitation followed by high-throughput DNA sequencing. The cohesin-associated gene score (CAGS) was quantified using the The Cancer Genome Atlas HCC cohort and single-sample gene set enrichment analysis. Distinct cohesin-associated gene patterns were identified by calculating the euclidean distance of each patient. We assessed the potential ability of the CAGS in predicting immune checkpoint blockade (ICB) treatment response using IMvigor210 and GSE78220 cohorts. STAG1 was upregulated in 3313 HCC tissue samples compared with 2692 normal liver tissue samples (standard mean difference = 0.54). A total of three cohesin-associated gene patterns were identified, where cluster 2 had a high TP53 mutated rate and a poor survival outcome. Low CAGS predicted a significant survival advantage but presaged poor immunotherapy response. Differentially expressed STAG1 co-expression genes were enriched in the mitotic cell cycle, lymphocyte activation, and blood vessel development. PDS5A and PDGFRA were predicted as the downstream transcriptional targets of STAG1. In summary, STAG1 is significantly upregulated in global HCC tissue samples and may participate in blood vessel development and the mitotic cell cycle. A cohesin-associated gene scoring system may have potential to predict the ICB response.

Maintaining Transcriptional Specificity Through Mitosis

Virtually all cell types have the same DNA, yet each type exhibits its own cell-specific pattern of gene expression. During the brief period of mitosis, the chromosomes exhibit changes in protein composition and modifications, a marked condensation, and a consequent reduction in transcription. Yet as cells exit mitosis, they reactivate their cell-specific programs with high fidelity. Initially, the field focused on the subset of transcription factors that are selectively retained in, and hence bookmark, chromatin in mitosis. However, recent studies show that many transcription factors can be retained in mitotic chromatin and that, surprisingly, such retention can be due to nonspecific chromatin binding. Here, we review the latest studies focusing on low-level transcription via promoters, rather than enhancers, as contributing to mitotic memory, as well as new insights into chromosome structure dynamics, histone modifications, cell cycle signaling, and nuclear envelope proteins that together ensure the fidelity of gene expression through a round of mitosis.

SETDB1 acts as a topological accessory to Cohesin via an H3K9me3-independent, genomic shunt for regulating cell fates

SETDB1 is a key regulator of lineage-specific genes and endogenous retroviral elements (ERVs) through its deposition of repressive H3K9me3 mark. Apart from its H3K9me3 regulatory role, SETDB1 has seldom been studied in terms of its other potential regulatory roles. To investigate this, a genomic survey of SETDB1 binding in mouse embryonic stem cells across multiple libraries was conducted, leading to the unexpected discovery of regions bereft of common repressive histone marks (H3K9me3, H3K27me3). These regions were enriched with the CTCF motif that is often associated with the topological regulator Cohesin. Further profiling of these non-H3K9me3 regions led to the discovery of a cluster of non-repeat loci that were co-bound by SETDB1 and Cohesin. These regions, which we named DiSCs (domains involving SETDB1 and Cohesin) were seen to be proximal to the gene promoters involved in embryonic stem cell pluripotency and lineage development. Importantly, it was found that SETDB1-Cohesin co-regulate target gene expression and genome topology at these DiSCs. Depletion of SETDB1 led to localized dysregulation of Cohesin binding thereby locally disrupting topological structures. Dysregulated gene expression trends revealed the importance of this cluster in ES cell maintenance as well as at gene 'islands' that drive differentiation to other lineages. The 'unearthing' of the DiSCs thus unravels a unique topological and transcriptional axis of control regulated chiefly by SETDB1.

Large-scale multi-omics analysis suggests specific roles for intragenic cohesin in transcriptional regulation

Cohesin, an essential protein complex for chromosome segregation, regulates transcription through a variety of mechanisms. It is not a trivial task to assign diverse cohesin functions. Moreover, the context-specific roles of cohesin-mediated interactions, especially on intragenic regions, have not been thoroughly investigated. Here we perform a comprehensive characterization of cohesin binding sites in several human cell types. We integrate epigenomic, transcriptomic and chromatin interaction data to explore the context-specific functions of intragenic cohesin related to gene activation. We identify a specific subset of cohesin binding sites, decreased intragenic cohesin sites (DICs), which are negatively correlated with transcriptional regulation. A subgroup of DICs is enriched with enhancer markers and RNA polymerase II, while the others are more correlated to chromatin architecture. DICs are observed in various cell types, including cells from patients with cohesinopathy. We also implement machine learning to our data and identified genomic features for isolating DICs from all cohesin sites. These results suggest a previously unidentified function of cohesin on intragenic regions for transcriptional regulation.

Alcohol induced increases in sperm Histone H3 lysine 4 trimethylation correlate with increased placental CTCF occupancy and altered developmental programming

Using a mouse model, studies by our group reveal that paternal preconception alcohol intake affects offspring fetal-placental growth, with long-lasting consequences on adult metabolism. Here, we tested the hypothesis that chronic preconception male alcohol exposure impacts histone enrichment in sperm and that these changes are associated with altered developmental programming in the placenta. Using chromatin immunoprecipitation, we find alcohol-induced increases in sperm histone H3 lysine 4 trimethylation (H3K4me3) that map to promoters and presumptive enhancer regions enriched in genes driving neurogenesis and craniofacial development. Given the colocalization of H3K4me3 with the chromatin binding factor CTCF across both sperm and embryos, we next examined CTCF localization in the placenta. We find global changes in CTCF binding within placentae derived from the male offspring of alcohol-exposed sires. Furthermore, altered CTCF localization correlates with dysregulated gene expression across multiple gene clusters; however, these transcriptional changes only occur in male offspring. Finally, we identified a correlation between genomic regions exhibiting alcohol-induced increases in sperm H3K4me3 and increased CTCF binding in male placentae. Collectively, our analysis demonstrates that the chromatin landscape of sperm is sensitive to chronic alcohol exposure and that a subset of these affected regions exhibits increased placental CTCF enrichment.

ESR1 mutant breast cancers show elevated basal cytokeratins and immune activation

Estrogen receptor alpha (ER/ESR1) is frequently mutated in endocrine resistant ER-positive (ER+) breast cancer and linked to ligand-independent growth and metastasis. Despite the distinct clinical features of ESR1 mutations, their role in intrinsic subtype switching remains largely unknown. Here we find that ESR1 mutant cells and clinical samples show a significant enrichment of basal subtype markers, and six basal cytokeratins (BCKs) are the most enriched genes. Induction of BCKs is independent of ER binding and instead associated with chromatin reprogramming centered around a progesterone receptor-orchestrated insulated neighborhood. BCK-high ER+ primary breast tumors exhibit a number of enriched immune pathways, shared with ESR1 mutant tumors. S100A8 and S100A9 are among the most induced immune mediators and involve in tumor-stroma paracrine crosstalk inferred by single-cell RNA-seq from metastatic tumors. Collectively, these observations demonstrate that ESR1 mutant tumors gain basal features associated with increased immune activation, encouraging additional studies of immune therapeutic vulnerabilities.

The glucocorticoid receptor associates with the cohesin loader NIPBL to promote long-range gene regulation

The cohesin complex is central to chromatin looping, but mechanisms by which these long-range chromatin interactions are formed and persist remain unclear. We demonstrate that interactions between a transcription factor (TF) and the cohesin loader NIPBL regulate enhancer-dependent gene activity. Using mass spectrometry, genome mapping, and single-molecule tracking methods, we demonstrate that the glucocorticoid (GC) receptor (GR) interacts with NIPBL and the cohesin complex at the chromatin level, promoting loop extrusion and long-range gene regulation. Real-time single-molecule experiments show that loss of cohesin markedly diminishes the concentration of TF molecules at specific nuclear confinement sites, increasing TF local concentration and promoting gene regulation. Last, patient-derived acute myeloid leukemia cells harboring cohesin mutations exhibit a reduced response to GCs, suggesting that the GR-NIPBL-cohesin interaction is defective in these patients, resulting in poor response to GC treatment.

Evaluating Expression of the STAG1 Gene as a Potential Breast Cancer Biomarker

STAG proteins, which are part of the cohesin complex and encoded by the STAG genes, are known as Irr1/Scc3 in yeast and as SA/STAG/stromalin in mammals. There are more variants as there are alternate splice sites, maybe three open reading frames (ORFs) code for three main proteins, including: SA1 (STAG1), SA2 (STAG2) and SA3 (STAG3). The cohesin protein complex has various essential roles in eukaryotic cell biology. This study compared the expression of the STAG1 gene in four different breast cancer cell lines, including: MCF-7, T-47D, MDA-MB-468, and MDA-MB-231 and normal breast tissue. RNA was extracted from these cell lines and mRNA was converted to cDNA, and then expression of the STAG1 gene was quantified by three sets of specific primer pairs using Real Time-quantitative PCR (RT-qPCR). The findings show significantly different over-expression of STAG1 in these cancer cell lines in comparison with the normal tissue, and the cell lines were different in their expression levels. In conclusion, the STAG1 gene can be postulated as a candidate breast cancer biomarker that needs to be further evaluated in breast tumor biopsies.

Master lineage transcription factors anchor trans mega transcriptional complexes at highly accessible enhancer sites to promote long-range chromatin clustering and transcription of distal target genes

The term 'super enhancers' (SE) has been widely used to describe stretches of closely localized enhancers that are occupied collectively by large numbers of transcription factors (TFs) and co-factors, and control the transcription of highly-expressed genes. Through integrated analysis of >600 DNase-seq, ChIP-seq, GRO-seq, STARR-seq, RNA-seq, Hi-C and ChIA-PET data in five human cancer cell lines, we identified a new class of autonomous SEs (aSEs) that are excluded from classic SE calls by the widely used Rank Ordering of Super-Enhancers (ROSE) method. TF footprint analysis revealed that compared to classic SEs and regular enhancers, aSEs are tightly bound by a dense array of master lineage TFs, which serve as anchors to recruit additional TFs and co-factors in trans. In addition, aSEs are preferentially enriched for Cohesins, which likely involve in stabilizing long-distance interactions between aSEs and their distal target genes. Finally, we showed that aSEs can be reliably predicted using a single DNase-seq data or combined with Mediator and/or P300 ChIP-seq. Overall, our study demonstrates that aSEs represent a unique class of functionally important enhancer elements that distally regulate the transcription of highly expressed genes.

Mouse Embryonic Fibroblast Adipogenic Potential: A Comprehensive Transcriptome Analysis

Our understanding of adipose tissue has progressed from an inert tissue for energy storage to be one of the largest endocrine organs regulating metabolic homoeostasis through its ability to synthesize and release various adipokines that regulate a myriad of pathways. The field of adipose tissue biology is growing due to this association with various chronic metabolic diseases. An important process in the regulation of adipose tissue biology is adipogenesis, which is the formation of new adipocytes. Investigating adipogenesis in vitro is currently a focus for identifying factors that might be utilized in clinically. A powerful tool for such work is high-throughput sequencing which can rapidly identify changes at gene expression level. Various cell models exist for studying adipogenesis and has been used in high-throughput studies, yet little is known about transcriptome profile that underlies adipogenesis in mouse embryonic fibroblasts. This study utilizes RNA-sequencing and computational analysis with DESeq2, gene ontology, protein–protein networks, and robust rank analysis to understand adipogenesis in mouse embryonic fibroblasts in-depth. Our analyses confirmed the requirement of mitotic clonal expansion prior to adipogenesis in this cell model and highlight the role of Cebpa and Cebpb in regulating adipogenesis through interactions of large numbers of genes.

Inner nuclear protein Matrin-3 coordinates cell differentiation by stabilizing chromatin architecture

Precise control of gene expression during differentiation relies on the interplay of chromatin and nuclear structure. Despite an established contribution of nuclear membrane proteins to developmental gene regulation, little is known regarding the role of inner nuclear proteins. Here we demonstrate that loss of the nuclear scaffolding protein Matrin-3 (Matr3) in erythroid cells leads to morphological and gene expression changes characteristic of accelerated maturation, as well as broad alterations in chromatin organization similar to those accompanying differentiation. Matr3 protein interacts with CTCF and the cohesin complex, and its loss perturbs their occupancy at a subset of sites. Destabilization of CTCF and cohesin binding correlates with altered transcription and accelerated differentiation. This association is conserved in embryonic stem cells. Our findings indicate Matr3 negatively affects cell fate transitions and demonstrate that a critical inner nuclear protein impacts occupancy of architectural factors, culminating in broad effects on chromatin organization and cell differentiation.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

GSK3-mediated phosphorylation of DEK3 regulates chromatin accessibility and stress tolerance in Arabidopsis

Chromatin dynamics enable the precise control of transcriptional programmes. The balance between restricting and opening of regulatory sequences on the DNA needs to be adjusted to prevailing conditions and is fine-tuned by chromatin remodelling proteins. DEK is an evolutionarily conserved chromatin architectural protein regulating important chromatin-related processes. However, the molecular link between DEK-induced chromatin reconfigurations and upstream signalling events remains unknown. Here, we show that ASKβ/AtSK31 is a salt stress-activated glycogen synthase kinase 3 (GSK3) from Arabidopsis thaliana that phosphorylates DEK3. This specific phosphorylation alters nuclear DEK3 protein complex composition and affects nucleosome occupancy and chromatin accessibility that is translated into changes in gene expression, contributing to salt stress tolerance. These findings reveal that DEK3 phosphorylation is critical for chromatin function and cellular stress response and provide a mechanistic example of how GSK3-based signalling is directly linked to chromatin, facilitating a transcriptional response.

Visualizing looping of two endogenous genomic loci using synthetic zinc-finger proteins with anti-FLAG and anti-HA frankenbodies in living cells

In eukaryotic nuclei, chromatin loops mediated through cohesin are critical structures that regulate gene expression and DNA replication. Here, we demonstrate a new method to see endogenous genomic loci using synthetic zinc-finger proteins harboring repeat epitope tags (ZF probes) for signal amplification via binding of tag-specific intracellular antibodies, or frankenbodies, fused with fluorescent proteins. We achieve this in two steps: First, we develop an anti-FLAG frankenbody that can bind FLAG-tagged proteins in diverse live-cell environments. The anti-FLAG frankenbody complements the anti-HA frankenbody, enabling two-color signal amplification from FLAG- and HA-tagged proteins. Second, we develop a pair of cell-permeable ZF probes that specifically bind two endogenous chromatin loci predicted to be involved in chromatin looping. By coupling our anti-FLAG and anti-HA frankenbodies with FLAG- and HA-tagged ZF probes, we simultaneously see the dynamics of the two loci in single living cells. This shows a close association between the two loci in the majority of cells, but the loci markedly separate from the triggered degradation of the cohesin subunit RAD21. Our ability to image two endogenous genomic loci simultaneously in single living cells provides a proof of principle that ZF probes coupled with frankenbodies are useful new tools for exploring genome dynamics in multiple colors.

Characterization of a functional endothelial super-enhancer that regulates ADAMTS18 and angiogenesis

Super-enhancers are clusters of enhancers associated with cell lineage. They can be powerful gene-regulators and may be useful in cell-type specific viral-vector development. Here, we have screened for endothelial super-enhancers and identified an enhancer from within a cluster that conferred 5–70-fold increase in transgene expression. Importantly, CRISPR/Cas9 deletion of enhancers demonstrated regulation of ADAMTS18, corresponding to evidence of chromatin contacts between these genomic regions. Cell division-related pathways were primarily affected by the enhancer deletions, which correlated with significant reduction in cell proliferation. Furthermore, we observed changes in angiogenesis-related genes consistent with the endothelial specificity of this SE. Indeed, deletion of the enhancers affected tube formation, resulting in reduced or shortened sprouts. The super-enhancer angiogenic role is at least partly due to its regulation of ADAMTS18, as siRNA knockdown of ADAMTS18 resulted in significantly shortened endothelial sprouts. Hence, functional characterization of a novel endothelial super-enhancer has revealed substantial downstream effects from single enhancer deletions and led to the discovery of the cis-target gene ADAMTS18 and its role in endothelial function.

BETting on a Transcriptional Deficit as the Main Cause for Cornelia de Lange Syndrome

Cornelia de Lange Syndrome (CdLS) is a human developmental syndrome with complex multisystem phenotypic features. It has been traditionally considered a cohesinopathy together with other phenotypically related diseases because of their association with mutations in subunits of the cohesin complex. Despite some overlap, the clinical manifestations of cohesinopathies vary considerably and, although their precise molecular mechanisms are not well defined yet, the potential pathomechanisms underlying these diverse developmental defects have been theoretically linked to alterations of the cohesin complex function. The cohesin complex plays a critical role in sister chromatid cohesion, but this function is not affected in CdLS. In the last decades, a non-cohesion-related function of this complex on transcriptional regulation has been well established and CdLS pathoetiology has been recently associated to gene expression deregulation. Up to 70% of CdLS cases are linked to mutations in the cohesin-loading factor NIPBL, which has been shown to play a prominent function on chromatin architecture and transcriptional regulation. Therefore, it has been suggested that CdLS can be considered a transcriptomopathy. Actually, CdLS-like phenotypes have been associated to mutations in chromatin-associated proteins, as KMT2A, AFF4, EP300, TAF6, SETD5, SMARCB1, MAU2, ZMYND11, MED13L, PHIP, ARID1B, NAA10, BRD4 or ANKRD11, most of which have no known direct association with cohesin. In the case of BRD4, a critical highly investigated transcriptional coregulator, an interaction with NIPBL has been recently revealed, providing evidence on their cooperation in transcriptional regulation of developmentally important genes. This new finding reinforces the notion of an altered gene expression program during development as the major etiological basis for CdLS. In this review, we intend to integrate the recent available evidence on the molecular mechanisms underlying the clinical manifestations of CdLS, highlighting data that favors a transcription-centered framework, which support the idea that CdLS could be conceptualized as a transcriptomopathy.

LIM homeodomain proteins and associated partners: Then and now

The field of molecular embryology started around 1990 by identifying new genes and analyzing their functions in early vertebrate embryogenesis. Those genes encode transcription factors, signaling molecules, their regulators, etc. Most of those genes are relatively highly expressed in specific regions or exhibit dramatic phenotypes when ectopically expressed or mutated. This review focuses on one of those genes, Lim1/Lhx1, which encodes a transcription factor. Lim1/Lhx1 is a member of the LIM homeodomain (LIM-HD) protein family, and its intimate partner, Ldb1/NLI, binds to two tandem LIM domains of LIM-HDs. The most ancient LIM-HD protein and its partnership with Ldb1 were innovated in the metazoan ancestor by gene fusion combining LIM domains and a homeodomain and by creating the LIM domain-interacting domain (LID) in ancestral Ldb, respectively. The LIM domain has multiple interacting interphases, and Ldb1 has a dimerization domain (DD), the LID, and other interacting domains that bind to Ssbp2/3/4 and the boundary factor, CTCF. By means of these domains, LIM-HD-Ldb1 functions as a hub protein complex, enabling more intricate and elaborate gene regulation. The common, ancestral role of LIM-HD proteins is neuron cell-type specification. Additionally, Lim1/Lhx1 serves crucial roles in the gastrula organizer and in kidney development. Recent studies using Xenopus embryos have revealed Lim1/Lhx1 functions and regulatory mechanisms during development and regeneration, providing insight into evolutionary developmental biology, functional genomics, gene regulatory networks, and regenerative medicine. In this review, we also discuss recent progress at unraveling participation of Ldb1, Ssbp, and CTCF in enhanceosomes, long-distance enhancer-promoter interactions, and trans-interactions between chromosomes.

Enhancer RNAs Mediate Estrogen-Induced Decommissioning of Selective Enhancers by Recruiting ERα and Its Cofactor

The function of enhancer RNAs (eRNAs) in transcriptional regulation remains obscure. By analyzing the genome-wide nascent transcript profiles in breast cancer cells, we identify a special group of eRNAs that are essential for estrogen-induced transcriptional repression. Using eRNAs of TM4SF1 and EFEMP1 as the paradigms, we find that these RNA molecules not only stabilize promoter-enhancer interactions but also recruit liganded estrogen receptor α (ERα) to particular enhancer regions, facilitate the formation of a functional transcriptional complex, and cause gene silencing. Interestingly, ERα is shown to directly bind with eRNAs by its DNA-binding domain. These eRNAs help with the formation of a specific ERα-centered transcriptional complex and promote the association of the histone demethylase KDM2A, which dismisses RNA polymerase II from designated enhancers and suppresses the transcription of target genes. Our work demonstrates a complete mechanism underlying the action of eRNAs in modulating and refining the locus-specific transcriptional program.

Yang et al. identified a group of eRNAs that are essential for estrogen-induced transcriptional repression by assisting with the chromatin recruitment of ERα through binding to its DNA-binding domain and facilitating the interaction of ERα with its cofactors, which leads to the dismissal of RNA polymerase II.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

At the Crossroad of Gene Regulation and Genome Organization: Potential Roles for ATP-Dependent Chromatin Remodelers in the Regulation of CTCF-Mediated 3D Architecture

Simple Summary

The way DNA is packaged in the nucleus of a cell is important for when and how genes are expressed. There are many levels of packaging, and new techniques have revealed that long-range interactions are important for both promoting and restricting the transcription of genes. Some long-range interactions are mediated by physical loops in the genome where, like a rubber band, the ring-shaped cohesin complex loops sections of DNA bound by CCCTC-binding factor (CTCF). Both cohesin and CTCF act on DNA, and increasing evidence indicates that their function is inhibited by nucleosomes bound to the DNA. In this review, we summarize the current knowledge of how individual chromatin remodelers, which utilize ATP to move nucleosomes on DNA, facilitate or inhibit cohesin/CTCF-dependent looping interactions.

Abstract

In higher order organisms, the genome is assembled into a protein-dense structure called chromatin. Chromatin is spatially organized in the nucleus through hierarchical folding, which is tightly regulated both in cycling cells and quiescent cells. Assembly and folding are not one-time events in a cell’s lifetime; rather, they are subject to dynamic shifts to allow changes in transcription, DNA replication, or DNA damage repair. Chromatin is regulated at many levels, and recent tools have permitted the elucidation of specific factors involved in the maintenance and regulation of the three-dimensional (3D) genome organization. In this review/perspective, we aim to cover the potential, but relatively unelucidated, crosstalk between 3D genome architecture and the ATP-dependent chromatin remodelers with a specific focus on how the architectural proteins CTCF and cohesin are regulated by chromatin remodeling.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

A cohesive look at leukemogenesis: The cohesin complex and other driving mutations in AML

Acute myeloid leukemia (AML) affects tens of thousands of patients a year, yet survival rates are as low as 25% in certain populations. This poor survival rate is partially due to the vast genetic diversity of the disease. Rarely do 2 patients with AML have the same mutational profile, which makes the development of targeted therapies particularly challenging. However, a set of recurrent mutations in chromatin modifiers have been identified in many patients, including mutations in the cohesin complex, which have been identified in up to 20% of cases. Interestingly, the canonical function of the cohesin complex in establishing sister chromatid cohesin during mitosis is unlikely to be the affected role in leukemogenesis. Instead, the cohesin complex's role in DNA looping and gene regulation likely facilitates disease. The epigenetic mechanisms by which cohesin complex mutations promote leukemia are not completely elucidated, but alterations of enhancer-promoter interactions and differential histone modifications have been shown to drive oncogenic gene expression changes. Such changes commonly include HoxA upregulation, which may represent a common pathway that could be therapeutically targeted. As cohesin mutations rarely occur alone, examining the impact of common co-occurring mutations, including those in NPM1, the core-binding factor complex, FLT3, and ASXL1, will yield additional insight. While further study of these mutational interactions is required, current research suggests that the use of combinatorial genetics could be the key to uncovering new targets, allowing for the treatment of AML patients based on their individual genetic profiles.

Tracking telomere fusions through crisis reveals conflict between DNA transcription and the DNA damage response

Identifying attributes that distinguish pre-malignant from senescent cells provides opportunities for targeted disease eradication and revival of anti-tumour immunity. We modelled a telomere-driven crisis in four human fibroblast lines, sampling at multiple time points to delineate genomic rearrangements and transcriptome developments that characterize the transition from dynamic proliferation into replicative crisis. Progression through crisis was associated with abundant intra-chromosomal telomere fusions with increasing asymmetry and reduced microhomology usage, suggesting shifts in DNA repair capacity. Eroded telomeres also fused with genomic loci actively engaged in transcription, with particular enrichment in long genes. Both gross copy number alterations and transcriptional responses to crisis likely underpin the elevated frequencies of telomere fusion with chromosomes 9, 16, 17, 19 and most exceptionally, chromosome 12. Juxtaposition of crisis-regulated genes with loci undergoing de novo recombination exposes the collusive contributions of cellular stress responses to the evolving cancer genome.

ESR1 ChIP-Seq Identifies Distinct Ligand-Free ESR1 Genomic Binding Sites in Human Hepatocytes and Liver Tissue

The estrogen receptor alpha (ESR1) is an important gene transcriptional regulator, known to mediate the effects of estrogen. Canonically, ESR1 is activated by its ligand estrogen. However, the role of unliganded ESR1 in transcriptional regulation has been gaining attention. We have recently shown that ligand-free ESR1 is a key regulator of several cytochrome P450 (CYP) genes in the liver, however ligand-free ESR1 has not been characterized genome-wide in the human liver. To address this, ESR1 ChIP-Seq was conducted in human liver samples and in hepatocytes with or without 17beta-estradiol (E2) treatment. We identified both ligand-dependent and ligand-independent binding sites throughout the genome. These two ESR1 binding categories showed different genomic localization, pathway enrichment, and cofactor colocalization, indicating different ESR1 regulatory function depending on ligand availability. By analyzing existing ESR1 data from additional human cell lines, we uncovered a potential ligand-independent ESR1 activity, namely its co-enrichment with the zinc finger protein 143 (ZNF143). Furthermore, we identified ESR1 binding sites near many gene loci related to drug therapy, including the CYPs. Overall, this study shows distinct ligand-free and ligand-bound ESR1 chromatin binding profiles in the liver and suggests the potential broad influence of ESR1 in drug metabolism and drug therapy.

The relationship between genome structure and function

Precise patterns of gene expression in metazoans are controlled by three classes of regulatory elements: promoters, enhancers and boundary elements. During differentiation and development, these elements form specific interactions in dynamic higher-order chromatin structures. However, the relationship between genome structure and its function in gene regulation is not completely understood. Here we review recent progress in this field and discuss whether genome structure plays an instructive role in regulating gene expression or is a reflection of the activity of the regulatory elements of the genome.

Cohesin subunit RAD21: From biology to disease

RAD21 (also known as KIAA0078, NXP1, HR21, Mcd1, Scc1, and hereafter called RAD21), an essential gene, encodes a DNA double-strand break (DSB) repair protein that is evolutionarily conserved in all eukaryotes from budding yeast to humans. RAD21 protein is a structural component of the highly conserved cohesin complex consisting of RAD21, SMC1a, SMC3, and SCC3 [STAG1 (SA1) and STAG2 (SA2) in metazoans] proteins, involved in sister chromatid cohesion. This function is essential for proper chromosome segregation, post-replicative DNA repair, and prevention of inappropriate recombination between repetitive regions. In interphase, cohesin also functions in the control of gene expression by binding to numerous sites within the genome. In addition to playing roles in the normal cell cycle and DNA DSB repair, RAD21 is also linked to the apoptotic pathways. Germline heterozygous or homozygous missense mutations in RAD21 have been associated with human genetic disorders, including developmental diseases such as Cornelia de Lange syndrome (CdLS) and chronic intestinal pseudo-obstruction (CIPO) called Mungan syndrome, respectively, and collectively termed as cohesinopathies. Somatic mutations and amplification of the RAD21 have also been widely reported in both human solid and hematopoietic tumors. Considering the role of RAD21 in a broad range of cellular processes that are hot spots in neoplasm, it is not surprising that the deregulation of RAD21 has been increasingly evident in human cancers. Herein, we review the biology of RAD21 and the cellular processes that this important protein regulates and discuss the significance of RAD21 deregulation in cancer and cohesinopathies.

Cohesin-Dependent and -Independent Mechanisms Mediate Chromosomal Contacts between Promoters and Enhancers

It is currently assumed that 3D chromosomal organization plays a central role in transcriptional control. However, depletion of cohesin and CTCF affects the steady-state levels of only a minority of transcripts. Here, we use high-resolution Capture Hi-C to interrogate the dynamics of chromosomal contacts of all annotated human gene promoters upon degradation of cohesin and CTCF. We show that a majority of promoter-anchored contacts are lost in these conditions, but many contacts with distinct properties are maintained, and some new ones are gained. The rewiring of contacts between promoters and active enhancers upon cohesin degradation associates with rapid changes in target gene transcription as detected by SLAM sequencing (SLAM-seq). These results provide a mechanistic explanation for the limited, but consistent, effects of cohesin and CTCF depletion on steady-state transcription and suggest the existence of both cohesin-dependent and -independent mechanisms of enhancer-promoter pairing.

Impact of 3D genome organization, guided by cohesin and CTCF looping, on sex-biased chromatin interactions and gene expression in mouse liver

Several thousand sex-differential distal enhancers have been identified in mouse liver; however, their links to sex-biased genes and the impact of any sex-differences in nuclear organization and chromatin interactions are unknown. To address these issues, we first characterized 1847 mouse liver genomic regions showing significant sex differential occupancy by cohesin and CTCF, two key 3D nuclear organizing factors. These sex-differential binding sites were primarily distal to sex-biased genes but rarely generated sex-differential TAD (topologically associating domain) or intra-TAD loop anchors, and were sometimes found in TADs without sex-biased genes. A substantial subset of sex-biased cohesin-non-CTCF binding sites, but not sex-biased cohesin-and-CTCF binding sites, overlapped sex-biased enhancers. Cohesin depletion reduced the expression of male-biased genes with distal, but not proximal, sex-biased enhancers by >10-fold, implicating cohesin in long-range enhancer interactions regulating sex-biased genes. Using circularized chromosome conformation capture-based sequencing (4C-seq), we showed that sex differences in distal sex-biased enhancer-promoter interactions are common. Intra-TAD loops with sex-independent cohesin-and-CTCF anchors conferred sex specificity to chromatin interactions indirectly, by insulating sex-biased enhancer-promoter contacts and by bringing sex-biased genes into closer proximity to sex-biased enhancers. Furthermore, sex-differential chromatin interactions involving sex-biased gene promoters, enhancers, and lncRNAs were associated with sex-biased binding of cohesin and/or CTCF. These studies elucidate how 3D genome organization impacts sex-biased gene expression in a non-reproductive tissue through both direct and indirect effects of cohesin and CTCF looping on distal enhancer interactions with sex-differentially expressed genes.

Epigenomics-based identification of oestrogen-regulated long noncoding RNAs in ER+ breast cancer

Breast cancer is one of the most prevalent cancers in women worldwide. Through the regulation of many coding and non-coding target genes, oestrogen (E₂ or 17β-oestradiol) and its nuclear receptor ERα play important roles in breast cancer development and progression. Despite the astounding advances in our understanding of oestrogen-regulated coding genes over the past decades, our knowledge on oestrogen-regulated non-coding targets has just begun to expand. Here we leverage epigenomic approaches to systematically analyse oestrogen-regulated long non-coding RNAs (lncRNAs). Similar to the coding targets of ERα, the transcription of oestrogen-regulated lncRNAs correlates with the activation status of ERα enhancers, measured by eRNA production, chromatin accessibility, and the occupancy of the enhancer regulatory components including P300, MED1, and ARID1B. Our 3D chromatin architecture analyses suggest that lncRNAs and their neighbouring E₂-resonsive coding genes, exemplified by LINC00160 and RUNX1, might be regulated as a 3D structural unit resulted from enhancer-promoter interactions. Finally, we evaluated the expression levels of LINC00160 and RUNX1 in various types of breast cancer and found that their expression positively correlated with the survival rate in ER+ breast cancer patients, implying that the oestrogen-regulated LINC00160 and its neighbouring RUNX1 might represent potential biomarkers for ER+ breast cancers.

Downexpression of HSD17B6 correlates with clinical prognosis and tumor immune infiltrates in hepatocellular carcinoma

Background

Hydroxysteroid 17-Beta Dehydrogenase 6 (HSD17B6), a key protein involved in synthetizing dihydrotestosterone, is abundant in the liver. Previous studies have suggested a role for dihydrotestosterone in modulating progress of various malignancies, and HSD17B6 dysfunction was associated with lung cancer and prostate cancer. However, little is known about the detailed role of HSD17B6 in hepatocellular carcinoma (HCC).

Methods

Clinical implication and survival data related to HSD17B6 expression in patients with HCC were obtained through TCGA, ICGC, ONCOMINE, GEO and HPA databases. Survival analysis plots were drawn with Kaplan–Meier Plotter. The ChIP-seq data were obtained from Cistrome DB. Protein–Protein Interaction and gene functional enrichment analyses were performed in STRING database. The correlations between HSD17B6 and tumor immune infiltrates was investigated via TIMER and xCell. The proliferation, migration and invasion of liver cancer cells transfected with HSD17B6 were evaluated by the CCK8 assay, wound healing test and transwell assay respectively. Expression of HSD17B6, TGFB1 and PD-L1 were assessed by quantitative RT-PCR.

Results

HSD17B6 expression was lower in HCC compared to normal liver and correlated with tumor stage and grade. Lower expression of HSD17B6 was associated with worse OS, PFS, RFS and DSS in HCC patients. HNF4A bound to enhancer and promoter regions of HSD17B6 gene, activating its transcription, and DNA methylation of HSD17B6 promoter negatively controlled the expression. HSD17B6 and its interaction partners were involved in androgen metabolism and biosynthesis in liver. HSD17B6 inhibited tumor cell proliferation, migration and invasion in liver cancer cells and low expression of HSD17B6 correlated with high immune cells infiltration, relative reduction of immune responses and multiple immune checkpoint genes expression in HCC, probably by regulating the expression of TGFB1.

Conclusions

This study indicate that HSD17B6 could be a new biomarker for the prognosis of HCC and an important negative regulator of immune responses in HCC.

Dynamics of the 4D genome during in vivo lineage specification and differentiation

Mammalian gene expression patterns are controlled by regulatory elements, which interact within topologically associating domains (TADs). The relationship between activation of regulatory elements, formation of structural chromatin interactions and gene expression during development is unclear. Here, we present Tiled-C, a low-input chromosome conformation capture (3C) technique. We use this approach to study chromatin architecture at high spatial and temporal resolution through in vivo mouse erythroid differentiation. Integrated analysis of chromatin accessibility and single-cell expression data shows that regulatory elements gradually become accessible within pre-existing TADs during early differentiation. This is followed by structural re-organization within the TAD and formation of specific contacts between enhancers and promoters. Our high-resolution data show that these enhancer-promoter interactions are not established prior to gene expression, but formed gradually during differentiation, concomitant with progressive upregulation of gene activity. Together, these results provide new insight into the close, interdependent relationship between chromatin architecture and gene regulation during development.

Estrogen induces dynamic ERα and RING1B recruitment to control gene and enhancer activities in luminal breast cancer

RING1B, a core Polycomb repressive complex 1 subunit, is a histone H2A ubiquitin ligase essential for development. RING1B is overexpressed in patients with luminal breast cancer (BC) and recruited to actively transcribed genes and enhancers co-occupied by the estrogen receptor α (ERα). Whether ERα-induced transcriptional programs are mediated by RING1B is not understood. We show that prolonged estrogen administration induces transcriptional output and chromatin landscape fluctuations. RING1B loss impairs full estrogen-mediated gene expression and chromatin accessibility for key BC transcription factors. These effects were mediated, in part, by RING1B enzymatic activity and nucleosome binding functions. RING1B is recruited in a cyclic manner to ERα, FOXA1, and GRHL2 cobound sites and regulates estrogen-induced enhancers and ERα recruitment. Last, ChIP exo revealed multiple binding events of these factors at single-nucleotide resolution, including RING1B occupancy approximately 10 base pairs around ERα bound sites. We propose RING1B as a key regulator of the dynamic, liganded-ERα transcriptional regulatory circuit in luminal BC.

Redundant and specific roles of cohesin STAG subunits in chromatin looping and transcriptional control

Cohesin is a ring-shaped multiprotein complex that is crucial for 3D genome organization and transcriptional regulation during differentiation and development. It also confers sister chromatid cohesion and facilitates DNA damage repair. Besides its core subunits SMC3, SMC1A, and RAD21, cohesin in somatic cells contains one of two orthologous STAG subunits, STAG1 or STAG2. How these variable subunits affect the function of the cohesin complex is still unclear. STAG1- and STAG2-cohesin were initially proposed to organize cohesion at telomeres and centromeres, respectively. Here, we uncover redundant and specific roles of STAG1 and STAG2 in gene regulation and chromatin looping using HCT116 cells with an auxin-inducible degron (AID) tag fused to either STAG1 or STAG2. Following rapid depletion of either subunit, we perform high-resolution Hi-C, gene expression, and sequential ChIP studies to show that STAG1 and STAG2 do not co-occupy individual binding sites and have distinct ways by which they affect looping and gene expression. These findings are further supported by single-molecule localizations via direct stochastic optical reconstruction microscopy (dSTORM) super-resolution imaging. Since somatic and congenital mutations of the STAG subunits are associated with cancer (STAG2) and intellectual disability syndromes with congenital abnormalities (STAG1 and STAG2), we verified STAG1-/STAG2-dependencies using human neural stem cells, hence highlighting their importance in particular disease contexts.