Fig. 6

Phylogenetic relationship and orthology of putative Diabrotica virgifera virgifera caspases. Maximum-likelihood (ML) analysis of a partial enzymatic domain alignment between putative D. v. virgifera caspases with seven known caspases from Drosophila melanogaster (Supplementary Fig. 5). The D. v. virgifera transcripts (DIAVI0NNNNN) and protein models (XP 0281NNNNN.1) from RefSeq accession GCF_003013835.1 for the draft genome assembly accession GCF_003013835.1 Dvir_v2 are integrated for homolog estimation. Caspase-encoding transcripts DIAVI027204 and DIAVI022989 up-regulated in both Cry3Bb1 and Cry34/35Ab1 exposed larvae (Table 7) are enclosed in boxes. Clades corresponding to effector and initiator classes are indicated by brackets. The ML of the constructed tree using the LG model of protein sequence evolution (Le and Gascuel 2008) [77] with a gamma shape parameter (G) = 2.7553 and 6.77% of sites defined as evolutionarily invariable (I) was minimized at a log likelihood of − 2224.28, and resulted in a tree with a total branch length of 10.157. Node support obtained using 1000 bootstrap pseudo-replications of the aligned protein sequence data