To gain further insights into the molecular mechanisms of attractive quality formation during AG development between PP and SP, gene ontology enrichment analysis was performed to determine which kind of GO term DEGs were mainly enriched in. Finally, 2,345 DEGs were annotated to 3,588 GO terms. The functions of 854 DEGs with higher expression in SP than in PP AG were assigned to biological process, cellular component, and molecular function. In biological process, SRP-dependent cotranslational protein targeting to membrane (41 DEGs, GO: 0006614), nuclear-transcribed mRNA catabolic process, nonsense-mediated decay (42 DEGs, GO: 0000184), cotranslational protein targeting to membrane (41 DEGs, GO: 0006613), protein targeting to ER (41 DEGs, GO: 0045047), and establishment of protein localization to endoplasmic reticulum (41 DEGs, GO: 0072599) were the most prominent terms. Cytosolic ribosome (44 DEGs, GO: 0022626) was the most prominent term within the cellular component followed by ribosomal subunit (46 DEGs, GO: 0044391), cytosolic large ribosomal subunit (27 DEGs, GO: 0022625), ribosome (49 DEGs, GO: 0005840), and extracellular exosome (96 DEGs, GO: 0070062). In molecular function, most of the annotated unique sequences were assigned to phosphatidylserine binding (13 DEGs, GO: 0001786), structural constituent of ribosome (48 DEGs, GO: 0003735), and modified amino acid binding (14 DEGs, GO: 0072341).