Cell type specificity of chromatin organization mediated by CTCF and cohesin

Effective modeling of the chromatin structure by coarse-grained methods

The interphase chromatin structure is extremely complex, precise and dynamic. Experimental methods can only show the frequency of interaction of the various parts of the chromatin. Therefore, it is extremely important to develop theoretical methods to predict the chromatin structure. In this publication, we implemented an extended version of the SBS model described by Barbieri et al. and created the ChroMC program that is easy to use and freely available (https://github.com/regulomics/chroMC) to other users. We also describe the necessary factors for the effective modeling of the chromatin structure in Drosophila melanogaster. We compared results of chromatin structure predictions using two methods: Monte Carlo and Molecular Dynamic. Our simulations suggest that incorporating black, non-reactive chromatin is necessary for successful prediction of chromatin structure, while the loop extrusion model with a long range attraction potential or Lennard-Jones (with local attraction force) as well as using Hi-C data as input are not essential for the basic structure reconstruction. We also proposed a new way to calculate the similarity of the properties of contact maps including the calculation of local similarity.Communicated by Ramaswamy H. Sarma.

Impact of Genetic Variations on Thromboembolic Risk in Saudis with Sickle Cell Disease

BACKGROUND

Sickle cell disease (SCD) is a Mendelian disease characterized by multigenic phenotypes. Previous reports indicated a higher rate of thromboembolic events (TEEs) in SCD patients. A number of candidate polymorphisms in certain genes (e.g., FVL, PRT, and MTHFR) were previously reported as risk factors for TEEs in different clinical conditions. This study aimed to genotype these genes and other loci predicted to underlie TEEs in SCD patients.

METHODOLOGY

A multi-center genome-wide association study (GWAS) involving Saudi SCD adult patients with a history of TEEs (n = 65) and control patients without TEE history (n = 285) was performed. Genotyping used the 10× Affymetrix Axiom array, which includes 683,030 markers. Fisher's exact test was used to generate p-values of TEE associations with each single-nucleotide polymorphism (SNP). The haplotype analysis software tool version 1.05, designed by the University of Göttingen, Germany, was used to identify the common inherited haplotypes.

RESULTS

No association was identified between the targeted single-nucleotide polymorphism rs1801133 in MTHFR and TEEs in SCD (p = 0.79). The allele frequency of rs6025 in FVL and rs1799963 in PRT in our cohort was extremely low (<0.01); thus, both variants were excluded from the analysis as no meaningful comparison was possible. In contrast, the GWAS analysis showed novel genome-wide associations (p < 5 × 10-8) with seven signals; five of them were located on Chr 11 (rs35390334, rs331532, rs317777, rs147062602, and rs372091), one SNP on Chr 20 (rs139341092), and another on Chr 9 (rs76076035). The other 34 SNPs located on known genes were also detected at a signal threshold of p < 5 × 10-6. Seven of the identified variants are located in olfactory receptor family 51 genes (OR51B5, OR51V1, OR51A1P, and OR51E2), and five variants were related to family 52 genes (OR52A5, OR52K1, OR52K2, and OR52T1P). The previously reported association between rs5006884-A in OR51B5 and fetal hemoglobin (HbF) levels was confirmed in our study, which showed significantly lower levels of HbF (p = 0.002) and less allele frequency (p = 0.003) in the TEE cases than in the controls. The assessment of the haplotype inheritance pattern involved the top ten significant markers with no LD (rs353988334, rs317777, rs14788626882, rs49188823, rs139349992, rs76076035, rs73395847, rs1368823, rs8888834548, and rs1455957). A haplotype analysis revealed significant associations between two haplotypes (a risk, TT-AA-del-AA-ins-CT-TT-CC-CC-AA, and a reverse protective, CC-GG-ins-GG-del-TT-CC-TT-GG-GG) and TEEs in SCD (p = 0.024, OR = 6.16, CI = 1.34-28.24, and p = 0.019, OR = 0.33, CI = 0.13-0.85, respectively).

CONCLUSIONS

Seven markers showed novel genome-wide associations; two of them were exonic variants (rs317777 in OLFM5P and rs147062602 in OR51B5), and less significant associations (p < 5 × 10-6) were identified for 34 other variants in known genes with TEEs in SCD. Moreover, two 10-SNP common haplotypes were determined with contradictory effects. Further replication of these findings is needed.

Understanding fundamental principles of enhancer biology at a model locus: Analysing the structure and function of an enhancer cluster at the α-globin locus

Despite ever-increasing accumulation of genomic data, the fundamental question of how individual genes are switched on during development, lineage-specification and differentiation is not fully answered. It is widely accepted that this involves the interaction between at least three fundamental regulatory elements: enhancers, promoters and insulators. Enhancers contain transcription factor binding sites which are bound by transcription factors (TFs) and co-factors expressed during cell fate decisions and maintain imposed patterns of activation, at least in part, via their epigenetic modification. This information is transferred from enhancers to their cognate promoters often by coming into close physical proximity to form a 'transcriptional hub' containing a high concentration of TFs and co-factors. The mechanisms underlying these stages of transcriptional activation are not fully explained. This review focuses on how enhancers and promoters are activated during differentiation and how multiple enhancers work together to regulate gene expression. We illustrate the currently understood principles of how mammalian enhancers work and how they may be perturbed in enhanceropathies using expression of the α-globin gene cluster during erythropoiesis, as a model.

RNA polymerase II pausing regulates chromatin organization in erythrocytes

Chicken erythrocytes are nucleated cells often considered to be transcriptionally inactive, although the epigenetic changes and chromatin remodeling that would mediate transcriptional repression and the extent of gene silencing during avian terminal erythroid differentiation are not fully understood. Here, we characterize the changes in gene expression, chromatin accessibility, genome organization and chromatin nuclear disposition during the terminal stages of erythropoiesis in chicken and uncover complex chromatin reorganization at different genomic scales. We observe a robust decrease in transcription in erythrocytes, but a set of genes maintains their expression, including genes involved in RNA polymerase II (Pol II) promoter-proximal pausing. Erythrocytes exhibit a reoriented nuclear architecture, with accessible chromatin positioned towards the nuclear periphery together with the paused RNA Pol II. In erythrocytes, chromatin domains are partially lost genome-wide, except at minidomains retained around paused promoters. Our results suggest that promoter-proximal pausing of RNA Pol II contributes to the transcriptional regulation of the erythroid genome and highlight the role of RNA polymerase in the maintenance of local chromatin organization.

High-throughput Pore-C reveals the single-allele topology and cell type-specificity of 3D genome folding

Canonical three-dimensional (3D) genome structures represent the ensemble average of pairwise chromatin interactions but not the single-allele topologies in populations of cells. Recently developed Pore-C can capture multiway chromatin contacts that reflect regional topologies of single chromosomes. By carrying out high-throughput Pore-C, we reveal extensive but regionally restricted clusters of single-allele topologies that aggregate into canonical 3D genome structures in two human cell types. We show that fragments in multi-contact reads generally coexist in the same TAD. In contrast, a concurrent significant proportion of multi-contact reads span multiple compartments of the same chromatin type over megabase distances. Synergistic chromatin looping between multiple sites in multi-contact reads is rare compared to pairwise interactions. Interestingly, the single-allele topology clusters are cell type-specific even inside highly conserved TADs in different types of cells. In summary, HiPore-C enables global characterization of single-allele topologies at an unprecedented depth to reveal elusive genome folding principles.

Characterization and perturbation of CTCF-mediated chromatin interactions for enhancing myogenic transdifferentiation

Expression of key transcription factors can induce transdifferentiation in somatic cells; however, this conversion is usually incomplete due to undefined intrinsic barriers. Here, we employ MyoD-induced transdifferentiation of fibroblasts as a model to illustrate the chromatin structures that impede the cell-fate transition. Focusing on the three-dimensional (3D) chromatin interactions, we show that MyoD directly establishes chromatin loops to activate myogenic transcriptional program. Similarly, dynamic changes of CTCF-mediated chromatin interactions are favorable for fibroblast-to-myoblast conversion. However, a substantial portion of CTCF-mediated chromatin interactions remain stable, and the associated genes are steady in expression and enriched for fibroblast function that may restrict cell-identity transformation. Temporal CTCF depletion can interrupt the resistant chromatin loops to enhance myogenic transdifferentiation in mice, pig, and chicken fibroblasts. Therefore, during induced transdifferentiation, the transcription factor can directly reorganize the 3D chromatin interactions, and perturbation of CTCF-mediated genome topology may resolve the limitations of cell fate transitions.

Interrogating cell type-specific cooperation of transcriptional regulators in 3D chromatin

Context-specific activities of transcription regulators (TRs) in the nucleus modulate spatiotemporal gene expression precisely. Using the largest ChIP-seq data and chromatin loops in the human K562 cell line, we initially interrogated TR cooperation in 3D chromatin via a graphical model and revealed many known and novel TRs manipulating context-specific pathways. To explore TR cooperation across broad tissue/cell types, we systematically leveraged large-scale open chromatin profiles, computational footprinting, and high-resolution chromatin interactions to investigate tissue/cell type-specific TR cooperation. We first delineated a landscape of TR cooperation across 40 human tissue/cell types. Network modularity analyses uncovered the commonality and specificity of TR cooperation in different conditions. We also demonstrated that TR cooperation information can better interpret the disease-causal variants identified by genome-wide association studies and recapitulate cell states during neural development. Our study characterizes shared and unique patterns of TR cooperation associated with the cell type specificity of gene regulation in 3D chromatin.

•

Computational inference of transcriptional regulator (TR) cooperation in 3D chromatin

•

A landscape of 3D TR cooperation across 40 human tissue/cell types

•

TR cooperation can better interpret the disease-causal variants identified by GWAS

•

Cooperation of certain TRs shapes context-specific gene regulation in cell development

CGGBP1-dependent CTCF-binding sites restrict ectopic transcription

Binding sites of the chromatin regulator protein CTCF function as important landmarks in the human genome. The recently characterized CTCF-binding sites at LINE-1 repeats depend on another repeat-regulatory protein CGGBP1. These CGGBP1-dependent CTCF-binding sites serve as potential barrier elements for epigenetic marks such as H3K9me3. Such CTCF-binding sites are associated with asymmetric H3K9me3 levels as well as RNA levels in their flanks. The functions of these CGGBP1-dependent CTCF-binding sites remain unknown. By performing targeted studies on candidate CGGBP1-dependent CTCF-binding sites cloned in an SV40 promoter-enhancer episomal system we show that these regions act as inhibitors of ectopic transcription from the SV40 promoter. CGGBP1-dependent CTCF-binding sites that recapitulate their genomic function of loss of CTCF binding upon CGGBP1 depletion and H3K9me3 asymmetry in immediate flanks are also the ones that show the strongest inhibition of ectopic transcription. By performing a series of strand-specific reverse transcription PCRs we demonstrate that this ectopic transcription results in the synthesis of RNA from the SV40 promoter in a direction opposite to the downstream reporter gene in a strand-specific manner. The unleashing of the bidirectionality of the SV40 promoter activity and a breach of the transcription barrier seems to depend on depletion of CGGBP1 and loss of CTCF binding proximal to the SV40 promoter. RNA-sequencing reveals that CGGBP1-regulated CTCF-binding sites act as barriers to transcription at multiple locations genome-wide. These findings suggest a role of CGGBP1-dependent binding sites in restricting ectopic transcription.

Dynamics of the host chromatin three-dimensional response to influenza virus infection

The spatial organization of chromatin is known to be highly dynamic in response to environmental stress. However, it remains unknown how chromatin dynamics contributes to or modulates the pathogenesis of immune and infectious diseases. Influenza virus is a single-stranded RNA virus, and transcription and replication of the virus genome occur in the nucleus. Since viral infection is generally associated with virus-driven hijack of the host cellular machineries, influenza virus may utilize and/or affect the nuclear system. In this review article, we focus on recent studies showing that the three-dimensional structure of chromatin changes with influenza virus infection, which affects the pathology of infection. Also, we discuss studies showing the roles of epigenetics in influenza virus infection. Understanding how this affects immune responses may lead to novel strategies to combat immune and infectious diseases.

Chromatin interaction neural network (ChINN): a machine learning-based method for predicting chromatin interactions from DNA sequences

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples.

Suv4-20h2 protects against influenza virus infection by suppression of chromatin loop formation

The spatial organization of chromatin is known to be highly dynamic in response to environmental stress. However, it remains unknown how chromatin dynamics contributes to or modulates disease pathogenesis. Here, we show that upon influenza virus infection, the H4K20me3 methyltransferase Suv4-20h2 binds the viral protein NP, which results in the inactivation of Suv4-20h2 and the dissociation of cohesin from Suv4-20h2. Inactivation of Suv4-20h2 by viral infection or genetic deletion allows the formation of an active chromatin loop at the HoxC8-HoxC6 loci coincident with cohesin loading. HoxC8 and HoxC6 proteins in turn enhance viral replication by inhibiting the Wnt-β-catenin mediated interferon response. Importantly, loss of Suv4-20h2 augments the pathology of influenza infection in vivo. Thus, Suv4-20h2 acts as a safeguard against influenza virus infection by suppressing cohesin-mediated loop formation.

•

H4K20me3 methyltransferase Suv4-20h2 suppresses influenza viral replication

•

Influenza virus NP protein binds to Suv4-20h2 and causes dissociation from cohesin

•

Suv4-20h2 inactivation generates cohesion-mediated loop formation at HoxC8 -HoxC6

•

HoxC8-HoxC6 enhance viral replication by suppressing Wnt/β-catenin signaling

A comparative analysis of chromatin accessibility in cattle, pig, and mouse tissues

Background

Although considerable progress has been made towards annotating the noncoding portion of the human and mouse genomes, regulatory elements in other species, such as livestock, remain poorly characterized. This lack of functional annotation poses a substantial roadblock to agricultural research and diminishes the value of these species as model organisms. As active regulatory elements are typically characterized by chromatin accessibility, we implemented the Assay for Transposase Accessible Chromatin (ATAC-seq) to annotate and characterize regulatory elements in pigs and cattle, given a set of eight adult tissues.

Results

Overall, 306,304 and 273,594 active regulatory elements were identified in pig and cattle, respectively. 71,478 porcine and 47,454 bovine regulatory elements were highly tissue-specific and were correspondingly enriched for binding motifs of known tissue-specific transcription factors. However, in every tissue the most prevalent accessible motif corresponded to the insulator CTCF, suggesting pervasive involvement in 3-D chromatin organization. Taking advantage of a similar dataset in mouse, open chromatin in pig, cattle, and mice were compared, revealing that the conservation of regulatory elements, in terms of sequence identity and accessibility, was consistent with evolutionary distance; whereas pig and cattle shared about 20% of accessible sites, mice and ungulates only had about 10% of accessible sites in common. Furthermore, conservation of accessibility was more prevalent at promoters than at intergenic regions.

Conclusions

The lack of conserved accessibility at distal elements is consistent with rapid evolution of enhancers, and further emphasizes the need to annotate regulatory elements in individual species, rather than inferring elements based on homology. This atlas of chromatin accessibility in cattle and pig constitutes a substantial step towards annotating livestock genomes and dissecting the regulatory link between genome and phenome.

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.

Genome-Wide Transcriptional Regulation of the Long Non-coding RNA Steroid Receptor RNA Activator in Human Erythroblasts

Erythropoiesis of human hematopoietic stem cells (HSCs) maintains generation of red blood cells throughout life. However, little is known how human erythropoiesis is regulated by long non-coding RNAs (lncRNAs). By using ChIRP-seq, we report here that the lncRNA steroid receptor RNA activator (SRA) occupies chromatin, and co-localizes with CTCF, H3K4me3, and H3K27me3 genome-wide in human erythroblast cell line K562. CTCF binding sites that are also occupied by SRA are enriched for either H3K4me3 or H3K27me3. Transcriptome-wide analyses reveal that SRA facilitates expression of erythroid-associated genes, while repressing leukocyte-associated genes in both K562 and CD36-positive primary human proerythroblasts derived from HSCs. We find that SRA-regulated genes are enriched by both CTCF and SRA bindings. Further, silencing of SRA decreases expression of the erythroid-specific markers TFRC and GYPA, and down-regulates expression of globin genes in both K562 and human proerythroblast cells. Taken together, our findings establish that the lncRNA SRA occupies chromatin, and promotes transcription of erythroid genes, therefore facilitating human erythroid transcriptional program.

CGGBP1-regulated cytosine methylation at CTCF-binding motifs resists stochasticity

Background

The human CGGBP1 binds to GC-rich regions and interspersed repeats, maintains homeostasis of stochastic cytosine methylation and determines DNA-binding of CTCF. Interdependence between regulation of cytosine methylation and CTCF occupancy by CGGBP1 remains unknown.

Results

By analyzing methylated DNA-sequencing data obtained from CGGBP1-depleted cells, we report that some transcription factor-binding sites, including CTCF, resist stochastic changes in cytosine methylation. By analysing CTCF-binding sites we show that cytosine methylation changes at CTCF motifs caused by CGGBP1 depletion resist stochastic changes. These CTCF-binding sites are positioned at locations where the spread of cytosine methylation in cis depends on the levels of CGGBP1.

Conclusion

Our findings suggest that CTCF occupancy and functions are determined by CGGBP1-regulated cytosine methylation patterns.

BICORN: An R package for integrative inference of de novo cis-regulatory modules

Genome-wide transcription factor (TF) binding signal analyses reveal co-localization of TF binding sites based on inferred cis-regulatory modules (CRMs). CRMs play a key role in understanding the cooperation of multiple TFs under specific conditions. However, the functions of CRMs and their effects on nearby gene transcription are highly dynamic and context-specific and therefore are challenging to characterize. BICORN (Bayesian Inference of COoperative Regulatory Network) builds a hierarchical Bayesian model and infers context-specific CRMs based on TF-gene binding events and gene expression data for a particular cell type. BICORN automatically searches for a list of candidate CRMs based on the input TF bindings at regulatory regions associated with genes of interest. Applying Gibbs sampling, BICORN iteratively estimates model parameters of CRMs, TF activities, and corresponding regulation on gene transcription, which it models as a sparse network of functional CRMs regulating target genes. The BICORN package is implemented in R (version 3.4 or later) and is publicly available on the CRAN server at https://cran.r-project.org/web/packages/BICORN/index.html.

MAGIC: A tool for predicting transcription factors and cofactors driving gene sets using ENCODE data

Transcriptomic profiling is an immensely powerful hypothesis generating tool. However, accurately predicting the transcription factors (TFs) and cofactors that drive transcriptomic differences between samples is challenging. A number of algorithms draw on ChIP-seq tracks to define TFs and cofactors behind gene changes. These approaches assign TFs and cofactors to genes via a binary designation of 'target', or 'non-target' followed by Fisher Exact Tests to assess enrichment of TFs and cofactors. ENCODE archives 2314 ChIP-seq tracks of 684 TFs and cofactors assayed across a 117 human cell lines under a multitude of growth and maintenance conditions. The algorithm presented herein, Mining Algorithm for GenetIc Controllers (MAGIC), uses ENCODE ChIP-seq data to look for statistical enrichment of TFs and cofactors in gene bodies and flanking regions in gene lists without an a priori binary classification of genes as targets or non-targets. When compared to other TF mining resources, MAGIC displayed favourable performance in predicting TFs and cofactors that drive gene changes in 4 settings: 1) A cell line expressing or lacking single TF, 2) Breast tumors divided along PAM50 designations 3) Whole brain samples from WT mice or mice lacking a single TF in a particular neuronal subtype 4) Single cell RNAseq analysis of neurons divided by Immediate Early Gene expression levels. In summary, MAGIC is a standalone application that produces meaningful predictions of TFs and cofactors in transcriptomic experiments.

