Deletions of chromosomal regulatory boundaries are associated with congenital disease

An integrative TAD catalog in lymphoblastoid cell lines discloses the functional impact of deletions and insertions in human genomes

The human genome is packaged within a three-dimensional (3D) nucleus and organized into structural units known as compartments, topologically associating domains (TADs), and loops. TAD boundaries, separating adjacent TADs, have been found to be well conserved across mammalian species and more evolutionarily constrained than TADs themselves. Recent studies show that structural variants (SVs) can modify 3D genomes through the disruption of TADs, which play an essential role in insulating genes from outside regulatory elements' aberrant regulation. However, how SV affects the 3D genome structure and their association among different aspects of gene regulation and candidate cis-regulatory elements (cCREs) have rarely been studied systematically. Here, we assess the impact of SVs intersecting with TAD boundaries by developing an integrative Hi-C analysis pipeline, which enables the generation of an in-depth catalog of TADs and TAD boundaries in human lymphoblastoid cell lines (LCLs) to fill the gap of limited resources. Our catalog contains 18,865 TADs, including 4596 sub-TADs, with 185 SVs (TAD-SVs) that alter chromatin architecture. By leveraging the ENCODE registry of cCREs in humans, we determine that 34 of 185 TAD-SVs intersect with cCREs and observe significant enrichment of TAD-SVs within cCREs. This study provides a database of TADs and TAD-SVs in the human genome that will facilitate future investigations of the impact of SVs on chromatin structure and gene regulation in health and disease.

Structural Variants at the LMNB1 Locus: Deciphering Pathomechanisms in Autosomal Dominant Adult-Onset Demyelinating Leukodystrophy

OBJECTIVE

We aimed to elucidate the pathogenic mechanisms underlying autosomal dominant adult-onset demyelinating leukodystrophy (ADLD), and to understand the genotype/phenotype correlation of structural variants (SVs) in the LMNB1 locus.

BACKGROUND

Since the discovery of 3D genome architectures and topologically associating domains (TADs), new pathomechanisms have been postulated for SVs, regardless of gene dosage changes. ADLD is a rare genetic disease associated with duplications (classical ADLD) or noncoding deletions (atypical ADLD) in the LMNB1 locus.

METHODS

High-throughput chromosome conformation capture, RNA sequencing, histopathological analyses of postmortem brain tissues, and clinical and neuroradiological investigations were performed.

RESULTS

We collected data from >20 families worldwide carrying SVs in the LMNB1 locus and reported strong clinical variability, even among patients carrying duplications of the entire LMNB1 gene, ranging from classical and atypical ADLD to asymptomatic carriers. We showed that patients with classic ADLD always carried intra-TAD duplications, resulting in a simple gene dose gain. Atypical ADLD was caused by LMNB1 forebrain-specific misexpression due to inter-TAD deletions or duplications. The inter-TAD duplication, which extends centromerically and crosses the 2 TAD boundaries, did not cause ADLD. Our results provide evidence that astrocytes are key players in ADLD pathology.

INTERPRETATION

Our study sheds light on the 3D genome and TAD structural changes associated with SVs in the LMNB1 locus, and shows that a duplication encompassing LMNB1 is not sufficient per se to diagnose ADLD, thereby strongly affecting genetic counseling. Our study supports breaking TADs as an emerging pathogenic mechanism that should be considered when studying brain diseases. ANN NEUROL 2024;96:855-870.

Polymer Physics Models Reveal Structural Folding Features of Single-Molecule Gene Chromatin Conformations

Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2 Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.

The probability of chromatin to be at the nuclear lamina has no systematic effect on its transcription level in fruit flies

BACKGROUND

Multiple studies have demonstrated a negative correlation between gene expression and positioning of genes at the nuclear envelope (NE) lined by nuclear lamina, but the exact relationship remains unclear, especially in light of the highly stochastic, transient nature of the gene association with the NE.

RESULTS

In this paper, we ask whether there is a causal, systematic, genome-wide relationship between the expression levels of the groups of genes in topologically associating domains (TADs) of Drosophila nuclei and the probabilities of TADs to be found at the NE. To investigate the nature of this possible relationship, we combine a coarse-grained dynamic model of the entire Drosophila nucleus with genome-wide gene expression data; we analyze the TAD averaged transcription levels of genes against the probabilities of individual TADs to be in contact with the NE in the control and lamins-depleted nuclei. Our findings demonstrate that, within the statistical error margin, the stochastic positioning of Drosophila melanogaster TADs at the NE does not, by itself, systematically affect the mean level of gene expression in these TADs, while the expected negative correlation is confirmed. The correlation is weak and disappears completely for TADs not containing lamina-associated domains (LADs) or TADs containing LADs, considered separately. Verifiable hypotheses regarding the underlying mechanism for the presence of the correlation without causality are discussed. These include the possibility that the epigenetic marks and affinity to the NE of a TAD are determined by various non-mutually exclusive mechanisms and remain relatively stable during interphase.

CONCLUSIONS

At the level of TADs, the probability of chromatin being in contact with the nuclear envelope has no systematic, causal effect on the transcription level in Drosophila. The conclusion is reached by combining model-derived time-evolution of TAD locations within the nucleus with their experimental gene expression levels.

[Epigenetics, principles and examples of applications]

Since the discovery of DNA as the support of genetic information, the challenge for generations of life scientists was to understand the mechanisms underlying the process that translate the sequence of a gene to a phenotype. In the 1950s, the concept of epigenetics was defined by the British biologist Conrad H. Waddington as the study of "epigenesis" that governs the biological processes involved in the development of any organism. The term epigenetics, now best defined as "above the DNA sequence" reflects the gene-environment interactions by which genes determine traits. Since, its first description, studies underlying the mechanisms involved in these processes has led to an increasing understanding of the regulation all genome transactions such as transcription, replication, repair and the biological pathways coordinated by these mechanisms. We will discuss here the main principles regulating epigenetic processes, their roles in physiology, their evolution over the life time and their implications in medicine.

POSTRE: a tool to predict the pathological effects of human structural variants

Understanding the pathological impact of non-coding genetic variation is a major challenge in medical genetics. Accumulating evidences indicate that a significant fraction of genetic alterations, including structural variants (SVs), can cause human disease by altering the function of non-coding regulatory elements, such as enhancers. In the case of SVs, described pathomechanisms include changes in enhancer dosage and long-range enhancer-gene communication. However, there is still a clear gap between the need to predict and interpret the medical impact of non-coding variants, and the existence of tools to properly perform these tasks. To reduce this gap, we have developed POSTRE (Prediction Of STRuctural variant Effects), a computational tool to predict the pathogenicity of SVs implicated in a broad range of human congenital disorders. By considering disease-relevant cellular contexts, POSTRE identifies SVs with either coding or long-range pathological consequences with high specificity and sensitivity. Furthermore, POSTRE not only identifies pathogenic SVs, but also predicts the disease-causative genes and the underlying pathological mechanism (e.g, gene deletion, enhancer disconnection, enhancer adoption, etc.). POSTRE is available at https://github.com/vicsanga/Postre.

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.

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.

