Computational prediction of CTCF/cohesin-based intra-TAD loops that insulate chromatin contacts and gene expression in mouse liver

SMC3 contributes to heart development by regulating super-enhancer associated genes

Abnormal cardiac development has been observed in individuals with Cornelia de Lange syndrome (CdLS) due to mutations in genes encoding members of the cohesin complex. However, the precise role of cohesin in heart development remains elusive. In this study, we aimed to elucidate the indispensable role of SMC3, a component of the cohesin complex, in cardiac development and its underlying mechanism. Our investigation revealed that CdLS patients with SMC3 mutations have high rates of congenital heart disease (CHD). We utilized heart-specific Smc3-knockout (SMC3-cKO) mice, which exhibit varying degrees of outflow tract (OFT) abnormalities, to further explore this relationship. Additionally, we identified 16 rare SMC3 variants with potential pathogenicity in individuals with isolated CHD. By employing single-nucleus RNA sequencing and chromosome conformation capture high-throughput genome-wide translocation sequencing, we revealed that Smc3 deletion downregulates the expression of key genes, including Ets2, in OFT cardiac muscle cells by specifically decreasing interactions between super-enhancers (SEs) and promoters. Notably, Ets2-SE-null mice also exhibit delayed OFT development in the heart. Our research revealed a novel role for SMC3 in heart development via the regulation of SE-associated genes, suggesting its potential relevance as a CHD-related gene and providing crucial insights into the molecular basis of cardiac development.

CatLearning: highly accurate gene expression prediction from histone mark

Histone modifications, known as histone marks, are pivotal in regulating gene expression within cells. The vast array of potential combinations of histone marks presents a considerable challenge in decoding the regulatory mechanisms solely through biological experimental approaches. To overcome this challenge, we have developed a method called CatLearning. It utilizes a modified convolutional neural network architecture with a specialized adaptation Residual Network to quantitatively interpret histone marks and predict gene expression. This architecture integrates long-range histone information up to 500Kb and learns chromatin interaction features without 3D information. By using only one histone mark, CatLearning achieves a high level of accuracy. Furthermore, CatLearning predicts gene expression by simulating changes in histone modifications at enhancers and throughout the genome. These findings help comprehend the architecture of histone marks and develop diagnostic and therapeutic targets for diseases with epigenetic changes.

Systematic decoding of cis gene regulation defines context-dependent control of the multi-gene costimulatory receptor locus in human T cells

Cis-regulatory elements (CREs) interact with trans regulators to orchestrate gene expression, but how transcriptional regulation is coordinated in multi-gene loci has not been experimentally defined. We sought to characterize the CREs controlling dynamic expression of the adjacent costimulatory genes CD28, CTLA4 and ICOS, encoding regulators of T cell-mediated immunity. Tiling CRISPR interference (CRISPRi) screens in primary human T cells, both conventional and regulatory subsets, uncovered gene-, cell subset- and stimulation-specific CREs. Integration with CRISPR knockout screens and assay for transposase-accessible chromatin with sequencing (ATAC-seq) profiling identified trans regulators influencing chromatin states at specific CRISPRi-responsive elements to control costimulatory gene expression. We then discovered a critical CCCTC-binding factor (CTCF) boundary that reinforces CRE interaction with CTLA4 while also preventing promiscuous activation of CD28. By systematically mapping CREs and associated trans regulators directly in primary human T cell subsets, this work overcomes longstanding experimental limitations to decode context-dependent gene regulatory programs in a complex, multi-gene locus critical to immune homeostasis.

Exploring the dynamic three-dimensional chromatin architecture and transcriptional landscape in goose liver tissues underlying metabolic adaptations induced by a high-fat diet

BACKGROUND

Goose, descendants of migratory ancestors, have undergone extensive selective breeding, resulting in their remarkable ability to accumulate fat in the liver and exhibit a high tolerance for significant energy intake. As a result, goose offers an excellent model for studying obesity, metabolic disorders, and liver diseases in mammals. Although the impact of the three-dimensional arrangement of chromatin within the cell nucleus on gene expression and transcriptional regulation is widely acknowledged, the precise functions of chromatin architecture reorganization during fat deposition in goose liver tissues still need to be fully comprehended.

RESULTS

In this study, geese exhibited more pronounced changes in the liver index and triglyceride (TG) content following the consumption of the high-fat diet (HFD) than mice without significant signs of inflammation. Additionally, we performed comprehensive analyses on 10 goose liver tissues (5 HFD, 5 normal), including generating high-resolution maps of chromatin architecture, conducting whole-genome gene expression profiling, and identifying H3K27ac peaks in the livers of geese and mice subjected to the HFD. Our results unveiled a multiscale restructuring of chromatin architecture, encompassing Compartment A/B, topologically associated domains, and interactions between promoters and enhancers. The dynamism of the three-dimensional genome architecture, prompted by the HFD, assumed a pivotal role in the transcriptional regulation of crucial genes. Furthermore, we identified genes that regulate chromatin conformation changes, contributing to the metabolic adaptation process of lipid deposition and hepatic fat changes in geese in response to excessive energy intake. Moreover, we conducted a cross-species analysis comparing geese and mice exposed to the HFD, revealing unique characteristics specific to the goose liver compared to a mouse. These chromatin conformation changes help elucidate the observed characteristics of fat deposition and hepatic fat regulation in geese under conditions of excessive energy intake.

CONCLUSIONS

We examined the dynamic modifications in three-dimensional chromatin architecture and gene expression induced by an HFD in goose liver tissues. We conducted a cross-species analysis comparing that of mice. Our results contribute significant insights into the chromatin architecture of goose liver tissues, offering a novel perspective for investigating mammal liver diseases.

Plasma growth hormone pulses induce male-biased pulsatile chromatin opening and epigenetic regulation in adult mouse liver

Sex differences in plasma growth hormone (GH) profiles, pulsatile in males and persistent in females, regulate sex differences in hepatic STAT5 activation linked to sex differences in gene expression and liver disease susceptibility, but little is understood about the fundamental underlying, GH pattern-dependent regulatory mechanisms. Here, DNase-I hypersensitivity site (DHS) analysis of liver chromatin accessibility in a cohort of 18 individual male mice established that the endogenous male rhythm of plasma GH pulse-stimulated liver STAT5 activation induces dynamic, repeated cycles of chromatin opening and closing at several thousand liver DHS and comprises a novel mechanism conferring male bias to liver chromatin accessibility. Strikingly, a single physiological replacement dose of GH given to hypophysectomized male mice restored, within 30 min, liver STAT5 activity and chromatin accessibility at 83% of the dynamic, pituitary hormone-dependent male-biased DHS. Sex-dependent transcription factor binding patterns and chromatin state analysis identified key genomic and epigenetic features distinguishing this dynamic, STAT5-driven mechanism of male-biased chromatin opening from a second GH-dependent mechanism operative at static male-biased DHS, which are constitutively open in male liver. Dynamic but not static male-biased DHS adopt a bivalent-like epigenetic state in female liver, as do static female-biased DHS in male liver, albeit using distinct repressive histone marks in each sex, namely, H3K9me3 at male-biased DHS in female liver and H3K27me3 at female-biased DHS in male liver. Moreover, sex-biased H3K36me3 marks are uniquely enriched at static sex-biased DHS, which may serve to keep these sex-dependent hepatocyte enhancers free of H3K27me3 repressive marks and thus constitutively open. Pulsatile chromatin opening stimulated by endogenous, physiological hormone pulses is thus one of two distinct GH-determined mechanisms for establishing widespread sex differences in hepatic chromatin accessibility and epigenetic regulation, both closely linked to sex-biased gene transcription and the sexual dimorphism of liver function.