Construction of a High-Density Genetic Linkage Map and QTL Mapping for Growth-Related Traits in Takifugu bimaculatus

Takifugu bimaculatus is a euryhaline species, distributed ranging from the southern Yellow Sea to the South China Sea. Their tolerance to a wide range of salinity and temperature, coupled with a desirable firm texture, makes T. bimaculatus a strong candidate for Takifugu aquaculture in subtropics areas. Due to the increasing demand in markets and emerging of the Takifugu aquaculture industry, close attention has been paid to improvement on the T. bimaculatus production. In aquaculture, the great effort has been put into marker-assisted selective breeding, and efficient improvement was realized. However, few genetic resources on T. bimaculatus are provided so far. Aiming at understanding the genetic basis underlying important economic growth traits, facilitating genetic improvement and enriching the genetic resource in T. bimaculatus, we constructed the first genetic linkage map for T. bimaculatus via double digestion restriction-site association DNA sequencing and conducted quantitative traits locus (QTL) mapping for growth-related traits. The map comprised 1976 single nucleotide polymorphism markers distributed on 22 linkage groups (LG), with a total genetic distance of 2039.74 cM. Based on the linkage map, a chromosome-level assembly was constructed whereby we carried out comparative genomics analysis, verifying the high accuracy on contigs ordering of the linkage map. On the other hand, 18 QTLs associated with growth traits were detected on LG6, LG7, LG8, LG10, LG20, and LG21 with phenotypical variance ranging from 15.1 to 56.4%. Candidate genes participating in cartilage development, fat accumulation, and other growth-related regulation activities were identified from these QTLs, including col11a1, foxa2, and thrap3. The linkage map provided a solid foundation for chromosomes assembly and refinement. QTLs reported here unraveled the genomic architecture of some growth traits, which will advance the investigation of aquaculture breeding efforts in T. bimaculatus.

