Making contacts on a nucleic acid polymer

Deregulated enhancer-promoter communication in cancer through altered nuclear architecture

Enhancers are critical regulators of gene expression. Structural variations in cancer genomes can lead to enhancer hijacking, where oncogenes are activated by mistargeted enhancer activity. Novel enhancer-promoter interactions may also arise through chromosomal rearrangements that create extrachromosomal DNA elements. Additionally, fusion proteins and other mutation-induced alterations in protein properties can lead to the aberrant assembly of proteins into large complexes on the size scale of 0.1-1 μm termed onco-condensates. Transcription factors and co-activators accumulate with cis-regulatory elements in these structures, driving oncogenic programs. Here, we review current evidence of how altered genome architecture and macromolecular assembly result in deregulated enhancer-promoter communication. We discuss emerging strategies to exploit these mechanisms for clinical applications.

Heterogeneous flexibility can contribute to chromatin segregation in the cell nucleus

The highly and slightly condensed forms of chromatin, heterochromatin and euchromatin, respectively, segregate in the cell nucleus. Heterochromatin is more abundant in the nucleus periphery. Here we study the mechanism of heterochromatin segregation by modeling interphase chromosomes as diblock ring copolymers confined in a rigid spherical shell using molecular dynamics simulations. In our model, heterochromatin and euchromatin are distinguished by their bending stiffnesses only, while an interaction potential between the spherical shell and chromatin is used to model lamin-associated proteins. Our simulations indicate that in the absence of attractive interactions between the nuclear shell and the chromatin, most heterochromatin segregates towards the nuclear interior due to the depletion of less flexible heterochromatin segments from the nuclear periphery. This inverted chromatin distribution,which is opposite to the conventional case with heterochromatin dominating at the periphery, is in accord with experimental observations in rod cells. This "inversion" is also found to be independent of the heterochromatin concentration and chromosome number. The chromatin distribution at the periphery found in vivo can be recovered by further increasing the bending stiffness of heterochromatin segments or by turning on attractive interactions between the nuclear shell and heterochromatin. Our results indicate that the bending stiffness of chromatin could be a contributor to chromosome organization along with differential effects of HP1α-driven phase segregation and of loop extruders and interactions with the nuclear envelope and topological constraints.

A single dose of cocaine rewires the 3D genome structure of midbrain dopamine neurons

Midbrain dopamine neurons (DNs) respond to a first exposure to addictive drugs and play key roles in chronic drug usage1-3. As the synaptic and transcriptional changes that follow an acute cocaine exposure are mostly resolved within a few days4,5, the molecular changes that encode the long-term cellular memory of the exposure within DNs remain unknown. To investigate whether a single cocaine exposure induces long-term changes in the 3D genome structure of DNs, we applied Genome Architecture Mapping and single nucleus transcriptomic analyses in the mouse midbrain. We found extensive rewiring of 3D genome architecture at 24 hours past exposure which remains or worsens by 14 days, outlasting transcriptional responses. The cocaine-induced chromatin rewiring occurs at all genomic scales and affects genes with major roles in cocaine-induced synaptic changes. A single cocaine exposure triggers extensive long-lasting changes in chromatin condensation in post-synaptic and post-transcriptional regulatory genes, for example the unfolding of Rbfox1 which becomes most prominent 14 days post exposure. Finally, structurally remodeled genes are most expressed in a specific DN sub-type characterized by low expression of the dopamine auto-receptor Drd2, a key feature of highly cocaine-sensitive cells. These results reveal an important role for long-lasting 3D genome remodelling in the cellular memory of a single cocaine exposure, providing new hypotheses for understanding the inception of drug addiction and 3D genome plasticity.

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.

Spatial organization of transcribed eukaryotic genes

Despite the well-established role of nuclear organization in the regulation of gene expression, little is known about the reverse: how transcription shapes the spatial organization of the genome. Owing to the small sizes of most previously studied genes and the limited resolution of microscopy, the structure and spatial arrangement of a single transcribed gene are still poorly understood. Here we study several long highly expressed genes and demonstrate that they form open-ended transcription loops with polymerases moving along the loops and carrying nascent RNAs. Transcription loops can span across micrometres, resembling lampbrush loops and polytene puffs. The extension and shape of transcription loops suggest their intrinsic stiffness, which we attribute to decoration with multiple voluminous nascent ribonucleoproteins. Our data contradict the model of transcription factories and suggest that although microscopically resolvable transcription loops are specific for long highly expressed genes, the mechanisms underlying their formation could represent a general aspect of eukaryotic transcription.

Genome reconstruction and haplotype phasing using chromosome conformation capture methodologies

Genomic analysis of individuals or organisms is predicated on the availability of high-quality reference and genotype information. With the rapidly dropping costs of high-throughput DNA sequencing, this is becoming readily available for diverse organisms and for increasingly large populations of individuals. Despite these advances, there are still aspects of genome sequencing that remain challenging for existing sequencing methods. This includes the generation of long-range contiguity during genome assembly, identification of structural variants in both germline and somatic tissues, the phasing of haplotypes in diploid organisms and the resolution of genome sequence for organisms derived from complex samples. These types of information are valuable for understanding the role of genome sequence and genetic variation on genome function, and numerous approaches have been developed to address them. Recently, chromosome conformation capture (3C) experiments, such as the Hi-C assay, have emerged as powerful tools to aid in these challenges for genome reconstruction. We will review the current use of Hi-C as a tool for aiding in genome sequencing, addressing the applications, strengths, limitations and potential future directions for the use of 3C data in genome analysis. We argue that unique features of Hi-C experiments make this data type a powerful tool to address challenges in genome sequencing, and that future integration of Hi-C data with alternative sequencing assays will facilitate the continuing revolution in genomic analysis and genome sequencing.

Using computational modelling to reveal mechanisms of epigenetic Polycomb control

The Polycomb system is essential for stable gene silencing in many organisms. This regulation is achieved in part through addition of the histone modifications H3K27me2/me3 by Polycomb Repressive Complex 2 (PRC2). These modifications are believed to be the causative epigenetic memory elements of PRC2-mediated silencing. As these marks are stored locally in the chromatin, PRC2-based memory is a cis-acting system. A key feature of stable epigenetic memory in cis is PRC2-mediated, self-reinforcing feedback from K27-methylated histones onto nearby histones in a read-write paradigm. However, it was not clear under what conditions such feedback can lead to stable memory, able, for example, to survive the perturbation of histone dilution at DNA replication. In this context, computational modelling has allowed a rigorous exploration of possible underlying memory mechanisms and has also greatly accelerated our understanding of switching between active and silenced states. Specifically, modelling has predicted that switching and memory at Polycomb loci is digital, with a locus being either active or inactive, rather than possessing intermediate, smoothly varying levels of activation. Here, we review recent advances in models of Polycomb control, focusing on models of epigenetic switching through nucleation and spreading of H3K27me2/me3. We also examine models that incorporate transcriptional feedback antagonism and those including bivalent chromatin states. With more quantitative experimental data on histone modification kinetics, as well as single-cell resolution data on transcription and protein levels for PRC2 targets, we anticipate an expanded need for modelling to help dissect increasingly interconnected and complex memory mechanisms.

