Hi-C metagenomics facilitate comparative genome analysis of bacteria and yeast from spontaneous beer and cider

Sequence-based analysis of fermented foods and beverages' microbiomes offers insights into their impact on taste and consumer health. High-throughput metagenomics provide detailed taxonomic and functional community profiling, but bacterial and yeast genome reconstruction and mobile genetic elements tracking are to be improved. We established a pipeline for exploring fermented foods microbiomes using metagenomics coupled with chromosome conformation capture (Hi-C metagenomics). The approach was applied to analyze a collection of spontaneously fermented beers and ciders (n = 12). The Hi-C reads were used to reconstruct the metagenome-assembled genomes (MAGs) of bacteria and yeasts facilitating subsequent comparative genomic analysis, assembly scaffolding and exploration of "plasmid-bacteria" links. For a subset of beverages, yeasts were isolated and characterized phenotypically. The reconstructed Hi-C MAGs primarily belonged to the Lactobacillaceae family in beers, along with Acetobacteraceae and Enterobacteriaceae in ciders, exhibiting improved quality compared to conventional metagenomic MAGs. Comparative genomic analysis of Lactobacillaceae Hi-C MAGs revealed clustering by niche and suggested genetic determinants of survival and probiotic potential. For Pediococcus damnosus, Hi-C-based networks of contigs enabled linking bacteria with plasmids. Analyzing phylogeny and accessory genes in the context of known reference genomes offered insights into the niche specialization of beer lactobacilli. The subspecies-level diversity of cider Tatumella spp. was disentangled using a Hi-C-based graph. We obtained highly complete yeast Hi-C MAGs primarily represented by Brettanomyces and Saccharomyces, with Hi-C-facilitated chromosome-level genome assembly for the former. Utilizing Hi-C metagenomics to unravel the genomic content of individual species can provide a deeper understanding of the ecological interactions within the food microbiome, aid in bioprospecting beneficial microorganisms, improving quality control and improving innovative fermented products.

Extensive long-range polycomb interactions and weak compartmentalization are hallmarks of human neuronal 3D genome

Chromatin architecture regulates gene expression and shapes cellular identity, particularly in neuronal cells. Specifically, polycomb group (PcG) proteins enable establishment and maintenance of neuronal cell type by reorganizing chromatin into repressive domains that limit the expression of fate-determining genes and sustain distinct gene expression patterns in neurons. Here, we map the 3D genome architecture in neuronal and non-neuronal cells isolated from the Wernicke's area of four human brains and comprehensively analyze neuron-specific aspects of chromatin organization. We find that genome segregation into active and inactive compartments is greatly reduced in neurons compared to other brain cells. Furthermore, neuronal Hi-C maps reveal strong long-range interactions, forming a specific network of PcG-mediated contacts in neurons that is nearly absent in other brain cells. These interacting loci contain developmental transcription factors with repressed expression in neurons and other mature brain cells. But only in neurons, they are rich in bivalent promoters occupied by H3K4me3 histone modification together with H3K27me3, which points to a possible functional role of PcG contacts in neurons. Importantly, other layers of chromatin organization also exhibit a distinct structure in neurons, characterized by an increase in short-range interactions and a decrease in long-range ones.

Trisomies Reorganize Human 3D Genome

Trisomy is the presence of one extra copy of an entire chromosome or its part in a cell nucleus. In humans, autosomal trisomies are associated with severe developmental abnormalities leading to embryonic lethality, miscarriage or pronounced deviations of various organs and systems at birth. Trisomies are characterized by alterations in gene expression level, not exclusively on the trisomic chromosome, but throughout the genome. Here, we applied the high-throughput chromosome conformation capture technique (Hi-C) to study chromatin 3D structure in human chorion cells carrying either additional chromosome 13 (Patau syndrome) or chromosome 16 and in cultured fibroblasts with extra chromosome 18 (Edwards syndrome). The presence of extra chromosomes results in systematic changes of contact frequencies between small and large chromosomes. Analyzing the behavior of individual chromosomes, we found that a limited number of chromosomes change their contact patterns stochastically in trisomic cells and that it could be associated with lamina-associated domains (LAD) and gene content. For trisomy 13 and 18, but not for trisomy 16, the proportion of compacted loci on a chromosome is correlated with LAD content. We also found that regions of the genome that become more compact in trisomic cells are enriched in housekeeping genes, indicating a possible decrease in chromatin accessibility and transcription level of these genes. These results provide a framework for understanding the mechanisms of pan-genome transcription dysregulation in trisomies in the context of chromatin spatial organization.

Unraveling the Polysaccharide Biosynthesis Potential of Ganoderma lucidum: A Chromosome-Level Assembly Using Hi-C Sequencing

