SRX13183627 GSM5696934_r1 SRP347099 PRJNA782319 SRS11112852 GSM5696934 GSM5696934 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183628 GSM5696935_r1 SRP347099 PRJNA782319 SRS11112853 GSM5696935 GSM5696935 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183629 GSM5696936_r1 SRP347099 PRJNA782319 SRS11112854 GSM5696936 GSM5696936 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183630 GSM5696937_r1 SRP347099 PRJNA782319 SRS11112855 GSM5696937 GSM5696937 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183631 GSM5696938_r1 SRP347099 PRJNA782319 SRS11112856 GSM5696938 GSM5696938 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183632 GSM5696939_r1 SRP347099 PRJNA782319 SRS11112857 GSM5696939 GSM5696939 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183633 GSM5696940_r1 SRP347099 PRJNA782319 SRS11112858 GSM5696940 GSM5696940 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183634 GSM5696941_r1 SRP347099 PRJNA782319 SRS11112859 GSM5696941 GSM5696941 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183635 GSM5696942_r1 SRP347099 PRJNA782319 SRS11112860 GSM5696942 GSM5696942 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183636 GSM5696943_r1 SRP347099 PRJNA782319 SRS11112861 GSM5696943 GSM5696943 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183637 GSM5696944_r1 SRP347099 PRJNA782319 SRS11112862 GSM5696944 GSM5696944 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000 SRX13183638 GSM5696945_r1 SRP347099 PRJNA782319 SRS11112863 GSM5696945 GSM5696945 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). The RNA integrity was assessed by Agilent 2100 with RIN number >7.0. Poly(A) RNA is purified from total RNA(5ug) using poly-T oligo-attached magnetic beads using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using divalent cations under high temperature. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA, which were next used to synthesise U-labeled second-stranded DNAs with E. coli DNA polymerase I, RNase H and dUTP. An A-base is then added to the blunt ends of each strand, then ligated to modified Illumina multiplex barcode adapters, which including custom Unique Molecular Identifiers for minimizing sequence-dependent bias and amplification noise according to (Shiroguchi et al. 2012), and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme treatment of the U-labeled second-stranded DNAs, The ligated products are amplified with PCR by the following conditions: initial denaturation at 95℃ for 3 min; 8cycles of denaturation at 98℃ for 15 sec, annealing at 60℃ for 15 sec, and extension at 72℃ for 30 sec; and then final extension at 72℃ for 5 min. The average insert size for the final cDNA library was 300 bp (±50 bp). At last, we performed the paired-end sequencing on an Illumina Novaseq TM 6000 following the vendor's recommended protocol. Illumina NovaSeq 6000