TADs as modular and dynamic units for gene regulation by hormones

Implications of the three-dimensional chromatin organization for genome evolution in a fungal plant pathogen

The spatial organization of eukaryotic genomes is linked to their biological functions, although it is not clear how this impacts the overall evolution of a genome. Here, we uncover the three-dimensional (3D) genome organization of the phytopathogen Verticillium dahliae, known to possess distinct genomic regions, designated adaptive genomic regions (AGRs), enriched in transposable elements and genes that mediate host infection. Short-range DNA interactions form clear topologically associating domains (TADs) with gene-rich boundaries that show reduced levels of gene expression and reduced genomic variation. Intriguingly, TADs are less clearly insulated in AGRs than in the core genome. At a global scale, the genome contains bipartite long-range interactions, particularly enriched for AGRs and more generally containing segmental duplications. Notably, the patterns observed for V. dahliae are also present in other Verticillium species. Thus, our analysis links 3D genome organization to evolutionary features conserved throughout the Verticillium genus.

Long-range gene regulation in hormone-dependent cancer

The human genome is organized into multiple structural layers, ranging from chromosome territories to progressively smaller substructures, such as topologically associating domains (TADs) and chromatin loops. These substructures, collectively referred to as long-range chromatin interactions (LRIs), have a significant role in regulating gene expression. TADs are regions of the genome that harbour groups of genes and regulatory elements that frequently interact with each other and are insulated from other regions, thereby preventing widespread uncontrolled DNA contacts. Chromatin loops formed within TADs through enhancer and promoter interactions are elastic, allowing transcriptional heterogeneity and stochasticity. Over the past decade, it has become evident that the 3D genome structure, also referred to as the chromatin architecture, is central to many transcriptional cellular decisions. In this Review, we delve into the intricate relationship between steroid receptors and LRIs, discussing how steroid receptors interact with and modulate these chromatin interactions. Genetic alterations in the many processes involved in organizing the nuclear architecture are often associated with the development of hormone-dependent cancers. A better understanding of the interplay between architectural proteins and hormone regulatory networks can ultimately be exploited to develop improved approaches for cancer treatment.

Domain Model of Eukaryotic Genome Organization: From DNA Loops Fixed on the Nuclear Matrix to TADs

The article reviews the development of ideas on the domain organization of eukaryotic genome, with special attention on the studies of DNA loops anchored to the nuclear matrix and their role in the emergence of the modern model of eukaryotic genome spatial organization. Critical analysis of results demonstrating that topologically associated chromatin domains are structural-functional blocks of the genome supports the notion that these blocks are fundamentally different from domains whose existence was proposed by the domain hypothesis of eukaryotic genome organization formulated in the 1980s. Based on the discussed evidence, it is concluded that the model postulating that eukaryotic genome is built from uniformly organized structural-functional blocks has proven to be untenable.

The Clock Takes Shape-24 h Dynamics in Genome Topology

Circadian rhythms orchestrate organismal physiology and behavior in order to anticipate daily changes in the environment. Virtually all cells have an internal rhythm that is synchronized every day by Zeitgebers (environmental cues). The synchrony between clocks within the animal enables the fitness and the health of organisms. Conversely, disruption of rhythms is linked to a variety of disorders: aging, cancer, metabolic diseases, and psychological disorders among others. At the cellular level, mammalian circadian rhythms are built on several layers of complexity. The transcriptional-translational feedback loop (TTFL) was the first to be described in the 90s. Thereafter oscillations in epigenetic marks highlighted the role of chromatin state in organizing the TTFL. More recently, studies on the 3D organization of the genome suggest that genome topology could be yet another layer of control on cellular circadian rhythms. The dynamic nature of genome topology over a solar day implies that the 3D mammalian genome has to be considered in the fourth dimension-in time. Whether oscillations in genome topology are a consequence of 24 h gene-expression or a driver of transcriptional cycles remains an open question. All said and done, circadian clock-gated phenomena such as gene expression, DNA damage response, cell metabolism and animal behavior-go hand in hand with 24 h rhythms in genome topology.

Co-Regulated Genes and Gene Clusters

There are many co-regulated genes in eukaryotic cells. The coordinated activation or repression of such genes occurs at specific stages of differentiation, or under the influence of external stimuli. As a rule, co-regulated genes are dispersed in the genome. However, there are also gene clusters, which contain paralogous genes that encode proteins with similar functions. In this aspect, they differ significantly from bacterial operons containing functionally linked genes that are not paralogs. In this review, we discuss the reasons for the existence of gene clusters in vertebrate cells and propose that clustering is necessary to ensure the possibility of selective activation of one of several similar genes.

The Role of Epigenomic Regulatory Pathways in the Gut-Brain Axis and Visceral Hyperalgesia

