WT-A-R1
RNAseq libraries were constructed by QIAquick PCR kit following manufacturer's protocols. Then highthroughput sequencing was performed using Illumina HiSeq 2000. The clean reads of 6 samples were separately mapped to the reference genome and reference genes by using Tophat1.3.2 (Trapnell et al., 2009) and Bowtie2, respectively. The FPKM value (fragments per kilobase of transcript per million mapped reads)) of each gene in each sample was calculated by RSEM V1.2.15 (Li & Dewey, 2011), which was used to compare the differences of gene expression. The differential expression of gene in each sample was analyzed by edgeR, P value and FDR (False Discovery Rate) were calculated. GO annotation and GO enrichment analysis of the differentially expressed genes were performed by topGO. KEGGPathway analysis was performed by KOBAS, and the significant enriched pathways was found by hypergeometric test.