## 3 Results
### Transcriptome Analysis
#### 3.2.2 Annotation of Unigenes
Meanwhile, 7% of unigenes sequences were between 30% and 60% (Figure 4c). Further analysis using the InterPro database, the distribution of enzyme classes identified in flower part transcripts of the three study species (Figure 5) and the protein class transferases dominated the enzyme class. In contrast, the enzyme hydrolases contributed to the second-highest distribution of protein classes. In _H. agauvanensis_, a total of 25,186 unigenes were classified into 4018 KEGG orthogonal groups and 413 KEGG pathways. The distribution of these KEGG orthogonal groups with different cell function categories is shown in Figure 6a. The protein families: genetic information process group (5044 unigenes) and the genetic information processing group (4785 unigenes) were the largest KEGG orthological groups. For _H. lacunosa_, a total of 24,703 unigenes were classified into 3706 KEGG orthological groups covering 413 KEGG pathways. The orthology group of the protein families: genetic information process group (5964 unigenes) and the genetic information processing group (4667 unigenes) were found  
Figure 5: Distribution of protein classes identified on flower part of three species _Hoya_; (**a**) _H. cagauvanensis_; (**b**) _H. lacunosa_, and (**c**) _H. coriacea_. Each area of the pie represents a value from the unigen. to be the largest orthogonal groups of KEGG (Figure 6b). While for _H. coriacea_, a total of 23,752 unigenes were classified into 3643 KEGG orthogonal groups and 408 KEGG pathways through the GhostKOALA program. The greatest orthogonal group of KEGG were discovered to be the protein families: genetic information process group (4873 unigenes) and the genetic information processing group (4682 unigenes) (Figure 6c). Based on the unigene InterPro database, samples were classified into functions according to the classification from the GO database using Blast2 GO software. A total of 98,203 GO _H. cagayanensis_, 335,81 GO _H. lacunosa_, and 56,914 GO _H. coriacea_ terms were matched in the unigene database. The unigenes distributions successfully matched with the GO database showed elevations in the biological processes, between levels 8-10 for _H. cagayanensis_; levels 7-9 for _H. lacunosa_; and levels 6-9 for _H. coriacea_. For molecular functional processes, the highest concentration levels were 5-6 for all species study. Meanwhile, in the cellular component, all species are highest in levels 5-6 (Figure 7). Figure 6: Classification of KEGG function on the flower part of three species _Hova_; (**a**) _H. cagayanensis_; (**b**) _H. lacunosa_, and (**c**) _H. coriacea_. Unigen functions were identified and categorized using GhostKOALA software over the entire eukaryotic database. Furthermore, for the category of BP in this study, the main term GO is 'organic substance metabolic process' for all species study; _H. cagayanensis_ (GO:0071074; 43,785 unigenes), _H. lacunosa_ (GO:0071704; 17,101 unigenes), and _H. coriacea_ (GO:0071704; 28,533 unigenes). Meanwhile, the category of MF was dominated by the term GO 'organic cyclic compound  
Figure 7: GO category classification in three _Hoya_ species **(a)**_H. cagayanensis;_ **(b)**_H. lacunosa;_ and **(c)**_H. coriaca_: Division of terms in GO from various levels across the categories of cellular component (CC), molecular function (MF) and biological process (BP). binding' with _H. cagayamensis_ (GO:0097159; 30,966 unigenes), _H. lacunosa_ (GO:0097159; 11,901 unigenes), and _H. coriacea_ (GO:0097159; 19,357 unigenes). Next, the CC category was monopolized by the term GO 'Intracellular anatomical structure' with _H. cagayamensis_ (GO:0005622; 36,653 unigenes), _H. lacunosa_ (GO:0005622; 15,785 unigenes), and _H. coriacea_ (GO:0005622; 29,351 unigenes) (Figure 8).