Search

We present HiGlass, an open source visualization tool built on web technologies that provides a rich interface for rapid, multiplex, and multiscale navigation of 2D genomic maps alongside 1D genomic tracks, allowing users to combine various data types, synchronize multiple visualization modalities, and share fully customizable views with others. We demonstrate its utility in exploring different experimental conditions, comparing the results of analyses, and creating interactive snapshots to share with collaborators and the broader public. HiGlass is accessible online at http://higlass.io and is also available as a containerized application that can be run on any platform., National Institutes of Health (U.S.) (U01 CA200059), National Institutes of Health (U.S.) (R00 HG007583), National Institutes of Health (U.S.) (U54 HG007963)

During mitosis, human chromosomes are linearly compacted about 1000-fold by loop-extruding motors. Recent experiments have shown that condensins extrude DNA loops but in a “one-sided” manner. This contrasts with existing models, which predict that symmetric, “two-sided” loop extrusion accounts for mitotic chromosome compaction. We explore whether one-sided extrusion, as it is currently seen in experiments, can compact chromosomes by developing a mean-field theoretical model for polymer compaction by motors that actively extrude loops and dynamically turnover. The model establishes a stringent upper bound of only about tenfold for compaction by strictly one-sided extrusion. We confirm this result with stochastic simulations. Thus, strictly one-sided extrusion as it has been observed so far cannot be the sole mechanism of chromosome compaction. However, as shown by the model, other two-sided or effectively two-sided mechanisms can achieve sufficient compaction.

Figures S1.â S7. Table S1. Supplementary material, and Supplementary references. (PDF 755 kb)

Mitotic chromosomes are among the most recognizable structures in the cell, yet for over a century their internal organization remains largely unsolved. We applied chromosome conformation capture methods, 5C and Hi-C, across the cell cycle and revealed two alternative three-dimensional folding states of the human genome. We show that the highly compartmentalized and cell-type-specific organization described previously for non-synchronous cells is restricted to interphase. In metaphase, we identify a homogenous folding state, which is locus-independent, common to all chromosomes, and consistent among cell types, suggesting a general principle of metaphase chromosome organization. Using polymer simulations, we find that metaphase Hi-C data are inconsistent with classic hierarchical models, and is instead best described by a linearly-organized longitudinally compressed array of consecutive chromatin loops., National Cancer Institute (U.S.) (Grant U54CA143874)

In this Communication we present statistical analysis of conservation profiles in families of homologous sequences for nine proteins whose folding nucleus was determined by protein engineering methods. We show that in all but one protein (AcP) folding nucleus residues are significantly more conserved than the rest of the protein. Two aspects of our study are especially important: 1) grouping of amino acids into classes according to their physical-chemical properties and 2) proper normalization of amino acid probabilities that reflects the fact that evolutionary pressure to conserve some amino acid types may itself affect concentration of various amino acid types in protein families. Neglect of any of those two factors may make physical and biological ``signals'' from conservation profiles disappear., Comment: submitted to Journal of Molecular Biology

Hematopoietic stem cells (HSCs) sustain long-term reconstitution of hematopoiesis in transplantation recipients, yet their role in the endogenous steady-state hematopoiesis remains unclear. In particular, recent studies suggested that HSCs provide a relatively minor contribution to immune cell development in adults. We directed transgene expression in a fraction of HSCs that maintained reconstituting activity during serial transplantations. Inducible genetic labeling showed that transgene-expressing HSCs gave rise to other phenotypic HSCs, confirming their top position in the differentiation hierarchy. The labeled HSCs rapidly contributed to committed progenitors of all lineages and to mature myeloid cells and lymphocytes, but not to B-1a cells or tissue macrophages. Importantly, labeled HSCs gave rise to more than two-thirds of all myeloid cells and platelets in adult mice, and this contribution could be accelerated by an induced interferon response. Thus, classically defined HSCs maintain immune cell development in the steady state and during systemic cytokine responses.

