Chromatin conformation analysis of primary patient tissue using a low input Hi-C method

Obtention of viable cell suspensions from breast cancer tumor biopsies for 3D chromatin conformation and single-cell transcriptome analysis

Molecular and cellular characterization of tumors is essential due to the complex and heterogeneous nature of cancer. In recent decades, many bioinformatic tools and experimental techniques have been developed to achieve personalized characterization of tumors. However, sample handling continues to be a major challenge as limitations such as prior treatments before sample acquisition, the amount of tissue obtained, transportation, or the inability to process fresh samples pose a hurdle for experimental strategies that require viable cell suspensions. Here, we present an optimized protocol that allows the recovery of highly viable cell suspensions from breast cancer primary tumor biopsies. Using these cell suspensions we have successfully characterized genome architecture through Hi-C. Also, we have evaluated single-cell gene expression and the tumor cellular microenvironment through single-cell RNAseq. Both technologies are key in the detailed and personalized molecular characterization of tumor samples. The protocol described here is a cost-effective alternative to obtain viable cell suspensions from biopsies simply and efficiently.

Integrated analyses highlight interactions between the three-dimensional genome and DNA, RNA and epigenomic alterations in metastatic prostate cancer

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

Biochemical properties of chromatin domains define genome compartmentalization

Chromatin three-dimensional (3D) organization inside the cell nucleus determines the separation of euchromatin and heterochromatin domains. Their segregation results in the definition of active and inactive chromatin compartments, whereby the local concentration of associated proteins, RNA and DNA results in the formation of distinct subnuclear structures. Thus, chromatin domains spatially confined in a specific 3D nuclear compartment are expected to share similar epigenetic features and biochemical properties, in terms of accessibility and solubility. Based on this rationale, we developed the 4f-SAMMY-seq to map euchromatin and heterochromatin based on their accessibility and solubility, starting from as little as 10 000 cells. Adopting a tailored bioinformatic data analysis approach we reconstruct also their 3D segregation in active and inactive chromatin compartments and sub-compartments, thus recapitulating the characteristic properties of distinct chromatin states. A key novelty of the new method is the capability to map both the linear segmentation of open and closed chromatin domains, as well as their compartmentalization in one single experiment.

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

scGrapHiC: deep learning-based graph deconvolution for Hi-C using single cell gene expression

SUMMARY

Single-cell Hi-C (scHi-C) protocol helps identify cell-type-specific chromatin interactions and sheds light on cell differentiation and disease progression. Despite providing crucial insights, scHi-C data is often underutilized due to the high cost and the complexity of the experimental protocol. We present a deep learning framework, scGrapHiC, that predicts pseudo-bulk scHi-C contact maps using pseudo-bulk scRNA-seq data. Specifically, scGrapHiC performs graph deconvolution to extract genome-wide single-cell interactions from a bulk Hi-C contact map using scRNA-seq as a guiding signal. Our evaluations show that scGrapHiC, trained on seven cell-type co-assay datasets, outperforms typical sequence encoder approaches. For example, scGrapHiC achieves a substantial improvement of 23.2% in recovering cell-type-specific Topologically Associating Domains over the baselines. It also generalizes to unseen embryo and brain tissue samples. scGrapHiC is a novel method to generate cell-type-specific scHi-C contact maps using widely available genomic signals that enables the study of cell-type-specific chromatin interactions.

AVAILABILITY AND IMPLEMENTATION

The GitHub link: https://github.com/rsinghlab/scGrapHiC contains the source code of scGrapHiC and associated scripts to preprocess publicly available datasets to produce the results and visualizations we have discuss in this manuscript.

Hi-C, a chromatin 3D structure technique advancing the functional genomics of immune cells

The functional performance of immune cells relies on a complex transcriptional regulatory network. The three-dimensional structure of chromatin can affect chromatin status and gene expression patterns, and plays an important regulatory role in gene transcription. Currently available techniques for studying chromatin spatial structure include chromatin conformation capture techniques and their derivatives, chromatin accessibility sequencing techniques, and others. Additionally, the recently emerged deep learning technology can be utilized as a tool to enhance the analysis of data. In this review, we elucidate the definition and significance of the three-dimensional chromatin structure, summarize the technologies available for studying it, and describe the research progress on the chromatin spatial structure of dendritic cells, macrophages, T cells, B cells, and neutrophils.

'Enhancing' skeletal muscle and stem cells in three-dimensions: genome regulation of skeletal muscle in development and disease

The noncoding genome imparts important regulatory control over gene expression. In particular, gene enhancers represent a critical layer of control that integrates developmental and differentiation signals outside the cell into transcriptional outputs inside the cell. Recently, there has been an explosion in genomic techniques to probe enhancer control, function, and regulation. How enhancers are regulated and integrate signals in stem cell development and differentiation is largely an open question. In this review, we focus on the role gene enhancers play in muscle stem cell specification, differentiation, and progression. We pay specific attention toward the identification of muscle-specific enhancers, the binding of transcription factors to these enhancers, and how enhancers communicate to their target genes via three-dimensional looping.

Unzipped chromosome-level genomes reveal allopolyploid nematode origin pattern as unreduced gamete hybridization

The formation and consequences of polyploidization in animals with clonal reproduction remain largely unknown. Clade I root-knot nematodes (RKNs), characterized by parthenogenesis and allopolyploidy, show a widespread geographical distribution and extensive agricultural destruction. Here, we generated 4 unzipped polyploid RKN genomes and identified a putative novel alternative telomeric element. Then we reconstructed 4 chromosome-level assemblies and resolved their genome structures as AAB for triploid and AABB for tetraploid. The phylogeny of subgenomes revealed polyploid RKN origin patterns as hybridization between haploid and unreduced gametes. We also observed extensive chromosomal fusions and homologous gene expression decrease after polyploidization, which might offset the disadvantages of clonal reproduction and increase fitness in polyploid RKNs. Our results reveal a rare pathway of polyploidization in parthenogenic polyploid animals and provide a large number of high-precision genetic resources that could be used for RKN prevention and control.

Multidimensional scaling methods can reconstruct genomic DNA loops using Hi-C data properties