Plasma Growth Hormone Pulses Induce Male-biased Pulsatile Chromatin Opening and Epigenetic Regulation in Adult Mouse Liver

Sex-differences in plasma growth hormone (GH) profiles, pulsatile in males and persistent in females, regulate sex differences in hepatic STAT5 activation linked to sex differences in gene expression and liver disease susceptibility, but little is understood about the fundamental underlying, GH pattern-dependent regulatory mechanisms. Here, DNase hypersensitivity site (DHS) analysis of liver chromatin accessibility in a cohort of 18 individual male mice established that the endogenous male rhythm of plasma GH pulse-stimulated liver STAT5 activation induces dynamic, repeated cycles of chromatin opening and closing at several thousand liver DHS and comprises a novel mechanism conferring male bias to liver chromatin accessibility. Strikingly, a single physiological replacement dose of GH given to hypophysectomized male mice restored, within 30 min, liver STAT5 activity and chromatin accessibility at 83% of the pituitary hormone-dependent dynamic male-biased DHS. Sex-dependent transcription factor binding patterns and chromatin state analysis identified key genomic and epigenetic features distinguishing this dynamic, STAT5-driven mechanism of male-biased chromatin opening from a second GH-dependent mechanism operative at static male-biased DHS, which are constitutively open in male liver. Dynamic but not static male-biased DHS adopt a bivalent-like epigenetic state in female liver, as do static female-biased DHS in male liver, albeit using distinct repressive histone marks in each sex, namely, H3K27me3 at female-biased DHS in male liver, and H3K9me3 at male-biased DHS in female liver. Moreover, sex-biased H3K36me3 marks are uniquely enriched at static sex-biased DHS, which may serve to keep these sex-dependent hepatocyte enhancers free of H3K27me3 repressive marks and thus constitutively open. Pulsatile chromatin opening stimulated by endogenous, physiological hormone pulses is thus one of two distinct GH-determined mechanisms for establishing widespread sex differences in hepatic chromatin accessibility and epigenetic regulation, both closely linked to sex-biased gene transcription and the sexual dimorphism of liver function.

Transposable elements in mammalian chromatin organization

Transposable elements (TEs) are mobile DNA elements that comprise almost 50% of mammalian genomic sequence. TEs are capable of making additional copies of themselves that integrate into new positions in host genomes. This unique property has had an important impact on mammalian genome evolution and on the regulation of gene expression because TE-derived sequences can function as cis-regulatory elements such as enhancers, promoters and silencers. Now, advances in our ability to identify and characterize TEs have revealed that TE-derived sequences also regulate gene expression by both maintaining and shaping 3D genome architecture. Studies are revealing how TEs contribute raw sequence that can give rise to the structures that shape chromatin organization, and thus gene expression, allowing for species-specific genome innovation and evolutionary novelty.

Integrative analysis of single-cell RNA-seq and ATAC-seq reveals heterogeneity of induced pluripotent stem cell-derived hepatic organoids

To gain deeper insights into transcriptomes and epigenomes of organoids, liver organoids from two states (expandable and more differentiated) were subjected to single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) analyses. Mitochondrial gene expression was higher in differentiated than in non-differentiated hepatocytes, with ATAC-seq peaks increasing near the mitochondrial control region. Differentiation of liver organoids resulted in the expression of transcription factors that act as enhancers and repressors. In addition, epigenetic mechanisms regulating the expression of alpha-fetoprotein (AFP) and albumin (ALB) differed in liver organoids and adult liver. Knockdown of PDX1, an essential transcription factor for pancreas development, led to the hepatic maturation of liver organoids through regulation of AFP and ALB expression. This integrative analysis of the transcriptomes and epigenomes of liver organoids at the single-cell level may contribute to a better understanding of the regulatory networks during liver development and the further development of mature in vitro human liver models.

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

Low-affinity CTCF binding drives transcriptional regulation whereas high-affinity binding encompasses architectural functions

CTCF is a DNA-binding protein which plays critical roles in chromatin structure organization and transcriptional regulation; however, little is known about the functional determinants of different CTCF-binding sites (CBS). Using a conditional mouse model, we have identified one set of CBSs that are lost upon CTCF depletion (lost CBSs) and another set that persists (retained CBSs). Retained CBSs are more similar to the consensus CTCF-binding sequence and usually span tandem CTCF peaks. Lost CBSs are enriched at enhancers and promoters and associate with active chromatin marks and higher transcriptional activity. In contrast, retained CBSs are enriched at TAD and loop boundaries. Integration of ChIP-seq and RNA-seq data has revealed that retained CBSs are located at the boundaries between distinct chromatin states, acting as chromatin barriers. Our results provide evidence that transient, lost CBSs are involved in transcriptional regulation, whereas retained CBSs are critical for establishing higher-order chromatin architecture.

Active enhancers strengthen insulation by RNA-mediated CTCF binding at chromatin domain boundaries

Vertebrate genomes are partitioned into chromatin domains or topologically associating domains (TADs), which are typically bound by head-to-head pairs of CTCF binding sites. Transcription at domain boundaries correlates with better insulation; however, it is not known whether the boundary transcripts themselves contribute to boundary function. Here we characterize boundary-associated RNAs genome-wide, focusing on the disease-relevant INK4a/ARF and MYC TAD. Using CTCF site deletions and boundary-associated RNA knockdowns, we observe that boundary-associated RNAs facilitate recruitment and clustering of CTCF at TAD borders. The resulting CTCF enrichment enhances TAD insulation, enhancer-promoter interactions, and TAD gene expression. Importantly, knockdown of boundary-associated RNAs results in loss of boundary insulation function. Using enhancer deletions and CRISPRi of promoters, we show that active TAD enhancers, but not promoters, induce boundary-associated RNA transcription, thus defining a novel class of regulatory enhancer RNAs.

Cohesin and CTCF complexes mediate contacts in chromatin loops depending on nucleosome positions

