Epigenetic signals that direct cell type-specific interferon beta response in mouse cells.

The antiviral response induced by type I interferon (IFN) via the JAK-STAT signaling cascade activates hundreds of IFN-stimulated genes (ISGs) across human and mouse tissues but varies between cell types. However, the links between the underlying epigenetic features and the ISG profile are not well understood. We mapped ISGs, binding sites of the STAT1 and STAT2 transcription factors, chromatin accessibility, and histone H3 lysine modification by acetylation (ac) and mono-/tri-methylation (me1, me3) in mouse embryonic stem cells and fibroblasts before and after IFNß treatment. A large fraction of ISGs and STAT-binding sites was cell type specific with promoter binding of a STAT1/2 complex being a key driver of ISGs. Furthermore, STAT1/2 binding to putative enhancers induced ISGs as inferred from a chromatin co-accessibility analysis. STAT1/2 binding was dependent on the chromatin context and positively correlated with preexisting H3K4me1 and H3K27ac marks in an open chromatin state, whereas the presence of H3K27me3 had an inhibitory effect. Thus, chromatin features present before stimulation represent an additional regulatory layer for the cell type-specific antiviral response.

TelNet - a database for human and yeast genes involved in telomere maintenance.

BACKGROUND: The ends of linear chromosomes, the telomeres, comprise repetitive DNA sequences in complex with proteins that protects them from being processed by the DNA repair machinery. Cancer cells need to counteract the shortening of telomere repeats during replication for their unlimited proliferation by reactivating the reverse transcriptase telomerase or by using the alternative lengthening of telomeres (ALT) pathway. The different telomere maintenance (TM) mechanisms appear to involve hundreds of proteins but their telomere repeat length related activities are only partly understood. Currently, a database that integrates information on TM relevant genes is missing. DESCRIPTION: To provide a resource for studies that dissect TM features, we here introduce the TelNet database at http://www.cancertelsys.org/telnet/ . It offers a comprehensive compilation of more than 2000 human and 1100 yeast genes linked to telomere maintenance. These genes were annotated in terms of TM mechanism, associated specific functions and orthologous genes, a TM significance score and information from peer-reviewed literature. This TM information can be retrieved via different search and view modes and evaluated for a set of genes as demonstrated for an exemplary application. CONCLUSION: TelNet supports the annotation of genes identified from bioinformatics analysis pipelines to reveal possible connections with TM networks. We anticipate that TelNet will be a helpful resource for researchers that study telomeres.

Investigation of the spatial structure and interactions of the genome at sub-kilobase-pair resolution using T2C.

Chromosome conformation capture (3C) and its derivatives (e.g., 4C, 5C and Hi-C) are used to analyze the 3D organization of genomes. We recently developed targeted chromatin capture (T2C), an inexpensive method for studying the 3D organization of genomes, interactomes and structural changes associated with gene regulation, the cell cycle, and cell survival and development. Here, we present the protocol for T2C based on capture, describing all experimental steps and bio-informatic tools in full detail. T2C offers high resolution, a large dynamic interaction frequency range and a high signal-to-noise ratio. Its resolution is determined by the resulting fragment size of the chosen restriction enzyme, which can lead to sub-kilobase-pair resolution. T2C's high coverage allows the identification of the interactome of each individual DNA fragment, which makes binning of reads (often used in other methods) basically unnecessary. Notably, T2C requires low sequencing efforts. T2C also allows multiplexing of samples for the direct comparison of multiple samples. It can be used to study topologically associating domains (TADs), determining their position, shape, boundaries, and intra- and inter-domain interactions, as well as the composition of aggregated loops, interactions between nucleosomes, individual transcription factor binding sites, and promoters and enhancers. T2C can be performed by any investigator with basic skills in molecular biology techniques in ∼7-8 d. Data analysis requires basic expertise in bioinformatics and in Linux and Python environments.

The detailed 3D multi-loop aggregate/rosette chromatin architecture and functional dynamic organization of the human and mouse genomes.

