Breaking TADs: How Alterations of Chromatin Domains Result in Disease

Three-dimensional chromatin reorganization regulates B cell development during ageing

The contribution of three-dimensional genome organization to physiological ageing is not well known. Here we show that large-scale chromatin reorganization distinguishes young and old bone marrow progenitor (pro-) B cells. These changes result in increased interactions at the compartment level and reduced interactions within topologically associated domains (TADs). The gene encoding Ebf1, a key B cell regulator, switches from compartment A to B with age. Genetically reducing Ebf1 recapitulates some features of old pro-B cells. TADs that are most reduced with age contain genes important for B cell development, including the immunoglobulin heavy chain (Igh) locus. Weaker intra-TAD interactions at Igh correlate with altered variable (V), diversity (D) and joining (J) gene recombination. Our observations implicate three-dimensional chromatin reorganization as a major driver of pro-B cell phenotypes that impair B lymphopoiesis with age.

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.

A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

A potassium-chloride co-transporter with altered genome architecture functions as a suppressor in glioma

Gliomas, the most lethal tumours in brain, have a poor prognosis despite accepting standard treatment. Limited benefits from current therapies can be attributed to genetic, epigenetic and microenvironmental cues that affect cell programming and drive tumour heterogeneity. Through the analysis of Hi-C data, we identified a potassium-chloride co-transporter SLC12A5 associated with disrupted topologically associating domain which was downregulated in tumour tissues. Multiple independent glioma cohorts were included to analyse the characterization of SLC12A5 and found it was significantly associated with pathological features, prognostic value, genomic alterations, transcriptional landscape and drug response. We constructed two SLC12A5 overexpression cell lines to verify the function of SLC12A5 that suppressed tumour cell proliferation and migration in vitro. In addition, SLC12A5 was also positively associated with GABAA receptor activity and negatively associated with pro-tumour immune signatures and immunotherapy response. Collectively, our study provides a comprehensive characterization of SLC12A5 in glioma and supports SLC12A5 as a potential suppressor of disease progression.

Rare germline structural variants increase risk for pediatric solid tumors

Pediatric solid tumors are rare malignancies that represent a leading cause of death by disease among children in developed countries. The early age-of-onset of these tumors suggests that germline genetic factors are involved, yet conventional germline testing for short coding variants in established predisposition genes only identifies pathogenic events in 10-15% of patients. Here, we examined the role of germline structural variants (SVs)-an underexplored form of germline variation-in pediatric extracranial solid tumors using germline genome sequencing of 1,766 affected children, their 943 unaffected relatives, and 6,665 adult controls. We discovered a sex-biased association between very large (>1 megabase) germline chromosomal abnormalities and a four-fold increased risk of solid tumors in male children. The overall impact of germline SVs was greatest in neuroblastoma, where we revealed burdens of ultra-rare SVs that cause loss-of-function of highly expressed, mutationally intolerant, neurodevelopmental genes, as well as noncoding SVs predicted to disrupt three-dimensional chromatin domains in neural crest-derived tissues. Collectively, our results implicate rare germline SVs as a predisposing factor to pediatric solid tumors that may guide future studies and clinical practice.

Enhancer-promoter specificity in gene transcription: molecular mechanisms and disease associations

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

Hi-C, a chromatin 3D structure technique advancing the functional genomics of immune cells

The functional performance of immune cells relies on a complex transcriptional regulatory network. The three-dimensional structure of chromatin can affect chromatin status and gene expression patterns, and plays an important regulatory role in gene transcription. Currently available techniques for studying chromatin spatial structure include chromatin conformation capture techniques and their derivatives, chromatin accessibility sequencing techniques, and others. Additionally, the recently emerged deep learning technology can be utilized as a tool to enhance the analysis of data. In this review, we elucidate the definition and significance of the three-dimensional chromatin structure, summarize the technologies available for studying it, and describe the research progress on the chromatin spatial structure of dendritic cells, macrophages, T cells, B cells, and neutrophils.

Identifying TAD-like domains on single-cell Hi-C data by graph embedding and changepoint detection

MOTIVATION

Topologically associating domains (TADs) are fundamental building blocks of 3D genome. TAD-like domains in single cells are regarded as the underlying genesis of TADs discovered in bulk cells. Understanding the organization of TAD-like domains helps to get deeper insights into their regulatory functions. Unfortunately, it remains a challenge to identify TAD-like domains on single-cell Hi-C data due to its ultra-sparsity.

RESULTS

We propose scKTLD, an in silico tool for the identification of TAD-like domains on single-cell Hi-C data. It takes Hi-C contact matrix as the adjacency matrix for a graph, embeds the graph structures into a low-dimensional space with the help of sparse matrix factorization followed by spectral propagation, and the TAD-like domains can be identified using a kernel-based changepoint detection in the embedding space. The results tell that our scKTLD is superior to the other methods on the sparse contact matrices, including downsampled bulk Hi-C data as well as simulated and experimental single-cell Hi-C data. Besides, we demonstrated the conservation of TAD-like domain boundaries at single-cell level apart from heterogeneity within and across cell types, and found that the boundaries with higher frequency across single cells are more enriched for architectural proteins and chromatin marks, and they preferentially occur at TAD boundaries in bulk cells, especially at those with higher hierarchical levels.

AVAILABILITY AND IMPLEMENTATION

scKTLD is freely available at https://github.com/lhqxinghun/scKTLD.

Computational methods for analysing multiscale 3D genome organization

