SRX831898 GSM1577738 SRP051737 SRS810210 GSM1577738 RNA-Seq TRANSCRIPTOMIC cDNA RNA was extracted with Aurum™ Total RNA Mini Kit (Bio-Rad) A modified dUTP method was used for preparation of RNA-Seq libraries[27]. Briefly, mRNA was purified from total RNA by Dynabeads® mRNA Purification Kit (Invitrogen) and reverse-transcribed with SuperScript® III (Invitrogen). The second strand was created by DNA Polymerase 1 (NEB) in the presence of RNase H and a dNTP mix containing dUTP instead dTTP. The cDNA was fragmented on a Covaris sonicator, and end repair and A-tailing were performed using the Quick Blunting kit and exo- Klenow fragment of DNA polymerase I, respectively (NEB). Illumina genomic adapters were ligated, and the sample was size-fractionated (~170–400 bp) on E-Gel EX 2% agarose (Invitrogen). After uracil-dna glycosylase (UDG) treatment (NEB) and a final PCR amplification step with Illumina primers (18 cycles), the cDNA libraries were size-fractionated (~175-400 bp) on a 2% agarose gel, quantified using a Qubit 2.0 fluorometer and dsDNA HS kit (Invitrogen), and sent to sequencing on a Illumina HiSeq 2500 sequencing system Illumina HiSeq 2500 GEO Sample GSM1577738 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1577738 gds 301577738 GEO Accession GSM1577738 SRX831899 GSM1577739 SRP051737 SRS810211 GSM1577739 RNA-Seq TRANSCRIPTOMIC cDNA RNA was extracted with Aurum™ Total RNA Mini Kit (Bio-Rad) A modified dUTP method was used for preparation of RNA-Seq libraries[27]. Briefly, mRNA was purified from total RNA by Dynabeads® mRNA Purification Kit (Invitrogen) and reverse-transcribed with SuperScript® III (Invitrogen). The second strand was created by DNA Polymerase 1 (NEB) in the presence of RNase H and a dNTP mix containing dUTP instead dTTP. The cDNA was fragmented on a Covaris sonicator, and end repair and A-tailing were performed using the Quick Blunting kit and exo- Klenow fragment of DNA polymerase I, respectively (NEB). Illumina genomic adapters were ligated, and the sample was size-fractionated (~170–400 bp) on E-Gel EX 2% agarose (Invitrogen). After uracil-dna glycosylase (UDG) treatment (NEB) and a final PCR amplification step with Illumina primers (18 cycles), the cDNA libraries were size-fractionated (~175-400 bp) on a 2% agarose gel, quantified using a Qubit 2.0 fluorometer and dsDNA HS kit (Invitrogen), and sent to sequencing on a Illumina HiSeq 2500 sequencing system Illumina HiSeq 2500 GEO Sample GSM1577739 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1577739 gds 301577739 GEO Accession GSM1577739 SRX831900 GSM1577740 SRP051737 SRS810212 GSM1577740 RNA-Seq TRANSCRIPTOMIC cDNA RNA was extracted with Aurum™ Total RNA Mini Kit (Bio-Rad) A modified dUTP method was used for preparation of RNA-Seq libraries[27]. Briefly, mRNA was purified from total RNA by Dynabeads® mRNA Purification Kit (Invitrogen) and reverse-transcribed with SuperScript® III (Invitrogen). The second strand was created by DNA Polymerase 1 (NEB) in the presence of RNase H and a dNTP mix containing dUTP instead dTTP. The cDNA was fragmented on a Covaris sonicator, and end repair and A-tailing were performed using the Quick Blunting kit and exo- Klenow fragment of DNA polymerase I, respectively (NEB). Illumina genomic adapters were ligated, and the sample was size-fractionated (~170–400 bp) on E-Gel EX 2% agarose (Invitrogen). After uracil-dna glycosylase (UDG) treatment (NEB) and a final PCR amplification step with Illumina primers (18 cycles), the cDNA libraries were size-fractionated (~175-400 bp) on a 2% agarose gel, quantified using a Qubit 2.0 fluorometer and dsDNA HS kit (Invitrogen), and sent to sequencing on a Illumina HiSeq 2500 sequencing system Illumina HiSeq 2500 GEO Sample GSM1577740 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1577740 gds 301577740 GEO Accession GSM1577740 SRX831901 GSM1577741 SRP051737 SRS810213 GSM1577741 RNA-Seq TRANSCRIPTOMIC cDNA RNA was extracted with Aurum™ Total RNA Mini Kit (Bio-Rad) A modified dUTP method was used for preparation of RNA-Seq libraries[27]. Briefly, mRNA was purified from total RNA by Dynabeads® mRNA Purification Kit (Invitrogen) and reverse-transcribed with SuperScript® III (Invitrogen). The second strand was created by DNA Polymerase 1 (NEB) in the presence of RNase H and a dNTP mix containing dUTP instead dTTP. The cDNA was fragmented on a Covaris sonicator, and end repair and A-tailing were performed using the Quick Blunting kit and exo- Klenow fragment of DNA polymerase I, respectively (NEB). Illumina genomic adapters were ligated, and the sample was size-fractionated (~170–400 bp) on E-Gel EX 2% agarose (Invitrogen). After uracil-dna glycosylase (UDG) treatment (NEB) and a final PCR amplification step with Illumina primers (18 cycles), the cDNA libraries were size-fractionated (~175-400 bp) on a 2% agarose gel, quantified using a Qubit 2.0 fluorometer and dsDNA HS kit (Invitrogen), and sent to sequencing on a Illumina HiSeq 2500 sequencing system Illumina HiSeq 2500 GEO Sample GSM1577741 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1577741 gds 301577741 GEO Accession GSM1577741 SRX831902 GSM1577742 SRP051737 SRS810214 GSM1577742 RNA-Seq TRANSCRIPTOMIC cDNA RNA was extracted with Aurum™ Total RNA Mini Kit (Bio-Rad) A modified dUTP method was used for preparation of RNA-Seq libraries[27]. Briefly, mRNA was purified from total RNA by Dynabeads® mRNA Purification Kit (Invitrogen) and reverse-transcribed with SuperScript® III (Invitrogen). The second strand was created by DNA Polymerase 1 (NEB) in the presence of RNase H and a dNTP mix containing dUTP instead dTTP. The cDNA was fragmented on a Covaris sonicator, and end repair and A-tailing were performed using the Quick Blunting kit and exo- Klenow fragment of DNA polymerase I, respectively (NEB). Illumina genomic adapters were ligated, and the sample was size-fractionated (~170–400 bp) on E-Gel EX 2% agarose (Invitrogen). After uracil-dna glycosylase (UDG) treatment (NEB) and a final PCR amplification step with Illumina primers (18 cycles), the cDNA libraries were size-fractionated (~175-400 bp) on a 2% agarose gel, quantified using a Qubit 2.0 fluorometer and dsDNA HS kit (Invitrogen), and sent to sequencing on a Illumina HiSeq 2500 sequencing system Illumina HiSeq 2500 GEO Sample GSM1577742 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1577742 gds 301577742 GEO Accession GSM1577742 SRX831903 GSM1577743 SRP051737 SRS810215 GSM1577743 RNA-Seq TRANSCRIPTOMIC cDNA RNA was extracted with Aurum™ Total RNA Mini Kit (Bio-Rad) A modified dUTP method was used for preparation of RNA-Seq libraries[27]. Briefly, mRNA was purified from total RNA by Dynabeads® mRNA Purification Kit (Invitrogen) and reverse-transcribed with SuperScript® III (Invitrogen). The second strand was created by DNA Polymerase 1 (NEB) in the presence of RNase H and a dNTP mix containing dUTP instead dTTP. The cDNA was fragmented on a Covaris sonicator, and end repair and A-tailing were performed using the Quick Blunting kit and exo- Klenow fragment of DNA polymerase I, respectively (NEB). Illumina genomic adapters were ligated, and the sample was size-fractionated (~170–400 bp) on E-Gel EX 2% agarose (Invitrogen). After uracil-dna glycosylase (UDG) treatment (NEB) and a final PCR amplification step with Illumina primers (18 cycles), the cDNA libraries were size-fractionated (~175-400 bp) on a 2% agarose gel, quantified using a Qubit 2.0 fluorometer and dsDNA HS kit (Invitrogen), and sent to sequencing on a Illumina HiSeq 2500 sequencing system Illumina HiSeq 2500 GEO Sample GSM1577743 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1577743 gds 301577743 GEO Accession GSM1577743 SRX831904 GSM1577744 SRP051737 SRS810216 GSM1577744 RNA-Seq TRANSCRIPTOMIC cDNA RNA was extracted with Aurum™ Total RNA Mini Kit (Bio-Rad) A modified dUTP method was used for preparation of RNA-Seq libraries[27]. Briefly, mRNA was purified from total RNA by Dynabeads® mRNA Purification Kit (Invitrogen) and reverse-transcribed with SuperScript® III (Invitrogen). The second strand was created by DNA Polymerase 1 (NEB) in the presence of RNase H and a dNTP mix containing dUTP instead dTTP. The cDNA was fragmented on a Covaris sonicator, and end repair and A-tailing were performed using the Quick Blunting kit and exo- Klenow fragment of DNA polymerase I, respectively (NEB). Illumina genomic adapters were ligated, and the sample was size-fractionated (~170–400 bp) on E-Gel EX 2% agarose (Invitrogen). After uracil-dna glycosylase (UDG) treatment (NEB) and a final PCR amplification step with Illumina primers (18 cycles), the cDNA libraries were size-fractionated (~175-400 bp) on a 2% agarose gel, quantified using a Qubit 2.0 fluorometer and dsDNA HS kit (Invitrogen), and sent to sequencing on a Illumina HiSeq 2500 sequencing system Illumina HiSeq 2500 GEO Sample GSM1577744 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1577744 gds 301577744 GEO Accession GSM1577744