3D Genome Research Articles

Why the estrous cycle matters for neuroscience

BackgroundOvarian hormone fluctuations over the rodent estrous cycle and the human menstrual cycle are known to significantly impact brain physiology and disease risk, yet this variable is largely ignored in preclinical neuroscience research, clinical studies, and psychiatric practice.MethodsTo assess the importance of the estrous cycle information for the analysis of sex differences in neuroscience research, we re-analyzed our previously published data with or without the estrous cycle information, giving a side-by-side comparison of the analyses of behavior, brain structure, gene expression, and 3D genome organization in female and male mice. We also examined and compared the variance of female and male groups across all neurobehavioral measures.ResultsWe show that accounting for the estrous cycle significantly increases the resolution of the neuroscience studies and allows for: (a) identification of masked sex differences; (b) mechanistic insight(s) into the identified sex differences, across different neurobehavioral outcomes, from behavior to molecular phenotypes. We confirm previous findings that female data from either mixed- or staged-female groups are, on average, not more variable than that of males. However, we show that female variability is not, at all, predictive of whether the estrous cycle plays an important role in regulating the outcome of interest.ConclusionsWe argue that “bringing back” the estrous cycle variable to the main stage is important in order to enhance the resolution and quality of the data, to advance the health of women and other menstruators, and to make research more gender-inclusive. We strongly encourage the neuroscience community to incorporate the estrous cycle information in their study design and data analysis, whenever possible, and we debunk some myths that tend to de-emphasize the importance and discourage the inclusion of this critically important biological variable.HighlightsOvarian hormone fluctuation impacts brain physiology and is a major psychiatric risk factor, yet this variable has been overlooked in neuroscience research and psychiatric practice.From rodent behavior to gene regulation, accounting for the estrous cycle increases the resolution of the neuroscience data, allowing identification and mechanistic insight(s) into sex differences.Female variability does not equal (and is not predictive of) the estrous cycle effect and should not be used as a proxy for the effects of ovarian hormones on the outcome of interest.Neuroscience researchers are advised to incorporate the estrous cycle information in their studies to foster more equitable, female- and gender-inclusive research.Studies of the ovarian cycle are especially important for improving women’s mental health.

Subtype-specific 3D genome alteration in acute myeloid leukaemia.

Acute myeloid leukaemia (AML) represents a set of heterogeneous myeloid malignancies, and hallmarks include mutations in epigenetic modifiers, transcription factors and kinases1-5. The extent to which mutations in AML drive alterations in chromatin 3D structure and contribute to myeloid transformation is unclear. Here we use Hi-C and whole-genome sequencing to analyse 25 samples from patients with AML and 7 samples from healthy donors. Recurrent and subtype-specific alterations in A/B compartments, topologically associating domains and chromatin loops were identified. RNA sequencing, ATAC with sequencing and CUT&Tag for CTCF, H3K27ac and H3K27me3 in the same AML samples also revealed extensive and recurrent AML-specific promoter-enhancer and promoter-silencer loops. We validated the role of repressive loops on their target genes by CRISPR deletion and interference. Structural variation-induced enhancer-hijacking and silencer-hijacking events were further identified in AML samples. Hijacked enhancers play a part in AML cell growth, as demonstrated by CRISPR screening, whereas hijacked silencers have a downregulating role, as evidenced by CRISPR-interference-mediated de-repression. Finally, whole-genome bisulfite sequencing of 20 AML and normal samples revealed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Treatment of AML cells with a DNA hypomethylating agent and triple knockdown of DNMT1, DNMT3A and DNMT3B enabled the manipulation of DNA methylation to revert 3D genome organization and gene expression. Overall, this study provides a resource for leukaemia studies and highlights the role of repressive loops and hijacked cis elements in human diseases.

Spatially resolved epigenomic profiling of single cells in complex tissues

The recent development of spatial omics methods has enabled single-cell profiling of the transcriptome and 3D genome organization with high spatial resolution. Expanding the repertoire of spatial omics tools, a spatially resolved single-cell epigenomics method will accelerate understanding of the spatial regulation of cell and tissue functions. Here, we report a method for spatially resolved epigenomic profiling of single cells using in situ tagmentation and transcription followed by multiplexed imaging. We demonstrated the ability to profile histone modifications marking active promoters, putative enhancers, and silent promoters in individual cells, and generated high-resolution spatial atlas of hundreds of active promoters and putative enhancers in embryonic and adult mouse brains. Our results suggested putative promoter-enhancer pairs and enhancer hubs regulating developmentally important genes. We envision this approach will be generally applicable to spatial profiling of epigenetic modifications and DNA-binding proteins, advancing our understanding of how gene expression is spatiotemporally regulated by the epigenome.

