Principles of Chromosome Architecture Revealed by Hi-C

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.

3D Genome Reconstruction from Partially Phased Hi-C Data

The 3-dimensional (3D) structure of the genome is of significant importance for many cellular processes. In this paper, we study the problem of reconstructing the 3D structure of chromosomes from Hi-C data of diploid organisms, which poses additional challenges compared to the better-studied haploid setting. With the help of techniques from algebraic geometry, we prove that a small amount of phased data is sufficient to ensure finite identifiability, both for noiseless and noisy data. In the light of these results, we propose a new 3D reconstruction method based on semidefinite programming, paired with numerical algebraic geometry and local optimization. The performance of this method is tested on several simulated datasets under different noise levels and with different amounts of phased data. We also apply it to a real dataset from mouse X chromosomes, and we are then able to recover previously known structural features.

Locus folding mechanisms determine modes of antigen receptor gene assembly

The dynamic folding of genomes regulates numerous biological processes, including antigen receptor (AgR) gene assembly. We show that, unlike other AgR loci, homotypic chromatin interactions and bidirectional chromosome looping both contribute to structuring Tcrb for efficient long-range V(D)J recombination. Inactivation of the CTCF binding element (CBE) or promoter at the most 5'Vβ segment (Trbv1) impaired loop extrusion originating locally and extending to DβJβ CBEs at the opposite end of Tcrb. Promoter or CBE mutation nearly eliminated Trbv1 contacts and decreased RAG endonuclease-mediated Trbv1 recombination. Importantly, Trbv1 rearrangement can proceed independent of substrate orientation, ruling out scanning by DβJβ-bound RAG as the sole mechanism of Vβ recombination, distinguishing it from Igh. Our data indicate that CBE-dependent generation of loops cooperates with promoter-mediated activation of chromatin to juxtapose Vβ and DβJβ segments for recombination through diffusion-based synapsis. Thus, the mechanisms that fold a genomic region can influence molecular processes occurring in that space, which may include recombination, repair, and transcriptional programming.

Hold out the genome: a roadmap to solving the cis-regulatory code

Gene expression is regulated by transcription factors that work together to read cis-regulatory DNA sequences. The 'cis-regulatory code' - how cells interpret DNA sequences to determine when, where and how much genes should be expressed - has proven to be exceedingly complex. Recently, advances in the scale and resolution of functional genomics assays and machine learning have enabled substantial progress towards deciphering this code. However, the cis-regulatory code will probably never be solved if models are trained only on genomic sequences; regions of homology can easily lead to overestimation of predictive performance, and our genome is too short and has insufficient sequence diversity to learn all relevant parameters. Fortunately, randomly synthesized DNA sequences enable testing a far larger sequence space than exists in our genomes, and designed DNA sequences enable targeted queries to maximally improve the models. As the same biochemical principles are used to interpret DNA regardless of its source, models trained on these synthetic data can predict genomic activity, often better than genome-trained models. Here we provide an outlook on the field, and propose a roadmap towards solving the cis-regulatory code by a combination of machine learning and massively parallel assays using synthetic DNA.

An Atlas of Promoter Chromatin Modifications and HiChIP Regulatory Interactions in Human Subcutaneous Adipose-Derived Stem Cells

The genome of human adipose-derived stem cells (ADSCs) from abdominal and gluteofemoral adipose tissue depots are maintained in depot-specific stable epigenetic conformations that influence cell-autonomous gene expression patterns and drive unique depot-specific functions. The traditional approach to explore tissue-specific transcriptional regulation has been to correlate differential gene expression to the nearest-neighbor linear-distance regulatory region defined by associated chromatin features including open chromatin status, histone modifications, and DNA methylation. This has provided important information; nonetheless, the approach is limited because of the known organization of eukaryotic chromatin into a topologically constrained three-dimensional network. This network positions distal regulatory elements in spatial proximity with gene promoters which are not predictable based on linear genomic distance. In this work, we capture long-range chromatin interactions using HiChIP to identify remote genomic regions that influence the differential regulation of depot-specific genes in ADSCs isolated from different adipose depots. By integrating these data with RNA-seq results and histone modifications identified by ChIP-seq, we uncovered distal regulatory elements that influence depot-specific gene expression in ADSCs. Interestingly, a subset of the HiChIP-defined chromatin loops also provide previously unknown connections between waist-to-hip ratio GWAS variants with genes that are known to significantly influence ADSC differentiation and adipocyte function.

A Lightweight Framework For Chromatin Loop Detection at the Single-Cell Level

Single-cell Hi-C (scHi-C) has made it possible to analyze chromatin organization at the single-cell level. However, scHi-C experiments generate inherently sparse data, which poses a challenge for loop calling methods. The existing approach performs significance tests across the imputed dense contact maps, leading to substantial computational overhead and loss of information at the single-cell level. To overcome this limitation, a lightweight framework called scGSLoop is proposed, which sets a new paradigm for scHi-C loop calling by adapting the training and inferencing strategies of graph-based deep learning to leverage the sequence features and 1D positional information of genomic loci. With this framework, sparsity is no longer a challenge, but rather an advantage that the model leverages to achieve unprecedented computational efficiency. Compared to existing methods, scGSLoop makes more accurate predictions and is able to identify more loops that have the potential to play regulatory roles in genome functioning. Moreover, scGSLoop preserves single-cell information by identifying a distinct group of loops for each individual cell, which not only enables an understanding of the variability of chromatin looping states between cells, but also allows scGSLoop to be extended for the investigation of multi-connected hubs and their underlying mechanisms.

The human long non-coding RNA LINC00941 and its modes of action in health and disease