Serpentine: a flexible 2D binning method for differential Hi-C analysis

Motivation

Hi-C contact maps reflect the relative contact frequencies between pairs of genomic loci, quantified through deep sequencing. Differential analyses of these maps enable downstream biological interpretations. However, the multi-fractal nature of the chromatin polymer inside the cellular envelope results in contact frequency values spanning several orders of magnitude: contacts between loci pairs separated by large genomic distances are much sparser than closer pairs. The same is true for poorly covered regions, such as repeated sequences. Both distant and poorly covered regions translate into low signal-to-noise ratios. There is no clear consensus to address this limitation.

Results

We present Serpentine, a fast, flexible procedure operating on raw data, which considers the contacts in each region of a contact map. Binning is performed only when necessary on noisy regions, preserving informative ones. This results in high-quality, low-noise contact maps that can be conveniently visualized for rigorous comparative analyses.

Availability and implementation

Serpentine is available on the PyPI repository and https://github.com/koszullab/serpentine; documentation and tutorials are provided at https://serpentine.readthedocs.io/en/latest/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Characterization of Gene Repression by Designed Transcription Activator-like Effector Dimer Proteins

Gene regulation by control of transcription initiation is a fundamental property of living cells. Much of our understanding of gene repression originated from studies of the Escherichia coli lac operon switch, in which DNA looping plays an essential role. To validate and generalize principles from lac for practical applications, we previously described artificial DNA looping driven by designed transcription activator-like effector dimer (TALED) proteins. Because TALE monomers bind the idealized symmetrical lac operator sequence in two orientations, our prior studies detected repression due to multiple DNA loops. We now quantitatively characterize gene repression in living E. coli by a collection of individual TALED loops with systematic loop length variation. Fitting of a thermodynamic model allows unequivocal demonstration of looping and comparison of the engineered TALED repression system with the natural lac repressor system.

instaGRAAL: chromosome-level quality scaffolding of genomes using a proximity ligation-based scaffolder

Hi-C exploits contact frequencies between pairs of loci to bridge and order contigs during genome assembly, resulting in chromosome-level assemblies. Because few robust programs are available for this type of data, we developed instaGRAAL, a complete overhaul of the GRAAL program, which has adapted the latter to allow efficient assembly of large genomes. instaGRAAL features a number of improvements over GRAAL, including a modular correction approach that optionally integrates independent data. We validate the program using data for two brown algae, and human, to generate near-complete assemblies with minimal human intervention.

Sensitive effect of linker histone binding mode and subtype on chromatin condensation

The complex role of linker histone (LH) on chromatin compaction regulation has been highlighted by recent discoveries of the effect of LH binding variability and isoforms on genome structure and function. Here we examine the effect of two LH variants and variable binding modes on the structure of chromatin fibers. Our mesoscale modeling considers oligonucleosomes with H1C and H1E, bound in three different on and off-dyad modes, and spanning different LH densities (0.5–1.6 per nucleosome), over a wide range of physiologically relevant nucleosome repeat lengths (NRLs). Our studies reveal an LH-variant and binding-mode dependent heterogeneous ensemble of fiber structures with variable packing ratios, sedimentation coefficients, and persistence lengths. For maximal compaction, besides dominantly interacting with parental DNA, LHs must have strong interactions with nonparental DNA and promote tail/nonparental core interactions. An off-dyad binding of H1E enables both; others compromise compaction for bendability. We also find that an increase of LH density beyond 1 is best accommodated in chromatosomes with one on-dyad and one off-dyad LH. We suggest that variable LH binding modes and concentrations are advantageous, allowing tunable levels of chromatin condensation and DNA accessibility/interactions. Thus, LHs add another level of epigenetic regulation of chromatin.

peakC: a flexible, non-parametric peak calling package for 4C and Capture-C data

It is becoming increasingly clear that chromosome organization plays an important role in gene regulation. High-resolution methods such as 4C, Capture-C and promoter capture Hi-C (PCHiC) enable the study of chromatin loops such as those formed between promoters and enhancers or CTCF/cohesin binding sites. An important aspect of 4C/Capture-C/PCHiC analyses is the reliable identification of chromatin loops, preferably not based on visual inspection of a DNA contact profile, but on reproducible statistical analysis that robustly scores interaction peaks in the non-uniform contact background. Here, we present peakC, an R package for the analysis of 4C/Capture-C/PCHiC data. We generated 4C data for 13 viewpoints in two tissues in at least triplicate to test our methods. We developed a non-parametric peak caller based on rank-products. Sampling analysis shows that not read depth but template quality is the most important determinant of success in 4C experiments. By performing peak calling on single experiments we show that the peak calling results are similar to the replicate experiments, but that false positive rates are significantly reduced by performing replicates. Our software is user-friendly and enables robust peak calling for one-vs-all chromosome capture experiments. peakC is available at: https://github.com/deWitLab/peakC.

Transcription-driven genome organization: a model for chromosome structure and the regulation of gene expression tested through simulations

Current models for the folding of the human genome see a hierarchy stretching down from chromosome territories, through A/B compartments and topologically-associating domains (TADs), to contact domains stabilized by cohesin and CTCF. However, molecular mechanisms underlying this folding, and the way folding affects transcriptional activity, remain obscure. Here we review physical principles driving proteins bound to long polymers into clusters surrounded by loops, and present a parsimonious yet comprehensive model for the way the organization determines function. We argue that clusters of active RNA polymerases and their transcription factors are major architectural features; then, contact domains, TADs and compartments just reflect one or more loops and clusters. We suggest tethering a gene close to a cluster containing appropriate factors—a transcription factory—increases the firing frequency, and offer solutions to many current puzzles concerning the actions of enhancers, super-enhancers, boundaries and eQTLs (expression quantitative trait loci). As a result, the activity of any gene is directly influenced by the activity of other transcription units around it in 3D space, and this is supported by Brownian-dynamics simulations of transcription factors binding to cognate sites on long polymers.

Protein-mediated looping of DNA under tension requires supercoiling