Functionally Annotating Regulatory Elements in the Equine Genome Using Histone Mark ChIP-Seq

One of the primary aims of the Functional Annotation of ANimal Genomes (FAANG) initiative is to characterize tissue-specific regulation within animal genomes. To this end, we used chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map four histone modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in eight prioritized tissues collected as part of the FAANG equine biobank from two thoroughbred mares. Data were generated according to optimized experimental parameters developed during quality control testing. To ensure that we obtained sufficient ChIP and successful peak-calling, data and peak-calls were assessed using six quality metrics, replicate comparisons, and site-specific evaluations. Tissue specificity was explored by identifying binding motifs within unique active regions, and motifs were further characterized by gene ontology (GO) and protein–protein interaction analyses. The histone marks identified in this study represent some of the first resources for tissue-specific regulation within the equine genome. As such, these publicly available annotation data can be used to advance equine studies investigating health, performance, reproduction, and other traits of economic interest in the horse.

Mapping cis-regulatory chromatin contacts in neural cells links neuropsychiatric disorder risk variants to target genes

Mutations in gene regulatory elements have been associated with a wide range of complex neuropsychiatric disorders. However, due to their cell-type specificity and difficulties in characterizing their regulatory targets, the ability to identify causal genetic variants has remained limited. To address these constraints, we perform an integrative analysis of chromatin interactions, open chromatin regions and transcriptomes using promoter capture Hi-C, assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and RNA sequencing, respectively, in four functionally distinct neural cell types: induced pluripotent stem cell (iPSC)-induced excitatory neurons and lower motor neurons, iPSC-derived hippocampal dentate gyrus-like neurons and primary astrocytes. We identify hundreds of thousands of long-range cis-interactions between promoters and distal promoter-interacting regions, enabling us to link regulatory elements to their target genes and reveal putative processes that are dysregulated in disease. Finally, we validate several promoter-interacting regions by using clustered regularly interspaced short palindromic repeats (CRISPR) techniques in human excitatory neurons, demonstrating that CDK5RAP3, STRAP and DRD2 are transcriptionally regulated by physically linked enhancers.

