NCAPD3 promotes diffuse large B-cell lymphoma progression through modulating SIRT1 expression in an H3K9 monomethylation-dependent manner

INTRODUCTION

Condensin, a family of structural maintenance of chromosome complexes, has been shown to regulate chromosome compaction and segregation during mitosis. NCAPD3, a HEAT-repeat subunit of condensin II, plays a dominant role in condensin-mediated chromosome dynamics but remains unexplored in lymphoma.

OBJECTIVES

The study aims to unravel the molecular function and mechanism of NCAPD3 in diffuse large B-cell lymphoma (DLBCL).

METHODS

The expression and clinical significance of NCAPD3 were assessed in public database and clinical specimens. Chromosome spreads, co-immunoprecipitation (co-IP), mass spectrometry (MS), and chromatin immunoprecipitation (ChIP) assays were conducted to untangle the role and mechanism of NCAPD3 in DLBCL.

RESULTS

NCAPD3 was highly expressed in DLBCL, correlated with poor prognosis. NCAPD3 deficiency impeded cell proliferation, induced apoptosis and increased the chemosensitivity. Instead, NCAPD3 overexpression facilitated cell proliferation. In vivo experiments further indicated targeting NCAPD3 suppressed tumor growth. Noteworthily, NCAPD3 deficiency disturbed the mitosis, triggering the formation of aneuploids. To reveal the function of NCAPD3 in DLBCL, chromosome spreads were conducted, presenting that chromosomes became compact upon NCAPD3 overexpression, instead, loose, twisted and lacking axial rigidity upon NCAPD3 absence. Meanwhile, the classical transcription-activated marker, H3K4 trimethylation, was found globally upregulated after NCAPD3 knockout, suggesting that NCAPD3 might participate in chromatin remodeling and transcription regulation. MS revealed NCAPD3 could interact with transcription factor, TFII I. Further co-IP and ChIP assays verified NCAPD3 could be anchored at the promoter of SIRT1 by TFII I and then supported the transcription of SIRT1 via recognizing H3K9 monomethylation (H3K9me1) on SIRT1 promoter. Function reversion assay verified the oncogenic role of NCAPD3 in DLBCL was partially mediated by SIRT1.

CONCLUSION

This study demonstrated that dysregulation of NCAPD3 could disturb chromosome compaction and segregation and regulate the transcription activity of SIRT1 in an H3K9me1-dependent manner, which provided novel insights into targeted strategy for DLBCL.

Variant to gene mapping for carpal tunnel syndrome risk loci implicates skeletal muscle regulatory elements

BACKGROUND

Carpal tunnel syndrome (CTS) is a common disorder caused by compression of the median nerve in the wrist, resulting in pain and numbness throughout the hand and forearm. While multiple behavioural and physiological factors influence CTS risk, a growing body of evidence supports a strong genetic contribution. Recent genome-wide association study (GWAS) efforts have reported 53 independent signals associated with CTS. While GWAS can identify genetic loci conferring risk, it does not determine which cell types drive the genetic aetiology of the trait, which variants are "causal" at a given signal, and which effector genes correspond to these non-coding variants. These obstacles limit interpretation of potential disease mechanisms.

METHODS

We analysed CTS GWAS findings in the context of chromatin conformation between gene promoters and accessible chromatin regions across cellular models of bone, skeletal muscle, adipocytes and neurons. We identified proxy variants in high LD with the lead CTS sentinel SNPs residing in promoter connected open chromatin in the skeletal muscle and bone contexts.

FINDINGS

We detected significant enrichment for heritability in skeletal muscle myotubes, as well as a weaker correlation in human mesenchymal stem cell-derived osteoblasts. In myotubes, our approach implicated 117 genes contacting 60 proxy variants corresponding to 20 of the 53 GWAS signals. In the osteoblast context we implicated 30 genes contacting 24 proxy variants coinciding with 12 signals, of which 19 genes shared. We subsequently prioritized BZW2 as a candidate effector gene in CTS and implicated it as novel gene that perturbs myocyte differentiation in vitro.

INTERPRETATION

Taken together our results suggest that the CTS genetic component influences the size, integrity, and organization of multiple tissues surrounding the carpal tunnel, in particular muscle and bone, to predispose the nerve to being compressed in this disease setting.

FUNDING

This work was supported by NIH Grant UM1 DK126194 (SFAG and WY), R01AG072705 (SFAG & KDH) and the Center for Spatial and Functional Genomics at CHOP (SFAG & ADW). SFAG is supported by the Daniel B. Burke Endowed Chair for Diabetes Research. WY is supported by the Perelman School of Medicine of the University of Pennsylvania.

Co-localization of clusters of TCR-regulated genes with TAD rearrangements

BACKGROUND

Gene expression has long been known to be influenced by the relative proximity of DNA regulatory elements. Topologically associating domains (TADs) are self-interacting genomic regions involved in regulating gene expression by controlling the proximity of these elements. Prior studies of TADs and their biological roles have revealed correlations between TAD changes and cellular differentiation. Here, we used Hi-C and RNA-seq data to correlate TCR-induced changes in TAD structure and gene expression in human CD4+ T cells.

RESULTS

We developed a pipeline, Differentially Expressed Gene Enrichment Finder (DEGEF), that identifies regions of differentially expressed gene enrichment. Using DEGEF, we found that TCR-regulated genes cluster non-uniformly across the genome and that these clusters preferentially localized in regions of TAD rearrangement. Interestingly, clusters of upregulated genes preferentially formed new Hi-C contacts compared to downregulated clusters, suggesting that TCR-activated CD4+ T cells may regulate genes by changing stimulatory contacts rather than inhibitory contacts.

CONCLUSIONS

