Network analysis identifies chromosome intermingling regions as regulatory hotspots for transcription

The effect of trisomic chromosomes on spatial genome organization and global transcription in embryonic stem cells

Aneuploidy frequently occurs in cancer and developmental diseases such as Down syndrome, with its functional consequences implicated in dosage effects on gene expression and global perturbation of stress response and cell proliferation pathways. However, how aneuploidy affects spatial genome organization remains less understood. In this study, we addressed this question by utilizing the previously established isogenic wild-type (WT) and trisomic mouse embryonic stem cells (mESCs). We employed a combination of Hi-C, RNA-seq, chromosome painting and nascent RNA imaging technologies to compare the spatial genome structures and gene transcription among these cells. We found that trisomy has little effect on spatial genome organization at the level of A/B compartment or topologically associating domain (TAD). Inter-chromosomal interactions are associated with chromosome regions with high gene density, active histone modifications and high transcription levels, which are confirmed by imaging. Imaging also revealed contracted chromosome volume and weakened transcriptional activity for trisomic chromosomes, suggesting potential implications for the transcriptional output of these chromosomes. Our data resources and findings may contribute to a better understanding of the consequences of aneuploidy from the angle of spatial genome organization.

Transcriptional changes are tightly coupled to chromatin reorganization during cellular aging

Human life expectancy is constantly increasing and aging has become a major risk factor for many diseases, although the underlying gene regulatory mechanisms are still unclear. Using transcriptomic and chromosomal conformation capture (Hi-C) data from human skin fibroblasts from individuals across different age groups, we identified a tight coupling between the changes in co-regulation and co-localization of genes. We obtained transcription factors, cofactors, and chromatin regulators that could drive the cellular aging process by developing a time-course prize-collecting Steiner tree algorithm. In particular, by combining RNA-Seq data from different age groups and protein-protein interaction data we determined the key transcription regulators and gene regulatory changes at different life stage transitions. We then mapped these transcription regulators to the 3D reorganization of chromatin in young and old skin fibroblasts. Collectively, we identified key transcription regulators whose target genes are spatially rearranged and correlate with changes in their expression, thereby providing potential targets for reverting cellular aging.

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

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

Mechanical forces and the 3D genome

Traditionally, the field of genomics has been studied from a biochemical perspective. Besides chemical influences, cells are subject to a variety of mechanical signals from their surrounding tissue microenvironment. These mechanical signals can not only cause changes to a cell's physical structure but can also lead to alterations in their genomes and gene expression programs. Understanding the mechanical control of genome organization and expression may provide a new perspective on gene regulation.

Mesoscale, long-time mixing of chromosomes and its connection to polymer dynamics

Chromosomes are arranged in distinct territories within the nucleus of animal cells. Recent experiments have shown that these territories overlap at their edges, suggesting partial mixing during interphase. Experiments that knock-down of condensin II proteins during interphase indicate increased chromosome mixing, which demonstrates control of the mixing. In this study, we use a generic polymer simulation to quantify the dynamics of chromosome mixing over time. We introduce the chromosome mixing index, which quantifies the mixing of distinct chromosomes in the nucleus. We find that the chromosome mixing index in a small confinement volume (as a model of the nucleus), increases as a power-law of the time, with the scaling exponent varying non-monotonically with self-interaction and volume fraction. By comparing the chromosome mixing index with both monomer subdiffusion due to (non-topological) intermingling of chromosomes as well as even slower reptation, we show that for relatively large volume fractions, the scaling exponent of the chromosome mixing index is related to Rouse dynamics for relatively weak chromosome attractions and to reptation for strong attractions. In addition, we extend our model to more realistically account for the situation of the Drosophila chromosome by including the heterogeneity of the polymers and their lengths to account for microphase separation of euchromatin and heterochromatin and their interactions with the nuclear lamina. We find that the interaction with the lamina further impedes chromosome mixing.

Integrative chromatin domain annotation through graph embedding of Hi-C data

MOTIVATION

The organization of the genome into domains plays a central role in gene expression and other cellular activities. Researchers identify genomic domains mainly through two views: 1D functional assays such as ChIP-seq, and chromatin conformation assays such as Hi-C. Fully understanding domains requires integrative modeling that combines these two views. However, the predominant form of integrative modeling uses segmentation and genome annotation (SAGA) along with the rigid assumption that loci in contact are more likely to share the same domain type, which is not necessarily true for epigenomic domain types and genome-wide chromatin interactions.