This paper proposes multidimensional scaling (MDS) applied to high-throughput chromosome conformation capture (Hi-C) data on genomic interactions to visualize DNA loops. Currently, the mechanisms underlying the regulation of gene expression are poorly understood, and where and when DNA loops are formed remains undetermined. Previous studies have focused on reproducing the entire three-dimensional structure of chromatin; however, identifying DNA loops using these data is time-consuming and difficult. MDS is an unsupervised method for reconstructing the original coordinates from a distance matrix. Here, MDS was applied to high-throughput chromosome conformation capture (Hi-C) data on genomic interactions to visualize DNA loops. Hi-C data were converted to distances by taking the inverse to reproduce loops via MDS, and the missing values were set to zero. Using the converted data, MDS was applied to the log-transformed genomic coordinate distances and this process successfully reproduced the DNA loops in the given structure. Consequently, the reconstructed DNA loops revealed significantly more DNA-transcription factor interactions involved in DNA loop formation than those obtained from previously applied methods. Furthermore, the reconstructed DNA loops were significantly consistent with chromatin immunoprecipitation followed by sequencing (ChIP-seq) peak positions. In conclusion, the proposed method is an improvement over previous methods for identifying DNA loops.

Improving comparative analyses of Hi-C data via contrastive self-supervised learning

Hi-C is a widely applied chromosome conformation capture (3C)-based technique, which has produced a large number of genomic contact maps with high sequencing depths for a wide range of cell types, enabling comprehensive analyses of the relationships between biological functionalities (e.g. gene regulation and expression) and the three-dimensional genome structure. Comparative analyses play significant roles in Hi-C data studies, which are designed to make comparisons between Hi-C contact maps, thus evaluating the consistency of replicate Hi-C experiments (i.e. reproducibility measurement) and detecting statistically differential interacting regions with biological significance (i.e. differential chromatin interaction detection). However, due to the complex and hierarchical nature of Hi-C contact maps, it remains challenging to conduct systematic and reliable comparative analyses of Hi-C data. Here, we proposed sslHiC, a contrastive self-supervised representation learning framework, for precisely modeling the multi-level features of chromosome conformation and automatically producing informative feature embeddings for genomic loci and their interactions to facilitate comparative analyses of Hi-C contact maps. Comprehensive computational experiments on both simulated and real datasets demonstrated that our method consistently outperformed the state-of-the-art baseline methods in providing reliable measurements of reproducibility and detecting differential interactions with biological meanings.

Chromosome conformation capture technologies as tools to detect structural variations and their repercussion in chromatin 3D configuration

3D genome organization regulates gene expression in different physiological and pathological contexts. Characterization of chromatin structure at different scales has provided information about how the genome organizes in the nuclear space, from chromosome territories, compartments of euchromatin and heterochromatin, topologically associated domains to punctual chromatin loops between genomic regulatory elements and gene promoters. In recent years, chromosome conformation capture technologies have also been used to characterize structural variations (SVs) de novo in pathological conditions. The study of SVs in cancer, has brought information about transcriptional misregulation that relates directly to the incidence and prognosis of the disease. For example, gene fusions have been discovered arising from chromosomal translocations that upregulate oncogenes expression, and other types of SVs have been described that alter large genomic regions encompassing many genes. However, studying SVs in 2D cannot capture all their regulatory implications in the genome. Recently, several bioinformatic tools have been developed to identify and classify SVs from chromosome conformation capture data and clarify how they impact chromatin structure in 3D, resulting in transcriptional misregulation. Here, we review recent literature concerning bioinformatic tools to characterize SVs from chromosome conformation capture technologies and exemplify their vast potential to rebuild the 3D landscape of genomes in cancer. The study of SVs from the 3D perspective can produce essential information about drivers, molecular targets, and disease evolution.

Characterization of Gene Regulatory Elements in Human Fetal Cortical Development: Enhancing Our Understanding of Neurodevelopmental Disorders and Evolution

The neocortex is the region that most distinguishes human brain from other mammals and primates [Annu Rev Genet. 2021 Nov;55(1):555-81]. Studying the development of human cortex is important in understanding the evolutionary changes occurring in humans relative to other primates, as well as in elucidating mechanisms underlying neurodevelopmental disorders. Cortical development is a highly regulated process, spatially and temporally coordinated by expression of essential transcriptional factors in response to signaling pathways [Neuron. 2019 Sep;103(6):980-1004]. Enhancers are the most well-understood cis-acting, non-protein-coding regulatory elements that regulate gene expression [Nat Rev Genet. 2014 Apr;15(4):272-86]. Importantly, given the conservation of both DNA sequence and molecular function of the majority of proteins across mammals [Genome Res. 2003 Dec;13(12):2507-18], enhancers [Science. 2015 Mar;347(6226):1155-9], which are far more divergent at the sequence level, likely account for the phenotypes that distinguish the human brain by changing the regulation of gene expression. In this review, we will revisit the conceptual framework of gene regulation during human brain development, as well as the evolution of technologies to study transcriptional regulation, with recent advances in genome biology that open a window allowing us to systematically characterize cis-regulatory elements in developing human brain [Hum Mol Genet. 2022 Oct;31(R1):R84-96]. We provide an update on work to characterize the suite of all enhancers in the developing human brain and the implications for understanding neuropsychiatric disorders. Finally, we discuss emerging therapeutic ideas that utilize our emerging knowledge of enhancer function.

Rewiring of the 3D genome during acquisition of carboplatin resistance in a triple-negative breast cancer patient-derived xenograft

Changes in the three-dimensional (3D) structure of the genome are an emerging hallmark of cancer. Cancer-associated copy number variants and single nucleotide polymorphisms promote rewiring of chromatin loops, disruption of topologically associating domains (TADs), active/inactive chromatin state switching, leading to oncogene expression and silencing of tumor suppressors. However, little is known about 3D changes during cancer progression to a chemotherapy-resistant state. We integrated chromatin conformation capture (Hi-C), RNA-seq, and whole-genome sequencing obtained from triple-negative breast cancer patient-derived xenograft primary tumors (UCD52) and carboplatin-resistant samples and found increased short-range (< 2 Mb) interactions, chromatin looping, formation of TAD, chromatin state switching into a more active state, and amplification of ATP-binding cassette transporters. Transcriptome changes suggested the role of long-noncoding RNAs in carboplatin resistance. Rewiring of the 3D genome was associated with TP53, TP63, BATF, FOS-JUN family of transcription factors and led to activation of aggressiveness-, metastasis- and other cancer-related pathways. Integrative analysis highlighted increased ribosome biogenesis and oxidative phosphorylation, suggesting the role of mitochondrial energy metabolism. Our results suggest that 3D genome remodeling may be a key mechanism underlying carboplatin resistance.