TADMaster: a comprehensive web-based tool for the analysis of topologically associated domains

Background

Chromosome conformation capture and its derivatives have provided substantial genetic data for understanding how chromatin self-organizes. These techniques have identified regions of high intrasequence interactions called topologically associated domains (TADs). TADs are structural and functional units that shape chromosomes and influence genomic expression. Many of these domains differ across cell development and can be impacted by diseases. Thus, analysis of the identified domains can provide insight into genome regulation. Hence, there are many approaches to identifying such domains across many cell lines. Despite the availability of multiple tools for TAD detection, TAD callers' speed, flexibility, result inconsistency, and reproducibility remain challenges in this research area.

Results

In this work, we developed a computational webserver called TADMaster that provides an analysis suite to directly evaluate the concordance level and robustness of two or more TAD data on any given genome region. The suite provides multiple visual and quantitative metrics to compare the identified domains' number, size, and various comparisons of shared domains, domain boundaries, and domain overlap.

Conclusions

TADMaster is an efficient and easy-to-use web application that provides a set of consensus and unique TADs to inform the choice of TADs. It can be accessed at http://tadmaster.io and is also available as a containerized application that can be deployed and run locally on any platform or operating system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-05020-2.

Disruption of the topologically associated domain at Xp21.2 is related to 46,XY gonadal dysgenesis

Background

Duplications at the Xp21.2 locus have previously been linked to 46,XY gonadal dysgenesis (GD), which is thought to result from gene dosage effects of NR0B1 (DAX1), but the exact disease mechanism remains unknown.

Methods

Patients with 46,XY GD were analysed by whole genome sequencing. Identified structural variants were confirmed by array CGH and analysed by high-throughput chromosome conformation capture (Hi-C).

Results

We identified two unrelated patients: one showing a complex rearrangement upstream of NR0B1 and a second harbouring a 1.2 Mb triplication, including NR0B1. Whole genome sequencing and Hi-C analysis revealed the rewiring of a topological-associated domain (TAD) boundary close to NR0B1 associated with neo-TAD formation and may cause enhancer hijacking and ectopic NR0B1 expression. Modelling of previous Xp21.2 structural variations associated with isolated GD support our hypothesis and predict similar neo-TAD formation as well as TAD fusion.

Conclusion

Here we present a general mechanism how deletions, duplications or inversions at the NR0B1 locus can lead to partial or complete GD by disrupting the cognate TAD in the vicinity of NR0B1. This model not only allows better diagnosis of GD with copy number variations (CNVs) at Xp21.2, but also gives deeper insight on how spatiotemporal activation of developmental genes can be disrupted by reorganised TADs causing impairment of gonadal development.

TADeus2: a web server facilitating the clinical diagnosis by pathogenicity assessment of structural variations disarranging 3D chromatin structure

In recent years great progress has been made in identification of structural variants (SV) in the human genome. However, the interpretation of SVs, especially located in non-coding DNA, remains challenging. One of the reasons stems in the lack of tools exclusively designed for clinical SVs evaluation acknowledging the 3D chromatin architecture. Therefore, we present TADeus2 a web server dedicated for a quick investigation of chromatin conformation changes, providing a visual framework for the interpretation of SVs affecting topologically associating domains (TADs). This tool provides a convenient visual inspection of SVs, both in a continuous genome view as well as from a rearrangement’s breakpoint perspective. Additionally, TADeus2 allows the user to assess the influence of analyzed SVs within flaking coding/non-coding regions based on the Hi-C matrix. Importantly, the SVs pathogenicity is quantified and ranked using TADA, ClassifyCNV tools and sampling-based P-value. TADeus2 is publicly available at https://tadeus2.mimuw.edu.pl.

InsuLock: A Weakly Supervised Learning Approach for Accurate Insulator Prediction, and Variant Impact Quantification

Mapping chromatin insulator loops is crucial to investigating genome evolution, elucidating critical biological functions, and ultimately quantifying variant impact in diseases. However, chromatin conformation profiling assays are usually expensive, time-consuming, and may report fuzzy insulator annotations with low resolution. Therefore, we propose a weakly supervised deep learning method, InsuLock, to address these challenges. Specifically, InsuLock first utilizes a Siamese neural network to predict the existence of insulators within a given region (up to 2000 bp). Then, it uses an object detection module for precise insulator boundary localization via gradient-weighted class activation mapping (~40 bp resolution). Finally, it quantifies variant impacts by comparing the insulator score differences between the wild-type and mutant alleles. We applied InsuLock on various bulk and single-cell datasets for performance testing and benchmarking. We showed that it outperformed existing methods with an AUROC of ~0.96 and condensed insulator annotations to ~2.5% of their original size while still demonstrating higher conservation scores and better motif enrichments. Finally, we utilized InsuLock to make cell-type-specific variant impacts from brain scATAC-seq data and identified a schizophrenia GWAS variant disrupting an insulator loop proximal to a known risk gene, indicating a possible new mechanism of action for the disease.

Towards a Better Vision of Retinoic Acid Signaling during Eye Development

Retinoic acid (RA) functions as an essential signal for development of the vertebrate eye by controlling the transcriptional regulatory activity of RA receptors (RARs). During eye development, the optic vesicles and later the retina generate RA as a metabolite of vitamin A (retinol). Retinol is first converted to retinaldehyde by retinol dehydrogenase 10 (RDH10) and then to RA by all three retinaldehyde dehydrogenases (ALDH1A1, ALDH1A2, and ALDH1A3). In early mouse embryos, RA diffuses to tissues throughout the optic placode, optic vesicle, and adjacent mesenchyme to stimulate folding of the optic vesicle to form the optic cup. RA later generated by the retina is needed for further morphogenesis of the optic cup and surrounding perioptic mesenchyme; loss of RA at this stage leads to microphthalmia and cornea plus eyelid defects. RA functions by binding to nuclear RARs at RA response elements (RAREs) that either activate or repress transcription of key genes. Binding of RA to RARs regulates recruitment of transcriptional coregulators such as nuclear receptor coactivator (NCOA) or nuclear receptor corepressor (NCOR), which in turn control binding of the generic coactivator p300 or the generic corepressor PRC2. No genes have been identified as direct targets of RA signaling during eye development, so future studies need to focus on identifying such genes and their RAREs. Studies designed to learn how RA normally controls eye development in vivo will provide basic knowledge valuable for determining how developmental eye defects occur and for improving strategies to treat eye defects.

