3D genome organization in health and disease: emerging opportunities in cancer translational medicine

scENCORE: leveraging single-cell epigenetic data to predict chromatin conformation using graph embedding

Dynamic compartmentalization of eukaryotic DNA into active and repressed states enables diverse transcriptional programs to arise from a single genetic blueprint, whereas its dysregulation can be strongly linked to a broad spectrum of diseases. While single-cell Hi-C experiments allow for chromosome conformation profiling across many cells, they are still expensive and not widely available for most labs. Here, we propose an alternate approach, scENCORE, to computationally reconstruct chromatin compartments from the more affordable and widely accessible single-cell epigenetic data. First, scENCORE constructs a long-range epigenetic correlation graph to mimic chromatin interaction frequencies, where nodes and edges represent genome bins and their correlations. Then, it learns the node embeddings to cluster genome regions into A/B compartments and aligns different graphs to quantify chromatin conformation changes across conditions. Benchmarking using cell-type-matched Hi-C experiments demonstrates that scENCORE can robustly reconstruct A/B compartments in a cell-type-specific manner. Furthermore, our chromatin confirmation switching studies highlight substantial compartment-switching events that may introduce substantial regulatory and transcriptional changes in psychiatric disease. In summary, scENCORE allows accurate and cost-effective A/B compartment reconstruction to delineate higher-order chromatin structure heterogeneity in complex tissues.

Three-dimensional nuclear architecture distinguishes thyroid cancer histotypes

Molecular markers can serve as diagnostic tools to support pathological analysis in thyroid neoplasms. However, because the same markers can be observed in some benign thyroid lesions, additional approaches are necessary to differentiate thyroid tumor subtypes, prevent overtreatment and tailor specific clinical management. This applies particularly to the recently described variant of thyroid cancer referred to as noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP). This variant has an estimated prevalence of 4.4% to 9.1% of all papillary thyroid carcinomas worldwide. We studied 60 thyroid lesions: 20 classical papillary thyroid carcinoma (CPTC), 20 follicular variant of PTC (FVPTC) and 20 NIFTP. We examined morphological and molecular features to identify parameters that can differentiate NIFTP from the other PTC subtypes. When blindly investigating the nuclear architecture of thyroid neoplasms, we observed that NIFTP has significantly longer telomeres than CPTC and FVPTC. Super-resolved 3D-structured illumination microscopy demonstrated that NIFTP is heterogeneous and that its nuclei contain more densely packed DNA and smaller interchromatin spaces than CPTC and FVPTC, a pattern that resembles normal thyroid tissue. These data are consistent with the observed indolent biological behavior and favorable prognosis associated with NIFTP, which lacks BRAFV600E mutations. Of note, next-generation thyroid oncopanel sequencing was unable to distinguish the thyroid cancer histotypes in our study cohort. In summary, our data suggest that 3D nuclear architecture can be a powerful analytical tool to diagnose and guide clinical management of NIFTP.

IChrom-Deep: An Attention-Based Deep Learning Model for Identifying Chromatin Interactions

Identification of chromatin interactions is crucial for advancing our knowledge of gene regulation. However, due to the limitations of high-throughput experimental techniques, there is an urgent need to develop computational methods for predicting chromatin interactions. In this study, we propose a novel attention-based deep learning model, termed IChrom-Deep, to identify chromatin interactions using sequence features and genomic features. The experimental results based on the datasets of three cell lines demonstrate that the IChrom-Deep achieves satisfactory performance and is superior to the previous methods. We also investigate the effect of DNA sequence and associated features and genomic features on chromatin interactions, and highlight the applicable scenarios of some features, such as sequence conservation and distance. Moreover, we identify a few genomic features that are extremely important across different cell lines, and IChrom-Deep achieves comparable performance with only these significant genomic features versus using all genomic features. It is believed that IChrom-Deep can serve as a useful tool for future studies that seek to identify chromatin interactions.

