CRAL_TRIO domain proteins are known to bind small lipophilic molecules such as retinal, inositol and Vitamin E and include such gene family members as PINTA, α-tocopherol transfer (ATT) proteins, retinoid binding proteins, and clavesins. In insects, very little is known about either the molecular evolution of this family of proteins or their ligand specificity. Here we characterize insect CRAL_TRIO domain proteins and present the first insect CRAL_TRIO protein phylogeny constructed by performing reciprocal BLAST searches of the reference genomes of Drosophila melanogaster, Anopheles gambiae, Apis mellifera, Tribolium castaneum, Bombyx mori, Manduca sexta and Danaus plexippus. We find several highly conserved amino acid residues in the CRAL_TRIO domain-containing genes across insects and a gene expansion resulting in more than twice as many gene family members in lepidopterans than in other surveyed insect species, but no lepidopteran homolog of the PINTA gene in Drosophila. In addition, we examined the expression pattern of CRAL_TRIO domain genes in Manduca sexta heads using RNA-Seq data. Of the 42 gene family members found in the M.sexta reference genome, we found 30 expressed in the head tissue with similar expression profiles between males and females. Our results suggest this gene family underwent a large expansion in lepidopteran, making the lepidopteran CRAL_TRIO domain family distinct from other holometabolous insect lineages.

Multi-omic approaches promise to supply the power to detect genes underlying disease and fitness-related phenotypes. Optimal use of the resulting profusion of data requires detailed investigation of individual candidate genes, a challenging proposition. Here, we combine transcriptomic and genomic data with molecular modelling of candidate enzymes to characterize the evolutionary history and function of the serine protease cocoonase. Heliconius butterflies possess the unique ability to feed on pollen; recent work has identified cocoonase as a candidate gene in pollen digestion. Cocoonase was first described in moths, where it aids in eclosure from the cocoon and is present as a single copy gene. In heliconiine butterflies it is duplicated and highly expressed in the mouthparts of adults. At least six copies of cocoonase are present in Heliconius melpomene and copy number varies across H. melpomene sub-populations. Most cocoonase genes are under purifying selection, however branch-site analyses suggest cocoonase 3 genes may have evolved under episodic diversifying selection. Molecular modelling of cocoonase proteins and examination of their predicted structures revealed that the active site region of each type has a similar structure to trypsin, with the same predicted substrate specificity across types. Variation among heliconiine cocoonases instead lies in the outward-facing residues involved in solvent interaction. Thus, the neofunctionalization of cocoonase duplicates appears to have resulted from the need for these serine proteases to operate in diverse biochemical environments. We suggest that cocoonase may have played a buffering role in feeding during the diversification of Heliconius across the neotropics by enabling these butterflies to digest protein from a range of biochemical milieux.