DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation

Epigenome modulation potentially provides a mechanism for organisms to adapt, within and between generations. However, neither the extent to which this occurs, nor the mechanisms involved are known. Here we investigate DNA methylation variation in Swedish Arabidopsis thaliana accessions grown at two different temperatures. Environmental effects were limited to transposons, where CHH methylation was found to increase with temperature. Genome-wide association studies (GWAS) revealed that the extensive CHH methylation variation was strongly associated with genetic variants in both cis and trans, including a major trans-association close to the DNA methyltransferase CMT2. Unlike CHH methylation, CpG gene body methylation (GBM) was not affected by growth temperature, but was instead correlated with the latitude of origin. Accessions from colder regions had higher levels of GBM for a significant fraction of the genome, and this was associated with increased transcription for the genes affected. GWAS revealed that this effect was largely due to trans-acting loci, many of which showed evidence of local adaptation. DOI: http://dx.doi.org/10.7554/eLife.05255.001

[1]  D. Weigel,et al.  Evolution of DNA Methylation Patterns in the Brassicaceae is Driven by Differences in Genome Organization , 2014, PLoS genetics.

[2]  Hao Wu,et al.  R/qtl: QTL Mapping in Experimental Crosses , 2003, Bioinform..

[3]  M. Matzke,et al.  A stepwise pathway for biogenesis of 24‐nt secondary siRNAs and spreading of DNA methylation , 2009, The EMBO journal.

[4]  R. O’Hara,et al.  QST–FST comparisons: evolutionary and ecological insights from genomic heterogeneity , 2013, Nature Reviews Genetics.

[5]  J. P. Jackson,et al.  Requirement of CHROMOMETHYLASE3 for Maintenance of CpXpG Methylation , 2001, Science.

[6]  D. Patel,et al.  Non-CG methylation patterns shape the epigenetic landscape in Arabidopsis , 2013, Nature Structural & Molecular Biology.

[7]  Robert J. Schmitz,et al.  Genomic Distribution of H 3 K 9 me 2 and DNA Methylation in a Maize Genome , 2014 .

[8]  Detlef Weigel,et al.  Recombination and linkage disequilibrium in Arabidopsis thaliana , 2007, Nature Genetics.

[9]  Huy Q. Dinh,et al.  Epigenetic Regulation of Repetitive Elements Is Attenuated by Prolonged Heat Stress in Arabidopsis[W][OA] , 2010, Plant Cell.

[10]  I. Hellmann,et al.  Massive genomic variation and strong selection in Arabidopsis thaliana lines from Sweden , 2013, Nature Genetics.

[11]  Arndt von Haeseler,et al.  NextGenMap: fast and accurate read mapping in highly polymorphic genomes , 2013, Bioinform..

[12]  Bjarni J. Vilhjálmsson,et al.  An efficient multi-locus mixed model approach for genome-wide association studies in structured populations , 2012, Nature Genetics.

[13]  Bjarni J. Vilhjálmsson,et al.  Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines , 2010 .

[14]  Julie A. Law,et al.  DDR complex facilitates global association of RNA Polymerase V to promoters and evolutionarily young transposons , 2012, Nature Structural &Molecular Biology.

[15]  S. Henikoff,et al.  Genome-wide analysis of Arabidopsis thaliana DNA methylation uncovers an interdependence between methylation and transcription , 2007, Nature Genetics.

[16]  S. Jacobsen,et al.  Comprehensive Analysis of Silencing Mutants Reveals Complex Regulation of the Arabidopsis Methylome , 2013, Cell.

[17]  P. Wigge,et al.  H2A.Z-Containing Nucleosomes Mediate the Thermosensory Response in Arabidopsis , 2010, Cell.

[18]  B. Gaut,et al.  Body-methylated genes in Arabidopsis thaliana are functionally important and evolve slowly. , 2012, Molecular biology and evolution.

[19]  P. Laird,et al.  Bis-SNP: Combined DNA methylation and SNP calling for Bisulfite-seq data , 2012, Genome Biology.

[20]  A. von Haeseler,et al.  Advanced Methylome Analysis after Bisulfite Deep Sequencing: An Example in Arabidopsis , 2012, PloS one.

[21]  Dario Strbenac,et al.  Repitools: an R package for the analysis of enrichment-based epigenomic data , 2010, Bioinform..

