Stochastic changes in cytosine methylation are a source of heritable epigenetic and phenotypic diversity in plants. Using the model plant Arabidopsis thaliana, we derive robust estimates of the rate at which methylation is spontaneously gained (forward epimutation) or lost (backward epimutation) at individual cytosines and construct a comprehensive picture of the epimutation landscape in this species. We demonstrate that the dynamic interplay between forward and backward epimutations is modulated by genomic context and show that subtle contextual differences have profoundly shaped patterns of methylation diversity in A. thaliana natural populations over evolutionary timescales. Theoretical arguments indicate that the epimutation rates reported here are high enough to rapidly uncouple genetic from epigenetic variation, but low enough for new epialleles to sustain long-term selection responses. Our results provide new insights into methylome evolution and its population-level consequences.epigenetics | epimutation | DNA methylation | evolution | Arabidopsis P lant genomes make extensive use of cytosine methylation to control the expression of transposable elements (TEs) and genes (1). Despite its tight regulation, methylation losses or gains at individual cytosines or clusters of cytosines can emerge spontaneously, in an event termed "epimutation" (2, 3). Many examples of segregating epimutations have been documented in experimental and wild populations of plants and in some cases contribute to heritable variation in phenotypes independently of DNA sequence variation (4, 5). These observations have led to much speculation about the role of DNA methylation in plant evolution (6-8), and its potential in breeding programs (9). In the model plant Arabidopsis thaliana, spontaneous methylation changes at CG dinucleotides accumulate in a rapid but nonlinear fashion over generations (2,3,10), thus pointing to high forward-backward epimutation rates (11). Precise estimates of these rates are necessary to be able to quantify the long-term dynamics of epigenetic variation under laboratory or natural conditions, and to understand the molecular mechanisms that drive methylome evolution (12-14). Here we combine theoretical modeling with high-resolution methylome analysis of multiple independent A. thaliana mutation accumulation (MA) lines (15), including measurements of methylation changes in continuous generations, to obtain robust estimates of forward and backward epimutation rates. ResultsWe joined whole-genome MethylC-seq (16) data from two earlier MA studies (2, 3) with extensive multigenerational MethylC-seq measurements from three additional MA lines (Fig. 1A and SI Appendix, Tables S1-S6). The first of these new MA lines (MA1 3) was propagated for 30 generations and includes measurements for 13 (nearly) consecutive generations (Fig. 1A). The other two MA lines (MA2 3) were propagated for 17 generations and were measured every four generations on average (Fig. 1A). These new data therefore allowed us to track epimutation...

BackgroundChromosome instability leads to aneuploidy, a state in which cells have abnormal numbers of chromosomes, and is found in two out of three cancers. In a chromosomal instable p53 deficient mouse model with accelerated lymphomagenesis, we previously observed whole chromosome copy number changes affecting all lymphoma cells. This suggests that chromosome instability is somehow suppressed in the aneuploid lymphomas or that selection for frequently lost/gained chromosomes out-competes the CIN-imposed mis-segregation.ResultsTo distinguish between these explanations and to examine karyotype dynamics in chromosome instable lymphoma, we use a newly developed single-cell whole genome sequencing (scWGS) platform that provides a complete and unbiased overview of copy number variations (CNV) in individual cells. To analyse these scWGS data, we develop AneuFinder, which allows annotation of copy number changes in a fully automated fashion and quantification of CNV heterogeneity between cells. Single-cell sequencing and AneuFinder analysis reveals high levels of copy number heterogeneity in chromosome instability-driven murine T-cell lymphoma samples, indicating ongoing chromosome instability. Application of this technology to human B cell leukaemias reveals different levels of karyotype heterogeneity in these cancers.ConclusionOur data show that even though aneuploid tumours select for particular and recurring chromosome combinations, single-cell analysis using AneuFinder reveals copy number heterogeneity. This suggests ongoing chromosome instability that other platforms fail to detect. As chromosome instability might drive tumour evolution, karyotype analysis using single-cell sequencing technology could become an essential tool for cancer treatment stratification.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-0971-7) contains supplementary material, which is available to authorized users.

No abstract