RESULTS

Here, we present an integrative approach that annotates domains using both 1D functional genomic signals and Hi-C measurements of genome-wide 3D interactions without the use of a pairwise prior. We do so by using a graph embedding to learn structural features corresponding to each genomic region, then inputting learned structural features along with functional genomic signals to a SAGA algorithm. We show that our domain types recapitulate well-known subcompartments with an additional granularity that distinguishes a combination of the spatial and functional states of the genomic regions. In particular, we identified a division of the previously identified A2 subcompartment such that the divided domain types have significantly varying expression levels.

AVAILABILITY AND IMPLEMENTATION

https://github.com/nedashokraneh/IChDA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Long-distance association of topological boundaries through nuclear condensates

The eukaryotic genome is partitioned into distinct topological domains separated by boundary elements. Emerging data support the concept that several well-established nuclear compartments are ribonucleoprotein condensates assembled through the physical process of phase separation. Here, based on our demonstration that chemical disruption of nuclear condensate assembly weakens the insulation properties of a specific subset (∼20%) of topologically associated domain (TAD) boundaries, we report that the disrupted boundaries are characterized by a high level of transcription and striking spatial clustering. These topological boundary regions tend to be spatially associated, even interchromosomally, segregate with nuclear speckles, and harbor a specific subset of "housekeeping" genes widely expressed in diverse cell types. These observations reveal a previously unappreciated mode of genome organization mediated by conserved boundary elements harboring highly and widely expressed transcription units and associated transcriptional condensates.

Chromatin Hubs: A biological and computational outlook

This review discusses our current understanding of chromatin biology and bioinformatics under the unifying concept of “chromatin hubs.” The first part reviews the biology of chromatin hubs, including chromatin–chromatin interaction hubs, chromatin hubs at the nuclear periphery, hubs around macromolecules such as RNA polymerase or lncRNAs, and hubs around nuclear bodies such as the nucleolus or nuclear speckles. The second part reviews existing computational methods, including enhancer–promoter interaction prediction, network analysis, chromatin domain callers, transcription factory predictors, and multi-way interaction analysis. We introduce an integrated model that makes sense of the existing evidence. Understanding chromatin hubs may allow us (i) to explain long-unsolved biological questions such as interaction specificity and redundancy of mechanisms, (ii) to develop more realistic kinetic and functional predictions, and (iii) to explain the etiology of genomic disease.

Egr1 is a 3D matrix-specific mediator of mechanosensitive stem cell lineage commitment

While extracellular matrix (ECM) mechanics strongly regulate stem cell commitment, the field's mechanistic understanding of this phenomenon largely derives from simplified two-dimensional (2D) culture substrates. Here, we found a 3D matrix-specific mechanoresponsive mechanism for neural stem cell (NSC) differentiation. NSC lineage commitment in 3D is maximally stiffness sensitive in the range of 0.1 to 1.2 kPa, a narrower and more brain-mimetic range than we had previously identified in 2D (0.75 to 75 kPa). Transcriptomics revealed stiffness-dependent up-regulation of early growth response 1 (Egr1) in 3D but not in 2D. Egr1 knockdown enhanced neurogenesis in stiff ECMs by driving β-catenin nuclear localization and activity in 3D, but not in 2D. Mechanical modeling and experimental studies under osmotic pressure indicate that stiff 3D ECMs are likely to stimulate Egr1 via increases in confining stress during cell volumetric growth. To our knowledge, Egr1 represents the first 3D-specific stem cell mechanoregulatory factor.

Cytogenomic Characterization of Giant Ring or Rod Marker Chromosome in Four Cases of Well-Differentiated and Dedifferentiated Liposarcoma

Chromosome and array comparative genomic hybridization (aCGH) analyses were performed on two cases of well-differentiated liposarcoma (WDLPS) and two cases of dedifferentiated liposarcoma (DDLPS). The results revealed the characteristic giant ring (GR) or giant rod marker (GRM) chromosomes in all four cases and amplification of numerous somatic copy number alterations (SCNAs) involving a core segment of 12q14.1q15 and other chromosomal regions in three cases. The levels of amplification for oncogenes OS9, CDK4, HMGA2, NUP107, MDM2, YEATS4, and FRS2 at the core segment or other SCNAs should be characterized to facilitate pathologic correlation and prognostic prediction. Further studies for the initial cellular crisis event affecting chromosome intermingling regions for cell-type specific gene regulation may reveal the underlying mutagenesis mechanism for GR and GRM in WDLPS and DDLPS.

