Close encounters of the 3C kind: long-range chromatin interactions and transcriptional regulation

Enhancer Function in the 3D Genome

In this review, we consider various aspects of enhancer functioning in the context of the 3D genome. Particular attention is paid to the mechanisms of enhancer-promoter communication and the significance of the spatial juxtaposition of enhancers and promoters in 3D nuclear space. A model of an activator chromatin compartment is substantiated, which provides the possibility of transferring activating factors from an enhancer to a promoter without establishing direct contact between these elements. The mechanisms of selective activation of individual promoters or promoter classes by enhancers are also discussed.

Investigating the genetic architecture of eye colour in a Canadian cohort

Eye color is highly variable in populations with European ancestry, ranging from low to high quantities of melanin in the iris. Polymorphisms in the HERC2/OCA2 locus have the largest effect on eye color in these populations, although other genomic regions also influence eye color. We performed genome-wide association studies of eye color in a Canadian cohort of European ancestry (N = 5,641) and investigated candidate causal variants. We uncovered several candidate causal signals in the HERC2/OCA2 region, whereas other loci likely harbor a single causal signal. We observed colocalization of eye color signals with the expression or methylation profiles of cultured primary melanocytes. Genetic correlations of eye and hair color suggest high genome-wide pleiotropy, but locus-level differences in the genetic architecture of both traits. Overall, we provide a better picture of the polymorphisms underpinning eye color variation, which may be a consequence of specific molecular processes in the iris melanocytes.

•

Genome-wide association studies of eye color in 5,641 participants

•

Multiple independent candidate causal variants were identified across HERC2/OCA2

•

Single candidate causal variants observed on or near IRF4, SLC24A4, TYR, and TYRP1

•

Colocalization of eye color signals with expression and methylation profiles

Allele-specific long-distance regulation dictates IL-32 isoform switching and mediates susceptibility to HIV-1

We integrated data obtained from HIV-1 genome-wide association studies with T cell-derived epigenome data and found that the noncoding intergenic variant rs4349147, which is statistically associated with HIV-1 acquisition, is located in a CD4⁺ T cell-specific deoxyribonuclease I hypersensitive region, suggesting regulatory potential for this variant. Deletion of the rs4349147 element in Jurkat cells strongly reduced expression of interleukin-32 (IL-32), approximately 10-kb upstream, and chromosome conformation capture assays identified a chromatin loop between rs4349147 and the IL-32 promoter validating its function as a long-distance enhancer. We generated single rs4349147-A or rs4349147-G allele clones and demonstrated that IL-32 enhancer activity and interaction with the IL-32 promoter are strongly allele dependent; rs4349147 -/A cells display reduced IL-32 expression and altered chromatin conformation as compared to rs4349147 G/- cells. Moreover, RNA sequencing demonstrated that rs4349147 G/- cells express a lower relative ratio of IL-32α to non-α isoforms than rs4349147 -/A cells and display increased expression of lymphocyte activation factors rendering them more prone to infection with HIV-1. In agreement, in primary CD4⁺ T cells, both treatment with recombinant IL-32γ (rIL-32γ) but not rIL-32α, and exogenous lentiviral overexpression of IL-32γ or IL-32β but not IL-32α resulted in a proinflammatory T cell cytokine environment concomitant with increased susceptibility to HIV infection. Our data demonstrate that rs4349147-G promotes transcription of non-IL-32α isoforms, generating a proinflammatory environment more conducive to HIV infection. This study provides a mechanistic link between a HIV-associated noncoding DNA variant and the expression of different IL-32 isoforms that display discrete anti-HIV properties.

Runx1 Structure and Function in Blood Cell Development

RUNX transcription factors belong to a highly conserved class of transcriptional regulators which play various roles in the development of the majority of metazoans. In this review we focus on the founding member of the family, RUNX1, and its role in the transcriptional control of blood cell development in mammals. We summarize data showing that RUNX1 functions both as activator and repressor within a chromatin environment, a feature that requires its interaction with multiple other transcription factors and co-factors. Furthermore, we outline how RUNX1 works together with other factors to reshape the epigenetic landscape and the three-dimensional structure of gene loci within the nucleus. Finally, we review how aberrant forms of RUNX1 deregulate blood cell development and cause hematopoietic malignancies.

Allele-specific transcriptional regulation of IRF4 in melanocytes is mediated by chromatin looping of the intronic rs12203592 enhancer to the IRF4 promoter

