Three-dimensional maps of all chromosomes in human male fibroblast nuclei and prometaphase rosettes

Major nuclear locales define nuclear genome organization and function beyond A and B compartments

Models of nuclear genome organization often propose a binary division into active versus inactive compartments yet typically overlook nuclear bodies. Here, we integrated analysis of sequencing and image-based data to compare genome organization in four human cell types relative to three different nuclear locales: the nuclear lamina, nuclear speckles, and nucleoli. Although gene expression correlates mostly with nuclear speckle proximity, DNA replication timing correlates with proximity to multiple nuclear locales. Speckle attachment regions emerge as DNA replication initiation zones whose replication timing and gene composition vary with their attachment frequency. Most facultative LADs retain a partially repressed state as iLADs, despite their positioning in the nuclear interior. Knock out of two lamina proteins, Lamin A and LBR, causes a shift of H3K9me3-enriched LADs from lamina to nucleolus, and a reciprocal relocation of H3K27me3-enriched partially repressed iLADs from nucleolus to lamina. Thus, these partially repressed iLADs appear to compete with LADs for nuclear lamina attachment with consequences for replication timing. The nuclear organization in adherent cells is polarized with nuclear bodies and genomic regions segregating both radially and relative to the equatorial plane. Together, our results underscore the importance of considering genome organization relative to nuclear locales for a more complete understanding of the spatial and functional organization of the human genome.

Organization and Dynamics of Chromosomes

How long thread-like eukaryotic chromosomes fit tidily in the small volume of the nucleus without significant entanglement is just beginning to be understood, thanks to major advances in experimental techniques. Several polymer models, which reproduce contact maps that measure the probabilities that two loci are in spatial contact, have predicted the 3D structures of interphase chromosomes. Data-driven approaches, using contact maps as input, predict that mitotic helical chromosomes are characterized by a switch in handedness, referred to as perversion. By using experimentally derived effective interactions between chromatin loci in simulations, structures of conventional and inverted nuclei have been accurately predicted. Polymer theory and simulations show that the dynamics of individual loci in chromatin exhibit subdiffusive behavior but the diffusion exponents are broadly distributed, which accords well with experiments. Although coarse-grained models are successful, many challenging problems remain, which require the creation of new experimental and computational tools to understand genome biology.

Beckwith-Wiedemann syndrome and large offspring syndrome involve alterations in methylome, transcriptome, and chromatin configuration

Beckwith-Wiedemann Syndrome (BWS) is the most common epigenetic overgrowth syndrome, caused by epigenetic alterations on chromosome 11p15. In ∼50% of patients with BWS, the imprinted region KvDMR1 (IC2) is hypomethylated. Nearly all children with BWS develop organ overgrowth and up to 28% develop cancer during childhood. The global epigenetic alterations beyond the 11p15 region in BWS are not currently known. Uncovering these alterations at the methylome, transcriptome, and chromatin architecture levels are necessary steps to improve the diagnosis and understanding of patients with BWS. Here we characterized the complete epigenetic profiles of BWS IC2 individuals together with the animal model of BWS, bovine large offspring syndrome (LOS). A novel finding of this research is the identification of two molecular subgroups of BWS IC2 individuals. Genome-wide alternations were detected for DNA methylation, transcript abundance, alternative splicing events of RNA, chromosome compartments, and topologically associating domains (TADs) in BWS and LOS, with shared alterations identified between species. Altered chromosome compartments and TADs were correlated with differentially expressed genes in BWS and LOS. Together, we highlight genes and genomic regions that have the potential to serve as targets for biomarker development to improve current molecular diagnostic methodologies for BWS.

BmNPV interacts with super-enhancer regions of the host chromatin to hijack cellular transcription machinery

Effective transcriptional activation relies on the spatial interaction between specific DNA elements. DNA interactions have also been observed between DNA viruses and their hosts, with limited understanding of the involved details. Baculovirus is a representative species of DNA virus and has been reported to interact with the host genome in our previous study. However, the biological significance of the baculovirus-host trans-species DNA interaction and its underlying mechanisms remain elusive. Here, using Bombyx mori nucleopolyhedrovirus (BmNPV) as the model virus, we combine epigenome, transcriptome, and biochemical assays to investigate the baculovirus-host DNA interaction. Our data show that BmNPV hijacks the transcriptional regulatory capacity of host super-enhancers (SEs) by physically interacting with these regions on the host genome. This results in the usurpation of the activating capacity of an SE-binding transcription factor GATA by the virus, thereby impairing the SE-induced specific transcriptional activation of the target antiviral genes. Moreover, the hijacked regulatory capacity is spread on BmNPV genome through cis-interaction of viral DNA, leading to enhanced viral gene expression. Overall, our results provide novel insights into the intricate interplay of viruses with host gene expression regulatory networks and broaden the vision in the mechanisms of viral exploitation on cellular machinery.

Nanoscale DNA tracing reveals the self-organization mechanism of mitotic chromosomes

How genomic DNA is folded during cell division to form the characteristic rod-shaped mitotic chromosomes essential for faithful genome inheritance is a long-standing open question in biology. Here, we use nanoscale DNA tracing in single dividing cells to directly visualize how the 3D fold of genomic DNA changes during mitosis at scales from single loops to entire chromosomes. Our structural analysis reveals a characteristic genome scaling minimum of 6-8 megabases in mitosis. Combined with data-driven modeling and molecular perturbations, we can show that very large and strongly overlapping loops formed by condensins are the fundamental structuring principle of mitotic chromosomes. These loops compact chromosomes locally and globally to the limit set by chromatin self-repulsion. The characteristic length, density, and increasingly overlapping structure of mitotic loops we observe in 3D fully explain how the rod-shaped mitotic chromosome structure emerges by self-organization during cell division.

Visualization of chromosomal reorganization induced by heterologous fusions in the mammalian nucleus

