3D Genome Research Articles

Spatial Chromosome Organization and Adaptation of Escherichia coli under Heat Stress.

The spatial organization of bacterial chromosomes is crucial for cellular functions. It remains unclear how bacterial chromosomes adapt to high-temperature stress. This study delves into the 3D genome architecture and transcriptomic responses of Escherichia coli under heat-stress conditions to unravel the intricate interplay between the chromosome structure and environmental cues. By examining the role of macrodomains, chromosome interaction domains (CIDs), and nucleoid-associated proteins (NAPs), this work unveils the dynamic changes in chromosome conformation and gene expression patterns induced by high-temperature stress. It was observed that, under heat stress, the short-range interaction frequency of the chromosomes decreased, while the long-range interaction frequency of the Ter macrodomain increased. Furthermore, two metrics, namely, Global Compactness (GC) and Local Compactness (LC), were devised to measure and compare the compactness of the chromosomes based on their 3D structure models. The findings in this work shed light on the molecular mechanisms underlying thermal adaptation and chromosomal organization in bacterial cells, offering valuable insights into the complex inter-relationships between environmental stimuli and genomic responses.

Reprogramming of 3D genome structure underlying HSPC development in zebrafish

BackgroundDevelopment of hematopoietic stem and progenitor cells (HSPC) is a multi-staged complex process that conserved between zebrafish and mammals. Understanding the mechanism underlying HSPC development is a holy grail of hematopoietic biology, which is helpful for HSPC clinical application. Chromatin conformation plays important roles in transcriptional regulation and cell fate decision; however, its dynamic and role in HSPC development is poorly investigated.MethodsWe performed chromatin structure and multi-omics dissection across different stages of HSPC developmental trajectory in zebrafish for the first time, including Hi-C, RNA-seq, ATAC-seq, H3K4me3 and H3K27ac ChIP-seq.ResultsThe chromatin organization of zebrafish HSPC resemble mammalian cells with similar hierarchical structure. We revealed the multi-scale reorganization of chromatin structure and its influence on transcriptional regulation and transition of cell fate during HSPC development. Nascent HSPC is featured by loose conformation with obscure structure at all layers. Notably, PU.1 was identified as a potential factor mediating formation of promoter-involved loops and regulating gene expression of HSPC.ConclusionsOur results provided a global view of chromatin structure dynamics associated with development of zebrafish HSPC and discovered key transcription factors involved in HSPC chromatin interactions, which will provide new insights into the epigenetic regulatory mechanisms underlying vertebrate HSPC fate decision.

FISHnet: Detecting chromatin domains in single-cell sequential Oligopaints imaging data.

Sequential Oligopaints DNA FISH is an imaging technique that measures higher-order genome folding at single-allele resolution via multiplexed, probe-based tracing. Currently there is a paucity of algorithms to identify 3D genome features in sequential Oligopaints data. Here, we present FISHnet, a graph theory method based on optimization of network modularity to detect chromatin domains and boundaries in pairwise distance matrices. FISHnet uncovers cell type-specific domain-like folding patterns on single alleles, thus enabling future studies aiming to elucidate the role for single-cell folding variation on genome function.

From enhancers to genome conformation: complex transcriptional control underlies expression of a single herpesviral gene.

Complex transcriptional control is a conserved feature of both eukaryotes and the viruses that infect them. Here, we illustrate this by combining high-density functional genomics, expression profiling, and viral-specific chromosome conformation capture to define with unprecedented detail the transcriptional regulation of a single gene, ORF68, from Kaposi's sarcoma-associated herpesvirus (KSHV). We first identified seven cis-regulatory regions by densely tiling the ~154 kb KSHV genome with CRISPRi. A parallel Cas9 nuclease screen indicated that three of these regions act as promoters of genes that regulate ORF68. RNA expression profiling demonstrated that three more of these regions act by either repressing or enhancing other distal viral genes involved in ORF68 transcriptional regulation. Finally, we tracked how the 3D structure of the viral genome changes during its lifecycle, revealing that these enhancing regulatory elements are physically closer to their targets when active, and that disrupting some elements caused large-scale changes to the 3D genome. These data enable us to construct a complete model revealing that the mechanistic diversity of this essential regulatory circuit matches that of human genes.