Our observations support a significant relationship between TAD rearrangements and changes in local gene expression. These findings indicate potentially important roles for TAD rearrangements in shaping their local regulatory environments and thus driving differential expression of nearby genes during CD4+ T cell activation. Moreover, they provide new insights into global mechanisms that regulate gene expression.

Dynamic chromatin conformation and accessibility changes mediate the spatial-specific gene regulatory network in Bombyx mori

Silk gland genes of Bombyx mori can have strict spatial expression patterns, which impact their functions and silk quality; however, our understanding of their regulation mechanisms is currently insufficient. To address this, the middle silk gland (MSG) and posterior silk gland (PSG) of the silkworm were investigated. Gene ontology annotation showed that spatially specific expressed genes were involved in the formation of H3k9me and chromatin topology. Chromatin conformation data generated by Hi-C showed that the topologically associated domain boundaries around FibL and Sericin1 genes were significantly different between MSG and PSG. Changes in chromatin conformation led to changes in chromatin activity, which significantly affected the expression of nearby genes in silkworm. Chromatin accessibility regions of MSG and PSG were analyzed using FAIRE-seq, and 1006 transcription factor motifs were identified in open chromatin regions. Furthermore, the spatial-specific expression patterns of silk gland genes were mainly associated with homeobox-contained transcription factors, such as POU-M2, which was specifically bound and relatively highly expressed in the MSG. The regulatory network mediated by POU-M2 regulated most of the spatial-specific expressed genes in MSG, such as ADH1. These results can aid in improving silk performance, optimizing silkworm breeding, and improving the gene spatial regulatory model research for insects.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

Mini-review: Gene regulatory network benefits from three-dimensional chromatin conformation and structural biology

Gene regulatory networks are now at the forefront of precision biology, which can help researchers better understand how genes and regulatory elements interact to control cellular gene expression, offering a more promising molecular mechanism in biological research. Interactions between the genes and regulatory elements involve different promoters, enhancers, transcription factors, silencers, insulators, and long-range regulatory elements, which occur at a ∼10 µm nucleus in a spatiotemporal manner. In this way, three-dimensional chromatin conformation and structural biology are critical for interpreting the biological effects and the gene regulatory networks. In the review, we have briefly summarized the latest processes in three-dimensional chromatin conformation, microscopic imaging, and bioinformatics, and we have presented the outlook and future directions for these three aspects.

BET bromodomain inhibitor JQ1 regulates spermatid development by changing chromatin conformation in mouse spermatogenesis

As a BET bromodomain inhibitor, JQ1 has been proven have efficacy against a number of different cancers. In terms of male reproduction, JQ1 may be used as a new type of contraceptive, since JQ1 treatment in male mice could lead to germ cell defects and a decrease of sperm motility, moreover, this effect is reversible. However, the mechanism of JQ1 acting on gene regulation in spermatogenesis remains unclear. Here, we performed single-cell RNA sequencing (scRNA-seq) on mouse testes treated with JQ1 or vehicle control to determine the transcriptional regulatory function of JQ1 in spermatogenesis at the single cell resolution. We confirmed that JQ1 treatment could increase the numbers of somatic cells and spermatocytes and decrease the numbers of spermatid cells. Gene Ontology (GO) analysis demonstrated that differentially expressed genes which were down-regulated after JQ1 injection were mainly enriched in "DNA conformation change" biological process in early developmental germ cells and "spermatid development" biological process in spermatid cells. ATAC-seq data further confirmed that JQ1 injection could change the open state of chromatin. In addition, JQ1 could change the numbers of accessible meiotic DNA double-stranded break sites and the types of transcription factor motif that functioned in pachytene spermatocytes and round spermatids. The multi-omics analysis revealed that JQ1 had the ability to regulate gene transcription by changing chromatin conformation in mouse spermatogenesis, which would potentiate the availability of JQ1 in male contraceptive.

Fine mapping with epigenetic information and 3D structure

Since 2005, thousands of genome-wide association studies (GWAS) have been published, identifying hundreds of thousands of genetic variants that increase risk of complex traits such as autoimmune diseases. This wealth of data has the potential to improve patient care, through personalized medicine and the identification of novel drug targets. However, the potential of GWAS for clinical translation has not been fully achieved yet, due to the fact that the functional interpretation of risk variants and the identification of causal variants and genes are challenging. The past decade has seen the development of great advances that are facilitating the overcoming of these limitations, by utilizing a plethora of genomics and epigenomics tools to map and characterize regulatory elements and chromatin interactions, which can be used to fine map GWAS loci, and advance our understanding of the biological mechanisms that cause disease.

Analysis of Hi-C Data for Discovery of Structural Variations in Cancer

Structural variations (SVs) are large genomic rearrangements that can be challenging to identify with current short read sequencing technology due to various confounding factors such as existence of genomic repeats and complex SV structures. Hi-C breakfinder is the first computational tool that utilizes the technology of high-throughput chromatin conformation capture assay (Hi-C) to systematically identify SVs, without being interfered by regular confounding factors. SVs change the spatial distance of genomic regions and cause discontinuous signals in Hi-C, which are difficult to analyze by routine informatics practice. Here we provide step-by-step guidance for how to identify SVs using Hi-C data and how to reconstruct Hi-C maps in the presence of SVs.

Monocytic THP-1 cells diverge significantly from their primary counterparts: a comparative examination of the chromosomal conformations and transcriptomes