Protein-mediated DNA looping is ubiquitous in chromatin organization and gene regulation, but to what extent supercoiling or nucleoid associated proteins promote looping is poorly understood. Using the lac repressor (LacI), a paradigmatic loop-mediating protein, we measured LacI-induced looping as a function of either supercoiling or the concentration of the HU protein, an abundant nucleoid protein in Escherichia coli. Negative supercoiling to physiological levels with magnetic tweezers easily drove the looping probability from 0 to 100% in single DNA molecules under slight tension that likely exists in vivo. In contrast, even saturating (micromolar) concentrations of HU could not raise the looping probability above 30% in similarly stretched DNA or 80% in DNA without tension. Negative supercoiling is required to induce significant looping of DNA under any appreciable tension.

Two Old Dogs, One New Trick: A Review of RNA Polymerase and Ribosome Interactions during Transcription-Translation Coupling

The coupling of transcription and translation is more than mere translation of an mRNA that is still being transcribed. The discovery of physical interactions between RNA polymerase and ribosomes has spurred renewed interest into this long-standing paradigm of bacterial molecular biology. Here, we provide a concise presentation of recent insights gained from super-resolution microscopy, biochemical, and structural work, including cryo-EM studies. Based on the presented data, we put forward a dynamic model for the interaction between RNA polymerase and ribosomes, in which the interactions are repeatedly formed and broken. Furthermore, we propose that long intervening nascent RNA will loop out and away during the forming the interactions between the RNA polymerase and ribosomes. By comparing the effect of the direct interactions between RNA polymerase and ribosomes with those that transcription factors NusG and RfaH mediate, we submit that two distinct modes of coupling exist: Factor-free and factor-mediated coupling. Finally, we provide a possible framework for transcription-translation coupling and elude to some open questions in the field.

Bacterial gene control by DNA looping using engineered dimeric transcription activator like effector (TALE) proteins

Genetic switches must alternate between states whose probabilities are dependent on regulatory signals. Classical examples of transcriptional control in bacteria depend on repressive DNA loops anchored by proteins whose structures are sensitive to small molecule inducers or co-repressors. We are interested in exploiting these natural principles to engineer artificial switches for transcriptional control of bacterial genes. Here, we implement designed homodimeric DNA looping proteins ('Transcription Activator-Like Effector Dimers'; TALEDs) for this purpose in living bacteria. Using well-studied FKBP dimerization domains, we build switches that mimic regulatory characteristics of classical Escherichia coli lactose, galactose and tryptophan operon promoters, including induction or co-repression by small molecules. Engineered DNA looping using TALEDs is thus a new approach to tuning gene expression in bacteria. Similar principles may also be applicable for gene control in eukaryotes.

Programmable DNA looping using engineered bivalent dCas9 complexes

DNA looping is a ubiquitous and critical feature of gene regulation. Although DNA looping can be efficiently detected, tools to readily manipulate DNA looping are limited. Here we develop CRISPR-based DNA looping reagents for creation of programmable DNA loops. Cleavage-defective Cas9 proteins of different specificity are linked by heterodimerization or translational fusion to create bivalent complexes able to link two separate DNA regions. After model-directed optimization, the reagents are validated using a quantitative DNA looping assay in E. coli. Looping efficiency is ~15% for a 4.7 kb loop, but is significantly improved by loop multiplexing with additional guides. Bivalent dCas9 complexes are also used to activate endogenous norVW genes by rewiring chromosomal DNA to bring distal enhancer elements to the gene promoters. Such reagents should allow manipulation of DNA looping in a variety of cell types, aiding understanding of endogenous loops and enabling creation of new regulatory connections.

DNA loops are a ubiquitious feature of gene regulation across the kingdoms of life. Here the authors design a Cas9-based dimerization system for inducing DNA loops in E. coli, allowing activation and rewiring of gene expression.

Efficient chromosomal-scale DNA looping in Escherichia coli using multiple DNA-looping elements

Genes are frequently regulated by interactions between proteins that bind to the DNA near the gene and proteins that bind to DNA sites located far away, with the intervening DNA looped out. But it is not understood how efficient looping can occur when the sites are very far apart. We develop a simple theoretical framework that relates looping efficiency to the energetic cost and benefit of looping, allowing prediction of the efficiency of single or multiple nested loops at different distances. Measurements of absolute loop efficiencies for Lac repressor and λ CI using gene expression reporters in Escherichia coli cells show that, as predicted by the model, long-range DNA looping between a pair of sites can be strongly enhanced by the use of nested DNA loops or by the use of additional protein-binding sequences. A combination of these approaches was able to generate efficient DNA looping at a 200 kb distance.

Random positioning of nucleosomes enhances heritable bistability

Chromosomal regions are often dynamically modified by histones, leading to the uncertainty of nucleosome positions. Experiments have provided evidence for this randomness, but it is unclear how it impacts epigenetic heritability. Here, by analyzing a mechanic model at the molecular level, which considers three representative types of nucleosomes (unmodified, methylated, and acetylated) and dynamic nucleosome modifications, we find that in contrast to the equidistance partition of nucleosomes, random partition can significantly enhance heritable bistability. Moreover, the more "chaotic" the nucleosome positions are, the better the heritable bistability is, in contrast to the previous view. In both cases of nucleosome positioning, heritable bistability occurs only when the total nucleosome number is beyond a threshold, and it depends strongly on the allocation rate that enzymes regulate transitions between different nucleosome types. Thus, we conclude that random positioning of nucleosomes is an unneglectable factor impacting heritable bistability. A point worth mentioning is that our model established on a master equation can easily be extended to include other more complex processes underlying dynamic nucleosome modifications.

Proximity ligation scaffolding and comparison of two Trichoderma reesei strains genomes

BACKGROUND

The presence of low complexity and repeated regions in genomes often results in difficulties to assemble sequencing data into full chromosomes. However, the availability of full genome scaffolds is essential to several investigations, regarding for instance the evolution of entire clades, the analysis of chromosome rearrangements, and is pivotal to sexual crossing studies. In non-conventional but industrially relevant model organisms, such as the ascomycete Trichoderma reesei, a complete genome assembly is seldom available.

RESULTS

The chromosome scaffolds of T. reesei QM6a and Rut-C30 strains have been generated using a contact genomic/proximity ligation genomic approach. The original reference assembly, encompassing dozens of scaffolds, was reorganized into two sets of seven chromosomes. Chromosomal contact data also allowed to characterize 10-40 kb, gene-free, AT-rich (76%) regions corresponding to the T. reesei centromeres. Large chromosomal rearrangements (LCR) in Rut-C30 were then characterized, in agreement with former studies, and the position of LCR breakpoints used to assess the likely chromosome structure of other T. reesei strains [QM9414, CBS999.97 (1-1, re), and QM9978]. In agreement with published results, we predict that the numerous chromosome rearrangements found in highly mutated industrial strains may limit the efficiency of sexual reproduction for their improvement.