Haplotype-specific chromatin looping reveals genetic interactions of regulatory regions modulating gene expression in 8p23.1

A major goal of genetics research is to elucidate mechanisms explaining how genetic variation contributes to phenotypic variation. The genetic variants identified in genome-wide association studies (GWASs) generally explain only a small proportion of heritability of phenotypic traits, the so-called missing heritability problem. Recent evidence suggests that additional common variants beyond lead GWAS variants contribute to phenotypic variation; however, their mechanistic underpinnings generally remain unexplored. Herein, we undertake a study of haplotype-specific mechanisms of gene regulation at 8p23.1 in the human genome, a region associated with a number of complex diseases. The FAM167A-BLK locus in this region has been consistently found in the genome-wide association studies (GWASs) of systemic lupus erythematosus (SLE) in all major ancestries. Our haplotype-specific chromatin interaction (Hi-C) experiments, allele-specific enhancer activity measurements, genetic analyses, and epigenome editing experiments revealed that: 1) haplotype-specific long-range chromatin interactions are prevalent in 8p23.1; 2) BLK promoter and cis-regulatory elements cooperatively interact with haplotype-specificity; 3) genetic variants at distal regulatory elements are allele-specific modifiers of the promoter variants at FAM167A-BLK; 4) the BLK promoter interacts with and, as an enhancer-like promoter, regulates FAM167A expression and 5) local allele-specific enhancer activities are influenced by global haplotype structure due to chromatin looping. Although systemic lupus erythematosus causal variants at the FAM167A-BLK locus are thought to reside in the BLK promoter region, our results reveal that genetic variants at distal regulatory elements modulate promoter activity, changing BLK and FAM167A gene expression and disease risk. Our results suggest that global haplotype-specific 3-dimensional chromatin looping architecture has a strong influence on local allelic BLK and FAM167A gene expression, providing mechanistic details for how regional variants controlling the BLK promoter may influence disease risk.

Monocytes acquire prostate cancer specific chromatin conformations upon indirect co-culture with prostate cancer cells

Background

Three-dimensional chromosome loop conformations are powerful regulators of gene expression. These chromosome conformations can be detected both in tumour and in circulating cells and have significant disease biomarker potential. We have recently detected specific chromosome conformations in circulating cells of patients with prostate cancer (PCa) which were similar to ones found in their primary tumours, however, the possibility of horizontal transfer of chromosome conformations was not studied previously.

Methods

Human monocytes (U937) were co-cultured in Boyden chambers through 0.4 uM membrane with or without PC-3 human PCa cells or their conditioned media and a custom DNA microarray for 900,000 chromosomal loops covering all coding loci and non-coding RNA genes was performed on each part of the co-culture system.

Results

We have detected 684 PC-3 cell-specific chromosome conformations across the whole genome that were absent in naïve monocytes but appeared in monocytes co-cultured with PC-3 cells or with PC-3-conditioned media. Comparing PC3-specific conformations to the ones we have previously detected in systemic circulation of high-risk PCa patients revealed 9 positive loops present in both settings.

Conclusions

Our results demonstrate for the first time a proof of concept for horizontal transfer of chromosome conformations without direct cell-cell contact. This carries high clinical relevance as we have previously observed chromatin conformations in circulating cells of patients with melanoma and PCa similar to ones in their primary tumours. These changes can be used as highly specific biomarkers for diagnosis and prognosis. Further studies are required to elucidate the specific mechanism of chromosome conformations transfer and its clinical significance in particular diseases.

Identification of mechanism of cancer-cell-specific reactivation of hTERT offers therapeutic opportunities for blocking telomerase specifically in human colorectal cancer