A three-dimensional genomics view for speciation research

Genomic information is folded in a three-dimensional (3D) structure, a rarely explored evolutionary driver of speciation. Technological advances now enable the study of 3D genome structures (3DGSs) across the Tree of Life. At the onset of 3D speciation genomics, we discuss the putative roles of 3DGSs in speciation.

Estrogen-induced chromatin looping changes identify a subset of functional regulatory elements.

Transcriptional enhancers can regulate individual or multiple genes through long-range three-dimensional (3D) genome interactions, and these interactions are commonly altered in cancer. Yet, the functional relationship between changes in 3D interactions associated with regulatory regions and differential gene expression appears context-dependent. In this study, we used HiChiP to capture changes in 3D genome interactions between active regulatory regions of endometrial cancer cells in response to estrogen treatment and uncovered significant differential long-range interactions that are strongly enriched for estrogen receptor α (ER) bound sites (ERBS). The ERBS anchoring differential loops with either a gene's promoter or distal regions were correlated with larger transcriptional responses to estrogen compared to ERBS not involved in differential interactions. To functionally test this observation, CRISPR-based Enhancer-i was used to deactivate specific ERBS, which revealed a wide range of effects on the transcriptional response to estrogen. However, these effects are only subtly and not significantly stronger for ERBS in differential loops. In addition, we observed an enrichment of 3D interactions between the promoters of estrogen up-regulated genes and found that looped promoters can work together cooperatively. Overall, our work suggests that changes in 3D genome structure upon estrogen treatment identify some functionally important regulatory regions; however, these changes aren't required for a transcriptional response to E2 in endometrial cancer cells.

Altered chromatin topologies caused by balanced chromosomal translocation lead to central iris hypoplasia

Despite the advent of genomic sequencing, molecular diagnosis remains unsolved in approximately half of patients with Mendelian disorders, largely due to unclarified functions of noncoding regions and the difficulty in identifying complex structural variations. In this study, we map a unique form of central iris hypoplasia in a large family to 6q15-q23.3 and 18p11.31-q12.1 using a genome-wide linkage scan. Long-read sequencing reveals a balanced translocation t(6;18)(q22.31;p11.22) with intergenic breakpoints. By performing Hi-C on induced pluripotent stem cells from a patient, we identify two chromatin topologically associating domains spanning across the breakpoints. These alterations lead the ectopic chromatin interactions between APCDD1 on chromosome 18 and enhancers on chromosome 6, resulting in upregulation of APCDD1. Notably, APCDD1 is specifically localized in the iris of human eyes. Our findings demonstrate that noncoding structural variations can lead to Mendelian diseases by disrupting the 3D genome structure and resulting in altered gene expression.

A Constitutive Heterochromatic Region Shapes Genome Organization and Impacts Gene Expression in Neurospora crassa.

Genome organization is essential for proper function, including gene expression. In metazoan genome organization, chromatin loops and Topologically Associated Domains (TADs) facilitate local gene clustering, while chromosomes form distinct nuclear territories characterized by compartmentalization of silent heterochromatin at the nuclear periphery and active euchromatin in the nucleus center. A similar hierarchical organization occurs in the fungus Neurospora crassa where its seven chromosomes form a Rabl conformation, where heterochromatic centromeres and telomeres independently cluster at the nuclear membrane, while interspersed heterochromatic loci in Neurospora aggregate across megabases of linear genomic distance for forming TAD-like structures. However, the role of individual heterochromatic loci in normal genome organization and function is unknown. Here, we examined the genome organization of a Neurospora strain harboring a ~47.4 kilobase facultative (temporarily silent) heterochromatic region deletion, as well as the genome organization of a strain deleted of a 110.6 kilobase permanently silent constitutive heterochromatic region. While the facultative heterochromatin deletion had little effect on local chromatin structure, the constitutive heterochromatin deletion altered local TAD-like structures, gene expression, and the predicted 3D genome structure by qualitatively repositioning genes into the nucleus center. Our work elucidates the role of individual heterochromatic regions for genome organization and function.