Recent progress in whole-genome mapping and imaging technologies has enabled the characterization of the spatial organization and folding of the genome in the nucleus. In parallel, advanced computational methods have been developed to leverage these mapping data to reveal multiscale three-dimensional (3D) genome features and to provide a more complete view of genome structure and its connections to genome functions such as transcription. Here, we discuss how recently developed computational tools, including machine-learning-based methods and integrative structure-modelling frameworks, have led to a systematic, multiscale delineation of the connections among different scales of 3D genome organization, genomic and epigenomic features, functional nuclear components and genome function. However, approaches that more comprehensively integrate a wide variety of genomic and imaging datasets are still needed to uncover the functional role of 3D genome structure in defining cellular phenotypes in health and disease.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

A chromatin-regulated biphasic circuit coordinates IL-1β-mediated inflammation

Inflammation is characterized by a biphasic cycle consisting initially of a proinflammatory phase that is subsequently resolved by anti-inflammatory processes. Interleukin-1β (IL-1β) is a master regulator of proinflammation and is encoded within the same topologically associating domain (TAD) as IL-37, which is an anti-inflammatory cytokine that opposes the function of IL-1β. Within this TAD, we identified a long noncoding RNA called AMANZI, which negatively regulates IL-1β expression and trained immunity through the induction of IL37 transcription. We found that the activation of IL37 occurs through the formation of a dynamic long-range chromatin contact that leads to the temporal delay of anti-inflammatory responses. The common variant rs16944 present in AMANZI augments this regulatory circuit, predisposing individuals to enhanced proinflammation or immunosuppression. Our work illuminates a chromatin-mediated biphasic circuit coordinating expression of IL-1β and IL-37, thereby regulating two functionally opposed states of inflammation from within a single TAD.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Genomic characterization of IDH-mutant astrocytoma progression to grade 4 in the treatment setting

As the progression of low-grade diffuse astrocytomas into grade 4 tumors significantly impacts patient prognosis, a better understanding of this process is of paramount importance for improved patient care. In this project, we analyzed matched IDH-mutant astrocytomas before and after progression to grade 4 from six patients (discovery cohort) with genome-wide sequencing, 21 additional patients with targeted sequencing, and 33 patients from Glioma Longitudinal AnalySiS cohort for validation. The Cancer Genome Atlas data from 595 diffuse gliomas provided supportive information. All patients in our discovery cohort received radiation, all but one underwent chemotherapy, and no patient received temozolomide (TMZ) before progression to grade 4 disease. One case in the discovery cohort exhibited a hypermutation signature associated with the inactivation of the MSH2 and DNMT3A genes. In other patients, the number of chromosomal rearrangements and deletions increased in grade 4 tumors. The cell cycle checkpoint gene CDKN2A, or less frequently RB1, was most commonly inactivated after receiving both chemo- and radiotherapy when compared to other treatment groups. Concomitant activating PDGFRA/MET alterations were detected in tumors that acquired a homozygous CDKN2A deletion. NRG3 gene was significantly downregulated and recurrently altered in progressed tumors. Its decreased expression was associated with poorer overall survival in both univariate and multivariate analysis. We also detected progression-related alterations in RAD51B and other DNA repair pathway genes associated with the promotion of error-prone DNA repair, potentially facilitating tumor progression. In our retrospective analysis of patient treatment and survival timelines (n = 75), the combination of postoperative radiation and chemotherapy (mainly TMZ) outperformed radiation, especially in the grade 3 tumor cohort, in which it was typically given after primary surgery. Our results provide further insight into the contribution of treatment and genetic alterations in cell cycle, growth factor signaling, and DNA repair-related genes to tumor evolution and progression.

SHIELD: a platform for high-throughput screening of barrier-type DNA elements in human cells

Chromatin boundary elements contribute to the partitioning of mammalian genomes into topological domains to regulate gene expression. Certain boundary elements are adopted as DNA insulators for safe and stable transgene expression in mammalian cells. These elements, however, are ill-defined and less characterized in the non-coding genome, partially due to the lack of a platform to readily evaluate boundary-associated activities of putative DNA sequences. Here we report SHIELD (Site-specific Heterochromatin Insertion of Elements at Lamina-associated Domains), a platform tailored for the high-throughput screening of barrier-type DNA elements in human cells. SHIELD takes advantage of the high specificity of serine integrase at heterochromatin, and exploits the natural heterochromatin spreading inside lamina-associated domains (LADs) for the discovery of potent barrier elements. We adopt SHIELD to evaluate the barrier activity of 1000 DNA elements in a high-throughput manner and identify 8 candidates with barrier activities comparable to the core region of cHS4 element in human HCT116 cells. We anticipate SHIELD could facilitate the discovery of novel barrier DNA elements from the non-coding genome in human cells.

Fundamental insights into the correlation between chromosome configuration and transcription

Eukaryotic chromosomes exhibit a hierarchical organization that spans a spectrum of length scales, ranging from sub-regions known as loops, which typically comprise hundreds of base pairs, to much larger chromosome territories that can encompass a few mega base pairs. Chromosome conformation capture experiments that involve high-throughput sequencing methods combined with microscopy techniques have enabled a new understanding of inter- and intra-chromosomal interactions with unprecedented details. This information also provides mechanistic insights on the relationship between genome architecture and gene expression. In this article, we review the recent findings on three-dimensional interactions among chromosomes at the compartment, topologically associating domain, and loop levels and the impact of these interactions on the transcription process. We also discuss current understanding of various biophysical processes involved in multi-layer structural organization of chromosomes. Then, we discuss the relationships between gene expression and genome structure from perturbative genome-wide association studies. Furthermore, for a better understanding of how chromosome architecture and function are linked, we emphasize the role of epigenetic modifications in the regulation of gene expression. Such an understanding of the relationship between genome architecture and gene expression can provide a new perspective on the range of potential future discoveries and therapeutic research.

