These are supplementary files and script files to a manuscript in which we show that DNA methylation variation plays a role in heritable phenotypic variation. DNA methylation variations were introduced into the parasite vector snail Biomphalaria glabrata with low genetic diversity by chemical treatment in F0 and followed over 3 generations using epigenetic recombinant inbred lines (epiRILs). We observed phenotypic variation in complex traits such as fecundity and susceptibility to infestation by Schistosoma mansoni and DNA methylation differences in F3. Both, increase and decrease of infestation success (up to 100% and down to 20% prevalence in epiRILs and from 86% to 94% in control RILs) indicated variation in complex resistance/compatibility trait. Average prevalence in control RILs was 84±5% but only 68±21 % in epiRILs. Fecundity also changed and was in average 47±7% in control RILs and 59±18% in epiRILs, being 12% higher in epiRILs. We found that the heritability h2 of the fecundity in the epiRILs was between 0.5 and 0.6 depending on the method used to estimate it. We developed a model for introducing epimutant offspring snails into resident susceptible populations. If genetic assimilation of the resistant phenotype occured in a small fraction of the introduced epimutant snails, we predict that the susceptible phenotype is replaced by the resistant phenotype after 50-70 generations. Produced at IHPE (http://ihpe.univ-perp.fr)

The work received funding through the Wellcome Trust strategic award 107475/Z/15/Z and through a Ph.D. grant to NL from the Region Occitanie (EPIPARA project) and the University of Perpignan Via Domitia.