Immortalized cell lines have long been used as model systems to systematically investigate biological processes under controlled and reproducible conditions, providing insights that have greatly advanced cellular biology and medical sciences. Recently, the widely used monocytic leukemia cell line, THP-1, was comprehensively examined to understand mechanistic relationships between the 3D chromatin structure and transcription during the trans-differentiation of monocytes to macrophages. To corroborate these observations in primary cells, we analyze in situ Hi-C and RNA-seq data of human primary monocytes and their differentiated macrophages in comparison to that obtained from the monocytic/macrophagic THP-1 cells. Surprisingly, we find significant differences between the primary cells and the THP-1 cells at all levels of chromatin structure, from loops to topologically associated domains to compartments. Importantly, the compartment-level differences correlate significantly with transcription: those genes that are in A-compartments in the primary cells but are in B-compartments in the THP-1 cells exhibit a higher level of expression in the primary cells than in the THP-1 cells, and vice versa. Overall, the genes in these different compartments are enriched for a wide range of pathways, and, at least in the case of the monocytic cells, their altered expression in certain pathways in the THP-1 cells argues for a less immune cell-like phenotype, suggesting that immortalization or prolonged culturing of THP-1 caused a divergence of these cells from their primary counterparts. It is thus essential to reexamine phenotypic details observed in cell lines with their primary counterparts so as to ensure a proper understanding of functional cell states in vivo.

Understanding three-dimensional chromatin organization in diploid genomes

The three-dimensional (3D) organization of chromatin in the nucleus of diploid eukaryotic organisms has fascinated biologists for many years. Despite major progress in chromatin conformation studies, current knowledge regarding the spatial organization of diploid (maternal and paternal) genomes is still limited. Recent advances in Hi-C technology and data processing approaches have enabled construction of diploid Hi-C contact maps. These maps greatly accelerated the pace of novel discoveries in haplotype-resolved 3D genome studies, revealing the role of allele biased chromatin conformation in transcriptional regulation. Here, we review emerging concepts and haplotype phasing strategies of Hi-C data in 3D diploid genome studies. We discuss new insights on homologous chromosomal organization and the interplay between allelic biased chromatin architecture and several nuclear functions, explaining how haplotype-resolved Hi-C technologies have been used to resolve important biological questions.

Genomic properties of variably methylated retrotransposons in mouse

BACKGROUND

Transposable elements (TEs) are enriched in cytosine methylation, preventing their mobility within the genome. We previously identified a genome-wide repertoire of candidate intracisternal A particle (IAP) TEs in mice that exhibit inter-individual variability in this methylation (VM-IAPs) with implications for genome function.

RESULTS

Here we validate these metastable epialleles and discover a novel class that exhibit tissue specificity (tsVM-IAPs) in addition to those with uniform methylation in all tissues (constitutive- or cVM-IAPs); both types have the potential to regulate genes in cis. Screening for variable methylation at other TEs shows that this phenomenon is largely limited to IAPs, which are amongst the youngest and most active endogenous retroviruses. We identify sequences enriched within cVM-IAPs, but determine that these are not sufficient to confer epigenetic variability. CTCF is enriched at VM-IAPs with binding inversely correlated with DNA methylation. We uncover dynamic physical interactions between cVM-IAPs with low methylation ranges and other genomic loci, suggesting that VM-IAPs have the potential for long-range regulation.

CONCLUSION

Our findings indicate that a recently evolved interplay between genetic sequence, CTCF binding, and DNA methylation at young TEs can result in inter-individual variability in transcriptional outcomes with implications for phenotypic variation.

Straightforward protocol for allele-specific chromatin conformation capture

A key advance in our understanding of gene regulation came with the finding that the genome undergoes three-dimensional nuclear folding in a genetically determined process. This 3D conformation directly influences the association between enhancers and their target promoters. This complex interplay has been proven to be essential for gene regulation, and genetic variants affecting this process have been associated to human diseases. The development of new technologies that quantify these DNA interactions represented a revolution in the field. High throughput techniques like HiC provide a general picture of chromatin topology. However, they often lack resolution to evidence subtle effects that single nucleotide polymorphisms exert over the contacts between cis-regulatory regions and target promoters. Here we propose a cost-efficient approach to perform allele-specific chromatin conformation analysis. As a proof of concept, we analyzed the impact of a common deletion mapping between SIRPB1 promoter and one of its downstream enhancers.

The matrix revolutions: towards the decoding of the plant chromatin three-dimensional reality

In recent years, we have witnessed a significant increase in studies addressing the three-dimensional (3D) chromatin organization of the plant nucleus. Important advances in chromatin conformation capture (3C)-derived and related techniques have allowed the exploration of the nuclear topology of plants with large and complex genomes, including various crops. In addition, the increase in their resolution has permitted the depiction of chromatin compartmentalization and interactions at the gene scale. These studies have revealed the highly complex mechanisms governing plant nuclear architecture and the remarkable knowledge gaps in this field. Here we discuss the state-of-the-art in plant chromosome architecture, including our knowledge of the hierarchical organization of the genome in 3D space and regarding other nuclear components. Furthermore, we highlight the existence in plants of topologically associated domain (TAD)-like structures that display striking differences from their mammalian counterparts, proposing the concept of ICONS-intergenic condensed spacers. Similarly, we explore recent advances in the study of chromatin loops and R-loops, and their implication in the regulation of gene activity. Finally, we address the impact that polyploidization has had on the chromatin topology of modern crops, and how this is related to phenomena such as subgenome dominance and biased gene retention in these organisms.

TADsplimer reveals splits and mergers of topologically associating domains for epigenetic regulation of transcription

We present TADsplimer, the first computational tool to systematically detect topologically associating domain (TAD) splits and mergers across the genome between Hi-C samples. TADsplimer recaptures splits and mergers of TADs with high accuracy in simulation analyses and defines hundreds of TAD splits and mergers between pairs of different cell types, such as endothelial cells and fibroblasts. Our work reveals a key role for TAD remodeling in epigenetic regulation of transcription and delivers the first tool for the community to perform dynamic analysis of TAD splits and mergers in numerous biological and disease models.

Alu retrotransposons modulate Nanog expression through dynamic changes in regional chromatin conformation via aryl hydrocarbon receptor