Comprehensive multimodal and multiomic profiling reveals epigenetic and transcriptional reprogramming in lung tumors.

Epigenomic mechanisms are critically involved in mediation of genetic and environmental factors that underlie cancer development. Histone modifications represent highly informative epigenomic marks that reveal activation and repression of gene activities and dysregulation of transcriptional control due to tumorigenesis. Here, we present a comprehensive epigenomic and transcriptomic mapping of 18 tumor and 20 non-neoplastic tissues from non-small cell lung adenocarcinoma patients. Our profiling covers 5 histone marks including activating (H3K4me3, H3K4me1, and H3K27ac) and repressive (H3K27me3 and H3K9me3) marks and the transcriptome using only 20 mg of tissue per sample, enabled by low-input omic technologies. Using advanced integrative bioinformatic analysis, we uncovered cancer-driving signaling cascade networks, changes in 3D genome modularity, and differential expression and functionalities of transcription factors and noncoding RNAs. Many of these identified genes and regulatory molecules showed no significant change in their expression or a single epigenomic modality, emphasizing the power of integrative multimodal and multiomic analysis using patient samples.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer.

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

Accelerating 3D genomics data analysis with Microcket

The three-dimensional (3D) organization of genome is fundamental to cell biology. To explore 3D genome, emerging high-throughput approaches have produced billions of sequencing reads, which is challenging and time-consuming to analyze. Here we present Microcket, a package for mapping and extracting interacting pairs from 3D genomics data, including Hi-C, Micro-C, and derivant protocols. Microcket utilizes a unique read-stitch strategy that takes advantage of the long read cycles in modern DNA sequencers; benchmark evaluations reveal that Microcket runs much faster than the current tools along with improved mapping efficiency, and thus shows high potential in accelerating and enhancing the biological investigations into 3D genome. Microcket is freely available at https://github.com/hellosunking/Microcket.

Pairtools: From sequencing data to chromosome contacts.

The field of 3D genome organization produces large amounts of sequencing data from Hi-C and a rapidly-expanding set of other chromosome conformation protocols (3C+). Massive and heterogeneous 3C+ data require high-performance and flexible processing of sequenced reads into contact pairs. To meet these challenges, we present pairtools-a flexible suite of tools for contact extraction from sequencing data. Pairtools provides modular command-line interface (CLI) tools that can be flexibly chained into data processing pipelines. The core operations provided by pairtools are parsing of.sam alignments into Hi-C pairs, sorting and removal of PCR duplicates. In addition, pairtools provides auxiliary tools for building feature-rich 3C+ pipelines, including contact pair manipulation, filtration, and quality control. Benchmarking pairtools against popular 3C+ data pipelines shows advantages of pairtools for high-performance and flexible 3C+ analysis. Finally, pairtools provides protocol-specific tools for restriction-based protocols, haplotype-resolved contacts, and single-cell Hi-C. The combination of CLI tools and tight integration with Python data analysis libraries makes pairtools a versatile foundation for a broad range of 3C+ pipelines.

SuPreMo: a computational tool for streamlining in silico perturbation using sequence-based predictive models.

The increasing development of sequence-based machine learning models has raised the demand for manipulating sequences for this application. However, existing approaches to edit and evaluate genome sequences using models have limitations, such as incompatibility with structural variants, challenges in identifying responsible sequence perturbations, and the need for vcf file inputs and phased data. To address these bottlenecks, we present Sequence Mutator for Predictive Models (SuPreMo), a scalable and comprehensive tool for performing and supporting in silico mutagenesis experiments. We then demonstrate how pairs of reference and perturbed sequences can be used with machine learning models to prioritize pathogenic variants or discover new functional sequences. SuPreMo was written in Python, and can be run using only one line of code to generate both sequences and 3D genome disruption scores. The codebase, instructions for installation and use, and tutorials are on the GitHub page: https://github.com/ketringjoni/SuPreMo.

