Topologically associating domains are stable units of replication-timing regulation

3D disorganization and rearrangement of genome provide insights into pathogenesis of NAFLD by integrated Hi-C, Nanopore, and RNA sequencing

The three-dimensional (3D) conformation of chromatin is integral to the precise regulation of gene expression. The 3D genome and genomic variations in non-alcoholic fatty liver disease (NAFLD) are largely unknown, despite their key roles in cellular function and physiological processes. High-throughput chromosome conformation capture (Hi-C), Nanopore sequencing, and RNA-sequencing (RNA-seq) assays were performed on the liver of normal and NAFLD mice. A high-resolution 3D chromatin interaction map was generated to examine different 3D genome hierarchies including A/B compartments, topologically associated domains (TADs), and chromatin loops by Hi-C, and whole genome sequencing identifying structural variations (SVs) and copy number variations (CNVs) by Nanopore sequencing. We identified variations in thousands of regions across the genome with respect to 3D chromatin organization and genomic rearrangements, between normal and NAFLD mice, and revealed gene dysregulation frequently accompanied by these variations. Candidate target genes were identified in NAFLD, impacted by genetic rearrangements and spatial organization disruption. Our data provide a high-resolution 3D genome interaction resource for NAFLD investigations, revealed the relationship among genetic rearrangements, spatial organization disruption, and gene regulation, and identified candidate genes associated with these variations implicated in the pathogenesis of NAFLD. The newly findings offer insights into novel mechanisms of NAFLD pathogenesis and can provide a new conceptual framework for NAFLD therapy.

Visualizing Live Chromatin Dynamics through CRISPR-Based Imaging Techniques

The three-dimensional organization of chromatin and its time-dependent changes greatly affect virtually every cellular function, especially DNA replication, genome maintenance, transcription regulation, and cell differentiation. Sequencing-based techniques such as ChIP-seq, ATAC-seq, and Hi-C provide abundant information on how genomic elements are coupled with regulatory proteins and functionally organized into hierarchical domains through their interactions. However, visualizing the time-dependent changes of such organization in individual cells remains challenging. Recent developments of CRISPR systems for site-specific fluorescent labeling of genomic loci have provided promising strategies for visualizing chromatin dynamics in live cells. However, there are several limiting factors, including background signals, off-target binding of CRISPR, and rapid photobleaching of the fluorophores, requiring a large number of target-bound CRISPR complexes to reliably distinguish the target-specific foci from the background. Various modifications have been engineered into the CRISPR system to enhance the signal-to-background ratio and signal longevity to detect target foci more reliably and efficiently, and to reduce the required target size. In this review, we comprehensively compare the performances of recently developed CRISPR designs for improved visualization of genomic loci in terms of the reliability of target detection, the ability to detect small repeat loci, and the allowed time of live tracking. Longer observation of genomic loci allows the detailed identification of the dynamic characteristics of chromatin. The diffusion properties of chromatin found in recent studies are reviewed, which provide suggestions for the underlying biological processes.

Hierarchical chromatin features reveal the toxin production in Bungarus multicinctus

BACKGROUND

Bungarus multicinctus, from which a classical Chinese medicine is produced, is known as the most venomous land snake in the world, but the chromatin organization and transcription factor activity during venom replenishment progress have not been explored yet. This study aimed to determine the roles of chromatin structure in toxin activity via bioinformatics and experimental validation.

METHODS

Chromosome conformation capture (Hi-C) analysis was used to examine interactions among chromosomes and identify different scales of chromatin during envenomation in B. multicinctus. Correlations between epigenetic modifications and chromatin structure were verified through ChIP-seq analysis. RNA-seq was used to validate the influence of variation in chromatin structure and gene expression levels on venom production and regulation.

RESULTS

Our results suggested that intra-chromosomal interactions are more intense than inter-chromosomal interactions among the control group, 3-day group of venom glands and muscles. Through this, we found that compartmental transition was correlated with chromatin interactions. Interestingly, the up-regulated genes in more compartmental switch regions reflect the function of toxin activity. Topologically associated domain (TAD) boundaries enriched with histone modifications are associated with different distributions of genes and the expression levels. Toxin-coding genes in the same loop are highly expressed, implying that the importance of epigenetic regulation during envenomination. On a smaller scale, the epigenetic markers affect transcriptional regulation by controlling the recruitment/inhibition of transcription initiation complexes.

CONCLUSIONS

Chromatin structure and epigenetic modifications could play a vital status role in the mechanisms of venom regulation in B. multicinctus.

DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation

DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, reduced representation bisulfite sequencing (RRBS) profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, this study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

Spatial genome architecture and the emergence of malignancy

Human chromosomes are large spatially and hierarchically structured entities, the integrity of which needs to be preserved throughout the lifespan of the cell and in conjunction with cell cycle progression. Preservation of chromosomal structure is important for proper deployment of cell type-specific gene expression programs. Thus, aberrations in the integrity and structure of chromosomes will predictably lead to disease, including cancer. Here, we provide an updated standpoint with respect to chromatin misfolding and the emergence of various cancer types. We discuss recent studies implicating the disruption of topologically associating domains, switching between active and inactive compartments, rewiring of promoter-enhancer interactions in malignancy as well as the effects of single nucleotide polymorphisms in non-coding regions involved in long-range regulatory interactions. In light of these findings, we argue that chromosome conformation studies may now also be useful for patient diagnosis and drug target discovery.

Topological isolation of developmental regulators in mammalian genomes

Precise control of mammalian gene expression is facilitated through epigenetic mechanisms and nuclear organization. In particular, insulated chromosome structures are important for regulatory control, but the phenotypic consequences of their boundary disruption on developmental processes are complex and remain insufficiently understood. Here, we generated deeply sequenced Hi-C data for human pluripotent stem cells (hPSCs) that allowed us to identify CTCF loop domains that have highly conserved boundary CTCF sites and show a notable enrichment of individual developmental regulators. Importantly, perturbation of such a boundary in hPSCs interfered with proper differentiation through deregulated distal enhancer-promoter activity. Finally, we found that germline variations affecting such boundaries are subject to purifying selection and are underrepresented in the human population. Taken together, our findings highlight the importance of developmental gene isolation through chromosomal folding structures as a mechanism to ensure their proper expression.

The phenotypic consequence of 3D genome boundary disruption on developmental processes remains insufficiently understood. Here, the authors show that perturbation of a SOX17 boundary in human pluripotent stem cells interferes with proper differentiation and that germline variations affecting such boundaries are subject to selection, resulting in underrepresentation in the human population.