Transcriptional reactivation of hTERT is the limiting step in tumorigenesis. While mutations in hTERT promoter present in 19% of cancers are recognized as key drivers of hTERT reactivation, mechanisms by which wildtype hTERT (WT-hTERT) promoter is reactivated, in majority of human cancers, remain unknown. Using primary colorectal cancers (CRC) we identified Tert INTeracting region 2 (T-INT2), the critical chromatin region essential for reactivating WT-hTERT promoter in CRCs. Elevated β-catenin and JunD level in CRC facilitates chromatin interaction between hTERT promoter and T-INT2 that is necessary to turn on hTERTexpression. Pharmacological screens uncovered salinomycin, which inhibits JunD mediated hTERT-T-INT2 interaction that is required for the formation of a stable transcription complex on the hTERT promoter. Our results showed for the first time how known CRC alterations, such as APC, lead to WT-hTERT promoter reactivation during stepwise-tumorigenesis and provide a new perspective for developing cancer-specific drugs.

MYC overexpression leads to increased chromatin interactions at superenhancers and MYC binding sites

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC. MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter–enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein–protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

Biophysics is reshaping our perception of the epigenome: from DNA-level to high-throughput studies

Epigenetic research holds great promise to advance our understanding of biomarkers and regulatory processes in health and disease. An increasing number of new approaches, ranging from molecular to biophysical analyses, enable identifying epigenetic changes on the level of a single gene or the whole epigenome. The aim of this review is to highlight how the field is shifting from completely molecular-biology-driven solutions to multidisciplinary strategies including more reliance on biophysical analysis tools. Biophysics not only offers technical advancements in imaging or structure analysis but also helps to explore regulatory interactions. New computational methods are also being developed to meet the demand of growing data volumes and their processing. Therefore, it is important to capture these new directions in epigenetics from a biophysical perspective and discuss current challenges as well as multiple applications of biophysical methods and tools. Specifically, we gradually introduce different biophysical research methods by first considering the DNA-level information and eventually higher-order chromatin structures. Moreover, we aim to highlight that the incorporation of bioinformatics, machine learning, and artificial intelligence into biophysical analysis allows gaining new insights into complex epigenetic processes. The gained understanding has already proven useful in translational and clinical research providing better patient stratification options or new therapeutic insights. Together, this offers a better readiness to transform bench-top experiments into industrial high-throughput applications with a possibility to employ developed methods in clinical practice and diagnostics.

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

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

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Application of Hi-C and other omics data analysis in human cancer and cell differentiation research

With the development of 3C (chromosome conformation capture) and its derivative technology Hi-C (High-throughput chromosome conformation capture) research, the study of the spatial structure of the genomic sequence in the nucleus helps researchers understand the functions of biological processes such as gene transcription, replication, repair, and regulation. In this paper, we first introduce the research background and purpose of Hi-C data visualization analysis. After that, we discuss the Hi-C data analysis methods from genome 3D structure, A/B compartment, TADs (topologically associated domain), and loop detection. We also discuss how to apply genome visualization technologies to the identification of chromosome feature structures. We continue with a review of correlation analysis differences among multi-omics data, and how to apply Hi-C and other omics data analysis into cancer and cell differentiation research. Finally, we summarize the various problems in joint analyses based on Hi-C and other multi-omics data. We believe this review can help researchers better understand the progress and applications of 3D genome technology.

H3K27me3-rich genomic regions can function as silencers to repress gene expression via chromatin interactions

The mechanisms underlying gene repression and silencers are poorly understood. Here we investigate the hypothesis that H3K27me3-rich regions of the genome, defined from clusters of H3K27me3 peaks, may be used to identify silencers that can regulate gene expression via proximity or looping. We find that H3K27me3-rich regions are associated with chromatin interactions and interact preferentially with each other. H3K27me3-rich regions component removal at interaction anchors by CRISPR leads to upregulation of interacting target genes, altered H3K27me3 and H3K27ac levels at interacting regions, and altered chromatin interactions. Chromatin interactions did not change at regions with high H3K27me3, but regions with low H3K27me3 and high H3K27ac levels showed changes in chromatin interactions. Cells with H3K27me3-rich regions knockout also show changes in phenotype associated with cell identity, and altered xenograft tumor growth. Finally, we observe that H3K27me3-rich regions-associated genes and long-range chromatin interactions are susceptible to H3K27me3 depletion. Our results characterize H3K27me3-rich regions and their mechanisms of functioning via looping.

