Gene expression levels were directly computed from the read alignment results for each replicate in the different groups. Standard RPKM values31 (reads per kilobase of exon model per million mapped reads) were extracted for each gene covered by the sequencing data and each replicate used in this study.