Ganoderma lucidum exhibits the ability to synthesize a diverse range of biologically active molecules with significant pharmaceutical potential, including xylomannan and fucogalactan, which have demonstrated antitumor activity. However, there exists considerable intra-species variability in the capacity to produce these metabolites at high concentrations, likely reflecting the high genomic diversity observed from a limited number of strains sequenced to date. We employed high-throughput shotgun sequencing to obtain the complete genome sequence of G. lucidum strain 5.1, which is distinguished by its remarkable xylomannan synthesis capabilities. Through the utilization of semi-automatic reordering based on conformation capture (Hi-C) data, we substantially enhanced the assembly process, resulting in the generation of 12 chromosome-level scaffolds with a cumulative length of 39 Mbp. By employing both de novo and homology-based approaches, we performed comprehensive annotation of the genome, thereby identifying a diverse repertoire of genes likely involved in polysaccharide biosynthesis. The genome sequence generated in this study serves as a valuable resource for elucidating the molecular mechanisms underlying the medicinal potential of Ganoderma species, discovering novel pharmaceutically valuable compounds, and elucidating the ecological mechanisms of the species. Furthermore, the chromosome contact map obtained for the first time for this species extends our understanding of 3D fungal genomics and provides insights into the functional and structural organization within the fungal kingdom.

Extrusion fountains are hallmarks of chromosome organization emerging upon zygotic genome activation

The first activation of gene expression during development (zygotic genome activation, ZGA) is accompanied by massive changes in chromosome organization. The connection between these two processes remains unknown. Using Hi-C for zebrafish embryos, we found that chromosome folding starts by establishing "fountains", novel elements of chromosome organization, emerging selectively at enhancers upon ZGA. Using polymer simulations, we demonstrate that fountains can emerge as sites of targeted cohesin loading and require two-sided, yet desynchronized, loop extrusion. Specific loss of fountains upon loss of pioneer transcription factors that drive ZGA reveals a causal connection between enhancer activity and fountain formation. Finally, we show that fountains emerge in early Medaka and Xenopus embryos; moreover, we found cohesin-dependent fountain pattern on enhancers of mouse embryonic stem cells. Taken together, fountains are the first enhancer-specific elements of chromosome organization; they constitute starting points of chromosome folding during early development, likely serving as sites of targeted cohesin loading.

Enhancer Function in the 3D Genome

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

RNA-DNA interactomes of three prokaryotes uncovered by proximity ligation

Proximity ligation approaches, which are widely used to study the spatial organization of the genome, also make it possible to reveal patterns of RNA-DNA interactions. Here, we use RedC, an RNA-DNA proximity ligation approach, to assess the distribution of major RNA types along the genomes of E. coli, B. subtilis, and thermophilic archaeon T. adornatum. We find that (i) messenger RNAs preferentially interact with their cognate genes and the genes located downstream in the same operon, which is consistent with polycistronic transcription; (ii) ribosomal RNAs preferentially interact with active protein-coding genes in both bacteria and archaea, indicating co-transcriptional translation; and (iii) 6S noncoding RNA, a negative regulator of bacterial transcription, is depleted from active genes in E. coli and B. subtilis. We conclude that the RedC data provide a rich resource for studying both transcription dynamics and the function of noncoding RNAs in microbial organisms.

3D genome alterations and editing in pathology

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

HiConfidence: a novel approach uncovering the biological signal in Hi-C data affected by technical biases

The chromatin interaction assays, particularly Hi-C, enable detailed studies of genome architecture in multiple organisms and model systems, resulting in a deeper understanding of gene expression regulation mechanisms mediated by epigenetics. However, the analysis and interpretation of Hi-C data remain challenging due to technical biases, limiting direct comparisons of datasets obtained in different experiments and laboratories. As a result, removing biases from Hi-C-generated chromatin contact matrices is a critical data analysis step. Our novel approach, HiConfidence, eliminates biases from the Hi-C data by weighing chromatin contacts according to their consistency between replicates so that low-quality replicates do not substantially influence the result. The algorithm is effective for the analysis of global changes in chromatin structures such as compartments and topologically associating domains. We apply the HiConfidence approach to several Hi-C datasets with significant technical biases, that could not be analyzed effectively using existing methods, and obtain meaningful biological conclusions. In particular, HiConfidence aids in the study of how changes in histone acetylation pattern affect chromatin organization in Drosophila melanogaster S2 cells. The method is freely available at GitHub: https://github.com/victorykobets/HiConfidence.

Segregation of α- and β-Globin Gene Cluster in Vertebrate Evolution: Chance or Necessity?

The review is devoted to the patterns of evolution of α- and β-globin gene domains. A hypothesis is presented according to which segregation of the ancestral cluster of α/β-globin genes in Amniota occurred due to the performance by α-globins and β-globins of non-canonical functions not related to oxygen transport.

Genome-Directed Cell Nucleus Assembly