BACKGROUND: The dynamic three-dimensional chromatin architecture of genomes and its co-evolutionary connection to its function-the storage, expression, and replication of genetic information-is still one of the central issues in biology. Here, we describe the much debated 3D architecture of the human and mouse genomes from the nucleosomal to the megabase pair level by a novel approach combining selective high-throughput high-resolution chromosomal interaction capture (T2C), polymer simulations, and scaling analysis of the 3D architecture and the DNA sequence. RESULTS: The genome is compacted into a chromatin quasi-fibre with ~5 ± 1 nucleosomes/11 nm, folded into stable ~30-100 kbp loops forming stable loop aggregates/rosettes connected by similar sized linkers. Minor but significant variations in the architecture are seen between cell types and functional states. The architecture and the DNA sequence show very similar fine-structured multi-scaling behaviour confirming their co-evolution and the above. CONCLUSIONS: This architecture, its dynamics, and accessibility, balance stability and flexibility ensuring genome integrity and variation enabling gene expression/regulation by self-organization of (in)active units already in proximity. Our results agree with the heuristics of the field and allow "architectural sequencing" at a genome mechanics level to understand the inseparable systems genomic properties.

Affinity, stoichiometry and cooperativity of heterochromatin protein 1 (HP1) binding to nucleosomal arrays.

Heterochromatin protein 1 (HP1) participates in establishing and maintaining heterochromatin via its histone-modification-dependent chromatin interactions. In recent papers HP1 binding to nucleosomal arrays was measured in vitro and interpreted in terms of nearest-neighbour cooperative binding. This mode of chromatin interaction could lead to the spreading of HP1 along the nucleosome chain. Here, we reanalysed previous data by representing the nucleosome chain as a 1D binding lattice and showed how the experimental HP1 binding isotherms can be explained by a simpler model without cooperative interactions between neighboring HP1 dimers. Based on these calculations and spatial models of dinucleosomes and nucleosome chains, we propose that binding stoichiometry depends on the nucleosome repeat length (NRL) rather than protein interactions between HP1 dimers. According to our calculations, more open nucleosome arrays with long DNA linkers are characterized by a larger number of binding sites in comparison to chains with a short NRL. Furthermore, we demonstrate by Monte Carlo simulations that the NRL dependent folding of the nucleosome chain can induce allosteric changes of HP1 binding sites. Thus, HP1 chromatin interactions can be modulated by the change of binding stoichiometry and the type of binding to condensed (methylated) and non-condensed (unmethylated) nucleosome arrays in the absence of direct interactions between HP1 dimers.

Changing chromatin fiber conformation by nucleosome repositioning.

Chromatin conformation is dynamic and heterogeneous with respect to nucleosome positions, which can be changed by chromatin remodeling complexes in the cell. These molecular machines hydrolyze ATP to translocate or evict nucleosomes, and establish loci with regularly and more irregularly spaced nucleosomes as well as nucleosome-depleted regions. The impact of nucleosome repositioning on the three-dimensional chromatin structure is only poorly understood. Here, we address this issue by using a coarse-grained computer model of arrays of 101 nucleosomes considering several chromatin fiber models with and without linker histones, respectively. We investigated the folding of the chain in dependence of the position of the central nucleosome by changing the length of the adjacent linker DNA in basepair steps. We found in our simulations that these translocations had a strong effect on the shape and properties of chromatin fibers: i), Fiber curvature and flexibility at the center were largely increased and long-range contacts between distant nucleosomes on the chain were promoted. ii), The highest destabilization of the fiber conformation occurred for a nucleosome shifted by two basepairs from regular spacing, whereas effects of linker DNA changes of ?10 bp in phase with the helical twist of DNA were minimal. iii), A fiber conformation can stabilize a regular spacing of nucleosomes inasmuch as favorable stacking interactions between nucleosomes are facilitated. This can oppose nucleosome translocations and increase the energetic costs for chromatin remodeling. Our computational modeling framework makes it possible to describe the conformational heterogeneity of chromatin in terms of nucleosome positions, and thus advances theoretical models toward a better understanding of how genome compaction and access are regulated within the cell.

Targeted Chromatin Capture (T2C): a novel high resolution high throughput method to detect genomic interactions and regulatory elements.

