Despite major progress in dissecting the molecular pathways that control DNA methylation patterns in plants, little is known about the mechanisms that shape plant methylomes over evolutionary time. Drawing on recent intra- and interspecific epigenomic studies, we show that methylome evolution over long timescales is largely a byproduct of genomic changes. By contrast, methylome evolution over short timescales appears to be driven mainly by spontaneous epimutational events. We argue that novel methods based on analyses of the methylation site frequency spectrum (mSFS) of natural populations can provide deeper insights into the evolutionary forces that act at each timescale.

Transposable elements (TEs) have long been known to be major contributors to plant evolution, adaptation and crop domestication. Stress-induced TE mobilization is of particular interest because it may result in novel gene regulatory pathways responding to stresses and thereby contribute to stress adaptation. Here, we investigated the genomic impacts of stress induced TE mobilization in wild type Arabidopsis plants. We find that the heat-stress responsive ONSEN TE displays an insertion site preference that is associated with specific chromatin states, especially those rich in H2A.Z histone variant and H3K27me3 histone mark. In order to better understand how novel ONSEN insertions affect the plant's response to heat stress, we carried out an in-depth transcriptomic analysis. We find that in addition to simple gene knockouts, ONSEN can produce a plethora of gene expression changes such as: constitutive activation of gene expression, alternative splicing, acquisition of heat-responsiveness, exonisation and genesis of novel non-coding and antisense RNAs. This report shows how the mobilization of a single TE-family can lead to a rapid rise of its copy number increasing the host's genome size and contribute to a broad range of transcriptomic novelty on which natural selection can then act.

Among more than 20,000 species of hermaphroditic trematodes, Schistosomatidae are unusual since they have evolved gonochorism. In schistosomes, sex is determined by a female heterogametic system, but phenotypic sexual dimorphism appears only after infection of the vertebrate definitive host. The completion of gonad maturation occurs even later, after pairing. To date, the molecular mechanisms that trigger the sexual differentiation in these species remain unknown, and in vivo studies on the developing schistosomulum stages are lacking. To study the molecular basis of sex determination and sexual differentiation in schistosomes, we investigated the whole transcriptome of the human parasite Schistosoma mansoni in a stage- and sex-comparative manner. We performed a RNA-seq on males and females for five developmental stages: cercariae larvae, three in vivo schistosomulum stages and adults. We detected 7,168 genes differentially expressed between sexes in at least one of the developmental stages, and 4,065 of them were functionally annotated. Transcriptome data were completed with H3K27me3 histone modification analysis using ChIP-Seq before (in cercariae) and after (in adults) the phenotypic sexual dimorphism appearance. In this paper we present (i) candidate determinants of the sexual differentiation, (ii) sex-biased players of the interaction with the vertebrate host, and (iii) different dynamic of the H3K27me3 histone mark between sexes as an illustration of sex-biased epigenetic landscapes. Our work presents evidence that sexual differentiation in S. mansoni is accompanied by distinct male and female transcriptional landscapes of known players of the host-parasite crosstalk, genetic determinants and epigenetic regulators. Our results suggest that such combination could lead to the optimized sexual dimorphism of this parasitic species. As S. mansoni is pathogenic for humans, this study represents a promising source of therapeutic targets, providing not only data on the parasite development in interaction with its vertebrate host, but also new insights on its reproductive function. Parasitic flatworms include more than 20,000 species that are classically hermaphrodites. Among them, the roughly hundred species of Schistosomatidae are intriguing because they are gonochoric. Schistosomes are responsible of the second most important parasitic disease worldwide, and eggs are the main cause of the inflammatory symptoms. Thus, studying the sexual reproduction mechanisms of schistosomes is of particular interest for drug development. Schistosome’s sex is genetically determined by the presence of sex chromosomes: ZZ in males or ZW in females. There is, however, no phenotypic dimorphism in the larval stages: sexual dimorphism appears only in the vertebrate host. In order to understand the molecular mechanisms underlying phenotypic sexual dimorphism, we performed a transcriptome analysis (RNA-Seq) in five different stages of the parasite lifecycle as well as a chromatin status analysis (ChIP-Seq) in the non-differentiated stage cercariae and in the adult differentiated stage, for males and females separately. Our work presents evidence that sexual differentiation in S. mansoni is accompanied by distinct male and female transcriptional landscapes of known players of the host-parasite crosstalk, developmental pathways and epigenetic regulators. Our sex-comparative approach provides therefore new potential therapeutic targets to affect development and sexual reproduction of parasite.