3D genome organization links non-coding disease-associated variants to genes.

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.

Folding Features and Dynamics of 3D Genome Architecture in Plant Fungal Pathogens.

The folding and dynamics of three-dimensional (3D) genome organization are fundamental for eukaryotes executing genome functions but have been largely unexplored in nonmodel fungi. Using high-throughput sequencing coupled with chromosome conformation capture (Hi-C) data, we generated two chromosome-level assemblies for Puccinia striiformis f. sp. tritici, a fungus causing stripe rust disease on wheat, for studying 3D genome architectures of plant pathogenic fungi. The chromatin organization of the fungus followed a combination of the fractal globule model and the equilibrium globule model. Surprisingly, chromosome compartmentalization was not detected. Dynamics of 3D genome organization during two developmental stages of P. striiformis f. sp. tritici indicated that regulation of gene activities might be independent of the changes of genome organization. In addition, chromatin conformation conservation was found to be independent of genome sequence synteny conservation among different fungi. These results highlighted the distinct folding principles of fungal 3D genomes. Our findings should be an important step toward a holistic understanding of the principles and functions of genome architecture across different eukaryotic kingdoms. IMPORTANCE Previously, our understanding of 3D genome architecture has mainly come from model mammals, insects, and plants. However, the organization and regulatory functions of 3D genomes in fungi are largely unknown. In this study, we comprehensively investigated P. striiformis f. sp. tritici, a plant fungal pathogen, and revealed distinct features of the 3D genome, comparing it with the universal folding feature of 3D genomes in higher eukaryotic organisms. We further suggested that there might be distinct regulatory mechanisms of gene expression that are independent of chromatin organization changes during the developmental stages of this rust fungus. Moreover, we showed that the evolutionary pattern of 3D genomes in this fungus is also different from the cases in mammalian genomes. In addition, the genome assembly pipeline and the generated two chromosome-level genomes will be valuable resources. These results highlighted the unexplored distinct features of 3D genome organization in fungi. Therefore, our study provided complementary knowledge to holistically understand the organization and functions of 3D genomes across different eukaryotes.

DLoopCaller: A deep learning approach for predicting genome-wide chromatin loops by integrating accessible chromatin landscapes.

In recent years, major advances have been made in various chromosome conformation capture technologies to further satisfy the needs of researchers for high-quality, high-resolution contact interactions. Discriminating the loops from genome-wide contact interactions is crucial for dissecting three-dimensional(3D) genome structure and function. Here, we present a deep learning method to predict genome-wide chromatin loops, called DLoopCaller, by combining accessible chromatin landscapes and raw Hi-C contact maps. Some available orthogonal data ChIA-PET/HiChIP and Capture Hi-C were used to generate positive samples with a wider contact matrix which provides the possibility to find more potential genome-wide chromatin loops. The experimental results demonstrate that DLoopCaller effectively improves the accuracy of predicting genome-wide chromatin loops compared to the state-of-the-art method Peakachu. Moreover, compared to two of most popular loop callers, such as HiCCUPS and Fit-Hi-C, DLoopCaller identifies some unique interactions. We conclude that a combination of chromatin landscapes on the one-dimensional genome contributes to understanding the 3D genome organization, and the identified chromatin loops reveal cell-type specificity and transcription factor motif co-enrichment across different cell lines and species.

Genome organization in cardiomyocytes expressing mutated A-type lamins.

Cardiomyopathy is a myocardial disorder, in which the heart muscle is structurally and functionally abnormal, often leading to heart failure. Dilated cardiomyopathy is characterized by a compromised left ventricular function and contributes significantly to the heart failure epidemic, which represents a staggering clinical and public health problem worldwide. Gene mutations have been identified in 35% of patients with dilated cardiomyopathy. Pathogenic variants in LMNA, encoding nuclear A-type lamins, are one of the major causative causes of dilated cardiomyopathy (i.e. CardioLaminopathy). A-type lamins are type V intermediate filament proteins, which are the main components of the nuclear lamina. The nuclear lamina is connected to the cytoskeleton on one side, and to the chromatin on the other side. Among the models proposed to explain how CardioLaminopathy arises, the “chromatin model” posits an effect of mutated A-type lamins on the 3D genome organization and thus on the transcription activity of tissue-specific genes. Chromatin contacts with the nuclear lamina via specific genomic regions called lamina-associated domains lamina-associated domains. These LADs play a role in the chromatin organization and gene expression regulation. This review focuses on the identification of LADs and chromatin remodeling in cardiac muscle cells expressing mutated A-type lamins and discusses the methods and relevance of these findings in disease.