Simple Summary

Speckles and other nuclear bodies, the nucleolus and perinucleolar zone, transcription/replication factories and the lamina-associated compartment, serve as a structural basis for various genomic functions. In turn, genome activity and specific chromatin 3D organization directly impact the integrity of intranuclear assemblies, initiating/facilitating their formation and dictating their composition. Thus, the large-scale nucleus structure and genome activity mutually influence each other. The cell nucleus is frequently considered a compartment in which the genome is placed to protect it from external forces. Here, we discuss the evidence demonstrating that the cell nucleus should be considered, rather, as structure built around the folded genome. Decondensing chromosomes provide a scaffold for the assembly of the nuclear envelope after mitosis, whereas genome activity directs the assembly of various nuclear compartments, including nucleolus, speckles and transcription factories.

Abstract

The cell nucleus is frequently considered a cage in which the genome is placed to protect it from various external factors. Inside the nucleus, many functional compartments have been identified that are directly or indirectly involved in implementing genomic DNA’s genetic functions. For many years, it was assumed that these compartments are assembled on a proteinaceous scaffold (nuclear matrix), which provides a structural milieu for nuclear compartmentalization and genome folding while simultaneously offering some rigidity to the cell nucleus. The results of research in recent years have made it possible to consider the cell nucleus from a different angle. From the “box” in which the genome is placed, the nucleus has become a kind of mobile exoskeleton, which is formed around the packaged genome, under the influence of transcription and other processes directly related to the genome activity. In this review, we summarize the main arguments in favor of this point of view by analyzing the mechanisms that mediate cell nucleus assembly and support its resistance to mechanical stresses.

SETDB1 fuels the lung cancer phenotype by modulating epigenome, 3D genome organization and chromatin mechanical properties

Imbalance in the finely orchestrated system of chromatin-modifying enzymes is a hallmark of many pathologies such as cancers, since causing the affection of the epigenome and transcriptional reprogramming. Here, we demonstrate that a loss-of-function mutation (LOF) of the major histone lysine methyltransferase SETDB1 possessing oncogenic activity in lung cancer cells leads to broad changes in the overall architecture and mechanical properties of the nucleus through genome-wide redistribution of heterochromatin, which perturbs chromatin spatial compartmentalization. Together with the enforced activation of the epithelial expression program, cytoskeleton remodeling, reduced proliferation rate and restricted cellular migration, this leads to the reversed oncogenic potential of lung adenocarcinoma cells. These results emphasize an essential role of chromatin architecture in the determination of oncogenic programs and illustrate a relationship between gene expression, epigenome, 3D genome and nuclear mechanics.

Pentad: a tool for distance-dependent analysis of Hi-C interactions within and between chromatin compartments

Background

Understanding the role of various factors in 3D genome organization is essential to determine their impact on shaping large-scale chromatin units such as euchromatin (A) and heterochromatin (B) compartments. At this level, chromatin compaction is extensively modulated when transcription and epigenetic profiles change upon cell differentiation and response to various external impacts. However, detailed analysis of chromatin contact patterns within and between compartments is complicated because of a lack of suitable computational methods.

Results

We developed a tool, Pentad, to perform calculation, visualisation and quantitative analysis of the average chromatin compartment from the Hi-C matrices in cis, trans, and specified genomic distances. As we demonstrated by applying Pentad to publicly available Hi-C datasets, it helps to reliably detect redistribution of contact frequency in the chromatin compartments and assess alterations in the compartment strength.

Conclusions

Pentad is a simple tool for the analysis of changes in chromatin compartmentalization in various biological conditions. Pentad is freely available at https://github.com/magnitov/pentad.

Hi-C Metagenomics in the ICU: Exploring Clinically Relevant Features of Gut Microbiome in Chronically Critically Ill Patients

Gut microbiome in critically ill patients shows profound dysbiosis. The most vulnerable is the subgroup of chronically critically ill (CCI) patients – those suffering from long-term dependence on support systems in intensive care units. It is important to investigate their microbiome as a potential reservoir of opportunistic taxa causing co-infections and a morbidity factor. We explored dynamics of microbiome composition in the CCI patients by combining “shotgun” metagenomics with chromosome conformation capture (Hi-C). Stool samples were collected at 2 time points from 2 patients with severe brain injury with different outcomes within a 1–2-week interval. The metagenome-assembled genomes (MAGs) were reconstructed based on the Hi-C data using a novel hicSPAdes method (along with the bin3c method for comparison), as well as independently of the Hi-C using MetaBAT2. The resistomes of the samples were derived using a novel assembly graph-based approach. Links of bacteria to antibiotic resistance genes, plasmids and viruses were analyzed using Hi-C-based networks. The gut community structure was enriched in opportunistic microorganisms. The binning using hicSPAdes was superior to the conventional WGS-based binning as well as to the bin3c in terms of the number, completeness and contamination of the reconstructed MAGs. Using Klebsiella pneumoniae as an example, we showed how chromosome conformation capture can aid comparative genomic analysis of clinically important pathogens. Diverse associations of resistome with antimicrobial therapy from the level of assembly graphs to gene content were discovered. Analysis of Hi-C networks suggested multiple “host-plasmid” and “host-phage” links. Hi-C metagenomics is a promising technique for investigating clinical microbiome samples. It provides a community composition profile with increased details on bacterial gene content and mobile genetic elements compared to conventional metagenomics. The ability of Hi-C binning to encompass the MAG’s plasmid content facilitates metagenomic evaluation of virulence and drug resistance dynamics in clinically relevant opportunistic pathogens. These findings will help to identify the targets for developing cost-effective and rapid tests for assessing microbiome-related health risks.