Chromosomes are spatially organized and functionally folded into a specific macro-structure in the nucleus. Recently, we and others created haploid cells with chromosome fusions. However, there is still lack of an effective strategy for precisely investigating how the genome copes with fusions. Here, we developed a down-sampling method to convert the populational Hi-C dataset into single cell-like Khimaira Matrix (K-matrix). K-matrix preserves not only the most prominent functional genomic features but also cell-to-cell variations. K-matrix-originated genome 3D models display spatial approach of fused chromosomes and minor global structure alterations. Combined with a layered positional decomposition analysis, our models indicate slight re-adjustment of chromosome distributions accordingly with an increasing tendency following more fusions involved. Nevertheless, the radial distribution of the A/B compartment is not affected dramatically. By contrast, natural populations harboring Rb fusions display significant alterations of chromosome radial location. Overall, K-matrix-originated models enable visualization of chromosomal reorganization with high resolution.

Long-range genomic loci stochastically assemble into combinatorial forms of chromosome skeleton

One fundamental yet open question is how eukaryotic chromosomes fold into segregated territories, a process essential for gene transcription and cell fate. Through analyzing Hi-C and chromatin-tracing DNA-FISH data, we identify long-range chromo skeleton loop structures that span over 100 Mb, extending beyond the reach of loop extrusion models. Spatial density analyses point to assembly formation independent of major nuclear structures. A subset of genomic loci serves as nucleation centers, driving loop clustering. These complexes are highly stable, as shown by live-cell imaging with sequence-specific fluorescent labeling, and biophysical model analyses reveal a multivalent binding mechanism. Our findings suggest a redundant, distributed cluster mechanism that ensures robustness across cell types and against mutations, guiding both chromosome compaction and the formation of smaller-scale chromosomal structures.

Ipsilateral restriction of chromosome movement along a centrosome, and apical-basal axis during the cell cycle

Little is known about how distance between homologous chromosomes are controlled during the cell cycle. Here, we show that the distribution of centromere components display two discrete clusters placed to either side of the centrosome and apical/basal axis from prophase to G1 interphase. 4-Dimensional live cell imaging analysis of centromere and centrosome tracking reveals that centromeres oscillate largely within one cluster, but do not cross over to the other cluster. We propose a model of an axis-dependent ipsilateral restriction of chromosome oscillations throughout mitosis.

Chromosome architecture and low cohesion bias acrocentric chromosomes towards aneuploidy during mammalian meiosis

Aneuploidy in eggs is a leading cause of miscarriages or viable developmental syndromes. Aneuploidy rates differ between individual chromosomes. For instance, chromosome 21 frequently missegregates, resulting in Down Syndrome. What causes chromosome-specific aneuploidy in meiosis is unclear. Chromosome 21 belongs to the class of acrocentric chromosomes, whose centromeres are located close to the chromosome end, resulting in one long and one short chromosome arm. We demonstrate that acrocentric chromosomes are generally more often aneuploid than metacentric chromosomes in porcine eggs. Kinetochores of acrocentric chromosomes are often partially covered by the short chromosome arm during meiosis I in human and porcine oocytes and orient less efficiently toward the spindle poles. These partially covered kinetochores are more likely to be incorrectly attached to the spindle. Additionally, sister chromatids of acrocentric chromosomes are held together by lower levels of cohesin, making them more vulnerable to age-dependent cohesin loss. Chromosome architecture and low cohesion therefore bias acrocentric chromosomes toward aneuploidy during mammalian meiosis.

Cytogenetic bands and sharp peaks of Alu underlie large-scale segmental regulation of nuclear genome architecture

Cytogenetic bands reflect genomic organization in large blocks of DNA with similar properties. Because banding patterns are invariant, this organization may often be assumed unimportant for genome regulation. Results here challenge that view. Findings here suggest cytogenetic bands reflect a visible framework upon which regulated genome architecture is built. Given Alu and L1 densities differ in cytogenetic bands, we examined their distribution after X-chromosome inactivation or formation of senescent-associated heterochromatin foci (SAHFs). Alu-rich regions remain outside both SAHFs and the Barr Body (BB), affirming that the BB is not the whole chromosome but a condensed, L1-rich core. Hi-C analysis of senescent cells demonstrates large (~10 Mb) G-bands remodel as a contiguous unit, gaining distal intrachromosomal interactions as syntenic G-bands coalesce into SAHFs. Striking peaks of Alu within R-bands strongly resist condensation. Thus, large-scale segmental genome architectur relates to dark versus light cytogenetic bands and Alu-peaks, implicating both in chromatin regulation.

Nanoscale 3D DNA tracing reveals the mechanism of self-organization of mitotic chromosomes

How genomic DNA is folded during cell division to form the characteristic rod-shaped mitotic chromosomes essential for faithful genome inheritance is a long-standing open question in biology. Here, we use nanoscale DNA-tracing in single dividing cells to directly visualize how the 3D fold of genomic DNA changes during mitosis, at scales from single loops to entire chromosomes. Our structural analysis reveals a characteristic genome scaling minimum at 6-8 Mbp in mitosis. Combined with data-driven modeling and molecular perturbations, we can show that very large and strongly overlapping loops formed by Condensins are the fundamental structuring principle of mitotic chromosomes. These loops compact chromosomes locally and globally to the limit set by chromatin self-repulsion. The characteristic length, density and increasingly overlapping structure of mitotic loops we observe in 3D, fully explain how the rod-shaped mitotic chromosome structure emerges by self-organization during cell division.

Phase Separation Mediated Sub-Nuclear Compartmentalization of Androgen Receptors

