Deciphering multi-way interactions in the human genome.

Chromatin architecture, a key regulator of gene expression, can be inferred using chromatin contact data from chromosome conformation capture, or Hi-C. However, classical Hi-C does not preserve multi-way contacts. Here we use long sequencing reads to map genome-wide multi-way contacts and investigate higher order chromatin organization in the human genome. We use hypergraph theory for data representation and analysis, and quantify higher order structures in neonatal fibroblasts, biopsied adult fibroblasts, and B lymphocytes. By integrating multi-way contacts with chromatin accessibility, gene expression, and transcription factor binding, we introduce a data-driven method to identify cell type-specific transcription clusters. We provide transcription factor-mediated functional building blocks for cell identity that serve as a global signature for cell types.

Rearrangement of T Cell genome architecture regulates GVHD.

WAPL, cohesin's DNA release factor, regulates three-dimensional (3D) chromatin architecture. The 3D chromatin structure and its relevance to mature T cell functions is not well understood. We show that in vivo lymphopenic expansion, and alloantigen-driven proliferation, alters the 3D structure and function of the genome in mature T cells. Conditional deletion of WAPL, cohesin's DNA release factor, in T cells reduced long-range genomic interactions and altered chromatin A/B compartments and interactions within topologically associating domains (TADs) of the chromatin in T cells at baseline. WAPL deficiency in T cells reduced loop extensions, changed expression of cell cycling genes and reduced proliferation following in vitro and in vivo stimulation, and reduced severity of graft-versus-host disease (GVHD) following experimental allogeneic hematopoietic stem cell transplantation. These data collectively characterize 3D genomic architecture of T cells in vivo and demonstrate biological and clinical implications for its disruption by cohesin release factor WAPL.

Valproic acid-induced changes of 4D nuclear morphology in astrocyte cells.

Histone deacetylase inhibitors, such as valproic acid (VPA), have important clinical therapeutic and cellular reprogramming applications. They induce chromatin reorganization that is associated with altered cellular morphology. However, there is a lack of comprehensive characterization of VPA-induced changes of nuclear size and shape. Here, we quantify 3D nuclear morphology of primary human astrocyte cells treated with VPA over time (hence, 4D). We compared volumetric and surface-based representations and identified seven features that jointly discriminate between normal and treated cells with 85% accuracy on day 7. From day 3, treated nuclei were more elongated and flattened and then continued to morphologically diverge from controls over time, becoming larger and more irregular. On day 7, most of the size and shape descriptors demonstrated significant differences between treated and untreated cells, including a 24% increase in volume and 6% reduction in extent (shape regularity) for treated nuclei. Overall, we show that 4D morphometry can capture how chromatin reorganization modulates the size and shape of the nucleus over time. These nuclear structural alterations may serve as a biomarker for histone (de-)acetylation events and provide insights into mechanisms of astrocytes-to-neurons reprogramming.

Publisher Correction: 3D Shape Modeling for Cell Nuclear Morphological Analysis and Classification.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Hypothesis: Caco-2 cell rotational 3D mechanogenomic turing patterns have clinical implications to colon crypts.

Colon crypts are recognized as a mechanical and biochemical Turing patterning model. Colon epithelial Caco-2 cell monolayer demonstrated 2D Turing patterns via force analysis of apical tight junction live cell imaging which illuminated actomyosin meshwork linking the actomyosin network of individual cells. Actomyosin forces act in a mechanobiological manner that alters cell/nucleus/tissue morphology. We observed the rotational motion of the nucleus in Caco-2 cells that appears to be driven by actomyosin during the formation of a differentiated confluent epithelium. Single- to multi-cell ring/torus-shaped genomes were observed prior to complex fractal Turing patterns extending from a rotating torus centre in a spiral pattern consistent with a gene morphogen motif. These features may contribute to the well-described differentiation from stem cells at the crypt base to the luminal colon epithelium along the crypt axis. This observation may be useful to study the role of mechanogenomic processes and the underlying molecular mechanisms as determinants of cellular and tissue architecture in space and time, which is the focal point of the 4D nucleome initiative. Mathematical and bioengineer modelling of gene circuits and cell shapes may provide a powerful algorithm that will contribute to future precision medicine relevant to a number of common medical disorders.

Genome Architecture Mediates Transcriptional Control of Human Myogenic Reprogramming.

Genome architecture has emerged as a critical element of transcriptional regulation, although its role in the control of cell identity is not well understood. Here we use transcription factor (TF)-mediated reprogramming to examine the interplay between genome architecture and transcriptional programs that transition cells into the myogenic identity. We recently developed new methods for evaluating the topological features of genome architecture based on network centrality. Through integrated analysis of these features of genome architecture and transcriptome dynamics during myogenic reprogramming of human fibroblasts we find that significant architectural reorganization precedes activation of a myogenic transcriptional program. This interplay sets the stage for a critical transition observed at several genomic scales reflecting definitive adoption of the myogenic phenotype. Subsequently, TFs within the myogenic transcriptional program participate in entrainment of biological rhythms. These findings reveal a role for topological features of genome architecture in the initiation of transcriptional programs during TF-mediated human cellular reprogramming.

3D Shape Modeling for Cell Nuclear Morphological Analysis and Classification.

