Inter-chromosomal contact networks provide insights into Mammalian chromatin organization

Systematic identification of interchromosomal interaction networks supports the existence of specialized RNA factories

Most studies of genome organization have focused on intrachromosomal (cis) contacts because they harbor key features such as DNA loops and topologically associating domains. Interchromosomal (trans) contacts have received much less attention, and tools for interrogating potential biologically relevant trans structures are lacking. Here, we develop a computational framework that uses Hi-C data to identify sets of loci that jointly interact in trans This method, trans-C, initiates probabilistic random walks with restarts from a set of seed loci to traverse an input Hi-C contact network, thereby identifying sets of trans-contacting loci. We validate trans-C in three increasingly complex models of established trans contacts: the Plasmodium falciparum var genes, the mouse olfactory receptor "Greek islands," and the human RBM20 cardiac splicing factory. We then apply trans-C to systematically test the hypothesis that genes coregulated by the same trans-acting element (i.e., a transcription or splicing factor) colocalize in three dimensions to form "RNA factories" that maximize the efficiency and accuracy of RNA biogenesis. We find that many loci with multiple binding sites of the same DNA-binding proteins interact with one another in trans, especially those bound by factors with intrinsically disordered domains. Similarly, clustered binding of a subset of RNA-binding proteins correlates with trans interaction of the encoding loci. We observe that these trans-interacting loci are close to nuclear speckles. These findings support the existence of trans- interacting chromatin domains (TIDs) driven by RNA biogenesis. Trans-C provides an efficient computational framework for studying these and other types of trans interactions, empowering studies of a poorly understood aspect of genome architecture.

Extrachromosomal DNA Associates with Nuclear Condensates and Reorganizes Chromatin Structures to Enhance Oncogenic Transcription

Extrachromosomal, circular DNA (ecDNA) is a prevalent oncogenic alteration in cancer genomes, often associated with aggressive tumor behavior and poor patient outcome. While previous studies proposed a chromatin-based mobile enhancer model for ecDNA-driven oncogenesis, its precise mechanism and impact remains unclear across diverse cancer types. Our study, utilizing advanced multi-omics profiling, epigenetic editing, and imaging approaches in three cancer models, reveals that ecDNA hubs are an integrated part of nuclear condensates and exhibit cancer-type specific chromatin connectivity. Epigenetic silencing of the ecDNA-specific regulatory modules or chemically disrupting liquid-liquid phase separation breaks down ecDNA hubs, displaces MED1 co-activator binding, inhibits oncogenic transcription, and promotes cell death. These findings substantiate the trans -activator function of ecDNA and underscore a structural mechanism driving oncogenesis. This refined understanding expands our views of oncogene regulation and opens potential avenues for novel therapeutic strategies in cancer treatment.

Potential roles of inter-chromosomal interactions in cell fate determination

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

Human Coronavirus Infection Reorganizes Spatial Genomic Architecture in Permissive Lung Cells

Chromatin conformation capture followed by next-generation sequencing in combination with large-scale polymer simulations (4DHiC) produces detailed information on genomic loci interactions, allowing for the interrogation of 3D spatial genomic structures. Here, Hi-C data was acquired from the infection of fetal lung fibroblast (MRC5) cells with α-coronavirus 229E (CoV229E). Experimental Hi-C contact maps were used to determine viral-induced changes in genomic architecture over a 48-hour time period following viral infection, revealing substantial alterations in contacts within chromosomes and in contacts between different chromosomes. To gain further structural insight and quantify the underlying changes, we applied the 4DHiC polymer simulation method to reconstruct the 3D genomic structures and dynamics corresponding to the Hi-C maps. The models successfully reproduced experimental Hi-C data, including the changes in contacts induced by viral infection. Our 3D spatial simulations uncovered widespread chromatin restructuring, including increased chromosome compactness and A-B compartment mixing arising from infection. Our model also suggests increased spatial accessibility to regions containing interferon-stimulated genes upon infection with CoV229E, followed by chromatin restructuring at later time points, potentially inducing the migration of chromatin into more compact regions. This is consistent with previously observed suppression of gene expression. Our spatial genomics study provides a mechanistic structural basis for changes in chromosome architecture induced by coronavirus infection in lung cells.