Transcriptional repression of Nanog is an important hallmark of stem cell differentiation. Chromatin modifications have been linked to the epigenetic profile of the Nanog gene, but whether chromatin organization actually plays a causal role in Nanog regulation is still unclear. Here, we report that the formation of a chromatin loop in the Nanog locus is concomitant to its transcriptional downregulation during human NTERA-2 cell differentiation. We found that two Alu elements flanking the Nanog gene were bound by the aryl hydrocarbon receptor (AhR) and the insulator protein CTCF during cell differentiation. Such binding altered the profile of repressive histone modifications near Nanog likely leading to gene insulation through the formation of a chromatin loop between the two Alu elements. Using a dCAS9-guided proteomic screening, we found that interaction of the histone methyltransferase PRMT1 and the chromatin assembly factor CHAF1B with the Alu elements flanking Nanog was required for chromatin loop formation and Nanog repression. Therefore, our results uncover a chromatin-driven, retrotransposon-regulated mechanism for the control of Nanog expression during cell differentiation.

Integrative prediction of gene expression with chromatin accessibility and conformation data

BACKGROUND

Enhancers play a fundamental role in orchestrating cell state and development. Although several methods have been developed to identify enhancers, linking them to their target genes is still an open problem. Several theories have been proposed on the functional mechanisms of enhancers, which triggered the development of various methods to infer promoter-enhancer interactions (PEIs). The advancement of high-throughput techniques describing the three-dimensional organization of the chromatin, paved the way to pinpoint long-range PEIs. Here we investigated whether including PEIs in computational models for the prediction of gene expression improves performance and interpretability.

RESULTS

We have extended our [Formula: see text] framework to include DNA contacts deduced from chromatin conformation capture experiments and compared various methods to determine PEIs using predictive modelling of gene expression from chromatin accessibility data and predicted transcription factor (TF) motif data. We designed a novel machine learning approach that allows the prioritization of TFs binding to distal loop and promoter regions with respect to their importance for gene expression regulation. Our analysis revealed a set of core TFs that are part of enhancer-promoter loops involving YY1 in different cell lines.

CONCLUSION

We present a novel approach that can be used to prioritize TFs involved in distal and promoter-proximal regulatory events by integrating chromatin accessibility, conformation, and gene expression data. We show that the integration of chromatin conformation data can improve gene expression prediction and aids model interpretability.

DNase I Chromatin Accessibility Analysis

Chromatin consists of compacted DNA in complex with proteins and contributes to the organization of DNA and its stability. Furthermore, chromatin plays key roles in regulating cellular processes such as DNA replication, transcription, DNA repair, and mitosis. Chromatin assumes more compact (inaccessible) or decondensed (accessible) conformations depending on the function that is being supported in the genome, either locally or globally. The activity of nucleases has been used previously to assess the accessibility of specific genomic regions in vitro, such as origins of replication at varying points in the cell cycle. Here, we provide an assay to determine the accessibility of specific human genomic regions (example used herein: Lamin B2 origin of DNA replication) by measuring the effect of DNase I nuclease on qPCR signal from the studied site. This assay provides a powerful method to interrogate the molecular mechanisms that regulate chromatin accessibility, and how these processes affect various cellular functions involving the human genome that require manipulation of chromatin conformation.

Visualising three-dimensional genome organisation in two dimensions

The three-dimensional organisation of the genome plays a crucial role in developmental gene regulation. In recent years, techniques to investigate this organisation have become more accessible to labs worldwide due to improvements in protocols and decreases in the cost of high-throughput sequencing. However, the resulting datasets are complex and can be challenging to analyse and interpret. Here, we provide a guide to visualisation approaches that can aid the interpretation of such datasets and the communication of biological results.

Fine mapping chromatin contacts in capture Hi-C data

Background

Hi-C and capture Hi-C (CHi-C) are used to map physical contacts between chromatin regions in cell nuclei using high-throughput sequencing. Analysis typically proceeds considering the evidence for contacts between each possible pair of fragments independent from other pairs. This can produce long runs of fragments which appear to all make contact with the same baited fragment of interest.

Results

We hypothesised that these long runs could result from a smaller subset of direct contacts and propose a new method, based on a Bayesian sparse variable selection approach, which attempts to fine map these direct contacts. Our model is conceptually novel, exploiting the spatial pattern of counts in CHi-C data. Although we use only the CHi-C count data in fitting the model, we show that the fragments prioritised display biological properties that would be expected of true contacts: for bait fragments corresponding to gene promoters, we identify contact fragments with active chromatin and contacts that correspond to edges found in previously defined enhancer-target networks; conversely, for intergenic bait fragments, we identify contact fragments corresponding to promoters for genes expressed in that cell type. We show that long runs of apparently co-contacting fragments can typically be explained using a subset of direct contacts consisting of <10% of the number in the full run, suggesting that greater resolution can be extracted from existing datasets.

Conclusions

Our results appear largely complementary to those from a per-fragment analytical approach, suggesting that they provide an additional level of interpretation that may be used to increase resolution for mapping direct contacts in CHi-C experiments.

Electronic supplementary material

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

Migration through a small pore disrupts inactive chromatin organization in neutrophil-like cells

BACKGROUND

Mammalian cells are flexible and can rapidly change shape when they contract, adhere, or migrate. The nucleus must be stiff enough to withstand cytoskeletal forces, but flexible enough to remodel as the cell changes shape. This is particularly important for cells migrating through confined spaces, where the nuclear shape must change in order to fit through a constriction. This occurs many times in the life cycle of a neutrophil, which must protect its chromatin from damage and disruption associated with migration. Here we characterized the effects of constricted migration in neutrophil-like cells.

RESULTS