BACKGROUND: Significant efforts have recently been put into the investigation of the spatial organization and the chromatin-interaction networks of genomes. Chromosome conformation capture (3C) technology and its derivatives are important tools used in this effort. However, many of these have limitations, such as being limited to one viewpoint, expensive with moderate to low resolution, and/or requiring a large sequencing effort. Techniques like Hi-C provide a genome-wide analysis. However, it requires massive sequencing effort with considerable costs. Here we describe a new technique termed Targeted Chromatin Capture (T2C), to interrogate large selected regions of the genome. T2C provides an unbiased view of the spatial organization of selected loci at superior resolution (single restriction fragment resolution, from 2 to 6 kbp) at much lower costs than Hi-C due to the lower sequencing effort. RESULTS: We applied T2C on well-known model regions, the mouse ß-globin locus and the human H19/IGF2 locus. In both cases we identified all known chromatin interactions. Furthermore, we compared the human H19/IGF2 locus data obtained from different chromatin conformation capturing methods with T2C data. We observed the same compartmentalization of the locus, but at a much higher resolution (single restriction fragments vs. the common 40 kbp bins) and higher coverage. Moreover, we compared the ß-globin locus in two different biological samples (mouse primary erythroid cells and mouse fetal brain), where it is either actively transcribed or not, to identify possible transcriptional dependent interactions. We identified the known interactions in the ß-globin locus and the same topological domains in both mouse primary erythroid cells and in mouse fetal brain with the latter having fewer interactions probably due to the inactivity of the locus. Furthermore, we show that interactions due to the important chromatin proteins, Ldb1 and Ctcf, in both tissues can be analyzed easily to reveal their role on transcriptional interactions and genome folding. CONCLUSIONS: T2C is an efficient, easy, and affordable with high (restriction fragment) resolution tool to address both genome compartmentalization and chromatin-interaction networks for specific genomic regions at high resolution for both clinical and non-clinical research.

Dissecting DNA-histone interactions in the nucleosome by molecular dynamics simulations of DNA unwrapping.

The nucleosome complex of DNA wrapped around a histone protein octamer organizes the genome of eukaryotes and regulates the access of protein factors to the DNA. We performed molecular dynamics simulations of the nucleosome in explicit water to study the dynamics of its histone-DNA interactions. A high-resolution histone-DNA interaction map was derived that revealed a five-nucleotide periodicity, in which the two DNA strands of the double helix made alternating contacts. On the 100-ns timescale, the histone tails mostly maintained their initial positions relative to the DNA, and the spontaneous unwrapping of DNA was limited to 1-2 basepairs. In steered molecular dynamics simulations, external forces were applied to the linker DNA to investigate the unwrapping pathway of the nucleosomal DNA. In comparison with a nucleosome without the unstructured N-terminal histone tails, the following findings were obtained: 1), Two main barriers during unwrapping were identified at DNA position ±70 and ±45 basepairs relative to the central DNA basepair at the dyad axis. 2), DNA interactions of the histone H3 N-terminus and the histone H2A C-terminus opposed the initiation of unwrapping. 3), The N-terminal tails of H2A, H2B, and H4 counteracted the unwrapping process at later stages and were essential determinants of nucleosome dynamics. Our detailed analysis of DNA-histone interactions revealed molecular mechanisms for modulating access to nucleosomal DNA via conformational rearrangements of its structure.

Force spectroscopy of chromatin fibers: extracting energetics and structural information from Monte Carlo simulations.

The folding of the nucleosome chain into a chromatin fiber is a central factor for controlling the DNA access of protein factors involved in transcription, DNA replication and repair. Force spectroscopy experiments with chromatin fibers are ideally suited to dissect the interactions that drive this process, and to probe the underlying fiber conformation. However, the interpretation of the experimental data is fraught with difficulties due to the complex interplay of the nucleosome geometry and the different energy terms involved. Here, we apply a Monte Carlo simulation approach to derive virtual chromatin fiber force spectroscopy curves. In the simulations, the effect of the nucleosome geometry, repeat length, nucleosome-nucleosome interaction potential, and the unwrapping of the DNA from the histone protein core on the shape of the force-extension curves was investigated. These simulations provide a framework for the evaluation of experimental data sets. We demonstrate how the relative contributions of DNA bending and twisting, nucleosome unstacking and unwrapping the nucleosomal DNA from the histone octamer can be dissected for a given fiber geometry.

Visualization, analysis, and design of COMBO-FISH probes in the grid-based GLOBE 3D genome platform.