Chromatin conformation changes in peripheral blood can detect prostate cancer and stratify disease risk groups

BACKGROUND

Current diagnostic blood tests for prostate cancer (PCa) are unreliable for the early stage disease, resulting in numerous unnecessary prostate biopsies in men with benign disease and false reassurance of negative biopsies in men with PCa. Predicting the risk of PCa is pivotal for making an informed decision on treatment options as the 5-year survival rate in the low-risk group is more than 95% and most men would benefit from surveillance rather than active treatment. Three-dimensional genome architecture and chromosome structures undergo early changes during tumourigenesis both in tumour and in circulating cells and can serve as a disease biomarker.

METHODS

In this prospective study we screened whole blood of newly diagnosed, treatment naïve PCa patients (n = 140) and cancer-free controls (n = 96) for the presence of 14,241 chromosomal loops in the loci of 425 genes.

RESULTS

We have detected specific chromosome conformation changes in the loci of ETS1, MAP3K14, SLC22A3 and CASP2 genes in peripheral blood from PCa patients yielding PCa detection with 80% sensitivity and 80% specificity. Further analysis between PCa risk groups yielded prognostic validation sets consisting of HSD3B2, VEGFC, APAF1, BMP6, ERG, MSR1, MUC1, ACAT1 and DAPK1 genes that achieved 80% sensitivity and 93% specificity stratifying high-risk category 3 vs low risk category 1 and 84% sensitivity and 89% specificity stratifying high risk category 3 vs intermediate risk category 2 disease.

CONCLUSIONS

Our results demonstrate specific chromosome conformations in the blood of PCa patients that allow PCa diagnosis and risk stratification with high sensitivity and specificity.

The Biological Significance of Multi-copy Regions and Their Impact on Variant Discovery

Identification of genetic variants via high-throughput sequencing (HTS) technologies has been essential for both fundamental and clinical studies. However, to what extent the genome sequence composition affects variant calling remains unclear. In this study, we identified 63,897 multi-copy sequences (MCSs) with a minimum length of 300 bp, each of which occurs at least twice in the human genome. The 151,749 genomic loci (multi-copy regions, or MCRs) harboring these MCSs account for 1.98% of the genome and are distributed unevenly across chromosomes. MCRs containing the same MCS tend to be located on the same chromosome. Gene Ontology (GO) analyses revealed that 3800 genes whose UTRs or exons overlap with MCRs are enriched for Golgi-related cellular component terms and various enzymatic activities in the GO biological function category. MCRs are also enriched for loci that are sensitive to neocarzinostatin-induced double-strand breaks. Moreover, genetic variants discovered by genome-wide association studies and recorded in dbSNP are significantly underrepresented in MCRs. Using simulated HTS datasets, we show that false variant discovery rates are significantly higher in MCRs than in other genomic regions. These results suggest that extra caution must be taken when identifying genetic variants in the MCRs via HTS technologies.

Model-based analysis of chromatin interactions from dCas9-Based CAPTURE-3C-seq

Deciphering long-range chromatin interactions is critical for understanding temporal and tissue-specific gene expression regulated by cis- and trans-acting factors. By combining the chromosome conformation capture (3C) and biotinylated dCas9 system, we previously established a method CAPTURE-3C-seq to unbiasedly identify high-resolution and locus-specific long-range DNA interactions. Here we present the statistical model and a flexible pipeline, C3S, for analysing CAPTURE-3C-seq or similar experimental data from raw sequencing reads to significantly interacting chromatin loci. C3S provides all steps for data processing, quality control and result illustration. It can automatically define the bin size based on the binding peak of the dCas9-targeted regions. Furthermore, it supports the analysis of intra- and inter-chromosomal interactions for different mammalian cell types. We successfully applied C3S across multiple datasets in human K562 cells and mouse embryonic stem cells (mESC) for detecting known and new chromatin interactions at multiple scales. Integrative and topological analysis of the interacted loci at the human β-globin gene cluster provides new insights into mechanisms in developmental gene regulation and network structure in local chromosomal architecture. Furthermore, computational results in mESCs reveal a role for chromatin interacting loops between enhancers and promoters in regulating alternative transcripts of the pluripotency gene OCT4.