Proper expression of genes requires communication with their regulatory elements that can be located elsewhere along the chromosome. The physics of chromatin fibers imposes a range of constraints on such communication. The molecular and biophysical mechanisms by which chromosomal communication is established, or prevented, have become a topic of intense study, and important roles for the spatial organization of chromosomes are being discovered. Here we present a view of the interphase 3D genome characterized by extensive physical compartmentalization and insulation on the one hand and facilitated long-range interactions on the other. We propose the existence of topological machines dedicated to set up and to exploit a 3D genome organization to both promote and censor communication along and between chromosomes., National Human Genome Research Institute (U.S.) (Grant R01 HG003143), National Human Genome Research Institute (U.S.) (Grant U54 HG007010), National Human Genome Research Institute (U.S.) (Grant U01 HG007910), National Cancer Institute (U.S.) (Grant U54 CA193419), National Institutes of Health (U.S.) (Grant U54 DK107980), National Institutes of Health (U.S.) (Grant U01 DA 040588), National Institute of General Medical Sciences (U.S.) (Grant R01 GM 112720), National Institute of Allergy and Infectious Diseases (U.S.) (Grant U01 R01 AI 117839)

Mammalian interphase chromosomes interact with the nuclear lamina (NL) through hundreds of large lamina-associated domains (LADs). We report a method to map NL contacts genome-wide in single human cells. Analysis of nearly 400 maps reveals a core architecture consisting of gene-poor LADs that contact the NL with high cell-to-cell consistency, interspersed by LADs with more variable NL interactions. The variable contacts tend to be cell-type specific and are more sensitive to changes in genome ploidy than the consistent contacts. Single-cell maps indicate that NL contacts involve multivalent interactions over hundreds of kilobases. Moreover, we observe extensive intra-chromosomal coordination of NL contacts, even over tens of megabases. Such coordinated loci exhibit preferential interactions as detected by Hi-C. Finally, the consistency of NL contacts is inversely linked to gene activity in single cells and correlates positively with the heterochromatic histone modification H3K9me3. These results highlight fundamental principles of single-cell chromatin organization., National Institutes of Health (U.S.) (Grant R01 GM114190), National Human Genome Research Institute (U.S.) (Grant R01 HG003143)

Predicting and constraining RNA virus evolution require understanding the molecular factors that define the mutational landscape accessible to these pathogens. RNA viruses typically have high mutation rates, resulting in frequent production of protein variants with compromised biophysical properties. Their evolution is necessarily constrained by the consequent challenge to protein folding and function. We hypothesized that host proteostasis mechanisms may be significant determinants of the fitness of viral protein variants, serving as a critical force shaping viral evolution. Here, we test that hypothesis by propagating influenza in host cells displaying chemically-controlled, divergent proteostasis environments. We find that both the nature of selection on the influenza genome and the accessibility of specific mutational trajectories are significantly impacted by host proteostasis. These findings provide new insights into features of host-pathogen interactions that shape viral evolution, and into the potential design of host proteostasis-targeted antiviral therapeutics that are refractory to resistance., National Institutes of Health (U.S.) (Award 1DP2GM119162), National Institutes of Health (U.S.) (Grant P30-ES002109)

DNA instability is an important contributor to cancer development. Previously, defects in the chromosome segregation and excessive DNA double strand breaks due to the replication or oxidative stresses were implicated in DNA instability in cancer. Here, we demonstrate that DNA instability can directly result from the oncogene-induced senescence signaling. Expression of the activated form of Her2 oncogene, NeuT, in immortalized breast epithelial cells led to downregulation of the major DNA repair factor histone H2AX and a number of other components of the HR and NHEJ double strand DNA breaks repair pathways. H2AX expression was regulated at the transcriptional level via a senescence pathway involving p21-mediated regulation of CDK and Rb1. The p21-dependent downregulation of H2AX was seen both in cell culture and the MMTV-neu mouse model of Her2-positive breast cancer. Importantly, downregulation of H2AX upon Her2/NeuT expression impaired repair of double strand DNA breaks. This impairment resulted in both increased DNA instability in the form of somatic copy number alterations, and in increased sensitivity to the chemotherapeutic drug doxorubicin. Overall, these findings indicate that the Her2/NeuT oncogene signaling directly potentiates DNA instability and increases sensitivity to DNA damaging treatments, National Institutes of Health (U.S.) (NIH CA081244), National Institutes of Health (U.S.) (NIH R21CA178563), National Institutes of Health (U.S.) (NIH 5U54CA143874)

