SRX5106021 GSM3505035 SRP172870 SRS4115666 GSM3505035 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505035 GSM3505035 GEO Accession GSM3505035 SRX5106022 GSM3505036 SRP172870 SRS4115667 GSM3505036 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505036 GSM3505036 GEO Accession GSM3505036 SRX5106023 GSM3505037 SRP172870 SRS4115670 GSM3505037 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505037 GSM3505037 GEO Accession GSM3505037 SRX5106024 GSM3505038 SRP172870 SRS4115668 GSM3505038 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505038 GSM3505038 GEO Accession GSM3505038 SRX5106025 GSM3505039 SRP172870 SRS4115669 GSM3505039 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505039 GSM3505039 GEO Accession GSM3505039 SRX5106026 GSM3505040 SRP172870 SRS4115671 GSM3505040 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505040 GSM3505040 GEO Accession GSM3505040 SRX5106027 GSM3505041 SRP172870 SRS4115672 GSM3505041 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505041 GSM3505041 GEO Accession GSM3505041 SRX5106028 GSM3505042 SRP172870 SRS4115673 GSM3505042 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505042 GSM3505042 GEO Accession GSM3505042 SRX5106029 GSM3505043 SRP172870 SRS4115674 GSM3505043 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505043 GSM3505043 GEO Accession GSM3505043 SRX5106030 GSM3505044 SRP172870 SRS4115675 GSM3505044 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505044 GSM3505044 GEO Accession GSM3505044 SRX5106031 GSM3505045 SRP172870 SRS4115676 GSM3505045 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505045 GSM3505045 GEO Accession GSM3505045 SRX5106032 GSM3505046 SRP172870 SRS4115677 GSM3505046 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The total RNA quantity and purity were analysis of Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number >7.0. Approximately 10 ug of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligoattached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed with Random Hexamers to create the cDNA library, and then synthesized the second strands with addition of buffer, dNTPs, RNaseH and DNA polymerase to be double-stranded DNA. Double-stranded DNA fragments with overhangs were repaired to have blunt ends by Klenow DNA polymerase and T4 DNA polymerase and then added A-overhang as well as adapters before purified by AMPure XP beads. Purified DNA fragments were amplified by PCR before purification again to generate the final library for sequencing. The paired-end sequencing was performed on an Illumina Hiseq4000 following the vendor's recommended protocol and generated reads of 2x150bp (150 PE) Illumina HiSeq 4000 gds 303505046 GSM3505046 GEO Accession GSM3505046