Motif-Raptor: a cell type-specific and transcription factor centric approach for post-GWAS prioritization of causal regulators

MOTIVATION

Genome-wide association studies (GWASs) have identified thousands of common trait-associated genetic variants but interpretation of their function remains challenging. These genetic variants can overlap the binding sites of transcription factors (TFs) and therefore could alter gene expression. However, we currently lack a systematic understanding on how this mechanism contributes to phenotype.

RESULTS

We present Motif-Raptor, a TF-centric computational tool that integrates sequence-based predictive models, chromatin accessibility, gene expression datasets and GWAS summary statistics to systematically investigate how TF function is affected by genetic variants. Given trait-associated non-coding variants, Motif-Raptor can recover relevant cell types and critical TFs to drive hypotheses regarding their mechanism of action. We tested Motif-Raptor on complex traits such as rheumatoid arthritis and red blood cell count and demonstrated its ability to prioritize relevant cell types, potential regulatory TFs and non-coding SNPs which have been previously characterized and validated.

AVAILABILITY AND IMPLEMENTATION

Motif-Raptor is freely available as a Python package at: https://github.com/pinellolab/MotifRaptor.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome maintenance during embryogenesis

Genome maintenance during embryogenesis is critical, because defects during this period can be perpetuated and thus have a long-term impact on individual's health and longevity. Nevertheless, genome instability is normal during certain aspects of embryonic development, indicating that there is a balance between the exigencies of timely cell proliferation and mutation prevention. In particular, early embryos possess unique cellular and molecular features that underscore the challenge of having an appropriate balance. Here, we discuss genome instability during embryonic development, the mechanisms used in various cell compartments to manage genomic stress and address outstanding questions regarding the balance between genome maintenance mechanisms in key cell types that are important for adulthood and progeny.

The nuclear periphery is a scaffold for tissue-specific enhancers

Nuclear architecture influences gene regulation and cell identity by controlling the three-dimensional organization of genes and their distal regulatory sequences, which may be far apart in linear space. The genome is functionally and spatially segregated in the eukaryotic nucleus with transcriptionally active regions in the nuclear interior separated from repressive regions, including those at the nuclear periphery. Here, we describe the identification of a novel type of nuclear peripheral chromatin domain that is enriched for tissue-specific transcriptional enhancers. Like other chromatin at the nuclear periphery, these regions are marked by H3K9me2. But unlike the nuclear peripheral Lamina-Associated Domains (LADs), these novel, enhancer-rich domains have limited Lamin B interaction. We therefore refer to them as H3K9me2-Only Domains (KODs). In mouse embryonic stem cells, KODs are found in Hi-C-defined A compartments and feature relatively accessible chromatin. KODs are characterized by low gene expression and enhancers located in these domains bear the histone marks of an inactive or poised state. These results indicate that KODs organize a subset of inactive, tissue-specific enhancers at the nuclear periphery. We hypothesize that KODs may play a role in facilitating and perhaps constraining the enhancer-promoter interactions underlying spatiotemporal regulation of gene expression programs in differentiation and development.

Transcription-coupled structural dynamics of topologically associating domains regulate replication origin efficiency

BACKGROUND

Metazoan cells only utilize a small subset of the potential DNA replication origins to duplicate the whole genome in each cell cycle. Origin choice is linked to cell growth, differentiation, and replication stress. Although various genetic and epigenetic signatures have been linked to the replication efficiency of origins, there is no consensus on how the selection of origins is determined.

RESULTS

We apply dual-color stochastic optical reconstruction microscopy (STORM) super-resolution imaging to map the spatial distribution of origins within individual topologically associating domains (TADs). We find that multiple replication origins initiate separately at the spatial boundary of a TAD at the beginning of the S phase. Intriguingly, while both high-efficiency and low-efficiency origins are distributed homogeneously in the TAD during the G1 phase, high-efficiency origins relocate to the TAD periphery before the S phase. Origin relocalization is dependent on both transcription and CTCF-mediated chromatin structure. Further, we observe that the replication machinery protein PCNA forms immobile clusters around TADs at the G1/S transition, explaining why origins at the TAD periphery are preferentially fired.

CONCLUSION

Our work reveals a new origin selection mechanism that the replication efficiency of origins is determined by their physical distribution in the chromatin domain, which undergoes a transcription-dependent structural re-organization process. Our model explains the complex links between replication origin efficiency and many genetic and epigenetic signatures that mark active transcription. The coordination between DNA replication, transcription, and chromatin organization inside individual TADs also provides new insights into the biological functions of sub-domain chromatin structural dynamics.

Multiscale modeling of genome organization with maximum entropy optimization

Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.

Orphan CpG islands amplify poised enhancer regulatory activity and determine target gene responsiveness

CpG islands (CGIs) represent a widespread feature of vertebrate genomes, being associated with ~70% of all gene promoters. CGIs control transcription initiation by conferring nearby promoters with unique chromatin properties. In addition, there are thousands of distal or orphan CGIs (oCGIs) whose functional relevance is barely known. Here we show that oCGIs are an essential component of poised enhancers that augment their long-range regulatory activity and control the responsiveness of their target genes. Using a knock-in strategy in mouse embryonic stem cells, we introduced poised enhancers with or without oCGIs within topologically associating domains harboring genes with different types of promoters. Analysis of the resulting cell lines revealed that oCGIs act as tethering elements that promote the physical and functional communication between poised enhancers and distally located genes, particularly those with large CGI clusters in their promoters. Therefore, by acting as genetic determinants of gene-enhancer compatibility, CGIs can contribute to gene expression control under both physiological and potentially pathological conditions.

Defining genome architecture at base-pair resolution

In higher eukaryotes, many genes are regulated by enhancers that are 104-106 base pairs (bp) away from the promoter. Enhancers contain transcription-factor-binding sites (which are typically around 7-22 bp), and physical contact between the promoters and enhancers is thought to be required to modulate gene expression. Although chromatin architecture has been mapped extensively at resolutions of 1 kilobase and above; it has not been possible to define physical contacts at the scale of the proteins that determine gene expression. Here we define these interactions in detail using a chromosome conformation capture method (Micro-Capture-C) that enables the physical contacts between different classes of regulatory elements to be determined at base-pair resolution. We find that highly punctate contacts occur between enhancers, promoters and CCCTC-binding factor (CTCF) sites and we show that transcription factors have an important role in the maintenance of the contacts between enhancers and promoters. Our data show that interactions between CTCF sites are increased when active promoters and enhancers are located within the intervening chromatin. This supports a model in which chromatin loop extrusion1 is dependent on cohesin loading at active promoters and enhancers, which explains the formation of tissue-specific chromatin domains without changes in CTCF binding.

