Expansion and contractions of orthologous gene families were determined using CAFÉ 2.2 (Computational Analysis of gene Family Evolution). The program uses a birth and death process to model gene gain and loss over a phylogeny. Large changes in gene family size in a phylogeny were tested by calculating p-values on each branch using the Viterbi method with a randomly generated likelihood distribution. This method calculates exact p-values for transitions between the parent and child family sizes for all branches of the phylogenetic tree. Enrichment of Gene Ontology terms for G. elata expanded gene families were summarized and visualized using REVIGO (small list, similarity (0.5), SimRel similarity measure).