Suppression of liquid-liquid phase separation by 1,6-hexanediol partially compromises the 3D genome organization in living cells

Liquid-liquid phase separation (LLPS) contributes to the spatial and functional segregation of molecular processes within the cell nucleus. However, the role played by LLPS in chromatin folding in living cells remains unclear. Here, using stochastic optical reconstruction microscopy (STORM) and Hi-C techniques, we studied the effects of 1,6-hexanediol (1,6-HD)-mediated LLPS disruption/modulation on higher-order chromatin organization in living cells. We found that 1,6-HD treatment caused the enlargement of nucleosome clutches and their more uniform distribution in the nuclear space. At a megabase-scale, chromatin underwent moderate but irreversible perturbations that resulted in the partial mixing of A and B compartments. The removal of 1,6-HD from the culture medium did not allow chromatin to acquire initial configurations, and resulted in more compact repressed chromatin than in untreated cells. 1,6-HD treatment also weakened enhancer-promoter interactions and TAD insulation but did not considerably affect CTCF-dependent loops. Our results suggest that 1,6-HD-sensitive LLPS plays a limited role in chromatin spatial organization by constraining its folding patterns and facilitating compartmentalization at different levels.

The two waves in single-cell 3D genomics

For decades, biochemical methods for the analysis of genome structure and function provided cell-population-averaged data that allowed general principles and tendencies to be disclosed. Microscopy-based studies, which immanently involve single-cell analysis, did not provide sufficient spatial resolution to investigate the particularly small details of 3D genome folding. Nevertheless, these studies demonstrated that mutual positions of chromosome territories within cell nuclei and individual genomic loci within chromosomal territories can vary significantly in individual cells. The development of new technologies in biochemistry and the advent of super-resolution microscopy in the last decade have made possible the full-scale study of 3D genome organization in individual cells. Maps of the 3D genome build based on C-data and super-resolution microscopy are highly consistent and, therefore, biologically relevant. The internal structures of individual chromosomes, loci, and topologically associating domains (TADs) are resolved as well as cell-cycle dynamics. 3D modeling allows one to investigate the physical mechanisms underlying genome folding. Finally, joint profiling of genome topology and epigenetic features will allow 3D genomics to handle complex cell-to-cell heterogeneity. In this review, we summarize the present state of studies into 3D genome organization in individual cells, analyze the technical problems of single-cell studies, and outline perspectives of 3D genomics.

LASCA: loop and significant contact annotation pipeline

Chromatin loops represent one of the major levels of hierarchical folding of the genome. Although the situation is evolving, current methods have various difficulties with the accurate mapping of loops even in mammalian Hi-C data, and most of them fail to identify chromatin loops in animal species with substantially different genome architecture. This paper presents the loop and significant contact annotation (LASCA) pipeline, which uses Weibull distribution-based modeling to effectively identify loops and enhancer–promoter interactions in Hi-C data from evolutionarily distant species: from yeast and worms to mammals. Available at: https://github.com/ArtemLuzhin/LASCA_pipeline.

Divide and Rule: Phase Separation in Eukaryotic Genome Functioning

The functioning of a cell at various organizational levels is determined by the interactions between macromolecules that promote cellular organelle formation and orchestrate metabolic pathways via the control of enzymatic activities. Although highly specific and relatively stable protein-protein, protein-DNA, and protein-RNA interactions are traditionally suggested as the drivers for cellular function realization, recent advances in the discovery of weak multivalent interactions have uncovered the role of so-called macromolecule condensates. These structures, which are highly divergent in size, composition, function, and cellular localization are predominantly formed by liquid-liquid phase separation (LLPS): a physical-chemical process where an initially homogenous solution turns into two distinct phases, one of which contains the major portion of the dissolved macromolecules and the other one containing the solvent. In a living cell, LLPS drives the formation of membrane-less organelles such as the nucleolus, nuclear bodies, and viral replication factories and facilitates the assembly of complex macromolecule aggregates possessing regulatory, structural, and enzymatic functions. Here, we discuss the role of LLPS in the spatial organization of eukaryotic chromatin and regulation of gene expression in normal and pathological conditions.

