Genomic impact of stress-induced transposable element mobility in Arabidopsis.

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.

Histone methylation changes are required for life cycle progression in the human parasite Schistosoma mansoni.

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.

METHimpute: imputation-guided construction of complete methylomes from WGBS data.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Frequency and mitotic heritability of epimutations in Schistosoma mansoni.

Schistosoma mansoni is a parasitic platyhelminth responsible for intestinal bilharzia. It has a complex life cycle, infecting a freshwater snail of the Biomphalaria genus, and then a mammalian host. Schistosoma mansoni adapts rapidly to new (allopatric) strains of its intermediate host. To study the importance of epimutations in this process, we infected sympatric and allopatric mollusc strains with parasite clones. ChIP-Seq was carried out on four histone modifications (H3K4me3, H3K27me3, H3K27ac and H4K20me1) in parallel with genomewide DNA resequencing (i) on parasite larvae shed by the infected snails and (ii) on adult worms that had developed from the larvae. No change in single nucleotide polymorphisms and no mobilization of transposable elements were observed, but 58-105 copy number variations (CNVs) within the parasite clones in different molluscs were detected. We also observed that the allopatric environment induces three types of chromatin structure changes: (i) host-induced changes on larvae epigenomes in 51 regions of the genome that are independent of the parasites' genetic background, (ii) spontaneous changes (not related to experimental condition or genotype of the parasite) at 64 locations and (iii) 64 chromatin structure differences dependent on the parasite genotype. Up to 45% of the spontaneous, but none of the host-induced chromatin structure changes were transmitted to adults. In our model, the environment induces epigenetic changes at specific loci but only spontaneous epimutations are mitotically heritable and have therefore the potential to contribute to transgenerational inheritance. We also show that CNVs are the only source of genetic variation and occur at the same order of magnitude as epimutations.

Combination of de novo assembly of massive sequencing reads with classical repeat prediction improves identification of repetitive sequences in Schistosoma mansoni.