Diagnostic Gene Panel Testing in (Non)-Syndromic Patients with Cleft Lip, Alveolus and/or Palate in the Netherlands

Objectives

Clefts of the lip, alveolus and/or palate (CLA/P) are the most common craniofacial congenital malformations in humans. These oral clefts can be divided into non-syndromic (isolated) and syndromic forms. Many cleft-related syndromes are clinically variable and genetically heterogeneous, making it challenging to distinguish syndromic from non-syndromic cases. Recognition of syndromic/genetic causes is important for personalized tailored care, identification of (unrecognized) comorbidities, and accurate genetic counseling. Therefore, next generation sequencing (NGS)-based targeted gene panel testing is increasingly implemented in diagnostics of CLA/P patients. In this retrospective study, we assess the yield of NGS gene panel testing in a cohort of CLA/P cases.

Methods

Whole exome sequencing (WES) followed by variant detection and interpretation in an a priori selected set of genes associated with CLA/P phenotypes was performed in 212 unrelated CLA/P patients after genetic counseling between 2015 and 2020. Medical records including family history and results of additional genetic tests were evaluated.

Results

In 24 CLA/P cases (11.3%), a pathogenic genetic variant was identified. Twenty out of these 24 had a genetic syndrome requiring specific monitoring and follow-up. Six of these 24 cases (25%) were presumed to be isolated CLA/P cases prior to testing, corresponding to 2.8% of the total cohort. In eight CLA/P cases (3.8%) without a diagnosis after NGS-based gene panel testing, a molecular diagnosis was established by additional genetic analyses (e.g., SNP array, single gene testing, trio WES).

Conclusion

This study illustrates NGS-based gene panel testing is a powerful diagnostic tool in the diagnostic workup of CLA/P patients. Also, in apparently isolated cases and non-familial cases, a genetic diagnosis can be identified. Early diagnosis facilitates personalized care for patients and accurate genetic counseling of their families.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Recent advances in chromosome capture techniques unraveling 3D genome architecture in germ cells, health, and disease

In eukaryotes, the genome does not emerge in a specific shape but rather as a hierarchial bundle within the nucleus. This multifaceted genome organization consists of multiresolution cellular structures, such as chromosome territories, compartments, and topologically associating domains, which are frequently defined by architecture, design proteins including CTCF and cohesin, and chromatin loops. This review briefly discusses the advances in understanding the basic rules of control, chromatin folding, and functional areas in early embryogenesis. With the use of chromosome capture techniques, the latest advancements in technologies for visualizing chromatin interactions come close to revealing 3D genome formation frameworks with incredible detail throughout all genomic levels, including at single-cell resolution. The possibility of detecting variations in chromatin architecture might open up new opportunities for disease diagnosis and prevention, infertility treatments, therapeutic approaches, desired exploration, and many other application scenarios.

Enhancer Function in the 3D Genome

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

Three-dimensional genome rewiring in loci with human accelerated regions

Human accelerated regions (HARs) are conserved genomic loci that evolved at an accelerated rate in the human lineage and may underlie human-specific traits. We generated HARs and chimpanzee accelerated regions with an automated pipeline and an alignment of 241 mammalian genomes. Combining deep learning with chromatin capture experiments in human and chimpanzee neural progenitor cells, we discovered a significant enrichment of HARs in topologically associating domains containing human-specific genomic variants that change three-dimensional (3D) genome organization. Differential gene expression between humans and chimpanzees at these loci suggests rewiring of regulatory interactions between HARs and neurodevelopmental genes. Thus, comparative genomics together with models of 3D genome folding revealed enhancer hijacking as an explanation for the rapid evolution of HARs.

Topologically associating domain boundaries are required for normal genome function

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes^1-3, but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11-80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening.

3D genome alterations and editing in pathology

The human genome is folded into a multi-level 3D structure that controls many nuclear functions including gene expression. Recently, alterations in 3D genome organization were associated with several genetic diseases and cancer. As a consequence, experimental approaches are now being developed to modify the global 3D genome organization and that of specific loci. Here, we discuss emerging experimental approaches of 3D genome editing that may prove useful in biomedicine.

The non-coding genome in Autism Spectrum Disorders

Autism Spectrum Disorders (ASD) are a group of neurodevelopmental disorders (NDDs) characterized by difficulties in social interaction and communication, repetitive behavior, and restricted interests. While ASD have been proven to have a strong genetic component, current research largely focuses on coding regions of the genome. However, non-coding DNA, which makes up for ∼99% of the human genome, has recently been recognized as an important contributor to the high heritability of ASD, and novel sequencing technologies have been a milestone in opening up new directions for the study of the gene regulatory networks embedded within the non-coding regions. Here, we summarize current progress on the contribution of non-coding alterations to the pathogenesis of ASD and provide an overview of existing methods allowing for the study of their functional relevance, discussing potential ways of unraveling ASD's "missing heritability".

Rewiring of the 3D genome during acquisition of carboplatin resistance in a triple-negative breast cancer patient-derived xenograft