Effects of forces on chromatin

Chromatin is a unique structure of DNA and histone proteins in the cell nucleus and the site of dynamic regulation of gene expression. Soluble factors are known to affect the chromatin structure and function via activating or inhibiting specific transcription factors. Forces on chromatin come from exogenous stresses on the cell surface and/or endogenous stresses, which are regulated by substrate mechanics, geometry, and topology. Forces on chromatin involve direct (via adhesion molecules, cytoskeleton, and the linker of nucleoskeleton and cytoskeleton complexes) and indirect (via diffusion and/or translocation processes) signaling pathways to modulate levels of chromatin folding and deformation to regulate transcription, which is controlled by histone modifications and depends on magnitude, direction, rate/frequency, duration, and modes of stresses. The rapid force transmission pathway activates multiple genes simultaneously, and the force may act like a "supertranscription factor." The indirect mechanotransduction pathways and the rapid force transmission pathway together exert sustained impacts on the chromatin, the nucleus, and cell functions.

Opportunity to improve livestock traits using 3D genomics

The advent of high-throughput chromosome conformation capture and sequencing (Hi-C) has enabled researchers to probe the 3D architecture of the mammalian genome in a genome-wide manner. Simultaneously, advances in epigenomic assays, such as chromatin immunoprecipitation and sequencing (ChIP-seq) and DNase-seq, have enabled researchers to study cis-regulatory interactions and chromatin accessibility across the same genome-wide scale. The use of these data has revealed many unique insights into gene regulation and disease pathomechanisms in several model organisms. With the advent of these high-throughput sequencing technologies, there has been an ever-increasing number of datasets available for study; however, this is often limited to model organisms. Livestock species play critical roles in the economies of developing and developed nations alike. Despite this, they are greatly underrepresented in the 3D genomics space; Hi-C and related technologies have the potential to revolutionise livestock breeding by enabling a more comprehensive understanding of how production traits are controlled. The growth in human and model organism Hi-C data has seen a surge in the availability of computational tools for use in 3D genomics, with some tools using machine learning techniques to predict features and improve dataset quality. In this review, we provide an overview of the 3D genome and discuss the status of 3D genomics in livestock before delving into advancing the field by drawing inspiration from research in human and mouse. We end by offering future directions for livestock research in the field of 3D genomics.

Rapid nucleus-scale reorganization of chromatin in neurons enables transcriptional adaptation for memory consolidation

The interphase nucleus is functionally organized in active and repressed territories defining the transcriptional status of the cell. However, it remains poorly understood how the nuclear architecture of neurons adapts in response to behaviorally relevant stimuli that trigger fast alterations in gene expression patterns. Imaging of fluorescently tagged nucleosomes revealed that pharmacological manipulation of neuronal activity in vitro and auditory cued fear conditioning in vivo induce nucleus-scale restructuring of chromatin within minutes. Furthermore, the acquisition of auditory fear memory is impaired after infusion of a drug into auditory cortex which blocks chromatin reorganization in vitro. We propose that active chromatin movements at the nucleus scale act together with local gene-specific modifications to enable transcriptional adaptations at fast time scales. Introducing a transgenic mouse line for photolabeling of histones, we extend the realm of systems available for imaging of chromatin dynamics to living animals.

Reciprocal regulation of cellular mechanics and metabolism

Metabolism and mechanics are intrinsically intertwined. External forces, sensed through the cytoskeleton or distortion of the cell and organelles, induce metabolic changes in the cell. The resulting changes in metabolism, in turn, feed back to regulate every level of cell biology, including the mechanical properties of cells and tissues. Here we examine the links between metabolism and mechanics, highlighting signalling pathways involved in the regulation and response to cellular mechanosensing. We consider how forces and metabolism regulate one another through nanoscale molecular sensors, micrometre-scale cytoskeletal networks, organelles and dynamic biomolecular condensates. Understanding this cross-talk will create diagnostic and therapeutic opportunities for metabolic disorders such as cancer, cardiovascular pathologies and obesity.

Chromatin Structure and Function in Mosquitoes