Long non-coding RNAs have gained attention in recent years as they were shown to play crucial roles in the regulation of cellular processes, but the understanding of the exact mechanisms is still incomplete in most cases. This is also true for long non-coding RNA LINC00941, which was recently found to be highly upregulated in various types of cancer influencing cell proliferation and metastasis. Initial studies could not elucidate the mode of action to understand the role and real impact of LINC00941 in tissue homeostasis and cancer development. However, recent analyses have demonstrated multiple potential modes of action of LINC00941 influencing the functionality of various cancer cell types. Correspondingly, LINC00941 was proposed to be involved in regulation of mRNA transcription and modulation of protein stability, respectively. In addition, several experimental approaches suggest a function of LINC00941 as competitive endogenous RNA, thus acting in a post-transcriptional regulatory fashion. This review summarizes our recent knowledge about the mechanisms of action of LINC00941 elucidated so far and discusses its putative role in miRNA sequestering processes. In addition, the functional role of LINC00941 in regulating human keratinocytes is discussed to also highlight its role in normal tissue homeostasis tissue aside from its involvement in cancer.

METTL3-stabilized super enhancers-lncRNA SUCLG2-AS1 mediates the formation of a long-range chromatin loop between enhancers and promoters of SOX2 in metastasis and radiosensitivity of nasopharyngeal carcinoma

BACKGROUND

Super enhancers (SE) play pivotal roles in cell identity and diseases occur including tumorigenesis. The depletion of SE-associated lncRNA transcripts, also known as super-lncRNA, causes the activity of SE to be dysregulated.

METHODS

We screened and identified an elevated metastasis-associated SE-lncRNA SUCLG2-AS1 in nasopharyngeal carcinoma (NPC) using RNA-sequencing, real-time quantitative polymerase chain reaction (RT-qPCR) and bioinformatics. Western blotting, RT-qPCR, methylated RNA immunoprecipitation (MeRIP), RNA immunoprecipitation, chromatin immunoprecipitation, RNA pull-down and 3C (chromosome conformation capture assays) were used for mechanistic studies.

RESULTS

SUCLG2-AS1 was correlated with a poor prognosis. SUCLG2-AS1 promotes NPC cell invasion and metastasis while repressing apoptosis and radiosensitivity in vitro and in vivo. Mechanistically, high SUCLG2-AS1 expression occurred in an m6A-dependent manner. SUCLG2-AS1 was found to be located in the SE region of SOX2, and it regulated the expression of SOX2 via long-range chromatin loop formation, which via mediating CTCF (transcription factor) occupied the SE and promoter region of SOX2, thus regulating the metastasis and radiosensitivity of NPC.

CONCLUSIONS

Taken together, our data suggest that SUCLG2-AS1 may serve as a novel intervention target for the clinical treatment of NPC.

Fundamental insights into the correlation between chromosome configuration and transcription

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

MS-ACGAN: A modified auxiliary classifier generative adversarial network for schizophrenia's samples augmentation based on microarray gene expression data

Artificial intelligence-based models and robust computational methods have expedited the data-to-knowledge trajectory in precision medicine. Although machine learning models have been widely applied in medical data analysis, some barriers are yet to be challenging, such as available biosample shortage, prohibitive costs, rare diseases, and ethical considerations. Transcriptomics, an omics approach that studies gene activities and provides gene expression data such as microarray and RNA-Sequences faces the difficulties of biospecimen collection, particularly for mental disorders, as some psychiatric patients avoid medical care. Microarray data suffers from the low number of available samples, making it challenging to apply machine learning models. However, adversarial generative network (GAN), the hottest paradigm in deep learning, has created unprecedented momentum in data augmentation and efficiently expands datasets. This paper proposes a novel model termed MS-ACGAN, where the generator feeds on a bordered Gaussian distribution. In machine learning, calibration is of utmost importance, which gives insight into model uncertainty and is considered a crucial step toward improving the robustness and reliability of models. Therefore, we apply calibration techniques to classifiers and focus on estimating their probabilities as accurately as possible. Additionally, we present our trustworthy outputs by harnessing confidence intervals that confine the point estimate limitations and report a range of expected values for performance metrics. Both concepts statistically describe the implemented model's reliability in this study. Furthermore, we employ two quantitative measures, GAN-train and GAN-test, to demonstrate that the artificial data generated by our robust approach remarkably resembles the original data characteristics.

Assignment of the somatic A/B compartments to chromatin domains in giant transcriptionally active lampbrush chromosomes

BACKGROUND

The three-dimensional configuration of the eukaryotic genome is an emerging area of research. Chromosome conformation capture outlined genome segregation into large scale A and B compartments corresponding mainly to transcriptionally active and repressive chromatin. It remains unknown how the compartmentalization of the genome changes in growing oocytes of animals with hypertranscriptional type of oogenesis. Such oocytes are characterized by highly elongated chromosomes, called lampbrush chromosomes, which acquire a typical chromomere-loop appearance, representing one of the classical model systems for exploring the structural and functional organization of chromatin domains.

RESULTS

Here, we compared the distribution of A/B compartments in chicken somatic cells with chromatin domains in lampbrush chromosomes. We found that in lampbrush chromosomes, the extended chromatin domains, restricted by compartment boundaries in somatic cells, disintegrate into individual chromomeres. Next, we performed FISH-mapping of the genomic loci, which belong to A or B chromatin compartments as well as to A/B compartment transition regions in embryonic fibroblasts on isolated lampbrush chromosomes. We found, that in chicken lampbrush chromosomes, clusters of dense compact chromomeres bearing short lateral loops and enriched with repressive epigenetic modifications generally correspond to constitutive B compartments in somatic cells. A compartments align with lampbrush chromosome segments with smaller, less compact chromomeres, longer lateral loops, and a higher transcriptional status. Clusters of small loose chromomeres with relatively long lateral loops show no obvious correspondence with either A or B compartment identity. Some genes belonging to facultative B (sub-) compartments can be tissue-specifically transcribed during oogenesis, forming distinct lateral loops.

CONCLUSIONS

Here, we established a correspondence between the A/B compartments in somatic interphase nucleus and chromatin segments in giant lampbrush chromosomes from diplotene stage oocytes. The chromomere-loop structure of the genomic regions corresponding to interphase A and B compartments reveals the difference in how they are organized at the level of chromatin domains. The results obtained also suggest that gene-poor regions tend to be packed into chromomeres.

High-quality chromosome-level de novo assembly of the Trifolium repens

BACKGROUND