Opportunity to improve livestock traits using 3D genomics

The advent of high-throughput chromosome conformation capture and sequencing (Hi-C) has enabled researchers to probe the 3D architecture of the mammalian genome in a genome-wide manner. Simultaneously, advances in epigenomic assays, such as chromatin immunoprecipitation and sequencing (ChIP-seq) and DNase-seq, have enabled researchers to study cis-regulatory interactions and chromatin accessibility across the same genome-wide scale. The use of these data has revealed many unique insights into gene regulation and disease pathomechanisms in several model organisms. With the advent of these high-throughput sequencing technologies, there has been an ever-increasing number of datasets available for study; however, this is often limited to model organisms. Livestock species play critical roles in the economies of developing and developed nations alike. Despite this, they are greatly underrepresented in the 3D genomics space; Hi-C and related technologies have the potential to revolutionise livestock breeding by enabling a more comprehensive understanding of how production traits are controlled. The growth in human and model organism Hi-C data has seen a surge in the availability of computational tools for use in 3D genomics, with some tools using machine learning techniques to predict features and improve dataset quality. In this review, we provide an overview of the 3D genome and discuss the status of 3D genomics in livestock before delving into advancing the field by drawing inspiration from research in human and mouse. We end by offering future directions for livestock research in the field of 3D genomics.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

Cytogenetic and Array-CGH Characterization of a Simple Case of Reciprocal t(3;10) Translocation Reveals a Hidden Deletion at 5q12

Chromosome deletions, including band 5q12, have rarely been reported and have been associated with a wide range of clinical manifestations, such as postnatal growth retardation, intellectual disability, hyperactivity, nonspecific ocular defects, facial dysmorphism, and epilepsy. In this study, we describe for the first time a child with growth retardation in which we identified a balanced t(3;10) translocation by conventional cytogenetic analysis in addition to an 8.6 Mb 5q12 deletion through array-CGH. Our results show that the phenotypic abnormalities of a case that had been interpreted as "balanced" by conventional cytogenetics are mainly due to a cryptic deletion, highlighting the need for molecular investigation in subjects with an abnormal phenotype before assuming the cause is an apparently simple cytogenetic rearrangement. Finally, we identify PDE4D and PIK3R1 genes as the two major candidates responsible for the clinical features expressed in our patient.

Purifying selection on noncoding deletions of human regulatory loci detected using their cellular pleiotropy

Genomic deletions provide a powerful loss-of-function model in noncoding regions to assess the role of purifying selection on genetic variation. Regulatory element function is characterized by nonuniform tissue and cell type activity, necessarily linking the study of fitness consequences from regulatory variants to their corresponding cellular activity. We generated a callset of deletions from genomes in the Alzheimer's Disease Neuroimaging Initiative (ADNI) and used deletions from The 1000 Genomes Project Consortium (1000GP) in order to examine whether purifying selection preserves noncoding sites of chromatin accessibility marked by DNase I hypersensitivity (DHS), histone modification (enhancer, transcribed, Polycomb-repressed, heterochromatin), and chromatin loop anchors. To examine this in a cellular activity-aware manner, we developed a statistical method, pleiotropy ratio score (PlyRS), which calculates a correlation-adjusted count of "cellular pleiotropy" for each noncoding base pair by analyzing shared regulatory annotations across tissues and cell types. By comparing real deletion PlyRS values to simulations in a length-matched framework and by using genomic covariates in analyses, we found that purifying selection acts to preserve both DHS and enhancer noncoding sites. However, we did not find evidence of purifying selection for noncoding transcribed, Polycomb-repressed, or heterochromatin sites beyond that of the noncoding background. Additionally, we found evidence that purifying selection is acting on chromatin loop integrity by preserving colocalized CTCF binding sites. At regions of DHS, enhancer, and CTCF within chromatin loop anchors, we found evidence that both sites of activity specific to a particular tissue or cell type and sites of cellularly pleiotropic activity are preserved by selection.

Epigenetics and Heart Development

Epigenetic control of gene expression during cardiac development and disease has been a topic of intense research in recent years. Advances in experimental methods to study DNA accessibility, transcription factor occupancy, and chromatin conformation capture technologies have helped identify regions of chromatin structure that play a role in regulating access of transcription factors to the promoter elements of genes, thereby modulating expression. These chromatin structures facilitate enhancer contacts across large genomic distances and function to insulate genes from cis-regulatory elements that lie outside the boundaries for the gene of interest. Changes in transcription factor occupancy due to changes in chromatin accessibility have been implicated in congenital heart disease. However, the factors controlling this process and their role in changing gene expression during development or disease remain unclear. In this review, we focus on recent advances in the understanding of epigenetic factors controlling cardiac morphogenesis and their role in diseases.

Novel MYT1 variants expose the complexity of oculo-auriculo-vertebral spectrum genetic mechanisms

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder characterized by anomalies mainly involving the structures derived from the first and second pharyngeal arches. The spectrum presents with heterogeneous clinical features and complex etiology with genetic factors not yet completely understood. To date, MYT1 is the most important gene unambiguously associated with the spectrum and with functional data confirmation. In this work, we aimed to identify new single nucleotide variants (SNVs) affecting MYT1 in a cohort of 73 Brazilian patients diagnosed with OAVS. In addition, we investigated copy number variations (CNVs) encompassing this gene or its cis-regulatory elements and compared the frequency of these events in patients versus a cohort of 455 Brazilian control individuals. A new SNV, predicted as likely deleterious, was identified in five unrelated patients with OAVS. All five patients presented hearing impairment and orbital asymmetry suggesting an association with the variant. CNVs near MYT1, located in its neighboring topologically associating domain (TAD), were found to be enriched in patients when compared to controls, indicating a possible involvement of this region with OAVS pathogenicity. Our findings highlight the genetic complexity of the spectrum that seems to involve more than one variant type and inheritance patterns.

Chromosomes in the genomic age. Preserving cytogenomic competence of diagnostic genome laboratories