The principles and function of chromatin and nuclear architecture have been extensively studied in model organisms, such as Drosophila melanogaster. However, little is known about the role of these epigenetic processes in transcriptional regulation in other insects including mosquitoes, which are major disease vectors and a worldwide threat for human health. Some of these life-threatening diseases are malaria, which is caused by protozoan parasites of the genus Plasmodium and transmitted by Anopheles mosquitoes; dengue fever, which is caused by an arbovirus mainly transmitted by Aedes aegypti; and West Nile fever, which is caused by an arbovirus transmitted by Culex spp. In this contribution, we review what is known about chromatin-associated mechanisms and the 3D genome structure in various mosquito vectors, including Anopheles, Aedes, and Culex spp. We also discuss the similarities between epigenetic mechanisms in mosquitoes and the model organism Drosophila melanogaster, and advocate that the field could benefit from the cross-application of state-of-the-art functional genomic technologies that are well-developed in the fruit fly. Uncovering the mosquito regulatory genome can lead to the discovery of unique regulatory networks associated with the parasitic life-style of these insects. It is also critical to understand the molecular interactions between the vectors and the pathogens that they transmit, which could hold the key to major breakthroughs on the fight against mosquito-borne diseases. Finally, it is clear that epigenetic mechanisms controlling mosquito environmental plasticity and evolvability are also of utmost importance, particularly in the current context of globalization and climate change.

Non-Canonical Functions of the Gamma-Tubulin Meshwork in the Regulation of the Nuclear Architecture

Simple Summary

The appearance of a cell is connected to its function. For example, the fusiform of smooth muscle cells is adapted to facilitate muscle contraction, the lobed nucleus in white blood cells assists with the migratory behavior of these immune cells, and the condensed nucleus in sperm aids in their swimming efficiency. Thus, changes in appearance have been used for decades by doctors as a diagnostic method for human cancers. Here, we summarize our knowledge of how a cell maintains the shape of the nuclear compartment. Specifically, we discuss the role of a novel protein meshwork, the gamma-tubulin meshwork, in the regulation of nuclear morphology and as a therapeutic target against cancer.

Abstract

The nuclear architecture describes the organization of the various compartments in the nucleus of eukaryotic cells, where a plethora of processes such as nucleocytoplasmic transport, gene expression, and assembly of ribosomal subunits occur in a dynamic manner. During the different phases of the cell cycle, in post-mitotic cells and after oncogenic transformation, rearrangements of the nuclear architecture take place, and, among other things, these alterations result in reorganization of the chromatin and changes in gene expression. A member of the tubulin family, γtubulin, was first identified as part of a multiprotein complex that allows nucleation of microtubules. However, more than a decade ago, γtubulin was also characterized as a nuclear protein that modulates several crucial processes that affect the architecture of the nucleus. This review presents the latest knowledge regarding changes that arise in the nuclear architecture of healthy cells and under pathological conditions and, more specifically, considers the particular involvement of γtubulin in the modulation of the biology of the nuclear compartment.

Association between fetal sex and maternal plasma microRNA responses to prenatal alcohol exposure: evidence from a birth outcome-stratified cohort

Most persons with fetal alcohol spectrum disorders (FASDs) remain undiagnosed or are diagnosed in later life. To address the need for earlier diagnosis, we previously assessed miRNAs in the blood plasma of pregnant women who were classified as unexposed to alcohol (UE), heavily exposed with affected infants (HEa), or heavily exposed with apparently unaffected infants (HEua). We reported that maternal miRNAs predicted FASD-related growth and psychomotor deficits in infants. Here, we assessed whether fetal sex influenced alterations in maternal circulating miRNAs following prenatal alcohol exposure (PAE). To overcome the loss of statistical power due to disaggregating maternal samples by fetal sex, we adapted a strategy of iterative bootstrap resampling with replacement to assess the stability of statistical parameter estimates. Bootstrap estimates of parametric and effect size tests identified male and female fetal sex-associated maternal miRNA responses to PAE that were not observed in the aggregated sample. Additionally, we observed, in HEa mothers of female, but not male fetuses, a network of co-secreted miRNAs whose expression was linked to miRNAs encoded on the X-chromosome. Interestingly, the number of significant miRNA correlations for the HEua group mothers with female fetuses was intermediate between HEa and UE mothers at mid-pregnancy, but more similar to UE mothers by the end of pregnancy. Collectively, these data show that fetal sex predicts maternal circulating miRNA adaptations, a critical consideration when adopting maternal miRNAs as diagnostic biomarkers. Moreover, a maternal co-secretion network, predominantly in pregnancies with female fetuses, emerged as an index of risk for adverse birth outcomes due to PAE.