Mammalian DNA Replication Timing

Immediately following the discovery of the structure of DNA and the semi-conservative replication of the parental DNA sequence into two new DNA strands, it became apparent that DNA replication is organized in a temporal and spatial fashion during the S phase of the cell cycle, correlated with the large-scale organization of chromatin in the nucleus. After many decades of limited progress, technological advances in genomics, genome engineering, and imaging have finally positioned the field to tackle mechanisms underpinning the temporal and spatial regulation of DNA replication and the causal relationships between DNA replication and other features of large-scale chromosome structure and function. In this review, we discuss these major recent discoveries as well as expectations for the coming decade.

Long noncoding RNA IL6-AS1 is highly expressed in chronic obstructive pulmonary disease and is associated with interleukin 6 by targeting miR-149-5p and early B-cell factor 1

Chronic obstructive pulmonary disease is a complex condition with multiple etiologies, including inflammation. We identified a novel long noncoding RNA (lncRNA), interleukin 6 antisense RNA 1 (IL6-AS1), which is upregulated in this disease and is associated with airway inflammation. We found that IL6-AS1 promotes the expression of inflammatory factors, especially interleukin (IL) 6. Mechanistically, cytoplasmic IL6-AS1 acts as an endogenous sponge by competitively binding to the microRNA miR-149-5p to stabilize IL-6 mRNA. Nuclear IL6-AS1 promotes IL-6 transcription by recruiting early B-cell factor 1 to the IL-6 promoter, which increases the methylation of the H3K4 histone and acetylation of the H3K27 histone. We propose a model of lncRNA expression in both the nucleus and cytoplasm that exerts similar effects through differing mechanisms, and IL6-AS1 probably increases inflammation via multiple pathways.

Tetranucleosome Interactions Drive Chromatin Folding

The multiscale organizational structure of chromatin in eukaryotic cells is instrumental to DNA transcription, replication, and repair. At mesoscopic length scales, nucleosomes pack in a manner that serves to regulate gene expression through condensation and expansion of the genome. The particular structures that arise and their respective thermodynamic stabilities, however, have yet to be fully resolved. In this study, we combine molecular modeling using the 1CPN mesoscale model of chromatin with nonlinear manifold learning to identify and characterize the structure and free energy of metastable states of short chromatin segments comprising between 4- and 16-nucleosomes. Our results reveal the formation of two previously characterized tetranucleosomal conformations, the "α-tetrahedron" and the "β-rhombus", which have been suggested to play an important role in the accessibility of DNA and, respectively, induce local chromatin compaction or elongation. The spontaneous formation of these motifs is potentially responsible for the slow nucleosome dynamics observed in experimental studies. Increases of the nucleosome repeat length are accompanied by more pronounced structural irregularity and flexibility and, ultimately, a dynamic liquid-like behavior that allows for frequent structural reorganization. Our findings indicate that tetranucleosome motifs are intrinsically stable structural states, driven by local internucleosomal interactions, and support a mechanistic picture of chromatin packing, dynamics, and accessibility that is strongly influenced by emergent local mesoscale structure.

SBTD: A Novel Method for Detecting Topological Associated Domains from Hi-C Data

The development of Hi-C technology has generated terabytes of chromatin interaction data, which bring possibilities for insight study of chromatin structure. Several studies revealed that mammalian chromosomes are folded into topological associated domains (TADs), which are conserved across cell types. Accurate detection of topological associated domains is now a vital process for revealing the relationship between the structure and function of genome organization. Unfortunately, the current TAD detection methods require massive computing resources, careful parameter adjustment and/or encounter inconsistent results. In this paper, we propose a novel method, Spectral-Based TAD Detector (SBTD), and evaluate its performance with a set of widely accepted statistical methods. We treat the chromatin interaction matrix as a graph and first introduce cosine similarity as a measure of the interaction patterns between bins. The results show that SBTD identifies higher quality TADs than the popular methods (DomainCaller, TopDom and SpectralTAD) and the internal bins of TADs identified by SBTD have higher correlation. Besides, The TADs identified by SBTD show a highly similar histone modification signal enrichment pattern at the boundary as reported in the previous literature. Finally, the motif enrichment analysis shows that compared with the background region, the DNA motifs of known insulator proteins are significantly enriched in the TAD boundary region identified by our method, which proves the high performance of our proposed method. Overall, SBTD is much more effective than existing methods with only one easy-to-adjust parameter, cluster number, for which we provide optimization guidelines.

G-quadruplex binding protein Rif1, a key regulator of replication timing

DNA replication is spatially and temporally regulated during S phase to execute efficient and coordinated duplication of entire genome. Various epigenomic mechanisms operate to regulate the timing and locations of replication. Among them, Rif1 plays a major role to shape the 'replication domains' that dictate which segments of the genome are replicated when and where in the nuclei. Rif1 achieves this task by generating higher-order chromatin architecture near nuclear membrane and by recruiting a protein phosphatase. Rif1 is a G4 binding protein, and G4 binding activity of Rif1 is essential for replication timing regulation in fission yeast. In this article, we first summarize strategies by which cells regulate their replication timing and then describe how Rif1 and its interaction with G4 contribute to regulation of chromatin architecture and replication timing.

Transcription shapes DNA replication initiation to preserve genome integrity

BACKGROUND

Early DNA replication occurs within actively transcribed chromatin compartments in mammalian cells, raising the immediate question of how early DNA replication coordinates with transcription to avoid collisions and DNA damage.

RESULTS

We develop a high-throughput nucleoside analog incorporation sequencing assay and identify thousands of early replication initiation zones in both mouse and human cells. The identified early replication initiation zones fall in open chromatin compartments and are mutually exclusive with transcription elongation. Of note, early replication initiation zones are mainly located in non-transcribed regions adjacent to transcribed regions. Mechanistically, we find that RNA polymerase II actively redistributes the chromatin-bound mini-chromosome maintenance complex (MCM), but not the origin recognition complex (ORC), to actively restrict early DNA replication initiation outside of transcribed regions. In support of this finding, we detect apparent MCM accumulation and DNA replication initiation in transcribed regions due to anchoring of nuclease-dead Cas9 at transcribed genes, which stalls RNA polymerase II. Finally, we find that the orchestration of early DNA replication initiation by transcription efficiently prevents gross DNA damage.

CONCLUSION