PML modulates H3.3 targeting to telomeric and centromeric repeats in mouse fibroblasts

Targeted deposition of histone variant H3.3 into chromatin is paramount for proper regulation of chromatin integrity, particularly in heterochromatic regions including repeats. We have recently shown that the promyelocytic leukemia (PML) protein prevents H3.3 from being deposited in large heterochromatic PML-associated domains (PADs). However, to what extent PML modulates H3.3 loading on chromatin in other areas of the genome remains unexplored. Here, we examined the impact of PML on targeting of H3.3 to genes and repeat regions that reside outside PADs. We show that loss of PML increases H3.3 deposition in subtelomeric, telomeric, pericentric and centromeric repeats in mouse embryonic fibroblasts, while other repeat classes are not affected. Expression of major satellite, minor satellite and telomeric non-coding transcripts is altered in Pml-null cells. In particular, telomeric Terra transcripts are strongly upregulated, in concordance with a marked reduction in H4K20me3 at these sites. Lastly, for most genes H3.3 enrichment or gene expression outcomes are independent of PML. Our data argue towards the importance of a PML-H3.3 axis in preserving a heterochromatin state at centromeres and telomeres.

Revealing transcription factor and histone modification co-localization and dynamics across cell lines by integrating ChIP-seq and RNA-seq data

Background

Interactions among transcription factors (TFs) and histone modifications (HMs) play an important role in the precise regulation of gene expression. The context specificity of those interactions and further its dynamics in normal and disease remains largely unknown. Recent development in genomics technology enables transcription profiling by RNA-seq and protein’s binding profiling by ChIP-seq. Integrative analysis of the two types of data allows us to investigate TFs and HMs interactions both from the genome co-localization and downstream target gene expression.

Results

We propose a integrative pipeline to explore the co-localization of 55 TFs and 11 HMs and its dynamics in human GM12878 and K562 by matched ChIP-seq and RNA-seq data from ENCODE. We classify TFs and HMs into three types based on their binding enrichment around transcription start site (TSS). Then a set of statistical indexes are proposed to characterize the TF-TF and TF-HM co-localizations. We found that Rad21, SMC3, and CTCF co-localized across five cell lines. High resolution Hi-C data in GM12878 shows that they associate most of the Hi-C peak loci with a specific CTCF-motif “anchor” and supports that CTCF, SMC3, and RAD2 co-localization serves important role in 3D chromatin structure. Meanwhile, 17 TF-TF pairs are highly dynamic between GM12878 and K562. We then build SVM models to correlate high and low expression level of target genes with TF binding and HM strength. We found that H3k9ac, H3k27ac, and three TFs (ELF1, TAF1, and POL2) are predictive with the accuracy about 85~92%.

Conclusion

We propose a pipeline to analyze the co-localization of TF and HM and their dynamics across cell lines from ChIP-seq, and investigate their regulatory potency by RNA-seq. The integrative analysis of two level data reveals new insight for the cooperation of TFs and HMs and is helpful in understanding cell line specificity of TF/HM interactions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5278-5) contains supplementary material, which is available to authorized users.

Predicting CTCF-mediated chromatin interactions by integrating genomic and epigenomic features