The androgen receptor (AR), a member of the nuclear steroid hormone receptor family of transcription factors, plays a crucial role not only in the development of the male phenotype but also in the development and growth of prostate cancer. While AR structure and AR interactions with coregulators and chromatin have been studied in detail, improving our understanding of AR function in gene transcription regulation, the spatio-temporal organization and the role of microscopically discernible AR foci in the nucleus are still underexplored. This review delves into the molecular mechanisms underlying AR foci formation, focusing on liquid-liquid phase separation and its role in spatially organizing ARs and their binding partners within the nucleus at transcription sites, as well as the influence of 3D-genome organization on AR-mediated gene transcription.

Spandrels of the cell nucleus

S.J. Gould and R. Lewontin in their famous "Spandrels paper" (1979) argued that many anatomical elements arise in evolution not due to their "current utility" but rather due to other "reasons for origin", such as other developmental processes, physical constraints and mechanical forces. Here, in the same spirit, we argue that a variety of molecular processes, physical constraints, and mechanical forces, alone or together, generate structures that are detectable in the cell nucleus, yet these structures themselves may not carry any specific function, being a mere reflection of processes that produced them.

Nucleolus and centromere Tyramide Signal Amplification-Seq reveals variable localization of heterochromatin in different cell types

Genome differential positioning within interphase nuclei remains poorly explored. We extended and validated Tyramide Signal Amplification (TSA)-seq to map genomic regions near nucleoli and pericentric heterochromatin in four human cell lines. Our study confirmed that smaller chromosomes localize closer to nucleoli but further deconvolved this by revealing a preference for chromosome arms below 36-46 Mbp in length. We identified two lamina associated domain subsets through their differential nuclear lamina versus nucleolar positioning in different cell lines which showed distinctive patterns of DNA replication timing and gene expression across all cell lines. Unexpectedly, active, nuclear speckle-associated genomic regions were found near typically repressive nuclear compartments, which is attributable to the close proximity of nuclear speckles and nucleoli in some cell types, and association of centromeres with nuclear speckles in human embryonic stem cells (hESCs). Our study points to a more complex and variable nuclear genome organization than suggested by current models, as revealed by our TSA-seq methodology.

Hi-C-guided many-polymer model to decipher 3D genome organization

We propose a high-throughput chromosome conformation capture data-based many-polymer model that allows us to generate an ensemble of multi-scale genome structures. We demonstrate the efficacy of our model by validating the generated structures against experimental measurements and employ them to address key questions regarding genome organization. Our model first confirms a significant correlation between chromosome size and nuclear positioning. Specifically, smaller chromosomes are distributed at the core region, whereas larger chromosomes are at the periphery, interacting with the nuclear envelope. The spatial distribution of A- and B-type compartments, which is nontrivial to infer from the corresponding high-throughput chromosome conformation capture maps alone, can also be elucidated using our model, accounting for an issue such as the effect of chromatin-lamina interaction on the compartmentalization of conventional and inverted nuclei. In accordance with imaging data, the overall shape of the 3D genome structures generated from our model displays significant variation. As a case study, we apply our method to the yellow fever mosquito genome, finding that the predicted morphology displays, on average, a more globular shape than the previously suggested spindle-like organization and that our prediction better aligns with the fluorescence in situ hybridization data. Our model has great potential to be extended to investigate many outstanding issues concerning 3D genome organization.

Generation of densely labeled oligonucleotides for the detection of small genomic elements

The genome contains numerous regulatory elements that may undergo complex interactions and contribute to the establishment, maintenance, and change of cellular identity. Three-dimensional genome organization can be explored with fluorescence in situ hybridization (FISH) at the single-cell level, but the detection of small genomic loci remains challenging. Here, we provide a rapid and simple protocol for the generation of bright FISH probes suited for the detection of small genomic elements. We systematically optimized probe design and synthesis, screened polymerases for their ability to incorporate dye-labeled nucleotides, and streamlined purification conditions to yield nanoscopy-compatible oligonucleotides with dyes in variable arrays (NOVA probes). With these probes, we detect genomic loci ranging from genome-wide repetitive regions down to non-repetitive loci below the kilobase scale. In conclusion, we introduce a simple workflow to generate densely labeled oligonucleotide pools that facilitate detection and nanoscopic measurements of small genomic elements in single cells.

Modelling DNA damage-repair and beyond

The paper presents a review of mechanistic modelling studies of DNA damage and DNA repair, and consequences to follow in mammalian cell nucleus. We hypothesize DNA deletions are consequences of repair of double strand breaks leading to the modifications of genome that play crucial role in long term development of genetic inheritance and diseases. The aim of the paper is to review formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double strand breaks and deletions in damaged human genome from endogenous and exogenous events. The model of the cell nucleus presented enables simulation of DNA damage at molecular level identifying the spectrum of damage induced in all chromosomal territories and loops. Our mechanistic modelling of DNA repair for double stand breaks (DSB), single strand breaks (SSB) and base damage (BD), shows the complexity of DNA damage is responsible for the longer repair times and the reason for the biphasic feature of mammalian cells repair curves. In the absence of experimentally determined data, the mechanistic model of repair predicts the in vivo rate constants for the proteins involved in the repair of DSB, SSB, and of BD.

A critical role for X-chromosome architecture in mammalian X-chromosome dosage compensation

To regulate gene expression, the macromolecular components of the mammalian interphase nucleus are spatially organized into a myriad of functional compartments. Over the past decade, increasingly sophisticated genomics, microscopy, and functional approaches have probed this organization in unprecedented detail. These investigations have linked chromatin-associated noncoding RNAs to specific nuclear compartments and uncovered mechanisms by which these RNAs establish such domains. In this review, we focus on the long non-coding RNA Xist and summarize new evidence demonstrating the significance of chromatin reconfiguration in creating the inactive X-chromosome compartment. Differences in chromatin compaction correlate with distinct levels of gene repression on the X-chromosome, potentially explaining how human XIST can induce chromosome-wide dampening and silencing of gene expression at different stages of human development.

The Central Role of Cytogenetics in Radiation Biology