Hi-C chromosome conformation capture sequencing of avian genomes using the BGISEQ-500 platform

BACKGROUND

Hi-C experiments couple DNA-DNA proximity with next-generation sequencing to yield an unbiased description of genome-wide interactions. Previous methods describing Hi-C experiments have focused on the industry-standard Illumina sequencing. With new next-generation sequencing platforms such as BGISEQ-500 becoming more widely available, protocol adaptations to fit platform-specific requirements are useful to give increased choice to researchers who routinely generate sequencing data.

RESULTS

We describe an in situ Hi-C protocol adapted to be compatible with the BGISEQ-500 high-throughput sequencing platform. Using zebra finch (Taeniopygia guttata) as a biological sample, we demonstrate how Hi-C libraries can be constructed to generate informative data using the BGISEQ-500 platform, following circularization and DNA nanoball generation. Our protocol is a modification of an Illumina-compatible method, based around blunt-end ligations in library construction, using un-barcoded, distally overhanging double-stranded adapters, followed by amplification using indexed primers. The resulting libraries are ready for circularization and subsequent sequencing on the BGISEQ series of platforms and yield data similar to what can be expected using Illumina-compatible approaches.

CONCLUSIONS

Our straightforward modification to an Illumina-compatible in situHi-C protocol enables data generation on the BGISEQ series of platforms, thus expanding the options available for researchers who wish to utilize the powerful Hi-C techniques in their research.

IGAP-integrative genome analysis pipeline reveals new gene regulatory model associated with nonspecific TF-DNA binding affinity

The human genome is regulated in a multi-dimensional way. While biophysical factors like Non-specific Transcription factor Binding Affinity (nTBA) act at DNA sequence level, other factors act above sequence levels such as histone modifications and 3-D chromosomal interactions. This multidimensionality of regulation requires many of these factors for a proper understanding of the regulatory landscape of the human genome. Here, we propose a new biophysical model for estimating nTBA. Integration of nTBA with chromatin modifications and chromosomal interactions, using a new Integrative Genome Analysis Pipeline (IGAP), reveals additive effects of nTBA to regulatory DNA sequences and identifies three types of genomic zones in the human genome (Inactive Genomic Zones, Poised Genomic Zones, and Active Genomic Zones). It also unveils a novel long distance gene regulatory model: chromosomal interactions reduce the physical distance between the high occupancy target (HOT) regions that results in high nTBA to DNA in the area, which in turn attract TFs to such regions with higher binding potential. These findings will help to elucidate the three-dimensional diffusion process that TFs use during their search for the right targets.

Predicting cell lineages using autoencoders and optimal transport

Lineage tracing involves the identification of all ancestors and descendants of a given cell, and is an important tool for studying biological processes such as development and disease progression. However, in many settings, controlled time-course experiments are not feasible, for example when working with tissue samples from patients. Here we present ImageAEOT, a computational pipeline based on autoencoders and optimal transport for predicting the lineages of cells using time-labeled datasets from different stages of a cellular process. Given a single-cell image from one of the stages, ImageAEOT generates an artificial lineage of this cell based on the population characteristics of the other stages. These lineages can be used to connect subpopulations of cells through the different stages and identify image-based features and biomarkers underlying the biological process. To validate our method, we apply ImageAEOT to a benchmark task based on nuclear and chromatin images during the activation of fibroblasts by tumor cells in engineered 3D tissues. We further validate ImageAEOT on chromatin images of various breast cancer cell lines and human tissue samples, thereby linking alterations in chromatin condensation patterns to different stages of tumor progression. Our results demonstrate the promise of computational methods based on autoencoding and optimal transport principles for lineage tracing in settings where existing experimental strategies cannot be used.

Many key biological processes, such as development and disease progression, require analyzing lineages of cells backwards as well as forwards in time. However, current single-cell experiments tend to be destructive to cells, so that a single lineage can only be measured at one point in time. In this work, we introduce a computational framework for predicting the lineage of cells from a single snapshot in time based on measurements of other cells at other time points. The method generates these lineages by computing the most plausible path for a population of cells to transition from one time point to the next, assuming that a cell is more likely to transition to similar cells compared to dissimilar cells. We validate our computational method on imaging data of fibroblasts and cancer cells, though our method could also be applied to other modalities of single-cell data such as genomics and transcriptomics as well as multi-modal single-cell datasets.