Dynamic antagonism between key repressive pathways maintains the placental epigenome

DNA and Histone 3 Lysine 27 methylation typically function as repressive modifications and operate within distinct genomic compartments. In mammals, the majority of the genome is kept in a DNA methylated state, whereas the Polycomb repressive complexes regulate the unmethylated CpG-rich promoters of developmental genes. In contrast to this general framework, the extra-embryonic lineages display non-canonical, globally intermediate DNA methylation levels, including disruption of local Polycomb domains. Here, to better understand this unusual landscape's molecular properties, we genetically and chemically perturbed major epigenetic pathways in mouse trophoblast stem cells. We find that the extra-embryonic epigenome reflects ongoing and dynamic de novo methyltransferase recruitment, which is continuously antagonized by Polycomb to maintain intermediate, locally disordered methylation. Despite its disorganized molecular appearance, our data point to a highly controlled equilibrium between counteracting repressors within extra-embryonic cells, one that can seemingly persist indefinitely without bistable features typically seen for embryonic forms of epigenetic regulation.

Influence of human peripheral blood samples preprocessing on the quality of Hi-C libraries

The genome-wide variant of the chromatin conformation capture technique (Hi-C) is a powerful tool for revealing patterns of genome spatial organization, as well as for understanding the effects of their disturbance on disease development. In addition, Hi-C can be used to detect chromosomal rearrangements, including balanced translocations and inversions. The use of the Hi-C method for the detection of chromosomal rearrangements is becoming more widespread. Modern high-throughput methods of genome analysis can effectively reveal point mutations and unbalanced chromosomal rearrangements. However, their sensitivity for determining translocations and inversions remains rather low. The storage of whole blood samples can affect the amount and integrity of genomic DNA, and it can distort the results of subsequent analyses if the storage was not under proper conditions. The Hi-C method is extremely demanding on the input material. The necessary condition for successfully applying Hi-C and obtaining high-quality data is the preservation of the spatial chromatin organization within the nucleus. The purpose of this study was to determine the optimal storage conditions of blood samples for subsequent Hi-C analysis. We selected 10 different conditions for blood storage and sample processing. For each condition, we prepared and sequenced Hi-C libraries. The quality of the obtained data was compared. As a result of the work, we formulated the requirements for the storage and processing of samples to obtain high-quality Hi-C data. We have established the minimum volume of blood sufficient for conducting Hi-C analysis. In addition, we have identified the most suitable methods for isolation of peripheral blood mononuclear cells and their long-term storage. The main requirement we have formulated is not to freeze whole blood.

Subtype-specific 3D genome alteration in acute myeloid leukaemia

Acute myeloid leukaemia (AML) represents a set of heterogeneous myeloid malignancies, and hallmarks include mutations in epigenetic modifiers, transcription factors and kinases1-5. The extent to which mutations in AML drive alterations in chromatin 3D structure and contribute to myeloid transformation is unclear. Here we use Hi-C and whole-genome sequencing to analyse 25 samples from patients with AML and 7 samples from healthy donors. Recurrent and subtype-specific alterations in A/B compartments, topologically associating domains and chromatin loops were identified. RNA sequencing, ATAC with sequencing and CUT&Tag for CTCF, H3K27ac and H3K27me3 in the same AML samples also revealed extensive and recurrent AML-specific promoter-enhancer and promoter-silencer loops. We validated the role of repressive loops on their target genes by CRISPR deletion and interference. Structural variation-induced enhancer-hijacking and silencer-hijacking events were further identified in AML samples. Hijacked enhancers play a part in AML cell growth, as demonstrated by CRISPR screening, whereas hijacked silencers have a downregulating role, as evidenced by CRISPR-interference-mediated de-repression. Finally, whole-genome bisulfite sequencing of 20 AML and normal samples revealed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Treatment of AML cells with a DNA hypomethylating agent and triple knockdown of DNMT1, DNMT3A and DNMT3B enabled the manipulation of DNA methylation to revert 3D genome organization and gene expression. Overall, this study provides a resource for leukaemia studies and highlights the role of repressive loops and hijacked cis elements in human diseases.

Comparison and critical assessment of single-cell Hi-C protocols

Advances in single-cell sequencing technologies make it possible to study the genome architecture in single cells. The rapid growth of the field has been fueled by the development of innovative single-cell Hi-C protocols. However, the protocols vary considerably in their efficiency, bias, scale and costs, and their relative advantages for different applications are unclear.

Here, we compare the two most commonly used single-cell Hi-C protocols. We use long-read sequencing to analyze molecular products of the Hi-C assay and show that whole-genome amplification step results in increased number of artifacts, larger coverage biases, and increased amount of noise compared to PCR-based amplification. Our comparison provides guidance for researchers studying chromatin architecture in individual cells.

High-resolution simulations of chromatin folding at genomic rearrangements in malignant B cells provide mechanistic insights into proto-oncogene deregulation