The majority of significant single-nucleotide polymorphisms (SNPs) identified with genome-wide association studies are located in non-coding regions of the genome; it is therefore possible that they are involved in transcriptional regulation of a nearby gene rather than affecting an encoded protein's function. Previously, it was demonstrated that the SNP rs12203592, located in intron 4 of the IRF4 gene, is strongly associated with human skin pigmentation and modulates an enhancer element that controls expression of IRF4. In our study, we investigated the allele-specific effect of rs12203592 on IRF4 expression in epidermal skin samples and in melanocytic cells from donors of different skin color. We focused on the characteristics and activity of the enhancer, and on long-range chromatin interactions in melanocytic cells homozygous and heterozygous for rs12203592. We found that, irrespective of the trans-activating environment, IRF4 transcription is strongly correlated with the allelic status of rs12203592, the activity of the rs12203592 enhancer and that the chromatin features depend on the rs12203592 genotype. Furthermore, we demonstrate that the rs12203592 enhancer physically interacts with the IRF4 promoter through an allele-dependent chromatin loop, and suggest that subsequent allele-specific activation of IRF4 transcription is stabilized by another allele-specific loop from the rs12203592 enhancer to an additional regulatory element in IRF4. We conclude that the non-coding SNP rs12203592 is located in a regulatory region and affects a wide range of enhancer characteristics, resulting into modulation of the enhancer's activity, its interaction with the IRF4 promoter and subsequent allele-specific transcription of IRF4. Our findings provide another example of a non-coding SNP affecting skin color by modulating enhancer-mediated transcriptional regulation.

Genetic variation in regulatory DNA elements: the case of OCA2 transcriptional regulation

Mutations within the OCA2 gene or the complete absence of the OCA2 protein leads to oculocutaneous albinism type 2. The OCA2 protein plays a central role in melanosome biogenesis, and it is a strong determinant of the eumelanin content in melanocytes. Transcript levels of the OCA2 gene are strongly correlated with pigmentation intensities. Recent studies demonstrated that the transcriptional level of OCA2 is to a large extent determined by the noncoding SNP rs12913832 located 21.5 kb upstream of the OCA2 gene promoter. In this review, we discuss current hypotheses and the available data on the mechanism of OCA2 transcriptional regulation and how this is influenced by genetic variation. Finally, we will explore how future epigenetic studies can be used to advance our insight into the functional biology that connects genetic variation to human pigmentation.

In search of the determinants of enhancer-promoter interaction specificity

Although it was originally believed that enhancers activate only the nearest promoter, recent global analyses enabled by high-throughput technology suggest that the network of enhancer-promoter interactions is far more complex. The mechanisms that determine the specificity of enhancer-promoter interactions are still poorly understood, but they are thought to include biochemical compatibility, constraints imposed by the three-dimensional architecture of chromosomes, insulator elements, and possibly the effects of local chromatin composition. In this review, we assess the current insights into these determinants, and highlight the functional genomic approaches that will lead the way towards better mechanistic understanding.

In search of the determinants of enhancer-promoter interaction specificity

Although it was originally believed that enhancers activate only the nearest promoter, recent global analyses enabled by high-throughput technology suggest that the network of enhancer-promoter interactions is far more complex. The mechanisms that determine the specificity of enhancer-promoter interactions are still poorly understood, but they are thought to include biochemical compatibility, constraints imposed by the three-dimensional architecture of chromosomes, insulator elements, and possibly the effects of local chromatin composition. In this review, we assess the current insights into these determinants, and highlight the functional genomic approaches that will lead the way towards better mechanistic understanding.

The clustering of CpG islands may constitute an important determinant of the 3D organization of interphase chromosomes

We used the 4C-Seq technique to characterize the genome-wide patterns of spatial contacts of several CpG islands located on chromosome 14 in cultured chicken lymphoid and erythroid cells. We observed a clear tendency for the spatial clustering of CpG islands present on the same and different chromosomes, regardless of the presence or absence of promoters within these CpG islands. Accordingly, we observed preferential spatial contacts between Sp1 binding motifs and other GC-rich genomic elements, including the DNA sequence motifs capable of forming G-quadruplexes. However, an anchor placed in a gene/CpG island-poor area formed spatial contacts with other gene/CpG island-poor areas on chromosome 14 and other chromosomes. These results corroborate the two-compartment model of the spatial organization of interphase chromosomes and suggest that the clustering of CpG islands constitutes an important determinant of the 3D organization of the eukaryotic genome in the cell nucleus. Using the ChIP-Seq technique, we mapped the genome-wide CTCF deposition sites in the chicken lymphoid and erythroid cells that were used for the 4C analysis. We observed a good correlation between the density of CTCF deposition sites and the level of 4C signals for the anchors located in CpG islands but not for an anchor located in a gene desert. It is thus possible that CTCF contributes to the clustering of CpG islands observed in our experiments.