Benchmark of software tools for prokaryotic chromosomal interaction domain identification

Motivation

The application of genome-wide chromosome conformation capture (3C) methods to prokaryotes provided insights into the spatial organization of their genomes and identified patterns conserved across the tree of life, such as chromatin compartments and contact domains. Prokaryotic genomes vary in GC content and the density of restriction sites along the chromosome, suggesting that these properties should be considered when planning experiments and choosing appropriate software for data processing. Diverse algorithms are available for the analysis of eukaryotic chromatin contact maps, but their potential application to prokaryotic data has not yet been evaluated.

Results

Here, we present a comparative analysis of domain calling algorithms using available single-microbe experimental data. We evaluated the algorithms’ intra-dataset reproducibility, concordance with other tools and sensitivity to coverage and resolution of contact maps. Using RNA-seq as an example, we showed how orthogonal biological data can be utilized to validate the reliability and significance of annotated domains. We also suggest that in silico simulations of contact maps can be used to choose optimal restriction enzymes and estimate theoretical map resolutions before the experiment. Our results provide guidelines for researchers investigating microbes and microbial communities using high-throughput 3C assays such as Hi-C and 3C-seq.

Availability and implementation

The code of the analysis is available at https://github.com/magnitov/prokaryotic_cids.

Supplementary information

Supplementary data are available at Bioinformatics online.

Studying RNA-DNA interactome by Red-C identifies noncoding RNAs associated with various chromatin types and reveals transcription dynamics

Non-coding RNAs (ncRNAs) participate in various biological processes, including regulating transcription and sustaining genome 3D organization. Here, we present a method termed Red-C that exploits proximity ligation to identify contacts with the genome for all RNA molecules present in the nucleus. Using Red-C, we uncovered the RNA-DNA interactome of human K562 cells and identified hundreds of ncRNAs enriched in active or repressed chromatin, including previously undescribed RNAs. Analysis of the RNA-DNA interactome also allowed us to trace the kinetics of messenger RNA production. Our data support the model of co-transcriptional intron splicing, but not the hypothesis of the circularization of actively transcribed genes.

Non-backtracking walks reveal compartments in sparse chromatin interaction networks

Chromatin communities stabilized by protein machinery play essential role in gene regulation and refine global polymeric folding of the chromatin fiber. However, treatment of these communities in the framework of the classical network theory (stochastic block model, SBM) does not take into account intrinsic linear connectivity of the chromatin loci. Here we propose the polymer block model, paving the way for community detection in polymer networks. On the basis of this new model we modify the non-backtracking flow operator and suggest the first protocol for annotation of compartmental domains in sparse single cell Hi-C matrices. In particular, we prove that our approach corresponds to the maximum entropy principle. The benchmark analyses demonstrates that the spectrum of the polymer non-backtracking operator resolves the true compartmental structure up to the theoretical detectability threshold, while all commonly used operators fail above it. We test various operators on real data and conclude that the sizes of the non-backtracking single cell domains are most close to the sizes of compartments from the population data. Moreover, the found domains clearly segregate in the gene density and correlate with the population compartmental mask, corroborating biological significance of our annotation of the chromatin compartmental domains in single cells Hi-C matrices.

[3D Genomics]

The development of new research methods significantly changed our views on the role that the 3D organization of the genome plays in its functional activity. It was found that the genome is subdivided into structural-functional units that restrict the area of enhancer action at the level of spatial organization. Spatial reconfiguration of an extended genomic fragment was identified as a potential mechanism that activates or represses various genes. Accordingly, a distorted spatial organization of the genome often causes various diseases, including cancer. All these observations contributed to the emergence of 3D genomics as a new avenue of research. The review summarizes the most important discoveries in the field of 3D genomics and discusses the directions of its further development.

C-TALE, a new cost-effective method for targeted enrichment of Hi-C/3C-seq libraries

Studies performed using Hi-C and other high-throughput whole-genome C-methods have demonstrated that 3D organization of eukaryotic genomes is functionally relevant. Unfortunately, ultra-deep sequencing of Hi-C libraries necessary to detect loop structures in large vertebrate genomes remains rather expensive. However, many studies are in fact aimed at determining the fine-scale 3D structure of comparatively small genomic regions up to several Mb in length. Such studies typically focus on the spatial structure of domains of coregulated genes, molecular mechanisms of loop formation, and interrogation of functional significance of GWAS-revealed polymorphisms. Therefore, a handful of molecular techniques based on Hi-C have been developed to address such issues. These techniques commonly rely on in-solution hybridization of Hi-C/3C-seq libraries with pools of biotinylated baits covering the region of interest, followed by deep sequencing of the enriched library. Here, we describe a new protocol of this kind, C-TALE (Chromatin TArget Ligation Enrichment). Preparation of hybridization probes from bacterial artificial chromosomes and an additional round of enrichment make C-TALE a cost-effective alternative to existing many-versus-all C-methods.