Total RNA sequencing identified that migration of neutrophil-like cells through 5- or 14-μm pores was associated with changes in the transcript levels of inflammation and chemotaxis-related genes when compared to unmigrated cells. Differentially expressed transcripts specific to migration with constriction were enriched for groups of genes associated with cytoskeletal remodeling. Hi-C was used to capture the genome organization in control and migrated cells. Limited switching was observed between the active (A) and inactive (B) compartments after migration. However, global depletion of short-range contacts was observed following migration with constriction compared to migration without constriction. Regions with disrupted contacts, TADs, and compartments were enriched for inactive chromatin.

CONCLUSION

Short-range genome organization is preferentially altered in inactive chromatin, possibly protecting transcriptionally active contacts from the disruptive effects of migration with constriction. This is consistent with current hypotheses implicating heterochromatin as the mechanoresponsive form of chromatin. Further investigation concerning the contribution of heterochromatin to stiffness, flexibility, and protection of nuclear function will be important for understanding cell migration in relation to human health and disease.

From association to mechanism in complex disease genetics: the role of the 3D genome

Genome-wide association studies (GWAS) and fine mapping studies in autoimmune diseases have identified thousands of genetic variants, the majority of which are located in non-protein-coding enhancer regions. Enhancers function within the context of the three-dimensional (3D) genome to form long-range DNA looping events with target gene promoters that spatially and temporally regulate gene expression. Investigating the functional significance of GWAS variants in the context of the 3D genome is essential for mechanistic understanding of these variants and how they influence disease pathology by altering DNA looping between enhancers and the target gene promoters they regulate. In this review, we discuss the functional complexity of the 3D genome and the technological approaches used to characterize DNA looping events. We then highlight examples from the literature that illustrate how functional mapping of the 3D genome can assist in defining mechanisms that influence pathogenic gene expression. We conclude by highlighting future advances necessary to fully integrate 3D genome analyses into the functional workup of GWAS variants in the continuing effort to improve the health of patients with autoimmune diseases.

Control of directionality of chromatin folding for the inter- and intra-domain contacts at the Tfap2c-Bmp7 locus

BACKGROUND

Contact domains of chromatin serve as a fundamental unit to regulate action of enhancers for target genes. Looping between a pair of CCCTC-binding factor (CTCF)-binding sites in convergent orientations underlies the formation of contact domains, while those in divergent orientations establish domain boundaries. However, every CTCF site is not necessarily engaged in loop or boundary structures, leaving functions of CTCF in varied genomic contexts still elusive. The locus containing Tfap2c and Bmp7 encompasses two contact domains separated by a region between the two genes, termed transition zone (TZ), characterized by two arrays of CTCF sites in divergent configuration. In this study, we created deletion and inversion alleles of these and other regions across the locus and investigated how they impinge on the conformation.

RESULTS

Deletion of the whole two CTCF arrays with the CRISPR/Cas9 system resulted in impairment of blocking of chromatin contacts by the TZ, as assessed by the circular chromatin conformation capture assay (4C-seq). Deletion and inversion of either of the two arrays similarly, but less pronouncedly, led to reduction in the blocking activity. Thus, the divergent configuration provides the TZ with the strong boundary activity. Uniquely, we show the TZ harbors a 50-kb region within one of the two arrays that contacts broadly with the both flanking intervals, regardless of the presence or orientation of the other CTCF array. Further, we show the boundary CTCF array has little impact on intra-domain folding; instead, locally associating CTCF sites greatly affect it.

CONCLUSIONS

Our results show that the TZ not only separates the two domains, but also bears a wide interval that shows isotropic behavior of chromatin folding, indicating a potentially complex nature of actual boundaries in the genome. We also show that CTCF-binding sites inside a domain greatly contribute to the intra-domain folding of chromatin. Thus, the study reveals diverse and context-dependent roles of CTCF in organizing chromatin conformation at different levels.

Finding local genome rearrangements

Background

The double cut and join (DCJ) model of genome rearrangement is well studied due to its mathematical simplicity and power to account for the many events that transform gene order. These studies have mostly been devoted to the understanding of minimum length scenarios transforming one genome into another. In this paper we search instead for rearrangement scenarios that minimize the number of rearrangements whose breakpoints are unlikely due to some biological criteria. One such criterion has recently become accessible due to the advent of the Hi-C experiment, facilitating the study of 3D spacial distance between breakpoint regions.

Results

We establish a link between the minimum number of unlikely rearrangements required by a scenario and the problem of finding a maximum edge-disjoint cycle packing on a certain transformed version of the adjacency graph. This link leads to a 3/2-approximation as well as an exact integer linear programming formulation for our problem, which we prove to be NP-complete. We also present experimental results on fruit flies, showing that Hi-C data is informative when used as a criterion for rearrangements.

Conclusions

A new variant of the weighted DCJ distance problem is addressed that ignores scenario length in its objective function. A solution to this problem provides a lower bound on the number of unlikely moves necessary when transforming one gene order into another. This lower bound aids in the study of rearrangement scenarios with respect to chromatin structure, and could eventually be used in the design of a fixed parameter algorithm with a more general objective function.

Chromosome contacts in activated T cells identify autoimmune disease candidate genes

Background

Autoimmune disease-associated variants are preferentially found in regulatory regions in immune cells, particularly CD4⁺ T cells. Linking such regulatory regions to gene promoters in disease-relevant cell contexts facilitates identification of candidate disease genes.

Results

Within 4 h, activation of CD4⁺ T cells invokes changes in histone modifications and enhancer RNA transcription that correspond to altered expression of the interacting genes identified by promoter capture Hi-C. By integrating promoter capture Hi-C data with genetic associations for five autoimmune diseases, we prioritised 245 candidate genes with a median distance from peak signal to prioritised gene of 153 kb. Just under half (108/245) prioritised genes related to activation-sensitive interactions. This included IL2RA, where allele-specific expression analyses were consistent with its interaction-mediated regulation, illustrating the utility of the approach.