Specific subfamilies of transposable elements contribute to different domains of T lymphocyte enhancers

Transposable elements (TEs) compose nearly half of mammalian genomes and provide building blocks for cis-regulatory elements. Using high-throughput sequencing, we show that 84 TE subfamilies are overrepresented, and distributed in a lineage-specific fashion in core and boundary domains of CD8⁺ T cell enhancers. Endogenous retroviruses are most significantly enriched in core domains with accessible chromatin, and bear recognition motifs for immune-related transcription factors. In contrast, short interspersed elements (SINEs) are preferentially overrepresented in nucleosome-containing boundaries. A substantial proportion of these SINEs harbor a high density of the enhancer-specific histone mark H3K4me1 and carry sequences that match enhancer boundary nucleotide composition. Motifs with regulatory features are better preserved within enhancer-enriched TE copies compared to their subfamily equivalents located in gene deserts. TE-rich and TE-poor enhancers associate with both shared and unique gene groups and are enriched in overlapping functions related to lymphocyte and leukocyte biology. The majority of T cell enhancers are shared with other immune lineages and are accessible in common hematopoietic progenitors. A higher proportion of immune tissue-specific enhancers are TE-rich compared to enhancers specific to other tissues, correlating with higher TE occurrence in immune gene-associated genomic regions. Our results suggest that during evolution, TEs abundant in these regions and carrying motifs potentially beneficial for enhancer architecture and immune functions were particularly frequently incorporated by evolving enhancers. Their putative selection and regulatory cooption may have accelerated the evolution of immune regulatory networks.

Regulation of epigenetic state by non-histone chromatin proteins and transcription factors: Implications in disease

Besides the fundamental components of the chromatin, DNA and octameric histone, the non-histone chromatin proteins and non-coding RNA play a critical role in the organization of functional chromatin domains. The non-histone chromatin proteins therefore regulate the transcriptional outcome in both physiological and pathophysiological state as well. They also help to maintain the epigenetic state of the genome indirectly. Several transcription factors and histone interacting factors also contribute in the maintenance of the epigenetic states, especially acetylation by the induction of autoacetylation ability of p300/CBP. Alterations of KAT activity have been found to be causally related to disease manifestation, and thus could be potential therapeutic target.

Alterations in Chromatin Folding Patterns in Cancer Variant-Enriched Loci

In this study, we focus on the following question: do genomic regions enriched in cancer variant mutations have significantly different chromatin folding patterns? We utilize publicly available Hi-C data to characterize chromatin folding patterns in healthy (GM12878) and cancer (K562) cells based on status of A/B compartmentalization and random vs non-random chromatin physical interactions. We then perform statistical testing to assess if chromatin folding patterns in cancer variant-enriched loci are significantly different from non-enriched loci. Our results indicate that loci with cancer variant status have significantly altered (FDR < 0.05) chromatin folding patterns.

Human Enhancers Harboring Specific Sequence Composition, Activity, and Genome Organization Are Linked to the Immune Response

The FANTOM5 consortium recently characterized 65,423 human enhancers from 1829 cell and tissue samples using the Cap Analysis of Gene Expression technology. We showed that the guanine and cytosine content at enhancer regions distinguishes two classes of enhancers harboring distinct DNA structural properties at flanking regions. A functional analysis of their predicted gene targets highlighted one class of enhancers as significantly enriched for associations with immune response genes. Moreover, these enhancers were specifically enriched for regulatory motifs recognized by transcription factors involved in immune response. We observed that enhancers enriched for links to immune response genes were more cell-type specific, preferentially activated upon bacterial infection, and with specific response activity. Looking at chromatin capture data, we found that the two classes of enhancers were lying in distinct topologically associating domains and chromatin loops. Our results suggest that specific nucleotide compositions encode for classes of enhancers that are functionally distinct and specifically organized in the human genome.