Participation of clinical genetic laboratories in External Quality Assessment schemes (EQAs) is a powerful method to ascertain if any improvement or additional training is required in the diagnostic service. Here, we provide evidence from recent EQAs that the competence in recognizing and interpreting cytogenetic aberrations is variable and could impact patient management. We identify several trends that could affect cytogenomic competence. Firstly, as a result of the age distribution among clinical laboratory geneticists (CLGs) registered at the European Board of Medical Genetics, about 25-30% of those with experience in cytogenetics will retire during the next decade. At the same time, there are about twice as many molecular geneticists to cytogeneticists among the younger CLGs. Secondly, when surveying training programs for CLG, we observed that not all programs guarantee that candidates gather sufficient experience in clinical cytogenomics. Thirdly, we acknowledge that whole genome sequencing (WGS) has a great attraction to biomedical scientists that wish to enter a training program for CLG. This, with a larger number of positions available, makes a choice for specialization in molecular genetics logical. However, current WGS technology cannot provide a diagnosis in all cases. Understanding the etiology of chromosomal rearrangements is essential for appropriate follow-up and for ascertaining recurrence risks. We define the minimal knowledge a CLG should have about cytogenomics in a world dominated by WGS, and discuss how laboratory directors and boards of professional organizations in clinical genetics can uphold cytogenomic competence by providing adequate CLG training programs and attracting sufficient numbers of trainees.

The haplolethality paradox of the wupA gene in Drosophila

Haplolethals (HL) are regions of diploid genomes that in one dose are fatal for the organism. Their biological meaning is obscure because heterozygous loss-of-function mutations result in dominant lethality (DL) and, consequently, should be under strong negative selection. We report an in depth study of the HL associated to the gene wings up A (wupA). It encodes 13 transcripts (A-M) that yield 11 protein isoforms (A-K) of Troponin I (TnI). They are functionally diverse in their control of muscle contraction, cell polarity and cell proliferation. Isoform K transfers to the nucleus where it increases transcription of the cell proliferation related genes CDK2, CDK4, Rap and Rab5. The nuclear translocation of isoform K is prevented by the co-expression of A or B isoforms, which illustrates isoform interactions. The corresponding DL mutations are, either DNA rearrangements clustered towards the gene 3’ end, thus affecting the genomic organization of all transcripts, or CRISPR-induced mutations in one of the two ATG sites which eliminate a subset of wupA products. The joint elimination of isoforms C, F, G and H, however, do not cause DL phenotypes. Genetically driven expression of single isoforms rescue neither DL nor any of the mutants known in the gene, suggesting that normal function requires properly regulated expression of specific combinations, rather than single, TnI isoforms. We conclude that the wupA associated HL results from the combined haploinsufficiency of a large set of TnI isoforms. The qualitative and quantitative normal expression of which, requires the chromosomal integrity of the wupA genomic region. Since all fly TnI isoforms are encoded in the same gene, its HL condition becomes unavoidable. These wupA features are comparable to those of dpp, the only other HL studied to some extent, and reveal a scenario of strict dosage dependence with implications for gene expression regulation and splitting.

Most species contain two copies of their genetic endowment, each received from each progenitor. If one of the duplicated genes is non-functional, the remaining copy may supply enough product as to cover the requirements for normal function or, alternatively, may reflect the insufficiency through a visible phenotype. In rare occasions, however, having one copy is so deleterious that causes lethality. These so called “haplolethal regions”, exist across species and represent an evolutionary paradox since they should have been subject to intense negative selection. The inherent difficulties to study haplolethals have precluded their study so far. Here, we analyzed the case of one of the five haplolethal regions of Drosophila, the one associated to the Troponin I encoding gene wupA, by measuring the transcriptional effects of mutations and chromosomal rearrangements affecting this gene. The data show that this haplolethality results from the combined insufficiency of a large number of Troponin I isoforms, which are functionally specialized, show interference and require the integrity of the native chromatin structure for their quantitatively regulated expression. These features unveil novel aspects of gene expression and, possibly, on evolutionary gene splitting. Also, haplolethals underscore the biological significance of protein dosage, in particular for functionally related products.

A TAD Skeptic: Is 3D Genome Topology Conserved?

The notion that topologically associating domains (TADs) are highly conserved across species is prevalent in the field of 3D genomics. However, what exactly is meant by 'highly conserved' and what are the actual comparative data that support this notion? To address these questions, we performed a historical review of the relevant literature and retraced numerous citation chains to reveal the primary data that were used as the basis for the widely accepted conclusion that TADs are highly conserved across evolution. A thorough review of the available evidence suggests the answer may be more complex than what is commonly presented.

Dbx2 regulation in limbs suggests interTAD sharing of enhancers

BACKGROUND

During tetrapod limb development, the HOXA13 and HOXD13 transcription factors are critical for the emergence and organization of the autopod, the most distal aspect where digits will develop. Since previous work had suggested that the Dbx2 gene is a target of these factors, we set up to analyze in detail this potential regulatory interaction.

RESULTS

We show that HOX13 proteins bind to mammalian-specific sequences at the vicinity of the Dbx2 locus that have enhancer activity in developing digits. However, the functional inactivation of the DBX2 protein did not elicit any particular phenotype related to Hox genes inactivation in digits, suggesting either redundant or compensatory mechanisms. We report that the neighboring Nell2 and Ano6 genes are also expressed in distal limb buds and are in part controlled by the same Dbx2 enhancers despite being localized into two different topologically associating domains (TADs) flanking the Dbx2 locus.

CONCLUSIONS

We conclude that Hoxa13 and Hoxd genes cooperatively activate Dbx2 expression in developing digits through binding to mammalian specific regulatory sequences in the Dbx2 neighborhood. Furthermore, these enhancers can overcome TAD boundaries in either direction to co-regulate a set of genes located in distinct chromatin domains.

Electron transfer flavoprotein and its role in mitochondrial energy metabolism in health and disease