Changes in the three-dimensional (3D) structure of the genome are an emerging hallmark of cancer. Cancer-associated copy number variants and single nucleotide polymorphisms promote rewiring of chromatin loops, disruption of topologically associating domains (TADs), active/inactive chromatin state switching, leading to oncogene expression and silencing of tumor suppressors. However, little is known about 3D changes during cancer progression to a chemotherapy-resistant state. We integrated chromatin conformation capture (Hi-C), RNA-seq, and whole-genome sequencing obtained from triple-negative breast cancer patient-derived xenograft primary tumors (UCD52) and carboplatin-resistant samples and found increased short-range (< 2 Mb) interactions, chromatin looping, formation of TAD, chromatin state switching into a more active state, and amplification of ATP-binding cassette transporters. Transcriptome changes suggested the role of long-noncoding RNAs in carboplatin resistance. Rewiring of the 3D genome was associated with TP53, TP63, BATF, FOS-JUN family of transcription factors and led to activation of aggressiveness-, metastasis- and other cancer-related pathways. Integrative analysis highlighted increased ribosome biogenesis and oxidative phosphorylation, suggesting the role of mitochondrial energy metabolism. Our results suggest that 3D genome remodeling may be a key mechanism underlying carboplatin resistance.

Oncogenic epidermal growth factor receptor signal-induced histone deacetylation suppresses chemokine gene expression in human lung adenocarcinoma

Epidermal growth factor receptor (EGFR)-mutated (mt) lung adenocarcinoma (LA) is refractory to immune checkpoint inhibitors (ICIs). However, the mechanisms have not been fully elucidated. CD8⁺ T cell infiltration was significantly lower in EGFR-mt than in EGFR-wild-type LA, which was associated with suppression of chemokine expression. Since this T cell-deserted tumor microenvironment may lead to the refractoriness of ICIs against EGFR-mt LA, we investigated the mechanism by focusing on the regulation of chemokine expression. The expression of C-X-C motif ligand (CXCL) 9, 10 and 11, which constitute a gene cluster on chromosome 4, was suppressed under EGFR signaling. The assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) revealed open chromatin peaks near this gene cluster following EGFR-tyrosine kinase inhibitor (TKI) treatment. The histone deacetylase (HDAC) inhibitor recovered the expression of CXCL9, 10 and 11 in EGFR-mt LA. Nuclear HDAC activity, as well as histone H3 deacetylation, were dependent on oncogenic EGFR signaling. Furthermore, the Cleavage Under Targets and Tagmentation (CUT & Tag) assay revealed a histone H3K27 acetylation peak at 15 kb upstream of CXCL11 after treatment with EGFR-TKI, which corresponded to one of the open chromatin peaks detected by ATAC-seq. The data suggest that EGFR-HDAC axis mediates silencing of the chemokine gene cluster through chromatin conformational change, which might be relevant to the ICI resistance by creating T cell-deserted tumor microenvironment. Targeting this axis may develop a new therapeutic strategy to overcome the ICI resistance of EGFR-mt LA.

Breakpoint analysis for cytogenetically balanced translocation revealed unexpected complex structural abnormalities and suggested the position effect for MEF2C

Many disease-causing genes have been identified by determining the breakpoints of balanced chromosomal translocations. Recent progress in genomic analysis has accelerated the analysis of chromosomal translocation-breakpoints at the nucleotide level. Using a long-read whole-genome sequence, we analyzed the breakpoints of the cytogenetically balanced chromosomal translocation t(5;15)(q21;26.3), which was confirmed to be of de novo origin, in a patient with a neurodevelopmental disorder. The results showed complex rearrangements with seven fragments consisting of five breakpoint-junctions (BJs). Four of the five BJs showed microhomologies of 1-3-bp, and only one BJ displayed a signature of blunt-end ligation, indicating chromothripsis as the underlying mechanism. Although the BJs did not disrupt any disease-causing gene, the clinical features of the patient were compatible with MEF2C haploinsufficiency syndrome. Complex rearrangements were located approximately 2.5-Mb downstream of MEF2C. Therefore, position effects were considered the mechanism of the occurrence of MEF2C haploinsufficiency syndrome.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

A de novo transcription-dependent TAD boundary underpins critical multiway interactions during antibody class switch recombination

Interactions between transcription and cohesin-mediated loop extrusion can influence 3D chromatin architecture. However, their relevance in biology is unclear. Here, we report a direct role for such interactions in the mechanism of antibody class switch recombination (CSR) at the murine immunoglobulin heavy chain locus (Igh). Using Tri-C to measure higher-order multiway interactions on single alleles, we find that the juxtaposition (synapsis) of transcriptionally active donor and acceptor Igh switch (S) sequences, an essential step in CSR, occurs via the interaction of loop extrusion complexes with a de novo topologically associating domain (TAD) boundary formed via transcriptional activity across S regions. Surprisingly, synapsis occurs predominantly in proximity to the 3' CTCF-binding element (3'CBE) rather than the Igh super-enhancer, suggesting a two-step mechanism whereby transcription of S regions is not topologically coupled to synapsis, as has been previously proposed. Altogether, these insights advance our understanding of how 3D chromatin architecture regulates CSR.

Chromosome conformation capture approaches to investigate 3D genome architecture in Ankylosing Spondylitis

