SRX764932 GSM1551497 SRP050100 SRS749259 GSM1551497 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA samples were prepared using the Total RNA Purication Kit, TRK1001 (LC Science, Houston, TX), treated with RNase-free DNaseI following the manufacturer's procedure to remove contaminated genomic DNA. The quantity and purity of total RNA were analysed by Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA).Equal quantities of high-quantity RNA of each sample were blended for cDNA library construction. The poly (A) messgener RNA was isolated from the total RNA samples with oligo (dT) attached magnetic beads (Invitrogen). The mRNA was frangmented into short pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed to the first-strand cDNAs by random hexamer primers. Then the second-strand cDNAs were synthesized to construct the final cDNA library according to the instruction of mRNA-Seq sample preparation kit (Illumina, San Diego, USA) and the paired-end cDNA with a length of 300 bp (±50 bp). After that, the cDNA libraries were sequenced on the Illumina HiSeqTM 2000 platform at the LC-BIO (Hangzhou, China) following the manufacturer’s recommendations. Illumina HiSeq 2000 GEO Sample GSM1551497 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1551497 gds 301551497 GEO Accession GSM1551497 SRX764933 GSM1551498 SRP050100 SRS749260 GSM1551498 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA samples were prepared using the Total RNA Purication Kit, TRK1001 (LC Science, Houston, TX), treated with RNase-free DNaseI following the manufacturer's procedure to remove contaminated genomic DNA. The quantity and purity of total RNA were analysed by Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA).Equal quantities of high-quantity RNA of each sample were blended for cDNA library construction. The poly (A) messgener RNA was isolated from the total RNA samples with oligo (dT) attached magnetic beads (Invitrogen). The mRNA was frangmented into short pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed to the first-strand cDNAs by random hexamer primers. Then the second-strand cDNAs were synthesized to construct the final cDNA library according to the instruction of mRNA-Seq sample preparation kit (Illumina, San Diego, USA) and the paired-end cDNA with a length of 300 bp (±50 bp). After that, the cDNA libraries were sequenced on the Illumina HiSeqTM 2000 platform at the LC-BIO (Hangzhou, China) following the manufacturer’s recommendations. Illumina HiSeq 2000 GEO Sample GSM1551498 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1551498 gds 301551498 GEO Accession GSM1551498 SRX764934 GSM1551499 SRP050100 SRS749261 GSM1551499 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA samples were prepared using the Total RNA Purication Kit, TRK1001 (LC Science, Houston, TX), treated with RNase-free DNaseI following the manufacturer's procedure to remove contaminated genomic DNA. The quantity and purity of total RNA were analysed by Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA).Equal quantities of high-quantity RNA of each sample were blended for cDNA library construction. The poly (A) messgener RNA was isolated from the total RNA samples with oligo (dT) attached magnetic beads (Invitrogen). The mRNA was frangmented into short pieces using divalent cations under elevated temperature. The cleaved RNA fragments were reverse-transcribed to the first-strand cDNAs by random hexamer primers. Then the second-strand cDNAs were synthesized to construct the final cDNA library according to the instruction of mRNA-Seq sample preparation kit (Illumina, San Diego, USA) and the paired-end cDNA with a length of 300 bp (±50 bp). After that, the cDNA libraries were sequenced on the Illumina HiSeqTM 2000 platform at the LC-BIO (Hangzhou, China) following the manufacturer’s recommendations. Illumina HiSeq 2000 GEO Sample GSM1551499 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1551499 gds 301551499 GEO Accession GSM1551499