Genomic rearrangements are known to result in proto-oncogene deregulation in many cancers, but the link to 3D genome structure remains poorly understood. Here, we used the highly predictive heteromorphic polymer (HiP-HoP) model to predict chromatin conformations at the proto-oncogene CCND1 in healthy and malignant B cells. After confirming that the model gives good predictions of Hi-C data for the nonmalignant human B cell-derived cell line GM12878, we generated predictions for two cancer cell lines, U266 and Z-138. These possess genome rearrangements involving CCND1 and the immunoglobulin heavy locus (IGH), which we mapped using targeted genome sequencing. Our simulations showed that a rearrangement in U266 cells where a single IGH super-enhancer is inserted next to CCND1 leaves the local topologically associated domain (TAD) structure intact. We also observed extensive changes in enhancer-promoter interactions within the TAD, suggesting that it is the downstream chromatin remodeling which gives rise to the oncogene activation, rather than the presence of the inserted super-enhancer DNA sequence per se. Simulations of the IGH-CCND1 reciprocal translocation in Z-138 cells revealed that an oncogenic fusion TAD is created, encompassing CCND1 and the IGH super-enhancers. We predicted how the structure and expression of CCND1 changes in these different cell lines, validating this using qPCR and fluorescence in situ hybridization microscopy. Our work demonstrates the power of polymer simulations to predict differences in chromatin interactions and gene expression for different translocation breakpoints.

Integrating 3D genomic and epigenomic data to enhance target gene discovery and drug repurposing in transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) are popular approaches to test for association between imputed gene expression levels and traits of interest. Here, we propose an integrative method PUMICE (Prediction Using Models Informed by Chromatin conformations and Epigenomics) to integrate 3D genomic and epigenomic data with expression quantitative trait loci (eQTL) to more accurately predict gene expressions. PUMICE helps define and prioritize regions that harbor cis-regulatory variants, which outperforms competing methods. We further describe an extension to our method PUMICE +, which jointly combines TWAS results from single- and multi-tissue models. Across 79 traits, PUMICE + identifies 22% more independent novel genes and increases median chi-square statistics values at known loci by 35% compared to the second-best method, as well as achieves the narrowest credible interval size. Lastly, we perform computational drug repurposing and confirm that PUMICE + outperforms other TWAS methods.

Retinoic acid induced meiosis initiation in female germline stem cells by remodelling three-dimensional chromatin structure

Objectives

This study aimed to clarify the regulation and mechanism of meiotic initiation in FGSC development.

Materials and Methods

FGSCs were induced to differentiate into meiosis in differentiation medium. RNA sequencing was performed to analysis the difference of transcription level. High‐through chromosome conformation capture sequencing (Hi‐C) was performed to analysis changes of three‐dimensional chromatin structure. Chromosome conformation capture further confirmed a spatial chromatin loop. ChIP‐qPCR and dual luciferase reporter were used to test the interaction between Stimulated by retinoic acid gene 8 (STRA8) protein and Trip13 promoter.

Results

Compared with FGSCs, the average diameter of STRA8‐positive germ cells increased from 13 μm to 16.8 μm. Furthermore, there were 4788 differentially expressed genes between the two cell stages; Meiosis and chromatin structure‐associated terms were significantly enriched. Additionally, Hi‐C results showed that FGSCs underwent A/B compartment switching (switch rate was 29.81%), the number of topologically associating domains (TADs) increasing, the average size of TADs decreasing, and chromatin loop changes at genome region of Trip13 from undifferentiated stage to meiosis‐initiation stage. Furthermore, we validated that Trip13 promoter contacted distal enhancer to form spatial chromatin loop and STRA8 could bind Trip13 promoter to promote gene expression.

Conclusion

FGSCs underwent chromatin structure remodelling from undifferentiated stage to meiosis‐initiation stage, which facilitated STRA8 binding to Trip13 promoter and promoting its expression.

HiCuT: An efficient and low input method to identify protein-directed chromatin interactions

3D genome organization regulates gene expression, and disruption of these long-range (>20kB) DNA-protein interactions results in pathogenic phenotypes. Chromosome conformation methods in conjunction with chromatin immunoprecipitation were used to decipher protein-directed chromatin interactions. However, these methods required abundant starting material (>500,000 cells), sizable number of sequencing reads (>100 million reads), and elaborate data processing methods to reduce background noise, which limited their use in primary cells. Hi-C Coupled chromatin cleavage and Tagmentation (HiCuT) is a new transposase-assisted tagmentation method that generates high-resolution protein directed long-range chromatin interactions as efficiently as existing methods, HiChIP and ChIA-PET, despite using 100,000 cells (5-fold less) and 12 million sequencing reads (8-fold fewer). Moreover, HiCuT generates high resolution fragment libraries with low background signal that are easily interpreted with minimal computational processing. We used HiCuT in human primary skin cells to link previously identified single nucleotide polymorphisms (SNPs) in skin disease to candidate genes and to identify functionally relevant transcription factors in an unbiased manner. HiCuT broadens the capacity for genomic profiling in systems previously unmeasurable, including primary cells, human tissue samples, and rare cell populations, and may be a useful tool for all investigators studying human genetics and personalized epigenomics.

DNA is precisely packaged and organized within a nucleus to regulate gene expression. Altering this structure results in disease and developmental abnormalities. Current methods to probe 3D genome organization require a minimum of 500,000 cells, 100 million sequencing reads, and elaborate computational genomics skills. This limits general adoption and prevents use in small populations of cells, including primary tissues. We report a new method called Hi-C Coupled chromatin cleavage and Tagmentation (HiCuT). This method couples current 3D genome methods in conjunction with target specific tagmentation, an enzyme-assisted approach to cut and tag DNA. We benchmarked HiCuT against existing methods and found similar efficiency and specificity in assessing target specific 3D genome organization despite reducing assay requirements to 100,000 cells (5-fold less) and 12 million sequencing reads (8-fold fewer). Strikingly, HiCuT data exhibited reduced background noise and required minimal computational processing. Taken together, HiCuT broadens the capacity for 3D genome profiling in cell populations previously unmeasurable, including primary cells and human tissues, reduces sequencing costs, and lowers the need for computational expertise. HiCuT will benefit all investigators studying gene regulation and disease pathophysiology.

Parallel analysis of transcription, integration, and sequence of single HIV-1 proviruses

HIV-1-infected cells that persist despite antiretroviral therapy (ART) are frequently considered "transcriptionally silent," but active viral gene expression may occur in some cells, challenging the concept of viral latency. Applying an assay for profiling the transcriptional activity and the chromosomal locations of individual proviruses, we describe a global genomic and epigenetic map of transcriptionally active and silent proviral species and evaluate their longitudinal evolution in persons receiving suppressive ART. Using genome-wide epigenetic reference data, we show that proviral transcriptional activity is associated with activating epigenetic chromatin features in linear proximity of integration sites and in their inter- and intrachromosomal contact regions. Transcriptionally active proviruses were actively selected against during prolonged ART; however, this pattern was violated by large clones of virally infected cells that may outcompete negative selection forces through elevated intrinsic proliferative activity. Our results suggest that transcriptionally active proviruses are dynamically evolving under selection pressure by host factors.