RNA polymerase II redistributes MCM complexes, but not the ORC, to prevent early DNA replication from initiating within transcribed regions. This RNA polymerase II-driven MCM redistribution spatially separates transcription and early DNA replication events and avoids the transcription-replication initiation collision, thereby providing a critical regulatory mechanism to preserve genome stability.

Integrative Analysis of Genome, 3D Genome, and Transcriptome Alterations of Clinical Lung Cancer Samples

Genomic studies of cancer cell alterations, such as mutations, copy number variations (CNVs), and translocations, greatly promote our understanding of the genesis and development of cancer. However, the 3D genome architecture of cancers remains less studied due to the complexity of cancer genomes and technical difficulties. To explore the 3D genome structure in clinical lung cancer, we performed Hi-C experiments using paired normal and tumor cells harvested from patients with lung cancer, combining with RNA-seq analysis. We demonstrated the feasibility of studying 3D genome of clinical lung cancer samples with a small number of cells (1 × 10⁴), compared the genome architecture between clinical samples and cell lines of lung cancer, and identified conserved and changed spatial chromatin structures between normal and cancer samples. We also showed that Hi-C data can be used to infer CNVs and point mutations in cancer. By integrating those different types of cancer alterations, we showed significant associations between CNVs, 3D genome, and gene expression. We propose that 3D genome mediates the effects of cancer genomic alterations on gene expression through altering regulatory chromatin structures. Our study highlights the importance of analyzing 3D genomes of clinical cancer samples in addition to cancer cell lines and provides an integrative genomic analysis pipeline for future larger-scale studies in lung cancer and other cancers.

How DNA loop extrusion mediated by cohesin enables V(D)J recombination

'Structural maintenance of chromosomes' (SMC) complexes are required for the folding of genomic DNA into loops. Theoretical considerations and single-molecule experiments performed with the SMC complexes cohesin and condensin indicate that DNA folding occurs via loop extrusion. Recent work indicates that this process is essential for the assembly of antigen receptor genes by V(D)J recombination in developing B and T cells of the vertebrate immune system. Here, I review how recent studies of the mouse immunoglobulin heavy chain locus Igh have provided evidence for this hypothesis and how the formation of chromatin loops by cohesin and regulation of this process by CTCF and Wapl might ensure that all variable gene segments in this locus (VH segments) participate in recombination with a re-arranged DJH segment, to ensure generation of a maximally diverse repertoire of B-cell receptors and antibodies.

Igh Locus Polymorphism May Dictate Topological Chromatin Conformation and V Gene Usage in the Ig Repertoire

Vast repertoires of unique antigen receptors are created in developing B and T lymphocytes. The antigen receptor loci contain many variable (V), diversity (D) and joining (J) gene segments that are arrayed across very large genomic expanses and are joined to form variable-region exons of expressed immunoglobulins and T cell receptors. This process creates the potential for an organism to respond to large numbers of different pathogens. Here, we consider the possibility that genetic polymorphisms with alterations in a vast array of regulatory elements in the immunoglobulin heavy chain (IgH) locus lead to changes in locus topology and impact immune-repertoire formation.

Nuclear organisation and replication timing are coupled through RIF1-PP1 interaction

Three-dimensional genome organisation and replication timing are known to be correlated, however, it remains unknown whether nuclear architecture overall plays an instructive role in the replication-timing programme and, if so, how. Here we demonstrate that RIF1 is a molecular hub that co-regulates both processes. Both nuclear organisation and replication timing depend upon the interaction between RIF1 and PP1. However, whereas nuclear architecture requires the full complement of RIF1 and its interaction with PP1, replication timing is not sensitive to RIF1 dosage. The role of RIF1 in replication timing also extends beyond its interaction with PP1. Availing of this separation-of-function approach, we have therefore identified in RIF1 dual function the molecular bases of the co-dependency of the replication-timing programme and nuclear architecture.

Hi-C as a molecular rangefinder to examine genomic rearrangements

The mammalian genome is highly packed into the nucleus. Over the past decade, the development of Hi-C has contributed significantly to our understanding of the three-dimensional (3D) chromatin structure, uncovering the principles and functions of higher-order chromatin organizations. Recent studies have repositioned its property in spatial proximity measurement to address challenging problems in genome analyses including genome assembly, haplotype phasing, and the detection of genomic rearrangements. In particular, the power of Hi-C in detecting large-scale structural variations (SVs) in the cancer genome has been demonstrated, which is challenging to be addressed solely with short-read-based whole-genome sequencing analyses. In this review, we first provide a comprehensive view of Hi-C as an intuitive and effective SV detection tool. Then, we introduce recently developed bioinformatics tools utilizing Hi-C to investigate genomic rearrangements. Finally, we discuss the potential application of single-cell Hi-C to address the heterogeneity of genomic rearrangements and sub-population identification in the cancer genome.

New gene discoveries in skeletal diseases with short stature

In the last decade, the widespread use of massively parallel sequencing has considerably boosted the number of novel gene discoveries in monogenic skeletal diseases with short stature. Defects in genes playing a role in the maintenance and function of the growth plate, the site of longitudinal bone growth, are a well-known cause of skeletal diseases with short stature. However, several genes involved in extracellular matrix composition or maintenance as well as genes partaking in various biological processes have also been characterized. This review aims to describe the latest genetic findings in spondyloepiphyseal dysplasias, spondyloepimetaphyseal dysplasias, and some monogenic forms of isolated short stature. Some examples of novel genetic mechanisms leading to skeletal conditions with short stature will be described. Strategies on how to successfully characterize novel skeletal phenotypes with short stature and genetic approaches to detect and validate novel gene-disease correlations will be discussed in detail. In summary, we review the latest gene discoveries underlying skeletal diseases with short stature and emphasize the importance of characterizing novel molecular mechanisms for genetic counseling, for an optimal management of the disease, and for therapeutic innovations.

Chromatin tracing and multiplexed imaging of nucleome architectures (MINA) and RNAs in single mammalian cells and tissue

The genome is hierarchically organized into several 3D architectures, including chromatin loops, domains, compartments and regions associated with nuclear lamina and nucleoli. Changes in these architectures have been associated with normal development, aging and a wide range of diseases. Despite its critical importance, understanding how the genome is spatially organized in single cells, how organization varies in different cell types in mammalian tissue and how organization affects gene expression remains a major challenge. Previous approaches have been limited by a lack of capacity to directly trace chromatin folding in 3D and to simultaneously measure genomic organization in relation to other nuclear components and gene expression in the same single cells. We have developed an image-based 3D genomics technique termed 'chromatin tracing', which enables direct 3D tracing of chromatin folding along individual chromosomes in single cells. More recently, we also developed multiplexed imaging of nucleome architectures (MINA), which enables simultaneous measurements of multiscale chromatin folding, associations of genomic regions with nuclear lamina and nucleoli and copy numbers of numerous RNA species in the same single cells in mammalian tissue. Here, we provide detailed protocols for chromatin tracing in cell lines and MINA in mammalian tissue, which take 3-4 d for experimental work and 2-3 d for data analysis. We expect these developments to be broadly applicable and to affect many lines of research on 3D genomics by depicting multiscale genomic architectures associated with gene expression, in different types of cells and tissue undergoing different biological processes.