CRISPR/Cas9-generated mouse model of Duchenne muscular dystrophy recapitulating a newly identified large 430 kb deletion in the human DMD gene

Exon skipping is a promising strategy for Duchenne muscular dystrophy (DMD) disease-modifying therapy. To make this approach safe, ensuring that excluding one or more exons will restore the reading frame and that the resulting protein will retain critical functions of the full-length dystrophin protein is necessary. However, in vivo testing of the consequences of skipping exons that encode the N-terminal actin-binding domain (ABD) has been confounded by the absence of a relevant animal model. We created a mouse model of the disease recapitulating a novel human mutation, a large de novo deletion of exons 8-34 of the DMD gene, found in a Russian DMD patient. This mutation was achieved by deleting exons 8-34 of the X-linked mouse D md gene using CRISPR/Cas9 genome editing, which led to a reading frame shift and the absence of functional dystrophin production. Male mice carrying this deletion display several important signs of muscular dystrophy, including a gradual age-dependent decrease in muscle strength, increased creatine kinase, muscle fibrosis and central nucleation. The degrees of these changes are comparable to those observed in mdx mice, a standard laboratory model of DMD. This new model of DMD will be useful for validating therapies based on skipping exons that encode the N-terminal ABD and for improving our understanding of the role of the N-terminal domain and central rod domain in the biological function of dystrophin. Simultaneous skipping of exons 6 and 7 should restore the gene reading frame and lead to the production of a protein that might retain functionality despite the partial deletion of the ABD.

The Nuclear Lamina as an Organizer of Chromosome Architecture

The nuclear lamina (NL) is a meshwork of lamins and lamin-associated proteins adjoining the inner side of the nuclear envelope. In early embryonic cells, the NL mainly suppresses background transcription, whereas, in differentiated cell types, its disruption affects gene expression more severely. Normally, the NL serves as a backbone for multiple chromatin anchoring sites, thus shaping the spatial organization of chromosomes in the interphase nucleus. However, upon cell senescence, aging, or in some types of terminally differentiated cells and lamin-associated diseases, the loss of NL-chromatin tethering causes drastic alterations in chromosome architecture. Here, we provide an overview of the recent advances in the field of NL-chromatin interactions, focusing on their impact on chromatin positioning, compaction, repression, and spatial organization.

Quantitative differences in TAD border strength underly the TAD hierarchy in Drosophila chromosomes

Chromosomes in many organisms, including Drosophila and mammals, are folded into topologically associating domains (TADs). Increasing evidence suggests that TAD folding is hierarchical, wherein subdomains combine to form larger superdomains, instead of a sequence of nonoverlapping domains. Here, we studied the hierarchical structure of TADs in Drosophila. We show that the boundaries of TADs of different hierarchical levels are characterized by the presence of different portions of active chromatin, but do not vary in the binding of architectural proteins, such as CCCTC binding factor or cohesin. The apparent hierarchy of TADs in Drosophila chromosomes is not likely to have functional importance but rather reflects various options of long-range chromatin folding directed by the distribution of active and inactive chromatin segments and may represent population average.

Genetic and Epigenetic Mechanisms of β-Globin Gene Switching

Vertebrates have multiple forms of hemoglobin that differ in the composition of their polypeptide chains. During ontogenesis, the composition of these subunits changes. Genes encoding different α- and β-polypeptide chains are located in two multigene clusters on different chromosomes. Each cluster contains several genes that are expressed at different stages of ontogenesis. The phenomenon of stage-specific transcription of globin genes is referred to as globin gene switching. Mechanisms of expression switching, stage-specific activation, and repression of transcription of α- and β-globin genes are of interest from both theoretical and practical points of view. Alteration of balanced expression of globin genes, which usually occurs due to damage to adult β-globin genes, leads to development of severe diseases - hemoglobinopathies. In most cases, reactivation of the fetal hemoglobin gene in patients with β-thalassemia and sickle cell disease can reduce negative consequences of irreversible alterations of expression of the β-globin genes. This review focuses on the current state of research on genetic and epigenetic mechanisms underlying stage-specific switching of β-globin genes.

Gene functioning and storage within a folded genome

