The libraries were prepared by Illumina TruSeq™ RNA sample prep kit technology. The complete mRNA was obtained by the enrichment of magnetic beads with oligo (dT), and then the mRNA was randomly broken by adding fragmentation buffer. The small fragments of 300 bp were screened out by the magnetic beads. First-strand cDNA was synthesized using random hexamer primer and reverse transcriptase, and then the second strand was synthesized. End repair mix was added to make the sticky end of the cDNA into a flat end, and then a “A” base was added to the 3′ end.
The libraries were sequenced on an Illumina NovaSeq 6000 platform according to manufactures’ instructions. After sequencing, software seqprep and sickle were used to remove the linker sequence, low quality reads, any anonymous nucleotides greater than 10%, and adapter to obtain clean data. The clean data was assembly using trinity software (Grabherr et al., 2011 (link)). The software transrate and busco were used to evaluate and optimize the assembly results. Finally, we obtained the unigenes database by aggregating the assembled sequences.