Total RNA was isolated from plant tissue using the BIOCUBE Plant RNA prep kit (Biocube System Inc., Korea) according to the manufacturers instructions. The quality of the extracted total RNA was measured using Bioanalyzer 2100 (Agilent, SantaClara, CA, USA). Total RNA with an RNA Integrity Number (RIN)value greater than or equal to 7 was used for library preparation. RNA sequencing libraries were prepared using the TruSeq Stranded mRNA LTS ample Prep Kit according to the manufacturers instructions (Illumina, SanDiego, CA, USA). A total of nine libraries with the irrespective indices were paired-ends equenced using the NovaSeq 6000 system (Illumina).