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.     

Identifying multi-locus chromatin contacts in human cells using tethered multiple 3C

Background

Several recently developed experimental methods, each an extension of the chromatin conformation capture (3C) assay, have enabled the genome-wide profiling of chromatin contacts between pairs of genomic loci in 3D. Especially in complex eukaryotes, data generated by these methods, coupled with other genome-wide datasets, demonstrated that non-random chromatin folding correlates strongly with cellular processes such as gene expression and DNA replication.

Results

We describe a genome architecture assay, tethered multiple 3C (TM3C), that maps genome-wide chromatin contacts via a simple protocol of restriction enzyme digestion and religation of fragments upon agarose gel beads followed by paired-end sequencing. In addition to identifying contacts between pairs of loci, TM3C enables identification of contacts among more than two loci simultaneously. We use TM3C to assay the genome architectures of two human cell lines: KBM7, a near-haploid chronic leukemia cell line, and NHEK, a normal diploid human epidermal keratinocyte cell line. We confirm that the contact frequency maps produced by TM3C exhibit features characteristic of existing genome architecture datasets, including the expected scaling of contact probabilities with genomic distance, megabase scale chromosomal compartments and sub-megabase scale topological domains. We also confirm that TM3C captures several known cell type-specific contacts, ploidy shifts and translocations, such as Philadelphia chromosome formation (Ph+) in KBM7. We confirm a subset of the triple contacts involving the IGF2-H19 imprinting control region (ICR) using PCR analysis for KBM7 cells. Our genome-wide analysis of pairwise and triple contacts demonstrates their preference for linking open chromatin regions to each other and for linking regions with higher numbers of DNase hypersensitive sites (DHSs) to each other. For near-haploid KBM7 cells, we infer whole genome 3D models that exhibit clustering of small chromosomes with each other and large chromosomes with each other, consistent with previous studies of the genome architectures of other human cell lines.

Conclusion

TM3C is a simple protocol for ascertaining genome architecture and can be used to identify simultaneous contacts among three or four loci. Application of TM3C to a near-haploid human cell line revealed large-scale features of chromosomal organization and multi-way chromatin contacts that preferentially link regions of open chromatin.

Electronic supplementary material

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

Acetylation of Chromatin-Associated Histone H3 Lysine 56 Inhibits the Development of Encysted Artemia Embryos

Background

As a response to harsh environments, the crustacean artemia produces diapause gastrula embryos (cysts), in which cell division and embryonic development are totally arrested. This dormant state can last for very long periods but be terminated by specific environmental stimuli. Thus, artemia is an ideal model organism in which to study cell cycle arrest and embryonic development.

Principal Finding

Our study focuses on the roles of H3K56ac in the arrest of cell cycle and development during artemia diapause formation and termination. We found that the level of H3K56ac on chromatin increased during diapause formation, and decreased upon diapause termination, remaining basal level throughout subsequent embryonic development. In both HeLa cells and artemia, blocking the deacetylation with nicotinamide, a histone deacetylase inhibitor, increased the level of H3K56ac on chromatin and induced an artificial cell cycle arrest. Furthermore, we found that this arrest of the cell cycle and development was induced by H3K56ac and dephosphorylation of the checkpoint protein, retinoblastoma protein.

Conclusions/Significance

These results have revealed the dynamic change in H3K56ac on chromatin during artemia diapause formation and termination. Thus, our findings provide insight into the regulation of cell division during arrest of artemia embryonic development and provide further insight into the functions of H3K56ac.     

The Distribution of GYR- and YLP-Like Motifs in Drosophila Suggests a General Role in Cuticle Assembly and Other Protein-Protein Interactions

Background

Arthropod cuticle is composed predominantly of a self-assembling matrix of chitin and protein. Genes encoding structural cuticular proteins are remarkably abundant in arthropod genomes, yet there has been no systematic survey of conserved motifs across cuticular protein families.

Methodology/Principal Findings

Two short sequence motifs with conserved tyrosines were identified in Drosophila cuticular proteins that were similar to the GYR and YLP Interpro domains. These motifs were found in members of the CPR, Tweedle, CPF/CPFL, and (in Anopheles gambiae) CPLCG cuticular protein families, and the Dusky/Miniature family of cuticle-associated proteins. Tweedle proteins have a characteristic motif architecture that is shared with the Drosophila protein GCR1 and its orthologs in other species, suggesting that GCR1 is also cuticular. A resilin repeat, which has been shown to confer elasticity, matched one of the motifs; a number of other Drosophila proteins of unknown function exhibit a motif architecture similar to that of resilin. The motifs were also present in some proteins of the peritrophic matrix and the eggshell, suggesting molecular convergence among distinct extracellular matrices. More surprisingly, gene regulation, development, and proteolysis were statistically over-represented ontology terms for all non-cuticular matches in Drosophila. Searches against other arthropod genomes indicate that the motifs are taxonomically widespread.

Conclusions

This survey suggests a more general definition for GYR and YLP motifs and reveals their contribution to several types of extracellular matrix. They may define sites of protein interaction with DNA or other proteins, based on ontology analysis. These results can help guide experimental studies on the biochemistry of cuticle assembly.     

Optimising Controlled Human Malaria Infection Studies Using Cryopreserved P. falciparum Parasites Administered by Needle and Syringe

Background

Controlled human malaria infection (CHMI) studies have become a routine tool to evaluate efficacy of candidate anti-malarial drugs and vaccines. To date, CHMI trials have mostly been conducted using the bite of infected mosquitoes, restricting the number of trial sites that can perform CHMI studies. Aseptic, cryopreserved P. falciparum sporozoites (PfSPZ Challenge) provide a potentially more accurate, reproducible and practical alternative, allowing a known number of sporozoites to be administered simply by injection.