The role of 3D genome organization in disease: From compartments to single nucleotides

Since the advent of the chromosome conformation capture technology, our understanding of the human genome 3D organization has grown rapidly and we now know that human interphase chromosomes are folded into multiple layers of hierarchical structures and each layer can play a critical role in transcriptional regulation. Alterations in any one of these finely-tuned layers can lead to unwanted cascade of molecular events and ultimately drive the manifestation of diseases and phenotypes. Here we discuss, starting from chromosome level organization going down to single nucleotide changes, recent studies linking diseases or phenotypes to changes in the 3D genome architecture.

GrapHi-C: graph-based visualization of Hi-C datasets

Objectives

Hi-C is a proximity-based ligation reaction used to detect regions of the genome that are close in 3D space (or “interacting”). Typically, results from Hi-C experiments (contact maps) are visualized as heatmaps or Circos plots. While informative, these visualizations do not directly represent genomic structure and folding, making the interpretation of the underlying 3D genomic organization obscured. Our objective was to generate a graph-based contact map representation that leads to a more intuitive structural visualization.

Results

Normalized contact maps were converted into undirected graphs where each vertex represented a genomic region and each edge represented a detected (intra- and inter-chromosomal) or known (linear) interaction between two regions. Each edge was weighted by the inverse of the linear distance (Hi-C experimental resolution) or the interaction frequency from the contact map. Graphs were generated based on this representation scheme for contact maps from existing fission yeast datasets. Originally, these datasets were used to (1) identify specific principles influencing fission yeast genome organization and (2) uncover changes in fission yeast genome organization during the cell cycle. When compared to the equivalent heatmaps and/or Circos plots, the graph-based visualizations more intuitively depicted the changes in genome organization described in the original studies.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3507-2) contains supplementary material, which is available to authorized users.

Probing instructions for expression regulation in gene nucleotide compositions

Gene expression is orchestrated by distinct regulatory regions to ensure a wide variety of cell types and functions. A challenge is to identify which regulatory regions are active, what are their associated features and how they work together in each cell type. Several approaches have tackled this problem by modeling gene expression based on epigenetic marks, with the ultimate goal of identifying driving regions and associated genomic variations that are clinically relevant in particular in precision medicine. However, these models rely on experimental data, which are limited to specific samples (even often to cell lines) and cannot be generated for all regulators and all patients. In addition, we show here that, although these approaches are accurate in predicting gene expression, inference of TF combinations from this type of models is not straightforward. Furthermore these methods are not designed to capture regulation instructions present at the sequence level, before the binding of regulators or the opening of the chromatin. Here, we probe sequence-level instructions for gene expression and develop a method to explain mRNA levels based solely on nucleotide features. Our method positions nucleotide composition as a critical component of gene expression. Moreover, our approach, able to rank regulatory regions according to their contribution, unveils a strong influence of the gene body sequence, in particular introns. We further provide evidence that the contribution of nucleotide content can be linked to co-regulations associated with genome 3D architecture and to associations of genes within topologically associated domains.

[Computer methods of analysis of chromosome contacts in the cell nucleus based on sequencing technology data]