SummaryTumor-propagating cells in acute leukemia maintain a stem/progenitor-like immature phenotype and proliferative capacity. Acute myeloid leukemia (AML) and acute T-lymphoblastic leukemia (T-ALL) originate from different lineages through distinct oncogenic events such as MLL fusions and Notch signaling, respectively. We found that Zfx, a transcription factor that controls hematopoietic stem cell self-renewal, controls the initiation and maintenance of AML caused by MLL-AF9 fusion and of T-ALL caused by Notch1 activation. In both leukemia types, Zfx prevents differentiation and activates gene sets characteristic of immature cells of the respective lineages. In addition, endogenous Zfx contributes to gene induction and transformation by Myc overexpression in myeloid progenitors. Key Zfx target genes include the mitochondrial enzymes Ptpmt1 and Idh2, whose overexpression partially rescues the propagation of Zfx-deficient AML. These results show that distinct leukemia types maintain their undifferentiated phenotype and self-renewal by exploiting a common stem-cell-related genetic regulator.

The cell fate decision between interferon-producing plasmacytoid DC (pDC) and antigen-presenting classical DC (cDC) is controlled by the E protein transcription factor TCF4 (E2-2). We report that TCF4 comprises two transcriptional isoforms, both of which are required for optimal pDC development in vitro. The long Tcf4 isoform is expressed specifically in pDCs, and its deletion in mice impaired pDCs development and led to the expansion of non-canonical CD8+ cDCs. The expression of Tcf4 commenced in progenitors and was further upregulated in pDCs, correlating with stage-specific activity of multiple enhancer elements. A conserved enhancer downstream of Tcf4 was required for its upregulation during pDC differentiation, revealing a positive feedback loop. The expression of Tcf4 and the resulting pDC differentiation were selectively sensitive to the inhibition of enhancer-binding BET protein activity. Thus, lineage-specifying function of E proteins is facilitated by lineage-specific isoform expression and by BET-dependent feedback regulation through distal regulatory elements.

During cell division, chromosomes are compacted in length by more than a 100-fold. A wide range of experiments demonstrated that in their compacted state, mammalian chromosomes form arrays of closely stacked consecutive ∼100 kb loops. The mechanism underlying the active process of chromosome compaction into a stack of loops is unknown. Here we test the hypothesis that chromosomes are compacted by enzymatic machines that actively extrude chromatin loops. When such loop-extruding factors (LEF) bind to chromosomes, they p rogressively bridge sites that are further away along the chromosome, thus extruding a loop. We demonstrate that collective action of LEFs leads to formation of a dynamic array of consecutive loops. Simulations and an analytically solved model identify two distinct steady states: a sparse state, where loops are highly dynamic but provide little compaction; and a dense state, where there are more stable loops and dramatic chromosome compaction. We find that human chromosomes operate at the border of the dense steady state. Our analysis also shows how the macroscopic characteristics of the loop array are determined by the microscopic properties of LEFs and their abundance. When the number of LEFs are used that match experimentally based estimates, the model can quantitatively reproduce the average loop length, the degree of compaction, and the general loop-array morphology of compact human chromosomes. Our study demonstrates that efficient chromosome compaction can be achieved solely by an active loop-extrusion process., National Institutes of Health (U.S.) (Grant GM114190), National Institutes of Health (U.S.) (Grant R01HG003143)

Peptide-recognition modules (PRMs) are used throughout biology to mediate protein–protein interactions, and many PRMs are members of large protein domain families. Recent genome-wide measurements describe networks of peptide–PRM interactions. In these networks, very similar PRMs recognize distinct sets of peptides, raising the question of how peptide-recognition specificity is achieved using similar protein domains. The analysis of individual protein complex structures often gives answers that are not easily applicable to other members of the same PRM family. Bioinformatics-based approaches, one the other hand, may be difficult to interpret physically. Here we integrate structural information with a large, quantitative data set of SH2 domain–peptide interactions to study the physical origin of domain–peptide specificity. We develop an energy model, inspired by protein folding, based on interactions between the amino-acid positions in the domain and peptide. We use this model to successfully predict which SH2 domains and peptides interact and uncover the positions in each that are important for specificity. The energy model is general enough that it can be applied to other members of the SH2 family or to new peptides, and the cross-validation results suggest that these energy calculations will be useful for predicting binding interactions. It can also be adapted to study other PRM families, predict optimal peptides for a given SH2 domain, or study other biological interactions, e.g. protein–DNA interactions., National Institutes of Health. National Centers for Biomedical Computing (Informatics for Integrating Biology and the Bedside), National Institutes of Health (U.S.) (grant U54LM008748)