CONCLUSIONS

The GRAAL program allowed us to generate the karyotype of the Rut-C30 strain, and from there to predict chromosome structure for most T. reesei strains for which sequence is available. This method that exploits proximity ligation sequencing approach is a fast, cheap, and straightforward way to characterize both chromosome structure and centromere sequences and is likely to represent a popular convenient alternative to expensive and work-intensive resequencing projects.

The second decade of 3C technologies: detailed insights into nuclear organization

The relevance of three-dimensional (3D) genome organization for transcriptional regulation and thereby for cellular fate at large is now widely accepted. Our understanding of the fascinating architecture underlying this function is based on microscopy studies as well as the chromosome conformation capture (3C) methods, which entered the stage at the beginning of the millennium. The first decade of 3C methods rendered unprecedented insights into genome topology. Here, we provide an update of developments and discoveries made over the more recent years. As we discuss, established and newly developed experimental and computational methods enabled identification of novel, functionally important chromosome structures. Regulatory and architectural chromatin loops throughout the genome are being cataloged and compared between cell types, revealing tissue invariant and developmentally dynamic loops. Architectural proteins shaping the genome were disclosed, and their mode of action is being uncovered. We explain how more detailed insights into the 3D genome increase our understanding of transcriptional regulation in development and misregulation in disease. Finally, to help researchers in choosing the approach best tailored for their specific research question, we explain the differences and commonalities between the various 3C-derived methods.

Scaffolding bacterial genomes and probing host-virus interactions in gut microbiome by proximity ligation (chromosome capture) assay

The biochemical activities of microbial communities, or microbiomes, are essential parts of environmental and animal ecosystems. The dynamics, balance, and effects of these communities are strongly influenced by phages present in the population. Being able to characterize bacterium-phage relationships is therefore essential to investigate these ecosystems to the full extent of their complexity. However, this task is currently limited by (i) the ability to characterize complete bacterial and viral genomes from a complex mix of species and (ii) the difficulty to assign phage sequences to their bacterial hosts. We show that both limitations can be circumvented using meta3C, an experimental and computational approach that exploits the physical contacts between DNA molecules to infer their proximity. In a single experiment, dozens of bacterial and phage genomes present in a complex mouse gut microbiota were assembled and scaffolded de novo. The phage genomes were then assigned to their putative bacterial hosts according to the physical contacts between the different DNA molecules, opening new perspectives for a comprehensive picture of the genomic structure of the gut flora. Therefore, this work holds far-reaching implications for human health studies aiming to bridge the virome to the microbiome.

Cap-dependent translation is mediated by 'RNA looping' rather than 'ribosome scanning'

The 40S ribosomal subunit cannot directly recognize the start codon of eukaryotic mRNAs. Instead, it recognizes the start codon after its association with the 5′-cap structure via translation initiation factors. Base-by-base inspection of the 5′UTR by a scanning ribosome is the generally accepted hypothesis of start codon selection. As part of an effort to confirm the underlying mechanism of start codon selection by the 40S ribosome, we investigated the role of eIF4G, which participates in the recruitment of 40S ribosomes to various translation enhancers, such as 5′-cap structure, poly(A) tail, and several internal ribosome entry sites. We found that an artificial translation factor composed of recombinant eIF4G fused with MS2 greatly enhanced translation of an upstream reporter gene when it was tethered to the 3′UTR. These data suggest that the 40S ribosome recruited to a translation enhancer can find the start codon by looping of the intervening RNA segment. The ‘RNA-looping’ hypothesis of translation start codon recognition was further supported by an analysis of the effect of 5′UTR length on translation efficiency and the mathematically predicted probability of RNA-loop–mediated interactions between the start codon and the 40S ribosome associated at the 5′-end.

Spatial genome organization and cognition

Nonrandom chromosomal conformations, including promoter-enhancer loopings that bypass kilobases or megabases of linear genome, provide a crucial layer of transcriptional regulation and move vast amounts of non-coding sequence into the physical proximity of genes that are important for neurodevelopment, cognition and behaviour. Activity-regulated changes in the neuronal '3D genome' could govern transcriptional mechanisms associated with learning and plasticity, and loop-bound intergenic and intronic non-coding sequences have been implicated in psychiatric and adult-onset neurodegenerative disease. Recent studies have begun to clarify the roles of spatial genome organization in normal and abnormal cognition.

The Hierarchy of Transcriptional Activation: From Enhancer to Promoter

Regulatory elements (enhancers) that are remote from promoters play a critical role in the spatial, temporal, and physiological control of gene expression. Studies on specific loci, together with genome-wide approaches, suggest that there may be many common mechanisms involved in enhancer-promoter communication. Here, we discuss the multiprotein complexes that are recruited to enhancers and the hierarchy of events taking place between regulatory elements and promoters.

Entropy gives rise to topologically associating domains

We investigate chromosome organization within the nucleus using polymer models whose formulation is closely guided by experiments in live yeast cells. We employ bead-spring chromosome models together with loop formation within the chains and the presence of nuclear bodies to quantify the extent to which these mechanisms shape the topological landscape in the interphase nucleus. By investigating the genome as a dynamical system, we show that domains of high chromosomal interactions can arise solely from the polymeric nature of the chromosome arms due to entropic interactions and nuclear confinement. In this view, the role of bio-chemical related processes is to modulate and extend the duration of the interacting domains.

Effect of Interaction between Chromatin Loops on Cell-to-Cell Variability in Gene Expression

According to recent experimental evidence, the interaction between chromatin loops, which can be characterized by three factors-connection pattern, distance between regulatory elements, and communication form, play an important role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effect that addresses the question of how changes in these factors affect variability at the expression level in a systematic rather than case-by-case fashion. Here we make such an effort, based on a mechanic model that maps three fundamental patterns for two interacting DNA loops into a 4-state model of stochastic transcription. We first show that in contrast to side-by-side loops, nested loops enhance mRNA expression and reduce expression noise whereas alternating loops have just opposite effects. Then, we compare effects of facilitated tracking and direct looping on gene expression. We find that the former performs better than the latter in controlling mean expression and in tuning expression noise, but this control or tuning is distance-dependent, remarkable for moderate loop lengths, and there is a limit loop length such that the difference in effect between two communication forms almost disappears. Our analysis and results justify the facilitated chromatin-looping hypothesis.

A p53-bound enhancer region controls a long intergenic noncoding RNA required for p53 stress response