Radiation cytogenetics has a rich history seldom appreciated by those outside the field. Early radiobiology was dominated by physics and biophysical concepts that borrowed heavily from the study of radiation-induced chromosome aberrations. From such studies, quantitative relationships between biological effect and changes in absorbed dose, dose rate and ionization density were codified into key concepts of radiobiological theory that have persisted for nearly a century. This review aims to provide a historical perspective of some of these concepts, including evidence supporting the contention that chromosome aberrations underlie development of many, if not most, of the biological effects of concern for humans exposed to ionizing radiations including cancer induction, on the one hand, and tumor eradication on the other. The significance of discoveries originating from these studies has widened and extended far beyond their original scope. Chromosome structural rearrangements viewed in mitotic cells were first attributed to the production of breaks by the radiations during interphase, followed by the rejoining or mis-rejoining among ends of other nearby breaks. These relatively modest beginnings eventually led to the discovery and characterization of DNA repair of double-strand breaks by non-homologous end joining, whose importance to various biological processes is now widely appreciated. Two examples, among many, are V(D)J recombination and speciation. Rapid technological advancements in cytogenetics, the burgeoning fields of molecular radiobiology and third-generation sequencing served as a point of confluence between the old and new. As a result, the emergent field of "cytogenomics" now becomes uniquely positioned for the purpose of more fully understanding mechanisms underlying the biological effects of ionizing radiation exposure.

Live chromosome identifying and tracking reveals size-based spatial pathway of meiotic errors in oocytes

Meiotic errors of relatively small chromosomes in oocytes result in egg aneuploidies that cause miscarriages and congenital diseases. Unlike somatic cells, which preferentially mis-segregate larger chromosomes, aged oocytes preferentially mis-segregate smaller chromosomes through unclear processes. Here, we provide a comprehensive three-dimensional chromosome identifying-and-tracking dataset throughout meiosis I in live mouse oocytes. This analysis reveals a prometaphase pathway that actively moves smaller chromosomes to the inner region of the metaphase plate. In the inner region, chromosomes are pulled by stronger bipolar microtubule forces, which facilitates premature chromosome separation, a major cause of segregation errors in aged oocytes. This study reveals a spatial pathway that facilitates aneuploidy of small chromosomes preferentially in aged eggs and implicates the role of the M phase in creating a chromosome size-based spatial arrangement.

Genetic/epigenetic effects in NF1 microdeletion syndrome: beyond the haploinsufficiency, looking at the contribution of not deleted genes

NF1 microdeletion syndrome, accounting for 5-11% of NF1 patients, is caused by a deletion in the NF1 region and it is generally characterized by a severe phenotype. Although 70% of NF1 microdeletion patients presents the same 1.4 Mb type-I deletion, some patients may show additional clinical features. Therefore, the contribution of several pathogenic mechanisms, besides haploinsufficiency of some genes within the deletion interval, is expected and needs to be defined. We investigated an altered expression of deletion flanking genes by qPCR in patients with type-1 NF1 deletion, compared to healthy donors, possibly contributing to the clinical traits of NF1 microdeletion syndrome. In addition, the 1.4-Mb deletion leads to changes in the 3D chromatin structure in the 17q11.2 region. Specifically, this deletion alters DNA-DNA interactions in the regions flanking the breakpoints, as demonstrated by our 4C-seq analysis. This alteration likely causes position effect on the expression of deletion flanking genes.Interestingly, 4C-seq analysis revealed that in microdeletion patients, an interaction was established between the RHOT1 promoter and the SLC6A4 gene, which showed increased expression. We performed NGS on putative modifier genes, and identified two "likely pathogenic" rare variants in RAS pathway, possibly contributing to incidental phenotypic features.This study provides new insights into understanding the pathogenesis of NF1 microdeletion syndrome and suggests a novel pathomechanism that contributes to the expression phenotype in addition to haploinsufficiency of genes located within the deletion.This is a pivotal approach that can be applied to unravel microdeletion syndromes, improving precision medicine, prognosis and patients' follow-up.

Chromatin image-driven modelling

The challenge of modelling the spatial conformation of chromatin remains an open problem. While multiple data-driven approaches have been proposed, each has limitations. This work introduces two image-driven modelling methods based on the Molecular Dynamics Flexible Fitting (MDFF) approach: the force method and the correlational method. Both methods have already been used successfully in protein modelling. We propose a novel way to employ them for building chromatin models directly from 3D images. This approach is termed image-driven modelling. Additionally, we introduce the initial structure generator, a tool designed to generate optimal starting structures for the proposed algorithms. The methods are versatile and can be applied to various data types, with minor modifications to accommodate new generation imaging techniques.

Potential roles of inter-chromosomal interactions in cell fate determination

Mammalian genomic DNA is packed in a small nucleus, and its folding and organization in the nucleus are critical for gene regulation and cell fate determination. In interphase, chromosomes are compartmentalized into certain nuclear spaces and territories that are considered incompatible with each other. The regulation of gene expression is influenced by the epigenetic characteristics of topologically associated domains and A/B compartments within chromosomes (intrachromosomal). Previously, interactions among chromosomes detected via chromosome conformation capture-based methods were considered noise or artificial errors. However, recent studies based on newly developed ligation-independent methods have shown that inter-chromosomal interactions play important roles in gene regulation. This review summarizes the recent understanding of spatial genomic organization in mammalian interphase nuclei and discusses the potential mechanisms that determine cell identity. In addition, this review highlights the potential role of inter-chromosomal interactions in early mouse development.

Ipsilateral restriction of chromosome movement along a centrosome, and apical-basal axis during the cell cycle

Little is known about how distance between homologous chromosomes are controlled during the cell cycle. Here, we show that the distribution of centromere components display two discrete clusters placed to either side of the centrosome and apical/basal axis from prophase to G1 interphase. 4-Dimensional live cell imaging analysis of centromere and centrosome tracking reveals that centromeres oscillate largely within one cluster, but do not cross over to the other cluster. We propose a model of an axis-dependent ipsilateral restriction of chromosome oscillations throughout mitosis.