The spatial organization of the eukaryotic genome plays an important role in regulating transcriptional activity. In the nucleus, chromatin forms loops that assemble into fundamental units called topologically associating domains that facilitate or inhibit long-range contacts. These loops are formed and held together by the ring-shaped cohesin protein complex, and this can involve binding of CCCTC-binding factor (CTCF). High-resolution conformation capture experiments provide the frequency at which two DNA fragments physically associate in three-dimensional space. However, technical limitations of this approach, such as low throughput, low resolution, or noise in contact maps, make data interpretation and identification of chromatin intraloop contacts, e.g., between distal regulatory elements and their target genes, challenging. Herein, an existing coarse-grained model of chromatin at single-nucleosome resolution was extended by integrating potentials describing CTCF and cohesin. We performed replica-exchange Monte Carlo simulations with regularly spaced nucleosomes and experimentally determined nucleosome positions in the presence of cohesin-CTCF, as well as depleted systems as controls. In fully extruded loops caused by the presence of cohesin and CTCF, the number of contacts within the formed loops was increased. The number and types of these contacts were impacted by the nucleosome distribution and loop size. Microloops were observed within cohesin-mediated loops due to thermal fluctuations without additional influence of other factors, and the number, size, and shape of microloops were determined by nucleosome distribution and loop size. Nucleosome positions directly affect the spatial structure and contact probability within a loop, with presumed consequences for transcriptional activity.

Nuclear envelope, chromatin organizers, histones, and DNA: The many achilles heels exploited across cancers

In eukaryotic cells, the genome is organized in the form of chromatin composed of DNA and histones that organize and regulate gene expression. The dysregulation of chromatin remodeling, including the aberrant incorporation of histone variants and their consequent post-translational modifications, is prevalent across cancers. Additionally, nuclear envelope proteins are often deregulated in cancers, which impacts the 3D organization of the genome. Altered nuclear morphology, genome organization, and gene expression are defining features of cancers. With advances in single-cell sequencing, imaging technologies, and high-end data mining approaches, we are now at the forefront of designing appropriate small molecules to selectively inhibit the growth and proliferation of cancer cells in a genome- and epigenome-specific manner. Here, we review recent advances and the emerging significance of aberrations in nuclear envelope proteins, histone variants, and oncohistones in deregulating chromatin organization and gene expression in oncogenesis.

Integrating extrusion complex-associated pattern to predict cell type-specific long-range chromatin loops

The chromatin loop plays a critical role in the study of gene expression and disease. Supervised learning-based algorithms to predict the chromatin loops require large priori information to satisfy the model construction, while the prediction sensitivity of unsupervised learning-based algorithms is still unsatisfactory. Therefore, we propose an unsupervised algorithm, Ecomap-loop. It takes advantage of extrusion complex-associated patterns, including CTCF, RAD21, and SMC enrichments, as well as the orientation distribution of CTCF motif of loops to build feature matrices; then the eigen decomposition model is employed to obtain the cell type-specific loops. We compare the performance of Ecomap-loop with the state-of-the-art unsupervised algorithm using Hi-C, ChIA-PET, expression quantitative trait locus (eQTL), and CRISPR interference (CRISPRi) screen data; the results show that Ecomap-loop achieves the best in four cell types. In addition, the functional analysis reveals the ability of Ecomap-loop to predict active functionality-related and cell type-specific loops.

The HASTER lncRNA promoter is a cis-acting transcriptional stabilizer of HNF1A

The biological purpose of long non-coding RNAs (lncRNAs) is poorly understood. Haploinsufficient mutations in HNF1A homeobox A (HNF1A), encoding a homeodomain transcription factor, cause diabetes mellitus. Here, we examine HASTER, the promoter of an lncRNA antisense to HNF1A. Using mouse and human models, we show that HASTER maintains cell-specific physiological HNF1A concentrations through positive and negative feedback loops. Pancreatic β cells from Haster mutant mice consequently showed variegated HNF1A silencing or overexpression, resulting in hyperglycaemia. HASTER-dependent negative feedback was essential to prevent HNF1A binding to inappropriate genomic regions. We demonstrate that the HASTER promoter DNA, rather than the lncRNA, modulates HNF1A promoter-enhancer interactions in cis and thereby regulates HNF1A transcription. Our studies expose a cis-regulatory element that is unlike classic enhancers or silencers, it stabilizes the transcription of its target gene and ensures the fidelity of a cell-specific transcription factor program. They also show that disruption of a mammalian lncRNA promoter can cause diabetes mellitus.

Postmitotic differentiation of human monocytes requires cohesin-structured chromatin

Cohesin is a major structural component of mammalian genomes and is required to maintain loop structures. While acute depletion in short-term culture models suggests a limited importance of cohesin for steady-state transcriptional circuits, long-term studies are hampered by essential functions of cohesin during replication. Here, we study genome architecture in a postmitotic differentiation setting, the differentiation of human blood monocytes (MO). We profile and compare epigenetic, transcriptome and 3D conformation landscapes during MO differentiation (either into dendritic cells or macrophages) across the genome and detect numerous architectural changes, ranging from higher level compartments down to chromatin loops. Changes in loop structures correlate with cohesin-binding, as well as epigenetic and transcriptional changes during differentiation. Functional studies show that the siRNA-mediated depletion of cohesin (and to a lesser extent also CTCF) markedly disturbs loop structures and dysregulates genes and enhancers that are primarily regulated during normal MO differentiation. In addition, gene activation programs in cohesin-depleted MO-derived macrophages are disturbed. Our findings implicate an essential function of cohesin in controlling long-term, differentiation- and activation-associated gene expression programs.

How chromatin structure and gene accessibility changes during monocyte differentiation is not clearly defined. Here the authors characterize the chromatin changes during macrophage or dendritic cell maturation from monocytes and the dependence of this upon cohesin and CTCF.

Mammalian PERIOD2 regulates H2A.Z incorporation in chromatin to orchestrate circadian negative feedback

Mammalian circadian oscillators are built on a feedback loop in which the activity of the transcription factor CLOCK-BMAL1 is repressed by the PER-CRY complex. Here, we show that murine Per^-/- fibroblasts display aberrant nucleosome occupancy around transcription start sites (TSSs) and at promoter-proximal and distal CTCF sites due to impaired histone H2A.Z deposition. Knocking out H2A.Z mimicked the Per null chromatin state and disrupted cellular rhythms. We found that endogenous mPER2 complexes retained CTCF as well as the specific H2A.Z-deposition chaperone YL1-a component of the ATP-dependent remodeler SRCAP and p400-TIP60 complex. While depleting YL1 or mutating chaperone-binding sites on H2A.Z lengthened the circadian period, H2A.Z deletion abrogated BMAL1 chromatin recruitment and promoted its proteasomal degradation. We propose that a PER2-mediated H2A.Z deposition pathway (1) compacts CLOCK-BMAL1 binding sites to establish negative feedback, (2) organizes circadian chromatin landscapes using CTCF and (3) bookmarks genomic loci for BMAL1 binding to impinge on the positive arm of the subsequent cycle.

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.

Chromatin Structure and Dynamics: Focus on Neuronal Differentiation and Pathological Implication