The CCCTC-binding zinc-finger protein (CTCF)-mediated network of long-range chromatin interactions is important for genome organization and function. Although this network has been considered largely invariant, we find that it exhibits extensive cell-type-specific interactions that contribute to cell identity. Here, we present Lollipop, a machine-learning framework, which predicts CTCF-mediated long-range interactions using genomic and epigenomic features. Using ChIA-PET data as benchmark, we demonstrate that Lollipop accurately predicts CTCF-mediated chromatin interactions both within and across cell types, and outperforms other methods based only on CTCF motif orientation. Predictions are confirmed computationally and experimentally by Chromatin Conformation Capture (3C). Moreover, our approach identifies other determinants of CTCF-mediated chromatin wiring, such as gene expression within the loops. Our study contributes to a better understanding about the underlying principles of CTCF-mediated chromatin interactions and their impact on gene expression.

CTCF mediates long-range chromatin interactions which are important for genome organization and function. Here, the authors demonstrate that CTCF-mediated interactome exhibits extensive plasticity and present Lollipop, a machine-learning framework which predicts CTCF-mediated long-range interactions using genomic and epigenomic features.

Towards a map of cis-regulatory sequences in the human genome

Accumulating evidence indicates that transcription factor (TF) binding sites, or cis-regulatory elements (CREs), and their clusters termed cis-regulatory modules (CRMs) play a more important role than do gene-coding sequences in specifying complex traits in humans, including the susceptibility to common complex diseases. To fully characterize their roles in deriving the complex traits/diseases, it is necessary to annotate all CREs and CRMs encoded in the human genome. However, the current annotations of CREs and CRMs in the human genome are still very limited and mostly coarse-grained, as they often lack the detailed information of CREs in CRMs. Here, we integrated 620 TF ChIP-seq datasets produced by the ENCODE project for 168 TFs in 79 different cell/tissue types and predicted an unprecedentedly completely map of CREs in CRMs in the human genome at single nucleotide resolution. The map includes 305 912 CRMs containing a total of 1 178 913 CREs belonging to 736 unique TF binding motifs. The predicted CREs and CRMs tend to be subject to either purifying selection or positive selection, thus are likely to be functional. Based on the results, we also examined the status of available ChIP-seq datasets for predicting the entire regulatory genome of humans.

Pan-cancer screen for mutations in non-coding elements with conservation and cancer specificity reveals correlations with expression and survival

Cancer develops by accumulation of somatic driver mutations, which impact cellular function. Mutations in non-coding regulatory regions can now be studied genome-wide and further characterized by correlation with gene expression and clinical outcome to identify driver candidates. Using a new two-stage procedure, called ncDriver, we first screened 507 ICGC whole-genomes from 10 cancer types for non-coding elements, in which mutations are both recurrent and have elevated conservation or cancer specificity. This identified 160 significant non-coding elements, including the TERT promoter, a well-known non-coding driver element, as well as elements associated with known cancer genes and regulatory genes (e.g., PAX5, TOX3, PCF11, MAPRE3). However, in some significant elements, mutations appear to stem from localized mutational processes rather than recurrent positive selection in some cases. To further characterize the driver potential of the identified elements and shortlist candidates, we identified elements where presence of mutations correlated significantly with expression levels (e.g., TERT and CDH10) and survival (e.g., CDH9 and CDH10) in an independent set of 505 TCGA whole-genome samples. In a larger pan-cancer set of 4128 TCGA exomes with expression profiling, we identified mutational correlation with expression for additional elements (e.g., near GATA3, CDC6, ZNF217, and CTCF transcription factor binding sites). Survival analysis further pointed to MIR122, a known marker of poor prognosis in liver cancer. In conclusion, the screen for significant mutation patterns coupled with correlative mutational analysis identified new individual driver candidates and suggest that some non-coding mutations recurrently affect expression and play a role in cancer development.

Mutations in the “non-coding” part of the genome have been identified that could be involved in driving cancer development. Jakob Pedersen, Henrik Hornshøj and colleagues from Aarhus University Hospital in Denmark and MIT in the United States developed a two-stage procedure to identify elements that could be driving cancer development in the part of DNA that does not code for proteins. They conducted statisical analyses on catalogs of tumor genomes to identify recurrent mutations. They then evaluated how specific these mutations were to different cancer types, their predicted functional impact, and their association with gene expression and patient survival. The analyses identified mutations in the non-coding part of cancer genomes that could be driving tumor development, but further analyses on larger sample sets need to be conducted to validate the results, which could provide a basis for biomarker discovery and precision medical treatment.

CTCF, Cohesin, and Chromatin in Human Cancer

It is becoming increasingly clear that eukaryotic genomes are subjected to higher-order chromatin organization by the CCCTC-binding factor/cohesin complex. Their dynamic interactions in three dimensions within the nucleus regulate gene transcription by changing the chromatin architecture. Such spatial genomic organization is functionally important for the spatial disposition of chromosomes to control cell fate during development and differentiation. Thus, the dysregulation of proper long-range chromatin interactions may influence the development of tumorigenesis and cancer progression.

Characteristics of allelic gene expression in human brain cells from single-cell RNA-seq data analysis

Background

Monoallelic expression of autosomal genes has been implicated in human psychiatric disorders. However, there is a paucity of allelic expression studies in human brain cells at the single cell and genome wide levels.

Results

In this report, we reanalyzed a previously published single-cell RNA-seq dataset from several postmortem human brains and observed pervasive monoallelic expression in individual cells, largely in a random manner. Examining single nucleotide variants with a predicted functional disruption, we found that the “damaged” alleles were overall expressed in fewer brain cells than their counterparts, and at a lower level in cells where their expression was detected. We also identified many brain cell type-specific monoallelically expressed genes. Interestingly, many of these cell type-specific monoallelically expressed genes were enriched for functions important for those brain cell types. In addition, function analysis showed that genes displaying monoallelic expression and correlated expression across neuronal cells from different individual brains were implicated in the regulation of synaptic function.

Conclusions

Our findings suggest that monoallelic gene expression is prevalent in human brain cells, which may play a role in generating cellular identity and neuronal diversity and thus increasing the complexity and diversity of brain cell functions.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4261-x) contains supplementary material, which is available to authorized users.

SIRT1 activates the expression of fetal hemoglobin genes

High fetal hemoglobin (HbF, α₂ γ₂ ) levels ameliorate the clinical manifestations of sickle cell disease and β thalassemia. The mechanisms that repress HbF expression and silence γ-globin genes in adults are incompletely characterized and only a single HbF inducer, hydroxyurea, is approved for treatment, and only in patients with sickle cell disease. We identified SIRT1, a protein deacetylase, as a new inducer of γ-globin. SIRT1 knockdown decreased, while SIRT1 ectopic expression upregulated γ-globin gene (HBG) expression in primary human erythroid cells and in K562 cells. The small molecule SIRT1 activators SRT2104 and SRT1720 enhanced HBG expression in cord blood human erythroblasts and reactivated silenced HBG in adult human erythroblasts. Furthermore, SIRT1 binds in the β-globin gene cluster locus control region (LCR) and HBG promoters, promotes the looping of the LCR to HBG promoter, and increases the binding of RNA polymerase II and H4K16Ac in the HBG promoter. SIRT1 suppressed the expression of the HBG suppressors BCL11A, KLF1, HDAC1 and HDAC2. Lastly, SIRT1 did not change the proliferation of human erythroid progenitor cells or the expression of differentiation marker CD235a. These data suggest that SIRT1 activates HBG expression through facilitating LCR looping to the HBG promoter, inhibiting the expression of transcriptional suppressors of HBG, and indirectly increasing histone acetylation in the HBG promoter. SIRT1 is a potential therapeutic target for γ-globin gene induction, and small molecule SIRT1 activators might serve as a lead compound for the development of new HbF inducers.