Non-random spatial organization of telomeres varies during the cell cycle and requires LAP2 and BAF

Spatial genome organization within the nucleus influences major biological processes and is impacted by the configuration of linear chromosomes. Here, we applied 3D spatial statistics and modeling on high-resolution telomere and centromere 3D-structured illumination microscopy images in cancer cells. We found a multi-scale organization of telomeres that dynamically evolved from a mixed clustered-and-regular distribution in early G1 to a purely regular distribution as cells progressed through the cell cycle. In parallel, our analysis revealed two pools of peripheral and internal telomeres, the proportions of which were inverted during the cell cycle. We then conducted a targeted screen using MadID to identify the molecular pathways driving or maintaining telomere anchoring to the nuclear envelope observed in early G1. Lamina-associated polypeptide (LAP) proteins were found transiently localized to telomeres in anaphase, a stage where LAP2α initiates the reformation of the nuclear envelope, and impacted telomere redistribution in the next interphase together with their partner barrier-to-autointegration factor (BAF).

The scales, mechanisms, and dynamics of the genome architecture

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

Ipsilateral restriction of chromosome movement along a centrosome, and apical-basal axis during the cell cycle

Little is known about how distance between homologous chromosomes are controlled during the cell cycle. Here, we show that the distribution of centromere components display two discrete clusters placed to either side of the centrosome and apical/basal axis from prophase to G 1 interphase. 4-Dimensional live cell imaging analysis of centromere and centrosome tracking reveals that centromeres oscillate largely within one cluster, but do not cross over to the other cluster. We propose a model of an axis-dependent ipsilateral restriction of chromosome oscillations throughout mitosis.

Nuclear morphology is shaped by loop-extrusion programs

It is well established that neutrophils adopt malleable polymorphonuclear shapes to migrate through narrow interstitial tissue spaces1-3. However, how polymorphonuclear structures are assembled remains unknown4. Here we show that in neutrophil progenitors, halting loop extrusion-a motor-powered process that generates DNA loops by pulling in chromatin5-leads to the assembly of polymorphonuclear genomes. Specifically, we found that in mononuclear neutrophil progenitors, acute depletion of the loop-extrusion loading factor nipped-B-like protein (NIPBL) induced the assembly of horseshoe, banded, ringed and hypersegmented nuclear structures and led to a reduction in nuclear volume, mirroring what is observed during the differentiation of neutrophils. Depletion of NIPBL also induced cell-cycle arrest, activated a neutrophil-specific gene program and conditioned a loss of interactions across topologically associating domains to generate a chromatin architecture that resembled that of differentiated neutrophils. Removing NIPBL resulted in enrichment for mega-loops and interchromosomal hubs that contain genes associated with neutrophil-specific enhancer repertoires and an inflammatory gene program. On the basis of these observations, we propose that in neutrophil progenitors, loop-extrusion programs produce lineage-specific chromatin architectures that permit the packing of chromosomes into geometrically confined lobular structures. Our data also provide a blueprint for the assembly of polymorphonuclear structures, and point to the possibility of engineering de novo nuclear shapes to facilitate the migration of effector cells in densely populated tumorigenic environments.

Boundary stacking interactions enable cross-TAD enhancer-promoter communication during limb development

Although promoters and their enhancers are frequently contained within a topologically associating domain (TAD), some developmentally important genes have their promoter and enhancers within different TADs. Hypotheses about molecular mechanisms enabling cross-TAD interactions remain to be assessed. To test these hypotheses, we used optical reconstruction of chromatin architecture to characterize the conformations of the Pitx1 locus on single chromosomes in developing mouse limbs. Our data support a model in which neighboring boundaries are stacked as a result of loop extrusion, bringing boundary-proximal cis-elements into contact. This stacking interaction also contributes to the appearance of architectural stripes in the population average maps. Through molecular dynamics simulations, we found that increasing boundary strengths facilitates the formation of the stacked boundary conformation, counter-intuitively facilitating border bypass. This work provides a revised view of the TAD borders' function, both facilitating and preventing cis-regulatory interactions, and introduces a framework to distinguish border-crossing from border-respecting enhancer-promoter pairs.

Epigenetic regulatory layers in the 3D nucleus

Nearly 7 decades have elapsed since Francis Crick introduced the central dogma of molecular biology, as part of his ideas on protein synthesis, setting the fundamental rules of sequence information transfer from DNA to RNAs and proteins. We have since learned that gene expression is finely tuned in time and space, due to the activities of RNAs and proteins on regulatory DNA elements, and through cell-type-specific three-dimensional conformations of the genome. Here, we review major advances in genome biology and discuss a set of ideas on gene regulation and highlight how various biomolecular assemblies lead to the formation of structural and regulatory features within the nucleus, with roles in transcriptional control. We conclude by suggesting further developments that will help capture the complex, dynamic, and often spatially restricted events that govern gene expression in mammalian cells.

Differences in Alu vs L1-rich chromosome bands underpin architectural reorganization of the inactive-X chromosome and SAHFs