Electron transfer flavoprotein (ETF) is an enzyme with orthologs from bacteria to humans. Human ETF is nuclear encoded by two separate genes, ETFA and ETFB, respectively. After translation, the two subunits are imported to the mitochondrial matrix space and assemble into a heterodimer containing one FAD and one AMP as cofactors. ETF functions as a hub taking up electrons from at least 14 flavoenzymes, feeding them into the respiratory chain. This represents a major source of reducing power for the electron transport chain from fatty acid oxidation and amino acid degradation. Transfer of electrons from the donor enzymes to ETF occurs by direct transfer between the enzyme bound flavins, a process that is tightly regulated by the polypeptide chain and by protein:protein interactions. ETF, in turn relays electrons to the iron sulfur cluster of the inner membrane protein ETF:QO, from where they travel via the FAD in ETF:QO to ubiquinone, entering the respiratory chain at the level of complex III. ETF recognizes its dehydrogenase partners via a recognition loop that anchors the protein on its partner followed by dynamic movements of the ETF flavin domain that bring redox cofactors in close proximity, thus promoting electron transfer. Genetic mutations in the ETFA or ETFB genes cause the Mendelian disorder multiple acyl-CoA dehydrogenase deficiency (MADD; OMIM #231680). We here review the knowledge on human ETF and investigations of the effects of disease-associated missense mutations in this protein that have promoted the understanding of the essential role that ETF plays in cellular metabolism and human disease.

Independence of chromatin conformation and gene regulation during Drosophila dorsoventral patterning

The relationship between chromatin organization and gene regulation remains unclear. While disruption of chromatin domains and domain boundaries can lead to misexpression of developmental genes, acute depletion of regulators of genome organization has a relatively small effect on gene expression. It is therefore uncertain whether gene expression and chromatin state drive chromatin organization or whether changes in chromatin organization facilitate cell-type-specific activation of gene expression. Here, using the dorsoventral patterning of the Drosophila melanogaster embryo as a model system, we provide evidence for the independence of chromatin organization and dorsoventral gene expression. We define tissue-specific enhancers and link them to expression patterns using single-cell RNA-seq. Surprisingly, despite tissue-specific chromatin states and gene expression, chromatin organization is largely maintained across tissues. Our results indicate that tissue-specific chromatin conformation is not necessary for tissue-specific gene expression but rather acts as a scaffold facilitating gene expression when enhancers become active.

A Guide to Using ClinTAD for Interpretation of DNA Copy Number Variants in the Context of Topologically Associated Domains

DNA copy number variants (CNVs) are routinely evaluated as part of clinical diagnosis in both the prenatal and postnatal genetic settings. Current guidelines for interpreting the potential clinical significance of these CNVs, typically identified by chromosomal microarray, focus entirely on genes localized within the CNV region. However, recent work has suggested that some CNVs can actually produce clinical impacts by influencing transcription of genes outside the CNV region. These alterations of transcription appear to occur by disrupting the composition of DNA topologically associated domains (TADs), which strongly influence contacts between gene promoters and their associated enhancers. Here we present a set of detailed protocols for the use of the free software tool ClinTAD (https://www.clintad.com). This decision-support software allows for prediction as to whether a given CNV may potentially disrupt a TAD boundary, and offers phenotype matching to genes near, but not within the CNV region, whose expression could be influenced by altered TAD architecture and that have phenotypic impacts related to that reported in a given patient. Our protocols here provide specific examples of how to implement these tools. In addition, the software has the capability to impact genomic research by evaluating multiple cases in parallel. We propose that this decision-support tool can benefit and improve genetic diagnosis. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Evaluating a single case using ClinTAD Basic Protocol 2: Evaluating a single case with multiple variants using ClinTAD Basic Protocol 3: Evaluating multiple cases using ClinTAD Basic Protocol 4: Creating tracks with custom data.

Lymphedema distichiasis syndrome may be caused by FOXC2 promoter-enhancer dissociation and disruption of a topological associated domain

Lymphedema distichiasis syndrome (LDS) is a rare autosomal dominant condition characterized by lower limb lymphedema, distichiasis, and variable additional features. LDS is usually caused by heterozygous sequence variants in the FOXC2 gene located at 16q24, but in one previous instance LDS has resulted from a balanced reciprocal translocation with a breakpoint at 16q24, 120 kb distal to the FOXC2 gene suggesting a position effect. Here, we describe a second family with LDS caused by a translocation involving 16q24. The family were ascertained after detection of a paternally inherited balanced reciprocal translocation t(16;22)(q24;q13.1) in a pregnancy complicated by severe fetal hydrops. There was a past history of multiple miscarriages in the father's family, and a personal and family history of lymphedema and distichiasis, consistent with the diagnosis of LDS. Using whole genome amplified DNA from single sperm of the male proband, bead array analysis demonstrated that the FOXC2 gene was intact and the chromosome 16 breakpoint mapped to the same region 120Kb distal to the FOXC2 gene. This case highlights the clinical consequences that can arise from a translocation of genomic material without dosage imbalance, and that it is increasingly feasible to predict and characterize possible effects with improved access to molecular techniques.

Prioritizing long range interactions in noncoding regions using GWAS and deletions perturbed TADs

Genome-wide association studies (GWAS) have contributed significantly to predisposing the disease etiology by associating single nucleotide polymorphisms (SNPs) with complex diseases. However, most GWAS-SNPs are in the noncoding regions that may affect distal genes via long range enhancer-promoter interactions. Thus, the common practice on GWAS discoveries cannot fully reveal the molecular mechanisms underpinning complex diseases. It is known that perturbations of topological associated domains (TADs) lead to long range interactions which underlie disease etiology. To identify the probable long range interactions in noncoding regions via GWAS and TADs perturbed by deletions, we integrated datasets from GWAS-SNPs, enhancers, TADs, and deletions. After ranking and clustering, we prioritized 201,132 high confident pairs of GWAS-SNPs and target genes. In this study, we performed a systematic inference on noncoding regions via GWAS-SNPs and deletion-perturbed TADs to boost GWAS discovery power. The high confident pairs of GWAS-SNPs and target genes (SE-Gs) provide the promising candidates to understand the molecular mechanisms underlying complex diseases with emphasis on the three-dimensional genome.

Rare or Overlooked? Structural Disruption of Regulatory Domains in Human Neurocristopathies

In the last few years, the role of non-coding regulatory elements and their involvement in human disease have received great attention. Among the non-coding regulatory sequences, enhancers are particularly important for the proper establishment of cell type-specific gene-expression programs. Furthermore, the disruption of enhancers can lead to human disease through two main mechanisms: (i) Mutations or copy number variants can directly alter the enhancer sequences and thereby affect expression of their target genes; (ii) structural variants can provoke changes in 3-D chromatin organization that alter neither the enhancers nor their target genes, but rather the physical communication between them. In this review, these pathomechanisms are mostly discussed in the context of neurocristopathies, congenital disorders caused by defects that occur during neural crest development. We highlight why, due to its contribution to multiple tissues and organs, the neural crest represents an important, yet understudied, cell type involved in multiple congenital disorders. Moreover, we discuss currently available resources and experimental models for the study of human neurocristopathies. Last, we provide some practical guidelines that can be followed when investigating human neurocristopathies caused by structural variants. Importantly, these guidelines can be useful not only to uncover the etiology of human neurocristopathies, but also of other human congenital disorders in which enhancer disruption is involved.

The preconception environment and sperm epigenetics

BACKGROUND

Infertility is a common reproductive disorder, with male factor infertility accounting for approximately half of all cases. Taking a paternal perceptive, recent research has shown that sperm epigenetics, such as changes in DNA methylation, histone modification, chromatin structure, and noncoding RNA expression, can impact reproductive and offspring health. Importantly, environmental conditions during the preconception period has been demonstrated to shape sperm epigenetics.

OBJECTIVES

To provide an overview on epigenetic modifications that regulate normal gene expression and epigenetic remodeling that occurs during spermatogenesis, and to discuss the epigenetic alterations that may occur to the paternal germline as a consequence of preconception environmental conditions and exposures.

MATERIALS AND METHODS

We examined published literature available on databases (PubMed, Google Scholar, ScienceDirect) focusing on adult male preconception environmental exposures and sperm epigenetics in epidemiologic studies and animal models.

RESULTS

The preconception period is a sensitive developmental window in which a variety of exposures such as toxicants, nutrition, drugs, stress, and exercise, affects sperm epigenetics.

DISCUSSION AND CONCLUSION

Understanding the environmental legacy of the sperm epigenome during spermatogenesis will enhance our understanding of reproductive health and improve reproductive success and offspring well-being.

Hi-C Identifies Complex Genomic Rearrangements and TAD-Shuffling in Developmental Diseases

Genome-wide analysis methods, such as array comparative genomic hybridization (CGH) and whole-genome sequencing (WGS), have greatly advanced the identification of structural variants (SVs) in the human genome. However, even with standard high-throughput sequencing techniques, complex rearrangements with multiple breakpoints are often difficult to resolve, and predicting their effects on gene expression and phenotype remains a challenge. Here, we address these problems by using high-throughput chromosome conformation capture (Hi-C) generated from cultured cells of nine individuals with developmental disorders (DDs). Three individuals had previously been identified as harboring duplications at the SOX9 locus and six had been identified with translocations. Hi-C resolved the positions of the duplications and was instructive in interpreting their distinct pathogenic effects, including the formation of new topologically associating domains (neo-TADs). Hi-C was very sensitive in detecting translocations, and it revealed previously unrecognized complex rearrangements at the breakpoints. In several cases, we observed the formation of fused-TADs promoting ectopic enhancer-promoter interactions that were likely to be involved in the disease pathology. In summary, we show that Hi-C is a sensible method for the detection of complex SVs in a clinical setting. The results help interpret the possible pathogenic effects of the SVs in individuals with DDs.

TADCompare: An R Package for Differential and Temporal Analysis of Topologically Associated Domains

Recent research using chromatin conformation capture technologies, such as Hi-C, has demonstrated the importance of topologically associated domains (TADs) and smaller chromatin loops, collectively referred hereafter as "interacting domains." Many such domains change during development or disease, and exhibit cell- and condition-specific differences. Quantification of the dynamic behavior of interacting domains will help to better understand genome regulation. Methods for comparing interacting domains between cells and conditions are highly limited. We developed TADCompare, a method for differential analysis of boundaries of interacting domains between two or more Hi-C datasets. TADCompare is based on a spectral clustering-derived measure called the eigenvector gap, which enables a loci-by-loci comparison of boundary differences. Using this measure, we introduce methods for identifying differential and consensus boundaries of interacting domains and tracking boundary changes over time. We further propose a novel framework for the systematic classification of boundary changes. Colocalization- and gene enrichment analysis of different types of boundary changes demonstrated distinct biological functionality associated with them. TADCompare is available on https://github.com/dozmorovlab/TADCompare and Bioconductor (submitted).

Disruption of chromatin folding domains by somatic genomic rearrangements in human cancer

Chromatin is folded into successive layers to organize linear DNA. Genes within the same topologically associating domains (TADs) demonstrate similar expression and histone-modification profiles, and boundaries separating different domains have important roles in reinforcing the stability of these features. Indeed, domain disruptions in human cancers can lead to misregulation of gene expression. However, the frequency of domain disruptions in human cancers remains unclear. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), which aggregated whole-genome sequencing data from 2,658 cancers across 38 tumor types, we analyzed 288,457 somatic structural variations (SVs) to understand the distributions and effects of SVs across TADs. Notably, SVs can lead to the fusion of discrete TADs, and complex rearrangements markedly change chromatin folding maps in the cancer genomes. Notably, only 14% of the boundary deletions resulted in a change in expression in nearby genes of more than twofold.

Clustered CTCF binding is an evolutionary mechanism to maintain topologically associating domains

BACKGROUND

CTCF binding contributes to the establishment of a higher-order genome structure by demarcating the boundaries of large-scale topologically associating domains (TADs). However, despite the importance and conservation of TADs, the role of CTCF binding in their evolution and stability remains elusive.

RESULTS

We carry out an experimental and computational study that exploits the natural genetic variation across five closely related species to assess how CTCF binding patterns stably fixed by evolution in each species contribute to the establishment and evolutionary dynamics of TAD boundaries. We perform CTCF ChIP-seq in multiple mouse species to create genome-wide binding profiles and associate them with TAD boundaries. Our analyses reveal that CTCF binding is maintained at TAD boundaries by a balance of selective constraints and dynamic evolutionary processes. Regardless of their conservation across species, CTCF binding sites at TAD boundaries are subject to stronger sequence and functional constraints compared to other CTCF sites. TAD boundaries frequently harbor dynamically evolving clusters containing both evolutionarily old and young CTCF sites as a result of the repeated acquisition of new species-specific sites close to conserved ones. The overwhelming majority of clustered CTCF sites colocalize with cohesin and are significantly closer to gene transcription start sites than nonclustered CTCF sites, suggesting that CTCF clusters particularly contribute to cohesin stabilization and transcriptional regulation.

CONCLUSIONS

Dynamic conservation of CTCF site clusters is an apparently important feature of CTCF binding evolution that is critical to the functional stability of a higher-order chromatin structure.

A novel 14q13.1-21.1 deletion identified by CNV-Seq in a patient with brain-lung-thyroid syndrome, tooth agenesis and immunodeficiency

Background

Chromosome 14q11-q22 deletion syndrome (OMIM 613457) is a rare genomic disorder. The phenotype heterogeneity depends on the deletion size, breakpoints and genes deleted. Critical genes like FOXG1, NKX2–1, PAX9 were identified.

Case presentation

We performed whole exome sequencing (WES) and copy number variation sequencing (CNV-seq) for a patient with mild speech and motor developmental delay, short stature, recurrent pulmonary infections, tooth agenesis and triad of brain-lung-thyroid syndrome. By using CNV-seq, we identified a 3.1 Mb de novo interstitial deletion of the 14q13.2q21.1 region encompassing 17 OMIM genes including NKX2–1, PAX9 and NFKBIA. Our patient’s phenotype is consistent with other published 14q13 deletion patients.

Conclusion

Our results showed the combination of WES and CNV-seq is an effective diagnostic strategy for patients with genetic or genomic disorders. After reviewing published patients, we also proposed a new critical region for 14q13 deletion syndrome with is a more benign disorder compared to 14q11-q22 deletion syndrome.

The distributions of protein coding genes within chromatin domains in relation to human disease

Background

Our understanding of the nuclear chromatin structure has increased hugely during the last years mainly as a consequence of the advances in chromatin conformation capture methods like Hi-C. The unprecedented resolution of genome-wide interaction maps shows functional consequences that extend the initial thought of an efficient DNA packaging mechanism: gene regulation, DNA repair, chromosomal translocations and evolutionary rearrangements seem to be only the peak of the iceberg. One key concept emerging from this research is the topologically associating domains (TADs) whose functional role in gene regulation and their association with disease is not fully untangled.

Results

We report that the lower the number of protein coding genes inside TADs, the higher the tendency of those genes to be associated with disease (p-value = 4 × \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-54}$$\end{document}10-54). Moreover, housekeeping genes are less associated with disease than other genes. Accordingly, they are depleted in TADs containing less than three protein coding genes (p-value = 3.9 × \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-34}$$\end{document}10-34). We observed that TADs with higher ratios of enhancers versus genes contained higher numbers of disease-associated genes. We interpret these results as an indication that sharing enhancers among genes reduces their involvement in disease. Larger TADs would have more chances to accommodate many genes and select for enhancer sharing along evolution.