The genome architecture in cell nuclei plays an important role in modern microscopy for the monitoring of medical diagnosis and therapy since changes of function and dynamics of genes are interlinked with changing geometrical parameters. The planning of corresponding diagnostic experiments and their imaging is a complex and often interactive IT intensive challenge and thus makes high-performance grids a necessity. To detect genetic changes we recently developed a new form of fluorescence in situ hybridization (FISH) - COMBinatorial Oligonucleotide FISH (COMBO-FISH) - which labels small nucleotide sequences clustering at a desired genomic location. To achieve a unique hybridization spot other side clusters have to be excluded. Therefore, we have designed an interactive pipeline using the grid-based GLOBE 3D Genome Viewer and Platform to design and display different labelling variants of candidate probe sets. Thus, we have created a grid-based virtual "paper" tool for easy interactive calculation, analysis, management, and representation for COMBO-FISH probe design with many an advantage: Since all the calculations and analysis run in a grid, one can instantly and with great visual ease locate duplications of gene subsequences to guide the elimination of side clustering sequences during the probe design process, as well as get at least an impression of the 3D architectural embedding of the respective chromosome region, which is of major importance to estimate the hybridization probe dynamics. Beyond, even several people at different locations could work on the same process in a team wise manner. Consequently, we present how a complex interactive process can profit from grid infrastructure technology using our unique GLOBE 3D Genome Platform gateway towards a real interactive curative diagnosis planning and therapy monitoring.

Parallel high-performance grid computing: capabilities and opportunities of a novel demanding service and business class allowing highest resource efficiency.

Especially in the life-science and the health-care sectors the huge IT requirements are imminent due to the large and complex systems to be analysed and simulated. Grid infrastructures play here a rapidly increasing role for research, diagnostics, and treatment, since they provide the necessary large-scale resources efficiently. Whereas grids were first used for huge number crunching of trivially parallelizable problems, increasingly parallel high-performance computing is required. Here, we show for the prime example of molecular dynamic simulations how the presence of large grid clusters including very fast network interconnects within grid infrastructures allows now parallel high-performance grid computing efficiently and thus combines the benefits of dedicated super-computing centres and grid infrastructures. The demands for this service class are the highest since the user group has very heterogeneous requirements: i) two to many thousands of CPUs, ii) different memory architectures, iii) huge storage capabilities, and iv) fast communication via network interconnects, are all needed in different combinations and must be considered in a highly dedicated manner to reach highest performance efficiency. Beyond, advanced and dedicated i) interaction with users, ii) the management of jobs, iii) accounting, and iv) billing, not only combines classic with parallel high-performance grid usage, but more importantly is also able to increase the efficiency of IT resource providers. Consequently, the mere "yes-we-can" becomes a huge opportunity like e.g. the life-science and health-care sectors as well as grid infrastructures by reaching higher level of resource efficiency.

Exploring the conformational space of chromatin fibers and their stability by numerical dynamic phase diagrams.

The three-dimensional structure of chromatin affects DNA accessibility and is therefore a key regulator of gene expression. However, the path of the DNA between consecutive nucleosomes, and the resulting chromatin fiber organization remain controversial. The conformational space available for the folding of the nucleosome chain has been analytically described by phase diagrams with a two-angle model, which describes the chain trajectory by a DNA entry-exit angle at the nucleosome and a torsion angle between consecutive nucleosomes. Here, a novel type of numerical phase diagrams is introduced that relates the geometric phase space to the energy associated with a given chromatin conformation. The resulting phase diagrams revealed differences in the energy landscape that reflect the probability of a given conformation to form in thermal equilibrium. Furthermore, we investigated the effects of entropy and additional degrees of freedom in the dynamic phase diagrams by performing Monte Carlo simulations of the initial chain trajectories. Using our approach, we were able to demonstrate that conformations that initially were geometrically impossible could evolve into energetically favorable states in thermal equilibrium due to DNA bending and torsion. In addition, dynamic phase diagrams were applied to identify chromatin fibers that reflect certain experimentally determined features.

The GLOBE 3D Genome Platform - towards a novel system-biological paper tool to integrate the huge complexity of genome organization and function.