The linear DNA sequence of mammalian chromosomes is organized in large blocks of DNA with similar sequence properties, producing a pattern of dark and light staining bands on mitotic chromosomes. Cytogenetic banding is essentially invariant between people and cell-types and thus may be assumed unrelated to genome regulation. We investigate whether large blocks of Alu-rich R-bands and L1-rich G-bands provide a framework upon which functional genome architecture is built. We examine two models of large-scale chromatin condensation: X-chromosome inactivation and formation of senescence-associated heterochromatin foci (SAHFs). XIST RNA triggers gene silencing but also formation of the condensed Barr Body (BB), thought to reflect cumulative gene silencing. However, we find Alu-rich regions are depleted from the L1-rich BB, supporting it is a dense core but not the entire chromosome. Alu-rich bands are also gene-rich, affirming our earlier findings that genes localize at the outer periphery of the BB. SAHFs similarly form within each territory by coalescence of syntenic L1 regions depleted for highly Alu-rich DNA. Analysis of senescent cell Hi-C data also shows large contiguous blocks of G-band and R-band DNA remodel as a segmental unit. Entire dark-bands gain distal intrachromosomal interactions as L1-rich regions form the SAHF. Most striking is that sharp Alu peaks within R-bands resist these changes in condensation. We further show that Chr19, which is exceptionally Alu rich, fails to form a SAHF. Collective results show regulation of genome architecture corresponding to large blocks of DNA and demonstrate resistance of segments with high Alu to chromosome condensation.

Profiling Chromosome Topological Features by Super-Resolution 3D Structured Illumination Microscopy

Three-dimensional structured illumination microscopy (3D-SIM) and fluorescence in situ hybridization on three-dimensional preserved cells (3D-FISH) have proven to be robust and efficient methodologies for analyzing nuclear architecture and profiling the genome's topological features. These methods have allowed the simultaneous visualization and evaluation of several target structures at super-resolution. In this chapter, we focus on the application of 3D-SIM for the visualization of 3D-FISH preparations of chromosomes in interphase, known as Chromosome Territories (CTs). We provide a workflow and detailed guidelines for sample preparation, image acquisition, and image analysis to obtain quantitative measurements for profiling chromosome topological features. In parallel, we address a practical example of these protocols in the profiling of CTs 9 and 22 involved in the translocation t(9;22) in Chronic Myeloid Leukemia (CML). The profiling of chromosome topological features described in this chapter allowed us to characterize a large-scale topological disruption of CTs 9 and 22 that correlates directly with patients' response to treatment and as a possible potential change in the inheritance systems. These findings open new insights into how the genome structure is associated with the response to cancer treatments, highlighting the importance of microscopy in analyzing the topological features of the genome.

Reconstructing diploid 3D chromatin structures from single cell Hi-C data with a polymer-based approach

Detailed understanding of the 3D structure of chromatin is a key ingredient to investigate a variety of processes inside the cell. Since direct methods to experimentally ascertain these structures lack the desired spatial fidelity, computational inference methods based on single cell Hi-C data have gained significant interest. Here, we develop a progressive simulation protocol to iteratively improve the resolution of predicted interphase structures by maximum-likelihood association of ambiguous Hi-C contacts using lower-resolution predictions. Compared to state-of-the-art methods, our procedure is not limited to haploid cell data and allows us to reach a resolution of up to 5,000 base pairs per bead. High resolution chromatin models grant access to a multitude of structural phenomena. Exemplarily, we verify the formation of chromosome territories and holes near aggregated chromocenters as well as the inversion of the CpG content for rod photoreceptor cells.

Adhesome Receptor Clustering is Accompanied by the Colocalization of the Associated Genes in the Cell Nucleus

Proteins on the cell membrane cluster to respond to extracellular signals; for example, adhesion proteins cluster to enhance extracellular matrix sensing; or T-cell receptors cluster to enhance antigen sensing. Importantly, the maturation of such receptor clusters requires transcriptional control to adapt and reinforce the extracellular signal sensing. However, it has been unclear how such efficient clustering mechanisms are encoded at the level of the genes that code for these receptor proteins. Using the adhesome as an example, we show that genes that code for adhesome receptor proteins are spatially co-localized and co-regulated within the cell nucleus. Towards this, we use Hi-C maps combined with RNA-seq data of adherent cells to map the correspondence between adhesome receptor proteins and their associated genes. Interestingly, we find that the transcription factors that regulate these genes are also co-localized with the adhesome gene loci, thereby potentially facilitating a transcriptional reinforcement of the extracellular matrix sensing machinery. Collectively, our results highlight an important layer of transcriptional control of cellular signal sensing.

Subgraph extraction and graph representation learning for single cell Hi-C imputation and clustering

Single-cell Hi-C (scHi-C) technology enables the investigation of 3D chromatin structure variability across individual cells. However, the analysis of scHi-C data is challenged by a large number of missing values. Here, we present a scHi-C data imputation model HiC-SGL, based on Subgraph extraction and graph representation learning. HiC-SGL can also learn informative low-dimensional embeddings of cells. We demonstrate that our method surpasses existing methods in terms of imputation accuracy and clustering performance by various metrics.

Corticosteroid-induced chromatin loop dynamics at the FKBP5 gene

FKBP5 is a 115-kb-long glucocorticoid-inducible gene implicated in psychiatric disorders. To investigate the complexities of chromatin interaction frequencies at the FKBP5 topologically associated domain (TAD), we deployed 15 one-to-all chromatin capture viewpoints near gene promoters, enhancers, introns, and CTCF-loop anchors. This revealed a "one-TAD-one-gene" structure encompassing the FKBP5 promoter and its enhancers. The FKBP5 promoter and its two glucocorticoid-stimulated enhancers roam the entire TAD while displaying subtle cell type-specific interactomes. The FKBP5 TAD consists of two nested CTCF loops that are coordinated by one CTCF site in the eighth intron of FKBP5 and another beyond its polyadenylation site, 61 kb further. Loop extension correlates with transcription increases through the intronic CTCF site. This is efficiently compensated for, since the short loop is restored even under high transcription regimes. The boundaries of the FKBP5 TAD consist of divergent CTCF site patterns, harbor multiple smaller genes, and are resilient to glucocorticoid stimulation. Interestingly, both FKBP5 TAD boundaries harbor H3K27me3-marked heterochromatin blocks that may reinforce them. We propose that cis-acting genetic and epigenetic polymorphisms underlying FKBP5 expression variation are likely to reside within a 240-kb region that consists of the FKBP5 TAD, its left sub-TAD, and both its boundaries.