Offspring production of haploid spermatid-like cells derived from mouse female germline stem cells with chromatin condensation

Background

During male meiosis, the Y chromosome can form perfect pairing with the X chromosome. However, it is unclear whether mammalian Female germline stem cells (FGSCs) without a Y chromosome can transdifferentiate into functional haploid spermatid-like cells (SLCs).

Results

We found that spermatogenesis was restarted by transplanting FGSCs into Kit^w/wv mutant testes. Complete meiosis and formation of SLCs was induced in vitro by testicular cells of Kit^w/wv mutant mice, cytokines and retinoic acid. Healthy offspring were produced by sperm and SLCs derived from the in vivo and in vitro transdifferentiation of FGSCs, respectively. Furthermore, high-throughput chromosome conformation capture sequencing(Hi-C-seq) and “bivalent” (H3K4me3-H3K27me3) micro chromatin immunoprecipitation sequencing (μChIP-seq) experiments showed that stimulated by retinoic acid gene 8 (STRA8)/protamine 1 (PRM1)-positive transdifferentiated germ cells (tGCs) and male germ cells (mGCs) display similar chromatin dynamics and chromatin condensation during in vitro spermatogenesis.

Conclusion

This study demonstrates that sperm can be produced from FGSCs without a Y chromosome. This suggests a strategy for dairy cattle breeding to produce only female offspring with a high-quality genetic background.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-021-00697-z.

CENP-A: A Histone H3 Variant with Key Roles in Centromere Architecture in Healthy and Diseased States

Centromeres are key architectural components of chromosomes. Here, we examine their construction, maintenance, and functionality. Focusing on the mammalian centromere- specific histone H3 variant, CENP-A, we highlight its coevolution with both centromeric DNA and its chaperone, HJURP. We then consider CENP-A de novo deposition and the importance of centromeric DNA recently uncovered with the added value from new ultra-long-read sequencing. We next review how to ensure the maintenance of CENP-A at the centromere throughout the cell cycle. Finally, we discuss the impact of disrupting CENP-A regulation on cancer and cell fate.

Isolation of Boechera stricta Developing Embryos for Hi-C

The possibility of analyzing chromatin topology in developing plant embryos is hampered by inaccessibility of the embryo sac, deeply embedded in the maternal seed tissue, following double fertilization. Here we describe a protocol to isolate, purify, and prepare developing Boechera stricta embryos for chromosome conformation capture-based methods as in situ Hi-C experiments. Early globular embryos can be isolated by air-pressure microaspiration, and subsequently washed to eliminate residual cells from the endosperm and maternal seed coat, allowing for pure sampling of selected stages of embryogenesis. This protocol allows for the possibility of comparing genome topology during plant embryonic differentiation since early until late embryo development stages.

Acute Protein Depletion Strategies to Functionally Dissect the 3D Genome

The organization of the genome inside the nucleus facilitates many nuclear processes. Because the nuclear genome is highly dynamic and often regulated by essential proteins, rapid depletion strategies are necessary to perform loss-of-function analyses. Fortunately, in recent years, various methods have been developed to manipulate the cellular levels of a protein directly and acutely. Here, we describe different methods that have been developed to rapidly deplete proteins from cells, with a focus on auxin inducible degron and dTAG methods, as these are most commonly used in 3D genome organization studies. We outline best practices for designing a knockin strategy, as well as generation and validation of knockin cell lines. Acute depletion strategies have been transformative for the study of the 3D genome and will be important tools for delineating the processes and factors that determine organization of the genome inside the nucleus.

In Situ HiC

In situ HiC uses the relative frequency of DNA-DNA ligation events to reconstruct the three-dimensional architecture of a genome. As such, restriction enzyme digested ends of genomic DNA within fixed nuclei are tagged with biotinylated dNTPs. DNA-DNA ligation events generated via proximity ligation are then captured, amplified and next generation sequenced to determine their linear genomic position, but also their three-dimensional relationship. Here, we describe these steps in detail.

Reorganization of the 3D genome pinpoints non-coding drivers of primary prostate tumors

Prostate cancer is a heterogeneous disease whose progression is linked to genome instability. However, the impact of this instability on the non-coding genome and its three-dimensional organization to aid progression is unclear. Using primary benign and tumor tissue, we find a high concordance in higher order three-dimensional genome organization. This concordance argues for constraints to the topology of prostate tumor genomes. Nonetheless, we identified changes in focal chromatin interactions, typical of loops bridging non-coding cis-regulatory elements, and showed how structural variants can induce these changes to guide cis-regulatory element hijacking. Such events resulted in opposing differential expression of genes found at antipodes of rearrangements. Collectively, these results argue that changes to focal chromatin interactions, as opposed to higher order genome organization, allow for aberrant gene regulation and are repeatedly mediated by structural variants in primary prostate cancer.

Offspring production of ovarian organoids derived from spermatogonial stem cells by defined factors with chromatin reorganization

Introduction

Fate determination of germline stem cells remains poorly understood at the chromatin structure level.

Objectives

Our research hopes to develop successful offspring production of ovarian organoids derived from spermatogonial stem cells (SSCs) by defined factors.

Methods

The offspring production from oocytes transdifferentiated from mouse SSCs with tracking of transplanted SSCs in vivo, single cell whole exome sequencing, and in 3D cell culture reconstitution of the process of oogenesis derived from SSCs. The defined factors were screened with ovarian organoids. We uncovered extensive chromatin reorganization during SSC conversion into induced germline stem cells (iGSCs) using high throughput chromosome conformation.

Results