Potential roles of condensin in genome organization and beyond in fission yeast

The genome is highly organized hierarchically by the function of structural maintenance of chromosomes (SMC) complex proteins such as condensin and cohesin from bacteria to humans. Although the roles of SMC complex proteins have been well characterized, their specialized roles in nuclear processes remain unclear. Condensin and cohesin have distinct binding sites and mediate long-range and short-range genomic associations, respectively, to form cell cycle-specific genome organization. Condensin can be recruited to highly expressed genes as well as dispersed repeat genetic elements, such as Pol III-transcribed genes, LTR retrotransposon, and rDNA repeat. In particular, mitotic transcription factors Ace2 and Ams2 recruit condensin to their target genes, forming centromeric clustering during mitosis. Condensin is potentially involved in various chromosomal processes such as the mobility of chromosomes, chromosome territories, DNA reannealing, and transcription factories. The current knowledge of condensin in fission yeast summarized in this review can help us understand how condensin mediates genome organization and participates in chromosomal processes in other organisms.

Replication Stress, Genomic Instability, and Replication Timing: A Complex Relationship

The replication-timing program constitutes a key element of the organization and coordination of numerous nuclear processes in eukaryotes. This program is established at a crucial moment in the cell cycle and occurs simultaneously with the organization of the genome, thus indicating the vital significance of this process. With recent technological achievements of high-throughput approaches, a very strong link has been confirmed between replication timing, transcriptional activity, the epigenetic and mutational landscape, and the 3D organization of the genome. There is also a clear relationship between replication stress, replication timing, and genomic instability, but the extent to which they are mutually linked to each other is unclear. Recent evidence has shown that replication timing is affected in cancer cells, although the cause and consequence of this effect remain unknown. However, in-depth studies remain to be performed to characterize the molecular mechanisms of replication-timing regulation and clearly identify different cis- and trans-acting factors. The results of these studies will potentially facilitate the discovery of new therapeutic pathways, particularly for personalized medicine, or new biomarkers. This review focuses on the complex relationship between replication timing, replication stress, and genomic instability.

Single-cell chromatin accessibility landscape identifies tissue repair program in human regulatory T cells

Murine regulatory T (Treg) cells in tissues promote tissue homeostasis and regeneration. We sought to identify features that characterize human Treg cells with these functions in healthy tissues. Single-cell chromatin accessibility profiles of murine and human tissue Treg cells defined a conserved, microbiota-independent tissue-repair Treg signature with a prevailing footprint of the transcription factor BATF. This signature, combined with gene expression profiling and TCR fate mapping, identified a population of tissue-like Treg cells in human peripheral blood that expressed BATF, chemokine receptor CCR8 and HLA-DR. Human BATF⁺CCR8⁺ Treg cells from normal skin and adipose tissue shared features with nonlymphoid T follicular helper-like (Tfh-like) cells, and induction of a Tfh-like differentiation program in naive human Treg cells partially recapitulated tissue Treg regenerative characteristics, including wound healing potential. Human BATF⁺CCR8⁺ Treg cells from healthy tissue share features with tumor-resident Treg cells, highlighting the importance of understanding the context-specific functions of these cells.

Utilization of Host Cell Chromosome Conformation by Viral Pathogens: Knowing When to Hold and When to Fold

Though viruses have their own genomes, many depend on the nuclear environment of their hosts for replication and survival. A substantial body of work has therefore been devoted to understanding how viral and eukaryotic genomes interact. Recent advances in chromosome conformation capture technologies have provided unprecedented opportunities to visualize how mammalian genomes are organized and, by extension, how packaging of nuclear DNA impacts cellular processes. Recent studies have indicated that some viruses, upon entry into host cell nuclei, produce factors that alter host chromatin topology, and thus, impact the 3D organization of the host genome. Additionally, a variety of distinct viruses utilize host genome architectural factors to advance various aspects of their life cycles. Indeed, human gammaherpesviruses, known for establishing long-term reservoirs of latent infection in B lymphocytes, utilize 3D principles of genome folding to package their DNA and establish latency in host cells. This manipulation of host epigenetic machinery by latent viral genomes is etiologically linked to the onset of B cell oncogenesis. Small DNA viruses, by contrast, are tethered to distinct cellular sites that support virus expression and replication. Here, we briefly review the recent findings on how viruses and host genomes spatially communicate, and how this impacts virus-induced pathology.

A unified framework for inferring the multi-scale organization of chromatin domains from Hi-C

Chromosomes are giant chain molecules organized into an ensemble of three-dimensional structures characterized with its genomic state and the corresponding biological functions. Despite the strong cell-to-cell heterogeneity, the cell-type specific pattern demonstrated in high-throughput chromosome conformation capture (Hi-C) data hints at a valuable link between structure and function, which makes inference of chromatin domains (CDs) from the pattern of Hi-C a central problem in genome research. Here we present a unified method for analyzing Hi-C data to determine spatial organization of CDs over multiple genomic scales. By applying statistical physics-based clustering analysis to a polymer physics model of the chromosome, our method identifies the CDs that best represent the global pattern of correlation manifested in Hi-C. The multi-scale intra-chromosomal structures compared across different cell types uncover the principles underlying the multi-scale organization of chromatin chain: (i) Sub-TADs, TADs, and meta-TADs constitute a robust hierarchical structure. (ii) The assemblies of compartments and TAD-based domains are governed by different organizational principles. (iii) Sub-TADs are the common building blocks of chromosome architecture. Our physically principled interpretation and analysis of Hi-C not only offer an accurate and quantitative view of multi-scale chromatin organization but also help decipher its connections with genome function.

An array of square blocks and checkerboard patterns characteristic to Hi-C data reflects the multi-scale organization of the chromatin chain. Deciphering the structures of chromatin domains from Hi-C and associating them with genome function are open problems of great importance in genome research. However, most existing methods are specialized in finding domains at different scales, making it difficult to integrate the solutions. Here we develop a unified framework for modeling and inferring domain structures over multiple scales, based on a physical model of the chromosome that reflects its nature as a three-dimensional object. Our method efficiently explores the space of domain solutions at different genomic scales, and systematically infers the chromatin domains over multiple scales from Hi-C data by employing a single tuning parameter. Our principled interpretation of Hi-C not only offers a quantitative view of multi-scale chromatin organization but also helps understand its connections with genome function.