Conclusions

Our systematic experimental framework offers an alternative approach to candidate causal gene identification for variants with cell state-specific functional effects, with achievable sample sizes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1285-0) contains supplementary material, which is available to authorized users.

Heterogeneous pattern of DNA methylation in developmentally important genes correlates with its chromatin conformation

Background

DNA methylation is a major epigenetic modification, playing a crucial role in the development and differentiation of higher organisms. DNA methylation is also known to regulate transcription by gene repression. Various developmental genes such as c-mos, HoxB5, Sox11, and Sry show tissue-specific gene expression that was shown to be regulated by promoter DNA methylation. The aim of the present study is to investigate the establishment of chromatin marks (active or repressive) in relation to heterogeneous methylation in the promoter regions of these developmentally important genes.

Results

Chromatin-immunoprecipitation (ChIP) assays were performed to immuno-precipitate chromatin by antibodies against both active (H3K4me3) and repressive (H3K9me3) chromatin regions. The analysis of ChIP results showed that both the percentage input and fold enrichment of activated chromatin was higher in tissues expressing the respective genes as compared to the tissues not expressing the same set of genes. This was true for all the genes selected for the study (c-mos, HoxB5, Sox11, and Sry). These findings illustrate that inconsistent DNA methylation patterns (sporadic, mosaic and heterogeneous) may also influence gene regulation, thereby resulting in the modulation of chromatin conformation.

Conclusions

These findings illustrate that various patterns of DNA methylation (asynchronous, mosaic and heterogeneous) correlates with chromatin modification, resulting in the gene regulation.

Identification of a RAI1-associated disease network through integration of exome sequencing, transcriptomics, and 3D genomics

Background

Smith-Magenis syndrome (SMS) is a developmental disability/multiple congenital anomaly disorder resulting from haploinsufficiency of RAI1. It is characterized by distinctive facial features, brachydactyly, sleep disturbances, and stereotypic behaviors.

Methods

We investigated a cohort of 15 individuals with a clinical suspicion of SMS who showed neither deletion in the SMS critical region nor damaging variants in RAI1 using whole exome sequencing. A combination of network analysis (co-expression and biomedical text mining), transcriptomics, and circularized chromatin conformation capture (4C-seq) was applied to verify whether modified genes are part of the same disease network as known SMS-causing genes.

Results

Potentially deleterious variants were identified in nine of these individuals using whole-exome sequencing. Eight of these changes affect KMT2D, ZEB2, MAP2K2, GLDC, CASK, MECP2, KDM5C, and POGZ, known to be associated with Kabuki syndrome 1, Mowat-Wilson syndrome, cardiofaciocutaneous syndrome, glycine encephalopathy, mental retardation and microcephaly with pontine and cerebellar hypoplasia, X-linked mental retardation 13, X-linked mental retardation Claes-Jensen type, and White-Sutton syndrome, respectively. The ninth individual carries a de novo variant in JAKMIP1, a regulator of neuronal translation that was recently found deleted in a patient with autism spectrum disorder. Analyses of co-expression and biomedical text mining suggest that these pathologies and SMS are part of the same disease network. Further support for this hypothesis was obtained from transcriptome profiling that showed that the expression levels of both Zeb2 and Map2k2 are perturbed in Rai1^–/– mice. As an orthogonal approach to potentially contributory disease gene variants, we used chromatin conformation capture to reveal chromatin contacts between RAI1 and the loci flanking ZEB2 and GLDC, as well as between RAI1 and human orthologs of the genes that show perturbed expression in our Rai1^–/– mouse model.

Conclusions

These holistic studies of RAI1 and its interactions allow insights into SMS and other disorders associated with intellectual disability and behavioral abnormalities. Our findings support a pan-genomic approach to the molecular diagnosis of a distinctive disorder.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0359-z) contains supplementary material, which is available to authorized users.

Plant programmed cell death from a chromatin point of view

Programmed cell death (PCD) is a ubiquitous genetically regulated process consisting of the activation of finely controlled signalling pathways that lead to cellular suicide. PCD can be part of a developmental programme (dPCD) or be triggered by environmental conditions (ePCD). In plant cells, as in animal cells, extensive chromatin condensation and degradation of the nuclear DNA are among the most conspicuous features of cells undergoing PCD. Changes in chromatin condensation could either reflect the structural changes required for internucleosomal fragmentation of nuclear DNA or relate to large-scale chromatin rearrangements associated with a major transcriptional switch occurring during cell death. The aim of this review is to give an update on plant PCD processes from a chromatin point of view. The first part will be dedicated to chromatin conformational changes associated with cell death observed in various developmental and physiological conditions, whereas the second part will be devoted to histone dynamics and DNA modifications associated with critical changes in genome expression during the cell death process.

Roles and regulation of bone morphogenetic protein-7 in kidney development and diseases

The gene encoding bone morphogenetic protein-7 (BMP7) is expressed in the developing kidney in embryos and also in the mature organ in adults. During kidney development, expression of BMP7 is essential to determine the final number of nephrons in and proper size of the organ. The secreted BMP7 acts on the nephron progenitor cells to exert its dual functions: To maintain and expand the progenitor population and to provide them with competence to respond to differentiation cues, each relying on distinct signaling pathways. Intriguingly, in the adult organ, BMP7 has been implicated in protection against and regeneration from injury. Exogenous administration of recombinant BMP7 to animal models of kidney diseases has shown promising effects in counteracting inflammation, apoptosis and fibrosis evoked upon injury. Although the expression pattern of BMP7 has been well described, the mechanisms by which it is regulated have remained elusive and the processes by which the secretion sites of BMP7 impinge upon its functions in kidney development and diseases have not yet been assessed. Understanding the regulatory mechanisms will pave the way towards gaining better insight into the roles of BMP7, and to achieving desired control of the gene expression as a therapeutic strategy for kidney diseases.