We demonstrate successful production of offspring from oocytes transdifferentiated from mouse spermatogonial stem cells (SSCs). Furthermore, we demonstrate direct induction of germline stem cells (iGSCs) differentiated into functional oocytes by transduction of H19, Stella, and Zfp57 and inactivation of Plzf in SSCs after screening with ovarian organoids. We uncovered extensive chromatin reorganization during SSC conversion into iGSCs, which was highly similar to female germline stem cells. We observed that although topologically associating domains were stable during SSC conversion, chromatin interactions changed in a striking manner, altering 35% of inactive and active chromosomal compartments throughout the genome.

Conclusion

We demonstrate successful offspring production of ovarian organoids derived from SSCs by defined factors with chromatin reorganization. These findings have important implications in various areas including mammalian gametogenesis, genetic and epigenetic reprogramming, biotechnology, and medicine.

A look into retinal organoids: methods, analytical techniques, and applications

Inherited retinal diseases (IRDs) cause progressive loss of light-sensitive photoreceptors in the eye and can lead to blindness. Gene-based therapies for IRDs have shown remarkable progress in the past decade, but the vast majority of forms remain untreatable. In the era of personalised medicine, induced pluripotent stem cells (iPSCs) emerge as a valuable system for cell replacement and to model IRD because they retain the specific patient genome and can differentiate into any adult cell type. Three-dimensional (3D) iPSCs-derived retina-like tissue called retinal organoid contains all major retina-specific cell types: amacrine, bipolar, horizontal, retinal ganglion cells, Müller glia, as well as rod and cone photoreceptors. Here, we describe the main applications of retinal organoids and provide a comprehensive overview of the state-of-art analysis methods that apply to this model system. Finally, we will discuss the outlook for improvements that would bring the cellular model a step closer to become an established system in research and treatment development of IRDs.

Profiling of 3D Genome Organization in Nasopharyngeal Cancer Needle Biopsy Patient Samples by a Modified Hi-C Approach

Nasopharyngeal cancer (NPC), a cancer derived from epithelial cells in the nasopharynx, is a cancer common in China, Southeast Asia, and Africa. The three-dimensional (3D) genome organization of nasopharyngeal cancer is poorly understood. A major challenge in understanding the 3D genome organization of cancer samples is the lack of a method for the characterization of chromatin interactions in solid cancer needle biopsy samples. Here, we developed Biop-C, a modified in situ Hi-C method using solid cancer needle biopsy samples. We applied Biop-C to characterize three nasopharyngeal cancer solid cancer needle biopsy patient samples. We identified topologically associated domains (TADs), chromatin interaction loops, and frequently interacting regions (FIREs) at key oncogenes in nasopharyngeal cancer from the Biop-C heatmaps. We observed that the genomic features are shared at some important oncogenes, but the patients also display extensive heterogeneity at certain genomic loci. On analyzing the super enhancer landscape in nasopharyngeal cancer cell lines, we found that the super enhancers are associated with FIREs and can be linked to distal genes via chromatin loops in NPC. Taken together, our results demonstrate the utility of our Biop-C method in investigating 3D genome organization in solid cancers.

tagHi-C Reveals 3D Chromatin Architecture Dynamics during Mouse Hematopoiesis

Spatiotemporal chromatin reorganization during hematopoietic differentiation has not been comprehensively characterized, mainly because of the large numbers of starting cells required for current chromatin conformation capture approaches. Here, we introduce a low-input tagmentation-based Hi-C (tagHi-C) method to capture the chromatin structures of hundreds of cells. Using tagHi-C, we are able to map the spatiotemporal dynamics of chromatin structure in ten primary hematopoietic stem, progenitor, and differentiated cell populations from mouse bone marrow. Our results reveal that changes in compartment dynamics and the Rabl configuration occur during hematopoietic cell differentiation. We identify gene-body-associating domains (GADs) as general structures for highly expressed genes. Moreover, we extend the body of knowledge regarding genes influenced by genome-wide association study (GWAS) loci through spatial chromatin looping. Our study provides the tagHi-C method for studying the three-dimensional (3D) genome of a small number of cells and maps the comprehensive 3D chromatin landscape of bone marrow hematopoietic cells.

Four-Dimensional Chromosome Structure Prediction

Chromatin conformation plays an important role in a variety of genomic processes, including genome replication, gene expression, and gene methylation. Hi-C data is frequently used to analyze structural features of chromatin, such as AB compartments, topologically associated domains, and 3D structural models. Recently, the genomics community has displayed growing interest in chromatin dynamics. Here, we present 4DMax, a novel method, which uses time-series Hi-C data to predict dynamic chromosome conformation. Using both synthetic data and real time-series Hi-C data from processes, such as induced pluripotent stem cell reprogramming and cardiomyocyte differentiation, we construct smooth four-dimensional models of individual chromosomes. These predicted 4D models effectively interpolate chromatin position across time, permitting prediction of unknown Hi-C contact maps at intermittent time points. Furthermore, 4DMax correctly recovers higher order features of chromatin, such as AB compartments and topologically associated domains, even at time points where Hi-C data is not made available to the algorithm. Contact map predictions made using 4DMax outperform naïve numerical interpolation in 87.7% of predictions on the induced pluripotent stem cell dataset. A/B compartment profiles derived from 4DMax interpolation showed higher similarity to ground truth than at least one profile generated from a neighboring time point in 100% of induced pluripotent stem cell experiments. Use of 4DMax may alleviate the cost of expensive Hi-C experiments by interpolating intermediary time points while also providing valuable visualization of dynamic chromatin changes.

Spatial genome architecture and the emergence of malignancy

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

3D genome organization during lymphocyte development and activation

Chromosomes have a complex three-dimensional (3D) architecture comprising A/B compartments, topologically associating domains and promoter-enhancer interactions. At all these levels, the 3D genome has functional consequences for gene transcription and therefore for cellular identity. The development and activation of lymphocytes involves strict control of gene expression by transcription factors (TFs) operating in a three-dimensionally organized chromatin landscape. As lymphocytes are indispensable for tissue homeostasis and pathogen defense, and aberrant lymphocyte activity is involved in a wide range of human morbidities, acquiring an in-depth understanding of the molecular mechanisms that control lymphocyte identity is highly relevant. Here we review current knowledge of the interplay between 3D genome organization and transcriptional control during B and T lymphocyte development and antigen-dependent activation, placing special emphasis on the role of TFs.