White clover (Trifolium repens L.), an excellent perennial legume forage, is an allotetraploid native to southeastern Europe and southern Asia. It has high nutritional, ecological, genetic breeding, and medicinal values and exhibits excellent resistance to cold, drought, trample, and weed infestation. Thus, white clover is widely planted in Europe, America, and China; however, the lack of reference genome limits its breeding and cultivation. This study generated a white clover de novo genome assembly at the chromosomal level and annotated its components.

RESULTS

The PacBio third-generation Hi-Fi assembly and sequencing methods generated a 1096 Mb genome size of T. repens, with contigs of N50 = 14 Mb, scaffolds of N50 = 65 Mb, and BUSCO value of 98.5%. The newly assembled genome has better continuity and integrity than the previously reported white clover reference genome; thus provides important resources for the molecular breeding and evolution of white clover and other forage. Additionally, we annotated 90,128 high-confidence gene models from the genome. White clover was closely related to Trifolium pratense and Trifolium medium but distantly related to Glycine max, Vigna radiata, Medicago truncatula, and Cicer arietinum. The expansion, contraction, and GO functional enrichment analysis of the gene families showed that T. repens gene families were associated with biological processes, molecular function, cellular components, and environmental resistance, which explained its excellent agronomic traits.

CONCLUSIONS

This study reports a high-quality de novo assembly of white clover genome obtained at the chromosomal level using PacBio Hi-Fi sequencing, a third-generation sequencing. The generated high-quality genome assembly of white clover provides a key basis for accelerating the research and molecular breeding of this important forage crop. The genome is also valuable for future studies on legume forage biology, evolution, and genome-wide mapping of quantitative trait loci associated with the relevant agronomic traits.

Regulation of loop extrusion on the interphase genome

In the human cell nucleus, dynamically organized chromatin is the substrate for gene regulation, DNA replication, and repair. A central mechanism of DNA loop formation is an ATPase motor cohesin-mediated loop extrusion. The cohesin complexes load and unload onto the chromosome under the control of other regulators that physically interact and affect motor activity. Regulation of the dynamic loading cycle of cohesin influences not only the chromatin structure but also genome-associated human disorders and aging. This review focuses on the recently spotlighted genome organizing factors and the mechanism by which their dynamic interactions shape the genome architecture in interphase.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

Double-strand break repair and mis-repair in 3D

DNA double-strand breaks (DSBs) are lesions that arise frequently from exposure to damaging agents as well as from ongoing physiological DNA transactions. Mis-repair of DSBs leads to rearrangements and structural variations in chromosomes, including insertions, deletions, and translocations implicated in disease. The DNA damage response (DDR) limits pathologic mutations and large-scale chromosome rearrangements. DSB repair initiates in 2D at DNA lesions with the stepwise recruitment of repair proteins and local chromatin remodeling which facilitates break accessibility. More complex structures are then formed via protein assembly into nanodomains and via genome folding into chromatin loops. Subsequently, 3D reorganization of DSBs is guided by clustering forces which drive the assembly of repair domains harboring multiple lesions. These domains are further stabilized and insulated into condensates via liquid-liquid phase-separation. Here, we discuss the benefits and risks associated with this 3D reorganization of the broken genome.

Next generation sequencing technologies to explore the diversity of germplasm resources: achievements and trends in tomato

Tomato is one of the major vegetable crops grown worldwide and a model species for genetic and biological research. Progress in genomic technologies made possible the development of forefront methods for high-scale sequencing, providing comprehensive insight into the genetic architecture of germplasm resources. This review revisits next-generation sequencing strategies and applications to investigate the diversity of tomato, describing the common platforms used for SNP genotyping of large collections, de novo sequencing, and whole genome resequencing. Significant findings in evolutionary history are outlined, thus discussing how genomics has provided new hints about the processes behind domestication. Finally, achievement and perspectives on pan-genome construction and graphical pan-genome development toward precise mining of the natural variation to be exploited for breeding purposes are presented.

Sample-multiplexing approaches for single-cell sequencing

Single-cell sequencing is widely used in biological and medical studies. However, its application with multiple samples is hindered by inefficient sample processing, high experimental costs, ambiguous identification of true single cells, and technical batch effects. Here, we introduce sample-multiplexing approaches for single-cell sequencing in transcriptomics, epigenomics, genomics, and multiomics. In single-cell transcriptomics, sample multiplexing uses variants of native or artificial features as sample markers, enabling sample pooling and decoding. Such features include: (1) natural genetic variation, (2) nucleotide-barcode anchoring on cellular or nuclear membranes, (3) nucleotide-barcode internalization to the cytoplasm or nucleus, (4) vector-based barcode expression in cells, and (5) nucleotide-barcode incorporation during library construction. Other single-cell omics methods are based on similar concepts, particularly single-cell combinatorial indexing. These methods overcome current challenges, while enabling super-loading of single cells. Finally, selection guidelines are presented that can accelerate technological application.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

Chromosome-level genome assembly of the black widow spider Latrodectus elegans illuminates composition and evolution of venom and silk proteins

BACKGROUND

The black widow spider has both extraordinarily neurotoxic venom and three-dimensional cobwebs composed of diverse types of silk. However, a high-quality reference genome for the black widow spider was still unavailable, which hindered deep understanding and application of the valuable biomass.

FINDINGS

We assembled the Latrodectus elegans genome, including a genome size of 1.57 Gb with contig N50 of 4.34 Mb and scaffold N50 of 114.31 Mb. Hi-C scaffolding assigned 98.08% of the genome to 14 pseudo-chromosomes, and with BUSCO, completeness analysis revealed that 98.4% of the core eukaryotic genes were completely present in this genome. Annotation of this genome identified that repetitive sequences account for 506.09 Mb (32.30%) and 20,167 protein-coding genes, and specifically, we identified 55 toxin genes and 26 spidroins and provide preliminary analysis of their composition and evolution.

CONCLUSIONS

We present the first chromosome-level genome assembly of a black widow spider and provide substantial toxin and spidroin gene resources. These high-qualified genomic data add valuable resources from a representative spider group and contribute to deep exploration of spider genome evolution, especially in terms of the important issues on the diversification of venom and web-weaving pattern. The sequence data are also firsthand templates for further application of the spider biomass.