Conclusions

Genes associated with human disease do not distribute randomly over the TADs. Our observations suggest general rules that confer functional stability to TADs, adding more evidence to the role of TADs as regulatory units.

Prioritization of genes driving congenital phenotypes of patients with de novo genomic structural variants

BACKGROUND

Genomic structural variants (SVs) can affect many genes and regulatory elements. Therefore, the molecular mechanisms driving the phenotypes of patients carrying de novo SVs are frequently unknown.

METHODS

We applied a combination of systematic experimental and bioinformatic methods to improve the molecular diagnosis of 39 patients with multiple congenital abnormalities and/or intellectual disability harboring apparent de novo SVs, most with an inconclusive diagnosis after regular genetic testing.

RESULTS

In 7 of these cases (18%), whole-genome sequencing analysis revealed disease-relevant complexities of the SVs missed in routine microarray-based analyses. We developed a computational tool to predict the effects on genes directly affected by SVs and on genes indirectly affected likely due to the changes in chromatin organization and impact on regulatory mechanisms. By combining these functional predictions with extensive phenotype information, candidate driver genes were identified in 16/39 (41%) patients. In 8 cases, evidence was found for the involvement of multiple candidate drivers contributing to different parts of the phenotypes. Subsequently, we applied this computational method to two cohorts containing a total of 379 patients with previously detected and classified de novo SVs and identified candidate driver genes in 189 cases (50%), including 40 cases whose SVs were previously not classified as pathogenic. Pathogenic position effects were predicted in 28% of all studied cases with balanced SVs and in 11% of the cases with copy number variants.