Quantitative analysis of morphological changes in a cell nucleus is important for the understanding of nuclear architecture and its relationship with pathological conditions such as cancer. However, dimensionality of imaging data, together with a great variability of nuclear shapes, presents challenges for 3D morphological analysis. Thus, there is a compelling need for robust 3D nuclear morphometric techniques to carry out population-wide analysis. We propose a new approach that combines modeling, analysis, and interpretation of morphometric characteristics of cell nuclei and nucleoli in 3D. We used robust surface reconstruction that allows accurate approximation of 3D object boundary. Then, we computed geometric morphological measures characterizing the form of cell nuclei and nucleoli. Using these features, we compared over 450 nuclei with about 1,000 nucleoli of epithelial and mesenchymal prostate cancer cells, as well as 1,000 nuclei with over 2,000 nucleoli from serum-starved and proliferating fibroblast cells. Classification of sets of 9 and 15 cells achieved accuracy of 95.4% and 98%, respectively, for prostate cancer cells, and 95% and 98% for fibroblast cells. To our knowledge, this is the first attempt to combine these methods for 3D nuclear shape modeling and morphometry into a highly parallel pipeline workflow for morphometric analysis of thousands of nuclei and nucleoli in 3D.

Chronic stress and intestinal barrier dysfunction: Glucocorticoid receptor and transcription repressor HES1 regulate tight junction protein Claudin-1 promoter.

Chronic stress and elevated glucocorticoid hormone are associated with decreases in the intestinal epithelial tight junction protein claudin-1 (CLDN1). Human/rat CLDN1 promoters contain glucocorticoid response elements (GREs) and adjacent transcription repressor HES1 binding N-boxes. Notch signaling target HES1 expression was high and glucocorticoid receptor (NR3C1) low at the crypt base and the pattern reversed at the crypt apex. Chronic stress reduced overall rat colon HES1 and NR3C1 that was associated with CLDN1 downregulation. Chromatin-immunoprecipitation experiments showed that HES1 and NR3C1 bind to the CLDN1 promoter in rat colon crypts. The binding of NR3C1 but not HES1 to CLDN1 promoter significantly decreased in chronically stressed animals, which was prevented by the NR3C1 antagonist RU486. We employed the 21-day Caco-2/BBe cell model to replicate cell differentiation along the crypt axis. HES1 siRNA treatment early in differentiation increased CLDN1. In contrast, stress levels of cortisol decreased CLDN1 in late differentiation stage but not in the early stage. HES1 was high, whereas NR3C1 and CLDN1 were low in the early stage which reversed in the late stage, e.g. HES1/NR3C1 binding to CLDN1 promoter demonstrates a dynamic and reciprocal pattern. These results suggest that chronic stress impairs colon epithelium homeostasis and barrier function via different mechanisms along the crypt axis.

Insight into dynamic genome imaging: Canonical framework identification and high-throughput analysis.

The human genome is dynamic in structure, complicating researcher's attempts at fully understanding it. Time series "Fluorescent in situ Hybridization" (FISH) imaging has increased our ability to observe genome structure, but due to cell type and experimental variability this data is often noisy and difficult to analyze. Furthermore, computational analysis techniques are needed for homolog discrimination and canonical framework detection, in the case of time-series images. In this paper we introduce novel ideas for nucleus imaging analysis, present findings extracted using dynamic genome imaging, and propose an objective algorithm for high-throughput, time-series FISH imaging. While a canonical framework could not be detected beyond statistical significance in the analyzed dataset, a mathematical framework for detection has been outlined with extension to 3D image analysis.

Functional organization of the human 4D Nucleome.

The 4D organization of the interphase nucleus, or the 4D Nucleome (4DN), reflects a dynamical interaction between 3D genome structure and function and its relationship to phenotype. We present initial analyses of the human 4DN, capturing genome-wide structure using chromosome conformation capture and 3D imaging, and function using RNA-sequencing. We introduce a quantitative index that measures underlying topological stability of a genomic region. Our results show that structural features of genomic regions correlate with function with surprising persistence over time. Furthermore, constructing genome-wide gene-level contact maps aided in identifying gene pairs with high potential for coregulation and colocalization in a manner consistent with expression via transcription factories. We additionally use 2D phase planes to visualize patterns in 4DN data. Finally, we evaluated gene pairs within a circadian gene module using 3D imaging, and found periodicity in the movement of clock circadian regulator and period circadian clock 2 relative to each other that followed a circadian rhythm and entrained with their expression.

Periodicity of nuclear morphology in human fibroblasts.

MOTIVATION: Morphology of the cell nucleus has been used as a key indicator of disease state and prognosis, but typically without quantitative rigor. It is also not well understood how nuclear morphology varies with time across different genetic backgrounds in healthy cells. To help answer these questions we measured the size and shape of nuclei in cell-cycle-synchronized primary human fibroblasts from 6 different individuals at 32 time points over a 75 hour period. RESULTS: The nucleus was modeled as an ellipsoid and its dynamics analyzed. Shape and volume changed significantly over this time. Two prominent frequencies were found in the 6 individuals: a 17 hour period consistent with the cell cycle and a 26 hour period. Our findings suggest that the shape of the nucleus changes over time and thus any time-invariant shape property may provide a misleading characterization of cellular populations at different phases of the cell cycle. The proposed methodology provides a general method to analyze morphological change using multiple time points even for non-live-cell experiments.

The cytocaud: a hair cell pathology in the waltzing Guinea pig.

The waltzing guinea pig displays severe inner ear dysfunction that involves both an auditory and a vestibular manifestation. The aim of this study was to characterize a pathological tail-like extension of the vestibular hair cells, the cytocaud. Our data suggest that nearly all type I hair cells in the waltzing guinea pig have cytocauds, which appear as membrane-bound tails containing mitochondria and cytoplasm that proceed in a basal direction toward the basement membrane. The extensions either attach to the basement membrane or penetrate it, and further proceed into the extracellular matrix. A core made of a thick and long (30 microm) actin-rich structure supports the slender long process. The actin core has cross-links that are periodically placed along the length of the cytocaud. Our data suggest that the cytocauds in vestibular hair cells of the waltzing guinea pig are highly organized structures associated with a failure to detach from the basement membrane.