Systematic identification of inter-chromosomal interaction networks supports the existence of RNA factories

Most studies of genome organization have focused on intra-chromosomal (cis) contacts because they harbor key features such as DNA loops and topologically associating domains. Inter-chromosomal (trans) contacts have received much less attention, and tools for interrogating potential biologically relevant trans structures are lacking. Here, we develop a computational framework to identify sets of loci that jointly interact in trans from Hi-C data. This method, trans-C, initiates probabilistic random walks with restarts from a set of seed loci to traverse an input Hi-C contact network, thereby identifying sets of trans-contacting loci. We validate trans-C in three increasingly complex models of established trans contacts: the Plasmodium falciparum var genes, the mouse olfactory receptor "Greek islands", and the human RBM20 cardiac splicing factory. We then apply trans-C to systematically test the hypothesis that genes co-regulated by the same trans-acting element (i.e., a transcription or splicing factor) co-localize in three dimensions to form "RNA factories" that maximize the efficiency and accuracy of RNA biogenesis. We find that many loci with multiple binding sites of the same transcription factor interact with one another in trans, especially those bound by transcription factors with intrinsically disordered domains. Similarly, clustered binding of a subset of RNA binding proteins correlates with trans interaction of the encoding loci. These findings support the existence of trans interacting chromatin domains (TIDs) driven by RNA biogenesis. Trans-C provides an efficient computational framework for studying these and other types of trans interactions, empowering studies of a poorly understood aspect of genome architecture.

Mapping the semi-nested community structure of 3D chromosome contact networks

Mammalian DNA folds into 3D structures that facilitate and regulate genetic processes such as transcription, DNA repair, and epigenetics. Several insights derive from chromosome capture methods, such as Hi-C, which allow researchers to construct contact maps depicting 3D interactions among all DNA segment pairs. These maps show a complex cross-scale organization spanning megabase-pair compartments to short-ranged DNA loops. To better understand the organizing principles, several groups analyzed Hi-C data assuming a Russian-doll-like nested hierarchy where DNA regions of similar sizes merge into larger and larger structures. Apart from being a simple and appealing description, this model explains, e.g., the omnipresent chequerboard pattern seen in Hi-C maps, known as A/B compartments, and foreshadows the co-localization of some functionally similar DNA regions. However, while successful, this model is incompatible with the two competing mechanisms that seem to shape a significant part of the chromosomes' 3D organization: loop extrusion and phase separation. This paper aims to map out the chromosome's actual folding hierarchy from empirical data. To this end, we take advantage of Hi-C experiments and treat the measured DNA-DNA interactions as a weighted network. From such a network, we extract 3D communities using the generalized Louvain algorithm. This algorithm has a resolution parameter that allows us to scan seamlessly through the community size spectrum, from A/B compartments to topologically associated domains (TADs). By constructing a hierarchical tree connecting these communities, we find that chromosomes are more complex than a perfect hierarchy. Analyzing how communities nest relative to a simple folding model, we found that chromosomes exhibit a significant portion of nested and non-nested community pairs alongside considerable randomness. In addition, by examining nesting and chromatin types, we discovered that nested parts are often associated with active chromatin. These results highlight that cross-scale relationships will be essential components in models aiming to reach a deep understanding of the causal mechanisms of chromosome folding.

Satellite DNAs and human sex chromosome variation

Satellite DNAs are present on every chromosome in the cell and are typically enriched in repetitive, heterochromatic parts of the human genome. Sex chromosomes represent a unique genomic and epigenetic context. In this review, we first report what is known about satellite DNA biology on human X and Y chromosomes, including repeat content and organization, as well as satellite variation in typical euploid individuals. Then, we review sex chromosome aneuploidies that are among the most common types of aneuploidies in the general population, and are better tolerated than autosomal aneuploidies. This is demonstrated also by the fact that aging is associated with the loss of the X, and especially the Y chromosome. In addition, supernumerary sex chromosomes enable us to study general processes in a cell, such as analyzing heterochromatin dosage (i.e. additional Barr bodies and long heterochromatin arrays on Yq) and their downstream consequences. Finally, genomic and epigenetic organization and regulation of satellite DNA could influence chromosome stability and lead to aneuploidy. In this review, we argue that the complete annotation of satellite DNA on sex chromosomes in human, and especially in centromeric regions, will aid in explaining the prevalence and the consequences of sex chromosome aneuploidies.