Ankylosing Spondylitis (AS) is a chronic inflammatory arthritis of the spine exhibiting a strong genetic background. The mechanistic and functional understanding of the AS-associated genomic loci, identified with Genome Wide Association Studies (GWAS), remains challenging. Chromosome conformation capture (3C) and derivatives are recent techniques which are of great help in elucidating the spatial genome organization and of enormous support in uncover a mechanistic explanation for disease-associated genetic variants. The perturbation of three-dimensional (3D) genome hierarchy may lead to a plethora of human diseases, including rheumatological disorders. Here we illustrate the latest approaches and related findings on the field of genome organization, highlighting how the instability of 3D genome conformation may be among the causes of rheumatological disease phenotypes. We suggest a new perspective on the inclusive potential of a 3C approach to inform GWAS results in rheumatic diseases. 3D genome organization may ultimately lead to a more precise and comprehensive functional interpretation of AS association, which is the starting point for emerging and more specific therapies.

Interphase chromosomes of the Aedes aegypti mosquito are liquid crystalline and can sense mechanical cues

We use data-driven physical simulations to study the three-dimensional architecture of the Aedes aegypti genome. Hi-C maps exhibit both a broad diagonal and compartmentalization with telomeres and centromeres clustering together. Physical modeling reveals that these observations correspond to an ensemble of 3D chromosomal structures that are folded over and partially condensed. Clustering of the centromeres and telomeres near the nuclear lamina appears to be a necessary condition for the formation of the observed structures. Further analysis of the mechanical properties of the genome reveals that the chromosomes of Aedes aegypti, by virtue of their atypical structural organization, are highly sensitive to the deformation of the nuclei. This last finding provides a possible physical mechanism linking mechanical cues to gene regulation.

Enhancer architecture and chromatin accessibility constrain phenotypic space during Drosophila development

Developmental enhancers bind transcription factors and dictate patterns of gene expression during development. Their molecular evolution can underlie phenotypical evolution, but the contributions of the evolutionary pathways involved remain little understood. Here, using mutation libraries in Drosophila melanogaster embryos, we observed that most point mutations in developmental enhancers led to changes in gene expression levels but rarely resulted in novel expression outside of the native pattern. In contrast, random sequences, often acting as developmental enhancers, drove expression across a range of cell types; random sequences including motifs for transcription factors with pioneer activity acted as enhancers even more frequently. Our findings suggest that the phenotypic landscapes of developmental enhancers are constrained by enhancer architecture and chromatin accessibility. We propose that the evolution of existing enhancers is limited in its capacity to generate novel phenotypes, whereas the activity of de novo elements is a primary source of phenotypic novelty.

A Comparison of Topologically Associating Domain Callers Based on Hi-C Data

Topologically associating domains (TADs) are local chromatin interaction domains, which have been shown to play an important role in gene expression regulation. TADs were originally discovered in the investigation of 3D genome organization based on High-throughput Chromosome Conformation Capture (Hi-C) data. Continuous considerable efforts have been dedicated to developing methods for detecting TADs from Hi-C data. Different computational methods for TADs identification vary in their assumptions and criteria in calling TADs. As a consequence, the TADs called by these methods differ in their similarities and biological features they are enriched in. In this work, we performed a systematic comparison of twenty-six TAD callers. We first compared the TADs and gaps between adjacent TADs across different methods, resolutions, and sequencing depths. We then assessed the quality of TADs and TAD boundaries according to three criteria: the decay of contact frequencies over the genomic distance, enrichment and depletion of regulatory elements around TAD boundaries, and reproducibility of TADs and TAD boundaries in replicate samples. Last, due to the lack of a gold standard of TADs, we also evaluated the performance of the methods on synthetic datasets. We discussed the key principles of TAD callers, and pinpointed current situation in the detection of TADs. We provide a concise, comprehensive, and systematic framework for evaluating the performance of TAD callers, and expect our work will provide useful guidance in choosing suitable approaches for the detection and evaluation of TADs.

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.

Clinical and genomic delineation of the new proximal 19p13.3 microduplication syndrome

A small but growing body of scientific literature is emerging about clinical findings in patients with 19p13.3 microdeletion or duplication. Recently, a proximal 19p13.3 microduplication syndrome was described, associated with growth delay, microcephaly, psychomotor delay and dysmorphic features. The aim of our study was to better characterize the syndrome associated with duplications in the proximal 19p13.3 region (prox 19p13.3 dup), and to propose a comprehensive analysis of the underlying genomic mechanism. We report the largest cohort of patients with prox 19p13.3 dup through a collaborative study. We collected 24 new patients with terminal or interstitial 19p13.3 duplication characterized by array-based Comparative Genomic Hybridization (aCGH). We performed mapping, phenotype-genotype correlations analysis, critical region delineation and explored three-dimensional chromatin interactions by analyzing Topologically Associating Domains (TADs). We define a new 377 kb critical region (CR 1) in chr19: 3,116,922-3,494,377, GRCh37, different from the previously described critical region (CR 2). The new 377 kb CR 1 includes a TAD boundary and two enhancers whose common target is PIAS4. We hypothesize that duplications of CR 1 are responsible for tridimensional structural abnormalities by TAD disruption and misregulation of genes essentials for the control of head circumference during development, by breaking down the interactions between enhancers and the corresponding targeted gene.

The Role of Long Noncoding RNA (lncRNAs) Biomarkers in Renal Cell Carcinoma

