Fig. 1

A pipeline summarizing the guilt-by-association and functional annotation approaches used in this study. The 344 PUFs measured via LC–MS/MS were subjected to co-expression and differential expression analyses. Structural modeling with SwissModel was used to determine structural templates which best fit a PUFs protein sequence. Domain and function prediction were performed using InterProScan, eggNOG-mapper, and PANNZER2. Phylogenetic coevolution analysis was used to test for coevolution. Gene trees were generated from a homology search in BLAST, followed by alignment of the top 200 hits and tree construction using FastTree. Regulatory information based on shared operons was extracted from the DOOR database