Comparative 3D genome architecture in vertebrates

BACKGROUND

The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species.

RESULTS

This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture.

CONCLUSIONS

Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

Neuron-specific chromosomal megadomain organization is adaptive to recent retrotransposon expansions

Regulatory mechanisms associated with repeat-rich sequences and chromosomal conformations in mature neurons remain unexplored. Here, we map cell-type specific chromatin domain organization in adult mouse cerebral cortex and report strong enrichment of Endogenous Retrovirus 2 (ERV2) repeat sequences in the neuron-specific heterochromatic B₂^NeuN+ megabase-scaling subcompartment. Single molecule long-read sequencing and comparative Hi-C chromosomal contact mapping in wild-derived SPRET/EiJ (Mus spretus) and laboratory inbred C57BL/6J (Mus musculus) reveal neuronal reconfigurations tracking recent ERV2 expansions in the murine germline, with significantly higher B₂^NeuN+ contact frequencies at sites with ongoing insertions in Mus musculus. Neuronal ablation of the retrotransposon silencer Kmt1e/Setdb1 triggers B₂^NeuN+ disintegration and rewiring with open chromatin domains enriched for cellular stress response genes, along with severe neuroinflammation and proviral assembly with infiltration of dendrites . We conclude that neuronal megabase-scale chromosomal architectures include an evolutionarily adaptive heterochromatic organization which, upon perturbation, results in transcriptional dysregulation and unleashes ERV2 proviruses with strong neuronal tropism.

The role of sex chromosomes and sex hormones in vocal learning systems

Vocal learning is the ability to imitate and modify sounds through auditory experience, a rare trait found in only a few lineages of mammals and birds. It is a critical component of human spoken language, allowing us to verbally transmit speech repertoires and knowledge across generations. In many vocal learning species, the vocal learning trait is sexually dimorphic, where it is either limited to males or present in both sexes to different degrees. In humans, recent findings have revealed subtle sexual dimorphism in vocal learning/spoken language brain regions and some associated disorders. For songbirds, where the neural mechanisms of vocal learning have been well studied, vocal learning appears to have been present in both sexes at the origin of the lineage and was then independently lost in females of some subsequent lineages. This loss is associated with an interplay between sex chromosomes and sex steroid hormones. Even in species with little dimorphism, like humans, sex chromosomes and hormones still have some influence on learned vocalizations. Here we present a brief synthesis of these studies, in the context of sex determination broadly, and identify areas of needed investigation to further understand how sex chromosomes and sex steroid hormones help establish sexually dimorphic neural structures for vocal learning.

Oncogenic extrachromosomal DNA functions as mobile enhancers to globally amplify chromosomal transcription

Extrachromosomal, circular DNA (ecDNA) is emerging as a prevalent yet less characterized oncogenic alteration in cancer genomes. We leverage ChIA-PET and ChIA-Drop chromatin interaction assays to characterize genome-wide ecDNA-mediated chromatin contacts that impact transcriptional programs in cancers. ecDNAs in glioblastoma patient-derived neurosphere and prostate cancer cell cultures are marked by widespread intra-ecDNA and genome-wide chromosomal interactions. ecDNA-chromatin contact foci are characterized by broad and high-level H3K27ac signals converging predominantly on chromosomal genes of increased expression levels. Prostate cancer cells harboring synthetic ecDNA circles composed of characterized enhancers result in the genome-wide activation of chromosomal gene transcription. Deciphering the chromosomal targets of ecDNAs at single-molecule resolution reveals an association with actively expressed oncogenes spatially clustered within ecDNA-directed interaction networks. Our results suggest that ecDNA can function as mobile transcriptional enhancers to promote tumor progression and manifest a potential synthetic aneuploidy mechanism of transcription control in cancer.