Chromatin structure is an essential regulator of gene expression. Its state of compaction contributes to the regulation of genetic programs, in particular during differentiation. Epigenetic processes, which include post-translational modifications of histones, DNA methylation and implication of non-coding RNA, are powerful regulators of gene expression. Neurogenesis and neuronal differentiation are spatio-temporally regulated events that allow the formation of the central nervous system components. Here, we review the chromatin structure and post-translational histone modifications associated with neuronal differentiation. Studying the impact of histone modifications on neuronal differentiation improves our understanding of the pathophysiological mechanisms of chromatinopathies and opens up new therapeutic avenues. In addition, we will discuss techniques for the analysis of histone modifications on a genome-wide scale and the pathologies associated with the dysregulation of the epigenetic machinery.

Three-Dimensional Genome Organization in Breast and Gynecological Cancers: How Chromatin Folding Influences Tumorigenic Transcriptional Programs

A growing body of research on the transcriptome and cancer genome has demonstrated that many gynecological tumor-specific gene mutations are located in cis-regulatory elements. Through chromosomal looping, cis-regulatory elements interact which each other to control gene expression by bringing distant regulatory elements, such as enhancers and insulators, into close proximity with promoters. It is well known that chromatin connections may be disrupted in cancer cells, promoting transcriptional dysregulation and the expression of abnormal tumor suppressor genes and oncogenes. In this review, we examine the roles of alterations in 3D chromatin interactions. This includes changes in CTCF protein function, cancer-risk single nucleotide polymorphisms, viral integration, and hormonal response as part of the mechanisms that lead to the acquisition of enhancers or super-enhancers. The translocation of existing enhancers, as well as enhancer loss or acquisition of insulator elements that interact with gene promoters, is also revised. Remarkably, similar processes that modify 3D chromatin contacts in gene promoters may also influence the expression of non-coding RNAs, such as long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), which have emerged as key regulators of gene expression in a variety of cancers, including gynecological malignancies.

HiC1Dmetrics: framework to extract various one-dimensional features from chromosome structure data

Eukaryotic genomes are organized in a three-dimensional spatial structure. In this regard, the development of chromosome conformation capture methods has enabled studies of chromosome organization on a genomic scale. Hi-C, the high-throughput chromosome conformation capture method, can reveal a population-averaged, hierarchical chromatin structure. The typical Hi-C analysis uses a two-dimensional (2D) contact matrix that indicates contact frequencies between all possible genomic position pairs. Oftentimes, however, such a 2D matrix is not amenable to handling quantitative comparisons, visualizations and integrations across multiple datasets. Although several one-dimensional (1D) metrics have been proposed to depict structural information in Hi-C data, their effectiveness is still underappreciated. Here, we first review the currently available 1D metrics for individual Hi-C samples or two-sample comparisons and then discuss their validity and suitable analysis scenarios. We also propose several new 1D metrics to identify additional unique features of chromosome structures. We highlight that the 1D metrics are reproducible and robust for comparing and visualizing multiple Hi-C samples. Moreover, we show that 1D metrics can be easily combined with epigenome tracks to annotate chromatin states in greater details. We develop a new framework, called HiC1Dmetrics, to summarize all 1D metrics discussed in this study. HiC1Dmetrics is open-source (github.com/wangjk321/HiC1Dmetrics) and can be accessed from both command-line and web-based interfaces. Our tool constitutes a useful resource for the community of chromosome-organization researchers.

Harnessing natural variation to identify cis regulators of sex-biased gene expression in a multi-strain mouse liver model

Sex differences in gene expression are widespread in the liver, where many autosomal factors act in tandem with growth hormone signaling to regulate individual variability of sex differences in liver metabolism and disease. Here, we compare hepatic transcriptomic and epigenetic profiles of mouse strains C57BL/6J and CAST/EiJ, representing two subspecies separated by 0.5-1 million years of evolution, to elucidate the actions of genetic factors regulating liver sex differences. We identify 144 protein coding genes and 78 lncRNAs showing strain-conserved sex bias; many have gene ontologies relevant to liver function, are more highly liver-specific and show greater sex bias, and are more proximally regulated than genes whose sex bias is strain-dependent. The strain-conserved genes include key growth hormone-dependent transcriptional regulators of liver sex bias; however, three other transcription factors, Trim24, Tox, and Zfp809, lose their sex-biased expression in CAST/EiJ mouse liver. To elucidate the observed strain specificities in expression, we characterized the strain-dependence of sex-biased chromatin opening and enhancer marks at cis regulatory elements (CREs) within expression quantitative trait loci (eQTL) regulating liver sex-biased genes. Strikingly, 208 of 286 eQTLs with strain-specific, sex-differential effects on expression were associated with a complete gain, loss, or reversal of the sex differences in expression between strains. Moreover, 166 of the 286 eQTLs were linked to the strain-dependent gain or loss of localized sex-biased CREs. Remarkably, a subset of these CREs apparently lacked strain-specific genetic variants yet showed coordinated, strain-dependent sex-biased epigenetic regulation. Thus, we directly link hundreds of strain-specific genetic variants to the high variability in CRE activity and expression of sex-biased genes and uncover underlying genetically-determined epigenetic states controlling liver sex bias in genetically diverse mouse populations.

Cohesin: behind dynamic genome topology and gene expression reprogramming

Beyond its originally discovered role tethering replicated sister chromatids, cohesin has emerged as a master regulator of gene expression. Recent advances in chromatin topology resolution and single-cell studies have revealed that cohesin has a pivotal role regulating highly dynamic chromatin interactions linked to transcription control. The dynamic association of cohesin with chromatin and its capacity to perform loop extrusion contribute to the heterogeneity of chromatin contacts. Additionally, different cohesin subcomplexes, with specific properties and regulation, control gene expression across the cell cycle and during developmental cell commitment. Here, we discuss the most recent literature in the field to highlight the role of cohesin in gene expression regulation during transcriptional shifts and its relationship with human diseases.

Chromatin insulation dynamics in glioblastoma: challenges and future perspectives of precision oncology

Glioblastoma (GBM) is the most aggressive primary brain tumor, having a poor prognosis and a median overall survival of less than two years. Over the last decade, numerous findings regarding the distinct molecular and genetic profiles of GBM have led to the emergence of several therapeutic approaches. Unfortunately, none of them has proven to be effective against GBM progression and recurrence. Epigenetic mechanisms underlying GBM tumor biology, including histone modifications, DNA methylation, and chromatin architecture, have become an attractive target for novel drug discovery strategies. Alterations on chromatin insulator elements (IEs) might lead to aberrant chromatin remodeling via DNA loop formation, causing oncogene reactivation in several types of cancer, including GBM. Importantly, it is shown that mutations affecting the isocitrate dehydrogenase (IDH) 1 and 2 genes, one of the most frequent genetic alterations in gliomas, lead to genome-wide DNA hypermethylation and the consequent IE dysfunction. The relevance of IEs has also been observed in a small population of cancer stem cells known as glioma stem cells (GSCs), which are thought to participate in GBM tumor initiation and drug resistance. Recent studies revealed that epigenomic alterations, specifically chromatin insulation and DNA loop formation, play a crucial role in establishing and maintaining the GSC transcriptional program. This review focuses on the relevance of IEs in GBM biology and their implementation as a potential theranostic target to stratify GBM patients and develop novel therapeutic approaches. We will also discuss the state-of-the-art emerging technologies using big data analysis and how they will settle the bases on future diagnosis and treatment strategies in GBM patients.