The gut-brain axis (GBA) is broadly accepted to describe the bidirectional circuit that links the gastrointestinal tract with the central nervous system (CNS). Interest in the GBA has grown dramatically over past two decades along with advances in our understanding of the importance of the axis in the pathophysiology of numerous common clinical disorders including mood disorders, neurodegenerative disease, diabetes mellitus, non-alcohol fatty liver disease (NAFLD) and enhanced abdominal pain (visceral hyperalgesia). Paralleling the growing interest in the GBA, there have been seminal developments in our understanding of how environmental factors such as psychological stress and other extrinsic factors alter gene expression, primarily via epigenomic regulatory mechanisms. This process has been driven by advances in next-generation multi-omics methods and bioinformatics. Recent reviews address various components of GBA, but the role of epigenomic regulatory pathways in chronic stress-associated visceral hyperalgesia in relevant regions of the GBA including the amygdala, spinal cord, primary afferent (nociceptive) neurons, and the intestinal barrier has not been addressed. Rapidly developing evidence suggests that intestinal epithelial barrier dysfunction and microbial dysbiosis play a potentially significant role in chronic stress-associated visceral hyperalgesia in nociceptive neurons innervating the lower intestine via downregulation in intestinal epithelial cell tight junction protein expression and increase in paracellular permeability. These observations support an important role for the regulatory epigenome in the development of future diagnostics and therapeutic interventions in clinical disorders affecting the GBA.

TAD boundary and strength prediction by integrating sequence and epigenetic profile information

Topologically associated domains (TADs) are one of the important higher order chromatin structures with various sizes in the eukaryotic genomes. TAD boundaries, as the flanking regions between adjacent domains, can restrict the interactions of regulatory elements, including enhancers and promoters, and are generally dynamic and variable in different cells. However, the influence of sequence and epigenetic profile-based features in the identification of TAD boundaries is largely unknown. In this work, we proposed a method called pTADS (prediction of TAD boundary and strength), to predict TAD boundaries and boundary strength across multiple cell lines with DNA sequence and epigenetic profile information. The performance was assessed in seven cell lines and three TAD calling methods. The results demonstrate that the TAD boundary can be well predicted by the selected shared features across multiple cell lines. Especially, the model can be transferable to predict the TAD boundary from one cell line to other cell lines. The boundary strength can be characterized by boundary score with good performance. The predicted TAD boundary and TAD boundary strength are further confirmed by three Hi-C contact matrix-based methods across multiple cell lines. The codes and datasets are available at https://github.com/chrom3DEpi/pTADS.

Independence of chromatin conformation and gene regulation during Drosophila dorsoventral patterning

The relationship between chromatin organization and gene regulation remains unclear. While disruption of chromatin domains and domain boundaries can lead to misexpression of developmental genes, acute depletion of regulators of genome organization has a relatively small effect on gene expression. It is therefore uncertain whether gene expression and chromatin state drive chromatin organization or whether changes in chromatin organization facilitate cell-type-specific activation of gene expression. Here, using the dorsoventral patterning of the Drosophila melanogaster embryo as a model system, we provide evidence for the independence of chromatin organization and dorsoventral gene expression. We define tissue-specific enhancers and link them to expression patterns using single-cell RNA-seq. Surprisingly, despite tissue-specific chromatin states and gene expression, chromatin organization is largely maintained across tissues. Our results indicate that tissue-specific chromatin conformation is not necessary for tissue-specific gene expression but rather acts as a scaffold facilitating gene expression when enhancers become active.

Impact of CAR Agonist Ligand TCPOBOP on Mouse Liver Chromatin Accessibility

Activation of the nuclear receptor and transcription factor CAR (Nr1i3) by its specific agonist ligand TCPOBOP (1, 4-bis[2-(3, 5-dichloropyridyloxy)]benzene) dysregulates hundreds of genes in mouse liver and is linked to male-biased hepatocarcinogenesis. To elucidate the genomic organization of CAR-induced gene responses, we investigated the distribution of TCPOBOP-responsive RefSeq coding and long noncoding RNA (lncRNA) genes across the megabase-scale topologically associating domains (TADs) that segment the genome, and which provide a structural framework that functionally constrains enhancer-promoter interactions. We show that a subset of TCPOBOP-responsive genes cluster within TADs, and that TCPOBOP-induced genes and TCPOBOP-repressed genes are often found in different TADs. Further, using DNase-seq and DNase hypersensitivity site (DHS) analysis, we identified several thousand genomic regions (ΔDHS) where short-term exposure to TCPOBOP induces localized changes (increases or decreases) in mouse liver chromatin accessibility, many of which cluster in TADs together with TCPOBOP-responsive genes. Sites of chromatin opening were highly enriched nearby genes induced by TCPOBOP and chromatin closing was highly enriched nearby genes repressed by TCPOBOP, consistent with TCPOBOP-responsive ΔDHS serving as enhancers and promoters that positively regulate CAR-responsive genes. Gene expression changes lagged behind chromatin opening or closing for a subset of TCPOBOP-responsive ΔDHS. ΔDHS that were specifically responsive to TCPOBOP in male liver were significantly enriched for genomic regions with a basal male bias in chromatin accessibility; however, the male-biased response of hepatocellular carcinoma-related genes to TCPOBOP was not associated with a correspondingly male-biased ΔDHS response. These studies elucidate the genome-wide organization of CAR-responsive genes and of the thousands of associated genomic sites where TCPOBOP exposure induces both rapid and persistent changes in chromatin accessibility.