RNA Biogenesis Instructs Functional Inter-Chromosomal Genome Architecture

Three-dimensional (3D) genome organization has emerged as an important layer of gene regulation in development and disease. The functional properties of chromatin folding within individual chromosomes (i.e., intra-chromosomal or in cis) have been studied extensively. On the other hand, interactions across different chromosomes (i.e., inter-chromosomal or in trans) have received less attention, being often regarded as background noise or technical artifacts. This viewpoint has been challenged by emerging evidence of functional relationships between specific trans chromatin interactions and epigenetic control, transcription, and splicing. Therefore, it is an intriguing possibility that the key processes involved in the biogenesis of RNAs may both shape and be in turn influenced by inter-chromosomal genome architecture. Here I present the rationale behind this hypothesis, and discuss a potential experimental framework aimed at its formal testing. I present a specific example in the cardiac myocyte, a well-studied post-mitotic cell whose development and response to stress are associated with marked rearrangements of chromatin topology both in cis and in trans. I argue that RNA polymerase II clusters (i.e., transcription factories) and foci of the cardiac-specific splicing regulator RBM20 (i.e., splicing factories) exemplify the existence of trans-interacting chromatin domains (TIDs) with important roles in cellular homeostasis. Overall, I propose that inter-molecular 3D proximity between co-regulated nucleic acids may be a pervasive functional mechanism in biology.

Nuclear gene proximity and protein interactions shape transcript covariations in mammalian single cells

Single-cell RNA sequencing studies on gene co-expression patterns could yield important regulatory and functional insights, but have so far been limited by the confounding effects of differentiation and cell cycle. We apply a tailored experimental design that eliminates these confounders, and report thousands of intrinsically covarying gene pairs in mouse embryonic stem cells. These covariations form a network with biological properties, outlining known and novel gene interactions. We provide the first evidence that miRNAs naturally induce transcriptome-wide covariations and compare the relative importance of nuclear organization, transcriptional and post-transcriptional regulation in defining covariations. We find that nuclear organization has the greatest impact, and that genes encoding for physically interacting proteins specifically tend to covary, suggesting importance for protein complex formation. Our results lend support to the concept of post-transcriptional RNA operons, but we further present evidence that nuclear proximity of genes may provide substantial functional regulation in mammalian single cells.

Gene expression covariation can be studied by single-cell RNA sequencing. Here the authors analyze intrinsically covarying gene pairs by eliminating the confounding effects in single-cell experiments and observe covariation of proximal genes and miRNA-induced covariation of target mRNAs.

Dynamic evolution of great ape Y chromosomes

The male-specific Y chromosome harbors genes important for sperm production. Because Y is repetitive, its DNA sequence was deciphered for only a few species, and its evolution remains elusive. Here we compared the Y chromosomes of great apes (human, chimpanzee, bonobo, gorilla, and orangutan) and found that many of their repetitive sequences and multicopy genes were likely already present in their common ancestor. Y repeats had increased intrachromosomal contacts, which might facilitate preservation of genes and gene regulatory elements. Chimpanzee and bonobo, experiencing high sperm competition, underwent many DNA changes and gene losses on the Y. Our research is significant for understanding the role of the Y chromosome in reproduction of nonhuman great apes, all of which are endangered.

The mammalian male-specific Y chromosome plays a critical role in sex determination and male fertility. However, because of its repetitive and haploid nature, it is frequently absent from genome assemblies and remains enigmatic. The Y chromosomes of great apes represent a particular puzzle: their gene content is more similar between human and gorilla than between human and chimpanzee, even though human and chimpanzee share a more recent common ancestor. To solve this puzzle, here we constructed a dataset including Ys from all extant great ape genera. We generated assemblies of bonobo and orangutan Ys from short and long sequencing reads and aligned them with the publicly available human, chimpanzee, and gorilla Y assemblies. Analyzing this dataset, we found that the genus Pan, which includes chimpanzee and bonobo, experienced accelerated substitution rates. Pan also exhibited elevated gene death rates. These observations are consistent with high levels of sperm competition in Pan. Furthermore, we inferred that the great ape common ancestor already possessed multicopy sequences homologous to most human and chimpanzee palindromes. Nonetheless, each species also acquired distinct ampliconic sequences. We also detected increased chromatin contacts between and within palindromes (from Hi-C data), likely facilitating gene conversion and structural rearrangements. Our results highlight the dynamic mode of Y chromosome evolution and open avenues for studies of male-specific dispersal in endangered great ape species.