Disruption of NIPBL/Scc2 in Cornelia de Lange Syndrome provokes cohesin genome-wide redistribution with an impact in the transcriptome

Cornelia de Lange syndrome (CdLS) is a rare disease affecting multiple organs and systems during development. Mutations in the cohesin loader, NIPBL/Scc2, were first described and are the most frequent in clinically diagnosed CdLS patients. The molecular mechanisms driving CdLS phenotypes are not understood. In addition to its canonical role in sister chromatid cohesion, cohesin is implicated in the spatial organization of the genome. Here, we investigate the transcriptome of CdLS patient-derived primary fibroblasts and observe the downregulation of genes involved in development and system skeletal organization, providing a link to the developmental alterations and limb abnormalities characteristic of CdLS patients. Genome-wide distribution studies demonstrate a global reduction of NIPBL at the NIPBL-associated high GC content regions in CdLS-derived cells. In addition, cohesin accumulates at NIPBL-occupied sites at CpG islands potentially due to reduced cohesin translocation along chromosomes, and fewer cohesin peaks colocalize with CTCF.

Neonatal thyroxine activation modifies epigenetic programming of the liver

The type 2 deiodinase (D2) in the neonatal liver accelerates local thyroid hormone triiodothyronine (T3) production and expression of T3-responsive genes. Here we show that this surge in T3 permanently modifies hepatic gene expression. Liver-specific Dio2 inactivation (Alb-D2KO) transiently increases H3K9me3 levels during post-natal days 1-5 (P1-P5), and results in methylation of 1,508 DNA sites (H-sites) in the adult mouse liver. These sites are associated with 1,551 areas of reduced chromatin accessibility (RCA) within core promoters and 2,426 within intergenic regions, with reduction in the expression of 1,363 genes. There is strong spatial correlation between density of H-sites and RCA sites. Chromosome conformation capture (Hi-C) data reveals a set of 81 repressed genes with a promoter RCA in contact with an intergenic RCA ~300 Kbp apart, within the same topologically associating domain (χ2 = 777; p < 0.00001). These data explain how the systemic hormone T3 acts locally during development to define future expression of hepatic genes.

Predicting pathogenic non-coding SVs disrupting the 3D genome in 1646 whole cancer genomes using multiple instance learning

Over the past years, large consortia have been established to fuel the sequencing of whole genomes of many cancer patients. Despite the increased abundance in tools to study the impact of SNVs, non-coding SVs have been largely ignored in these data. Here, we introduce svMIL2, an improved version of our Multiple Instance Learning-based method to study the effect of somatic non-coding SVs disrupting boundaries of TADs and CTCF loops in 1646 cancer genomes. We demonstrate that svMIL2 predicts pathogenic non-coding SVs with an average AUC of 0.86 across 12 cancer types, and identifies non-coding SVs affecting well-known driver genes. The disruption of active (super) enhancers in open chromatin regions appears to be a common mechanism by which non-coding SVs exert their pathogenicity. Finally, our results reveal that the contribution of pathogenic non-coding SVs as opposed to driver SNVs may highly vary between cancers, with notably high numbers of genes being disrupted by pathogenic non-coding SVs in ovarian and pancreatic cancer. Taken together, our machine learning method offers a potent way to prioritize putatively pathogenic non-coding SVs and leverage non-coding SVs to identify driver genes. Moreover, our analysis of 1646 cancer genomes demonstrates the importance of including non-coding SVs in cancer diagnostics.

Global analysis of expression, maturation and subcellular localization of mouse liver transcriptome identifies novel sex-biased and TCPOBOP-responsive long non-coding RNAs

BACKGROUND

While nuclear transcription and RNA processing and localization are well established for protein coding genes (PCGs), these processes are poorly understood for long non-coding (lnc)RNAs. Here, we characterize global patterns of transcript expression, maturation and localization for mouse liver RNA, including more than 15,000 lncRNAs. PolyA-selected liver RNA was isolated and sequenced from four subcellular fractions (chromatin, nucleoplasm, total nucleus, and cytoplasm), and from the chromatin-bound fraction without polyA selection.

RESULTS

Transcript processing, determined from normalized intronic to exonic sequence read density ratios, progressively increased for PCG transcripts in going from the chromatin-bound fraction to the nucleoplasm and then on to the cytoplasm. Transcript maturation was similar for lncRNAs in the chromatin fraction, but was significantly lower in the nucleoplasm and cytoplasm. LncRNA transcripts were 11-fold more likely to be significantly enriched in the nucleus than cytoplasm, and 100-fold more likely to be significantly chromatin-bound than nucleoplasmic. Sequencing chromatin-bound RNA greatly increased the sensitivity for detecting lowly expressed lncRNAs and enabled us to discover and localize hundreds of novel regulated liver lncRNAs, including lncRNAs showing sex-biased expression or responsiveness to TCPOBOP a xenobiotic agonist ligand of constitutive androstane receptor (Nr1i3).

CONCLUSIONS

Integration of our findings with prior studies and lncRNA annotations identified candidate regulatory lncRNAs for a variety of hepatic functions based on gene co-localization within topologically associating domains or transcription divergent or antisense to PCGs associated with pathways linked to hepatic physiology and disease.

The Dm-Myb Oncoprotein Contributes to Insulator Function and Stabilizes Repressive H3K27me3 PcG Domains

Drosophila Myb (Dm-Myb) encodes a protein that plays a key role in regulation of mitotic phase genes. Here, we further refine its role in the context of a developing tissue as a potentiator of gene expression required for proper RNA polymerase II (RNA Pol II) function and efficient H3K4 methylation at promoters. In contrast to its role in gene activation, Myb is also required for repression of many genes, although no specific mechanism for this role has been proposed. We now reveal a critical role for Myb in contributing to insulator function, in part by promoting binding of insulator proteins BEAF-32 and CP190 and stabilizing H3K27me3 Polycomb-group (PcG) domains. In the absence of Myb, H3K27me3 is markedly reduced throughout the genome, leading to H3K4me3 spreading and gene derepression. Finally, Myb is enriched at boundaries that demarcate chromatin environments, including chromatin loop anchors. These results reveal functions of Myb that extend beyond transcriptional regulation.

Myb has been considered a transcriptional activator of primarily M phase genes. Here, Santana et al. show that Myb also contributes to insulator function, in part by promoting binding of insulator factors, and is required to stabilize repressive domains in the genome.

STAT5 Regulation of Sex-Dependent Hepatic CpG Methylation at Distal Regulatory Elements Mapping to Sex-Biased Genes