Genome-wide chromatin studies identified the tumor suppressor p53 as both a promoter and an enhancer-binding transcription factor. As an enhancer factor, p53 can induce local production of enhancer RNAs, as well as transcriptional activation of distal neighboring genes. Beyond the regulation of protein-coding genes, p53 has the capacity to regulate long intergenic noncoding RNA molecules (lincRNAs); however, their importance to the p53 tumor suppressive function remains poorly characterized. Here, we identified and characterized a novel p53-bound intronic enhancer that controls the expression of its host, the lincRNA00475 (linc-475). We demonstrate the requirement of linc-475 for the proper induction of a p53-dependent cell cycle inhibitory response. We further confirm the functional importance of linc-475 in the maintenance of CDKN1A/p21 levels, a cell cycle inhibitor and a major p53 target gene, following p53 activation. Interestingly, loss of linc-475 reduced the binding of both p53 and RNA polymerase II (RNAPII) to the promoter of p21, attenuating its transcription rate following p53 activation. Altogether, our data suggest a direct role of p53-bound enhancer domains in the activation of lincRNAs required for an efficient p53 transcriptional response.

Invariant TAD Boundaries Constrain Cell-Type-Specific Looping Interactions between Promoters and Distal Elements around the CFTR Locus

Three-dimensional genome structure plays an important role in gene regulation. Globally, chromosomes are organized into active and inactive compartments while, at the gene level, looping interactions connect promoters to regulatory elements. Topologically associating domains (TADs), typically several hundred kilobases in size, form an intermediate level of organization. Major questions include how TADs are formed and how they are related to looping interactions between genes and regulatory elements. Here we performed a focused 5C analysis of a 2.8 Mb chromosome 7 region surrounding CFTR in a panel of cell types. We find that the same TAD boundaries are present in all cell types, indicating that TADs represent a universal chromosome architecture. Furthermore, we find that these TAD boundaries are present irrespective of the expression and looping of genes located between them. In contrast, looping interactions between promoters and regulatory elements are cell-type specific and occur mostly within TADs. This is exemplified by the CFTR promoter that in different cell types interacts with distinct sets of distal cell-type-specific regulatory elements that are all located within the same TAD. Finally, we find that long-range associations between loci located in different TADs are also detected, but these display much lower interaction frequencies than looping interactions within TADs. Interestingly, interactions between TADs are also highly cell-type-specific and often involve loci clustered around TAD boundaries. These data point to key roles of invariant TAD boundaries in constraining as well as mediating cell-type-specific long-range interactions and gene regulation.

Identifying and Reducing Systematic Errors in Chromosome Conformation Capture Data

Chromosome conformation capture (3C)-based techniques have recently been used to uncover the mystic genomic architecture in the nucleus. These techniques yield indirect data on the distances between genomic loci in the form of contact frequencies that must be normalized to remove various errors. This normalization process determines the quality of data analysis. In this study, we describe two systematic errors that result from the heterogeneous local density of restriction sites and different local chromatin states, methods to identify and remove those artifacts, and three previously described sources of systematic errors in 3C-based data: fragment length, mappability, and local DNA composition. To explain the effect of systematic errors on the results, we used three different published data sets to show the dependence of the results on restriction enzymes and experimental methods. Comparison of the results from different restriction enzymes shows a higher correlation after removing systematic errors. In contrast, using different methods with the same restriction enzymes shows a lower correlation after removing systematic errors. Notably, the improved correlation of the latter case caused by systematic errors indicates that a higher correlation between results does not ensure the validity of the normalization methods. Finally, we suggest a method to analyze random error and provide guidance for the maximum reproducibility of contact frequency maps.

SCRaMbLE generates designed combinatorial stochastic diversity in synthetic chromosomes

Synthetic chromosome rearrangement and modification by loxP-mediated evolution (SCRaMbLE) generates combinatorial genomic diversity through rearrangements at designed recombinase sites. We applied SCRaMbLE to yeast synthetic chromosome arm synIXR (43 recombinase sites) and then used a computational pipeline to infer or unscramble the sequence of recombinations that created the observed genomes. Deep sequencing of 64 synIXR SCRaMbLE strains revealed 156 deletions, 89 inversions, 94 duplications, and 55 additional complex rearrangements; several duplications are consistent with a double rolling circle mechanism. Every SCRaMbLE strain was unique, validating the capability of SCRaMbLE to explore a diverse space of genomes. Rearrangements occurred exclusively at designed loxPsym sites, with no significant evidence for ectopic rearrangements or mutations involving synthetic regions, the 99% nonsynthetic nuclear genome, or the mitochondrial genome. Deletion frequencies identified genes required for viability or fast growth. Replacement of 3′ UTR by non-UTR sequence had surprisingly little effect on fitness. SCRaMbLE generates genome diversity in designated regions, reveals fitness constraints, and should scale to simultaneous evolution of multiple synthetic chromosomes.

Distinct polymer physics principles govern chromatin dynamics in mouse and Drosophila topological domains

Background

In higher eukaryotes, the genome is partitioned into large "Topologically Associating Domains" (TADs) in which the chromatin displays favoured long-range contacts. While a crumpled/fractal globule organization has received experimental supports at higher-order levels, the organization principles that govern chromatin dynamics within these TADs remain unclear. Using simple polymer models, we previously showed that, in mouse liver cells, gene-rich domains tend to adopt a statistical helix shape when no significant locus-specific interaction takes place.

Results

Here, we use data from diverse 3C-derived methods to explore chromatin dynamics within mouse and Drosophila TADs. In mouse Embryonic Stem Cells (mESC), that possess large TADs (median size of 840 kb), we show that the statistical helix model, but not globule models, is relevant not only in gene-rich TADs, but also in gene-poor and gene-desert TADs. Interestingly, this statistical helix organization is considerably relaxed in mESC compared to liver cells, indicating that the impact of the constraints responsible for this organization is weaker in pluripotent cells. Finally, depletion of histone H1 in mESC alters local chromatin flexibility but not the statistical helix organization. In Drosophila, which possesses TADs of smaller sizes (median size of 70 kb), we show that, while chromatin compaction and flexibility are finely tuned according to the epigenetic landscape, chromatin dynamics within TADs is generally compatible with an unconstrained polymer configuration.

Conclusions

Models issued from polymer physics can accurately describe the organization principles governing chromatin dynamics in both mouse and Drosophila TADs. However, constraints applied on this dynamics within mammalian TADs have a peculiar impact resulting in a statistical helix organization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1786-8) contains supplementary material, which is available to authorized users.

Why the activity of a gene depends on its neighbors

Sixty years ago, the position of a gene on a chromosome was seen to be a major determinant of gene activity; however, position effects are rarely central to current discussions of gene expression. We describe a comprehensive and simplifying view of how position in 1D sequence and 3D nuclear space underlies expression. We suggest that apparently-different regulatory motifs including enhancers, silencers, insulators, barriers, and boundaries act similarly - they are active promoters that tether target genes close to, or distant from, appropriate transcription sites or 'factories'. We also suggest that any active transcription unit regulates the firing of its neighbors - and thus can be categorized as one or other type of motif; this is consistent with expression quantitative trait loci (eQTLs) being widely dispersed.

