Gene functional annotation was performed by sequence comparison with public databases. The unigenes of the de novo assembly were searched against the NCBI nonredundant (NR), Swiss-Prot, KEGG, and COG protein databases by using BLASTX with a cut-off E value of 10−5.