Robust detection of translocations in lymphoma FFPE samples using targeted locus capture-based sequencing

In routine diagnostic pathology, cancer biopsies are preserved by formalin-fixed, paraffin-embedding (FFPE) procedures for examination of (intra-) cellular morphology. Such procedures inadvertently induce DNA fragmentation, which compromises sequencing-based analyses of chromosomal rearrangements. Yet, rearrangements drive many types of hematolymphoid malignancies and solid tumors, and their manifestation is instructive for diagnosis, prognosis, and treatment. Here, we present FFPE-targeted locus capture (FFPE-TLC) for targeted sequencing of proximity-ligation products formed in FFPE tissue blocks, and PLIER, a computational framework that allows automated identification and characterization of rearrangements involving selected, clinically relevant, loci. FFPE-TLC, blindly applied to 149 lymphoma and control FFPE samples, identifies the known and previously uncharacterized rearrangement partners. It outperforms fluorescence in situ hybridization (FISH) in sensitivity and specificity, and shows clear advantages over standard capture-NGS methods, finding rearrangements involving repetitive sequences which they typically miss. FFPE-TLC is therefore a powerful clinical diagnostics tool for accurate targeted rearrangement detection in FFPE specimens.

Chromatin architecture transitions from zebrafish sperm through early embryogenesis

Chromatin architecture mapping in 3D formats has increased our understanding of how regulatory sequences and gene expression are connected and regulated in a genome. The 3D chromatin genome shows extensive remodeling during embryonic development, and although the cleavage-stage embryos of most species lack structure before zygotic genome activation (pre-ZGA), zebrafish has been reported to have structure. Here, we aimed to determine the chromosomal architecture in paternal/sperm zebrafish gamete cells to discern whether it either resembles or informs early pre-ZGA zebrafish embryo chromatin architecture. First, we assessed the higher-order architecture through advanced low-cell in situ Hi-C. The structure of zebrafish sperm, packaged by histones, lacks topological associated domains and instead displays “hinge-like” domains of ∼150 kb that repeat every 1–2 Mbs, suggesting a condensed repeating structure resembling mitotic chromosomes. The pre-ZGA embryos lacked chromosomal structure, in contrast to prior work, and only developed structure post-ZGA. During post-ZGA, we find chromatin architecture beginning to form at small contact domains of a median length of ∼90 kb. These small contact domains are established at enhancers, including super-enhancers, and chemical inhibition of Ep300a (p300) and Crebbpa (CBP) activity, lowering histone H3K27ac, but not transcription inhibition, diminishes these contacts. Together, this study reveals hinge-like domains in histone-packaged zebrafish sperm chromatin and determines that the initial formation of high-order chromatin architecture in zebrafish embryos occurs after ZGA primarily at enhancers bearing high H3K27ac.

Deciphering the Complexity of 3D Chromatin Organization Driving Lymphopoiesis and Lymphoid Malignancies

Proper lymphopoiesis and immune responses depend on the spatiotemporal control of multiple processes, including gene expression, DNA recombination and cell fate decisions. High-order 3D chromatin organization is increasingly appreciated as an important regulator of these processes and dysregulation of genomic architecture has been linked to various immune disorders, including lymphoid malignancies. In this review, we present the general principles of the 3D chromatin topology and its dynamic reorganization during various steps of B and T lymphocyte development and activation. We also discuss functional interconnections between architectural, epigenetic and transcriptional changes and introduce major key players of genomic organization in B/T lymphocytes. Finally, we present how alterations in architectural factors and/or 3D genome organization are linked to dysregulation of the lymphopoietic transcriptional program and ultimately to hematological malignancies.

Comparison of the Hi-C, GAM and SPRITE methods using polymer models of chromatin

Hi-C, split-pool recognition of interactions by tag extension (SPRITE) and genome architecture mapping (GAM) are powerful technologies utilized to probe chromatin interactions genome wide, but how faithfully they capture three-dimensional (3D) contacts and how they perform relative to each other is unclear, as no benchmark exists. Here, we compare these methods in silico in a simplified, yet controlled, framework against known 3D structures of polymer models of murine and human loci, which can recapitulate Hi-C, GAM and SPRITE experiments and multiplexed fluorescence in situ hybridization (FISH) single-molecule conformations. We find that in silico Hi-C, GAM and SPRITE bulk data are faithful to the reference 3D structures whereas single-cell data reflect strong variability among single molecules. The minimal number of cells required in replicate experiments to return statistically similar contacts is different across the technologies, being lowest in SPRITE and highest in GAM under the same conditions. Noise-to-signal levels follow an inverse power law with detection efficiency and grow with genomic distance differently among the three methods, being lowest in GAM for genomic separations >1 Mb.

A cookbook for DNase Hi-C

Background

The Hi-C technique is widely employed to study the 3-dimensional chromatin architecture and to assemble genomes. The conventional in situ Hi-C protocol employs restriction enzymes to digest chromatin, which results in nonuniform genomic coverage. Using sequence-agnostic restriction enzymes, such as DNAse I, could help to overcome this limitation.

Results

In this study, we compare different DNAse Hi-C protocols and identify the critical steps that significantly affect the efficiency of the protocol. In particular, we show that the SDS quenching strategy strongly affects subsequent chromatin digestion. The presence of biotinylated oligonucleotide adapters may lead to ligase reaction by-products, which can be avoided by rational design of the adapter sequences. Moreover, the use of nucleotide-exchange enzymes for biotin fill-in enables simultaneous labelling and repair of DNA ends, similar to the conventional Hi-C protocol. These improvements simplify the protocol, making it less expensive and time-consuming.

Conclusions

We propose a new robust protocol for the preparation of DNAse Hi-C libraries from cultured human cells and blood samples supplemented with experimental controls and computational tools for the evaluation of library quality.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00389-5.

Freeze substitution Hi-C, a convenient and cost-effective method for capturing the natural 3D chromatin conformation from frozen samples