Network-based method for regions with statistically frequent interchromosomal interactions at single-cell resolution

BACKGROUND

Chromosome conformation capture-based methods, especially Hi-C, enable scientists to detect genome-wide chromatin interactions and study the spatial organization of chromatin, which plays important roles in gene expression regulation, DNA replication and repair etc. Thus, developing computational methods to unravel patterns behind the data becomes critical. Existing computational methods focus on intrachromosomal interactions and ignore interchromosomal interactions partly because there is no prior knowledge for interchromosomal interactions and the frequency of interchromosomal interactions is much lower while the search space is much larger. With the development of single-cell technologies, the advent of single-cell Hi-C makes interrogating the spatial structure of chromatin at single-cell resolution possible. It also brings a new type of frequency information, the number of single cells with chromatin interactions between two disjoint chromosome regions.

RESULTS

Considering the lack of computational methods on interchromosomal interactions and the unsurprisingly frequent intrachromosomal interactions along the diagonal of a chromatin contact map, we propose a computational method dedicated to analyzing interchromosomal interactions of single-cell Hi-C with this new frequency information. To the best of our knowledge, our proposed tool is the first to identify regions with statistically frequent interchromosomal interactions at single-cell resolution. We demonstrate that the tool utilizing networks and binomial statistical tests can identify interesting structural regions through visualization, comparison and enrichment analysis and it also supports different configurations to provide users with flexibility.

CONCLUSIONS

It will be a useful tool for analyzing single-cell Hi-C interchromosomal interactions.

Evidence of constraint in the 3D genome for trans-splicing in human cells

Fusion transcripts are commonly found in eukaryotes, and many aberrant fusions are associated with severe diseases, including cancer. One class of fusion transcripts is generated by joining separate transcripts through trans-splicing. However, the mechanism of trans-splicing in mammals remains largely elusive. Here we showed evidence to support an intuitive hypothesis that attributes trans-sphcing to the spatial proximity between premature transcripts. A novel trans-splicing detection tool (TSD) was developed to reliably identify intra-chromosomal trans-splicing events (iTSEs) from RNA-seq data. TSD can maintain a remarkable balance between sensitivity and accuracy, thus distinguishing it from most state-of-the-art tools. The accuracy of TSD was experimentally demonstrated by excluding potential false discovery from mosaic genome or template switching during PCR. We showed that iTSEs identified by TSD were frequently found between genomic regulatory elements, which are known to be more prone to interact with each other. Moreover, iTSE sites may be more physically adjacent to each other than random control in the tested human lymphoblastoid cell line according to Hi-C data. Our results suggest that trans-splicing and 3D genome architecture may be coupled in mammals and that our pipeline, TSD, may facilitate investigations of trans-splicing on a systematic and accurate level previously thought impossible.

Three-dimensional Organization of Polytene Chromosomes in Somatic and Germline Tissues of Malaria Mosquitoes

