The quality of RNA-seq reads were checked using Trimmomatic (Bolger et al., 2014) and the low quality of reads were lter out. Then, Hisat 2 was used to align clean RNA-seq reads to reference genome. The RSEM software (Li and Dewey, 2011) was used to calculate gene expression values. The DESeq2 package (Anders and Huber, 2010) was used for identify differential expressed genes. Total RNA was extracted for qRT-PCR experiment according to the protocol of RNAprep pure plant kit (TIANGEN, China).
The reverse transcription from RNA to cDNA was performed using cDNA Synthesis Kit (TIANGEN, China). Peanut Actin gene was used as the internal control and relative gene expression was calculated by the 2-ΔΔCt method.