Growth hormone-activated STAT5b is an essential regulator of sex-differential gene expression in mouse liver; however, its impact on hepatic gene expression and epigenetic responses is poorly understood. Here, we found a substantial, albeit incomplete loss of liver sex bias in hepatocyte-specific STAT5a/STAT5b (collectively, STAT5)-deficient mouse liver. In male liver, many male-biased genes were downregulated in direct association with the loss of STAT5 binding; many female-biased genes, which show low STAT5 binding, were derepressed, indicating an indirect mechanism for repression by STAT5. Extensive changes in CpG methylation were seen in STAT5-deficient liver, where sex differences were abolished at 88% of ∼1,500 sex-differentially methylated regions, largely due to increased DNA methylation upon STAT5 loss. STAT5-dependent CpG hypomethylation was rarely found at proximal promoters of STAT5-dependent genes. Rather, STAT5 primarily regulated the methylation of distal enhancers, where STAT5 deficiency induced widespread hypermethylation at genomic regions enriched for accessible chromatin, enhancer histone marks (histone H3 lysine 4 monomethylation [H3K4me1] and histone H3 lysine 27 acetylation [H3K27ac]), STAT5 binding, and DNA motifs for STAT5 and other transcription factors implicated in liver sex differences. Thus, the sex-dependent binding of STAT5 to liver chromatin is closely linked to the sex-dependent demethylation of distal regulatory elements linked to STAT5-dependent genes important for liver sex bias.

Genome-wide prediction of topoisomerase IIβ binding by architectural factors and chromatin accessibility

DNA topoisomerase II-β (TOP2B) is fundamental to remove topological problems linked to DNA metabolism and 3D chromatin architecture, but its cut-and-reseal catalytic mechanism can accidentally cause DNA double-strand breaks (DSBs) that can seriously compromise genome integrity. Understanding the factors that determine the genome-wide distribution of TOP2B is therefore not only essential for a complete knowledge of genome dynamics and organization, but also for the implications of TOP2-induced DSBs in the origin of oncogenic translocations and other types of chromosomal rearrangements. Here, we conduct a machine-learning approach for the prediction of TOP2B binding using publicly available sequencing data. We achieve highly accurate predictions, with accessible chromatin and architectural factors being the most informative features. Strikingly, TOP2B is sufficiently explained by only three features: DNase I hypersensitivity, CTCF and cohesin binding, for which genome-wide data are widely available. Based on this, we develop a predictive model for TOP2B genome-wide binding that can be used across cell lines and species, and generate virtual probability tracks that accurately mirror experimental ChIP-seq data. Our results deepen our knowledge on how the accessibility and 3D organization of chromatin determine TOP2B function, and constitute a proof of principle regarding the in silico prediction of sequence-independent chromatin-binding factors.

Hexavalent chromium promotes differential binding of CTCF to its cognate sites in Euchromatin

Hexavalent chromium compounds are well-established respiratory carcinogens to which humans are commonly exposed in industrial and occupational settings. In addition, natural and anthropogenic sources of these compounds contribute to the exposure of global populations through multiple routes, including dermal, ingestion and inhalation that elevate the risk of cancer by largely unresolved mechanisms. Cr(VI) has genotoxic properties that include ternary adduct formation with DNA, increases in DNA damage, mostly by double-strand break formation, and altered transcriptional responses. Our previous work using ATAC-seq showed that CTCF motifs were enriched in Cr(VI)-dependent differentially accessible chromatin, suggesting that CTCF, a key transcription factor responsible for the regulation of the transcriptome, might be a chromium target. To test this hypothesis, we investigated the effect of Cr(VI) treatment on the binding of CTCF to its cognate sites and ensuing changes in transcription-related histone modifications. Differentially bound CTCF sites were enriched by Cr(VI) treatment within transcription-related regions, specifically transcription start sites and upstream genic regions. Functional annotation of the affected genes highlighted biological processes previously associated with Cr(VI) exposure. Notably, we found that differentially bound CTCF sites proximal to the promoters of this subset of genes were frequently associated with the active histone marks H3K27ac, H3K4me3, and H3K36me3, in agreement with the concept that Cr(VI) targets CTCF in euchromatic regions of the genome. Our results support the conclusion that Cr(VI) treatment promotes the differential binding of CTCF to its cognate sites in genes near transcription-active boundaries, targeting these genes for dysregulation.

DeepYY1: a deep learning approach to identify YY1-mediated chromatin loops

The protein Yin Yang 1 (YY1) could form dimers that facilitate the interaction between active enhancers and promoter-proximal elements. YY1-mediated enhancer-promoter interaction is the general feature of mammalian gene control. Recently, some computational methods have been developed to characterize the interactions between DNA elements by elucidating important features of chromatin folding; however, no computational methods have been developed for identifying the YY1-mediated chromatin loops. In this study, we developed a deep learning algorithm named DeepYY1 based on word2vec to determine whether a pair of YY1 motifs would form a loop. The proposed models showed a high prediction performance (AUCs$\ge$0.93) on both training datasets and testing datasets in different cell types, demonstrating that DeepYY1 has an excellent performance in the identification of the YY1-mediated chromatin loops. Our study also suggested that sequences play an important role in the formation of YY1-mediated chromatin loops. Furthermore, we briefly discussed the distribution of the replication origin site in the loops. Finally, a user-friendly web server was established, and it can be freely accessed at http://lin-group.cn/server/DeepYY1.

Genetic factors contributing to extensive variability of sex-specific hepatic gene expression in Diversity Outbred mice

Sex-specific transcription characterizes hundreds of genes in mouse liver, many implicated in sex-differential drug and lipid metabolism and disease susceptibility. While the regulation of liver sex differences by growth hormone-activated STAT5 is well established, little is known about autosomal genetic factors regulating the sex-specific liver transcriptome. Here we show, using genotyping and expression data from a large population of Diversity Outbred mice, that genetic factors work in tandem with growth hormone to control the individual variability of hundreds of sex-biased genes, including many long non-coding RNA genes. Significant associations between single nucleotide polymorphisms and sex-specific gene expression were identified as expression quantitative trait loci (eQTLs), many of which showed strong sex-dependent associations. Remarkably, autosomal genetic modifiers of sex-specific genes were found to account for more than 200 instances of gain or loss of sex-specificity across eight Diversity Outbred mouse founder strains. Sex-biased STAT5 binding sites and open chromatin regions with strain-specific variants were significantly enriched at eQTL regions regulating correspondingly sex-specific genes, supporting the proposed functional regulatory nature of the eQTL regions identified. Binding of the male-biased, growth hormone-regulated repressor BCL6 was most highly enriched at trans-eQTL regions controlling female-specific genes. Co-regulated gene clusters defined by overlapping eQTLs included sets of highly correlated genes from different chromosomes, further supporting trans-eQTL action. These findings elucidate how an unexpectedly large number of autosomal factors work in tandem with growth hormone signaling pathways to regulate the individual variability associated with sex differences in liver metabolism and disease.

Topologically Associated Domains Delineate Susceptibility to Somatic Hypermutation