Molecular Mechanisms of Alu and LINE-1 Interspersed Repetitive Sequences Reveal Diseases of Visual System Dysfunction

BACKGROUND

Short interspersed nuclear elements (SINEs) and long interspersed nuclear elements (LINE-1s) are the abundant and well-characterized repetitive elements in the human genome.

METHODS

For this review, all relevant original research studies were assessed by searching electronic databases, including PubMed, Google Scholar, and Web of Science, by using relevant keywords. Accumulating evidence indicates that the disorder of gene expression regulated by these repetitive sequences is one of the causes of the diseases of visual system dysfunction, including retinal degenerations, glaucoma, retinitis punctata albescens, retinitis pigmentosa, geographic atrophy, and age-related macular degeneration, suggesting that SINEs and LINE-1s may have great potential implications in ophthalmology.

RESULTS

Alu elements belonging to the SINEs are present in more than one million copies, comprising 10% of the human genome.

CONCLUSION

This study offers recent advances in Alu and LINE-1 mechanisms in the development of eye diseases. The current study could advance our knowledge of the roles of SINEs and LINE-1s in the developing process of eye diseases, suggesting new diagnostic biomarkers, therapeutic strategies, and significant points for future studies.

Nucleolus and centromere TSA-Seq reveals variable localization of heterochromatin in different cell types

Genome differential positioning within interphase nuclei remains poorly explored. We extended and validated TSA-seq to map genomic regions near nucleoli and pericentric heterochromatin in four human cell lines. Our study confirmed that smaller chromosomes localize closer to nucleoli but further deconvolved this by revealing a preference for chromosome arms below 36-46 Mbp in length. We identified two lamina associated domain subsets through their differential nuclear lamina versus nucleolar positioning in different cell lines which showed distinctive patterns of DNA replication timing and gene expression across all cell lines. Unexpectedly, active, nuclear speckle-associated genomic regions were found near typically repressive nuclear compartments, which is attributable to the close proximity of nuclear speckles and nucleoli in some cell types, and association of centromeres with nuclear speckles in hESCs. Our study points to a more complex and variable nuclear genome organization than suggested by current models, as revealed by our TSA-seq methodology.

The Landscape of MYB/MYBL1- and Peri-MYB/MYBL1-Associated Rearrangements in Adenoid Cystic Carcinoma

Approximately 60% of adenoid cystic carcinoma (AdCC) cases are positive for MYB::NFIB or MYBL1::NFIB, whereas MYB/MYBL1 oncoprotein, a key driver of AdCC, is overexpressed in most cases. Juxtaposition of superenhancer regions in NFIB and other genes into the MYB/MYBL1 locus is an attractive oncogenic hypothesis for AdCC cases, either negative or positive for MYB/MYBL1::NFIB. However, evidence supporting this hypothesis is insufficient. We examined 160 salivary AdCC cases for rearrangements in MYB/MYBL1 loci and peri-MYB/MYBL1 areas (centromeric and telomeric areas of 10 Mb each) using formalin-fixed, paraffin-embedded tumor sections. For the detection of the rearrangements, we employed conventional fluorescence in situ hybridization split and fusion assays and a 5 Mb fluorescence in situ hybridization split assay. The latter is a novel assay that enabled us to detect any possible splits within a 5 Mb distance of a chromosome. We found MYB/MYBL1- and peri-MYB/MYBL1-associated rearrangements in 149/160 patients (93%). AdCC cases positive for rearrangements in MYB, MYBL1, the peri-MYB area, and the peri-MYBL1 area numbered 105 (66%), 20 (13%), 19 (12%), and 5 (3%), respectively. In 24 peri-MYB/MYBL1 rearrangement-positive cases, 14 (58%) were found to have a juxtaposition of the NFIB or RAD51B locus into the MYB/MYBL1 loci. On comparing with a tumor group positive for MYB::NFIB, a hallmark of AdCC, other genetically classified tumor groups had similar features of overexpression of the MYB transcript and MYB oncoprotein as detected by semiquantitative RT-qPCR and immunohistochemistry, respectively. In addition, clinicopathological and prognostic features were similar among these groups. Our study suggests that peri-MYB/MYBL1 rearrangements may be a frequent event in AdCC and may result in biological and clinicopathological consequences comparable to MYB/MYBL1 rearrangements. The landscape of MYB/MYBL1 and peri-MYB/MYBL1 rearrangements shown here strongly suggests that juxtaposition of superenhancers into MYB/MYBL1 or peri-MYB/MYBL1 loci is an alteration that acts as a key driver for AdCC oncogenesis and may unify MYB/MYBL1 rearrangement-positive and negative cases.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

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.

Three-dimensional genome structure and function

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

Region-based epigenetic clock design improves RRBS-based age prediction

Recent studies suggest that epigenetic rejuvenation can be achieved using drugs that mimic calorie restriction and techniques such as reprogramming-induced rejuvenation. To effectively test rejuvenation in vivo, mouse models are the safest alternative. However, we have found that the recent epigenetic clocks developed for mouse reduced-representation bisulphite sequencing (RRBS) data have significantly poor performance when applied to external datasets. We show that the sites captured and the coverage of key CpGs required for age prediction vary greatly between datasets, which likely contributes to the lack of transferability in RRBS clocks. To mitigate these coverage issues in RRBS-based age prediction, we present two novel design strategies that use average methylation over large regions rather than individual CpGs, whereby regions are defined by sliding windows (e.g. 5 kb), or density-based clustering of CpGs. We observe improved correlation and error in our regional blood clocks (RegBCs) compared to published individual-CpG-based techniques when applied to external datasets. The RegBCs are also more robust when applied to low coverage data and detect a negative age acceleration in mice undergoing calorie restriction. Our RegBCs offer a proof of principle that age prediction of RRBS datasets can be improved by accounting for multiple CpGs over a region, which negates the lack of read depth currently hindering individual-CpG-based approaches.

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

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