The genome of the parasitic platyhelminth Schistosoma mansoni is composed of approximately 40% of repetitive sequences of which roughly 20% correspond to transposable elements. When the genome sequence became available, conventional repeat prediction programs were used to find these repeats, but only a fraction could be identified. To exhaustively characterize the repeats we applied a new massive sequencing based strategy: we re-sequenced the genome by next generation sequencing, aligned the sequencing reads to the genome and assembled all multiple-hit reads into contigs corresponding to the repetitive part of the genome. We present here, for the first time, this de novo repeat assembly strategy and we confirm that such assembly is feasible. We identified and annotated 4,143 new repeats in the S. mansoni genome. At least one third of the repeats are transcribed. This strategy allowed us also to identify 14 new microsatellite markers, which can be used for pedigree studies. Annotations and the combined (previously known and new) 5,420 repeat sequences (corresponding to 47% of the genome) are available for download (http://methdb.univ-perp.fr/downloads/).

DNA methylation dynamics during stress response in woodland strawberry (Fragaria vesca).

Environmental stresses can result in a wide range of physiological and molecular responses in plants. These responses can also impact epigenetic information in genomes, especially at the level of DNA methylation (5-methylcytosine). DNA methylation is the hallmark heritable epigenetic modification and plays a key role in silencing transposable elements (TEs). Although DNA methylation is an essential epigenetic mechanism, fundamental aspects of its contribution to stress responses and adaptation remain obscure. We investigated epigenome dynamics of wild strawberry (Fragaria vesca) in response to variable ecologically relevant environmental conditions at the DNA methylation level. F. vesca methylome responded with great plasticity to ecologically relevant abiotic and hormonal stresses. Thermal stress resulted in substantial genome-wide loss of DNA methylation. Notably, all tested stress conditions resulted in marked hot spots of differential DNA methylation near centromeric or pericentromeric regions, particularly in the non-symmetrical DNA methylation context. Additionally, we identified differentially methylated regions (DMRs) within promoter regions of transcription factor (TF) superfamilies involved in plant stress-response and assessed the effects of these changes on gene expression. These findings improve our understanding on stress-response at the epigenome level by highlighting the correlation between DNA methylation, TEs and gene expression regulation in plants subjected to a broad range of environmental stresses.

The Complete Genome of the "Flavescence Dorée" Phytoplasma Reveals Characteristics of Low Genome Plasticity.

Members of the genus 'Candidatus Phytoplasma' are obligate intracellular bacteria restricted to phloem sieve elements and are able to colonize several tissues and the hemolymph in their insect vectors. The current unfeasibility of axenic culture and the low complexity of genomic sequences are obstacles in assembling complete chromosomes. Here, a method combining pathogen DNA enrichment from infected insects and dual deep-sequencing technologies was used to obtain the complete genome of a phytoplasma causing Grapevine Flavescence dorée. The de novo assembly generated a circular chromosome of 654,223 bp containing 506 protein-coding genes. Quality assessment of the draft showed a high degree of completeness. Comparative analysis with other phytoplasmas revealed the absence of potential mobile units and a reduced amount of putative phage-derived segments, suggesting a low genome plasticity. Phylogenetic analyses identified Candidatus Phytoplasma ziziphi as the closest fully sequenced relative. The "Flavescence dorée" phytoplasma strain CH genome also encoded for several putative effector proteins potentially playing a role in pathogen virulence. The availability of this genome provides the basis for the study of the pathogenicity mechanisms and evolution of the Flavescence dorée phytoplasma.

Mutations in Fanconi anemia genes and the risk of esophageal cancer.

The incidence of esophageal squamous cell carcinoma (ESCC) is very high in northeastern Iran. Previously, we reported a strong familial component of ESCC among Turkmens, who constitute approximately one-half of the population of this region. We hypothesized that the genes which cause Fanconi anemia might be candidate genes for ESCC. We sequenced the entire coding regions of 12 Fanconi anemia genes in the germline DNA of 190 Turkmen cases of ESCC. We identified three heterozygous insertion/deletion mutations: one in FANCD2 (p.Val1233del), one in FANCE (p.Val311SerfsX2), and one in FANCL (p.Thr367AsnfsX13). All three patients had a strong family history of ESCC. In addition, four patients (out of 746 tested) were homozygous for the FANCA p.Ser858Arg mutation, compared to none of 1,373 matched controls (OR = 16.7, 95% CI = 6.2-44.2, P = 0.01). The p. Lys3326X mutation in BRCA2 (also known as Fanconi anemia gene FANCD1) was present in 27 of 746 ESCC cases and in 16 of 1,373 controls (OR = 3.38, 95% CI = 1.97-6.91, P = 0.0002). In summary, both heterozygous and homozygous mutations in several Fanconi anemia-predisposing genes are associated with an increased risk of ESCC in Iran.

The tropical coral Pocillopora acuta displays an unusual chromatin structure and shows histone H3 clipping plasticity upon bleaching.

Background: Pocillopora acuta is a hermatypic coral with strong ecological importance. Anthropogenic disturbances and global warming are major threats that can induce coral bleaching, the disruption of the mutualistic symbiosis between the coral host and its endosymbiotic algae. Previous works have shown that somaclonal colonies display different levels of survival depending on the environmental conditions they previously faced. Epigenetic mechanisms are good candidates to explain this phenomenon. However, almost no work had been published on the P. acuta epigenome, especially on histone modifications. In this study, we aim at providing the first insight into chromatin structure of this species. Methods: We aligned the amino acid sequence of P. acuta core histones with histone sequences from various phyla. We developed a centri-filtration on sucrose gradient to separate chromatin from the host and the symbiont. The presence of histone H3 protein and specific histone modifications were then detected by western blot performed on histone extraction done from bleached and healthy corals. Finally, micrococcal nuclease (MNase) digestions were undertaken to study nucleosomal organization. Results: The centri-filtration enabled coral chromatin isolation with less than 2% of contamination by endosymbiont material. Histone sequences alignments with other species show that P. acuta displays on average ~90% of sequence similarities with mice and ~96% with other corals. H3 detection by western blot showed that H3 is clipped in healthy corals while it appeared to be intact in bleached corals. MNase treatment failed to provide the usual mononucleosomal digestion, a feature shared with some cnidarian, but not all; suggesting an unusual chromatin structure. Conclusions: These results provide a first insight into the chromatin, nucleosome and histone structure of P. acuta. The unusual patterns highlighted in this study and partly shared with other cnidarian will need to be further studied to better understand its role in corals.

Divergent DNA Methylation Signatures of Juvenile Seedlings, Grafts and Adult Apple Trees.

The vast majority of previous studies on epigenetics in plants have centered on the study of inheritance of DNA methylation patterns in annual plants. In contrast, perennial plants may have the ability to accumulate changes in DNA methylation patterns over numerous years. However, currently little is known about long-lived perennial and clonally reproducing plants that may have evolved different DNA methylation inheritance mechanisms as compared to annual plants. To study the transmission of DNA methylation patterns in a perennial plant, we used apple (Malus domestica) as a model plant. First, we investigated the inheritance of DNA methylation patterns during sexual reproduction in apple by comparing DNA methylation patterns of mature trees to juvenile seedlings resulting from selfing. While we did not observe a drastic genome-wide change in DNA methylation levels, we found clear variations in DNA methylation patterns localized in regions enriched for genes involved in photosynthesis. Using transcriptomics, we also observed that genes involved in this pathway were overexpressed in seedlings. To assess how DNA methylation patterns are transmitted during clonal propagation we then compared global DNA methylation of a newly grafted tree to its mature donor tree. We identified significant, albeit weak DNA methylation changes resulting from grafting. Overall, we found that a majority of DNA methylation patterns from the mature donor tree are transmitted to newly grafted plants, however with detectable specific local differences. Both the epigenomic and transcriptomic data indicate that grafted plants are at an intermediate phase between an adult tree and seedling and inherit part of the epigenomic history of their donor tree.

Correction to: Measuring Histone Modifications in the Human Parasite Schistosoma mansoni.

Correction to: Chapter 9 in: David J. Timson (ed.), Schistosoma mansoni: Methods and Protocols, Methods in Molecular Biology, vol. 2151.

Measuring Histone Modifications in the Human Parasite Schistosoma mansoni.

DNA-binding proteins play critical roles in many major processes such as development and sexual biology of Schistosoma mansoni and are important for the pathogenesis of schistosomiasis. Chromatin immunoprecipitation (ChIP) experiments followed by sequencing (ChIP-seq) are useful to characterize the association of genomic regions with posttranslational chemical modifications of histone proteins. Challenges in the standard ChIP protocol have motivated recent enhancements in this approach, such as reducing the number of cells required and increasing the resolution. In this chapter, we describe the latest advances made by our group in the ChIP methods to improve the standard ChIP protocol to reduce the number of input cells required and to increase the resolution and robustness of ChIP in S. mansoni.

Identification and characterisation of hypomethylated DNA loci controlling quantitative resistance in Arabidopsis.

Variation in DNA methylation enables plants to inherit traits independently of changes to DNA sequence. Here, we have screened an Arabidopsis population of epigenetic recombinant inbred lines (epiRILs) for resistance against Hyaloperonospora arabidopsidis (Hpa). These lines share the same genetic background, but show variation in heritable patterns of DNA methylation. We identified four epigenetic quantitative trait loci (epiQTLs) that provide quantitative resistance without reducing plant growth or resistance to other (a)biotic stresses. Phenotypic characterisation and RNA-sequencing analysis revealed that Hpa-resistant epiRILs are primed to activate defence responses at the relatively early stages of infection. Collectively, our results show that hypomethylation at selected pericentromeric regions is sufficient to provide quantitative disease resistance, which is associated with genome-wide priming of defence-related genes. Based on comparisons of global gene expression and DNA methylation between the wild-type and resistant epiRILs, we discuss mechanisms by which the pericentromeric epiQTLs could regulate the defence-related transcriptome.

Dosage Compensation throughout the Schistosoma mansoni Lifecycle: Specific Chromatin Landscape of the Z Chromosome.

Differentiated sex chromosomes are accompanied by a difference in gene dose between X/Z-specific and autosomal genes. At the transcriptomic level, these sex-linked genes can lead to expression imbalance, or gene dosage can be compensated by epigenetic mechanisms and results into expression level equalization. Schistosoma mansoni has been previously described as a ZW species (i.e., female heterogamety, in opposition to XY male heterogametic species) with a partial dosage compensation, but underlying mechanisms are still unexplored. Here, we combine transcriptomic (RNA-Seq) and epigenetic data (ChIP-Seq against H3K4me3, H3K27me3, and H4K20me1 histone marks) in free larval cercariae and intravertebrate parasitic stages. For the first time, we describe differences in dosage compensation status in ZW females, depending on the parasitic status: free cercariae display global dosage compensation, whereas intravertebrate stages show a partial dosage compensation. We also highlight regional differences of gene expression along the Z chromosome in cercariae, but not in the intravertebrate stages. Finally, we feature a consistent permissive chromatin landscape of the Z chromosome in both sexes and stages. We argue that dosage compensation in schistosomes is characterized by chromatin remodeling mechanisms in the Z-specific region.

Atypical methylmalonic aciduria: frequency of mutations in the methylmalonyl CoA epimerase gene (MCEE).

Methylmalonic aciduria is known to result from defects in the enzyme methylmalonyl CoA mutase (MCM) (mut complementation group) and from defects in the synthesis of the MCM cofactor adenosylcobalamin (cblA, cblB, cblC, cblD, and cblF groups). Two patients who excrete methylmalonic acid have recently been shown to have a homozygous nonsense mutation in the gene coding for methylmalonyl CoA epimerase (MCEE). To further understand the cause of methylmalonic acid excretion, the MCEE gene was sequenced in 229 patients with elevations of methylmalonic acid excretion for which no cause was known. Mutations in MCEE were detected in five patients: two patients homozygous for c.139C>T, p.R47X, one patient homozygous for c.178A>C, p.K60Q, and two patients heterozygous for c.427C>T, p.R143C. Fusion of fibroblast lines from two patients homozygous for c.139C>T, p.R47X did not result in correction of [(14)C]propionate incorporation toward control values while the defect in these fibroblasts was complemented by mut, cblA, and cblB fibroblasts. Infection with wild-type MCEE cDNA resulted in correction of the biochemical phenotype in cells from both patients.

Methylome evolution in plants.

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.

Epigenetic origin of adaptive phenotypic variants in the human blood fluke Schistosoma mansoni.

BACKGROUND: Adaptive evolution is not possible without the generation of phenotypic variants. The origin of these variations has been a central topic in evolutionary biology. Up to now, it was commonly accepted that standing genetic variation is the only cause of phenotypic variants. However, epigenetic information is emerging as a complementary source of heritable phenotypic variation that contributes to evolution. The relative importance of genetics and epigenetics in generating heritable phenotypic variation is nevertheless a matter of debate. RESULTS: We used a host-parasite system to address this question. The human blood fluke Schistosoma mansoni can adapt rapidly to new intermediate snail hosts. The interaction between parasite and mollusk is characterized by a compatibility polymorphism illustrating the evolutionary dynamics in this system. The principal molecular marker for compatibility (infection success) is the expression pattern of a group of polymorphic mucins (SmPoMuc) in the parasite. We show here that chromatin structure changes as the SmPoMuc promoters are the cause for SmPoMuc transcription polymorphism leading to phenotypic novelty and increase in infection success, i.e., fitness. CONCLUSION: We establish that epigenetic changes can be the major if not only cause of adaptive phenotypic variants in Schistosoma mansoni, suggesting that epimutations can provide material for adaptive evolution in the absence of genetic variation in other systems. In addition, our results indicate that epidrugs can be used to control parasite development but also parasite evolution.

Sex-Biased Transcriptome of Schistosoma mansoni: Host-Parasite Interaction, Genetic Determinants and Epigenetic Regulators Are Associated with Sexual Differentiation.

BACKGROUND: 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. METHODOLOGY/ PRINCIPAL FINDINGS: 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. CONCLUSIONS/ SIGNIFICANCE: 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.

Macrophage scavenger receptor 1 999C>T (R293X) mutation and risk of prostate cancer.

BACKGROUND: Variants in the gene encoding the macrophage scavenger receptor 1 (MSR1(4)) protein have been identified in men with prostate cancer, and several small studies have suggested that the 999C>T (R293X) protein-truncating mutation may be associated with an increased risk for this disease. METHODS: Using large case-control, cohort, and prostate cancer family studies conducted in several Western countries, we tested for the 999C>T mutation in 2,943 men with invasive prostate carcinoma, including 401 males from multiple-case families, 1,982 cases unselected for age, and 575 men diagnosed before the age of 56 years, and in 2,870 male controls. Risk ratios were estimated by unconditional logistic regression adjusting for country and by a modified segregation analysis. A meta-analysis was conducted pooling our data with published data. RESULTS: The prevalence of MSR1*999C>T mutation carriers was 0.027 (SE, 0.003) in cases and 0.022 (SE, 0.002) in controls, and did not differ by country, ethnicity, or source. The adjusted risk ratio for prostate cancer associated with being a 999C>T carrier was 1.31 [95% confidence interval (CI), 0.93-1.84; P = 0.16]. The modified segregation analysis estimated the risk ratio to be 1.20 (95% CI, 0.87-1.66; P = 0.16). The risk ratio estimated from the meta-analysis was 1.34 (95% CI, 0.94-1.89; P = 0.10). CONCLUSION: Our large-scale analysis of case and controls from several countries found no evidence that the 999C>T mutation is associated with increased risk of prostate cancer. The meta-analysis suggests it is unlikely that this mutation confers more than a 2-fold increased risk.