Inferring 3D chromatin structure using a multiscale approach based on quaternions

Background

The knowledge of the spatial organisation of the chromatin fibre in cell nuclei helps researchers to understand the nuclear machinery that regulates dna activity. Recent experimental techniques of the type Chromosome Conformation Capture (3c, or similar) provide high-resolution, high-throughput data consisting in the number of times any possible pair of dna fragments is found to be in contact, in a certain population of cells. As these data carry information on the structure of the chromatin fibre, several attempts have been made to use them to obtain high-resolution 3d reconstructions of entire chromosomes, or even an entire genome. The techniques proposed treat the data in different ways, possibly exploiting physical-geometric chromatin models. One popular strategy is to transform contact data into Euclidean distances between pairs of fragments, and then solve a classical distance-to-geometry problem.

Results

We developed and tested a reconstruction technique that does not require translating contacts into distances, thus avoiding a number of related drawbacks. Also, we introduce a geometrical chromatin chain model that allows us to include sound biochemical and biological constraints in the problem. This model can be scaled at different genomic resolutions, where the structures of the coarser models are influenced by the reconstructions at finer resolutions. The search in the solution space is then performed by a classical simulated annealing, where the model is evolved efficiently through quaternion operators. The presence of appropriate constraints permits the less reliable data to be overlooked, so the result is a set of plausible chromatin configurations compatible with both the data and the prior knowledge.

Conclusions

To test our method, we obtained a number of 3d chromatin configurations from hi-c data available in the literature for the long arm of human chromosome 1, and validated their features against known properties of gene density and transcriptional activity. Our results are compatible with biological features not introduced a priori in the problem: structurally different regions in our reconstructions highly correlate with functionally different regions as known from literature and genomic repositories.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0667-0) contains supplementary material, which is available to authorized users.

Real sequence effects on the search dynamics of transcription factors on DNA

Recent experiments show that transcription factors (TFs) indeed use the facilitated diffusion mechanism to locate their target sequences on DNA in living bacteria cells: TFs alternate between sliding motion along DNA and relocation events through the cytoplasm. From simulations and theoretical analysis we study the TF-sliding motion for a large section of the DNA-sequence of a common E. coli strain, based on the two-state TF-model with a fast-sliding search state and a recognition state enabling target detection. For the probability to detect the target before dissociating from DNA the TF-search times self-consistently depend heavily on whether or not an auxiliary operator (an accessible sequence similar to the main operator) is present in the genome section. Importantly, within our model the extent to which the interconversion rates between search and recognition states depend on the underlying nucleotide sequence is varied. A moderate dependence maximises the capability to distinguish between the main operator and similar sequences. Moreover, these auxiliary operators serve as starting points for DNA looping with the main operator, yielding a spectrum of target detection times spanning several orders of magnitude. Auxiliary operators are shown to act as funnels facilitating target detection by TFs.

Contact genomics: scaffolding and phasing (meta)genomes using chromosome 3D physical signatures

High-throughput DNA sequencing technologies are fuelling an accelerating trend to assemble de novo or resequence the genomes of numerous species as well as to complete unfinished assemblies. While current DNA sequencing technologies remain limited to reading stretches of a few hundreds or thousands of base pairs, experimental and computational methods are continuously improving with the goal of assembling entire genomes from large numbers of short DNA sequences. However, the algorithms that piece together DNA strands face important limitations due, notably, to the presence of repeated sequences or of multiple haplotypes within one genome, thus leaving many assemblies incomplete. Recently, the realization that the physical contacts experienced by a portion of a DNA molecule could be used as a robust and quantitative assay to determine its genomic position has led to the emerging field of contact genomics, which promises to revolutionize current genome assembly approaches by exploiting the flexible polymer properties of chromosomes. Here we review the current applications of contact genomics to genome scaffolding, haplotyping and metagenomic assembly, then outline the future developments we envision.

Translation initiation mediated by RNA looping

Eukaryotic translation initiation commences at the initiation codon near the 5' end of mRNA by a 40S ribosomal subunit, and the recruitment of a 40S ribosome to an mRNA is facilitated by translation initiation factors interacting with the m(7)G cap and/or poly(A) tail. The 40S ribosome recruited to an mRNA is then transferred to the AUG initiation codon with the help of translation initiation factors. To understand the mechanism by which the ribosome finds an initiation codon, we investigated the role of eIF4G in finding the translational initiation codon. An artificial polypeptide eIF4G fused with MS2 was localized downstream of the reporter gene through MS2-binding sites inserted in the 3' UTR of the mRNA. Translation of the reporter was greatly enhanced by the eIF4G-MS2 fusion protein regardless of the presence of a cap structure. Moreover, eIF4G-MS2 tethered at the 3' UTR enhanced translation of the second cistron of a dicistronic mRNA. The encephalomyocarditis virus internal ribosome entry site, a natural translational-enhancing element facilitating translation through an interaction with eIF4G, positioned downstream of a reporter gene, also enhanced translation of the upstream gene in a cap-independent manner. Finally, we mathematically modeled the effect of distance between the cap structure and initiation codon on the translation efficiency of mRNAs. The most plausible explanation for translational enhancement by the translational-enhancing sites is recognition of the initiation codon by the ribosome bound to the ribosome-recruiting sites through "RNA looping." The RNA looping hypothesis provides a logical explanation for augmentation of translation by enhancing elements located upstream and/or downstream of a protein-coding region.

Metagenomic chromosome conformation capture (meta3C) unveils the diversity of chromosome organization in microorganisms

Genomic analyses of microbial populations in their natural environment remain limited by the difficulty to assemble full genomes of individual species. Consequently, the chromosome organization of microorganisms has been investigated in a few model species, but the extent to which the features described can be generalized to other taxa remains unknown. Using controlled mixes of bacterial and yeast species, we developed meta3C, a metagenomic chromosome conformation capture approach that allows characterizing individual genomes and their average organization within a mix of organisms. Not only can meta3C be applied to species already sequenced, but a single meta3C library can be used for assembling, scaffolding and characterizing the tridimensional organization of unknown genomes. By applying meta3C to a semi-complex environmental sample, we confirmed its promising potential. Overall, this first meta3C study highlights the remarkable diversity of microorganisms chromosome organization, while providing an elegant and integrated approach to metagenomic analysis.

DOI:http://dx.doi.org/10.7554/eLife.03318.001