Decoding the spatial chromatin organization and dynamic epigenetic landscapes of macrophage cells during differentiation and immune activation

Immunocytes dynamically reprogram their gene expression profiles during differentiation and immunoresponse. However, the underlying mechanism remains elusive. Here, we develop a single-cell Hi-C method and systematically delineate the 3D genome and dynamic epigenetic atlas of macrophages during these processes. We propose “degree of disorder” to measure genome organizational patterns inside topologically-associated domains, which is correlated with the chromatin epigenetic states, gene expression, and chromatin structure variability in individual cells. Furthermore, we identify that NF-κB initiates systematic chromatin conformation reorganization upon Mycobacterium tuberculosis infection. The integrated Hi-C, eQTL, and GWAS analysis depicts the atlas of the long-range target genes of mycobacterial disease susceptible loci. Among these, the SNP rs1873613 is located in the anchor of a dynamic chromatin loop with LRRK2, whose inhibitor AdoCbl could be an anti-tuberculosis drug candidate. Our study provides comprehensive resources for the 3D genome structure of immunocytes and sheds insights into the order of genome organization and the coordinated gene transcription during immunoresponse.

Ultrafast and interpretable single-cell 3D genome analysis with Fast-Higashi.

Single-cell Hi-C (scHi-C) technologies can probe three-dimensional (3D) genome structures in individual cells. However, existing scHi-C analysis methods are hindered by the data quality and complex 3D genome patterns. The lack of computational scalability and interpretability poses further challenges for large-scale analysis. Here, we introduce Fast-Higashi, an ultrafast and interpretable method based on tensor decomposition and partial random walk with restart, enabling joint identification of cell identities and chromatin meta-interactions from sparse scHi-C data. Extensive evaluations demonstrate the advantage of Fast-Higashi over existing methods, leading to improved delineation of rare cell types and continuous developmental trajectories. Fast-Higashi can directly identify 3D genome features that define distinct cell types and help elucidate cell-type-specific connections between genome structure and function. Moreover, Fast-Higashi can generalize to incorporate other single-cell omics data. Fast-Higashi provides a highly efficient and interpretable scHi-C analysis solution that is applicable to a broad range of biological contexts.

Right on target: Chromatin jets arise from targeted cohesin loading in wild-type cells

Uncovering an informative feature of 3D genome structure, Guo etal. (2022) describe chromatin jets in quiescent murine thymocytes: 1-2 Mb structures formed by targeted cohesin loading at narrow accessible chromatin regions and visible as prominent off-diagonal stripes on contact maps.

Establishment of 3D chromatin structure after fertilization and the metabolic switch at the morula-to-blastocyst transition require CTCF.

The eukaryotic genome is organized in 3D at different scales. This structure is driven and maintained by different chromatin states and by architectural factors, such as the zinc finger protein CTCF. Zygotic genome structure is established de novo after fertilization, but its impact during the first stages of mammalian development is unclear. We show that deletion of Ctcf in mouse embryos impairs the establishment of chromatin structure, but the first cell fate decision is unperturbed and embryos are viable until the late blastocyst. Furthermore, maternal CTCF is not necessary for development. Gene expression changes in metabolic and protein homeostasis programs that occur during the morula-to-blastocyst transition depend on CTCF. However, these changes do not correlate with disruption of chromatin but with binding of CTCF to the promoter of downregulated genes. Our results show that CTCF regulates both 3D genome organization and transcription during mouse preimplantation development, but as independent processes.

Bioinformatics Pipelines for Identification of Super-Enhancers and 3D Chromatin Contacts.

Precise regulation of gene expression is integral in development. Emerging studies have highlighted that super-enhancers (SEs), which are clusters of multiple enhancers, play critical roles in regulating cell type-specific gene expression via 3D chromatin, thereby defining the cellular identities of given cells. Here we provide optimized bioinformatics pipelines to identify SEs and 3D chromatin contacts. Our pipelines encompass the processing of chromatin immunoprecipitation sequencing (ChIP-seq) data to identify SEs and the processing of genome-wide chromosome conformation capture (Hi-C) data. We can then infer long-range chromatin contacts between SEs and other genomic regions. This integrative computational approach, which can be applied to CUT&RUN and CUT&Tag, alternative technologies to ChIP-seq, allows us to identify genomic locations of SEs and their 3D genome configuration, whereby multiple SEs act in concert. We show an analysis of mouse spermatogenesis as an example of this application.