How many roads lead to cohesinopathies?

Genetic mapping studies reveal that mutations in cohesion pathways are responsible for multispectrum developmental abnormalities termed cohesinopathies. These include Roberts syndrome (RBS), Cornelia de Lange Syndrome (CdLS), and Warsaw Breakage Syndrome (WABS). The cohesinopathies are characterized by overlapping phenotypes ranging from craniofacial deformities, limb defects, and mental retardation. Though these syndromes share a similar suite of phenotypes and arise due to mutations in a common cohesion pathway, the underlying mechanisms are currently believed to be distinct. Defects in mitotic failure and apoptosis i.e. trans DNA tethering events are believed to be the underlying cause of RBS, whereas the underlying cause of CdLS is largely modeled as occurring through defects in transcriptional processes i.e. cis DNA tethering events. Here, we review recent findings described primarily in zebrafish, paired with additional studies in other model systems, including human patient cells, which challenge the notion that cohesinopathies represent separate syndromes. We highlight numerous studies that illustrate the utility of zebrafish to provide novel insights into the phenotypes, genes affected and the possible mechanisms underlying cohesinopathies. We propose that transcriptional deregulation is the predominant mechanism through which cohesinopathies arise. Developmental Dynamics 246:881-888, 2017. © 2017 Wiley Periodicals, Inc.

Comparative analysis of three-dimensional chromosomal architecture identifies a novel fetal hemoglobin regulatory element

In this study, Huang et al. compared the chromosomal architectures of fetal and adult human erythroblasts and found that, globally, chromatin structures and compartments A/B are highly similar at both developmental stages. Their results uncover a new critical regulatory region as a potential target for therapeutic genome editing for hemoglobinopathies and highlight the power of chromosome conformation analysis in discovering new cis control elements.

Chromatin structure is tightly intertwined with transcription regulation. Here we compared the chromosomal architectures of fetal and adult human erythroblasts and found that, globally, chromatin structures and compartments A/B are highly similar at both developmental stages. At a finer scale, we detected distinct folding patterns at the developmentally controlled β-globin locus. Specifically, new fetal stage-specific contacts were uncovered between a region separating the fetal (γ) and adult (δ and β) globin genes (encompassing the HBBP1 and BGLT3 noncoding genes) and two distal chromosomal sites (HS5 and 3′HS1) that flank the locus. In contrast, in adult cells, the HBBP1–BGLT3 region contacts the embryonic ε-globin gene, physically separating the fetal globin genes from the enhancer (locus control region [LCR]). Deletion of the HBBP1 region in adult cells alters contact landscapes in ways more closely resembling those of fetal cells, including increased LCR–γ-globin contacts. These changes are accompanied by strong increases in γ-globin transcription. Notably, the effects of HBBP1 removal on chromatin architecture and gene expression closely mimic those of deleting the fetal globin repressor BCL11A, implicating BCL11A in the function of the HBBP1 region. Our results uncover a new critical regulatory region as a potential target for therapeutic genome editing for hemoglobinopathies and highlight the power of chromosome conformation analysis in discovering new cis control elements.

Analyses of germline variants associated with ovarian cancer survival identify functional candidates at the 1q22 and 19p12 outcome loci

We previously identified associations with ovarian cancer outcome at five genetic loci. To identify putatively causal genetic variants and target genes, we prioritized two ovarian outcome loci (1q22 and 19p12) for further study. Bioinformatic and functional genetic analyses indicated that MEF2D and ZNF100 are targets of candidate outcome variants at 1q22 and 19p12, respectively. At 19p12, the chromatin interaction of a putative regulatory element with the ZNF100 promoter region correlated with candidate outcome variants. At 1q22, putative regulatory elements enhanced MEF2D promoter activity and haplotypes containing candidate outcome variants modulated these effects. In a public dataset, MEF2D and ZNF100 expression were both associated with ovarian cancer progression-free or overall survival time. In an extended set of 6,162 epithelial ovarian cancer patients, we found that functional candidates at the 1q22 and 19p12 loci, as well as other regional variants, were nominally associated with patient outcome; however, no associations reached our threshold for statistical significance (p<1×10-5). Larger patient numbers will be needed to convincingly identify any true associations at these loci.

Regulatory dynamics of 11p13 suggest a role for EHF in modifying CF lung disease severity

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene cause cystic fibrosis (CF), but are not good predictors of lung phenotype. Genome-wide association studies (GWAS) previously identified additional genomic sites associated with CF lung disease severity. One of these, at chromosome 11p13, is an intergenic region between Ets homologous factor (EHF) and Apaf-1 interacting protein (APIP). Our goal was to determine the functional significance of this region, which being intergenic is probably regulatory. To identify cis-acting elements, we used DNase-seq and H3K4me1 and H3K27Ac ChIP-seq to map open and active chromatin respectively, in lung epithelial cells. Two elements showed strong enhancer activity for the promoters of EHF and the 5' adjacent gene E47 like ETS transcription factor 5 (ELF5) in reporter gene assays. No enhancers of the APIP promoter were found. Circular chromosome conformation capture (4C-seq) identified direct physical interactions of elements within 11p13. This confirmed the enhancer-promoter associations, identified additional interacting elements and defined topologically associating domain (TAD) boundaries, enriched for CCCTC-binding factor (CTCF). No strong interactions were observed with the APIP promoter, which lies outside the main TAD encompassing the GWAS signal. These results focus attention on the role of EHF in modifying CF lung disease severity.

Spatial organization of genome architecture in neuronal development and disease

Although mammalian genomes encode genetic information in their linear sequences, their fundamental function with regard to gene expression depends on the higher-order structure of chromosomes. Current techniques for the evaluation of chromosomal structure have revealed that genomes are arranged at several hierarchical levels in three-dimensional space. The spatial organization of genomes involves the formation of chromatin loops that bypass a wide range of genomic distances, providing a connection between enhancers and chromosomal domains. Furthermore, they form chromatin domains that are arranged into chromosome territories in the three-dimensional space of the cell nucleus. Recent studies have shown that the spatial organization of the genome is essential for normal brain development and function. Activity-dependent alterations in the spatial organization of the genome can regulate transcriptional activity related to neuronal plasticity. Disruptions in the higher-order chromatin architecture have been implicated in neuropsychiatric disorders, such as cognitive dysfunction and anxiety. Here, we discuss the growing interest in the role of genome organization in brain development and neurological disorders.