Microbial communities play vital roles in the environment and sustain animal and plant life. Marine microbes are part of the ocean's food chain; soil microbes support the turnover of major nutrients and facilitate plant growth; and the microbial communities residing in the human gut support digestion and the immune system, among other roles. These communities are very complex systems, often containing 1000s of different species engaged in co-dependent relationships, and are therefore very difficult to study.

The entire DNA sequence of an organism constitutes its genome, and much of this genetic information is stored in large structures called chromosomes. Examining the genome of a species can provide important clues about its lifestyle and how it evolved. To do this, DNA is extracted from cells and is then usually cut into smaller fragments, amplified, and sequenced. The small stretches of sequence obtained, called reads, are finally assembled, yielding ideally the complete genome of the organism under study.

Metagenomics attempts to interpret the combined genome of all the different species in a microbial community and has been instrumental in deciphering how the different species interact with each other. Metagenomics involves sequencing stretches of the community's DNA and matching these pieces to individual species to ultimately assemble whole genomes. While this may be a relatively straightforward task for communities that contain only a handful of members, the metagenomes derived from complex microbial communities are huge, fragmented, and incomplete. This often makes it very difficult or even nearly impossible to match the inferred DNA stretches to individual species.

A method called chromosome conformation capture (or ‘3C’ for short) can reveal the physical contacts between different regions of a chromosome and between the different chromosomes of a cell. How often each of these chromosomal contacts occurs provides a kind of physical signature to each genome and each individual chromosome within it.

Marbouty et al. took advantage of these interactions to develop a technique that combines metagenomics and chromosome conformation capture—called meta3C—that can analyze the DNA of many different species mixed together. Testing meta3C on artificial mixtures of a few species of yeast or bacteria showed that meta3C can separate the genomes of the different species without any prior knowledge of the composition of the mix. In a single experiment, meta3C can identify individual chromosomes, match each of them to its species of origin, and reveal the three-dimensional structure of each genome in the mix. Further tests showed that meta3C can also interpret more complex communities where the number and types of the species present are not known.

Meta3C holds great promise for understanding how microbial communities work and how the genomes of the species within a community are organized. However, further developments of the technique will be required to investigate communities as diverse as those present in most natural environments.

DOI:http://dx.doi.org/10.7554/eLife.03318.002

High-quality genome (re)assembly using chromosomal contact data

Closing gaps in draft genome assemblies can be costly and time-consuming, and published genomes are therefore often left ‘unfinished.’ Here we show that genome-wide chromosome conformation capture (3C) data can be used to overcome these limitations, and present a computational approach rooted in polymer physics that determines the most likely genome structure using chromosomal contact data. This algorithm—named GRAAL—generates high-quality assemblies of genomes in which repeated and duplicated regions are accurately represented and offers a direct probabilistic interpretation of the computed structures. We first validated GRAAL on the reference genome of Saccharomyces cerevisiae, as well as other yeast isolates, where GRAAL recovered both known and unknown complex chromosomal structural variations. We then applied GRAAL to the finishing of the assembly of Trichoderma reesei and obtained a number of contigs congruent with the know karyotype of this species. Finally, we showed that GRAAL can accurately reconstruct human chromosomes from either fragments generated in silico or contigs obtained from de novo assembly. In all these applications, GRAAL compared favourably to recently published programmes implementing related approaches.

The correct assembly of genomes from sequencing data remains a challenge due to difficulties in correctly assigning the location of repeated DNA elements. Here the authors describe GRAAL, an algorithm that utilizes genome-wide chromosome contact data within a probabilistic framework to produce accurate genome assemblies.

DNA supercoiling: a regulatory signal for the λ repressor

Topoisomerases, polymerases, and the chirality introduced by the binding of histones or nucleoid-associated proteins affect DNA supercoiling in vivo. However, supercoiling is not just a by-product of DNA metabolism. Supercoiling is an indicator of cell health, it modifies the accessibility of chromatin, and coordinates the transcription of genes. This suggests that regulatory, protein-mediated loops in DNA may sense supercoiling of the genome in which they are embedded. The λ repressor (CI) maintains the quiescent (lysogenic) transcriptome of bacteriophage λ in infected Escherichia coli. CI-mediated looping prevents overexpression of the repressor protein to preserve sensitivity to conditions that trigger virulence (lysis). Experiments were performed to assess how well the CI-mediated DNA loop traps superhelicity and determine whether supercoiling enhances CI-mediated DNA looping. CI oligomers partitioned plasmids into topological domains and prevented the passage of supercoiling between them. Furthermore, in single DNA molecules stretched and twisted with magnetic tweezers, levels of superhelical density confined in CI-mediated DNA loops ranged from -15% or +11%. Finally, in DNA under tensions that may occur in vivo, supercoiling lowered the free energy of loop formation and was essential for DNA looping. Supercoiling-enhanced looping can influence the maintenance of lysogeny in the λ repressor system; it can encode sensitivity to the energy level of the cell and creates independent topological domains of distinct superhelical density.

Quantitation of interactions between two DNA loops demonstrates loop domain insulation in E. coli cells

Eukaryotic gene regulation involves complex patterns of long-range DNA-looping interactions between enhancers and promoters, but how these specific interactions are achieved is poorly understood. Models that posit other DNA loops--that aid or inhibit enhancer-promoter contact--are difficult to test or quantitate rigorously in eukaryotic cells. Here, we use the well-characterized DNA-looping proteins Lac repressor and phage λ CI to measure interactions between pairs of long DNA loops in E. coli cells in the three possible topological arrangements. We find that side-by-side loops do not affect each other. Nested loops assist each other's formation consistent with their distance-shortening effect. In contrast, alternating loops, where one looping element is placed within the other DNA loop, inhibit each other's formation, thus providing clear support for the loop domain model for insulation. Modeling shows that combining loop assistance and loop interference can provide strong specificity in long-range interactions.

Specificity, propagation, and memory of pericentric heterochromatin

The cell establishes heritable patterns of active and silenced chromatin via interacting factors that set, remove, and read epigenetic marks. To understand how the underlying networks operate, we have dissected transcriptional silencing in pericentric heterochromatin (PCH) of mouse fibroblasts. We assembled a quantitative map for the abundance and interactions of 16 factors related to PCH in living cells and found that stably bound complexes of the histone methyltransferase SUV39H1/2 demarcate the PCH state. From the experimental data, we developed a predictive mathematical model that explains how chromatin-bound SUV39H1/2 complexes act as nucleation sites and propagate a spatially confined PCH domain with elevated histone H3 lysine 9 trimethylation levels via chromatin dynamics. This “nucleation and looping” mechanism is particularly robust toward transient perturbations and stably maintains the PCH state. These features make it an attractive model for establishing functional epigenetic domains throughout the genome based on the localized immobilization of chromatin-modifying enzymes.

Targeted sequencing by proximity ligation for comprehensive variant detection and local haplotyping