Chromatin interactions functionally affect genome architecture and gene regulation, but to date, only fresh samples must be used in High-through chromosome conformation capture (Hi-C) to keep natural chromatin conformation intact. This requirement has impeded the advancement of 3D genome research by limiting sample collection and storage options for researchers and severely limiting the number of samples that can be processed in a short time. Here, we develop a freeze substitution Hi-C (FS-Hi-C) technique that overcomes the need for fresh samples. FS-Hi-C can be used with samples stored in liquid nitrogen (LN2): the water in a vitreous form in the sample cells is replaced with ethanol via automated freeze substitution. After confirming that the FS step preserves the natural chromosome conformation during sample thawing, we tested the performance of FS-Hi-C with Drosophila melanogaster and Gossypium hirsutum. Beyond allowing the use of frozen samples and confirming that FS-Hi-C delivers robust data for generating contact heat maps and delineating A/B compartments and topologically associating domains, we found that FS-Hi-C outperforms the in situ Hi-C in terms of library quality, reproducibility, and valid interactions. Thus, FS-Hi-C will probably extend the application of 3D genome structure analysis to the vast number of experimental contexts in biological and medical research for which Hi-C methods have been unfeasible to date.

Structural variant detection in cancer genomes: computational challenges and perspectives for precision oncology

Cancer is generally characterized by acquired genomic aberrations in a broad spectrum of types and sizes, ranging from single nucleotide variants to structural variants (SVs). At least 30% of cancers have a known pathogenic SV used in diagnosis or treatment stratification. However, research into the role of SVs in cancer has been limited due to difficulties in detection. Biological and computational challenges confound SV detection in cancer samples, including intratumor heterogeneity, polyploidy, and distinguishing tumor-specific SVs from germline and somatic variants present in healthy cells. Classification of tumor-specific SVs is challenging due to inconsistencies in detected breakpoints, derived variant types and biological complexity of some rearrangements. Full-spectrum SV detection with high recall and precision requires integration of multiple algorithms and sequencing technologies to rescue variants that are difficult to resolve through individual methods. Here, we explore current strategies for integrating SV callsets and to enable the use of tumor-specific SVs in precision oncology.

Three-dimensional interactions between enhancers and promoters during intestinal differentiation depend upon HNF4

Cells in renewing tissues exhibit dramatic transcriptional changes as they differentiate. The contribution of chromatin looping to tissue renewal is incompletely understood. Enhancer-promoter interactions could be relatively stable as cells transition from progenitor to differentiated states; alternatively, chromatin looping could be as dynamic as the gene expression from their loci. The intestinal epithelium is the most rapidly renewing mammalian tissue. Proliferative cells in crypts of Lieberkühn sustain a stream of differentiated cells that are continually shed into the lumen. We apply chromosome conformation capture combined with chromatin immunoprecipitation (HiChIP) and sequencing to measure enhancer-promoter interactions in progenitor and differentiated cells of the intestinal epithelium. Despite dynamic gene regulation across the differentiation axis, we find that enhancer-promoter interactions are relatively stable. Functionally, we find HNF4 transcription factors are required for chromatin looping at target genes. Depletion of HNF4 disrupts local chromatin looping, histone modifications, and target gene expression. This study provides insights into transcriptional regulatory mechanisms governing homeostasis in renewing tissues.

Chen et al. provide a survey of enhancer-promoter 3D looping in the intestinal epithelium by HiChIP, in vivo. They find that enhancer-promoter interactions are highly dependent upon the key intestinal transcription factor HNF4. Their findings provide insights into transcriptional regulatory mechanisms governing homeostasis in renewing tissues.

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.

FAN-C: a feature-rich framework for the analysis and visualisation of chromosome conformation capture data

Chromosome conformation capture data, particularly from high-throughput approaches such as Hi-C, are typically very complex to analyse. Existing analysis tools are often single-purpose, or limited in compatibility to a small number of data formats, frequently making Hi-C analyses tedious and time-consuming. Here, we present FAN-C, an easy-to-use command-line tool and powerful Python API with a broad feature set covering matrix generation, analysis, and visualisation for C-like data ( https://github.com/vaquerizaslab/fanc ). Due to its compatibility with the most prevalent Hi-C storage formats, FAN-C can be used in combination with a large number of existing analysis tools, thus greatly simplifying Hi-C matrix analysis.

4D nucleome modeling

The intrinsic dynamic nature of chromosomes is emerging as a fundamental component in regulating DNA transcription, replication, and damage-repair among other nuclear functions. With this increased awareness, reinforced over the last ten years, many new experimental techniques, mainly based on microscopy and chromosome conformation capture, have been introduced to study the genome in space and time. Owing to the increasing complexity of these cutting-edge techniques, computational approaches have become of paramount importance to interpret, contextualize, and complement such experiments with new insights. Hence, it is becoming crucial for experimental biologists to have a clear understanding of the diverse theoretical modeling approaches available and the biological information each of them can provide.

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.

The Transition from Quiescent to Activated States in Human Hematopoietic Stem Cells Is Governed by Dynamic 3D Genome Reorganization

Lifelong blood production requires long-term hematopoietic stem cells (LT-HSCs), marked by stemness states involving quiescence and self-renewal, to transition into activated short-term HSCs (ST-HSCs) with reduced stemness. As few transcriptional changes underlie this transition, we used single-cell and bulk assay for transposase-accessible chromatin sequencing (ATAC-seq) on human HSCs and hematopoietic stem and progenitor cell (HSPC) subsets to uncover chromatin accessibility signatures, one including LT-HSCs (LT/HSPC signature) and another excluding LT-HSCs (activated HSPC [Act/HSPC] signature). These signatures inversely correlated during early hematopoietic commitment and differentiation. The Act/HSPC signature contains CCCTC-binding factor (CTCF) binding sites mediating 351 chromatin interactions engaged in ST-HSCs, but not LT-HSCs, enclosing multiple stemness pathway genes active in LT-HSCs and repressed in ST-HSCs. CTCF silencing derepressed stemness genes, restraining quiescent LT-HSCs from transitioning to activated ST-HSCs. Hence, 3D chromatin interactions centrally mediated by CTCF endow a gatekeeper function that governs the earliest fate transitions HSCs make by coordinating disparate stemness pathways linked to quiescence and self-renewal.