Detecting TAD-like domains from RNA-associated interactions

Topologically associated domains (TADs) are crucial chromatin structural units. Evidence has illustrated that RNA–chromatin and RNA–RNA spatial interactions, so-called RNA-associated interactions (RAIs), may be associated with TAD-like domains (TLDs). To decode hierarchical TLDs from RAIs, we proposed SuperTLD, a domain detection algorithm incorporating imputation. We applied SuperTLD on four RAI data sets and compared TLDs with the TADs identified from the corresponding Hi-C datasets. The TLDs and TADs share a moderate similarity of hierarchies ≥ 0.5312 and the finest structures ≥ 0.8295. Comparison between boundaries and domains further demonstrated the novelty of TLDs. Enrichment analysis of epigenetic characteristics illustrated that the novel TLDs exhibit an enriched CTCF by 0.6245 fold change and H3 histone marks enriched within domains. GO analysis on the TLD novel boundaries exhibited enriched diverse terms, revealing TLDs’ formation mechanism related closely to gene regulation.

Persistent spectral simplicial complex-based machine learning for chromosomal structural analysis in cellular differentiation

The three-dimensional (3D) chromosomal structure plays an essential role in all DNA-templated processes, including gene transcription, DNA replication and other cellular processes. Although developing chromosome conformation capture (3C) methods, such as Hi-C, which can generate chromosomal contact data characterized genome-wide chromosomal structural properties, understanding 3D genomic nature-based on Hi-C data remains lacking. Here, we propose a persistent spectral simplicial complex (PerSpectSC) model to describe Hi-C data for the first time. Specifically, a filtration process is introduced to generate a series of nested simplicial complexes at different scales. For each of these simplicial complexes, its spectral information can be calculated from the corresponding Hodge Laplacian matrix. PerSpectSC model describes the persistence and variation of the spectral information of the nested simplicial complexes during the filtration process. Different from all previous models, our PerSpectSC-based features provide a quantitative global-scale characterization of chromosome structures and topology. Our descriptors can successfully classify cell types and also cellular differentiation stages for all the 24 types of chromosomes simultaneously. In particular, persistent minimum best characterizes cell types and Dim (1) persistent multiplicity best characterizes cellular differentiation. These results demonstrate the great potential of our PerSpectSC-based models in polymeric data analysis.

Plant In Situ Hi-C Experimental Protocol and Bioinformatic Analysis

Hi-C enables the characterization of the 0conformation of the genome in the three-dimensional nuclear space. This technique has revolutionized our ability to detect interactions between linearly distant genomic sites on a genome-wide scale. Here, we detail a protocol to carry out in situ Hi-C in plants and describe a straightforward bioinformatics pipeline for the analysis of such data, in particular for comparing samples from different organs or conditions.

Advances in profiling chromatin architecture shed light on the regulatory dynamics underlying brain disorders

Understanding the exquisitely complex nature of the three-dimensional organization of the genome and how it affects gene regulation remains a central question in biology. Recent advances in sequencing- and imaging-based approaches in decoding the three-dimensional chromatin landscape have enabled a systematic characterization of gene regulatory architecture. In this review, we outline how chromatin architecture provides a reference atlas to predict the functional consequences of non-coding variants associated with human traits and disease. High-throughput perturbation assays such as massively parallel reporter assays (MPRA) and CRISPR-based genome engineering in combination with a reference atlas opened an avenue for going beyond observational studies to experimentally validating the regulatory principles of the genome. We conclude by providing a suggested path forward by calling attention to barriers that can be addressed for a more complete understanding of the regulatory landscape of the human brain.

Replication Labeling Methods for Super-Resolution Imaging of Chromosome Territories and Chromatin Domains