Spatial organization of chromosome territories and interactions between interphase chromosomes themselves, as well as with the nuclear periphery, play important roles in epigenetic regulation of the genome function. However, the interplay between inter-chromosomal contacts and chromosome-nuclear envelope attachments in an organism’s development is not well-understood. To address this question, we conducted microscopic analyses of the three-dimensional chromosome organization in malaria mosquitoes. We employed multi-colored oligonucleotide painting probes, spaced 1 Mb apart along the euchromatin, to quantitatively study chromosome territories in larval salivary gland cells and adult ovarian nurse cells of Anopheles gambiae, An. coluzzii, and An. merus. We found that the X chromosome territory has a significantly smaller volume and is more compact than the autosomal arm territories. The number of inter-chromosomal, and the percentage of the chromosome–nuclear envelope, contacts were conserved among the species within the same cell type. However, the percentage of chromosome regions located at the nuclear periphery was typically higher, while the number of inter-chromosomal contacts was lower, in salivary gland cells than in ovarian nurse cells. The inverse correlation was considerably stronger for the autosomes. Consistent with previous theoretical arguments, our data indicate that, at the genome-wide level, there is an inverse relationship between chromosome-nuclear envelope attachments and chromosome–chromosome interactions, which is a key feature of the cell type-specific nuclear architecture.

Revealing Gene Function and Transcription Relationship by Reconstructing Gene-Level Chromatin Interaction

Chromatin is hierarchically organized in human interphase nuclei. Dynamic chromatin interactions are thought to influence gene transcription and cell fate determination. A consensus concept is that genes may form transcription factories within nucleus by spatially interaction. However, it is still not well known whether the function-related genes co-locate in three-dimensional (3D) space for co-transcription. Especially, there is a lack of visualization method that directly reflect the relationship between gene spatial interaction, gene function and co-transcription. In this study, we constructed three kinds of matrices based on gene ontology annotations, high-through chromosome conformation capture (Hi-C) data and RNA-seq data from twenty human tissues and cell lines. The comparative analysis for gene pairs revealed that 3D genome organization influences gene transcription predominantly at local scale. We found that the local genes within family clusters have similar transcription patterns. We also found that spatial reorganization of a histone gene cluster could control gene transcription. These observations suggest that function-related genes are close in space and activated or repressed together. Our work provided a framework for genome-wide studying the relationship between gene function, co-transcription and spatial interaction.

Identification of Genomic Alterations Through Multilevel DNA Structural Analysis

Current methods to identify genomic alterations using whole-genome sequencing (WGS) data are often limited to single nucleotide polymorphisms and insertions and deletions that are less than 10 bp in length. These limitations are largely due to challenges in accurately mapping short sequencing reads that significantly diverge from the reference genome. Newer sequencing-based methods have been developed to define and characterize larger DNA structural elements. This is achieved by enriching for and sequencing regions of the genome that contain a specific element, followed by identifying genomic regions with high densities of mapped short reads that designate the location of these elements. This process essentially aggregates short read data into larger structural units for further characterization. Here, we describe protocols for identifying various types of genomic alterations using differential analysis of these structural units. We focus on changes in DNA accessibility, protein-DNA interactions, and chromosomal contacts as measured by ATAC-Seq, ChIP-Seq, and Hi-C respectively. As many protocols have been published describing the generation and processing of these data, we focus on simple methods that can be used to identify mutations in these data, and can be executed by someone with limited computational expertise.

Nuclear organization mediates cancer-compromised genetic and epigenetic control

Nuclear organization is functionally linked to genetic and epigenetic regulation of gene expression for biological control and is modified in cancer. Nuclear organization supports cell growth and phenotypic properties of normal and cancer cells by facilitating physiologically responsive interactions of chromosomes, genes and regulatory complexes at dynamic three-dimensional microenvironments. We will review nuclear structure/function relationships that include: 1. Epigenetic bookmarking of genes by phenotypic transcription factors to control fidelity and plasticity of gene expression as cells enter and exit mitosis; 2. Contributions of chromatin remodeling to breast cancer nuclear morphology, metabolism and effectiveness of chemotherapy; 3. Relationships between fidelity of nuclear organization and metastasis of breast cancer to bone; 4. Dynamic modifications of higher-order inter- and intra-chromosomal interactions in breast cancer cells; 5. Coordinate control of cell growth and phenotype by tissue-specific transcription factors; 6. Oncofetal epigenetic control by bivalent histone modifications that are functionally related to sustaining the stem cell phenotype; and 7. Noncoding RNA-mediated regulation in the onset and progression of breast cancer. The discovery of components to nuclear organization that are functionally related to cancer and compromise gene expression have the potential for translation to innovative cancer diagnosis and targeted therapy.