Using 3D epigenomic maps of primary olfactory neuronal cells from living individuals to understand gene regulation

As part of PsychENCODE, we developed a three-dimensional (3D) epigenomic map of primary cultured neuronal cells derived from olfactory neuroepithelium (CNON). We mapped topologically associating domains and high-resolution chromatin interactions using Hi-C and identified regulatory elements using chromatin immunoprecipitation and nucleosome positioning assays. Using epigenomic datasets from biopsies of 63 living individuals, we found that epigenetic marks at distal regulatory elements are more variable than marks at proximal regulatory elements. By integrating genotype and metadata, we identified enhancers that have different levels corresponding to differences in genetic variation, gender, smoking, and schizophrenia. Motif searches revealed that many CNON enhancers are bound by neuronal-related transcription factors. Last, we combined 3D epigenomic maps and gene expression profiles to predict enhancer-target gene interactions on a genome-wide scale. This study not only provides a framework for understanding individual epigenetic variation using a primary cell model system but also contributes valuable data resources for epigenomic studies of neuronal epithelium.

Genome-wide map of proximity linkage to renin proximal promoter in rat

A challenge to understanding enhancer-gene relationships is that enhancers are not always sequentially close to the gene they regulate. Physical proximity mapping through sequencing can provide an unbiased view of the chromatin close to the proximal promoter of the renin gene ( Ren). Our objective was to determine genomic regions that physically interact with the renin proximal promoter, using two different genetic backgrounds, the Dahl salt sensitive and normotensive SS-13BN, which have been shown to have different regulation of plasma renin in vivo. The chromatin conformation capture method with sequencing focused at the Ren proximal promoter in rat-derived cardiac endothelial cells was used. Cells were fixed, chromatin close to the Ren promoter was captured, and fragments were sequenced. The clustering of mapped reads produced a genome-wide map of chromatin in contact with the Ren promoter. The largest number of contacts was found on chromosome 13, the chromosome with Ren, and contacts were found on all other chromosomes except chromosome X. These contacts were significantly enriched with genes positively correlated with Ren expression and with mapped quantitative trait loci associated with blood pressure, cardiovascular, and renal phenotypes. The results were reproducible in an independent biological replicate. The findings reported here represent the first map between a critical cardiovascular gene and physical interacting loci throughout the genome and will provide the basis for several new directions of research.

Signaling by Steroid Hormones in the 3D Nuclear Space

Initial studies showed that ligand-activated hormone receptors act by binding to the proximal promoters of individual target genes. Genome-wide studies have now revealed that regulation of transcription by steroid hormones mainly depends on binding of the receptors to distal regulatory elements. Those distal elements, either enhancers or silencers, act on the regulation of target genes by chromatin looping to the gene promoters. In the nucleus, this level of chromatin folding is integrated within dynamic higher orders of genome structures, which are organized in a non-random fashion. Terminally differentiated cells exhibit a tissue-specific three-dimensional (3D) organization of the genome that favors or restrains the activity of transcription factors and modulates the function of steroid hormone receptors, which are transiently activated upon hormone exposure. Conversely, integration of the hormones signal may require modifications of the 3D organization to allow appropriate transcriptional outcomes. In this review, we summarize the main levels of organization of the genome, review how they can modulate the response to steroids in a cell specific manner and discuss the role of receptors in shaping and rewiring the structure in response to hormone. Taking into account the dynamics of 3D genome organization will contribute to a better understanding of the pleiotropic effects of steroid hormones in normal and cancer cells.

Phytochemicals Targeting Estrogen Receptors: Beneficial Rather Than Adverse Effects?

In mammals, the effects of estrogen are mainly mediated by two different estrogen receptors, ERα and ERβ. These proteins are members of the nuclear receptor family, characterized by distinct structural and functional domains, and participate in the regulation of different biological processes, including cell growth, survival and differentiation. The two estrogen receptor (ER) subtypes are generated from two distinct genes and have partially distinct expression patterns. Their activities are modulated differently by a range of natural and synthetic ligands. Some of these ligands show agonistic or antagonistic effects depending on ER subtype and are described as selective ER modulators (SERMs). Accordingly, a few phytochemicals, called phytoestrogens, which are synthesized from plants and vegetables, show low estrogenic activity or anti-estrogenic activity with potentially anti-proliferative effects that offer nutraceutical or pharmacological advantages. These compounds may be used as hormonal substitutes or as complements in breast cancer treatments. In this review, we discuss and summarize the in vitro and in vivo effects of certain phytoestrogens and their potential roles in the interaction with estrogen receptors.

Genomic effects of glucocorticoids

Glucocorticoids and their receptor (GR) have been an important area of research because of their pleiotropic physiological functions and extensive use in the clinic. In addition, the association between GR and glucocorticoids, which is highly specific, leads to rapid nuclear translocation where GR associates with chromatin to regulate gene transcription. This simplified model system has been instrumental for studying the complexity of transcription regulation processes occurring at chromatin. In this review we discuss our current understanding of GR action that has been enhanced by recent developments in genome wide measurements of chromatin accessibility, histone marks, chromatin remodeling and 3D chromatin structure in various cell types responding to glucocorticoids.