SRX3892932 GSM3082818 SRP137711 SRS3130739 GSM3082818 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was assessed by electrophoresis on a denaturing agarose gel and quantified by NanoDrop spectrophotometer (NanoDrop, USA) and Agilent 2100 Bioanalyzer. The total RNA samples (3 mg) were treated with the epicentre Ribo-ZeroTM Kit to remove rRNA before constructing the RNA-seq libraries. Then the rRNA depleted RNA samples were fragmented and used to synthesize first- and second-strand complementary DNA (cDNA) with random hexamer primers, dNTPs and DNA Polymerase I. The cDNA fragments were cleaned and concentrated using AMPure XP beads, then the ends were repaired and modified with T4 DNA polymerase Klenow DNA polymerase to add A and adapter at the 3' end of the DNA fragments. The ligated cDNA products were purified with AMPure XP beads and treated with uracil DNA glycosylase to remove the second-strand cDNA. Purified first-strand cDNA was subjected to PCR amplification, and the libraries were quality controlled with a Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) and sequenced by Illumina HiSeq 3000 (LC Sciences, USA) on a 125 bp paired-end run. Illumina HiSeq 3000 gds 303082818 GSM3082818 GEO Accession GSM3082818 SRX3892933 GSM3082819 SRP137711 SRS3130744 GSM3082819 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was assessed by electrophoresis on a denaturing agarose gel and quantified by NanoDrop spectrophotometer (NanoDrop, USA) and Agilent 2100 Bioanalyzer. The total RNA samples (3 mg) were treated with the epicentre Ribo-ZeroTM Kit to remove rRNA before constructing the RNA-seq libraries. Then the rRNA depleted RNA samples were fragmented and used to synthesize first- and second-strand complementary DNA (cDNA) with random hexamer primers, dNTPs and DNA Polymerase I. The cDNA fragments were cleaned and concentrated using AMPure XP beads, then the ends were repaired and modified with T4 DNA polymerase Klenow DNA polymerase to add A and adapter at the 3' end of the DNA fragments. The ligated cDNA products were purified with AMPure XP beads and treated with uracil DNA glycosylase to remove the second-strand cDNA. Purified first-strand cDNA was subjected to PCR amplification, and the libraries were quality controlled with a Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) and sequenced by Illumina HiSeq 3000 (LC Sciences, USA) on a 125 bp paired-end run. Illumina HiSeq 3000 gds 303082819 GSM3082819 GEO Accession GSM3082819 SRX3892934 GSM3082820 SRP137711 SRS3130743 GSM3082820 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was assessed by electrophoresis on a denaturing agarose gel and quantified by NanoDrop spectrophotometer (NanoDrop, USA) and Agilent 2100 Bioanalyzer. The total RNA samples (3 mg) were treated with the epicentre Ribo-ZeroTM Kit to remove rRNA before constructing the RNA-seq libraries. Then the rRNA depleted RNA samples were fragmented and used to synthesize first- and second-strand complementary DNA (cDNA) with random hexamer primers, dNTPs and DNA Polymerase I. The cDNA fragments were cleaned and concentrated using AMPure XP beads, then the ends were repaired and modified with T4 DNA polymerase Klenow DNA polymerase to add A and adapter at the 3' end of the DNA fragments. The ligated cDNA products were purified with AMPure XP beads and treated with uracil DNA glycosylase to remove the second-strand cDNA. Purified first-strand cDNA was subjected to PCR amplification, and the libraries were quality controlled with a Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) and sequenced by Illumina HiSeq 3000 (LC Sciences, USA) on a 125 bp paired-end run. Illumina HiSeq 3000 gds 303082820 GSM3082820 GEO Accession GSM3082820 SRX3892935 GSM3082821 SRP137711 SRS3130742 GSM3082821 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was assessed by electrophoresis on a denaturing agarose gel and quantified by NanoDrop spectrophotometer (NanoDrop, USA) and Agilent 2100 Bioanalyzer. The total RNA samples (3 mg) were treated with the epicentre Ribo-ZeroTM Kit to remove rRNA before constructing the RNA-seq libraries. Then the rRNA depleted RNA samples were fragmented and used to synthesize first- and second-strand complementary DNA (cDNA) with random hexamer primers, dNTPs and DNA Polymerase I. The cDNA fragments were cleaned and concentrated using AMPure XP beads, then the ends were repaired and modified with T4 DNA polymerase Klenow DNA polymerase to add A and adapter at the 3' end of the DNA fragments. The ligated cDNA products were purified with AMPure XP beads and treated with uracil DNA glycosylase to remove the second-strand cDNA. Purified first-strand cDNA was subjected to PCR amplification, and the libraries were quality controlled with a Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) and sequenced by Illumina HiSeq 3000 (LC Sciences, USA) on a 125 bp paired-end run. Illumina HiSeq 3000 gds 303082821 GSM3082821 GEO Accession GSM3082821 SRX3892936 GSM3082822 SRP137711 SRS3130741 GSM3082822 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was assessed by electrophoresis on a denaturing agarose gel and quantified by NanoDrop spectrophotometer (NanoDrop, USA) and Agilent 2100 Bioanalyzer. The total RNA samples (3 mg) were treated with the epicentre Ribo-ZeroTM Kit to remove rRNA before constructing the RNA-seq libraries. Then the rRNA depleted RNA samples were fragmented and used to synthesize first- and second-strand complementary DNA (cDNA) with random hexamer primers, dNTPs and DNA Polymerase I. The cDNA fragments were cleaned and concentrated using AMPure XP beads, then the ends were repaired and modified with T4 DNA polymerase Klenow DNA polymerase to add A and adapter at the 3' end of the DNA fragments. The ligated cDNA products were purified with AMPure XP beads and treated with uracil DNA glycosylase to remove the second-strand cDNA. Purified first-strand cDNA was subjected to PCR amplification, and the libraries were quality controlled with a Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) and sequenced by Illumina HiSeq 3000 (LC Sciences, USA) on a 125 bp paired-end run. Illumina HiSeq 3000 gds 303082822 GSM3082822 GEO Accession GSM3082822 SRX3892937 GSM3082823 SRP137711 SRS3130740 GSM3082823 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was assessed by electrophoresis on a denaturing agarose gel and quantified by NanoDrop spectrophotometer (NanoDrop, USA) and Agilent 2100 Bioanalyzer. The total RNA samples (3 mg) were treated with the epicentre Ribo-ZeroTM Kit to remove rRNA before constructing the RNA-seq libraries. Then the rRNA depleted RNA samples were fragmented and used to synthesize first- and second-strand complementary DNA (cDNA) with random hexamer primers, dNTPs and DNA Polymerase I. The cDNA fragments were cleaned and concentrated using AMPure XP beads, then the ends were repaired and modified with T4 DNA polymerase Klenow DNA polymerase to add A and adapter at the 3' end of the DNA fragments. The ligated cDNA products were purified with AMPure XP beads and treated with uracil DNA glycosylase to remove the second-strand cDNA. Purified first-strand cDNA was subjected to PCR amplification, and the libraries were quality controlled with a Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) and sequenced by Illumina HiSeq 3000 (LC Sciences, USA) on a 125 bp paired-end run. Illumina HiSeq 3000 gds 303082823 GSM3082823 GEO Accession GSM3082823