Force-induced gene up-regulation does not follow the weak power law but depends on H3K9 demethylation

Mechanical forces play important roles in development, physiology, and diseases, but how force is transduced into gene transcription remains elusive. Here, we show that transcription of transgene DHFR or endogenous genes egr-1 and Cav1 is rapidly up-regulated in response to cyclic forces applied via integrins at low frequencies but not at 100 Hz. Gene up-regulation does not follow the weak power law with force frequency. Force-induced transcription up-regulation at the nuclear interior is associated with demethylation of histone H3 lysine-9 trimethylation (H3K9me3), whereas no transcription up-regulation near the nuclear periphery is associated with H3K9me3 that inhibits Pol II recruitment to the promoter site. H3K9me3 demethylation induces Pol II recruitment and increases force-induced transcription of egr-1 and Cav1 at the nuclear interior and activates mechano-nonresponsive gene FKBP5 near the nuclear periphery, whereas H3K9me3 hypermethylation has opposite effects. Our findings demonstrate that rapid up-regulation of endogenous mechanoresponsive genes depends on H3K9me3 demethylation.

Multivariate analysis reveals activation-primed fibroblast geometric states in engineered 3D tumor microenvironments

Fibroblasts are a heterogeneous group of cells comprising subpopulations that have been found to be activated in the stromal microenvironment that regulates tumor initiation and growth. The underlying mechanisms of such selective activation of fibroblasts are not understood. We propose that the intrinsic geometric heterogeneity of fibroblasts modulates the nuclear mechanotransduction of signals from the microenvironment, resulting in their selective activation. To test this, we developed an engineered 3D fibroblast tumor coculture system and used high resolution images to quantify multiple cell geometry sensitive nuclear morphological and chromatin organizational features. These features were then mapped to activation levels as measured by the nuclear abundance of transcription cofactor, megakaryoblastic leukemia, and protein levels of its target, αSMA. Importantly, our results indicate the presence of activation-“primed” cell geometries that present higher activation levels, which are further enhanced in the presence of stimuli from cancer cells. Further, we show that by enriching the population of activation-primed cell geometric states by either increasing matrix rigidity or micropatterning primed cell shapes, fibroblast activation levels can be increased. Collectively, our results reveal important cellular geometric states that select for fibroblast activation within the heterogenous tumor microenvironment.

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.

MOCHI enables discovery of heterogeneous interactome modules in 3D nucleome

The composition of the cell nucleus is highly heterogeneous, with different constituents forming complex interactomes. However, the global patterns of these interwoven heterogeneous interactomes remain poorly understood. Here we focus on two different interactomes, chromatin interaction network and gene regulatory network, as a proof of principle to identify heterogeneous interactome modules (HIMs), each of which represents a cluster of gene loci that is in spatial contact more frequently than expected and that is regulated by the same group of transcription factors. HIM integrates transcription factor binding and 3D genome structure to reflect “transcriptional niche” in the nucleus. We develop a new algorithm, MOCHI, to facilitate the discovery of HIMs based on network motif clustering in heterogeneous interactomes. By applying MOCHI to five different cell types, we found that HIMs have strong spatial preference within the nucleus and show distinct functional properties. Through integrative analysis, this work shows the utility of MOCHI to identify HIMs, which may provide new perspectives on the interplay between transcriptional regulation and 3D genome organization.

H3K9me3, H3K36me3, and H4K20me3 Expression Correlates with Patient Outcome in Esophageal Squamous Cell Carcinoma as Epigenetic Markers

BACKGROUND

Histone methylation, as an essential pattern of posttranslational modifications, contributes to multiple cancer-related biological processes. Dysregulation of histone methylation is now considered a biomarker for cancer prognosis.

AIMS

This study investigated and evaluated the potential role of four histone lysine trimethylation markers as biomarkers for esophageal squamous cell carcinoma (ESCC) prognosis.

METHODS