Continuing progress in super-resolution microscopy enables the study of sub-chromosomal chromatin organization in single cells with unprecedented detail. Here we describe refined methods for pulse-chase replication labeling of individual chromosome territories (CTs) and replication domain units in mammalian cell nuclei, with specific focus on their application to three-dimensional structured illumination microscopy (3D-SIM). We provide detailed protocols for highly efficient electroporation-based delivery or scratch loading of cell-impermeable fluorescent nucleotides for live-cell studies. Furthermore, we describe the application of (2'S)-2'-deoxy-2'-fluoro-5-ethynyluridine (F-ara-EdU) and 5-vinyl-2'-deoxyuridine (VdU) for the in situ detection of segregated chromosome territories and sister chromatids with minimized cytotoxic side effects.

Chromatin as self-returning walks: From population to single cell and back

With a growing understanding of the chromatin structure, many efforts remain focused on bridging the gap between what is suggested by population-averaged data and what is visualized for single cells. A popular approach to traversing these scales is to fit a polymer model to Hi-C contact data. However, Hi-C is an average of millions to billions of cells, and each cell may not contain all population-averaged contacts. Therefore, we employ a novel approach of summing individual chromosome trajectories-determined by our Self-Returning Random Walk model-to create populations of cells. We allow single cells to consist of disparate structures and reproduce a variety of experimentally relevant contact maps. We show that the amount of shared topology between cells, and their mechanism of formation, changes the population-averaged structure. Therefore, we present a modeling technique that, with few constraints and little oversight, can be used to understand which single-cell chromatin structures underlie population-averaged behavior.

HiC1Dmetrics: framework to extract various one-dimensional features from chromosome structure data

Eukaryotic genomes are organized in a three-dimensional spatial structure. In this regard, the development of chromosome conformation capture methods has enabled studies of chromosome organization on a genomic scale. Hi-C, the high-throughput chromosome conformation capture method, can reveal a population-averaged, hierarchical chromatin structure. The typical Hi-C analysis uses a two-dimensional (2D) contact matrix that indicates contact frequencies between all possible genomic position pairs. Oftentimes, however, such a 2D matrix is not amenable to handling quantitative comparisons, visualizations and integrations across multiple datasets. Although several one-dimensional (1D) metrics have been proposed to depict structural information in Hi-C data, their effectiveness is still underappreciated. Here, we first review the currently available 1D metrics for individual Hi-C samples or two-sample comparisons and then discuss their validity and suitable analysis scenarios. We also propose several new 1D metrics to identify additional unique features of chromosome structures. We highlight that the 1D metrics are reproducible and robust for comparing and visualizing multiple Hi-C samples. Moreover, we show that 1D metrics can be easily combined with epigenome tracks to annotate chromatin states in greater details. We develop a new framework, called HiC1Dmetrics, to summarize all 1D metrics discussed in this study. HiC1Dmetrics is open-source (github.com/wangjk321/HiC1Dmetrics) and can be accessed from both command-line and web-based interfaces. Our tool constitutes a useful resource for the community of chromosome-organization researchers.

Perspective of the GEMSTONE Consortium on Current and Future Approaches to Functional Validation for Skeletal Genetic Disease Using Cellular, Molecular and Animal-Modeling Techniques

The availability of large human datasets for genome-wide association studies (GWAS) and the advancement of sequencing technologies have boosted the identification of genetic variants in complex and rare diseases in the skeletal field. Yet, interpreting results from human association studies remains a challenge. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary. Multiple unknowns exist for putative causal genes, including cellular localization of the molecular function. Intermediate traits ("endophenotypes"), e.g. molecular quantitative trait loci (molQTLs), are needed to identify mechanisms of underlying associations. Furthermore, index variants often reside in non-coding regions of the genome, therefore challenging for interpretation. Knowledge of non-coding variance (e.g. ncRNAs), repetitive sequences, and regulatory interactions between enhancers and their target genes is central for understanding causal genes in skeletal conditions. Animal models with deep skeletal phenotyping and cell culture models have already facilitated fine mapping of some association signals, elucidated gene mechanisms, and revealed disease-relevant biology. However, to accelerate research towards bridging the current gap between association and causality in skeletal diseases, alternative in vivo platforms need to be used and developed in parallel with the current -omics and traditional in vivo resources. Therefore, we argue that as a field we need to establish resource-sharing standards to collectively address complex research questions. These standards will promote data integration from various -omics technologies and functional dissection of human complex traits. In this mission statement, we review the current available resources and as a group propose a consensus to facilitate resource sharing using existing and future resources. Such coordination efforts will maximize the acquisition of knowledge from different approaches and thus reduce redundancy and duplication of resources. These measures will help to understand the pathogenesis of osteoporosis and other skeletal diseases towards defining new and more efficient therapeutic targets.

ARID1A deficiency weakens BRG1-RAD21 interaction that jeopardizes chromatin compactness and drives liver cancer cell metastasis

ARID1A, encoding a subunit of SWI/SNF chromatin remodeling complex, is widely recognized as a tumor suppressor gene in multiple tumor types including liver cancer. Previous studies have demonstrated that ARID1A deficiency can cause liver cancer metastasis, possibly due to the altered chromatin organization, however the underlying mechanisms remain poorly understood. To address the effect of Arid1a deficiency on chromatin organization, we generated chromatin interaction matrices, and exploited the conformation changes upon Arid1a depletion in hepatocytes. Our results demonstrated that Arid1a deficiency induced A/B compartment switching, topologically associated domain (TAD) remodeling, and decrease of chromatin loops. Further mechanism studies revealed that ATPase BRG1 of SWI/SNF complex could physically interact with RAD21, a structural subunit of chromatin architectural element cohesin; whereas ARID1A deficiency significantly diminished the coupled BRG1-RAD21. Interestingly, the tumor-associated genes within the switched compartments were differentially expressed depending upon Arid1a depletion or not. As a consequence of ARID1A deficiency-induced conformational alteration, the dysregulation of some genes such as PMP22 and GSC, promoted the invasion capacity of liver cancer cells. This study provides an insight into liver cancer tumorigenesis and progression related to ARID1A mutations.

Dissecting and predicting different types of binding sites in nucleic acids based on structural information

The biological functions of DNA and RNA generally depend on their interactions with other molecules, such as small ligands, proteins and nucleic acids. However, our knowledge of the nucleic acid binding sites for different interaction partners is very limited, and identification of these critical binding regions is not a trivial work. Herein, we performed a comprehensive comparison between binding and nonbinding sites and among different categories of binding sites in these two nucleic acid classes. From the structural perspective, RNA may interact with ligands through forming binding pockets and contact proteins and nucleic acids using protruding surfaces, while DNA may adopt regions closer to the middle of the chain to make contacts with other molecules. Based on structural information, we established a feature-based ensemble learning classifier to identify the binding sites by fully using the interplay among different machine learning algorithms, feature spaces and sample spaces. Meanwhile, we designed a template-based classifier by exploiting structural conservation. The complementarity between the two classifiers motivated us to build an integrative framework for improving prediction performance. Moreover, we utilized a post-processing procedure based on the random walk algorithm to further correct the integrative predictions. Our unified prediction framework yielded promising results for different binding sites and outperformed existing methods.

Expanding studies of chromosome structure and function in the era of T2T genomics

The recent accomplishment of a truly complete human genome has afforded a new view of chromosome structure and function that was limited 30 years ago. Here, we discuss the expansion of knowledge from the early cytological studies of the genome to the current high-resolution genomic, epigenetic and functional maps that have been achieved by recent technology and computational advances. These studies have revealed unexpected complexities of genome organization and function and uncovered new views of fundamental chromosomal elements. Comprehensive genomic maps will enable accurate diagnosis of human diseases caused by altered chromosome structure and function, facilitate development of chromosome-based therapies and shape the future of preventative medicine and healthcare.

Chromatin interaction neural network (ChINN): a machine learning-based method for predicting chromatin interactions from DNA sequences

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples.

Disruption of NIPBL/Scc2 in Cornelia de Lange Syndrome provokes cohesin genome-wide redistribution with an impact in the transcriptome

Cornelia de Lange syndrome (CdLS) is a rare disease affecting multiple organs and systems during development. Mutations in the cohesin loader, NIPBL/Scc2, were first described and are the most frequent in clinically diagnosed CdLS patients. The molecular mechanisms driving CdLS phenotypes are not understood. In addition to its canonical role in sister chromatid cohesion, cohesin is implicated in the spatial organization of the genome. Here, we investigate the transcriptome of CdLS patient-derived primary fibroblasts and observe the downregulation of genes involved in development and system skeletal organization, providing a link to the developmental alterations and limb abnormalities characteristic of CdLS patients. Genome-wide distribution studies demonstrate a global reduction of NIPBL at the NIPBL-associated high GC content regions in CdLS-derived cells. In addition, cohesin accumulates at NIPBL-occupied sites at CpG islands potentially due to reduced cohesin translocation along chromosomes, and fewer cohesin peaks colocalize with CTCF.

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.

Integrative analysis of genomic and epigenomic data reveals underlying super-enhancer-mediated microRNA regulatory network for human bone mineral density

Bone mineral density (BMD) is a highly heritable complex trait and is a key indicator for diagnosis and treatment for osteoporosis. In the last decade, numerous susceptibility loci for BMD and fracture have been identified by genome wide association studies (GWAS); however, fine mapping of these loci is challengeable. Here, we proposed a new long-range fine-mapping approach that combined super-enhancers (SEs) and microRNAs (miRNAs) data, which were two important factors in control of cell identity and specific differentiation, with the GWAS summary datasets in cell-type-restricted way. Genome-wide SE-based analysis found that the BMD-related variants were significantly enriched in the osteoblast SE regions, indicative of potential long-range effects of such SNPs. With the SNP-mapped SEs (mSEs), 13 accessible long-range mSE-interacted miRNAs (mSE-miRNAs) were identified by integrating osteoblast Hi-C and ATAC-seq data, including three known bone-related miRNAs (miR-132-3p, miR-212-3p and miR-125b-5p). The putative targets of the two newly identified mSE-miRNAs (miR-548aj-3p and miR-190a-3p) were found largely enriched in osteogenic-related pathway and processes, suggesting that these mSE-miRNAs could be functional in the regulation of osteoblast differentiation. Furthermore, we identified 54 genes with the long-range 'mSE-miRNA' approach, and 24 of them were previously reported to be related to skeletal development. Besides, enrichment analysis found that these genes were specifically enriched in the post-transcriptional regulation and bone formation processes. This study provided a new insight into the approach of fine-mapping of GWAS loci. A tool was provided for the genome-wide SE-based analysis and the detection of long-range osteoblast-restricted mSE-miRNAs (https://github.com/Zheng-Lab-Westlake/Osteo-Fine-Mapp-SNP2SE2miRNA).

DNA-Topology Simplification by Topoisomerases

The topological properties of DNA molecules, supercoiling, knotting, and catenation, are intimately connected with essential biological processes, such as gene expression, replication, recombination, and chromosome segregation. Non-trivial DNA topologies present challenges to the molecular machines that process and maintain genomic information, for example, by creating unwanted DNA entanglements. At the same time, topological distortion can facilitate DNA-sequence recognition through localized duplex unwinding and longer-range loop-mediated interactions between the DNA sequences. Topoisomerases are a special class of essential enzymes that homeostatically manage DNA topology through the passage of DNA strands. The activities of these enzymes are generally investigated using circular DNA as a model system, in which case it is possible to directly assay the formation and relaxation of DNA supercoils and the formation/resolution of knots and catenanes. Some topoisomerases use ATP as an energy cofactor, whereas others act in an ATP-independent manner. The free energy of ATP hydrolysis can be used to drive negative and positive supercoiling or to specifically relax DNA topologies to levels below those that are expected at thermodynamic equilibrium. The latter activity, which is known as topology simplification, is thus far exclusively associated with type-II topoisomerases and it can be understood through insight into the detailed non-equilibrium behavior of type-II enzymes. We use a non-equilibrium topological-network approach, which stands in contrast to the equilibrium models that are conventionally used in the DNA-topology field, to gain insights into the rates that govern individual transitions between topological states. We anticipate that our quantitative approach will stimulate experimental work and the theoretical/computational modeling of topoisomerases and similar enzyme systems.

'Enhancing' red cell fate through epigenetic mechanisms

PURPOSE OF REVIEW

Transcription of erythroid-specific genes is regulated by the three-dimensional (3D) structure and composition of chromatin, which dynamically changes during erythroid differentiation. Chromatin organization and dynamics are regulated by several epigenetic mechanisms involving DNA (de-)methylation, posttranslational modifications (PTMs) of histones, chromatin-associated structural proteins, and higher-order structural changes and interactions. This review addresses examples of recent developments in several areas delineating the interface of chromatin regulation and erythroid-specific lineage transcription.

RECENT FINDINGS

We survey and discuss recent studies that focus on the erythroid chromatin landscape, erythroid enhancer-promotor interactions, super-enhancer functionality, the role of chromatin modifiers and epigenetic crosstalk, as well as the progress in mapping red blood cell (RBC) trait-associated genetic variants within cis-regulatory elements (CREs) identified in genome-wide association study (GWAS) efforts as a step toward determining their impact on erythroid-specific gene expression.

SUMMARY

As one of the best characterized and accessible cell differentiation systems, erythropoiesis has been at the forefront of studies aiming to conceptualize how chromatin dynamics regulate transcription. New emerging technologies that bring a significantly enhanced spatial and temporal resolution of chromatin structure, and allow investigation of small cell numbers, have advanced our understanding of chromatin dynamics during erythroid differentiation in vivo.

Functional impact of cancer-associated cohesin variants on gene expression and cellular identity

Cohesin is a ring-shaped protein complex that controls dynamic chromosome structure. Cohesin activity is important for a variety of biological processes, including formation of DNA loops that regulate gene expression. The precise mechanisms by which cohesin shapes local chromosome structure and gene expression are not fully understood. Recurrent mutations in cohesin complex members have been reported in various cancers, though it is not clear whether many cohesin sequence variants have phenotypes and contribute to disease. Here, we utilized CRISPR/Cas9 genome editing to introduce a variety of cohesin sequence variants into murine embryonic stem cells and investigate their molecular and cellular consequences. Some of the cohesin variants tested caused changes to transcription, including altered expression of gene encoding lineage-specifying developmental regulators. Altered gene expression was also observed at insulated neighborhoods, where cohesin-mediated DNA loops constrain potential interactions between genes and enhancers. Furthermore, some cohesin variants altered the proliferation rate and differentiation potential of murine embryonic stem cells. This study provides a functional comparison of cohesin variants found in cancer within an isogenic system, revealing the relative roles of various cohesin perturbations on gene expression and maintenance of cellular identity.

Reversible chromatin condensation by the DNA repair and demethylation factor thymine DNA glycosylase

Chromatin structures (and modulators thereof) play a central role in genome organization and function. Herein, we report that thymine DNA glycosylase (TDG), an essential enzyme involved in DNA repair and demethylation, has the capacity to alter chromatin structure directly through its physical interactions with DNA. Using chemically defined nucleosome arrays, we demonstrate that TDG induces decompaction of individual chromatin fibers upon binding and promotes self-association of nucleosome arrays into higher-order oligomeric structures (i.e. condensation). Chromatin condensation is mediated by TDG’s disordered polycationic N-terminal domain, whereas its C-terminal domain antagonizes this process. Furthermore, we demonstrate that TDG-mediated chromatin condensation is reversible by growth arrest and DNA damage 45 alpha (GADD45a), implying that TDG cooperates with its binding partners to dynamically control chromatin architecture. Finally, we show that chromatin condensation by TDG is sensitive to the methylation status of the underlying DNA. This new paradigm for TDG has specific implications for associated processes, such as DNA repair, DNA demethylation, and transcription, and general implications for the role of DNA modification ‘readers’ in controlling chromatin organization.

Independence of chromatin conformation and gene regulation during Drosophila dorsoventral patterning

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

Histone H1 loss drives lymphoma by disrupting 3D chromatin architecture

Linker histone H1 proteins bind to nucleosomes and facilitate chromatin compaction¹, although their biological functions are poorly understood. Mutations in the genes that encode H1 isoforms B-E (H1B, H1C, H1D and H1E; also known as H1-5, H1-2, H1-3 and H1-4, respectively) are highly recurrent in B cell lymphomas, but the pathogenic relevance of these mutations to cancer and the mechanisms that are involved are unknown. Here we show that lymphoma-associated H1 alleles are genetic driver mutations in lymphomas. Disruption of H1 function results in a profound architectural remodelling of the genome, which is characterized by large-scale yet focal shifts of chromatin from a compacted to a relaxed state. This decompaction drives distinct changes in epigenetic states, primarily owing to a gain of histone H3 dimethylation at lysine 36 (H3K36me2) and/or loss of repressive H3 trimethylation at lysine 27 (H3K27me3). These changes unlock the expression of stem cell genes that are normally silenced during early development. In mice, loss of H1c and H1e (also known as H1f2 and H1f4, respectively) conferred germinal centre B cells with enhanced fitness and self-renewal properties, ultimately leading to aggressive lymphomas with an increased repopulating potential. Collectively, our data indicate that H1 proteins are normally required to sequester early developmental genes into architecturally inaccessible genomic compartments. We also establish H1 as a bona fide tumour suppressor and show that mutations in H1 drive malignant transformation primarily through three-dimensional genome reorganization, which leads to epigenetic reprogramming and derepression of developmentally silenced genes.

Hi-C implementation of genome structure for in silico models of radiation-induced DNA damage

Developments in the genome organisation field has resulted in the recent methodology to infer spatial conformations of the genome directly from experimentally measured genome contacts (Hi-C data). This provides a detailed description of both intra- and inter-chromosomal arrangements. Chromosomal intermingling is an important driver for radiation-induced DNA mis-repair. Which is a key biological endpoint of relevance to the fields of cancer therapy (radiotherapy), public health (biodosimetry) and space travel. For the first time, we leverage these methods of inferring genome organisation and couple them to nano-dosimetric radiation track structure modelling to predict quantities and distribution of DNA damage within cell-type specific geometries. These nano-dosimetric simulations are highly dependent on geometry and are benefited from the inclusion of experimentally driven chromosome conformations. We show how the changes in Hi-C contract maps impact the inferred geometries resulting in significant differences in chromosomal intermingling. We demonstrate how these differences propagate through to significant changes in the distribution of DNA damage throughout the cell nucleus, suggesting implications for DNA repair fidelity and subsequent cell fate. We suggest that differences in the geometric clustering for the chromosomes between the cell-types are a plausible factor leading to changes in cellular radiosensitivity. Furthermore, we investigate changes in cell shape, such as flattening, and show that this greatly impacts the distribution of DNA damage. This should be considered when comparing in vitro results to in vivo systems. The effect may be especially important when attempting to translate radiosensitivity measurements at the experimental in vitro level to the patient or human level.

Profiling chromatin regulatory landscape: insights into the development of ChIP-seq and ATAC-seq

Chromatin regulatory landscape plays a critical role in many disease processes and embryo development. Epigenome sequencing technologies such as chromatin immunoprecipitation sequencing (ChIP-seq) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) have enabled us to dissect the pan-genomic regulatory landscape of cells and tissues in both time and space dimensions by detecting specific chromatin state and its corresponding transcription factors. Pioneered by the advancement of chromatin immunoprecipitation-chip (ChIP-chip) technology, abundant epigenome profiling technologies have become available such as ChIP-seq, DNase I hypersensitive site sequencing (DNase-seq), ATAC-seq and so on. The advent of single-cell sequencing has revolutionized the next-generation sequencing, applications in single-cell epigenetics are enriched rapidly. Epigenome sequencing technologies have evolved from low-throughput to high-throughput and from bulk sample to the single-cell scope, which unprecedentedly benefits scientists to interpret life from different angles. In this review, after briefly introducing the background knowledge of epigenome biology, we discuss the development of epigenome sequencing technologies, especially ChIP-seq & ATAC-seq and their current applications in scientific research. Finally, we provide insights into future applications and challenges.

Determinants and role of chromatin organization in acute leukemia

DNA is compacted into higher order structures that have major implications in gene regulation. These structures allow for long-range interactions of DNA elements, such as the association of promoters with their cognate enhancers. In recent years, mutations in genes that control these structures, including the cohesin-complex and the insulator-binding protein CTCF, have been found in a spectrum of hematologic disorders, and especially in acute leukemias. Cohesin and CTCF are critical for mediating looping and establishing boundaries within chromatin. Cells that harbor mutations in these genes display aberrant chromatin architecture and resulting differences in gene expression that contribute to leukemia initiation and progression. Here, we provide detailed discussion of the nature of 3D interactions and the way that they are disrupted in acute leukemia. Continued research in this area will provide new insights into the mechanisms of leukemogenesis and may shed light on novel treatment strategies.

LINC00941 promotes oral squamous cell carcinoma progression via activating CAPRIN2 and canonical WNT/β-catenin signaling pathway

Dysregulation of long non-coding RNAs (lncRNAs) has been implicated in many cancer developments. Previous studies showed that lncRNA LINC00941 was aberrantly expressed in oral squamous cell carcinoma (OSCC). However, its role in OSCC development remains elusive. In this study, we demonstrated that in OSCC cells, EP300 activates LINC00941 transcription through up-regulating its promoter H3K27ac modification. Up-regulated LINC00941 in turn activates CAPRIN2 expression by looping to CAPRIN2 promoter. Functional assays suggest that both LINC00941 and CAPRIN2 play pivotal roles in promoting OSCC cell proliferation and colony formation. In vivo assay further confirmed the role of LINC00941 in promoting OSCC cell tumour formation. Lastly, we showed that the role of LINC00941 and CAPRIN2 in OSCC progression was mediated through activating the canonical WNT/β-catenin signaling pathway. Thus, LINC00941/CAPRIN2/ WNT/β-catenin signaling pathway provides new therapeutic targets for OSCC treatment.

Insulator-based loops mediate the spreading of H3K27me3 over distant micro-domains repressing euchromatin genes

Background

Chromosomes are subdivided spatially to delimit long-range interactions into topologically associating domains (TADs). TADs are often flanked by chromatin insulators and transcription units that may participate in such demarcation. Remarkably, single-cell Drosophila TAD units correspond to dynamic heterochromatin nano-compartments that can self-assemble. The influence of insulators on such dynamic compartmentalization remains unclear. Moreover, to what extent heterochromatin domains are fully compartmentalized away from active genes remains unclear from Drosophila to human.

Results

Here, we identify H3K27me3 micro-domains genome-wide in Drosophila, which are attributed to the three-dimensional spreading of heterochromatin marks into euchromatin. Whereas depletion of insulator proteins increases H3K27me3 spreading locally, across heterochromatin borders, it concomitantly decreases H3K27me3 levels at distant micro-domains discrete sites. Quantifying long-range interactions suggests that random interactions between heterochromatin TADs and neighbor euchromatin cannot predict the presence of micro-domains, arguing against the hypothesis that they reflect defects in self-folding or in insulating repressive TADs. Rather, micro-domains are predicted by specific long-range interactions with the TAD borders bound by insulator proteins and co-factors required for looping. Accordingly, H3K27me3 spreading to distant sites is impaired by insulator mutants that compromise recruitment of looping co-factors. Both depletions and insulator mutants significantly reduce H3K27me3 micro-domains, deregulating the flanking genes.

Conclusions

Our data highlight a new regulatory mode of H3K27me3 by insulator-based long-range interactions controlling distant euchromatic genes.

Chromatin Hyperacetylation Impacts Chromosome Folding by Forming a Nuclear Subcompartment

Delineating how chromosomes fold at length scales beyond one megabase remains obscure relative to smaller-scale folding into TADs, loops, and nucleosomes. We find that rather than simply unfolding chromatin, histone hyperacetylation results in interactions between distant genomic loci separated by tens to hundreds of megabases, even in the absence of transcription. These hyperacetylated "megadomains" are formed by the BRD4-NUT fusion oncoprotein, interact both within and between chromosomes, and form a specific nuclear subcompartment that has elevated gene activity with respect to other subcompartments. Pharmacological degradation of BRD4-NUT results in collapse of megadomains and attenuation of the interactions between them. In contrast, these interactions persist and contacts between newly acetylated regions are formed after inhibiting RNA polymerase II initiation. Our structure-function approach thus reveals that broad chromatin domains of identical biochemical composition, independent of transcription, form nuclear subcompartments, and also indicates the potential of altering chromosome structure for treating human disease.

Genome-Wide Identification and Analysis of Enhancer-Regulated microRNAs Across 31 Human Cancers

Enhancers are cis-regulatory DNA elements that positively regulate the transcription of target genes in a tissue-specific manner, and dysregulation of target genes could lead to various diseases, such as cancer. Recent studies have shown that enhancers can regulate microRNAs (miRNAs) and participate in their biological synthesis. However, the network of enhancer-regulated miRNAs across multiple cancers is still unclear. Here, a total of 2,418 proximal enhancer-miRNA interactions and 1,280 distal enhancer-miRNA interactions were identified through the integration of genomic distance, co-expression, and 3D genome data in 31 cancers. The results showed that both proximal and distal interactions exhibited a significant cancer type-specific feature trend at the tissue level rather than at the single-cell level, and there was a noteworthy positive correlation between the expression of the miRNA and the number of enhancers regulating the same miRNA in most cancers. Furthermore, we found that there was a high correlation between the formation of enhancer-miRNA pairs and the expression of enhancer RNAs (eRNAs) whether in distal or proximal regulation. The characteristics analysis showed that miRes (enhancers that regulated miRNAs) and non-miRes presented significant differences in sequence conservation, guanine-cytosine (GC) content, and histone modification signatures. Notably, GC content, H3K4me1, and H3K36me3 were present differently between distal and proximal regulation, suggesting that they might participate in chromosome looping of enhancer-miRNA interactions. Finally, we introduced a case study, enhancer: chr1:1186391-1186507 ∼ miR-200a was highly relevant to the survival of thyroid cancer patients and a cis-eQTL SNP on the enhancer affected the expression of the TNFRSF18 gene as a tumor suppressor.