[22]  Pjotr Prins,et al.  R/qtl: high-throughput multiple QTL mapping , 2010, Bioinform..

[23]  Century-scale methylome stability in a recently diverged Arabidopsis thaliana lineage , 2014 .

[24]  J. Hermisson,et al.  Keeping It Local: Evidence for Positive Selection in Swedish Arabidopsis thaliana , 2014, Molecular biology and evolution.

[25]  H. Vaucheret,et al.  RNA Silencing Is Resistant to Low-Temperature in Grapevine , 2013, PloS one.

[26]  Francine E. Garrett-Bakelman,et al.  methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles , 2012, Genome Biology.

[27]  David S. Lapointe,et al.  ChIPpeakAnno: a Bioconductor package to annotate ChIP-seq and ChIP-chip data , 2010, BMC Bioinformatics.

[28]  Mats E. Pettersson,et al.  Natural CMT2 Variation Is Associated With Genome-Wide Methylation Changes and Temperature Seasonality , 2013, bioRxiv.

[29]  Simon Fraser,et al.  LDheatmap: An R Function for Graphical Display of Pairwise Linkage Disequilibria Between Single Nucleotide Polymorphisms , 2006 .

[30]  Wei Li,et al.  BSMAP: whole genome bisulfite sequence MAPping program , 2009, BMC Bioinformatics.

[31]  Steven R. Eichten,et al.  Inheritance Patterns and Stability of DNA Methylation Variation in Maize Near-Isogenic Lines , 2013, Genetics.

[32]  W. Huber,et al.  Differential expression analysis for sequence count data , 2010 .

[33]  Joy Bergelson,et al.  References and Notes Supporting Online Material Adaptation to Climate across the Arabidopsis Thaliana Genome , 2022 .

[34]  N. Patterson,et al.  Estimating and interpreting FST: The impact of rare variants , 2013, Genome research.

[35]  H. Kang,et al.  Variance component model to account for sample structure in genome-wide association studies , 2010, Nature Genetics.

[36]  Christoph Lippert,et al.  LIMIX: genetic analysis of multiple traits , 2014, bioRxiv.

[37]  Nathan M. Springer,et al.  Epigenetic and Genetic Influences on DNA Methylation Variation in Maize Populations[C][W] , 2013, Plant Cell.

[38]  Vipin T. Sreedharan,et al.  RNA‐Seq Read Alignments with PALMapper , 2010, Current protocols in bioinformatics.

[39]  W. Peacock,et al.  DNA METHYLATION IN PLANTS. , 1998, Annual review of plant physiology and plant molecular biology.

[40]  E. Dermitzakis,et al.  Passive and active DNA methylation and the interplay with genetic variation in gene regulation , 2013, eLife.

[41]  Matthew D. Schultz,et al.  Patterns of Population Epigenomic Diversity , 2013, Nature.

[42]  E. Bucher,et al.  An siRNA pathway prevents transgenerational retrotransposition in plants subjected to stress , 2011, Nature.

[43]  S. Henikoff,et al.  A DNA methyltransferase homolog with a chromodomain exists in multiple polymorphic forms in Arabidopsis. , 1998, Genetics.

[44]  Robert J. Schmitz,et al.  Genomic Distribution of H3K9me2 and DNA Methylation in a Maize Genome , 2014, PloS one.

[45]  Natural CMT2 variation is associated with genome-wide methylation changes and temperature seasonality , 2014 .

[46]  R. Doerge,et al.  Epigenetic Natural Variation in Arabidopsis thaliana , 2007, PLoS biology.

[47]  Nicole Lettner,et al.  How a Retrotransposon Exploits the Plant's Heat Stress Response for Its Activation , 2014, PLoS genetics.

[48]  张静,et al.  Banana Ovate family protein MaOFP1 and MADS-box protein MuMADS1 antagonistically regulated banana fruit ripening , 2015 .

[49]  B. Langmead,et al.  BSmooth: from whole genome bisulfite sequencing reads to differentially methylated regions , 2012, Genome Biology.

[50]  D. Coleman-Derr,et al.  The Arabidopsis Nucleosome Remodeler DDM1 Allows DNA Methyltransferases to Access H1-Containing Heterochromatin , 2013, Cell.

[51]  Guo-Bo Chen,et al.  Estimating heritability of complex traits from genome-wide association studies using IBS-based Haseman–Elston regression , 2014, Front. Genet..