Comparative genome analysis
To identify ortholog groups among multiple species, including B. mori , D. plexippus , P. xuthus and P. xylostella , OrthoFinder (Emms, & Kelly, 2015) was used. Each gene set corresponded to the genome assembly, which was used for BUSCO analysis (Table S5). Regarding to D. plexippus , P. xuthus and P. xylostella , the proteome data was obtained from Lepbase (http://lepbase.org) (Challis, Kumar, Dasmahapatra, Jiggins, & Blaxter, 2016). The proteome data of B. mori was obtained from SilkBase (http://silkbase.ab.a.u-tokyo.ac.jp/cgi-bin/download.cgi).