Dynamics of Viral and Host 3D Genome Structure upon Infection.

Eukaryotic chromatin is highly organized in the 3D nuclear space and dynamically regulated in response to environmental stimuli. This genomic organization is arranged in a hierarchical fashion to support various cellular functions, including transcriptional regulation of gene expression. Like other host cellular mechanisms, viral pathogens utilize and modulate host chromatin architecture and its regulatory machinery to control features of their life cycle, such as lytic versus latent status. Combined with previous research focusing on individual loci, recent global genomic studies employing conformational assays coupled with high-throughput sequencing technology have informed models for host and, in some cases, viral 3D chromosomal structure re-organization during infection and the contribution of these alterations to virus-mediated diseases. Here, we review recent discoveries and progress in host and viral chromatin structural dynamics during infection, focusing on a subset of DNA (human herpesviruses and HPV) as well as RNA (HIV, influenza virus and SARS-CoV-2) viruses. An understanding of how host and viral genomic structure affect gene expression in both contexts and ultimately viral pathogenesis can facilitate the development of novel therapeutic strategies.

The three-dimensional landscape of cortical chromatin accessibility in Alzheimer's disease.

To characterize the dysregulation of chromatin accessibility in Alzheimer's disease (AD), we generated 636 ATAC-seq libraries from neuronal and nonneuronal nuclei isolated from the superior temporal gyrus and entorhinal cortex of 153 AD cases and 56 controls. By analyzing a total of ~20 billion read pairs, we expanded the repertoire of known open chromatin regions (OCRs) in the human brain and identified cell-type-specific enhancer-promoter interactions. We show that interindividual variability in OCRs can be leveraged to identify cis-regulatory domains (CRDs) that capture the three-dimensional structure of the genome (3D genome). We identified AD-associated effects on chromatin accessibility, the 3D genome and transcription factor (TF) regulatory networks. For one of the most AD-perturbed TFs, USF2, we validated its regulatory effect on lysosomal genes. Overall, we applied a systematic approach to understanding the role of the 3D genome in AD. We provide all data as an online resource for widespread community-based analysis.

CohesinDB: a comprehensive database for decoding cohesin-related epigenomes, 3D genomes and transcriptomes in human cells.