Computational Methods for Assessing Chromatin Hierarchy

The hierarchical organization of chromatin is known to associate with diverse cellular functions; however, the precise mechanisms and the 3D structure remain to be determined. With recent advances in high-throughput next generation sequencing (NGS) techniques, genome-wide profiling of chromatin structures is made possible. Here, we provide a comprehensive overview of NGS-based methods for profiling "higher-order" and "primary-order" chromatin structures from both experimental and computational aspects. Experimental requirements and considerations specific for each method were highlighted. For computational analysis, we summarized a common analysis strategy for both levels of chromatin assessment, focusing on the characteristic computing steps and the tools. The recently developed single-cell level techniques based on Hi-C and ATAC-seq present great potential to reveal cell-to-cell variability in chromosome architecture. A brief discussion on these methods in terms of experimental and data analysis features is included. We also touch upon the biological relevance of chromatin organization and how the combination with other techniques uncovers the underlying mechanisms. We conclude with a summary and our prospects on necessary improvements of currently available methods in order to advance understanding of chromatin hierarchy. Our review brings together the analyses of both higher- and primary-order chromatin structures, and serves as a roadmap when choosing appropriate experimental and computational methods for assessing chromatin hierarchy.

Structure of the human chromosome interaction network

New Hi-C technologies have revealed that chromosomes have a complex network of spatial contacts in the cell nucleus of higher organisms, whose organisation is only partially understood. Here, we investigate the structure of such a network in human GM12878 cells, to derive a large scale picture of nuclear architecture. We find that the intensity of intra-chromosomal interactions is power-law distributed. Inter-chromosomal interactions are two orders of magnitude weaker and exponentially distributed, yet they are not randomly arranged along the genomic sequence. Intra-chromosomal contacts broadly occur between epigenomically homologous regions, whereas inter-chromosomal contacts are especially associated with regions rich in highly expressed genes. Overall, genomic contacts in the nucleus appear to be structured as a network of networks where a set of strongly individual chromosomal units, as envisaged in the 'chromosomal territory' scenario derived from microscopy, interact with each other via on average weaker, yet far from random and functionally important interactions.

FOXO3 longevity interactome on chromosome 6

FOXO3 has been implicated in longevity in multiple populations. By DNA sequencing in long‐lived individuals, we identified all single nucleotide polymorphisms (SNPs) in FOXO3 and showed 41 were associated with longevity. Thirteen of these had predicted alterations in transcription factor binding sites. Those SNPs appeared to be in physical contact, via RNA polymerase II binding chromatin looping, with sites in the FOXO3 promoter, and likely function together as a cis‐regulatory unit. The SNPs exhibited a high degree of LD in the Asian population, in which they define a specific longevity haplotype that is relatively common. The haplotype was less frequent in whites and virtually nonexistent in Africans. We identified distant contact points between FOXO3 and 46 neighboring genes, through long‐range physical contacts via CCCTC‐binding factor zinc finger protein (CTCF) binding sites, over a 7.3 Mb distance on chromosome 6q21. When activated by cellular stress, we visualized movement of FOXO3 toward neighboring genes. FOXO3 resides at the center of this early‐replicating and highly conserved syntenic region of chromosome 6. Thus, in addition to its role as a transcription factor regulating gene expression genomewide, FOXO3 may function at the genomic level to help regulate neighboring genes by virtue of its central location in chromatin conformation via topologically associated domains. We believe that the FOXO3 ‘interactome’ on chromosome 6 is a chromatin domain that defines an aging hub. A more thorough understanding of the functions of these neighboring genes may help elucidate the mechanisms through which FOXO3 variants promote longevity and healthy aging.

Superresolution imaging of nanoscale chromosome contacts