Renal cell carcinoma is one of the common cancers whose incidence and mortality are continuously growing worldwide. Initially, this type of tumour is usually asymptomatic. Due to the lack of reliable diagnostic markers, one-third of ccRCC patients already have distant metastases at the time of diagnosis. This underlines the importance of establishing biomarkers that would enable the prediction of the disease's course and the risk of metastasis. LncRNA, which modulates genes at the epigenetic, transcriptional, and post-transcriptional levels, appears promising. The actions of lncRNA involve sponging and sequestering target miRNAs, thus affecting numerous biological processes. Studies have confirmed the involvement of RNAs in various diseases, including RCC. In this review, we focused on MALAT1 (a marker of serious pathological changes and a factor in the promotion of tumorigenesis), RCAT1 (tumour promoter in RCC), DUXAP9 (a plausible marker of localized ccRCC), TCL6 (exerting tumour-suppressive effects in renal cancer), LINC00342 (acting as an oncogene), AGAP2 Antisense1 (plausible predictor of RCC progression), DLEU2 (factor promoting tumours growth via the regulation of epithelial-mesenchymal transition), NNT-AS1 (sponge of miR-22 contributing to tumour progression), LINC00460 (favouring ccRCC development and progression) and Lnc-LSG1 (a factor that may stimulate ccRCC metastasis).

A cohesin traffic pattern genetically linked to gene regulation

Cohesin-mediated loop extrusion has been shown to be blocked at specific cis-elements, including CTCF sites, producing patterns of loops and domain boundaries along chromosomes. Here we explore such cis-elements, and their role in gene regulation. We find that transcription termination sites of active genes form cohesin- and RNA polymerase II-dependent domain boundaries that do not accumulate cohesin. At these sites, cohesin is first stalled and then rapidly unloaded. Start sites of transcriptionally active genes form cohesin-bound boundaries, as shown before, but are cohesin-independent. Together with cohesin loading, possibly at enhancers, these sites create a pattern of cohesin traffic that guides enhancer-promoter interactions. Disrupting this traffic pattern, by removing CTCF, renders cells sensitive to knockout of genes involved in transcription initiation, such as the SAGA complexes, and RNA processing such DEAD/H-Box RNA helicases. Without CTCF, these factors are less efficiently recruited to active promoters.

MoDLE: high-performance stochastic modeling of DNA loop extrusion interactions

DNA loop extrusion emerges as a key process establishing genome structure and function. We introduce MoDLE, a computational tool for fast, stochastic modeling of molecular contacts from DNA loop extrusion capable of simulating realistic contact patterns genome wide in a few minutes. MoDLE accurately simulates contact maps in concordance with existing molecular dynamics approaches and with Micro-C data and does so orders of magnitude faster than existing approaches. MoDLE runs efficiently on machines ranging from laptops to high performance computing clusters and opens up for exploratory and predictive modeling of 3D genome structure in a wide range of settings.

Predictable and stable epimutations induced during clonal plant propagation with embryonic transcription factor

Clonal propagation is frequently used in commercial plant breeding and biotechnology programs because it minimizes genetic variation, yet it is not uncommon to observe clonal plants with stable phenotypic changes, a phenomenon known as somaclonal variation. Several studies have linked epigenetic modifications induced during regeneration with this newly acquired phenotypic variation. However, the factors that determine the extent of somaclonal variation and the molecular changes underpinning this process remain poorly understood. To address this gap in our knowledge, we compared clonally propagated Arabidopsis thaliana plants derived from somatic embryogenesis using two different embryonic transcription factors- RWP-RK DOMAIN-CONTAINING 4 (RKD4) or LEAFY COTYLEDON2 (LEC2) and from two epigenetically distinct founder tissues. We found that both the epi(genetic) status of the explant and the regeneration protocol employed play critical roles in shaping the molecular and phenotypic landscape of clonal plants. Phenotypic variation in regenerated plants can be largely explained by the inheritance of tissue-specific DNA methylation imprints, which are associated with specific transcriptional and metabolic changes in sexual progeny of clonal plants. For instance, regenerants were particularly affected by the inheritance of root-specific epigenetic imprints, which were associated with an increased accumulation of salicylic acid in leaves and accelerated plant senescence. Collectively, our data reveal specific pathways underpinning the phenotypic and molecular variation that arise and accumulate in clonal plant populations.

Lamin A/C-dependent chromatin architecture safeguards naïve pluripotency to prevent aberrant cardiovascular cell fate and function

Tight control of cell fate choices is crucial for normal development. Here we show that lamin A/C plays a key role in chromatin organization in embryonic stem cells (ESCs), which safeguards naïve pluripotency and ensures proper cell fate choices during cardiogenesis. We report changes in chromatin compaction and localization of cardiac genes in Lmna-/- ESCs resulting in precocious activation of a transcriptional program promoting cardiomyocyte versus endothelial cell fate. This is accompanied by premature cardiomyocyte differentiation, cell cycle withdrawal and abnormal contractility. Gata4 is activated by lamin A/C loss and Gata4 silencing or haploinsufficiency rescues the aberrant cardiovascular cell fate choices induced by lamin A/C deficiency. We uncover divergent functions of lamin A/C in naïve pluripotent stem cells and cardiomyocytes, which have distinct contributions to the transcriptional alterations of patients with LMNA-associated cardiomyopathy. We conclude that disruption of lamin A/C-dependent chromatin architecture in ESCs is a primary event in LMNA loss-of-function cardiomyopathy.

CASPIAN: A method to identify chromatin topological associated domains based on spatial density cluster