Chromatin structure can control gene expression and can define specific transcription states. For example, bivalent methylation of histone H3K4 and H3K27 is linked to poised transcription in vertebrate embryonic stem cells (ESC). It allows them to rapidly engage specific developmental pathways. We reasoned that non-vertebrate metazoans that encounter a similar developmental constraint (i.e. to quickly start development into a new phenotype) might use a similar system. Schistosomes are parasitic platyhelminthes that are characterized by passage through two hosts: a mollusk as intermediate host and humans or rodents as definitive host. During its development, the parasite undergoes drastic changes, most notable immediately after infection of the definitive host, i.e. during the transition from the free-swimming cercariae into adult worms. We used Chromatin Immunoprecipitation followed by massive parallel sequencing (ChIP-Seq) to analyze genome-wide chromatin structure of S. mansoni on the level of histone modifications (H3K4me3, H3K27me3, H3K9me3, and H3K9ac) in cercariae, schistosomula and adults (available at http://genome.univ-perp.fr). We saw striking differences in chromatin structure between the developmental stages, but most importantly we found that cercariae possess a specific combination of marks at the transcription start sites (TSS) that has similarities to a structure found in ESC. We demonstrate that in cercariae no transcription occurs, and we provide evidences that cercariae do not possess large numbers of canonical stem cells. We describe here a broad view on the epigenome of a metazoan parasite. Most notably, we find bivalent histone H3 methylation in cercariae. Methylation of H3K27 is removed during transformation into schistosomula (and stays absent in adults) and transcription is activated. In addition, shifts of H3K9 methylation and acetylation occur towards upstream and downstream of the transcriptional start site (TSS). We conclude that specific H3 modifications are a phylogenetically older and probably more general mechanism, i.e. not restricted to stem cells, to poise transcription. Since adult couples must form to cause the disease symptoms, changes in histone modifications appear to be crucial for pathogenesis and represent therefore a therapeutic target. The blood fluke Schistosoma mansoni causes intestinal bilharzia. The parasite has a complex life cycle in which a freshwater snail serves as intermediate host from which the human infecting larvae hatch. These larvae will actively seek skin contact, penetrate through the epithelium and start developing straight away into adult worms. Development from larvae into adults needs thorough adjustment of gene expression through repositioning or modification of proteins that are associated with DNA (the chromatin). We decided to compare the chromatin of human infective larvae (cercariae), the first developmental stage after infection of the vertebrate host (schistosomula) and adults of S. mansoni. We found that cercariae possess chromatin structures (modifications of histone H3) around the beginning of genes that are very different from schistosomula and adults. We conclude that this structure serves to keep gene transcription in a poised state, i.e. transcription is initiated and can start immediately when the blocking histone modification is removed. A similar type of histone modification was found in embryonic stem cells of vertebrates and our data indicate that it is either a more ancient and/or more general means to poise transcription than previously assumed. Since many parasites possess infective stages that develop rapidly within the host, this particular chromatin structure could be a therapeutic target for a new class of antiparasitic drugs.

Whole-genome bisulfite sequencing (WGBS) has become the standard method for interrogating plant methylomes at base resolution. However, deep WGBS measurements remain cost prohibitive for large, complex genomes and for population-level studies. As a result, most published plant methylomes are sequenced far below saturation, with a large proportion of cytosines having either missing data or insufficient coverage. Here we present METHimpute, a Hidden Markov Model (HMM) based imputation algorithm for the analysis of WGBS data. Unlike existing methods, METHimpute enables the construction of complete methylomes by inferring the methylation status and level of all cytosines in the genome regardless of coverage. Application of METHimpute to maize, rice and Arabidopsis shows that the algorithm infers cytosine-resolution methylomes with high accuracy from data as low as 6X, compared to data with 60X, thus making it a cost-effective solution for large-scale studies. METHimpute provides methylation status calls and levels for all cytosines in the genome regardless of coverage, thus yielding complete methylomes even with low-coverage WGBS datasets. The method has been extensively tested in plants, but should also be applicable to other species. An implementation is available on Bioconductor.