In mammals, genomic DNA that is roughly 2 m long is folded to fit the size of the cell nucleus that has a diameter of about 10 μm. The folding of genomic DNA is mediated via assembly of DNA-protein complex, chromatin. In addition to the reduction of genomic DNA linear dimensions, the assembly of chromatin allows to discriminate and to mark active (transcribed) and repressed (non-transcribed) genes. Consequently, epigenetic regulation of gene expression occurs at the level of DNA packaging in chromatin. Taking into account the increasing attention of scientific community toward epigenetic systems of gene regulation, it is very important to understand how DNA folding in chromatin is related to gene activity. For many years the hierarchical model of DNA folding was the most popular. It was assumed that nucleosome fiber (10-nm fiber) is folded into 30-nm fiber and further on into chromatin loops attached to a nuclear/chromosome scaffold. Recent studies have demonstrated that there is much less regularity in chromatin folding within the cell nucleus. The very existence of 30-nm chromatin fibers in living cells was questioned. On the other hand, it was found that chromosomes are partitioned into self-interacting spatial domains that restrict the area of enhancers action. Thus, TADs can be considered as structural-functional domains of the chromosomes. Here we discuss the modern view of DNA packaging within the cell nucleus in relation to the regulation of gene expression. Special attention is paid to the possible mechanisms of the chromatin fiber self-assembly into TADs. We discuss the model postulating that partitioning of the chromosome into TADs is determined by the distribution of active and inactive chromatin segments along the chromosome.

This article was specially invited by the editors and represents work by leading researchers.

Activation of the alpha-globin gene expression correlates with dramatic upregulation of nearby non-globin genes and changes in local and large-scale chromatin spatial structure

BACKGROUND

In homeotherms, the alpha-globin gene clusters are located within permanently open genome regions enriched in housekeeping genes. Terminal erythroid differentiation results in dramatic upregulation of alpha-globin genes making their expression comparable to the rRNA transcriptional output. Little is known about the influence of the erythroid-specific alpha-globin gene transcription outburst on adjacent, widely expressed genes and large-scale chromatin organization. Here, we have analyzed the total transcription output, the overall chromatin contact profile, and CTCF binding within the 2.7 Mb segment of chicken chromosome 14 harboring the alpha-globin gene cluster in cultured lymphoid cells and cultured erythroid cells before and after induction of terminal erythroid differentiation.

RESULTS

We found that, similarly to mammalian genome, the chicken genomes is organized in TADs and compartments. Full activation of the alpha-globin gene transcription in differentiated erythroid cells is correlated with upregulation of several adjacent housekeeping genes and the emergence of abundant intergenic transcription. An extended chromosome region encompassing the alpha-globin cluster becomes significantly decompacted in differentiated erythroid cells, and depleted in CTCF binding and CTCF-anchored chromatin loops, while the sub-TAD harboring alpha-globin gene cluster and the upstream major regulatory element (MRE) becomes highly enriched with chromatin interactions as compared to lymphoid and proliferating erythroid cells. The alpha-globin gene domain and the neighboring loci reside within the A-like chromatin compartment in both lymphoid and erythroid cells and become further segregated from the upstream gene desert upon terminal erythroid differentiation.

CONCLUSIONS

Our findings demonstrate that the effects of tissue-specific transcription activation are not restricted to the host genomic locus but affect the overall chromatin structure and transcriptional output of the encompassing topologically associating domain.

Single-nucleus Hi-C reveals unique chromatin reorganization at oocyte-to-zygote transition

Chromatin is reprogrammed after fertilization to produce a totipotent zygote with the potential to generate a new organism. The maternal genome inherited from the oocyte and the paternal genome provided by sperm coexist as separate haploid nuclei in the zygote. How these two epigenetically distinct genomes are spatially organized is poorly understood. Existing chromosome conformation capture-based methods are not applicable to oocytes and zygotes owing to a paucity of material. To study three-dimensional chromatin organization in rare cell types, we developed a single-nucleus Hi-C (high-resolution chromosome conformation capture) protocol that provides greater than tenfold more contacts per cell than the previous method. Here we show that chromatin architecture is uniquely reorganized during the oocyte-to-zygote transition in mice and is distinct in paternal and maternal nuclei within single-cell zygotes. Features of genomic organization including compartments, topologically associating domains (TADs) and loops are present in individual oocytes when averaged over the genome, but the presence of each feature at a locus varies between cells. At the sub-megabase level, we observed stochastic clusters of contacts that can occur across TAD boundaries but average into TADs. Notably, we found that TADs and loops, but not compartments, are present in zygotic maternal chromatin, suggesting that these are generated by different mechanisms. Our results demonstrate that the global chromatin organization of zygote nuclei is fundamentally different from that of other interphase cells. An understanding of this zygotic chromatin 'ground state' could potentially provide insights into reprogramming cells to a state of totipotency.

Topologically-associating domains: gene warehouses adapted to serve transcriptional regulation

Structural-functional domains have long been hypothesized to occur in eukaryotic chromosomes, but their existence still remains controversial. Here, we discuss the current state of studies of 3D genome folding and the relation of this folding to the functional organization of the genome.