Metazoan genomes are spatially organized at multiple scales, from packaging of DNA around individual nucleosomes to segregation of whole chromosomes into distinct territories1–5. At the intermediate scale of kilobases to megabases, which encompasses the sizes of genes, gene clusters and regulatory domains, the three-dimensional (3D) organization of DNA is implicated in multiple gene regulatory mechanisms2–4,6–8, but understanding this organization remains a challenge. At this scale, the genome is partitioned into domains of different epigenetic states that are essential for regulating gene expression9–11. Here, we investigate the 3D organization of chromatin in different epigenetic states using super-resolution imaging. We classified genomic domains in Drosophila cells into transcriptionally active, inactive, or Polycomb-repressed states and observed distinct chromatin organizations for each state. Remarkably, all three types of chromatin domains exhibit power-law scaling between their physical sizes in 3D and their domain lengths, but each type has a distinct scaling exponent. Polycomb-repressed chromatin shows the densest packing and most intriguing folding behaviour in which packing density increases with domain length. Distinct from the self-similar organization displayed by transcriptionally active and inactive chromatin, the Polycomb-repressed domains are characterized by a high degree of chromatin intermixing within the domain. Moreover, compared to inactive domains, Polycomb-repressed domains spatially exclude neighbouring active chromatin to a much stronger degree. Computational modelling and knockdown experiments suggest that reversible chromatin interactions mediated by Polycomb-group proteins plays an important role in these unique packaging properties of the repressed chromatin. Taken together, our super-resolution images reveal distinct chromatin packaging for different epigenetic states at the kilobase-to-megabase scale, a length scale that is directly relevant to genome regulation.

Topological constraints can affect both equilibrium and dynamic properties of polymer systems and can play a role in the organization of chromosomes. Despite many theoretical studies, the effects of topological constraints on the equilibrium state of a single compact polymer have not been systematically studied. Here we use simulations to address this longstanding problem. We find that sufficiently long unknotted polymers differ from knotted ones in the spatial and topological states of their subchains. The unknotted globule has subchains that are mostly unknotted and form asymptotically compact RG(s) ∼ s1/3 crumples. However, crumples display a high fractal dimension of the surface db = 2.8, forming excessive contacts and interpenetrating each other. We conclude that this topologically constrained equilibrium state resembles a conjectured crumpled globule [Grosberg et al., Journal de Physique, 1988, 49, 2095], but differs from its idealized hierarchy of self-similar, isolated and compact crumples., MIT-France Seed Fund, National Cancer Institute (U.S.). Physical Sciences-Oncology Center (MIT, (U54CA143874)), National Research University Higher School of Economics (Program for Basic Research)

How DNA is organized in three dimensions inside the cell nucleus and how this affects the ways in which cells access, read and interpret genetic information are among the longest standing questions in cell biology. Using newly developed molecular, genomic and computational approaches based on the chromosome conformation capture technology (such as 3C, 4C, 5C and Hi-C), the spatial organization of genomes is being explored at unprecedented resolution. Interpreting the increasingly large chromatin interaction data sets is now posing novel challenges. Here we describe several types of statistical and computational approaches that have recently been developed to analyse chromatin interaction data., National Institutes of Health (U.S.), National Human Genome Research Institute (U.S.) (HG003143), National Human Genome Research Institute (U.S.) (HG003143-06S1), W. M. Keck Foundation, Spain. Ministerio de Ciencia e Innovación (BFU2010-19310/BMC), Human Frontier Science Program (Strasbourg, France) (RGP0044/2011), European Union (BLUEPRINT project (EU FP7 grant agreement 282510)), National Cancer Institute (U.S.) (Physical Sciences Oncology Center at MIT, U54CA143874)