Despite developments in targeted gene sequencing and whole-genome analysis techniques, the robust detection of all genetic variation, including structural variants, in and around genes of interest and in an allele-specific manner remains a challenge. Here we present targeted locus amplification (TLA), a strategy to selectively amplify and sequence entire genes on the basis of the crosslinking of physically proximal sequences. We show that, unlike other targeted re-sequencing methods, TLA works without detailed prior locus information, as one or a few primer pairs are sufficient for sequencing tens to hundreds of kilobases of surrounding DNA. This enables robust detection of single nucleotide variants, structural variants and gene fusions in clinically relevant genes, including BRCA1 and BRCA2, and enables haplotyping. We show that TLA can also be used to uncover insertion sites and sequences of integrated transgenes and viruses. TLA therefore promises to be a useful method in genetic research and diagnostics when comprehensive or allele-specific genetic information is needed.

Mathematical model of alternative mechanism of telomere length maintenance

Biopolymer length regulation is a complex process that involves a large number of biological, chemical, and physical subprocesses acting simultaneously across multiple spatial and temporal scales. An illustrative example important for genomic stability is the length regulation of telomeres-nucleoprotein structures at the ends of linear chromosomes consisting of tandemly repeated DNA sequences and a specialized set of proteins. Maintenance of telomeres is often facilitated by the enzyme telomerase but, particularly in telomerase-free systems, the maintenance of chromosomal termini depends on alternative lengthening of telomeres (ALT) mechanisms mediated by recombination. Various linear and circular DNA structures were identified to participate in ALT, however, dynamics of the whole process is still poorly understood. We propose a chemical kinetics model of ALT with kinetic rates systematically derived from the biophysics of DNA diffusion and looping. The reaction system is reduced to a coagulation-fragmentation system by quasi-steady-state approximation. The detailed treatment of kinetic rates yields explicit formulas for expected size distributions of telomeres that demonstrate the key role played by the J factor, a quantitative measure of bending of polymers. The results are in agreement with experimental data and point out interesting phenomena: an appearance of very long telomeric circles if the total telomere density exceeds a critical value (excess mass) and a nonlinear response of the telomere size distributions to the amount of telomeric DNA in the system. The results can be of general importance for understanding dynamics of telomeres in telomerase-independent systems as this mode of telomere maintenance is similar to the situation in tumor cells lacking telomerase activity. Furthermore, due to its universality, the model may also serve as a prototype of an interaction between linear and circular DNA structures in various settings.

fourSig: a method for determining chromosomal interactions in 4C-Seq data

The ability to correlate chromosome conformation and gene expression gives a great deal of information regarding the strategies used by a cell to properly regulate gene activity. 4C-Seq is a relatively new and increasingly popular technology where the set of genomic interactions generated by a single point in the genome can be determined. 4C-Seq experiments generate large, complicated data sets and it is imperative that signal is properly distinguished from noise. Currently, there are a limited number of methods for analyzing 4C-Seq data. Here, we present a new method, fourSig, which in addition to being precise and simple to use also includes a new feature that prioritizes detected interactions. Our results demonstrate the efficacy of fourSig with previously published and novel 4C-Seq data sets and show that our significance prioritization correlates with the ability to reproducibly detect interactions among replicates.

Quantitation of the DNA tethering effect in long-range DNA looping in vivo and in vitro using the Lac and λ repressors

Efficient and specific interactions between proteins bound to the same DNA molecule can be dependent on the length of the DNA tether that connects them. Measurement of the strength of this DNA tethering effect has been largely confined to short separations between sites, and it is not clear how it contributes to long-range DNA looping interactions, such as occur over separations of tens to hundreds of kilobase pairs in vivo. Here, gene regulation experiments using the LacI and λ CI repressors, combined with mathematical modeling, were used to quantitate DNA tethering inside Escherichia coli cells over the 250- to 10,000-bp range. Although LacI and CI loop DNA in distinct ways, measurements of the tethering effect were very similar for both proteins. Tethering strength decreased with increasing separation, but even at 5- to 10-kb distances, was able to increase contact probability 10- to 20-fold and drive efficient looping. Tethering in vitro with the Lac repressor was measured for the same 600-to 3,200-bp DNAs using tethered particle motion, a single molecule technique, and was 5- to 45-fold weaker than in vivo over this range. Thus, the enhancement of looping seen previously in vivo at separations below 500 bp extends to large separations, underlining the need to understand how in vivo factors aid DNA looping. Our analysis also suggests how efficient and specific looping could be achieved over very long DNA separations, such as what occurs between enhancers and promoters in eukaryotic cells.

Physical origin of the contact frequency in chromosome conformation capture data

Physical proximity between each pair of genomic loci in a nucleus is measured as a form of contact frequency in chromosome conformation capture-based methods. Complexity of chromosome structure in interphase can be characterized by measuring a statistical property of physical distance between genomic loci according to genomic separation along single chromatids. To find a relationship between the physical distance and the contact frequency, we propose a polymer model derived from the Langevin equation. The model is derived by considering a structure of a chromosome as a trajectory of a particle, where each consecutive segment in the chromosome corresponds to a transient position in the trajectory over time. Using chromosome conformation capture data, we demonstrate the functional relationship between the two quantities. The physical distances derived from the mean contact frequencies by the model show a good correlation with those from experimental data. From the model, we present that the mean contact frequency curve can be divided into three components that arise from different physical origins and show that the contact frequency is proportional to the contact surface area, not to the volume of segments suggested by the fractal globule model. The model explains both a decaying pattern of the contact frequency and the biphasic relationship between the physical distance and the genomic length.

Establishing epigenetic domains via chromatin-bound histone modifiers

The eukaryotic nucleus harbors the DNA genome, which associates with histones and other chromosomal proteins into a complex referred to as chromatin. It provides an additional layer of so-called epigenetic information via histone modifications and DNA methylation on top of the DNA sequence that determines the cell's active gene expression program. The nucleus is devoid of internal organelles separated by membranes. Thus, free diffusive transport of proteins and RNA can occur throughout the space accessible for a given macromolecule. At the same time, chromatin is partitioned into different specialized structures such as nucleoli, chromosome territories, and heterochromatin domains that serve distinct functions. Here, we address the question of how the activity of chromatin-modifying enzymes is confined to chromatin subcompartments. We discuss mechanisms for establishing activity gradients of diffusive chromatin-modifying enzymes that could give rise to distinct chromatin domains within the cell nucleus. Interestingly, such gradients might directly result from immobilization of the enzymes on the flexible chromatin chain. Thus, locus-specific tethering of these enzymes to chromatin could have the potential to establish, maintain, or modulate epigenetic patterns of characteristic domain size.