Transcriptional Regulation of ROS Homeostasis by the ERR Subfamily of Nuclear Receptors

Reactive oxygen species (ROS) such as superoxide anion (O₂^•-) and hydrogen peroxide (H₂O₂) are generated endogenously by processes such as mitochondrial oxidative phosphorylation, or they may arise from exogenous sources like bacterial invasion. ROS can be beneficial (oxidative eustress) as signaling molecules but also harmful (oxidative distress) to cells when ROS levels become unregulated in response to physiological, pathological or pharmacological insults. Indeed, abnormal ROS levels have been shown to contribute to the etiology of a wide variety of diseases. Transcriptional control of metabolic genes is a crucial mechanism to coordinate ROS homeostasis. Therefore, a better understanding of how ROS metabolism is regulated by specific transcription factors can contribute to uncovering new therapeutic strategies. A large body of work has positioned the estrogen-related receptors (ERRs), transcription factors belonging to the nuclear receptor superfamily, as not only master regulators of cellular energy metabolism but, most recently, of ROS metabolism. Herein, we will review the role played by the ERRs as transcriptional regulators of ROS generation and antioxidant mechanisms and also as ROS sensors. We will assess how the control of ROS homeostasis by the ERRs can be linked to physiology and disease and the possible contribution of manipulating ERR activity in redox medicine.

Human ORC/MCM density is low in active genes and correlates with replication time but does not delimit initiation zones

Eukaryotic DNA replication initiates during S phase from origins that have been licensed in the preceding G1 phase. Here, we compare ChIP-seq profiles of the licensing factors Orc2, Orc3, Mcm3, and Mcm7 with gene expression, replication timing, and fork directionality profiles obtained by RNA-seq, Repli-seq, and OK-seq. Both, the origin recognition complex (ORC) and the minichromosome maintenance complex (MCM) are significantly and homogeneously depleted from transcribed genes, enriched at gene promoters, and more abundant in early- than in late-replicating domains. Surprisingly, after controlling these variables, no difference in ORC/MCM density is detected between initiation zones, termination zones, unidirectionally replicating regions, and randomly replicating regions. Therefore, ORC/MCM density correlates with replication timing but does not solely regulate the probability of replication initiation. Interestingly, H4K20me3, a histone modification proposed to facilitate late origin licensing, was enriched in late-replicating initiation zones and gene deserts of stochastic replication fork direction. We discuss potential mechanisms specifying when and where replication initiates in human cells.

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.

Control of ribosomal protein synthesis by the Microprocessor complex

Ribosome biogenesis in eukaryotes requires the coordinated production and assembly of 80 ribosomal proteins and four ribosomal RNAs (rRNAs), and its rate must be synchronized with cellular growth. Here, we showed that the Microprocessor complex, which mediates the first step of microRNA processing, potentiated the transcription of ribosomal protein genes by eliminating DNA/RNA hybrids known as R-loops. Nutrient deprivation triggered the nuclear export of Drosha, a key component of the Microprocessor complex, and its subsequent degradation by the E3 ubiquitin ligase Nedd4, thereby reducing ribosomal protein production and protein synthesis. In mouse erythroid progenitors, conditional deletion of Drosha led to the reduced production of ribosomal proteins, translational inhibition of the mRNA encoding the erythroid transcription factor Gata1, and impaired erythropoiesis. This phenotype mirrored the clinical presentation of human "ribosomopathies." Thus, the Microprocessor complex plays a pivotal role in synchronizing protein synthesis capacity with cellular growth rate and is a potential drug target for anemias caused by ribosomal insufficiency.

Involvement of TFAP2A in the activation of GSDMD gene promoter in hyperoxia-induced ALI

Oxygen therapy is a common treatment in neonatal intensive care units, but long-term continuous hyperoxia ventilation may induce acute lung injury (ALI). Gasdermin D (GSDMD)-mediated pyroptosis participates in various diseases including ALI, but the role of GSDMD in hyperoxia-induced ALI is yet understood. Here, we showed a significant increase in GSDMD after exposure to high oxygen. To elucidate the molecular mechanisms involved in GSDMD regulation, we identified the core promoter of GSDMD, -98 ~ -12 bp relative to the transcriptional start site (TSS). The results of mutational analysis, overexpression or siRNA interference, EMSA and ChIP demonstrated that E2F4 and TFAP2A positively regulate the transcriptional activity of the GSDMD by binding to its promoter. However, only TFAP2A showed a regulatory effect on the expression of GSDMD. Moreover, TFAP2A was increased in the lung tissues of rats exposed to hyperoxia and showed a strong linear correlation with GSDMD. Our results indicated that TFAP2A positively regulates the GSDMD expression via binding to the promoter region of GSDMD.

Chromosomal coordination and differential structure of asynchronous replicating regions

Stochastic asynchronous replication timing (AS-RT) is a phenomenon in which the time of replication of each allele is different, and the identity of the early allele varies between cells. By taking advantage of stable clonal pre-B cell populations derived from C57BL6/Castaneous mice, we have mapped the genome-wide AS-RT loci, independently of genetic differences. These regions are characterized by differential chromatin accessibility, mono-allelic expression and include new gene families involved in specifying cell identity. By combining population level mapping with single cell FISH, our data reveal the existence of a novel regulatory program that coordinates a fixed relationship between AS-RT regions on any given chromosome, with some loci set to replicate in a parallel and others set in the anti-parallel orientation. Our results show that AS-RT is a highly regulated epigenetic mark established during early embryogenesis that may be used for facilitating the programming of mono-allelic choice throughout development.

Most regions of the mammalian genome replicate both alleles in a synchronous manner, but some loci have been found to replicate asynchronously and the time of replication of each allele is different. Here the authors, by employing clonal mouse cells from a hybrid strain chart replication timing over the entire genome, using polymorphisms to distinguish between the paternal and maternal alleles.

CTCF-binding element regulates ESC differentiation via orchestrating long-range chromatin interaction between enhancers and HoxA