Heterochromatic 3D genome organization is directed by HP1a- and H3K9-dependent and independent mechanisms

Whether and how histone post-translational modifications and the proteins that bind them drive 3D genome organization remains unanswered. Here, we evaluate the contribution of H3K9-methylated constitutive heterochromatin to 3D genome organization in Drosophila tissues. We find that the predominant organizational feature of wild-type tissues is the segregation of euchromatic chromosome arms from heterochromatic pericentromeres. Reciprocal perturbation of HP1a⋅H3K9me binding, using a point mutation in the HP1a chromodomain or replacement of the replication-dependent histone H3 with H3K9R mutant histones, revealed that HP1a binding to methylated H3K9 in constitutive heterochromatin is required to limit contact frequency between pericentromeres and chromosome arms and regulate the distance between arm and pericentromeric regions. Surprisingly, the self-association of pericentromeric regions is largely preserved despite the loss of H3K9 methylation and HP1a occupancy. Thus, the HP1a⋅H3K9 interaction contributes to but does not solely drive the segregation of euchromatin and heterochromatin inside the nucleus.

Stratifying TAD boundaries pinpoints focal genomic regions of regulation, damage, and repair.

Advances in chromatin mapping have exposed the complex chromatin hierarchical organization in mammals, including topologically associating domains (TADs) and their substructures, yet the functional implications of this hierarchy in gene regulation and disease progression are not fully elucidated. Our study delves into the phenomenon of shared TAD boundaries, which are pivotal in maintaining the hierarchical chromatin structure and regulating gene activity. By integrating high-resolution Hi-C data, chromatin accessibility, and DNA double-strand breaks (DSBs) data from various cell lines, we systematically explore the complex regulatory landscape at high-level TAD boundaries. Our findings indicate that these boundaries are not only key architectural elements but also vibrant hubs, enriched with functionally crucial genes and complex transcription factor binding site-clustered regions. Moreover, they exhibit a pronounced enrichment of DSBs, suggesting a nuanced interplay between transcriptional regulation and genomic stability. Our research provides novel insights into the intricate relationship between the 3D genome structure, gene regulation, and DNA repair mechanisms, highlighting the role of shared TAD boundaries in maintaining genomic integrity and resilience against perturbations. The implications of our findings extend to understanding the complexities of genomic diseases and open new avenues for therapeutic interventions targeting the structural and functional integrity of TAD boundaries.

The role of physicochemical processes in the formation of the 3D genome and compartmentalization of the cell nucleus.

The review analyzes the role of physicochemical processes in the formation of the function-dependent architecture of the cell nucleus, built on the platform of a folded genome. The main attention is paid to various forms of the phase separation process, primarily the processes of liquid-liquid phase separation and polymer-polymer phase separation. The role of these processes in the formation of chromatin compartments and maintenance of three-dimensional genome architecture is discussed in detail. The relationship between genome activity and the creation of functional compartments in the cell nucleus is also analyzed.

The GWAS SNP rs80207740 modulates erythrocyte traits via allele-specific binding of IKZF1 and targeting XPO7 gene.

Genome-wide association studies have identified many single nucleotide polymorphisms (SNPs) associated with erythrocyte traits. However, the functional variants and their working mechanisms remain largely unknown. Here, we reported that the SNP of rs80207740, which was associated with red blood cell (RBC) volume and hemoglobin content across populations, conferred enhancer activity to XPO7 gene via allele-differentially binding to Ikaros family zinc finger 1 (IKZF1). We showed that the region around rs80207740 was an erythroid-specific enhancer using reporter assays, and that the G-allele further enhanced activity. 3D genome evidence showed that the enhancer interacted with the XPO7 promoter, and eQTL analysis suggested that the G-allele upregulated expression of XPO7. We further showed that the rs80207740-G allele facilitated the binding of transcription factor IKZF1 in EMSA and ChIP analyses. Knockdown of IKZF1 and GATA1 resulted in decreased expression of Xpo7 in both human and mouse erythroid cells. Finally, we constructed Xpo7 knockout mouse by CRISPR/Cas9 and observed anemic phenotype with reduced volume and hemoglobin content of RBC, consistent to the effect of rs80207740 on erythrocyte traits. Overall, our study demonstrated that rs80207740 modulated erythroid indices by regulating IKZF1 binding and Xpo7 expression.