PyMEGABASE: Predicting cell-type-specific structural annotations of chromosomes using the epigenome

The folding patterns of interphase genomes in higher eukaryotes, as obtained from DNA-proximity-ligation or Hi-C experiments, are used to classify loci into structural classes called compartments and subcompartments. These structurally annotated (sub)compartments are known to exhibit specific epigenomic characteristics and cell-type-specific variations. To explore the relationship between genome structure and the epigenome, we present PyMEGABASE (PYMB), a maximum-entropy-based neural network model that predicts (sub)compartment annotations of a locus based solely on the local epigenome, such as ChIP-Seq of histone post-translational modifications. PYMB builds upon our previous model while improving robustness, capability to handle diverse inputs and user-friendly implementation. We employed PYMB to predict subcompartments for over a hundred human cell types available in ENCODE, shedding light on the links between subcompartments, cell identity, and epigenomic signals. The fact that PYMB, trained on data for human cells, can accurately predict compartments in mice suggests that the model is learning underlying physicochemical principles transferable across cell types and species. Reliable at higher resolutions (up to 5 kbp), PYMB is used to investigate compartment-specific gene expression. Not only can PYMB generate (sub)compartment information without Hi-C experiments, but its predictions are also interpretable. Analyzing PYMB's trained parameters, we explore the importance of various epigenomic marks in each subcompartment prediction. Furthermore, the predictions of the model can be used as input for OpenMiChroM software, which has been calibrated to generate three-dimensional structures of the genome. Detailed documentation of PYMB is available at https://pymegabase.readthedocs.io, including an installation guide using pip or conda, and Jupyter/Colab notebook tutorials.

Genome-wide tracing to decipher nuclear organization

Nuclear organization impacts gene expression activity and cell phenotype. Our current understanding is mainly derived from ensemble-level sequencing studies that reflect the 3D genome structure of millions of cells. These approaches have provided invaluable details on the 3D organizations of the genome and their relation to other nuclear landmarks. However, they mostly lack the ability to provide multimodal information simultaneously at the single-cell level. In recent years, cutting-edge imaging technologies have risen to the challenge of simultaneously describing multiple components of the nuclear space at the single-cell level, paving the way for a deeper understanding of the genome structure-function relationship. This review will focus on the development and utilization of such technologies to gain a multi-component view of the nucleus at single-cell resolution, dissecting the complexity and heterogeneity of nuclear organization.

Photoselective sequencing: microscopically guided genomic measurements with subcellular resolution

In biological systems, spatial organization and function are interconnected. Here we present photoselective sequencing, a new method for genomic and epigenomic profiling within morphologically distinct regions. Starting with an intact biological specimen, photoselective sequencing uses targeted illumination to selectively unblock a photocaged fragment library, restricting the sequencing-based readout to microscopically identified spatial regions. We validate photoselective sequencing by measuring the chromatin accessibility profiles of fluorescently labeled cell types within the mouse brain and comparing with published data. Furthermore, by combining photoselective sequencing with a computational strategy for decomposing bulk accessibility profiles, we find that the oligodendrocyte-lineage-cell population is relatively enriched for oligodendrocyte-progenitor cells in the cortex versus the corpus callosum. Finally, we leverage photoselective sequencing at the subcellular scale to identify features of chromatin that are correlated with positioning at the nuclear periphery. These results collectively demonstrate that photoselective sequencing is a flexible and generalizable platform for exploring the interplay of spatial structures with genomic and epigenomic properties.

Women’s contribution in understanding how topoisomerases, supercoiling, and transcription control genome organization

One of the biggest paradoxes in biology is that human genome is roughly 2 m long, while the nucleus containing it is almost one million times smaller. To fit into the nucleus, DNA twists, bends and folds into several hierarchical levels of compaction. Still, DNA has to maintain a high degree of accessibility to be readily replicated and transcribed by proteins. How compaction and accessibility co-exist functionally in human cells is still a matter of debate. Here, we discuss how the torsional stress of the DNA helix acts as a buffer, regulating both chromatin compaction and accessibility. We will focus on chromatin supercoiling and on the emerging role of topoisomerases as pivotal regulators of genome organization. We will mainly highlight the major breakthrough studies led by women, with the intention of celebrating the work of this group that remains a minority within the scientific community.

Genome-wide identification of copy neutral loss of heterozygosity reveals its possible association with spatial positioning of chromosomes

Loss of heterozygosity (LOH) is a genetic alteration that results from the loss of one allele at a heterozygous locus. In particular, copy neutral LOH (CN-LOH) events are generated, for example, by mitotic homologous recombination after monoallelic defection or gene conversion, resulting in novel homozygous locus having two copies of the normal counterpart allele. This phenomenon can serve as a source of genome diversity and is associated with various diseases. To clarify the nature of the CN-LOH such as the frequency, genomic distribution and inheritance pattern, we made use of whole-genome sequencing data of the three-generation CEPH/Utah family cohort, with the pedigree consisting of grandparents, parents and offspring. We identified an average of 40.7 CN-LOH events per individual taking advantage of 285 healthy individuals from 33 families in the cohort. On average 65% of them were classified as gonosomal-mosaicism-associated CN-LOH, which exists in both germline and somatic cells. We also confirmed that the incidence of the CN-LOH has little to do with the parents' age and sex. Furthermore, through the analysis of the genomic region including the CN-LOH, we found that the chance of the occurrence of the CN-LOH tends to increase at the GC-rich locus and/or on the chromosome having a relatively close inter-homolog distance. We expect that these results provide significant insights into the association between genetic alteration and spatial position of chromosomes as well as the intrinsic genetic property of the CN-LOH.

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.