CONCLUSIONS

These results demonstrate an integrated computational and experimental approach to predict driver genes based on analyses of WGS data with phenotype association and chromatin organization datasets. These analyses nominate new pathogenic loci and have strong potential to improve the molecular diagnosis of patients with de novo SVs.

Chromatin organization modulates the origin of heritable structural variations in human genome

Structural variations (SVs) in the human genome originate from different mechanisms related to DNA repair, replication errors, and retrotransposition. Our analyses of 26 927 SVs from the 1000 Genomes Project revealed differential distributions and consequences of SVs of different origin, e.g. deletions from non-allelic homologous recombination (NAHR) are more prone to disrupt chromatin organization while processed pseudogenes can create accessible chromatin. Spontaneous double stranded breaks (DSBs) are the best predictor of enrichment of NAHR deletions in open chromatin. This evidence, along with strong physical interaction of NAHR breakpoints belonging to the same deletion suggests that majority of NAHR deletions are non-meiotic i.e. originate from errors during homology directed repair (HDR) of spontaneous DSBs. In turn, the origin of the spontaneous DSBs is associated with transcription factor binding in accessible chromatin revealing the vulnerability of functional, open chromatin. The chromatin itself is enriched with repeats, particularly fixed Alu elements that provide the homology required to maintain stability via HDR. Through co-localization of fixed Alus and NAHR deletions in open chromatin we hypothesize that old Alu expansion had a stabilizing role on the human genome.

Epigenetic Modifications in Stress Response Genes Associated With Childhood Trauma

Adverse childhood experiences (ACEs) may be referred to by other terms (e.g., early life adversity or stress and childhood trauma) and have a lifelong impact on mental and physical health. For example, childhood trauma has been associated with posttraumatic stress disorder (PTSD), anxiety, depression, bipolar disorder, diabetes, and cardiovascular disease. The heritability of ACE-related phenotypes such as PTSD, depression, and resilience is low to moderate, and, moreover, is very variable for a given phenotype, which implies that gene by environment interactions (such as through epigenetic modifications) may be involved in the onset of these phenotypes. Currently, there is increasing interest in the investigation of epigenetic contributions to ACE-induced differential health outcomes. Although there are a number of studies in this field, there are still research gaps. In this review, the basic concepts of epigenetic modifications (such as methylation) and the function of the hypothalamic-pituitary-adrenal (HPA) axis in the stress response are outlined. Examples of specific genes undergoing methylation in association with ACE-induced differential health outcomes are provided. Limitations in this field, e.g., uncertain clinical diagnosis, conceptual inconsistencies, and technical drawbacks, are reviewed, with suggestions for advances using new technologies and novel research directions. We thereby provide a platform on which the field of ACE-induced phenotypes in mental health may build.

Functional Consequences of Chromosomal Rearrangements on Gene Expression: Not So Deleterious After All?

Chromosomes are folded and organized into topologically associating domains (TADs) which provide a framework for the interaction of enhancers with the promoter of their target gene(s). Structural rearrangements observed during evolution or in disease contexts suggest that changes in genome organization strongly affect gene expression and can have drastic phenotypic effects. In this review, I will discuss how recent genomic engineering experiments reveal a more contrasted picture, suggesting that TADs are important but not always essential for gene expression regulation.

Visualising three-dimensional genome organisation in two dimensions

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

Transcriptionally active HERV-H retrotransposons demarcate topologically associating domains in human pluripotent stem cells