Somatic hypermutation (SHM) introduces point mutations into immunoglobulin (Ig) genes but also causes mutations in other parts of the genome. We have used lentiviral SHM reporter vectors to identify regions of the genome that are susceptible (“hot”) and resistant (“cold”) to SHM, revealing that SHM susceptibility and resistance are often properties of entire topologically associated domains (TADs). Comparison of hot and cold TADs reveals that while levels of transcription are equivalent, hot TADs are enriched for the cohesin loader NIPBL, super-enhancers, markers of paused/stalled RNA polymerase 2, and multiple important B cell transcription factors. We demonstrate that at least some hot TADs contain enhancers that possess SHM targeting activity and that insertion of a strong Ig SHM-targeting element into a cold TAD renders it hot. Our findings lead to a model for SHM susceptibility involving the cooperative action of cis-acting SHM targeting elements and the dynamic and architectural properties of TADs.

Senigl et al. show that genome susceptibility to somatic hypermutation (SHM) is confined within topologically associated domains (TADs) and is linked to markers of strong enhancers and stalled transcription and high levels of the cohesin loader NIPBL. Insertion of an ectopic SHM targeting element renders an entire TAD susceptible to SHM.

GATA-1-dependent histone H3K27 acetylation mediates erythroid cell-specific chromatin interaction between CTCF sites

CCCTC-binding factor (CTCF) sites interact with each other in the chromatin environment, establishing chromatin domains. Our previous study showed that interaction between CTCF sites is cell type-specific around the β-globin locus and is dependent on erythroid-specific activator GATA-1. To find out molecular mechanisms of the cell type-specific interaction, we directly inhibited GATA-1 binding to the β-globin enhancers by deleting its binding motifs and found that histone H3K27 acetylation (H3K27ac) was decreased at CTCF sites surrounding the β-globin locus, even though CTCF binding itself was maintained at the sites. Forced H3K27ac by Trichostatin A treatment or CBP/p300 KD affected the interactions between CTCF sites around the β-globin locus without changes in CTCF binding. Analysis of public ChIA-PET data revealed that H3K27ac is higher at CTCF sites forming short interactions than long interactions. GATA-1 was identified as a representative transcription factor that relates with genes present inside the short interactions in erythroid K562 cells. Depletion of GATA-1-reduced H3K27ac at CTCF sites near erythroid-specific enhancers. These results indicate that H3K27ac at CTCF sites is required for cell type-specific chromatin interactions between them. Tissue-specific activator GATA-1 appears to play a role in H3K27ac at CTCF sites in erythroid cells.

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.

In silico prediction of high-resolution Hi-C interaction matrices

The three-dimensional (3D) organization of the genome plays an important role in gene regulation bringing distal sequence elements in 3D proximity to genes hundreds of kilobases away. Hi-C is a powerful genome-wide technique to study 3D genome organization. Owing to experimental costs, high resolution Hi-C datasets are limited to a few cell lines. Computational prediction of Hi-C counts can offer a scalable and inexpensive approach to examine 3D genome organization across multiple cellular contexts. Here we present HiC-Reg, an approach to predict contact counts from one-dimensional regulatory signals. HiC-Reg predictions identify topologically associating domains and significant interactions that are enriched for CCCTC-binding factor (CTCF) bidirectional motifs and interactions identified from complementary sources. CTCF and chromatin marks, especially repressive and elongation marks, are most important for HiC-Reg’s predictive performance. Taken together, HiC-Reg provides a powerful framework to generate high-resolution profiles of contact counts that can be used to study individual locus level interactions and higher-order organizational units of the genome.

Existing computational approaches to predict long-range regulatory interactions do not fully exploit high-resolution Hi-C datasets. Here the authors present a Random Forests regression-based approach to predict high-resolution Hi-C counts using one-dimensional regulatory genomic signals.

Tales from topographic oceans: topologically associated domains and cancer

The 3D organization of the genome within the cell nucleus has come into sharp focus over the last decade. This has largely arisen because of the application of genomic approaches that have revealed numerous levels of genomic and chromatin interactions, including topologically associated domains (TADs). The current review examines how these domains were identified, are organized, how their boundaries arise and are regulated, and how genes within TADs are coordinately regulated. There are many examples of the disruption to TAD structure in cancer and the altered regulation, structure and function of TADs are discussed in the context of hormone responsive cancers, including breast, prostate and ovarian cancer. Finally, some aspects of the statistical insight and computational skills required to interrogate TAD organization are considered and future directions discussed.

DNA sequence-dependent chromatin architecture and nuclear hubs formation

In this study, by exploring chromatin conformation capture data, we show that the nuclear segregation of Topologically Associated Domains (TADs) is contributed by DNA sequence composition. GC-peaks and valleys of TADs strongly influence interchromosomal interactions and chromatin 3D structure. To gain insight on the compositional and functional constraints associated with chromatin interactions and TADs formation, we analysed intra-TAD and intra-loop GC variations. This led to the identification of clear GC-gradients, along which, the density of genes, super-enhancers, transcriptional activity, and CTCF binding sites occupancy co-vary non-randomly. Further, the analysis of DNA base composition of nucleolar aggregates and nuclear speckles showed strong sequence-dependant effects. We conjecture that dynamic DNA binding affinity and flexibility underlay the emergence of chromatin condensates, their growth is likely promoted in mechanically soft regions (GC-rich) of the lowest chromatin and nucleosome densities. As a practical perspective, the strong linear association between sequence composition and interchromosomal contacts can help define consensus chromatin interactions, which in turn may be used to study alternative states of chromatin architecture.

Comparison of computational methods for 3D genome analysis at single-cell Hi-C level

Hi-C is a high-throughput chromosome conformation capture technology that is becoming routine in the literature. Although the price of sequencing has been dropping dramatically, high-resolution Hi-C data are not always an option for many studies, such as in single cells. However, the performance of current computational methods based on Hi-C at the ultra-sparse data condition has yet to be fully assessed. Therefore, in this paper, after briefly surveying the primary computational methods for Hi-C data analysis, we assess the performance of representative methods on data normalization, identification of compartments, Topologically Associating Domains (TADs) and chromatin loops under the condition of ultra-low resolution. We showed that most state-of-the-art methods do not work properly for that condition. Then, we applied the three best-performing methods on real single-cell Hi-C data, and their performance indicates that compartments may be a statistical feature emerging from the cell population, while TADs and chromatin loops may dynamically exist in single cells.

Widespread epigenetic changes to the enhancer landscape of mouse liver induced by a specific xenobiotic agonist ligand of the nuclear receptor CAR