The study spatial chromosome structure and chromosome folding in the interphase cell nucleus is an important challenge of world science. Detection of eukaryotic genome regions that physically interact with each other could be done by modern sequencing technologies. A basic method of chromosome folding by total sequencing of contacting DNA fragments is HI-C. Long-range chromosomal interactions play an important role in gene transcription and regulation. The study of chromosome interactions, 3D (three-dimensional) genome structure and its effect on gene transcription allows revealing fundamental biological processes from a viewpoint of structural regulation and are important for cancer research. The technique of chromatin immunoprecipitation and subsequent sequencing (ChIP-seq) make possible to determine binding sites of transcription factors that regulate expression of eukaryotic genes; genome transcription factors binding maps have been. The ChIA-PET technology allows exploring not only target protein binding sites, but also pairs of such sites on proximally located and interacting with each other chromosomes co-located in three-dimensional space of the cell nucleus. Here we discuss the principles of the construction of genomic maps and matrices of chromosome contacts according to ChIA-PET and Hi-C data that capture the chromosome conformation and overview existing software for 3D genome analysis including in house programs of gene location analysis in topological domains.

Expanding the effects of ERG on chromatin landscapes and dysregulated transcription in prostate cancer

ERG overexpression in prostate cancers promotes the development of widespread changes in gene expression and chromatin landscapes, leading to redistribution of key transcription factors in prostate cancers positive for the TMPRSS2-ERG fusion gene. The overexpression of ERG is further assisted by the development of a super-enhancer in the ERG locus.

Super-Enhancers and Broad H3K4me3 Domains Form Complex Gene Regulatory Circuits Involving Chromatin Interactions

Stretched histone regions, such as super-enhancers and broad H3K4me3 domains, are associated with maintenance of cell identity and cancer. We connected super-enhancers and broad H3K4me3 domains in the K562 chronic myelogenous leukemia cell line as well as the MCF-7 breast cancer cell line with chromatin interactions. Super-enhancers and broad H3K4me3 domains showed higher association with chromatin interactions than their typical counterparts. Interestingly, we identified a subset of super-enhancers that overlap with broad H3K4me3 domains and show high association with cancer-associated genes including tumor suppressor genes. Besides cell lines, we could observe chromatin interactions by a Chromosome Conformation Capture (3C)-based method, in primary human samples. Several chromatin interactions involving super-enhancers and broad H3K4me3 domains are constitutive and can be found in both cancer and normal samples. Taken together, these results reveal a new layer of complexity in gene regulation by super-enhancers and broad H3K4me3 domains.

The Genome Conformation As an Integrator of Multi-Omic Data: The Example of Damage Spreading in Cancer

Publicly available multi-omic databases, in particular if associated with medical annotations, are rich resources with the potential to lead a rapid transition from high-throughput molecular biology experiments to better clinical outcomes for patients. In this work, we propose a model for multi-omic data integration (i.e., genetic variations, gene expression, genome conformation, and epigenetic patterns), which exploits a multi-layer network approach to analyse, visualize, and obtain insights from such biological information, in order to use achieved results at a macroscopic level. Using this representation, we can describe how driver and passenger mutations accumulate during the development of diseases providing, for example, a tool able to characterize the evolution of cancer. Indeed, our test case concerns the MCF-7 breast cancer cell line, before and after the stimulation with estrogen, since many datasets are available for this case study. In particular, the integration of data about cancer mutations, gene functional annotations, genome conformation, epigenetic patterns, gene expression, and metabolic pathways in our multi-layer representation will allow a better interpretation of the mechanisms behind a complex disease such as cancer. Thanks to this multi-layer approach, we focus on the interplay of chromatin conformation and cancer mutations in different pathways, such as metabolic processes, that are very important for tumor development. Working on this model, a variance analysis can be implemented to identify normal variations within each omics and to characterize, by contrast, variations that can be accounted to pathological samples compared to normal ones. This integrative model can be used to identify novel biomarkers and to provide innovative omic-based guidelines for treating many diseases, improving the efficacy of decision trees currently used in clinic.

4D nucleomes in single cells: what can computational modeling reveal about spatial chromatin conformation?

Genome-wide sequencing technologies enable investigations of the structural properties of the genome in various spatial dimensions. Here, we review computational techniques developed to model the three-dimensional genome in single cells versus ensembles of cells and assess their underlying assumptions. We further address approaches to study the spatio-temporal aspects of genome organization from single-cell data.