Chromatin architecture has been implicated in cell type-specific gene regulatory programs, yet how chromatin remodels during development remains to be fully elucidated. Here, by interrogating chromatin reorganization during human pluripotent stem cell (hPSC) differentiation, we discover a role for the primate-specific endogenous retrotransposon human endogenous retrovirus subfamily H (HERV-H) in creating topologically associating domains (TADs) in hPSCs. Deleting these HERV-H elements eliminates their corresponding TAD boundaries and reduces the transcription of upstream genes, while de novo insertion of HERV-H elements can introduce new TAD boundaries. The ability of HERV-H to create TAD boundaries depends on high transcription, as transcriptional repression of HERV-H elements prevents the formation of boundaries. This ability is not limited to hPSCs, as these actively transcribed HERV-H elements and their corresponding TAD boundaries also appear in pluripotent stem cells from other hominids but not in more distantly related species lacking HERV-H elements. Overall, our results provide direct evidence for retrotransposons in actively shaping cell type- and species-specific chromatin architecture.

Structural variation in the 3D genome

Structural and quantitative chromosomal rearrangements, collectively referred to as structural variation (SV), contribute to a large extent to the genetic diversity of the human genome and thus are of high relevance for cancer genetics, rare diseases and evolutionary genetics. Recent studies have shown that SVs can not only affect gene dosage but also modulate basic mechanisms of gene regulation. SVs can alter the copy number of regulatory elements or modify the 3D genome by disrupting higher-order chromatin organization such as topologically associating domains. As a result of these position effects, SVs can influence the expression of genes distant from the SV breakpoints, thereby causing disease. The impact of SVs on the 3D genome and on gene expression regulation has to be considered when interpreting the pathogenic potential of these variant types.

Understanding the 3D genome: Emerging impacts on human disease

Recent burst of new technologies that allow for quantitatively delineating chromatin structure has greatly expanded our understanding of how the genome is organized in the three-dimensional (3D) space of the nucleus. It is now clear that the hierarchical organization of the eukaryotic genome critically impacts nuclear activities such as transcription, replication, as well as cellular and developmental events such as cell cycle, cell fate decision and embryonic development. In this review, we discuss new insights into how the structural features of the 3D genome hierarchy are established and maintained, how this hierarchy undergoes dynamic rearrangement during normal development and how its perturbation will lead to human disease, highlighting the accumulating evidence that links the diverse 3D genome architecture components to a multitude of human diseases and the emerging mechanisms by which 3D genome derangement causes disease phenotypes.

Navigating the non-coding genome in heart development and Congenital Heart Disease

Congenital Heart Disease (CHD) is characterised by a wide range of cardiac defects, from mild to life-threatening, which occur in babies worldwide. To date, there is no cure to CHD, however, progress in surgery has reduced its mortality allowing children affected by CHD to reach adulthood. In an effort to understand its genetic basis, several studies involving whole-genome sequencing (WGS) of patients with CHD have been undertaken and generated a great wealth of information. The majority of putative causative mutations identified in WGS studies fall into the non-coding part of the genome. Unfortunately, due to the lack of understanding of the function of these non-coding mutations, it is challenging to establish a causal link between the non-coding mutation and the disease. Thus, here we review the state-of-the-art approaches to interpret non-coding mutations in the context of CHD and address the following questions: What are the non-coding sequences important for cardiac function? Which technologies are used to identify them? Which resources are available to analyse them? What mutations are expected in these non-coding sequences? Learning from developmental process, what is their expected role in CHD?

Location, location, location: nuclear structure regulates gene expression in neurons

Genome architecture plays a critical role in regulating the expression of genes that are essential for nervous system development. During neuronal differentiation, spatially and temporally regulated transcription allows neuronal migration, the growth of dendrites and axons, and at later stages, synaptic formation and the establishment of neuronal circuitry. Genome topology and relocation of gene loci within the nucleus are now regarded as key factors that contribute to transcriptional regulation. Here, we review recent work supporting the hypothesis that the dynamic organization of chromatin within the nucleus impacts gene activation in response to extrinsic signalling and during neuronal differentiation. The consequences of disruption of the genome architecture on neuronal health will be also discussed.

ClinTAD: a tool for copy number variant interpretation in the context of topologically associated domains

Standard clinical interpretation of DNA copy number variants (CNVs) identified by cytogenomic microarray involves examining protein-coding genes within the region and comparison to other CNVs. Emerging basic research suggests that CNVs can also exert a pathogenic effect through disruption of DNA structural elements such as topologically associated domains (TADs). To begin to integrate these discoveries with current practice, we developed ClinTAD, a free browser-based tool to assist with interpretation of CNVs in the context of TADs ( www.clintad.com ). We used ClinTAD to examine 209 randomly selected single-nucleotide polymorphism microarray cases with a total of 236 CNVs. We compared 118 CNVs classified as variants of uncertain clinical significance (VUS), where additional insight into pathogenicity of these CNVs would be of greatest utility, to 118 CNVs classified as benign. We found that a higher proportion of VUS had at least two genes in a nearby TAD related to a phenotype seen in the patient based on Human Phenotype Ontology (HPO) annotation. We present example cases demonstrating scenarios where ClinTAD may either increase or decrease clinical suspicion of pathogenicity for VUS, depending on disruption of TAD boundaries and HPO phenotype match. ClinTAD is an easy-to-use tool, based on emerging research in chromatin architecture, that can help inform CNV interpretation.

Chromatin features constrain structural variation across evolutionary timescales

Noncoding DNA sequences play crucial roles in gene regulation, including via three-dimensional genome organization where they define chromatin boundaries and segment the genome into a sequence of insulated neighborhoods. However, the relative importance of noncoding DNA elements, particularly in comparison with protein-coding DNA sequences, remains more poorly characterized. Here, we systematically test if chromatin boundary disruptions are under purifying selection. Our analyses uncover a genomewide depletion of structural variants that would have the potential to alter chromatin structure. This in turn has implications for predicting not only which variants are likely pathogenic in clinical genetics settings, but also which are likely key innovations in primate evolution, and argues for expanding the current gene-centric paradigm for interpreting structural variants.

The potential impact of structural variants includes not only the duplication or deletion of coding sequences, but also the perturbation of noncoding DNA regulatory elements and structural chromatin features, including topological domains (TADs). Structural variants disrupting TAD boundaries have been implicated both in cancer and developmental disease; this likely occurs via “enhancer hijacking,” whereby removal of the TAD boundary exposes enhancers to new target transcription start sites (TSSs). With this functional role, we hypothesized that boundaries would display evidence for negative selection. Here we demonstrate that the chromatin landscape constrains structural variation both within healthy humans and across primate evolution. In contrast, in patients with developmental delay, variants occur remarkably uniformly across genomic features, suggesting a potentially broad role for enhancer hijacking in human disease.