DNA G-quadruplex formation in response to remote downstream transcription activity: long-range sensing and signal transducing in DNA double helix

G-quadruplexes, four-stranded structures formed by Guanine-rich nucleic acids, are implicated in many physiological and pathological processes. G-quadruplex-forming sequences are abundant in genomic DNA, and G-quadruplexes have recently been shown to exist in the genome of mammalian cells. However, how G-quadruplexes are formed in the genomes remains largely unclear. Here, we show that G-quadruplex formation can be remotely induced by downstream transcription events that are thousands of base pairs away. The induced G-quadruplexes alter protein recognition and cause transcription termination at the local region. These results suggest that a G-quadruplex-forming sequence can serve as a sensor or receiver to sense remote DNA tracking activity in response to the propagation of mechanical torsion in a DNA double helix. We propose that the G-quadruplex formation may provide a mean for long-range sensing and communication between distal genomic locations to coordinate regulatory transactions in genomic DNA.

Disclosure of a structural milieu for the proximity ligation reveals the elusive nature of an active chromatin hub

The current progress in the study of the spatial organization of interphase chromosomes became possible owing to the development of the chromosome conformation capture (3C) protocol. The crucial step of this protocol is the proximity ligation—preferential ligation of DNA fragments assumed to be joined within nuclei by protein bridges and solubilized as a common complex after formaldehyde cross-linking and DNA cleavage. Here, we show that a substantial, and in some cases the major, part of DNA is not solubilized from cross-linked nuclei treated with restriction endonuclease(s) and sodium dodecyl sulphate and that this treatment neither causes lysis of the nucleus nor drastically affects its internal organization. Analysis of the ligation frequencies of the mouse β-globin gene domain DNA fragments demonstrated that the previously reported 3C signals were generated predominantly, if not exclusively, in the insoluble portion of the 3C material. The proximity ligation thus occurs within the cross-linked chromatin cage in non-lysed nuclei. The finding does not compromise the 3C protocol but allows the consideration of an active chromatin hub as a folded chromatin domain or a nuclear compartment rather than a rigid complex of regulatory elements.

Spatial localization of co-regulated genes exceeds genomic gene clustering in the Saccharomyces cerevisiae genome

While it has been long recognized that genes are not randomly positioned along the genome, the degree to which its 3D structure influences the arrangement of genes has remained elusive. In particular, several lines of evidence suggest that actively transcribed genes are spatially co-localized, forming transcription factories; however, a generalized systematic test has hitherto not been described. Here we reveal transcription factories using a rigorous definition of genomic structure based on Saccharomyces cerevisiae chromosome conformation capture data, coupled with an experimental design controlling for the primary gene order. We develop a data-driven method for the interpolation and the embedding of such datasets and introduce statistics that enable the comparison of the spatial and genomic densities of genes. Combining these, we report evidence that co-regulated genes are clustered in space, beyond their observed clustering in the context of gene order along the genome and show this phenomenon is significant for 64 out of 117 transcription factors. Furthermore, we show that those transcription factors with high spatially co-localized targets are expressed higher than those whose targets are not spatially clustered. Collectively, our results support the notion that, at a given time, the physical density of genes is intimately related to regulatory activity.

A gene loop containing the floral repressor FLC is disrupted in the early phase of vernalization

Gene activation in eukaryotes frequently involves interactions between chromosomal regions. We have investigated whether higher-order chromatin structures are involved in the regulation of the Arabidopsis floral repressor gene FLC, a target of several chromatin regulatory pathways. Here, we identify a gene loop involving the physical interaction of the 5' and 3' flanking regions of the FLC locus using chromosome conformation capture. The FLC loop is unaffected by mutations disrupting conserved chromatin regulatory pathways leading to very different expression states. However, the loop is disrupted during vernalization, the cold-induced, Polycomb-dependent epigenetic silencing of FLC. Loop disruption parallels timing of the cold-induced FLC transcriptional shut-down and upregulation of FLC antisense transcripts, but does not need a cold-induced PHD protein required for the epigenetic silencing. We suggest that gene loop disruption is an early step in the switch from an expressed to a Polycomb-silenced state.

LCR 5' hypersensitive site specificity for globin gene activation within the active chromatin hub