Proper expression of Homeobox A cluster genes (HoxA) is essential for embryonic stem cell (ESC) differentiation and individual development. However, mechanisms controlling precise spatiotemporal expression of HoxA during early ESC differentiation remain poorly understood. Herein, we identified a functional CTCF-binding element (CBE⁺⁴⁷) closest to the 3'-end of HoxA within the same topologically associated domain (TAD) in ESC. CRISPR-Cas9-mediated deletion of CBE⁺⁴⁷ significantly upregulated HoxA expression and enhanced early ESC differentiation induced by retinoic acid (RA) relative to wild-type cells. Mechanistic analysis by chromosome conformation capture assay (Capture-C) revealed that CBE⁺⁴⁷ deletion decreased interactions between adjacent enhancers, enabling formation of a relatively loose enhancer-enhancer interaction complex (EEIC), which overall increased interactions between that EEIC and central regions of HoxA chromatin. These findings indicate that CBE⁺⁴⁷ organizes chromatin interactions between its adjacent enhancers and HoxA. Furthermore, deletion of those adjacent enhancers synergistically inhibited HoxA activation, suggesting that these enhancers serve as an EEIC required for RA-induced HoxA activation. Collectively, these results provide new insight into RA-induced HoxA expression during early ESC differentiation, also highlight precise regulatory roles of the CTCF-binding element in orchestrating high-order chromatin structure.

Topologically associating domains and their role in the evolution of genome structure and function in Drosophila

Topologically associating domains (TADs) were recently identified as fundamental units of three-dimensional eukaryotic genomic organization, although our knowledge of the influence of TADs on genome evolution remains preliminary. To study the molecular evolution of TADs in Drosophila species, we constructed a new reference-grade genome assembly and accompanying high-resolution TAD map for D. pseudoobscura Comparison of D. pseudoobscura and D. melanogaster, which are separated by ∼49 million years of divergence, showed that ∼30%-40% of their genomes retain conserved TADs. Comparative genomic analysis of 17 Drosophila species revealed that chromosomal rearrangement breakpoints are enriched at TAD boundaries but depleted within TADs. Additionally, genes within conserved TADs show lower expression divergence than those located in nonconserved TADs. Furthermore, we found that a substantial proportion of long genes (>50 kbp) in D. melanogaster (42%) and D. pseudoobscura (26%) constitute their own TADs, implying transcript structure may be one of the deterministic factors for TAD formation. By using structural variants (SVs) identified from 14 D. melanogaster strains, its three closest sibling species from the D. simulans species complex, and two obscura clade species, we uncovered evidence of selection acting on SVs at TAD boundaries, but with the nature of selection differing between SV types. Deletions are depleted at TAD boundaries in both divergent and polymorphic SVs, suggesting purifying selection, whereas divergent tandem duplications are enriched at TAD boundaries relative to polymorphism, suggesting they are adaptive. Our findings highlight how important TADs are in shaping the acquisition and retention of structural mutations that fundamentally alter genome organization.

Human DDK rescues stalled forks and counteracts checkpoint inhibition at unfired origins to complete DNA replication

Eukaryotic genomes replicate via spatially and temporally regulated origin firing. Cyclin-dependent kinase (CDK) and Dbf4-dependent kinase (DDK) promote origin firing, whereas the S phase checkpoint limits firing to prevent nucleotide and RPA exhaustion. We used chemical genetics to interrogate human DDK with maximum precision, dissect its relationship with the S phase checkpoint, and identify DDK substrates. We show that DDK inhibition (DDKi) leads to graded suppression of origin firing and fork arrest. S phase checkpoint inhibition rescued origin firing in DDKi cells and DDK-depleted Xenopus egg extracts. DDKi also impairs RPA loading, nascent-strand protection, and fork restart. Via quantitative phosphoproteomics, we identify the BRCA1-associated (BRCA1-A) complex subunit MERIT40 and the cohesin accessory subunit PDS5B as DDK effectors in fork protection and restart. Phosphorylation neutralizes autoinhibition mediated by intrinsically disordered regions in both substrates. Our results reveal mechanisms through which DDK controls the duplication of large vertebrate genomes.

SuperTAD: robust detection of hierarchical topologically associated domains with optimized structural information

Topologically associating domains (TADs) are the organizational units of chromosome structures. TADs can contain TADs, thus forming a hierarchy. TAD hierarchies can be inferred from Hi-C data through coding trees. However, the current method for computing coding trees is not optimal. In this paper, we propose optimal algorithms for this computation. In comparison with seven state-of-art methods using two public datasets, from GM12878 and IMR90 cells, SuperTAD shows a significant enrichment of structural proteins around detected boundaries and histone modifications within TADs and displays a high consistency between various resolutions of identical Hi-C matrices.

Large organized chromatin lysine domains help distinguish primitive from differentiated cell populations

The human genome is partitioned into a collection of genomic features, inclusive of genes, transposable elements, lamina interacting regions, early replicating control elements and cis-regulatory elements, such as promoters, enhancers, and anchors of chromatin interactions. Uneven distribution of these features within chromosomes gives rise to clusters, such as topologically associating domains (TADs), lamina-associated domains, clusters of cis-regulatory elements or large organized chromatin lysine (K) domains (LOCKs). Here we show that LOCKs from diverse histone modifications discriminate primitive from differentiated cell types. Active LOCKs (H3K4me1, H3K4me3 and H3K27ac) cover a higher fraction of the genome in primitive compared to differentiated cell types while repressive LOCKs (H3K9me3, H3K27me3 and H3K36me3) do not. Active LOCKs in differentiated cells lie proximal to highly expressed genes while active LOCKs in primitive cells tend to be bivalent. Genes proximal to bivalent LOCKs are minimally expressed in primitive cells. Furthermore, bivalent LOCKs populate TAD boundaries and are preferentially bound by regulators of chromatin interactions, including CTCF, RAD21 and ZNF143. Together, our results argue that LOCKs discriminate primitive from differentiated cell populations.

3DIV update for 2021: a comprehensive resource of 3D genome and 3D cancer genome

Three-dimensional (3D) genome organization is tightly coupled with gene regulation in various biological processes and diseases. In cancer, various types of large-scale genomic rearrangements can disrupt the 3D genome, leading to oncogenic gene expression. However, unraveling the pathogenicity of the 3D cancer genome remains a challenge since closer examinations have been greatly limited due to the lack of appropriate tools specialized for disorganized higher-order chromatin structure. Here, we updated a 3D-genome Interaction Viewer and database named 3DIV by uniformly processing ∼230 billion raw Hi-C reads to expand our contents to the 3D cancer genome. The updates of 3DIV are listed as follows: (i) the collection of 401 samples including 220 cancer cell line/tumor Hi-C data, 153 normal cell line/tissue Hi-C data, and 28 promoter capture Hi-C data, (ii) the live interactive manipulation of the 3D cancer genome to simulate the impact of structural variations and (iii) the reconstruction of Hi-C contact maps by user-defined chromosome order to investigate the 3D genome of the complex genomic rearrangement. In summary, the updated 3DIV will be the most comprehensive resource to explore the gene regulatory effects of both the normal and cancer 3D genome. ‘3DIV’ is freely available at http://3div.kr.