Tissue arrays were made from 135 paraffin-embedded ESCC samples and examined for histone markers by immunohistochemistry, and 10 pairs of cancer and noncancerous mucosa tissues from ESCC patients were investigated with Western blot. Chi-squared test, Kaplan-Meier analysis with log-rank test, and Cox proportional hazard trend analyses were performed to assess the prognostic values of the markers.

RESULTS

Histone 3 lysine 4 trimethylation (H3K4me3), histone 3 lysine 9 trimethylation (H3K9me3), and histone 4 lysine 20 trimethylation (H4K20me3), but not histone 3 lysine 36 trimethylation (H3K36me3), showed stronger immunostaining signals in tumor tissues than in the corresponding adjacent non-neoplastic mucosa tissues. The expression patterns of H3K36me3, H3K9me3, and H4K20me3 correlated with tumor infiltrating depth, lymph node involvement, and pTNM stage. Low-scoring H3K9me3 and H4K20me3 predicted better prognosis, while H3K36me3 manifested the opposite trend. Poor prognosis occurred in ESCC patients with expression patterns of high levels of H3K9me3, high levels of H4K20me3, and low levels of H3K36me3 expression.

CONCLUSIONS

H3K9me3, H4K20me3, and H3K36me3 showed a close relationship with clinical features and were considered independent risk factors for survival of ESCC patients. The combination of H3K9me3, H4K20me3, and H3K36me3 expression, rather than the expression of a single histone marker, is believed to further enhance evaluations of ESCC prognosis and management.

Mapping the spectrum of 3D communities in human chromosome conformation capture data

Several experiments show that the three dimensional (3D) organization of chromosomes affects genetic processes such as transcription and gene regulation. To better understand this connection, researchers developed the Hi-C method that is able to detect the pairwise physical contacts of all chromosomal loci. The Hi-C data show that chromosomes are composed of 3D compartments that range over a variety of scales. However, it is challenging to systematically detect these cross-scale structures. Most studies have therefore designed methods for specific scales to study foremost topologically associated domains (TADs) and A/B compartments. To go beyond this limitation, we tailor a network community detection method that finds communities in compact fractal globule polymer systems. Our method allows us to continuously scan through all scales with a single resolution parameter. We found: (i) polymer segments belonging to the same 3D community do not have to be in consecutive order along the polymer chain. In other words, several TADs may belong to the same 3D community. (ii) CTCF proteins-a loop-stabilizing protein that is ascribed a big role in TAD formation-are well correlated with community borders only at one level of organization. (iii) TADs and A/B compartments are traditionally treated as two weakly related 3D structures and detected with different algorithms. With our method, we detect both by simply adjusting the resolution parameter. We therefore argue that they represent two specific levels of a continuous spectrum 3D communities, rather than seeing them as different structural entities.

Mechanical regulation of genome architecture and cell-fate decisions

Landmark experiments in vitro showed that somatic cells can be reprogrammed to stem-cells by the constitutive expression of particular transcription factors. However, in vivo cells naturally exhibit de-differentiation and trans-differentiation programs, thereby suggesting that the signals from the local mechanical microenvironment may be sufficient to induce stem-cell state transitions. In this review, we discuss recent evidence for the biophysical regulation of genome architecture and nuclear programs. We start by discussing the coupling between cellular architecture, genome organization and gene expression. We then review the role of biophysical factors in regulating genome architecture and cell-state transitions. Finally, we discuss the molecular basis of cell-state transitions during nuclear reprogramming. Collectively, we highlight the importance of the mechanical regulation of genome organization on cell-fate decisions.

Exploring chromatin hierarchical organization via Markov State Modelling