Simulating topological domains in human chromosomes with a fitting-free model

We discuss a polymer model for the 3D organization of human chromosomes. A chromosome is represented by a string of beads, with each bead being "colored" according to 1D bioinformatic data (e.g., chromatin state, histone modification, GC content). Individual spheres (representing bi- and multi-valent transcription factors) can bind reversibly and selectively to beads with the appropriate color. During molecular dynamics simulations, the factors bind, and the string spontaneously folds into loops, rosettes, and topologically-associating domains (TADs). This organization occurs in the absence of any specified interactions between distant DNA segments, or between transcription factors. A comparison with Hi-C data shows that simulations predict the location of most boundaries between TADs correctly. The model is "fitting-free" in the sense that it does not use Hi-C data as an input; consequently, one of its strengths is that it can - in principle - be used to predict the 3D organization of any region of interest, or whole chromosome, in a given organism, or cell line, in the absence of existing Hi-C data. We discuss how this simple model might be refined to include more transcription factors and binding sites, and to correctly predict contacts between convergent CTCF binding sites.

Erythroid specific activator GATA-1-dependent interactions between CTCF sites around the β-globin locus

CTCF sites (binding motifs for CCCTC-binding factor, an insulator protein) are located considerable distances apart on genomes but are closely positioned in organized chromatin. The close positioning of CTCF sites is often cell type or tissue specific. Here we analyzed chromatin organization in eight CTCF sites around the β-globin locus by 3C assay and explored the roles of erythroid specific transcription activator GATA-1 and KLF1 in it. It was found five CTCF sites convergent to the locus interact with each other in erythroid K562 cells but not in non-erythroid 293 cells. The interaction was decreased by depletion of GATA-1 or KLF1. It accompanied reductions of CTCF and Rad21 occupancies and loss of active chromatin structure at the CTCF sites. Furthermore Rad21 occupancy was reduced in the β-globin locus control region (LCR) hypersensitive sites (HSs) by the depletion of GATA-1 or KLF1. The role of GATA-1 in interaction between CTCF sites was revealed by its ectopic expression in 293 cells and by deletion of a GATA-1 site in the LCR HS2. These findings indicate that erythroid specific activator GATA-1 acts at CTCF sites around the β-globin locus to establish tissue-specific chromatin organization.

CTCF interacts with the lytic HSV-1 genome to promote viral transcription

CTCF is an essential chromatin regulator implicated in important nuclear processes including in nuclear organization and transcription. Herpes Simplex Virus-1 (HSV-1) is a ubiquitous human pathogen, which enters productive infection in human epithelial and many other cell types. CTCF is known to bind several sites in the HSV-1 genome during latency and reactivation, but its function has not been defined. Here, we report that CTCF interacts extensively with the HSV-1 DNA during lytic infection by ChIP-seq, and its knockdown results in the reduction of viral transcription, viral genome copy number and virus yield. CTCF knockdown led to increased H3K9me3 and H3K27me3, and a reduction of RNA pol II occupancy on viral genes. Importantly, ChIP-seq analysis revealed that there is a higher level of CTD Ser2P modified RNA Pol II near CTCF peaks relative to the Ser5P form in the viral genome. Consistent with this, CTCF knockdown reduced the Ser2P but increased Ser5P modified forms of RNA Pol II on viral genes. These results suggest that CTCF promotes HSV-1 lytic transcription by facilitating the elongation of RNA Pol II and preventing silenced chromatin on the viral genome.

In the loop: how chromatin topology links genome structure to function in mechanisms underlying learning and memory

Different aspects of learning, memory, and cognition are regulated by epigenetic mechanisms such as covalent DNA modifications and histone post-translational modifications. More recently, the modulation of chromatin architecture and nuclear organization is emerging as a key factor in dynamic transcriptional regulation of the post-mitotic neuron. For instance, neuronal activity induces relocalization of gene loci to 'transcription factories', and specific enhancer-promoter looping contacts allow for precise transcriptional regulation. Moreover, neuronal activity-dependent DNA double-strand break formation in the promoter of immediate early genes appears to overcome topological constraints on transcription. Together, these findings point to a critical role for genome topology in integrating dynamic environmental signals to define precise spatiotemporal gene expression programs supporting cognitive processes.

High-Resolution 4C Reveals Rapid p53-Dependent Chromatin Reorganization of the CDKN1A Locus in Response to Stress

A regulatory program involving hundreds of genes is coordinated by p53 to prevent carcinogenesis in response to stress. Given the importance of chromatin loops in gene regulation, we investigated whether DNA interactions participate in the p53 stress response. To shed light on this issue, we measured the binding dynamics of cohesin in response to stress. We reveal that cohesin is remodeled at specific loci during the stress response and that its binding within genes negatively correlates with transcription. At p53 target genes, stress-induced eviction of cohesin from gene bodies is concomitant to spatial reorganization of loci through the disruption of functional chromatin loops. These findings demonstrate that chromatin loops can be remodeled upon stress and contribute to the p53-driven stress response. Additionally, we also propose a mechanism whereby transcription-coupled eviction of cohesin from CDKN1A might act as a molecular switch to control spatial interactions between regulatory elements.

Transcription factors mediate condensin recruitment and global chromosomal organization in fission yeast

It is becoming clear that structural-maintenance-of-chromosomes (SMC) complexes such as condensin and cohesin are involved in three-dimensional genome organization, yet their exact roles in functional organization remain unclear. We used chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) to comprehensively identify genome-wide associations mediated by condensin and cohesin in fission yeast. We found that although cohesin and condensin often bind to the same loci, they direct different association networks and generate small and larger chromatin domains, respectively. Cohesin mediates associations between loci positioned within 100 kb of each other; condensin can drive longer-range associations. Moreover, condensin, but not cohesin, connects cell cycle-regulated genes bound by mitotic transcription factors. This study describes the different functions of condensin and cohesin in genome organization and how specific transcription factors function in condensin loading, cell cycle-dependent genome organization and mitotic chromosome organization to support faithful chromosome segregation.

Dynamic cohesin-mediated chromatin architecture controls epithelial-mesenchymal plasticity in cancer

