SRX3375258 GSM2844207 SRP124624 SRS2671686 GSM2844207 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844207 GSM2844207 GEO Accession GSM2844207 SRX3375259 GSM2844208 SRP124624 SRS2671687 GSM2844208 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844208 GSM2844208 GEO Accession GSM2844208 SRX3375260 GSM2844209 SRP124624 SRS2671692 GSM2844209 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844209 GSM2844209 GEO Accession GSM2844209 SRX3375261 GSM2844210 SRP124624 SRS2671689 GSM2844210 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844210 GSM2844210 GEO Accession GSM2844210 SRX3375262 GSM2844211 SRP124624 SRS2671693 GSM2844211 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844211 GSM2844211 GEO Accession GSM2844211 SRX3375263 GSM2844212 SRP124624 SRS2671688 GSM2844212 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844212 GSM2844212 GEO Accession GSM2844212 SRX3375264 GSM2844213 SRP124624 SRS2671691 GSM2844213 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844213 GSM2844213 GEO Accession GSM2844213 SRX3375265 GSM2844214 SRP124624 SRS2671695 GSM2844214 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844214 GSM2844214 GEO Accession GSM2844214 SRX3375266 GSM2844215 SRP124624 SRS2671696 GSM2844215 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844215 GSM2844215 GEO Accession GSM2844215 SRX3375267 GSM2844216 SRP124624 SRS2671690 GSM2844216 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844216 GSM2844216 GEO Accession GSM2844216 SRX3375268 GSM2844217 SRP124624 SRS2671694 GSM2844217 Bisulfite-Seq GENOMIC RANDOM Genomic DNA was sheared by DNA Shearing Instrument (30W%,3sON,9sOFF，20min) . The sheared DNAs size were between 400-500 bp and verified by 1.5% agrose gel electrophoresis. The quantity of the sheared DNAs was measured by Qubit 2.0(Invitrogen). For each sample, 4μg of genomic DNA was used for bisulfite-seq library preparation. Sheared DNAs were denaturated at 98℃ followed by bisulfite conversion. Then DNA was synthesized with DNA synthesis primer harboring tagging sequence. The DNA was purified and amplified and PCR products corresponding to 300-500 bps were purified, quantified and stored at -80 ℃ before sequencing. For high-throughput sequencing, the libraries were prepared following the manufacturer's instructions and applied to Illumina Hiseq 2000 system for 101 nt pair-end sequencing by Majorbio. Inc (Shanghai, China). Illumina HiSeq 2000 gds 302844217 GSM2844217 GEO Accession GSM2844217