Co-expression of a specific group of genes requires physical associations among these genes, which form functional chromosomal contacts. While DNA fluorescence in situ hybridization (FISH) pinpoints the localization of genes within the 3D nuclear architecture, direct evidence of physical chromosomal contacts is still lacking. Here, we report a method for the direct visualization of transcription-dependent chromosomal contacts formed in two distinct mechanical states of cells. We prepared open chromatin spreads from isolated nuclei, ensuring 2D rendering of chromosome organization. Superresolution imaging of these chromatin spreads resolved the nanoscale organization of genome contacts. We optimized our imaging method using chromatin spreads from serum+/− cells. We then showed direct visualization of functional gene clusters targeted by YAP (Yes-associated protein) and SRF (Serum response factor) transcription factors. In addition, we showed the association of NF-κB bound gene clusters induced by TNF-α addition. Furthermore, EpiTect ChIP qPCR results showed that these nanoscale clusters were enriched with corresponding transcription factors. Taken together, our method provides a robust platform to directly visualize and study specific genome-wide chromosomal contacts.

Long noncoding RNA TUG1 is downregulated in non-small cell lung cancer and can regulate CELF1 on binding to PRC2

Background

Long noncoding RNAs (lncRNAs) play crucial roles in tumorigenesis, and lncRNA taurine-upregulated gene 1 (TUG1) has been proven to be associated with several human cancers. However, the mechanisms of TUG1-involved regulation remain largely unknown.

Methods

We examined the expressions of TUG1 in a cohort of 89 patients with non-small cell lung cancer (NSCLC) to determine the association between TUG1 expression and clinical parameters. We used circular chromosome conformation capture (4C) coupled with next-generation sequencing to explore the genome regions that interact with TUG1 and the TUG1-mediated regulation.

Results

TUG1 was significantly downregulated, and the TUG1 downregulation correlated with sex (p = 0.006), smoking status (p = 0.016), and tumor differentiation grade (p = 0.001). Knockdown of TUG1 significantly promoted the proliferation of NSCLC cells. According to the bioinformatic analysis result of TUG1 4C sequencing data, 83 candidate genes and their interaction regions were identified. Among these candidate genes, CUGBP and Elav-like family member 1 (CELF1) are potential targets of TUG1 in-trans regulation. To confirm the interaction between TUG1 and CELF1, relative expressions of CELF1 were examined in TUG1 knockdown H520 cells; results showed that CELF1 was significantly upregulated in TUG1 knockdown H520 cells. RNA immunoprecipitation was then performed to examine whether TUG1 RNA was bound to PRC2, a TUG1-involved regulation mechanism reported in previous studies. The results demonstrated that TUG1 RNA was bound to enhancer of zeste protein 2/embryonic ectoderm development (EZH2/EED), which is essential for PRC2. Finally, our designed ChIP assay revealed that the EZH2/EED was bound to the promotor region of CELF1 within 992 bp upstream of the transcript start site.

Conclusion

TUG1 is downregulated in NSCLC. Using TUG1 4C sequencing and bioinformatic analysis, we found CELF1 to be a potential target of TUG1 RNA in in-trans regulation. Moreover, subsequent experiments showed that TUG1 RNA could bind to PRC2 in the promotor region of CELF1 and negatively regulate CELF1 expressions in H520 cells. Our results may facilitate developing new treatment modalities targeting TUG1/PRC2/CELF1 interactions in patients with NSCLC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2569-6) contains supplementary material, which is available to authorized users.

Exploring the mechanisms of genome-wide long-range interactions: interpreting chromosome organization

Developments in chromosome conformation capture (3C) technologies have revealed that the three-dimensional organization of a genome leads widely separated functional elements to reside in close proximity. However, the mechanisms responsible for mediating long-range interactions are still not completely known. In this review, we firstly evaluate and compare the current seven 3C-based methods, summarize their advantages and discuss their limitations to our current understanding of genome structure. Then, software packages available to perform the analysis of 3C-based data are described. Moreover, we review the insights into the two main mechanisms of long-range interactions, which regulate gene expression by bringing together promoters and distal regulatory elements and by creating structural domains that contain functionally related genes with similar expression landscape. At last, we summarize what is known about the mediating factors involved in stimulation/repression of long-range interactions, such as transcription factors and noncoding RNAs.