We propose a new computational method for exploring chromatin structural organization based on Markov State Modelling of Hi-C data represented as an interaction network between genomic loci. A Markov process describes the random walk of a traveling probe in the corresponding energy landscape, mimicking the motion of a biomolecule involved in chromatin function. By studying the metastability of the associated Markov State Model upon annealing, the hierarchical structure of individual chromosomes is observed, and corresponding set of structural partitions is identified at each level of hierarchy. Then, the notion of effective interaction between partitions is derived, delineating the overall topology and architecture of chromosomes. Mapping epigenetic data on the graphs of intra-chromosomal effective interactions helps in understanding how chromosome organization facilitates its function. A sketch of whole-genome interactions obtained from the analysis of 539 partitions from all 23 chromosomes, complemented by distributions of gene expression regulators and epigenetic factors, sheds light on the structure-function relationships in chromatin, delineating chromosomal territories, as well as structural partitions analogous to topologically associating domains and active / passive epigenomic compartments. In addition to the overall genome architecture shown by effective interactions, the affinity between partitions of different chromosomes was analyzed as an indicator of the degree of association between partitions in functionally relevant genomic interactions. The overall static picture of whole-genome interactions obtained with the method presented in this work provides a foundation for chromatin structural reconstruction, for the modelling of chromatin dynamics, and for exploring the regulation of genome function. The algorithms used in this study are implemented in a freely available Python package ChromaWalker (https://bitbucket.org/ZhenWahTan/chromawalker).

Automatic detection of genomic regions with informative epigenetic patterns

BACKGROUND

Epigenetic phenomena are crucial for explaining the phenotypic plasticity seen in the cells of different tissues, developmental stages and diseases, all holding the same DNA sequence. As technology is allowing to retrieve epigenetic information in a genome-wide fashion, massive epigenomic datasets are being accumulated in public repositories. New approaches are required to mine those data to extract useful knowledge. We present here an automatic approach for detecting genomic regions with epigenetic variation patterns across samples related to a grouping of these samples, as a way of detecting regions functionally associated to the phenomenon behind the classification.

RESULTS

We show that the regions automatically detected by the method in the whole human genome associated to three different classifications of a set of epigenomes (cancer vs. healthy, brain vs. other organs, and fetal vs. adult tissues) are enriched in genes associated to these processes.

CONCLUSIONS

The method is fully automatic and can exhaustively scan the whole human genome at any resolution using large collections of epigenomes as input, although it also produces good results with small datasets. Consequently, it will be valuable for obtaining functional information from the incoming epigenomic information as it continues to accumulate.

Modeling three-dimensional genomic organization in evolution and pathogenesis

The regulation of gene expression is mediated via the complex three-dimensional (3D) conformation of the genetic material and its interactions with various intracellular factors. Various experimental and computational approaches have been developed in recent years for understating the relation between the 3D conformation of the genome and the phenotypes of cells in normal condition and diseases. In this review, we will discuss novel approaches for analyzing and modeling the 3D genomic conformation, focusing on deciphering disease-causing mutations that affect gene expression. We conclude that as this is a very challenging mission, an important direction should involve the comparative analysis of various 3D models from various organisms or cells.

The 3D genome: From fundamental principles to disease and cancer

In higher eukaryotes, the three-dimensional (3D) organization of the genome is intimately related to numerous key biological functions including gene expression, DNA repair and DNA replication regulations. Alteration of this 3D organization is detrimental to the organism and can give rise to a broad range of diseases such as cancers. Here, we review recent advances in the field. We first describe how the genome is packed in 3D to form chromosome territories, compartments and domains. We also give an overview of the recent techniques that allow to map the genome in 3D up to the kilobase resolution. We then discuss potential mechanisms by which genome misfolding can affect proper gene expression by distal enhancers, and how the 3D genome influences the formation of genomic rearrangements.

Forces driving the three-dimensional folding of eukaryotic genomes

The last decade has radically renewed our understanding of higher order chromatin folding in the eukaryotic nucleus. As a result, most current models are in support of a mostly hierarchical and relatively stable folding of chromosomes dividing chromosomal territories into A‐ (active) and B‐ (inactive) compartments, which are then further partitioned into topologically associating domains (TADs), each of which is made up from multiple loops stabilized mainly by the CTCF and cohesin chromatin‐binding complexes. Nonetheless, the structure‐to‐function relationship of eukaryotic genomes is still not well understood. Here, we focus on recent work highlighting the biophysical and regulatory forces that contribute to the spatial organization of genomes, and we propose that the various conformations that chromatin assumes are not so much the result of a linear hierarchy, but rather of both converging and conflicting dynamic forces that act on it.

Nuclear Mechanopathology and Cancer Diagnosis

Abnormalities in nuclear and chromatin organization are hallmarks of many diseases including cancer. In this review, we highlight our understanding of how the cellular microenvironment regulates nuclear morphology and, with it, the spatial organization of chromosomes and genes, resulting in cell type-specific genomic programs. We also discuss the molecular basis for maintaining nuclear and genomic integrity and how alterations in nuclear mechanotransduction pathways result in various diseases. Finally, we highlight the importance of digital pathology based on nuclear morphometric features combined with single-cell genomics for early cancer diagnostics.