High-resolution single-cell 3D-models of chromatin ensembles during Drosophila embryogenesis

Single-cell chromatin studies provide insights into how chromatin structure relates to functions of individual cells. However, balancing high-resolution and genome wide-coverage remains challenging. We describe a computational method for the reconstruction of large 3D-ensembles of single-cell (sc) chromatin conformations from population Hi-C that we apply to study embryogenesis in Drosophila. With minimal assumptions of physical properties and without adjustable parameters, our method generates large ensembles of chromatin conformations via deep-sampling. Our method identifies specific interactions, which constitute 5-6% of Hi-C frequencies, but surprisingly are sufficient to drive chromatin folding, giving rise to the observed Hi-C patterns. Modeled sc-chromatins quantify chromatin heterogeneity, revealing significant changes during embryogenesis. Furthermore, >50% of modeled sc-chromatin maintain topologically associating domains (TADs) in early embryos, when no population TADs are perceptible. Domain boundaries become fixated during development, with strong preference at binding-sites of insulator-complexes upon the midblastula transition. Overall, high-resolution 3D-ensembles of sc-chromatin conformations enable further in-depth interpretation of population Hi-C, improving understanding of the structure-function relationship of genome organization.

A new duck genome reveals conserved and convergently evolved chromosome architectures of birds and mammals

BACKGROUND

Ducks have a typical avian karyotype that consists of macro- and microchromosomes, but a pair of much less differentiated ZW sex chromosomes compared to chickens. To elucidate the evolution of chromosome architectures between ducks and chickens, and between birds and mammals, we produced a nearly complete chromosomal assembly of a female Pekin duck by combining long-read sequencing and multiplatform scaffolding techniques.

RESULTS

A major improvement of genome assembly and annotation quality resulted from the successful resolution of lineage-specific propagated repeats that fragmented the previous Illumina-based assembly. We found that the duck topologically associated domains (TAD) are demarcated by putative binding sites of the insulator protein CTCF, housekeeping genes, or transitions of active/inactive chromatin compartments, indicating conserved mechanisms of spatial chromosome folding with mammals. There are extensive overlaps of TAD boundaries between duck and chicken, and also between the TAD boundaries and chromosome inversion breakpoints. This suggests strong natural selection pressure on maintaining regulatory domain integrity, or vulnerability of TAD boundaries to DNA double-strand breaks. The duck W chromosome retains 2.5-fold more genes relative to chicken. Similar to the independently evolved human Y chromosome, the duck W evolved massive dispersed palindromic structures, and a pattern of sequence divergence with the Z chromosome that reflects stepwise suppression of homologous recombination.

CONCLUSIONS

Our results provide novel insights into the conserved and convergently evolved chromosome features of birds and mammals, and also importantly add to the genomic resources for poultry studies.

Unified Analysis of Multiple ChIP-Seq Datasets

High-throughput sequencing technologies are increasingly used in molecular cell biology to assess genome-wide chromatin dynamics of proteins bound to DNA, through techniques such as chromatin immunoprecipitation sequencing (ChIP-seq). These techniques often rely on an analysis strategy based on identifying genomic regions with increased sequencing signal to infer the binding location or chemical modifications of proteins bound to DNA. Peak calling within individual samples has been well described, however relatively little attention has been devoted to the merging of replicate samples, and the cross-comparison of many samples. Here, we present a generalized strategy to enable the unification of ChIP-seq datasets, enabling enhanced cross-comparison of binding patterns. The strategy works by merging peak data between different (even unrelated) samples, and then using a local background to recalculate enrichment. This strategy redefines the peaks within each experiment, allowing for more accurate cross-comparison of datasets.

Close to the edge: Heterochromatin at the nucleolar and nuclear peripheries

Chromatin is a dynamic structure composed of DNA, RNA, and proteins, regulating storage and expression of the genetic material in the nucleus. Heterochromatin plays a crucial role in driving the three-dimensional arrangement of the interphase genome, and in preserving genome stability by maintaining a subset of the genome in a silent state. Spatial genome organization contributes to normal patterns of gene function and expression, and is therefore of broad interest. Mammalian heterochromatin, the focus of this review, mainly localizes at the nuclear periphery, forming Lamina-associated domains (LADs), and at the nucleolar periphery, forming Nucleolus-associated domains (NADs). Together, these regions comprise approximately one-half of mammalian genomes, and most but not all loci within these domains are stochastically placed at either of these two locations after exit from mitosis at each cell cycle. Excitement about the role of these heterochromatic domains in early development has recently been heightened by the discovery that LADs appear at some loci in the preimplantation mouse embryo prior to other chromosomal features like compartmental identity and topologically-associated domains (TADs). While LADs have been extensively studied and mapped during cellular differentiation and early embryonic development, NADs have been less thoroughly studied. Here, we summarize pioneering studies of NADs and LADs, more recent advances in our understanding of cis/trans-acting factors that mediate these localizations, and discuss the functional significance of these associations.

Two Separate Cases: Complex Chromosomal Abnormality Involving Three Chromosomes and Small Supernumerary Marker Chromosome in Patients with Impaired Reproductive Function

For medical genetic counseling, estimating the chance of a child being born with chromosome abnormality is crucially important. Cytogenetic diagnostics of parents with a balanced karyotype are a special case. Such chromosome rearrangements cannot be detected with comprehensive chromosome screening. In the current paper, we consider chromosome diagnostics in two cases of chromosome rearrangement in patients with balanced karyotype and provide the results of a detailed analysis of complex chromosomal rearrangement (CCR) involving three chromosomes and a small supernumerary marker chromosome (sSMC) in a patient with impaired reproductive function. The application of fluorescent in situ hybridization, microdissection, and multicolor banding allows for describing analyzed karyotypes in detail. In the case of a CCR, such as the one described here, the probability of gamete formation with a karyotype, showing a balance of chromosome regions, is extremely low. Recommendation for the family in genetic counseling should take into account the obtained result. In the case of an sSMC, it is critically important to identify the original chromosome from which the sSMC has been derived, even if the euchromatin material is absent. Finally, we present our view on the optimal strategy of identifying and describing sSMCs, namely the production of a microdissectional DNA probe from the sSMC combined with a consequent reverse painting.