Three-dimensional (3D) organization of genomes affects critical cellular processes such as transcription, replication, and deoxyribo nucleic acid (DNA) repair. While previous studies have investigated the natural role, the 3D organization plays in limiting a possible set of genomic rearrangements following DNA repair, the influence of specific organizational principles on this process, particularly over longer evolutionary time scales, remains relatively unexplored. In budding yeast S.cerevisiae, chromosomes are organized into a Rabl-like configuration, with clustered centromeres and telomeres tethered to the nuclear periphery. Hi-C data for S.cerevisiae show that a consequence of this Rabl-like organization is that regions equally distant from centromeres are more frequently in contact with each other, between arms of both the same and different chromosomes. Here, we detect rearrangement events in Saccharomyces species using an automatic approach, and observe increased rearrangement frequency between regions with higher contact frequencies. Together, our results underscore how specific principles of 3D chromosomal organization can influence evolutionary events., National Institutes of Health (U.S.) (Grant GM114190)

The binding of a transcription factor (TF) to a DNA operator site can initiate or repress the expression of a gene. Computational prediction of sites recognized by a TF has traditionally relied upon knowledge of several cognate sites, rather than an ab initio approach. Here, we examine the possibility of using structure-based energy calculations that require no knowledge of bound sites but rather start with the structure of a protein–DNA complex. We study the PurR Escherichia coli TF, and explore to which extent atomistic models of protein–DNA complexes can be used to distinguish between cognate and noncognate DNA sites. Particular emphasis is placed on systematic evaluation of this approach by comparing its performance with bioinformatic methods, by testing it against random decoys and sites of homologous TFs. We also examine a set of experimental mutations in both DNA and the protein. Using our explicit estimates of energy, we show that the specificity for PurR is dominated by direct protein–DNA interactions, and weakly influenced by bending of DNA., National Science Foundation (U.S.) (Grant DMR-08- 03315), Deutsche Forschungsgemeinschaft (DFG) (Grant VI237/1), NEC Research Support Fund, National Institutes of Health. National Centers for Biomedical Computing (Informatics for Integrating Biology and the Bedside), National Institutes of Health (U.S.) (3U54LM008748-04S1)

Strong experimental and theoretical evidence shows that transcription factors and other specific DNA-binding proteins find their sites using a two-mode search: alternating between 3D diffusion through the cell and 1D sliding along the DNA. We consider the role spatial effects in the mechanism on two different scales. First, we reconcile recent experimental findings by showing that the 3D diffusion of the transcription factor is often local, i.e. the transcription factor lands quite near its dissociation site. Second, we discriminate between two types of searches: global searches and local searches. We show that these searches differ significantly in average search time and the variability of search time. Using experimentally measured parameter values, we also show that 1D and 3D search is not optimally balanced, leading to much larger estimates of search time. Together, these results lead to a number of biological implications including suggestions of how prokaryotes and eukaryotes achieve rapid gene regulation and the relationship between the search mechanism and noise in gene expression., Comment: 16 pages, 4 figures

The classic model of eukaryotic gene expression requires direct spatial contact between a distal enhancer and a proximal promoter. Recent Chromosome Conformation Capture (3C) studies show that enhancers and promoters are embedded in a complex network of looping interactions. Here we use a polymer model of chromatin fiber to investigate whether, and to what extent, looping interactions between elements in the vicinity of an enhancer-promoter pair can influence their contact frequency. Our equilibrium polymer simulations show that a chromatin loop, formed by elements flanking either an enhancer or a promoter, suppresses enhancer-promoter interactions, working as an insulator. A loop formed by elements located in the region between an enhancer and a promoter, on the contrary, facilitates their interactions. We find that different mechanisms underlie insulation and facilitation; insulation occurs due to steric exclusion by the loop, and is a global effect, while facilitation occurs due to an effective shortening of the enhancer-promoter genomic distance, and is a local effect. Consistently, we find that these effects manifest quite differently for in silico 3C and microscopy. Our results show that looping interactions that do not directly involve an enhancer-promoter pair can nevertheless significantly modulate their interactions. This phenomenon is analogous to allosteric regulation in proteins, where a conformational change triggered by binding of a regulatory molecule to one site affects the state of another site., National Cancer Institute (U.S.). Physical Sciences-Oncology Center (U54-CA143874-04), Massachusetts Institute of Technology. Undergraduate Research Opportunities Program, National Science Foundation (U.S.). Massachusetts Institute of Technology. Program for Research in Mathematics, Engineering and Science for High School Students (PRIMES)