Models and algorithms for genome rearrangement with positional constraints

Background

Traditionally, the merit of a rearrangement scenario between two gene orders has been measured based on a parsimony criteria alone; two scenarios with the same number of rearrangements are considered equally good. In this paper, we acknowledge that each rearrangement has a certain likelihood of occurring based on biological constraints, e.g. physical proximity of the DNA segments implicated or repetitive sequences.

Results

We propose optimization problems with the objective of maximizing overall likelihood, by weighting the rearrangements. We study a binary weight function suitable to the representation of sets of genome positions that are most likely to have swapped adjacencies. We give a polynomial-time algorithm for the problem of finding a minimum weight double cut and join scenario among all minimum length scenarios. In the process we solve an optimization problem on colored noncrossing partitions, which is a generalization of the Maximum Independent Set problem on circle graphs.

Conclusions

We introduce a model for weighting genome rearrangements and show that under simple yet reasonable conditions, a fundamental distance can be computed in polynomial time. This is achieved by solving a generalization of the Maximum Independent Set problem on circle graphs. Several variants of the problem are also mentioned.

Put your 3D glasses on: plant chromatin is on show

The three-dimensional organization of the eukaryotic nucleus and its chromosomal conformation have emerged as important features in the complex network of mechanisms behind gene activity and genome connectivity dynamics, which can be evidenced in the regionalized chromosomal spatial distribution and the clustering of diverse genomic regions with similar expression patterns. The development of chromatin conformation capture (3C) techniques has permitted the elucidation of commonalities between the eukaryotic phyla, as well as important differences among them. The growing number of studies in the field performed in plants has shed light on the structural and regulatory features of these organisms. For instance, it has been proposed that plant chromatin can be arranged into different conformations such as Rabl, Rosette-like, and Bouquet, and that both short- and long-range chromatin interactions occur in Arabidopsis. In this review, we compile the current knowledge about chromosome architecture characteristics in plants, as well as the molecular events and elements (including long non-coding RNAs, histone and DNA modifications, chromatin remodeling complexes, and transcription factors) shaping the genome three-dimensional conformation. Furthermore, we discuss the developmental outputs of genome topology-mediated gene expression regulation. It is becoming increasingly clear that new tools and techniques with higher resolution need to be developed and implemented in Arabidopsis and other model plants in order to better understand chromosome architecture dynamics, from an integrative perspective with other fields of plant biology such as development, stress biology, and finally agriculture.

Development of 5' LTR DNA methylation of latent HIV-1 provirus in cell line models and in long-term-infected individuals

BACKGROUND