The DNaseI hypersensitive sites (HSs) of the human β-globin locus control region (LCR) may function as part of an LCR holocomplex within a larger active chromatin hub (ACH). Differential activation of the globin genes during development may be controlled in part by preferential interaction of each gene with specific individual HSs during globin gene switching, a change in conformation of the LCR holocomplex, or both. To distinguish between these possibilities, human β-globin locus yeast artificial chromosome (β-YAC) lines were produced in which the ε-globin gene was replaced with a second marked β-globin gene (β^m), coupled to an intact LCR, a 5′HS3 complete deletion (5′ΔHS3) or a 5′HS3 core deletion (5′ΔHS3c). The 5′ΔHS3c mice expressed β^m-globin throughout development; γ-globin was co-expressed in the embryonic yolk sac, but not in the fetal liver; and wild-type β-globin was co-expressed in adult mice. Although the 5′HS3 core was not required for β^m-globin expression, previous work showed that the 5′HS3 core is necessary for ε-globin expression during embryonic erythropoiesis. A similar phenotype was observed in 5′HS complete deletion mice, except β^m-globin expression was higher during primitive erythropoiesis and γ-globin expression continued into fetal definitive erythropoiesis. These data support a site specificity model of LCR HS-globin gene interaction.

Super-resolution imaging reveals three-dimensional folding dynamics of the β-globin locus upon gene activation

The chromatin architecture is constantly changing because of cellular processes such as proliferation, differentiation and changes in the expression profile during gene activation or silencing. Unravelling the changes that occur in the chromatin structure during these processes has been a topic of interest for many years. It is known that gene activation of large gene loci is thought to occur by means of an active looping mechanism. It was also shown for the β-globin locus that the gene promoter interacts with an active chromatin hub by means of an active looping mechanism. This means that the locus changes in three-dimensional (3D) nuclear volume and chromatin shape. As a means of visualizing and measuring these dynamic changes in chromatin structure of the β-globin locus, we used a 3D DNA-FISH method in combination with 3D image acquisition to volume render fluorescent signals into 3D objects. These 3D chromatin structures were geometrically analysed, and results prior to and after gene activation were quantitatively compared. Confocal and super-resolution imaging revealed that the inactive locus occurs in several different conformations. These conformations change in shape and surface structure upon cell differentiation into a more folded and rounded structure that has a substantially smaller size and volume. These physical measurements represent the first non-biochemical evidence that, upon gene activation, an actively transcribing chromatin hub is formed by means of additional chromatin looping.

Hypermethylation of the gene LARP2 for noninvasive prenatal diagnosis of β-thalassemia based on DNA methylation profile

In order to identify epigenetic markers of β-thalassemia, a genome-wide profiling method named differential methylation hybridization was used to search these differentially methylated genes. Unsupervised hierarchical clustering and molecular annotation system were used to analyze the data, and methylation-specific PCR and real-time PCR were used to confirm the differentially methylated genes. This system was validated by detecting 13 cases, 10 of which were homo-zygous β-thalassaemia. Totally 113 genes were identified as methlyation-enriched genes (ratio ≥ 2.0, P < 0.05) and 96 genes were identified as hypomethylated genes in both groups (ratio ≤ 0.5, P < 0.05). The promoter of the gene of La ribonucleoprotein domain family (LARP2) was significantly hypermethylated in β-thalassemia, and the expression of LARP2 was significantly lower in β-thalassemia. Hypermethylation of the LARP2 promoter was correlated with its lower expression in β-thalassemia and our chip-based DNA methylation detection system can provide earlier diagnosis of β-thalassemia using this epigenetic marker.

HERC2 rs12913832 modulates human pigmentation by attenuating chromatin-loop formation between a long-range enhancer and the OCA2 promoter

Pigmentation of skin, eye, and hair reflects some of the most evident common phenotypes in humans. Several candidate genes for human pigmentation are identified. The SNP rs12913832 has strong statistical association with human pigmentation. It is located within an intron of the nonpigment gene HERC2, 21 kb upstream of the pigment gene OCA2, and the region surrounding rs12913832 is highly conserved among animal species. However, the exact functional role of HERC2 rs12913832 in human pigmentation is unknown. Here we demonstrate that the HERC2 rs12913832 region functions as an enhancer regulating OCA2 transcription. In darkly pigmented human melanocytes carrying the rs12913832 T-allele, we detected binding of the transcription factors HLTF, LEF1, and MITF to the HERC2 rs12913832 enhancer, and a long-range chromatin loop between this enhancer and the OCA2 promoter that leads to elevated OCA2 expression. In contrast, in lightly pigmented melanocytes carrying the rs12913832 C-allele, chromatin-loop formation, transcription factor recruitment, and OCA2 expression are all reduced. Hence, we demonstrate that allelic variation of a common noncoding SNP located in a distal regulatory element not only disrupts the regulatory potential of this element but also affects its interaction with the relevant promoter. We provide the key mechanistic insight that allele-dependent differences in chromatin-loop formation (i.e., structural differences in the folding of gene loci) result in differences in allelic gene expression that affects common phenotypic traits. This concept is highly relevant for future studies aiming to unveil the functional basis of genetically determined phenotypes, including diseases.