Cohesin is a multifunctional protein responsible for transcriptional regulation and chromatin organization. Cohesin binds to chromatin at tens of thousands of distinct sites in a conserved or tissue-specific manner, whereas the function of cohesin varies greatly depending on the epigenetic properties of specific chromatin loci. Cohesin also extensively mediates cis-regulatory modules (CRMs) and chromatin loops. Even though next-generation sequencing technologies have provided a wealth of information on different aspects of cohesin, the integration and exploration of the resultant massive cohesin datasets are not straightforward. Here, we present CohesinDB (https://cohesindb.iqb.u-tokyo.ac.jp), a comprehensive multiomics cohesin database in human cells. CohesinDB includes 2043 epigenomics, transcriptomics and 3D genomics datasets from 530 studies involving 176 cell types. By integrating these large-scale data, CohesinDB summarizes three types of 'cohesin objects': 751 590 cohesin binding sites, 957 868 cohesin-related chromatin loops and 2229 500 cohesin-related CRMs. Each cohesin object is annotated with locus, cell type, classification, function, 3D genomics and cis-regulatory information. CohesinDB features a user-friendly interface for browsing, searching, analyzing, visualizing and downloading the desired information. CohesinDB contributes a valuable resource for all researchers studying cohesin, epigenomics, transcriptional regulation and chromatin organization.

ParticleChromo3D: a Particle Swarm Optimization algorithm for chromosome 3D structure prediction from Hi-C data

BackgroundThe three-dimensional (3D) structure of chromatin has a massive effect on its function. Because of this, it is desirable to have an understanding of the 3D structural organization of chromatin. To gain greater insight into the spatial organization of chromosomes and genomes and the functions they perform, chromosome conformation capture (3C) techniques, particularly Hi-C, have been developed. The Hi-C technology is widely used and well-known because of its ability to profile interactions for all read pairs in an entire genome. The advent of Hi-C has greatly expanded our understanding of the 3D genome, genome folding, gene regulation and has enabled the development of many 3D chromosome structure reconstruction methods.ResultsHere, we propose a novel approach for 3D chromosome and genome structure reconstruction from Hi-C data using Particle Swarm Optimization (PSO) approach called ParticleChromo3D. This algorithm begins with a grouping of candidate solution locations for each chromosome bin, according to the particle swarm algorithm, and then iterates its position towards a global best candidate solution. While moving towards the optimal global solution, each candidate solution or particle uses its own local best information and a randomizer to choose its path. Using several metrics to validate our results, we show that ParticleChromo3D produces a robust and rigorous representation of the 3D structure for input Hi-C data. We evaluated our algorithm on simulated and real Hi-C data in this work. Our results show that ParticleChromo3D is more accurate than most of the existing algorithms for 3D structure reconstruction.ConclusionsOur results also show that constructed ParticleChromo3D structures are very consistent, hence indicating that it will always arrive at the global solution at every iteration. The source code for ParticleChromo3D, the simulated and real Hi-C datasets, and the models generated for these datasets are available here: https://github.com/OluwadareLab/ParticleChromo3D

Three-dimensional genome organization in immune cell fate and function.

Immune cell development and activation demand the precise and coordinated control of transcriptional programmes. Three-dimensional (3D) organization of the genome has emerged as an important regulator of chromatin state, transcriptional activity and cell identity by facilitating or impeding long-range genomic interactions among regulatory elements and genes. Chromatin folding thus enables cell type-specific and stimulus-specific transcriptional responses to extracellular signals, which are essential for the control of immune cell fate, for inflammatory responses and for generating a diverse repertoire of antigen receptor specificities. Here, we review recent findings connecting 3D genome organization to the control of immune cell differentiation and function, and discuss how alterations in genome folding may lead to immune dysfunction and malignancy.

Integrative epigenomics in chronic lymphocytic leukaemia: Biological insights and clinical applications.

Chronic lymphocytic leukaemia (CLL) is not only characterised by driver genetic alterations but by extensive epigenetic changes. Over the last decade, epigenomic studies have described the DNA methylome, chromatin accessibility, histone modifications and the three-dimensional (3D) genome architecture of CLL. Beyond its regulatory role, the DNA methylome contains imprints of the cellular origin and proliferative history of CLL cells. These two aspects are strong independent prognostic factors. Integrative analyses of chromatin marks have uncovered novel regulatory elements and altered transcription factor networks as non-genetic means mediating gene deregulation in CLL. Additionally, CLL cells display a disease-specific pattern of 3D genome interactions. From the technological perspective, we are currently witnessing a transition from bulk omics to single-cell analyses. This review aims at summarising the major findings from the epigenomics field as well as providing a prospect of the present and future of single-cell analyses in CLL.

Prostate Cancer Transcriptomic Regulation by the Interplay of Germline Risk Alleles, Somatic Mutations, and 3D Genomic Architecture.

Many prostate cancer germline risk alleles are enriched in the noncoding regions of the genome and are hypothesized to regulate transcription. We present a 3D genomics framework to unravel risk SNP function and describe the widespread germline-somatic interplay in transcription control. This article is highlighted in the In This Issue feature, p. 2711.

Abstract PR008: The 3D epigenome of rhabdomyosarcoma

Abstract Rhabdomyosarcoma (RMS) is driven by transcription factors that set up super-enhancers: large clusters of transcriptional machinery accumulated at oncogenes. Super-enhancers are also the key vulnerability of RMS, as we have revealed by CRISPR screening and chemical genomics. However, the identification of super-enhancers is inconsistent, because they are extracted from 1-dimensional ChIP-seq data, but are in fact 3D genomic objects (“super-clusters”) that form biomolecular condensates (phase-separated droplets). We have now applied new 3D genomics approaches to enable proper discovery of these epigenomic features, also yielding proper enhancer-gene connectivity maps. Furthermore, super-clusters are co-built by transcription factors recruiting p300 (acetylation writer) and BRD4 (acetylation reader). We have generated preliminary evidence that directly perturbing/degrading p300 and BRD4 are the most effective and selective means to shut down super-clusters at all oncogenes driving this childhood cancer. We applied Absolute Quantification of Architecture (AQuA) HiChIP, and were able to uncover the drug-induced 3D folding abnormalities that selectively halt RNA-Polymerase 2 at super-cluster driven oncogenes in this sarcoma. Citation Format: Berkley E. Gryder, Issra Osman, Adam Durbin, Jun Qi. The 3D epigenome of rhabdomyosarcoma [abstract]. In: Proceedings of the AACR Special Conference: Sarcomas; 2022 May 9-12; Montreal, QC, Canada. Philadelphia (PA): AACR; Clin Cancer Res 2022;28(18_Suppl):Abstract nr PR008.