Genomes are tremendous co-evolutionary holistic systems for molecular storage, processing and fabrication of information. Their system-biological complexity remains, however, still largely mysterious, despite immense sequencing achievements and huge advances in the understanding of the general sequential, three-dimensional and regulatory organization. Here, we present the GLOBE 3D Genome Platform a completely novel grid based virtual "paper" tool and in fact the first system-biological genome browser integrating the holistic complexity of genomes in a single easy comprehensible platform: Based on a detailed study of biophysical and IT requirements, every architectural level from sequence to morphology of one or several genomes can be approached in a real and in a symbolic representation simultaneously and navigated by continuous scale-free zooming within a unique three-dimensional OpenGL and grid driven environment. In principle an unlimited number of multi-dimensional data sets can be visualized, customized in terms of arrangement, shape, colour, and texture etc. as well as accessed and annotated individually or in groups using internal or external data bases/facilities. Any information can be searched and correlated by importing or calculating simple relations in real-time using grid resources. A general correlation and application platform for more complex correlative analysis and a front-end for system-biological simulations both using again the huge capabilities of grid infrastructures is currently under development. Hence, the GLOBE 3D Genome Platform is an example of a grid based approach towards a virtual desktop for genomic work combining the three fundamental distributed resources: i) visual data representation, ii) data access and management, and iii) data analysis and creation. Thus, the GLOBE 3D Genome Platform is the novel system-biology oriented information system urgently needed to access, present, annotate, and to simulate the holistic genome complexity in a unique gateway towards a real understanding, educative presentation and curative manipulation planning of this tremendous evolutionary information grail - genomes.

The effect of internucleosomal interaction on folding of the chromatin fiber.

The folding of the nucleosome chain into a chromatin fiber modulates DNA accessibility and is therefore an important factor for the control of gene expression. The fiber conformation depends crucially on the interaction between individual nucleosomes. However, this parameter has not been accurately determined experimentally, and it is affected by posttranslational histone modifications and binding of chromosomal proteins. Here, the effect of different internucleosomal interaction strengths on the fiber conformation was investigated by Monte Carlo computer simulations. The fiber geometry was modeled to fit that of chicken erythrocyte chromatin, which has been examined in numerous experimental studies. In the Monte Carlo simulation, the nucleosome shape was described as an oblate spherocylinder, and a replica exchange protocol was developed to reach thermal equilibrium for a broad range of internucleosomal interaction energies. The simulations revealed the large impact of the nucleosome geometry and the nucleosome repeat length on the compaction of the chromatin fiber. At high internucleosomal interaction energies, a lateral self-association of distant fiber parts and an interdigitation of nucleosomes were apparent. These results identify key factors for the control of the compaction and higher order folding of the chromatin fiber.

Nucleosome geometry and internucleosomal interactions control the chromatin fiber conformation.

Based on model structures with atomic resolution, a coarse-grained model for the nucleosome geometry was implemented. The dependence of the chromatin fiber conformation on the spatial orientation of nucleosomes and the path and length of the linker DNA was systematically explored by Monte Carlo simulations. Two fiber types were analyzed in detail that represent nucleosome chains without and with linker histones, respectively: two-start helices with crossed-linker DNA (CL conformation) and interdigitated one-start helices (ID conformation) with different nucleosome tilt angles. The CL conformation was derived from a tetranucleosome crystal structure that was extended into a fiber. At thermal equilibrium, the fiber shape persisted but relaxed into a structure with a somewhat lower linear mass density of 3.1 +/- 0.1 nucleosomes/11 nm fiber. Stable ID fibers required local nucleosome tilt angles between 40 degrees and 60 degrees. For these configurations, much higher mass densities of up to 7.9 +/- 0.2 nucleosomes/11 nm fiber were obtained. A model is proposed, in which the transition between a CL and ID fiber is mediated by relatively small changes of the local nucleosome geometry. These were found to be in very good agreement with changes induced by linker histone H1 binding as predicted from the high resolution model structures.

Folding and organization of a contiguous chromosome region according to the gene distribution pattern in primary genomic sequence.

Specific mammalian genes functionally and dynamically associate together within the nucleus. Yet, how an array of many genes along the chromosome sequence can be spatially organized and folded together is unknown. We investigated the 3D structure of a well-annotated, highly conserved 4.3-Mb region on mouse chromosome 14 that contains four clusters of genes separated by gene "deserts." In nuclei, this region forms multiple, nonrandom "higher order" structures. These structures are based on the gene distribution pattern in primary sequence and are marked by preferential associations among multiple gene clusters. Associating gene clusters represent expressed chromatin, but their aggregation is not simply dependent on ongoing transcription. In chromosomes with aggregated gene clusters, gene deserts preferentially align with the nuclear periphery, providing evidence for chromosomal region architecture by specific associations with functional nuclear domains. Together, these data suggest dynamic, probabilistic 3D folding states for a contiguous megabase-scale chromosomal region, supporting the diverse activities of multiple genes and their conserved primary sequence organization.