RUNX1 regulates the CD34 gene in haematopoietic stem cells by mediating interactions with a distal regulatory element

The transcription factor RUNX1 is essential to establish the haematopoietic gene expression programme; however, the mechanism of how it activates transcription of haematopoietic stem cell (HSC) genes is still elusive. Here, we obtained novel insights into RUNX1 function by studying regulation of the human CD34 gene, which is expressed in HSCs. Using transgenic mice carrying human CD34 PAC constructs, we identified a novel downstream regulatory element (DRE), which is bound by RUNX1 and is necessary for human CD34 expression in long-term (LT)-HSCs. Conditional deletion of Runx1 in mice harbouring human CD34 promoter-DRE constructs abrogates human CD34 expression. We demonstrate by chromosome conformation capture assays in LT-HSCs that the DRE physically interacts with the human CD34 promoter. Targeted mutagenesis of RUNX binding sites leads to perturbation of this interaction and decreased human CD34 expression in LT-HSCs. Overall, our in vivo data provide novel evidence about the role of RUNX1 in mediating interactions between distal and proximal elements of the HSC gene CD34.

From nucleosome to chromosome: a dynamic organization of genetic information

Gene activity is controlled at different levels of chromatin organization, which involve genomic sequences, nucleosome structure, chromatin folding and chromosome arrangement. These levels are interconnected and influence each other. At the basic level nucleosomes generally occlude the DNA sequence from interacting with DNA-binding proteins. Evidently, nucleosome positioning is a major factor in gene control and chromatin organization. Understanding the biological rules that govern the deposition and removal of the nucleosomes to and from the chromatin fiber is the key to understanding gene regulation and chromatin organization. In this review we describe and discuss the relationship between the different levels of chromatin organization in plants and animals.

In the loop: long range chromatin interactions and gene regulation

Enhancers, silencer and insulators are DNA elements that play central roles in regulation of the genome that are crucial for development and differentiation. In metazoans, these elements are often separated from target genes by distances that can reach 100  Kb. How regulation can be accomplished over long distances has long been intriguing. Current data indicate that although the mechanisms by which these diverse regulatory elements affect gene transcription may vary, an underlying feature is the establishment of close contacts or chromatin loops. With the generalization of this principle, new questions emerge, such as how the close contacts are formed and stabilized and, importantly, how they contribute to the regulation of transcriptional output at target genes. This review will concentrate on examples where a functional role and a mechanistic understanding has been explored for loops formed between genes and their regulatory elements or among the elements themselves.

Impact of nuclear organization and dynamics on epigenetic regulation in the central nervous system: implications for neurological disease states

Epigenetic mechanisms that are highly responsive to interoceptive and environmental stimuli mediate the proper execution of complex genomic programs, such as cell type-specific gene transcription and posttranscriptional RNA processing, and are increasingly thought to be important for modulating the development, homeostasis, and plasticity of the central nervous system (CNS). These epigenetic processes include DNA methylation, histone modifications, and chromatin remodeling, all of which play roles in neural cellular diversity, connectivity, and plasticity. Further, large-scale transcriptomic analyses have revealed that the eukaryotic genome is pervasively transcribed, forming interleaved protein-coding RNAs and regulatory nonprotein-coding RNAs (ncRNAs), which act through a broad array of molecular mechanisms. Most of these ncRNAs are transcribed in a cell type- and developmental stage-specific manner in the CNS. A broad array of posttranscriptional processes, such as RNA editing and transport, can modulate the functions of both protein-coding RNAs and ncRNAs. Additional studies implicate nuclear organization and dynamics in mediating epigenetic regulation. The compartmentalization of DNA sequences and other molecular machinery into functional nuclear domains, such as transcription factories, Cajal bodies, promyelocytic leukemia nuclear bodies, nuclear speckles, and paraspeckles, some of which are found prominently in neural cells, is associated with regulation of transcriptional activity and posttranscriptional RNA processing. These observations suggest that genomic architecture and RNA biology in the CNS are much more complex and nuanced than previously appreciated. Increasing evidence now suggests that most, if not all, human CNS diseases are associated with either primary or secondary perturbations in one or more aspects of the epigenome. In this review, we provide an update of our emerging understanding of genomic architecture, RNA biology, and nuclear organization and highlight the interconnected roles that deregulation of these factors may play in diverse CNS disorders.