DNA methylation-based high-resolution mapping of long-distance chromosomal interactions in nucleosome-depleted regions.

3C-based methods have significantly advanced our understanding of 3D genome organization. However, it remains a formidable task to precisely capture long-range chromosomal interactions between individual loci, such as those between promoters and distal enhancers. Here, we present Methyltransferase Targeting-based chromosome Architecture Capture (MTAC), a method that maps the contacts between a target site (viewpoint) and the rest of the genome in budding yeast with high resolution and sensitivity. MTAC detects hundreds of intra- and inter-chromosomal interactions within nucleosome-depleted regions (NDRs) that cannot be captured by 4C, Hi-C, or Micro-C. By applying MTAC to various viewpoints, we find that (1) most long-distance chromosomal interactions detected by MTAC reflect tethering by the nuclear pore complexes (NPCs), (2) genes co-regulated by methionine assemble into inter-chromosomal clusters near NPCs upon activation, (3) mediated by condensin, the mating locus forms a highly specific interaction with the recombination enhancer (RE) in a mating-type specific manner, and (4) correlation of MTAC signals among NDRs reveal spatial mixing and segregation of the genome. Overall, these results demonstrate MTAC as a powerful tool to resolve fine-scale long-distance chromosomal interactions and provide insights into the 3D genome organization.

Differential 3D genome architecture and imprinted gene expression: cause or consequence?

Imprinted genes provide an attractive paradigm to unravel links between transcription and genome architecture. The parental allele-specific expression of these essential genes - which are clustered in chromosomal domains - is mediated by parental methylation imprints at key regulatory DNA sequences. Recent chromatin conformation capture (3C)-based studies show differential organization of topologically associating domains between the parental chromosomes at imprinted domains, in embryonic stem and differentiated cells. At several imprinted domains, differentially methylated regions show allelic binding of the insulator protein CTCF, and linked focal retention of cohesin, at the non-methylated allele only. This generates differential patterns of chromatin looping between the parental chromosomes, already in the early embryo, and thereby facilitates the allelic gene expression. Recent research evokes also the opposite scenario, in which allelic transcription contributes to the differential genome organization, similarly as reported for imprinted X chromosome inactivation. This may occur through epigenetic effects on CTCF binding, through structural effects of RNA Polymerase II, or through imprinted long non-coding RNAs that have chromatin repressive functions. The emerging picture is that epigenetically-controlled differential genome architecture precedes and facilitates imprinted gene expression during development, and that at some domains, conversely, the mono-allelic gene expression also influences genome architecture.

Multiscale 3D genome organization underlies duck fatty liver with no adipose inflammation or serious injury

Non-alcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease. Little is known about how gene expression and chromatin structure are regulated in NAFLD due to lack of suitable model. Ducks naturally develop fatty liver similar to serious human non-alcoholic fatty liver (NAFL) without adipose inflammation and liver fibrosis, thus serves as a good model for investigating molecular mechanisms of adipose metabolism and anti-inflammation. Here, we constructed a NAFLD model without adipose inflammation and liver fibrosis in ducks. By performing dynamic pathological and transcriptomic analyses, we identified critical genes involving in regulation of the NF-κB and MHCII signaling, which usually lead to adipose inflammation and liver fibrosis. We further generated dynamic three-dimensional chromatin maps during liver fatty formation and recovery. This showed that ducks enlarged hepatocyte cell nuclei to reduce inter-chromosomal interaction, decompress chromatin structure, and alter strength of intra-TAD and loop interactions during fatty liver formation. These changes partially contributed to the tight control the NF-κB and the MHCII signaling. Our analysis uncovers duck chromatin reorganization might be advantageous to maintain liver regenerative capacity and reduce adipose inflammation. These findings shed light on new strategies for NAFLD control.