Abstract
The chemical interactions of insects and host plants are shaping the evolution of chemosensory receptor gene families. However, the correlation between host range and chemoreceptor gene repertoire sizes is still elusive in Papilionidae. Here, we addressed the issue of whether host plant diversities are correlated with the expansions of odorant (ORs) or gustatory (GRs) receptors in six Papilio butterflies. By combining genomics, transcriptomics and bioinformatics approaches, 381 ORs and 328 GRs were annotated in the genomes of a generalist P. glaucus and five specialists, P. xuthus, P. polytes, P. memnon, P. machaon and P. dardanus. Orthologous ORs or GRs in Papilio had highly conserved gene structure. Five Papilio specialists exhibited a similar frequency of intron lengths for ORs or GRs, but which was different from those in the generalist. Phylogenetic analysis revealed 60 orthologous OR groups, 45 of which shared one-to-one relationships. Such a single gene in each butterfly also occurred in 26 GR groups. Intriguingly, bitter GRs had fewer introns than other GRs and clustered into a large clade. Focusing on the two chemoreceptor gene families in P. xuthus, most PxutORs (52/58) were expressed in antennae and 31 genes in reproductive tissues. Eleven out of 28 foretarsus-expressed PxutGRs were female-biased genes, as strong candidates for sensing oviposition stimulants. These results indicate that the host range may not shape the large-scale expansions of ORs and GRs in Papilio butterflies and identify important molecular targets involved in olfaction, oviposition or reproduction in P. xuthus.