Active chromatin and transcription play a key role in chromosome partitioning into topologically associating domains

Recent advances enabled by the Hi-C technique have unraveled many principles of chromosomal folding that were subsequently linked to disease and gene regulation. In particular, Hi-C revealed that chromosomes of animals are organized into topologically associating domains (TADs), evolutionary conserved compact chromatin domains that influence gene expression. Mechanisms that underlie partitioning of the genome into TADs remain poorly understood. To explore principles of TAD folding in Drosophila melanogaster, we performed Hi-C and poly(A)⁺ RNA-seq in four cell lines of various origins (S2, Kc167, DmBG3-c2, and OSC). Contrary to previous studies, we find that regions between TADs (i.e., the inter-TADs and TAD boundaries) in Drosophila are only weakly enriched with the insulator protein dCTCF, while another insulator protein Su(Hw) is preferentially present within TADs. However, Drosophila inter-TADs harbor active chromatin and constitutively transcribed (housekeeping) genes. Accordingly, we find that binding of insulator proteins dCTCF and Su(Hw) predicts TAD boundaries much worse than active chromatin marks do. Interestingly, inter-TADs correspond to decompacted inter-bands of polytene chromosomes, whereas TADs mostly correspond to densely packed bands. Collectively, our results suggest that TADs are condensed chromatin domains depleted in active chromatin marks, separated by regions of active chromatin. We propose the mechanism of TAD self-assembly based on the ability of nucleosomes from inactive chromatin to aggregate, and lack of this ability in acetylated nucleosomal arrays. Finally, we test this hypothesis by polymer simulations and find that TAD partitioning may be explained by different modes of inter-nucleosomal interactions for active and inactive chromatin.

Spatial organization of the chicken beta-globin gene domain in erythroid cells of embryonic and adult lineages

Background

The β-globin gene domains of vertebrate animals constitute popular models for studying the regulation of eukaryotic gene transcription. It has previously been shown that in the mouse the developmental switching of globin gene expression correlates with the reconfiguration of an active chromatin hub (ACH), a complex of promoters of transcribed genes with distant regulatory elements. Although it is likely that observations made in the mouse β-globin gene domain are also relevant for this locus in other species, the validity of this supposition still lacks direct experimental evidence. Here, we have studied the spatial organization of the chicken β-globin gene domain. This domain is of particular interest because it represents the perfect example of the so-called ‘strong’ tissue-specific gene domain flanked by insulators, which delimit the area of preferential sensitivity to DNase I in erythroid cells.

Results

Using chromosome conformation capture (3C), we have compared the spatial configuration of the β-globin gene domain in chicken red blood cells (RBCs) expressing embryonic (3-day-old RBCs) and adult (9-day-old RBCs) β-globin genes. In contrast to observations made in the mouse model, we found that in the chicken, the early embryonic β-globin gene, Ε, did not interact with the locus control region in RBCs of embryonic lineage (3-day RBCs), where this gene is actively transcribed. In contrast to the mouse model, a strong interaction of the promoter of another embryonic β-globin gene, ρ, with the promoter of the adult β-globin gene, β^A, was observed in RBCs from both 3-day and 9-day chicken embryos. Finally, we have demonstrated that insulators flanking the chicken β-globin gene domain from the upstream and from the downstream interact with each other, which places the area characterized by lineage-specific sensitivity to DNase I in a separate chromatin loop.

Conclusions

Taken together, our results strongly support the ACH model but show that within a domain of tissue-specific genes, the active status of a promoter does not necessarily correlate with the recruitment of this promoter to the ACH.

The inactivation of the π gene in chicken erythroblasts of adult lineage is not mediated by packaging of the embryonic part of the α-globin gene domain into a repressive heterochromatin-like structure

The developmental switch of globin gene expression is a characteristic feature of vertebrate organisms. The switch of β-globin expression is believed to depend on reconfiguration of the active chromatin hub, which contains transcribed genes and regulatory elements. Mechanisms controlling the switch of α-globin gene expression are less clear. Here, we studied the mode of chromatin packaging of the chicken α-globin gene domain in red blood cells (RBCs) of primitive and definite lineages and the spatial configuration of this domain in RBCs of primitive lineage. It has been demonstrated that RBCs of primitive lineage already contain the adult-type active chromatin hub but the embryonal α-type globin π gene is not recruited to this hub. Distribution of active and repressive histone modifications over the α-globin gene domain in RBCs of definite and primitive lineages does not corroborate the hypothesis that inactivation of the π gene in RBCs of adult lineage is mediated via formation of a local repressed chromatin domain. This conclusion is supported by the demonstration that in chicken erythroblasts of adult lineage, the embryonal and adult segments of the α-globin gene domain show similar elevated sensitivities to DNase I.