Human immunodeficiency virus type 1 (HIV-1) latency represents the major barrier to virus eradication in infected individuals because cells harboring latent HIV-1 provirus are not affected by current antiretroviral therapy (ART). We previously demonstrated that DNA methylation of HIV-1 long terminal repeat (5' LTR) restricts HIV-1 reactivation and, together with chromatin conformation, represents an important mechanism of HIV-1 latency maintenance. Here, we explored the new issue of temporal development of DNA methylation in latent HIV-1 5' LTR.

RESULTS

In the Jurkat CD4(+) T cell model of latency, we showed that the stimulation of host cells contributed to de novo DNA methylation of the latent HIV-1 5' LTR sequences. Consecutive stimulations of model CD4(+) T cell line with TNF-α and PMA or with SAHA contributed to the progressive accumulation of 5' LTR DNA methylation. Further, we showed that once established, the high DNA methylation level of the latent 5' LTR in the cell line model was a stable epigenetic mark. Finally, we explored the development of 5' LTR DNA methylation in the latent reservoir of HIV-1-infected individuals who were treated with ART. We detected low levels of 5' LTR DNA methylation in the resting CD4(+) T cells of the group of patients who were treated for up to 3 years. However, after long-term ART, we observed an accumulation of 5' LTR DNA methylation in the latent reservoir. Importantly, within the latent reservoir of some long-term-treated individuals, we uncovered populations of proviral molecules with a high density of 5' LTR CpG methylation.

CONCLUSIONS

Our data showed the presence of 5' LTR DNA methylation in the long-term reservoir of HIV-1-infected individuals and implied that the transient stimulation of cells harboring latent proviruses may contribute, at least in part, to the methylation of the HIV-1 promoter.

Kisspeptin Induces Dynamic Chromatin Modifications to Control GnRH Gene Expression

In vitro studies have demonstrated an increase in GnRH gene expression associated with an elevated secretory response to kisspeptin administration, suggesting that kisspeptin mediates GnRH expression at both the secretory and pretranslational levels. However, the kisspeptin-mediated intracellular mechanisms associated with the dynamic chromatin modifications modulating GnRH gene expression are unclear. The studies in this manuscript describe specific histone modifications on the enhancer and promoter of the mouse GnRH (mGnRH) gene induced by kisspeptin in GnRH neuronal cell lines (GT1-7 cells). ChIP assays followed by quantitative real-time PCR (qPCR) demonstrate that 15 and 45 min of 10(-9) M kisspeptin significantly increased histone 3 acetylation (H3Ac) at the kisspeptin response element (KsRE) contained between -3446 and -2806 bp of the mGnRH enhancer (GnRHen) in GT1-7 cells, while no changes were observed in the downstream neuron-specific element (NSE). Moreover, kisspeptin specifically induced acetylation of H3AcK14 and K27 and trimethylation of H3 lysine 4 at the KsRE (markers of active chromatin) and no changes in dimethylation of H3K9 (a marker associated with gene repression). Occupancy of RNA Pol II (RNAPII) and a differential carboxyl-terminal domain (CTD) phosphorylation pattern was observed. An interaction between the NSE and the KsRE via a chromatin loop in the mGnRH gene by kisspeptin was detected by the chromosome conformation capture assay (3C). In conclusion, these results demonstrate that kisspeptin induces histone acetylation/methylation and consequently enhances the formation of a chromatin loop in the mGnRH gene which results in known increase in kisspeptin-dependent mGnRH expression.

Structural organization of human replication timing domains

Recent analysis of genome-wide epigenetic modification data, mean replication timing (MRT) profiles and chromosome conformation data in mammals have provided increasing evidence that flexibility in replication origin usage is regulated locally by the epigenetic landscape and over larger genomic distances by the 3D chromatin architecture. Here, we review the recent results establishing some link between replication domains and chromatin structural domains in pluripotent and various differentiated cell types in human. We reconcile the originally proposed dichotomic picture of early and late constant timing regions that replicate by multiple rather synchronous origins in separated nuclear compartments of open and closed chromatins, with the U-shaped MRT domains bordered by "master" replication origins specified by a localized (∼200-300 kb) zone of open and transcriptionally active chromatin from which a replication wave likely initiates and propagates toward the domain center via a cascade of origin firing. We discuss the relationships between these MRT domains, topologically associated domains and lamina-associated domains. This review sheds a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and the determination of differentiation properties.

Stress induces cell dedifferentiation in plants. Biochim Biophys Acta. 2015;1849:378-384. [PMID: 25086338 DOI: 10.1016/j.bbagrm.2014.07.015]

Accumulating evidence lends support to the proposal that a major theme in plant responses to stresses is dedifferentiation, whereby mature cells acquire stem cell features (e.g. open chromatin conformation) prior to acquisition of a new cell fate. In this review, we discuss data addressing plant cell plasticity and provide evidence linking stress, dedifferentiation and a switch in cell fate. We emphasize the epigenetic modifications associated with stress-induced global changes in chromatin structure and conclude with the implications for genetic variation and for induced pluripotent stem cells in animals. It appears that stress is perceived as a signal that directs plant cells to undergo reprogramming (dedifferentiation) as a means for adaptation and in preparation for a stimulus-based acquisition of a new cell fate. This article is part of a Special Issue entitled: Stress as a fundamental theme in cell plasticity.

Reproducibility of 3D chromatin configuration reconstructions

It is widely recognized that the three-dimensional (3D) architecture of eukaryotic chromatin plays an important role in processes such as gene regulation and cancer-driving gene fusions. Observing or inferring this 3D structure at even modest resolutions had been problematic, since genomes are highly condensed and traditional assays are coarse. However, recently devised high-throughput molecular techniques have changed this situation. Notably, the development of a suite of chromatin conformation capture (CCC) assays has enabled elicitation of contacts-spatially close chromosomal loci-which have provided insights into chromatin architecture. Most analysis of CCC data has focused on the contact level, with less effort directed toward obtaining 3D reconstructions and evaluating the accuracy and reproducibility thereof. While questions of accuracy must be addressed experimentally, questions of reproducibility can be addressed statistically-the purpose of this paper. We use a constrained optimization technique to reconstruct chromatin configurations for a number of closely related yeast datasets and assess reproducibility using four metrics that measure the distance between 3D configurations. The first of these, Procrustes fitting, measures configuration closeness after applying reflection, rotation, translation, and scaling-based alignment of the structures. The others base comparisons on the within-configuration inter-point distance matrix. Inferential results for these metrics rely on suitable permutation approaches. Results indicate that distance matrix-based approaches are preferable to Procrustes analysis, not because of the metrics per se but rather on account of the ability to customize permutation schemes to handle within-chromosome contiguity. It has recently been emphasized that the use of constrained optimization approaches to 3D architecture reconstruction are prone to being trapped in local minima. Our methods of reproducibility assessment provide a means for comparing 3D reconstruction solutions so that we can discern between local and global optima by contrasting solutions under perturbed inputs.

Formation of distinct chromatin conformation signatures epigenetically regulate macrophage activation

Microbial-lipopolysacharide (LPS), interleukin 4 (IL-4) and interferon gamma (IFN-γ) polarise macrophages into "innate", "alternative" and "classical", activation states by selective gene regulation. Expression of MARCO, CD200, CD200R1 (innate), MRC1 (alternative) and H2-Eb1 (classical) selectively marks these distinct activation states. Epigenetic events drive such activation upon stimuli and here we study one such mechanism, chromatin conformation signatures implicated in long-range chromatin interactions that regulate transcriptional switch and gene expression. The EpiSwitch™ technology was used to identify and analyse potential markers bordering such conformational signatures for these genes and juxtaposition of markers was compared between resting and activated macrophages. LPS, IL-4 and IFN-γ selectively altered chromatin conformations of their responsive genes in wild type, but not in MyD88(-/-), IL-4R(-/-) and IFN-γR(-/-) macrophages. In addition, two distinct conformations were observed in CD200R1 after LPS and IFN-γ stimulation. In summary, signal-specific alterations in chromatin conformation provide biomarkers that identify and determine distinct gene expression programmes during macrophage activation.

Dynamic changes in chromatin conformation at the TNF transcription start site in T helper lymphocyte subsets

Tumor necrosis factor (TNF) is one of the key primary response genes in the immune system that can be activated by a variety of stimuli. Previous analysis of chromatin accessibility to DNaseI demonstrated open chromatin conformation of the TNF proximal promoter in T cells. Here, using chromatin probing with restriction enzyme EcoNI and micrococcal nuclease we show that in contrast to the proximal promoter, the TNF transcription start site remains in a closed chromatin configuration in primary T helper (Th) cells, but acquires an open state after activation or polarization under Th1 and Th17 conditions. We further demonstrate that transcription factor c-Jun plays a pivotal role in the maintenance of open chromatin conformation at the transcription start site of the TNF gene.