SRX3401307 GSM2859860 SRP125148 SRS2695452 GSM2859860 Bisulfite-Seq GENOMIC Reduced Representation Genomic DNA was isolated by phenol-chloroform extraction followed by ethanol precipitation. Two aliquots of DNA were separately digested with either TaqI or MspI restriction enzymes over night and subsequently pooled for all downstream processing. Digested DNA was end-repaired, A-tailed, and adapters containing methylated cytosines ligated on using the KAPA HyperPrep kit according to the manufacturers instructions. Adapter-ligated DNA was bisulfite converted over using the EpiTech Bisulfite kit (Qiagen, Inc). Bisulfite treated libraries were PCR amplified using indexed primers and the HiFi HotStart Uracil+ polymerase (KAPA Biosystems). Illumina HiSeq 2500 gds 302859860 GSM2859860 GEO Accession GSM2859860 SRX3401308 GSM2859861 SRP125148 SRS2695453 GSM2859861 Bisulfite-Seq GENOMIC Reduced Representation Genomic DNA was isolated by phenol-chloroform extraction followed by ethanol precipitation. Two aliquots of DNA were separately digested with either TaqI or MspI restriction enzymes over night and subsequently pooled for all downstream processing. Digested DNA was end-repaired, A-tailed, and adapters containing methylated cytosines ligated on using the KAPA HyperPrep kit according to the manufacturers instructions. Adapter-ligated DNA was bisulfite converted over using the EpiTech Bisulfite kit (Qiagen, Inc). Bisulfite treated libraries were PCR amplified using indexed primers and the HiFi HotStart Uracil+ polymerase (KAPA Biosystems). Illumina HiSeq 2500 gds 302859861 GSM2859861 GEO Accession GSM2859861 SRX3401309 GSM2859862 SRP125148 SRS2695454 GSM2859862 Bisulfite-Seq GENOMIC Reduced Representation Genomic DNA was isolated by phenol-chloroform extraction followed by ethanol precipitation. Two aliquots of DNA were separately digested with either TaqI or MspI restriction enzymes over night and subsequently pooled for all downstream processing. Digested DNA was end-repaired, A-tailed, and adapters containing methylated cytosines ligated on using the KAPA HyperPrep kit according to the manufacturers instructions. Adapter-ligated DNA was bisulfite converted over using the EpiTech Bisulfite kit (Qiagen, Inc). Bisulfite treated libraries were PCR amplified using indexed primers and the HiFi HotStart Uracil+ polymerase (KAPA Biosystems). Illumina HiSeq 2500 gds 302859862 GSM2859862 GEO Accession GSM2859862 SRX3401310 GSM2859863 SRP125148 SRS2695455 GSM2859863 Bisulfite-Seq GENOMIC Reduced Representation Genomic DNA was isolated by phenol-chloroform extraction followed by ethanol precipitation. Two aliquots of DNA were separately digested with either TaqI or MspI restriction enzymes over night and subsequently pooled for all downstream processing. Digested DNA was end-repaired, A-tailed, and adapters containing methylated cytosines ligated on using the KAPA HyperPrep kit according to the manufacturers instructions. Adapter-ligated DNA was bisulfite converted over using the EpiTech Bisulfite kit (Qiagen, Inc). Bisulfite treated libraries were PCR amplified using indexed primers and the HiFi HotStart Uracil+ polymerase (KAPA Biosystems). Illumina HiSeq 2500 gds 302859863 GSM2859863 GEO Accession GSM2859863 SRX3401311 GSM2859864 SRP125148 SRS2695456 GSM2859864 Bisulfite-Seq GENOMIC Reduced Representation Genomic DNA was isolated by phenol-chloroform extraction followed by ethanol precipitation. Two aliquots of DNA were separately digested with either TaqI or MspI restriction enzymes over night and subsequently pooled for all downstream processing. Digested DNA was end-repaired, A-tailed, and adapters containing methylated cytosines ligated on using the KAPA HyperPrep kit according to the manufacturers instructions. Adapter-ligated DNA was bisulfite converted over using the EpiTech Bisulfite kit (Qiagen, Inc). Bisulfite treated libraries were PCR amplified using indexed primers and the HiFi HotStart Uracil+ polymerase (KAPA Biosystems). Illumina HiSeq 2500 gds 302859864 GSM2859864 GEO Accession GSM2859864 SRX3401312 GSM2859865 SRP125148 SRS2695457 GSM2859865 Bisulfite-Seq GENOMIC Reduced Representation Genomic DNA was isolated by phenol-chloroform extraction followed by ethanol precipitation. Two aliquots of DNA were separately digested with either TaqI or MspI restriction enzymes over night and subsequently pooled for all downstream processing. Digested DNA was end-repaired, A-tailed, and adapters containing methylated cytosines ligated on using the KAPA HyperPrep kit according to the manufacturers instructions. Adapter-ligated DNA was bisulfite converted over using the EpiTech Bisulfite kit (Qiagen, Inc). Bisulfite treated libraries were PCR amplified using indexed primers and the HiFi HotStart Uracil+ polymerase (KAPA Biosystems). Illumina HiSeq 2500 gds 302859865 GSM2859865 GEO Accession GSM2859865