With the development of Hi-C technology, the detection of topologically associated domains (TADs) boundaries plays an important role in exploring the relationship between gene structure and expression. However, a method that can identify accurate TAD boundaries from the Hi-C contact matrix with different resolutions is currently lacking. We proposed a method named CASPIAN that can identify chromatin TAD boundaries based on the spatial density clustering algorithm. CASPIAN requires few parameters to call TADs. This method is realized using the hierarchical density-based clustering method HDBSCAN, where the distance of pairwise bins is calculated based on three distance metrics (Euclidean, Manhattan, and Chebyshev distance metric) to adapt to the characteristics of the Hi-C contact matrix generated from simulation experiments or normalized methods. Our results show that, same as standard methods (e.g., Insulation Score, TopDom), CASPIAN can enrich factors related to promoting the gene expression, such as CTCF, H3K4me1, H3K4me3, RAD21, POLR2A, and SMC3. We also calculated the approximate proportion of various factors anchored at the TAD boundaries to observe the distribution of these factors surrounding the TAD boundaries. In conclusion, CASPIAN is an easy method to explore the relationship between transcription factors and TAD boundaries. CASPIAN is available online (https://gitee.com/ghaiyan/caspian).

Cohesin is required for long-range enhancer action at the Shh locus

The regulatory landscapes of developmental genes in mammals can be complex, with enhancers spread over many hundreds of kilobases. It has been suggested that three-dimensional genome organization, particularly topologically associating domains formed by cohesin-mediated loop extrusion, is important for enhancers to act over such large genomic distances. By coupling acute protein degradation with synthetic activation by targeted transcription factor recruitment, here we show that cohesin, but not CTCF, is required for activation of the target gene Shh by distant enhancers in mouse embryonic stem cells. Cohesin is not required for activation directly at the promoter or by an enhancer located closer to the Shh gene. Our findings support the hypothesis that chromatin compaction via cohesin-mediated loop extrusion allows for genes to be activated by enhancers that are located many hundreds of kilobases away in the linear genome and suggests that cohesin is dispensable for enhancers located more proximally.

Manipulating complex chromatin folding via CRISPR-guided bioorthogonal chemistry

Precise manipulation of chromatin folding is important for understanding the relationship between the three-dimensional genome and nuclear function. Existing tools can reversibly establish individual chromatin loops but fail to manipulate two or more chromatin loops. Here, we engineer a powerful CRISPR system which can manipulate multiple chromatin contacts using bioorthogonal reactions, termed the bioorthogonal reaction-mediated programmable chromatin loop (BPCL) system. The multiinput BPCL system employs engineered single-guide RNAs recognized by discrete bioorthogonal adaptors to independently and dynamically control different chromatin loops formation without cross-talk in the same cell or to establish hubs of multiway chromatin contacts. We use the BPCL system to successfully juxtapose the pluripotency gene promoters to enhancers and activate their endogenous expression. BPCL enables us to independently engineer multiway chromatin contacts without cross-talk, which provides a way to precisely dissect the high complexity and dynamic nature of chromatin folding.

Understanding the function of regulatory DNA interactions in the interpretation of non-coding GWAS variants

Genome-wide association studies (GWAS) have identified a vast number of variants associated with various complex human diseases and traits. However, most of these GWAS variants reside in non-coding regions producing no proteins, making the interpretation of these variants a daunting challenge. Prior evidence indicates that a subset of non-coding variants detected within or near cis-regulatory elements (e.g., promoters, enhancers, silencers, and insulators) might play a key role in disease etiology by regulating gene expression. Advanced sequencing- and imaging-based technologies, together with powerful computational methods, enabling comprehensive characterization of regulatory DNA interactions, have substantially improved our understanding of the three-dimensional (3D) genome architecture. Recent literature witnesses plenty of examples where using chromosome conformation capture (3C)-based technologies successfully links non-coding variants to their target genes and prioritizes relevant tissues or cell types. These examples illustrate the critical capability of 3D genome organization in annotating non-coding GWAS variants. This review discusses how 3D genome organization information contributes to elucidating the potential roles of non-coding GWAS variants in disease etiology.

Interstitial deletion 4p15.32p16.1 and complex chromoplexy in a female proband with severe neurodevelopmental delay, growth failure and dysmorphism

Complex chromosomal rearrangements involve the restructuring of genetic material within a single chromosome or across multiple chromosomes. These events can cause serious human disease by disrupting coding DNA and gene regulatory elements via deletions, duplications, and structural rearrangements. Here we describe a 5-year-old female with severe developmental delay, dysmorphic features, multi-suture craniosynostosis, and growth failure found to have a complex series of balanced intra- and inter-chromosomal rearrangements involving chromosomes 4, 11, 13, and X. Initial clinical studies were performed by karyotype, chromosomal microarray, and FISH with research-based short-read genome sequencing coupled with sanger sequencing to precisely map her breakpoints to the base pair resolution to understand the molecular basis of her phenotype. Genome analysis revealed two pathogenic deletions at 4p16.1-p15.32 and 4q31.1, accounting for her developmental delay and dysmorphism. We identified over 60 breakpoints, many with blunt ends and limited homology, supporting a role for non-homologous end joining in restructuring and resolution of the seminal chromoplexy event. We propose that the complexity of our patient’s genomic rearrangements with a high number of breakpoints causes dysregulation of gene expression by three-dimensional chromatin interactions or topologically associating domains leading to growth failure and craniosynostosis. Our work supports an important role for genome sequencing in understanding the molecular basis of complex chromosomal rearrangements in human disease.

The relevance of chromatin architecture to genome rearrangements in Drosophila

DNA within chromosomes in the nucleus is non-randomly organized into chromosome territories, compartments and topologically associated domains (TADs). Chromosomal rearrangements have the potential to alter chromatin organization and modify gene expression leading to selection against these structural variants. Drosophila pseudoobscura has a wealth of naturally occurring gene arrangements that were generated by overlapping inversion mutations caused by two chromosomal breaks that rejoin the central region in reverse order. Unlike humans, Drosophila inversion heterozygotes do not have negative effects associated with crossing over during meiosis because males use achiasmate mechanisms for proper segregation, and aberrant recombinant meiotic products generated in females are lost in polar bodies. As a result, Drosophila populations are found to harbour extensive inversion polymorphisms. It is not clear, however, whether chromatin architecture constrains which inversions breakpoints persist in populations. We mapped the breakpoints of seven inversions in D. pseudoobscura to the TAD map to determine if persisting inversion breakpoints are more likely to occur at boundaries between TADs. Our results show that breakpoints occur at TAD boundaries more than expected by chance. Some breakpoints may alter gene expression within TADs supporting the hypothesis that position effects contribute to inversion establishment. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Analysis of three-dimensional chromatin packing domains by chromatin scanning transmission electron microscopy (ChromSTEM)

Chromatin organization over multiple length scales plays a critical role in the regulation of transcription. Deciphering the interplay of these processes requires high-resolution, three-dimensional, quantitative imaging of chromatin structure in vitro. Herein, we introduce ChromSTEM, a method that utilizes high-angle annular dark-field imaging and tomography in scanning transmission electron microscopy combined with DNA-specific staining for electron microscopy. We utilized ChromSTEM for an in-depth quantification of 3D chromatin conformation with high spatial resolution and contrast, allowing for characterization of higher-order chromatin structure almost down to the level of the DNA base pair. Employing mass scaling analysis on ChromSTEM mass density tomograms, we observed that chromatin forms spatially well-defined higher-order domains, around 80 nm in radius. Within domains, chromatin exhibits a polymeric fractal-like behavior and a radially decreasing mass-density from the center to the periphery. Unlike other nanoimaging and analysis techniques, we demonstrate that our unique combination of this high-resolution imaging technique with polymer physics-based analysis enables us to (i) investigate the chromatin conformation within packing domains and (ii) quantify statistical descriptors of chromatin structure that are relevant to transcription. We observe that packing domains have heterogeneous morphological properties even within the same cell line, underlying the potential role of statistical chromatin packing in regulating gene expression within eukaryotic nuclei.

Domain Model of Eukaryotic Genome Organization: From DNA Loops Fixed on the Nuclear Matrix to TADs

The article reviews the development of ideas on the domain organization of eukaryotic genome, with special attention on the studies of DNA loops anchored to the nuclear matrix and their role in the emergence of the modern model of eukaryotic genome spatial organization. Critical analysis of results demonstrating that topologically associated chromatin domains are structural-functional blocks of the genome supports the notion that these blocks are fundamentally different from domains whose existence was proposed by the domain hypothesis of eukaryotic genome organization formulated in the 1980s. Based on the discussed evidence, it is concluded that the model postulating that eukaryotic genome is built from uniformly organized structural-functional blocks has proven to be untenable.

Regulation of 3D Organization and Its Role in Cancer Biology

Three-dimensional (3D) genomics is the frontier field in the post-genomics era, its foremost content is the relationship between chromatin spatial conformation and regulation of gene transcription. Cancer biology is a complex system resulting from genetic alterations in key tumor oncogenes and suppressor genes for cell proliferation, DNA replication, cell differentiation, and homeostatic functions. Although scientific research in recent decades has revealed how the genome sequence is mutated in many cancers, high-order chromosomal structures involved in the development and fate of cancer cells represent a crucial but rarely explored aspect of cancer genomics. Hence, dissection of the 3D genome conformation of cancer helps understand the unique epigenetic patterns and gene regulation processes that distinguish cancer biology from normal physiological states. In recent years, research in tumor 3D genomics has grown quickly. With the rapid progress of 3D genomics technology, we can now better determine the relationship between cancer pathogenesis and the chromatin structure of cancer cells. It is becoming increasingly explicit that changes in 3D chromatin structure play a vital role in controlling oncogene transcription. This review focuses on the relationships between tumor gene expression regulation, tumor 3D chromatin structure, and cancer phenotypic plasticity. Furthermore, based on the functional consequences of spatial disorganization in the cancer genome, we look forward to the clinical application prospects of 3D genomic biomarkers.

TADfit is a multivariate linear regression model for profiling hierarchical chromatin domains on replicate Hi-C data

Topologically associating domains (TADs) are fundamental building blocks of three dimensional genome, and organized into complex hierarchies. Identifying hierarchical TADs on Hi-C data helps to understand the relationship between genome architectures and gene regulation. Herein we propose TADfit, a multivariate linear regression model for profiling hierarchical chromatin domains, which tries to fit the interaction frequencies in Hi-C contact matrix with and without replicates using all-possible hierarchical TADs, and the significant ones can be determined by the regression coefficients obtained with the help of an online learning solver called Follow-The-Regularized-Leader (FTRL). Beyond the existing methods, TADfit has an ability to handle multiple contact matrix replicates and find partially overlapping TADs on them, which helps to find the comprehensive underlying TADs across replicates from different experiments. The comparative results tell that TADfit has better accuracy and reproducibility, and the hierarchical TADs called by it exhibit a reasonable biological relevance.

TADfit is a computational method that can identify hierarchical or partially overlapping TADs from Hi-C data, in part by using information from multiple replicates to improve detection power.