BackgroundAlzheimer’s disease (AD) is a neurodegenerative disease of the brain and the most common form of dementia in the elderly. Aneuploidy, a state in which cells have an abnormal number of chromosomes, has been proposed to play a role in neurodegeneration in AD patients. Several studies using fluorescence in situ hybridization have shown that the brains of AD patients contain an increased number of aneuploid cells. However, because the reported rate of aneuploidy in neurons ranges widely, a more sensitive method is needed to establish a possible role of aneuploidy in AD pathology.ResultsIn the current study, we used a novel single-cell whole genome sequencing (scWGS) approach to assess aneuploidy in isolated neurons from the frontal cortex of normal control individuals (n = 6) and patients with AD (n = 10). The sensitivity and specificity of our method was shown by the presence of three copies of chromosome 21 in all analyzed neuronal nuclei of a Down’s syndrome sample (n = 36). Very low levels of aneuploidy were found in the brains from control individuals (n = 589) and AD patients (n = 893). In contrast to other studies, we observe no selective gain of chromosomes 17 or 21 in neurons of AD patients.ConclusionscWGS showed no evidence for common aneuploidy in normal and AD neurons. Therefore, our results do not support an important role for aneuploidy in neuronal cells in the pathogenesis of AD. This will need to be confirmed by future studies in larger cohorts.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-0976-2) contains supplementary material, which is available to authorized users.

The field of epigenomics has rapidly progressed from the study of individual reference epigenomes to surveying epigenomic variation in populations. Recent studies in a number of species, from yeast to humans, have begun to dissect the cis- and trans-regulatory genetic mechanisms that shape patterns of population epigenomic variation at the level of single epigenetic marks, as well as at the level of integrated chromatin state maps. We show that this information is paving the way towards a more complete understanding of the heritable basis underlying population epigenomic variation. We also highlight important conceptual challenges when interpreting results from these genetic studies, particularly in plants, in which epigenomic variation can be determined both by genetic and epigenetic inheritance.

Histones are highly covalently modified, but the functions of many of these modifications remain unknown. In particular, it is unclear how histone marks are coupled to cellular metabolism and how this coupling affects chromatin architecture. We identified histone H3 Lys14 (H3K14) as a site of propionylation and butyrylation in vivo and carried out the first systematic characterization of histone propionylation. We found that H3K14pr and H3K14bu are deposited by histone acetyltransferases, are preferentially enriched at promoters of active genes and are recognized by acylation-state-specific reader proteins. In agreement with these findings, propionyl-CoA was able to stimulate transcription in an in vitro transcription system. Notably, genome-wide H3 acylation profiles were redefined following changes to the metabolic state, and deletion of the metabolic enzyme propionyl-CoA carboxylase altered global histone propionylation levels. We propose that histone propionylation, acetylation and butyrylation may act in combination to promote high transcriptional output and to couple cellular metabolism with chromatin structure and function.

Epigenetic mechanisms and chromatin structure play an important role in development. Their impact is therefore expected to be strong in parasites with complex life cycles and multiple, strikingly different, developmental stages, i.e. developmental plasticity. Some studies have already described how the chromatin structure, through histone modifications, varies from a developmental stage to another in a few unicellular parasites. While H3K4me3 profiles remain relatively constant, H3K27 trimethylation and bivalent methylation show strong variation. Inhibitors (A366 and GSK343) of H3K27 histone methyltransferase activity in S. mansoni efficiently blocked miracidium to sporocyst transition indicating that H3K27 trimethylation is required for life cycle progression. As S. mansoni is a multicellular parasite that significantly affects both the health and economy of endemic areas, a better understanding of fluke developmental processes within the definitive host will likely highlight novel disease control strategies. Towards this goal, we also studied H4K20me1 in female cercariae and adults. In particular, we found that bivalent trimethylation of H3K4 and H3K27 at the transcription start site of genes is a landmark of the cercarial stage. In cercariae, H3K27me3 presence and strong enrichment in H4K20me1 over long regions (10–100 kb) is associated with development related genes. Here, we provide a broad overview of the chromatin structure of a metazoan parasite throughout its most important lifecycle stages. The five developmental stages studied here present distinct chromatin structures, indicating that histone methylation plays an important role during development. Hence, components of the histone methylation (and demethylation) machinery may provide suitable Schistosomiasis control targets.

Motivation Strand-seq is a specialized single-cell DNA sequencing technique centered around the directionality of single-stranded DNA. Computational tools for Strand-seq analyses must capture the strand-specific information embedded in these data. Results Here we introduce breakpointR, an R/Bioconductor package specifically tailored to process and interpret single-cell strand-specific sequencing data obtained from Strand-seq. We developed breakpointR to detect local changes in strand directionality of aligned Strand-seq data, to enable fine-mapping of sister chromatid exchanges, germline inversion and to support global haplotype assembly. Given the broad spectrum of Strand-seq applications we expect breakpointR to be an important addition to currently available tools and extend the accessibility of this novel sequencing technique. Availability and implementation R/Bioconductor package https://bioconductor.org/packages/breakpointR. Supplementary information Supplementary data are available at Bioinformatics online.