RRBS was performed as described previously51 (link),52 , starting with 100 ng of genomic DNA per sample. Custom-designed methylated and unmethylated oligonucleotides were added at a concentration of 0.1% to serve as spike-in controls for monitoring bisulfite conversion efficiency. After adaptor ligation, RRBS libraries were quantified by qPCR and pooled in combinations of six. For library enrichment, the number of PCR cycles was determined by qPCR and never exceeded 18 cycles. The library was purified twice using Agencourt AMPure XP beads (Beckman Coulter, A63880). Quality control for the final library was performed by measuring the DNA concentration with the Qubit dsDNA HS assay (ThermoFisher Scientific, Q32851) on Qubit 2.0 Fluorometer (ThermoFisher Scientific, Q32866) and by determining library fragment sizes with the Experion DNA 1K Analysis kit (Bio-Rad, 700-7107) on the Experion Automated Electrophoresis Station (Bio-Rad, 701-7000). Libraries were sequenced on Illumina HiSeq 2000/2500 machines.