Epithelial to mesenchymal transition (EMT) and mesenchymal to epithelial transition (MET) are important interconnected events in tumorigenesis controlled by complex genetic networks. However, the cues that activate EMT-initiating factors and the mechanisms that reversibly connect EMT/MET are not well understood. Here, we show that cohesin-mediated chromatin organization coordinates EMT/MET by regulating mesenchymal genes. We report that RAD21, a subunit of the cohesin complex, is expressed in epithelial breast cancer cells, whereas its expression is decreased in mesenchymal cancer. Depletion of RAD21 in epithelial cancer cells causes transcriptional activation of TGFB1 and ITGA5, inducing EMT. Reduced binding of RAD21 changes intrachromosomal chromatin interactions within the TGFB1 and ITGA5 loci, creating an active transcriptional environment. Similarly, stem cell-like cancer cells also show an open chromatin structure at both genes, which correlates with high expression levels and mesenchymal fate characteristics. Conversely, overexpression of RAD21 in mesenchymal cancer cells induces MET-specific expression patterns. These findings indicate that dynamic cohesin-mediated chromatin structures are responsible for the initiation and regulation of essential EMT-related cell fate changes in cancer.

The Hierarchy of Transcriptional Activation: From Enhancer to Promoter

Regulatory elements (enhancers) that are remote from promoters play a critical role in the spatial, temporal, and physiological control of gene expression. Studies on specific loci, together with genome-wide approaches, suggest that there may be many common mechanisms involved in enhancer-promoter communication. Here, we discuss the multiprotein complexes that are recruited to enhancers and the hierarchy of events taking place between regulatory elements and promoters.

Entropy gives rise to topologically associating domains

We investigate chromosome organization within the nucleus using polymer models whose formulation is closely guided by experiments in live yeast cells. We employ bead-spring chromosome models together with loop formation within the chains and the presence of nuclear bodies to quantify the extent to which these mechanisms shape the topological landscape in the interphase nucleus. By investigating the genome as a dynamical system, we show that domains of high chromosomal interactions can arise solely from the polymeric nature of the chromosome arms due to entropic interactions and nuclear confinement. In this view, the role of bio-chemical related processes is to modulate and extend the duration of the interacting domains.

The Chromatin Signature of Pluripotency: Establishment and Maintenance

The revolutionary discovery that somatic cells can be reprogrammed by a defined set transcription factors to induced pluripotent stem cells (iPSCs) changed dramatically the way we perceive cell fate determination. Importantly, iPSCs, similar to embryo-derived stem cells (ESCs), are characterized by a remarkable developmental plasticity and the capacity to self-renew "indefinitely" under appropriate culture conditions, opening new avenues for personalized therapy and disease modeling. Elucidating the molecular mechanisms that maintain, induce, or alter stem cell identity is crucial for a deeper understanding of cell fate determination and potential translational applications. Intense research over the last 10 years exploiting technological advances in epigenomics and genome editing has unraveled many of the mysteries of pluripotent identity enabling novel and efficient ways to manipulate it for biomedical purposes. In this review, we focus on the chromatin and epigenetic characteristics that distinguish stem cells from somatic cells and their dynamic changes during differentiation and reprogramming.

Intergenic GWAS SNPs are key components of the spatial and regulatory network for human growth

Meta-analysis of genome-wide association studies has resulted in the identification of hundreds of genetic variants associated with growth and stature. Determining how these genetic variants influence growth is important, but most are non-coding, and there is little understanding of how these variants contribute to adult height. To determine the mechanisms by which human variation contributes to growth, we combined spatial genomic connectivity (high-throughput conformation capture) with functional (gene expression, expression Quantitative Trait Loci) data to determine how non-genic loci associated with infant length, pubertal and adult height and contribute to gene regulatory networks. This approach identified intergenic single-nucleotide polymorphisms (SNPs) ∼85 kb upstream of FBXW11 that spatially connect with distant loci. These regulatory connections are reinforced by evidence of SNP-enhancer effects and altered expression in genes influencing the action of human growth hormone. Functional assays provided evidence for enhancer activity of the intergenic region near FBXW11 that harbors SNP rs12153391, which is associated with an expression Quantitative Trait Loci. Our results suggest that variants in this locus have genome-wide effects as key modifiers of growth (both overgrowth and short stature) acting through a regulatory network. We believe that the genes and pathways connected with this regulatory network are potential targets that could be investigated for diagnostic, prenatal and carrier testing for growth disorders. Finally, the regulatory networks we generated illustrate the power of using existing datasets to interrogate the contribution of intergenic SNPs to common syndromes/diseases.

CTCF and CohesinSA-1 Mark Active Promoters and Boundaries of Repressive Chromatin Domains in Primary Human Erythroid Cells

Background

CTCF and cohesin^SA-1 are regulatory proteins involved in a number of critical cellular processes including transcription, maintenance of chromatin domain architecture, and insulator function. To assess changes in the CTCF and cohesin^SA-1 interactomes during erythropoiesis, chromatin immunoprecipitation coupled with high throughput sequencing and mRNA transcriptome analyses via RNA-seq were performed in primary human hematopoietic stem and progenitor cells (HSPC) and primary human erythroid cells from single donors.

Results

Sites of CTCF and cohesin^SA-1 co-occupancy were enriched in gene promoters in HSPC and erythroid cells compared to single CTCF or cohesin sites. Cell type-specific CTCF sites in erythroid cells were linked to highly expressed genes, with the opposite pattern observed in HSPCs. Chromatin domains were identified by ChIP-seq with antibodies against trimethylated lysine 27 histone H3, a modification associated with repressive chromatin. Repressive chromatin domains increased in both number and size during hematopoiesis, with many more repressive domains in erythroid cells than HSPCs. CTCF and cohesin^SA-1 marked the boundaries of these repressive chromatin domains in a cell-type specific manner.

Conclusion

These genome wide data, changes in sites of protein occupancy, chromatin architecture, and related gene expression, support the hypothesis that CTCF and cohesin^SA-1 have multiple roles in the regulation of gene expression during erythropoiesis including transcriptional regulation at gene promoters and maintenance of chromatin architecture. These data from primary human erythroid cells provide a resource for studies of normal and perturbed erythropoiesis.

Decoding transcriptional enhancers: Evolving from annotation to functional interpretation

Deciphering the intricate molecular processes that orchestrate the spatial and temporal regulation of genes has become an increasingly major focus of biological research. The differential expression of genes by diverse cell types with a common genome is a hallmark of complex cellular functions, as well as the basis for multicellular life. Importantly, a more coherent understanding of gene regulation is critical for defining developmental processes, evolutionary principles and disease etiologies. Here we present our current understanding of gene regulation by focusing on the role of enhancer elements in these complex processes. Although functional genomic methods have provided considerable advances to our understanding of gene regulation, these assays, which are usually performed on a genome-wide scale, typically provide correlative observations that lack functional interpretation. Recent innovations in genome editing technologies have placed gene regulatory studies at an exciting crossroads, as systematic, functional evaluation of enhancers and other transcriptional regulatory elements can now be performed in a coordinated, high-throughput manner across the entire genome. This review provides insights on transcriptional enhancer function, their role in development and disease, and catalogues experimental tools commonly used to study these elements. Additionally, we discuss the crucial role of novel techniques in deciphering the complex gene regulatory landscape and how these studies will shape future research.