CAR (Nr1i3), a liver nuclear receptor and xenobiotic sensor, induces drug, steroid and lipid metabolism and dysregulates genes linked to hepatocellular carcinogenesis, but its impact on the liver epigenome is poorly understood. TCPOBOP, a halogenated xenochemical and highly specific CAR agonist ligand, induces localized chromatin opening or closing at several thousand mouse liver genomic regions, discovered as differential DNase-hypersensitive sites (ΔDHS). Active enhancer and promoter histone marks induced by TCPOBOP were enriched at opening DHS and TCPOBOP-inducible genes. Enrichment of CAR binding and CAR motifs was seen at opening DHS and their inducible drug/lipid metabolism gene targets, and at many constitutively open DHS located nearby. TCPOBOP-responsive cell cycle and DNA replication genes co-dependent on MET/EGFR signaling for induction were also enriched for CAR binding. A subset of opening DHS and many closing DHS mapping to TCPOBOP-responsive target genes did not bind CAR, indicating an indirect mechanism for their changes in chromatin accessibility. TCPOBOP-responsive DHS were also enriched for induced binding of RXRA, CEBPA and CEBPB, and for motifs for liver-enriched factors that may contribute to liver-specific transcriptional responses to TCPOBOP exposure. These studies elucidate the enhancer landscape of TCPOBOP-exposed liver and the widespread epigenetic changes that are induced by both direct and indirect mechanisms linked to CAR activation. The global maps of thousands of environmental chemical-induced epigenetic changes described here constitute a rich resource for further research on xenochemical effects on liver chromatin states and the epigenome.

Sex-Biased lncRNAs Inversely Correlate With Sex-Opposite Gene Coexpression Networks in Diversity Outbred Mouse Liver

Sex differences in liver gene expression are determined by pituitary growth hormone secretion patterns, which regulate sex-dependent liver transcription factors and establish sex-specific chromatin states. Hypophysectomy (hypox) identifies two major classes of liver sex-biased genes, defined by their sex-dependent positive or negative responses to pituitary hormone ablation. However, the mechanisms that underlie each hypox-response class are unknown. We sought to discover candidate, regulatory, long noncoding RNAs (lncRNAs) controlling responsiveness to hypox. We characterized gene structures and expression patterns for 15,558 mouse liver-expressed lncRNAs, including many sex-specific lncRNAs regulated during postnatal development or subject to circadian regulation. Using the high natural allelic variance of Diversity Outbred (DO) mice, we discovered tightly coexpressed clusters of sex-specific protein-coding genes (gene modules) in male and female DO liver. Remarkably, many gene modules were strongly enriched for sex-specific genes within a single hypox-response class, indicating that the genetic heterogeneity of DO mice encompasses responsiveness to hypox. Moreover, several distant gene modules were enriched for gene subsets of the same hypox-response class, highlighting the complex regulation of hypox-responsiveness. Finally, we identified eight sex-specific lncRNAs with strong negative regulatory potential, as indicated by their strong negative correlation of expression across DO mouse livers with that of protein-coding gene modules enriched for genes of the opposite sex bias and inverse hypox-response class. These findings reveal an important role for genetic factors in regulating responsiveness to hypox, and present testable hypotheses for the roles of sex-biased liver lncRNAs in controlling the sex-bias of liver gene expression.

DNA (de)methylation in embryonic stem cells controls CTCF-dependent chromatin boundaries

Coordinated changes of DNA (de)methylation, nucleosome positioning, and chromatin binding of the architectural protein CTCF play an important role for establishing cell-type-specific chromatin states during differentiation. To elucidate molecular mechanisms that link these processes, we studied the perturbed DNA modification landscape in mouse embryonic stem cells (ESCs) carrying a double knockout (DKO) of the Tet1 and Tet2 dioxygenases. These enzymes are responsible for the conversion of 5-methylcytosine (5mC) into its hydroxymethylated (5hmC), formylated (5fC), or carboxylated (5caC) forms. We determined changes in nucleosome positioning, CTCF binding, DNA methylation, and gene expression in DKO ESCs and developed biophysical models to predict differential CTCF binding. Methylation-sensitive nucleosome repositioning accounted for a significant portion of CTCF binding loss in DKO ESCs, whereas unmethylated and nucleosome-depleted CpG islands were enriched for CTCF sites that remained occupied. A number of CTCF sites also displayed direct correlations with the CpG modification state: CTCF was preferentially lost from sites that were marked with 5hmC in wild-type (WT) cells but not from 5fC-enriched sites. In addition, we found that some CTCF sites can act as bifurcation points defining the differential methylation landscape. CTCF loss from such sites, for example, at promoters, boundaries of chromatin loops, and topologically associated domains (TADs), was correlated with DNA methylation/demethylation spreading and can be linked to down-regulation of neighboring genes. Our results reveal a hierarchical interplay between cytosine modifications, nucleosome positions, and DNA sequence that determines differential CTCF binding and regulates gene expression.

LRH1-driven transcription factor circuitry for hepatocyte identity: Super-enhancer cistromic analysis

BACKGROUND

The injured liver loses normal function, with concomitant decrease of key identity genes. Super-enhancers contribute to mammalian cell identity. Here, we identified core transcription factors (TFs) that are active in hepatocytes, using genome-wide analysis and hierarchical ordering of super-enhancer distribution.

METHODS

Expression of core TFs was assessed in a cohort of patients with hepatitis or cirrhosis and animal models. Quantitative PCR, chromatin immunoprecipitation assays, and hydrodynamic gene delivery methods were used to assess gene regulation and hepatocyte viability. RNA-sequencing data were generated to investigate the role of LRH1 in hepatocyte protection from injury.

RESULTS

Network analysis of super-enhancer-associated gene interactions and expression arrays for cohorts of patients with hepatitis and cirrhosis enabled us to identify a super-enhancer-associated network, and LRH1, HNF4α, PPARα, and RXRα as core TFs. In mouse models, expression of core TFs was robustly inhibited by single and multiple challenge(s) with liver toxicant. RNA-seq analysis revealed changes in expression in the super-enhancer-associated genes sensitively biased toward repression by intoxication. LRH1 gene delivery prevented the loss of hepatic super-enhancer-associated signaling circuitry in toxicant-challenged mice, and protected the liver from injury, indicating the role of LRH1 in hepatocyte identity and viability. In hepatocytes, overexpression of each core TF promoted induction of other TFs.

CONCLUSION

Overall, this study identified LRH1-driven pathway as a circuitry responsible for hepatocyte identity by using cistromic analysis, improving our understanding of liver pathophysiology and identifying novel therapeutic targets.

The structural and functional roles of CTCF in the regulation of cell type-specific and human disease-associated super-enhancers

BACKGROUND

Super-enhancers play critical roles in cell-type specific gene controls and human disease progression. CCCTC-binding factor (CTCF), a transcriptional repressor that insulates the expression of neighboring genes and is involved in chromatin interactions, is frequently present in the boundary regions of or within super-enhancers. However, the structural and functional roles of CTCF in regulating super-enhancers remain elusive.

OBJECTIVE

To provide a comprehensive review describing the distinct chromatin features and functional roles of CTCF within super-enhancers.

METHODS

This review compares the various tools used to study the three-dimensional (3D) chromatin architecture of super-enhancers; summarizes the chromatin features of CTCF within cell-type specific super-enhancers and their in vivo biological activities, as determined by CRISPR/Cas9 genome editing; and describes the structural and functional activities of CTCF within human disease-associated super-enhancers.

CONCLUSION

This review provides fundamental insights into the regulatory mechanisms of super-enhancers and facilitates studies of tissue-specific developmental processes and human disease progression.