Methodology

We sought to assess the infectivity of PfSPZ Challenge administered in different dosing regimens to malaria-naive healthy adults (n = 18). Six participants received 2,500 sporozoites intradermally (ID), six received 2,500 sporozoites intramuscularly (IM) and six received 25,000 sporozoites IM.

Findings

Five out of six participants receiving 2,500 sporozoites ID, 3/6 participants receiving 2,500 sporozoites IM and 6/6 participants receiving 25,000 sporozoites IM were successfully infected. The median time to diagnosis was 13.2, 17.8 and 12.7 days for 2,500 sporozoites ID, 2,500 sporozoites IM and 25,000 sporozoites IM respectively (Kaplan Meier method; p = 0.024 log rank test).

Conclusions

2,500 sporozoites ID and 25,000 sporozoites IM have similar infectivities. Given the dose response in infectivity seen with IM administration, further work should evaluate increasing doses of PfSPZ Challenge IM to identify a dosing regimen that reliably infects 100% of participants.

Trial Registration

ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT01465048","term_id":"NCT01465048"}}NCT01465048     

A Two-Dimensional Simulation Model of the Bicoid Gradient in Drosophila

Background

Bicoid (Bcd) is a Drosophila morphogenetic protein responsible for patterning the anterior structures in embryos. Recent experimental studies have revealed important insights into the behavior of this morphogen gradient, making it necessary to develop a model that can recapitulate the biological features of the system, including its dynamic and scaling properties.

Methodology/Principal Findings

We present a biologically realistic 2-D model of the dynamics of the Bcd gradient in Drosophila embryos. This model is based on equilibrium binding of Bcd molecules to non-specific, low affinity DNA sites throughout the Drosophila genome. It considers both the diffusion media within which the Bcd gradient is formed and the dynamic and other relevant properties of bcd mRNA from which Bcd protein is produced. Our model recapitulates key features of the Bcd protein gradient observed experimentally, including its scaling properties and the stability of its nuclear concentrations during development. Our simulation model also allows us to evaluate the effects of other biological activities on Bcd gradient formation, including the dynamic redistribution of bcd mRNA in early embryos. Our simulation results suggest that, in our model, Bcd protein diffusion is important for the formation of an exponential gradient in embryos.

Conclusions/Significance

The 2-D model described in this report is a simple and versatile simulation procedure, providing a quantitative evaluation of the Bcd gradient system. Our results suggest an important role of Bcd binding to non-specific, low-affinity DNA sites in proper formation of the Bcd gradient in our model. They demonstrate that highly complex biological systems can be effectively modeled with relatively few parameters.     

Differential DNA methylation with age displays both common and dynamic features across human tissues that are influenced by CpG landscape

Background

DNA methylation is an epigenetic modification that changes with age in human tissues, although the mechanisms and specificity of this process are still poorly understood. We compared CpG methylation changes with age across 283 human blood, brain, kidney, and skeletal muscle samples using methylation arrays to identify tissue-specific age effects.

Results

We found age-associated CpGs (ageCGs) that are both tissue-specific and common across tissues. Tissue-specific ageCGs are frequently located outside CpG islands with decreased methylation, and common ageCGs show the opposite trend. AgeCGs are significantly associated with poorly expressed genes, but those with decreasing methylation are linked with higher tissue-specific expression levels compared with increasing methylation. Therefore, tissue-specific gene expression may protect against common age-dependent methylation. Distinguished from other tissues, skeletal muscle ageCGs are more associated with expression, enriched near genes related to myofiber contraction, and closer to muscle-specific CTCF binding sites. Kidney-specific ageCGs are more increasingly methylated compared to other tissues as measured by affiliation with kidney-specific expressed genes. Underlying chromatin features also mark common and tissue-specific age effects reflective of poised and active chromatin states, respectively. In contrast with decreasingly methylated ageCGs, increasingly methylated ageCGs are also generally further from CTCF binding sites and enriched within lamina associated domains.

Conclusions

Our data identified common and tissue-specific DNA methylation changes with age that are reflective of CpG landscape and suggests both common and unique alterations within human tissues. Our findings also indicate that a simple epigenetic drift model is insufficient to explain all age-related changes in DNA methylation.     

Loss of the Histone Pre-mRNA Processing Factor Stem-Loop Binding Protein in Drosophila Causes Genomic Instability and Impaired Cellular Proliferation

Background

Metazoan replication-dependent histone mRNAs terminate in a conserved stem-loop structure rather than a polyA tail. Formation of this unique mRNA 3′ end requires Stem-loop Binding Protein (SLBP), which directly binds histone pre-mRNA and stimulates 3′ end processing. The 3′ end stem-loop is necessary for all aspects of histone mRNA metabolism, including replication coupling, but its importance to organism fitness and genome maintenance in vivo have not been characterized.

Methodology/Principal Findings

In Drosophila, disruption of the Slbp gene prevents normal histone pre-mRNA processing and causes histone pre-mRNAs to utilize the canonical 3′ end processing pathway, resulting in polyadenylated histone mRNAs that are no longer properly regulated. Here we show that Slbp mutants display genomic instability, including loss of heterozygosity (LOH), increased presence of chromosome breaks, tetraploidy, and changes in position effect variegation (PEV). During imaginal disc growth, Slbp mutant cells show defects in S phase and proliferate more slowly than control cells.

Conclusions/Significance

These data are consistent with a model in which changing the 3′ end of histone mRNA disrupts normal replication-coupled histone mRNA biosynthesis and alters chromatin assembly, resulting in genomic instability, inhibition of cell proliferation, and impaired development.