SRX22342361 GSM7880313_r1 SRP469926 PRJNA1035046 SRS19391229 GSM7880313 GSM7880313 Bisulfite-Seq TRANSCRIPTOMIC RANDOM Total RNAs were extracted using TRIzol Reagent (Invitrogen, cat. NO 15596026）. DNA digestion was carried out after RNA extraction by DNaseI. 2 μg total RNAs were used for mRNA enrichment by KAPA mRNA Capture Kit (Roche, KK8441). The conversion of unmethylated cytosines to thymidines of enriched mRNA was carried out using the EZ RNA Methylation Kit (Zymo Research, Cat. No. R5001) . Sequencing libraries were prepared using KCTM Stranded mRNA Library Prep Kit for Illumina® (Catalog NO. DR08402, Wuhan Seqhealth Co., Ltd. China) following the manufacturer's instruction. PCR products corresponding to 200-500 bps were enriched, quantified and finally sequenced on Novaseq 6000 sequencer (Illumina) with PE150 model. Illumina NovaSeq 6000 SRX22342362 GSM7880314_r1 SRP469926 PRJNA1035046 SRS19391231 GSM7880314 GSM7880314 Bisulfite-Seq TRANSCRIPTOMIC RANDOM Total RNAs were extracted using TRIzol Reagent (Invitrogen, cat. NO 15596026）. DNA digestion was carried out after RNA extraction by DNaseI. 2 μg total RNAs were used for mRNA enrichment by KAPA mRNA Capture Kit (Roche, KK8441). The conversion of unmethylated cytosines to thymidines of enriched mRNA was carried out using the EZ RNA Methylation Kit (Zymo Research, Cat. No. R5001) . Sequencing libraries were prepared using KCTM Stranded mRNA Library Prep Kit for Illumina® (Catalog NO. DR08402, Wuhan Seqhealth Co., Ltd. China) following the manufacturer's instruction. PCR products corresponding to 200-500 bps were enriched, quantified and finally sequenced on Novaseq 6000 sequencer (Illumina) with PE150 model. Illumina NovaSeq 6000 SRX22342363 GSM7880315_r1 SRP469926 PRJNA1035046 SRS19391230 GSM7880315 GSM7880315 Bisulfite-Seq TRANSCRIPTOMIC RANDOM Total RNAs were extracted using TRIzol Reagent (Invitrogen, cat. NO 15596026）. DNA digestion was carried out after RNA extraction by DNaseI. 2 μg total RNAs were used for mRNA enrichment by KAPA mRNA Capture Kit (Roche, KK8441). The conversion of unmethylated cytosines to thymidines of enriched mRNA was carried out using the EZ RNA Methylation Kit (Zymo Research, Cat. No. R5001) . Sequencing libraries were prepared using KCTM Stranded mRNA Library Prep Kit for Illumina® (Catalog NO. DR08402, Wuhan Seqhealth Co., Ltd. China) following the manufacturer's instruction. PCR products corresponding to 200-500 bps were enriched, quantified and finally sequenced on Novaseq 6000 sequencer (Illumina) with PE150 model. Illumina NovaSeq 6000 SRX22342364 GSM7880316_r1 SRP469926 PRJNA1035046 SRS19391233 GSM7880316 GSM7880316 Bisulfite-Seq TRANSCRIPTOMIC RANDOM Total RNAs were extracted using TRIzol Reagent (Invitrogen, cat. NO 15596026）. DNA digestion was carried out after RNA extraction by DNaseI. 2 μg total RNAs were used for mRNA enrichment by KAPA mRNA Capture Kit (Roche, KK8441). The conversion of unmethylated cytosines to thymidines of enriched mRNA was carried out using the EZ RNA Methylation Kit (Zymo Research, Cat. No. R5001) . Sequencing libraries were prepared using KCTM Stranded mRNA Library Prep Kit for Illumina® (Catalog NO. DR08402, Wuhan Seqhealth Co., Ltd. China) following the manufacturer's instruction. PCR products corresponding to 200-500 bps were enriched, quantified and finally sequenced on Novaseq 6000 sequencer (Illumina) with PE150 model. Illumina NovaSeq 6000 SRX22342365 GSM7880317_r1 SRP469926 PRJNA1035046 SRS19391232 GSM7880317 GSM7880317 Bisulfite-Seq TRANSCRIPTOMIC RANDOM Total RNAs were extracted using TRIzol Reagent (Invitrogen, cat. NO 15596026）. DNA digestion was carried out after RNA extraction by DNaseI. 2 μg total RNAs were used for mRNA enrichment by KAPA mRNA Capture Kit (Roche, KK8441). The conversion of unmethylated cytosines to thymidines of enriched mRNA was carried out using the EZ RNA Methylation Kit (Zymo Research, Cat. No. R5001) . Sequencing libraries were prepared using KCTM Stranded mRNA Library Prep Kit for Illumina® (Catalog NO. DR08402, Wuhan Seqhealth Co., Ltd. China) following the manufacturer's instruction. PCR products corresponding to 200-500 bps were enriched, quantified and finally sequenced on Novaseq 6000 sequencer (Illumina) with PE150 model. Illumina NovaSeq 6000 SRX22342366 GSM7880318_r1 SRP469926 PRJNA1035046 SRS19391234 GSM7880318 GSM7880318 Bisulfite-Seq TRANSCRIPTOMIC RANDOM Total RNAs were extracted using TRIzol Reagent (Invitrogen, cat. NO 15596026）. DNA digestion was carried out after RNA extraction by DNaseI. 2 μg total RNAs were used for mRNA enrichment by KAPA mRNA Capture Kit (Roche, KK8441). The conversion of unmethylated cytosines to thymidines of enriched mRNA was carried out using the EZ RNA Methylation Kit (Zymo Research, Cat. No. R5001) . Sequencing libraries were prepared using KCTM Stranded mRNA Library Prep Kit for Illumina® (Catalog NO. DR08402, Wuhan Seqhealth Co., Ltd. China) following the manufacturer's instruction. PCR products corresponding to 200-500 bps were enriched, quantified and finally sequenced on Novaseq 6000 sequencer (Illumina) with PE150 model. Illumina NovaSeq 6000