Genomics of the relict species Baronia brevicornis sheds light on its demographic history and genome size evolution across swallowtail butterflies.

Relict species, like coelacanth, gingko, tuatara, are the remnants of formerly more ecologically and taxonomically diverse lineages. It raises the questions of why they are currently species-poor, have restrained ecology, and are often vulnerable to extinction. Estimating heterozygosity level and demographic history can guide our understanding of the evolutionary history and conservation status of relict species. However, few studies have focused on relict invertebrates compared to vertebrates. We sequenced the genome of Baronia brevicornis (Lepidoptera: Papilionidae), which is an endangered species, the sister species of all swallowtail butterflies, and is the oldest lineage of all extant butterflies. From a dried specimen, we were able to generate both long-read and short-read data and assembled a genome of 406 Mb for Baronia. We found a fairly high level of heterozygosity (0.58%) compared to other swallowtail butterflies, which contrasts with its endangered and relict status. Taking into account the high ratio of recombination over mutation, demographic analyses indicated a sharp decline of the effective population size initiated in the last million years. Moreover, the Baronia genome was used to study genome size variation in Papilionidae. Genome sizes are mostly explained by transposable elements activities, suggesting that large genomes appear to be a derived feature in swallowtail butterflies as transposable elements activity is recent and involves different transposable elements classes among species. This first Baronia genome provides a resource for assisting conservation in a flagship and relict insect species as well as for understanding swallowtail genome evolution.

Convergence in sympatric swallowtail butterflies reveals ecological interactions as a key driver of worldwide trait diversification.

Ecological interactions can promote phenotypic diversification in sympatric species. While competition can enhance trait divergence, other ecological interactions may promote convergence in sympatric species. Within butterflies, evolutionary convergences in wing color patterns have been reported between distantly related species, especially in females of palatable species, where mimetic color patterns are promoted by predator communities shared with defended species living in sympatry. Wing color patterns are also often involved in species recognition in butterflies, and divergence in this trait has been reported in closely related species living in sympatry as a result of reproductive character displacement. Here, we investigate the effect of sympatry between species on the convergence vs. divergence of their wing color patterns in relation to phylogenetic distance, focusing on the iconic swallowtail butterflies (family Papilionidae). We developed an unsupervised machine learning-based method to estimate phenotypic distances among wing color patterns of 337 species, enabling us to finely quantify morphological diversity at the global scale among species and allowing us to compute pairwise phenotypic distances between sympatric and allopatric species pairs. We found phenotypic convergence in sympatry, stronger among distantly related species, while divergence was weaker and restricted to closely related males. The convergence was stronger among females than males, suggesting that differential selective pressures acting on the two sexes drove sexual dimorphism. Our results highlight the significant effect of ecological interactions driven by predation pressures on trait diversification in Papilionidae and provide evidence for the interaction between phylogenetic proximity and ecological interactions in sympatry, acting on macroevolutionary patterns of phenotypic diversification.

A draft reference genome assembly of the Pipevine Swallowtail butterfly, Battus philenor hirsuta.

The California Pipevine Swallowtail Butterfly, Battus philenor hirsuta, and its host plant, the California Pipevine or Dutchman's Pipe, Aristolochia californica Torr., are an important California endemic species pair. While this species pair is an ideal system to study co-evolution, genomic resources for both are lacking. Here, we report a new, chromosome-level assembly of B. philenor hirsuta as part of the California Conservation Genomics Project (CCGP). Following the sequencing and assembly strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin proximity sequencing technology to produce a de novo assembled genome. Our genome assembly, the first for any species in the genus, contains 109 scaffolds spanning 443 mega base (Mb) pairs, with a contig N50 of 14.6 Mb, a scaffold N50 of 15.2 Mb, and BUSCO complete score of 98.9%. In combination with the forthcoming A. californica reference genome, the B. philenor hirsuta genome will be a powerful tool for documenting landscape genomic diversity and plant-insect co-evolution in a rapidly changing California landscape.

Lost and Found: Taxonomic Notes on Some Butterflies (Lepidoptera: Papilionoidea) from J.F. Zikán found on F. Dissmann's Collection.

J. F. Zikán described several new taxa of butterflies from the Neotropical region. The majority of the butterfly types of J. F. Zikán was deposited at the Instituto Oswaldo Cruz, Rio de Janeiro, Brazil (IOC). However, some type series were found to be incomplete. The present paper discusses and illustrate some type specimens of Papilionidae and Nymphalidae from J. F. Zikán found in Dissmann's collection, that now is deposited at Museu de Diversidade Biológica from University of Campinas - UNICAMP, Campinas, São Paulo, Brazil (ZUEC).

Complete mitochondrial genomes of Papilio nephelus chaon and Papilio epycides (Lepidoptera: Papilionidae: Papilioninae) and phylogenetic analysis.

The complete mitochondrial genome (mitogenome) sequences of Papilio nephelus chaon and Papilio epycides were sequenced by Illumina and analyzed in this study. They are 15,287 bp and 15,012 bp in size, respectively, and contains 13 protein-coding genes (PCGs), 22 tRNA genes (tRNAs), 2 rRNA genes (rRNAs), and 1 AT-rich control region (CR). The phylogenetic relationships of 56 species in the Papilionidae were inferred based on concatenated nucleotide sequences by using Maximum Likelihood with the selected best-fit model GTR + F+R6. The phylogenetic analysis showed that P. nephelus chaon and P. epycides were located in the genus Papilio. This study provides a basis for further study on mitogenome and phylogenetics of the Papilionidae.

A phylogenetic study to assess the link between biome specialization and diversification in swallowtail butterflies.

The resource-use hypothesis, proposed by E.S. Vrba, states that habitat fragmentation caused by climatic oscillations would affect particularly biome specialists (species inhabiting only one biome), which might show higher speciation and extinction rates than biome generalists. If true, lineages would accumulate biome-specialist species. This effect would be particularly exacerbated for biomes located at the periphery of the global climatic conditions, namely, biomes that have high/low precipitation and high/low temperature such as rainforest (warm-humid), desert (warm-dry), steppe (cold-dry) and tundra (cold-humid). Here, we test these hypotheses in swallowtail butterflies, a clade with more than 570 species, covering all the continents but Antarctica, and all climatic conditions. Swallowtail butterflies are among the most studied insects, and they are a model group for evolutionary biology and ecology studies. Continental macroecological rules are normally tested using vertebrates, this means that there are fewer examples exploring terrestrial invertebrate patterns at global scale. Here, we compiled a large Geographic Information System database on swallowtail butterflies' distribution maps and used the most complete time-calibrated phylogeny to quantify diversification rates (DRs). In this paper, we aim to answer the following questions: (1) Are there more biome-specialist swallowtail butterflies than biome generalists? (2) Is DR related to biome specialization? (3) If so, do swallowtail butterflies inhabiting extreme biomes show higher DRs? (4) What is the effect of species distribution area? Our results showed that swallowtail family presents a great number of biome specialists which showed substantially higher DRs compared to generalists. We also found that biome specialists are unevenly distributed across biomes. Overall, our results are consistent with the resource-use hypothesis, species climatic niche and biome fragmentation as key factors promoting isolation.

Evidence of attack deflection suggests adaptive evolution of wing tails in butterflies.

Predation is a powerful selective force shaping many behavioural and morphological traits in prey species. The deflection of predator attacks from vital parts of the prey usually involves the coordinated evolution of prey body shape and colour. Here, we test the deflection effect of hindwing (HW) tails in the swallowtail butterfly Iphiclides podalirius. In this species, HWs display long tails associated with a conspicuous colour pattern. By surveying the wings within a wild population of I. podalirius, we observed that wing damage was much more frequent on the tails. We then used a standardized behavioural assay employing dummy butterflies with real I. podalirius wings to study the location of attacks by great tits Parus major. Wing tails and conspicuous coloration of the HWs were struck more often than the rest of the body by birds. Finally, we characterized the mechanical properties of fresh wings and found that the tail vein was more fragile than the others, suggesting facilitated escape ability of butterflies attacked at this location. Our results clearly support the deflective effect of HW tails and suggest that predation is an important selective driver of the evolution of wing tails and colour pattern in butterflies.

The complete mitochondrial genome data of the Common Rose butterfly, Pachliopta aristolochiae (Lepidoptera, Papilionoidea, Papilionidae) from Malaysia.

Here, we present the complete mitochondrial genome of Pachliopta aristolochiae, a Common Rose butterfly from Malaysia. The sequence was generated using Illumina NovaSeq 6000 sequencing platform. The mitogenome is 15,235bp long, consisting of 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and two D-loop regions. The total base composition was (81.6%), with A (39.3%), T (42.3%), C (11.0%) and G (7.3%). The gene order of the three tRNAs was trnM-trnI-trnQ, which differs from the ancestral insect gene order trnI-trnQ-trnM. Phylogenetic tree analysis revealed that the sequenced Pachliopta aristolochiae in this data is closely related to Losaria neptunus (NC 037868), with highly supported ML and BI analysis. The data presented in this work can provide useful resources for other researchers to study deeper into the phylogenetic relationships of Lepidoptera and the diversification of the Pachliopta species. Also, as one of the bioindicator species, this data can be used to assess environmental changes in the terrestrial and aquatic ecosystem via enviromental DNA approahes. The mitogenome of Pachliopta aristolochiae is available in GenBank under the accession number MZ781228.

Chemical Identity of Cuticular Lipid Components in the Mimetic Swallowtail Butterfly Papilio polytes.

The swallowtail Papilio polytes shows Batesian and female-limited polymorphic mimicry. In Japan, P. polytes females have two different forms: the cyrus form is non-mimetic and resembles males, whereas the polytes form mimics Pachliopta aristolochiae and Byasa (Atrophaneura) alcinous as unpalatable models. During mating, P. polytes males use cuticular lipids to distinguish non-mimetic females from conspecific males and sympatric sister species. In this study, we investigated whether compositional differences in cuticular lipids exist between mimetic and non-mimetic females of P. polytes and between mimetic females and their model species. The mimetic and non-mimetic females had nearly identical cuticular lipid profiles, which differed from those of males. The two model species exhibited sexually dimorphic and species-specific cuticular lipid compositions, which were distinctly different from those of mimetic P. polytes females. These results strongly suggest that P. polytes females maintain the identity of cuticular lipid profiles regardless of the mimicry type, and this feature helps males recognize mimetic females as the correct mating partners.

The complete mitogenome of the Paranticopsis xenocles (Lepidoptera: Papilionidae: Papilioninae) and phylogenetic implications.

Paranticopsis xenocles Doubleday belongs to the Paranticopsis of Papilionidae and is mainly distributed in China mainland. Herein, we report the complete mitogenome of P. xenocles reconstructing from Illumina sequence data. The mitogenome is 15,187 bp in length and contains 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. The phylogenetic analysis indicated that P. xenocles were clustered within Paranticopsis. This study would provide useful genetic information for future studies on taxonomy, phylogeny, and evolution of Papilionidae species.

α-Linolenic acid in Papilio machaon larvae regurgitant induces a defensive response in Apiaceae.

Papilio machaon hippocrates C. Felder et R. Felder, 1864 (Papilionidae) larvae are pests of plants of the family Apiaceae. It is unclear whether Apiaceae plants show induced defensive responses against P. machaon hippocrates larvae, and if so, how these responses are induced. Comparison of the fatty acid (FA) composition of the leaves of host plants and the regurgitant of P. machaon hippocrates larvae by LC-MS revealed a great increase in α-linolenic acid (α-LA) in the regurgitant compared with the FAs contained in the leaves. However, specific FA amino acid conjugates, known as elicitor compounds, such as volicitin, were not detected. Sterile host plants (Saposhnikovia divaricata (Turcz.) Schischk., Apiaceae) were treated with α-LA to mimic the damage made by P. machaon hippocrates larvae. After α-LA treatment to leaves, induced defensive reactions, i.e., release of volatile compounds such as α- and ß-pinene and camphene (possible induced indirect defense) and the accumulation of specialized metabolites such as (R)-falcarinol and bergapten (possible induced direct defense) were observed. These findings highlight the role of α-LA in the interaction between P. machaon hippocrates larvae and Apiaceae host plants.

Notes on some Tanzanian butterfly specimens in the Suffert Collection: a case of patria/ falsa/.

A number of mislabelled butterfly specimens in the Suffert collection purporting to be from Lindi in southeast Tanzania belong to taxa restricted to the Usambara Mountains in northeast Tanzania thus suggesting patria falsa. The history of Suffert's collection and the possible reasons behind this mislabelling event are discussed in more detail. One of the mislabelled specimens, the holotype of Papilio illyris hamatus Joicey Talbot, 1918 described from Tanzania, is the senior name for Graphium illyris girardeaui Guilbot Plantrou, 1978, the latter name formally synonymised herein with the former. A brief biography of the little-known Ernst Suffert is presented.

The complete mitochondrial genome of Parnassius mercurius Grum-Grshimailo (Lepidoptera: Papilionidae: Parnassiinae).

The mitogenome of Parnassius mercurius Grum-Grshimailo is determined to be 15,372 bp in length, including 37 typical insect mitochondrial genes and an AT-rich region. All PCGs start with ATN, except for COI with CGA; 12 of 13 PCGs harbour the common stop codon TAA or TAG, whereas COII end with a single T. The lrRNA and srRNA genes are 1344 bp and 775 bp in length, respectively. The AT-rich region contains several features characteristic of the lepidopterans. Phylogenetic analysis shows that P. mercurius is the closest relative of P. epaphus and P. nomion lineage, rather than the P. jacquemontii.

Phylogenomic test of mitochondrial clues to archaic ancestors in a group of hybridizing swallowtail butterflies.

Genomics has revolutionized our understanding of hybridization and introgression, but most of the early evidence for these processes came from studies of mitochondrial introgression. To expand these evolutionary insights from mitochondrial patterns, we evaluate phylogenetic discordance across the nuclear genomes of a hybridizing system, the Papilio machaon group of swallowtail butterflies. This species group contains three hybrid lineages (P. brevicauda, P. joanae, and P. m. kahli) that are geographically disjunct across North America and have complete fixation of a mitochondrial lineage that is otherwise primarily found in P. m. hudsonianus, a boreal subspecies of the Holarctic P. machaon. Genome-wide nuclear markers place the three hybrid lineages as a monophyletic group that is sister to P. polyxenes/P. zelicaon rather than P. machaon, although ancient hybridization between a subspecies of P. machaon and the ancestor of these three lineages is also shown by their greater nuclear affinity to P. m. hudsonianus than to other subspecies of P. machaon. Individuals from contemporary hybrid swarms in Alberta, where mitochondrial DNA fixation has not occurred, were more intermediate between their respective parent species, demonstrating diversity in mito-nuclear discordance following hybrid interactions. Our new phylogenetic findings for the P. machaon species group also include: subspecific paraphyly within P. machaon itself across its Holarctic distribution; paraphyly of P. zelicaon relative to P. polyxenes; and more divergent placement of a Mediterranean species, P. hospiton. These results provide the first comprehensive genomic evaluation of relationships within this species group and provide insight into the evolutionary dynamics of hybridization and mitochondrial introgression.

Swallowtail butterflies (Lepidoptera: Papilionidae) species diversity and distribution in Africa: The Papilionidae collection at the National Museums of Kenya, Nairobi, Kenya.

BACKGROUND: Species data from the Museum collections have been shown to be of great value as a tool for prioritising conservation actions in Africa (Fjeldsa and Tushabe 2005). The National Museums of Kenya (NMK) have an entomology collection, housed in 4,000 drawers in cabinets that contain over 1.5 million specimens, including the largest butterfly collection in Africa (Arnett et al. 1997). Lampe and Striebing (2005) demonstrated how to digitise large insect collections in order to make their associated label data into databases that can be used for functions, such as creating distribution maps. The NMK's swallowtail butterflies' collection had not been digitised and thus there was a need to capture the label data to create a database that can be used for mapping the distribution of the species in Kenya and elsewhere. These data have addressed one of the most significant challenges to insect conservation i.e. the lack of baseline information concerning species diversity and distribution (Summerville and Crist 2003). These data have provided key historic papilionid species diversity and distribution data that can be used to monitor their populations, as butterflies are declining due to changes in land use, intensive agriculture and pestcide use, diseases and pest and climate change (Potts et al. 2016; Bongaarts 2019). The publication of the occurrence data records in GBIF has been undertaken, thus making the data available to a wider audience and promoting availability for use. NEW INFORMATION: The swallowtail butterflies collection at the National Museums of Kenya was digitised from 2017-2019 and this paper presents details of the Papilionid collection at the Zoology Department, NMK, Nairobi, Kenya.The collection holds 7,345 voucher specimens, consisting of three genera and 133 species. The collection covers the period between 1850 to 2019.The distribution of the swallowtail butterflies, housed at the NMK, covers East Africa with 88%, Central Africa (6%), Western Africa (4%) and Southern Africa (2%).

Flower orientation in Gloriosa superba (Colchicaceae) promotes cross-pollination via butterfly wings.

BACKGROUND AND AIMS: Complex modifications of angiosperm flowers often function for precise pollen placement on pollinators and to promote cross-pollination. We explore the functional significance of the unusually elaborate morphology of Gloriosa superba flowers, which are divided into one hermaphrodite meranthium and five male meranthia (functional pollination units of a single flower). METHODS: We used controlled pollination experiments, floral measurements, pollen load analyses and visitor observations in four populations of G. superba in South Africa to determine the breeding system, mechanism of pollination and role of flower in the promotion of cross-pollination. KEY RESULTS: We established that G. superba is self-compatible, but reliant on pollinators for seed production. Butterflies, in particular the pierid Eronia cleodora, were the primary pollinators (>90 % of visitors). Butterflies brush against the anthers and stigma during nectar feeding and pollen is carried on their ventral wing surfaces. Butterfly scales were positively correlated with the number of pollen grains on stigmas. We demonstrate that the styles were orientated towards clearings in the vegetation and we confirm that the highest proportion of initial visits was to hermaphrodite meranthia pointing towards clearings. CONCLUSIONS: The flower morphology of G. superba results in effective pollen transfer on the wings of butterfly visitors. The style-bearing hermaphrodite meranthium of the flowers orientates towards open spaces in the vegetation, thus increasing the probability that butterflies land first on the hermaphrodite meranthium. This novel aspect of flower orientation is interpreted as a mechanism that promotes cross-pollination.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution.

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies-swallowtails of the family Papilionidae-that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6-142.7 Ma) and Papilionidae around 71.4 Ma (49.8-103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

An Oviposition Stimulant for a Magnoliaceae-Feeding Swallowtail Butterfly, Graphium doson, from its Primary Host Plant, Michelia compressa.

Chemical examination of plant constituents responsible for oviposition by a Magnoliaceae-feeding butterfly, Graphium doson, was conducted using its major host plant, Michelia compressa. A methanol extract prepared from young leaves of the plant elicited a strong oviposition response from females. The methanolic extract was then separated by solvent partition into three fractions: CHCl3, i-BuOH, and aqueous fractions. Active substance(s) resided in both i-BuOH- and water-soluble fractions. Bioassay-guided further fractionation of the water-soluble substances by means of various chromatographic techniques led to the isolation of an oviposition stimulant. The stimulant was identified as D-(+)-pinitol on the basis of 13C NMR spectra and physicochemical properties. D-(+)-Pinitol singly exhibited a moderate oviposition-stimulatory activity at a dose of 150 µg/cm2. This compound was present also in another host plant, Magnolia grandiflora, in a sufficient amount to induce oviposition behavior of G. doson females. Certain cyclitols including D-(+)-pinitol have been reported to be involved in stimulation of oviposition by some Aristolochiaceae- and Rutaceae-feeding papilionid butterflies. A possible pathway of phytochemical-mediated host shifts in the Papilionidae, in which certain cyclitols could enact important mediators, is discussed in relation to the evolution of cyclitol biosynthesis in plants.

Cyanogenesis in Arthropods: From Chemical Warfare to Nuptial Gifts.

Chemical defences are key components in insect⁻plant interactions, as insects continuously learn to overcome plant defence systems by, e.g., detoxification, excretion or sequestration. Cyanogenic glucosides are natural products widespread in the plant kingdom, and also known to be present in arthropods. They are stabilised by a glucoside linkage, which is hydrolysed by the action of β-glucosidase enzymes, resulting in the release of toxic hydrogen cyanide and deterrent aldehydes or ketones. Such a binary system of components that are chemically inert when spatially separated provides an immediate defence against predators that cause tissue damage. Further roles in nitrogen metabolism and inter- and intraspecific communication has also been suggested for cyanogenic glucosides. In arthropods, cyanogenic glucosides are found in millipedes, centipedes, mites, beetles and bugs, and particularly within butterflies and moths. Cyanogenic glucosides may be even more widespread since many arthropod taxa have not yet been analysed for the presence of this class of natural products. In many instances, arthropods sequester cyanogenic glucosides or their precursors from food plants, thereby avoiding the demand for de novo biosynthesis and minimising the energy spent for defence. Nevertheless, several species of butterflies, moths and millipedes have been shown to biosynthesise cyanogenic glucosides de novo, and even more species have been hypothesised to do so. As for higher plant species, the specific steps in the pathway is catalysed by three enzymes, two cytochromes P450, a glycosyl transferase, and a general P450 oxidoreductase providing electrons to the P450s. The pathway for biosynthesis of cyanogenic glucosides in arthropods has most likely been assembled by recruitment of enzymes, which could most easily be adapted to acquire the required catalytic properties for manufacturing these compounds. The scattered phylogenetic distribution of cyanogenic glucosides in arthropods indicates that the ability to biosynthesise this class of natural products has evolved independently several times. This is corroborated by the characterised enzymes from the pathway in moths and millipedes. Since the biosynthetic pathway is hypothesised to have evolved convergently in plants as well, this would suggest that there is only one universal series of unique intermediates by which amino acids are efficiently converted into CNglcs in different Kingdoms of Life. For arthropods to handle ingestion of cyanogenic glucosides, an effective detoxification system is required. In butterflies and moths, hydrogen cyanide released from hydrolysis of cyanogenic glucosides is mainly detoxified by β-cyanoalanine synthase, while other arthropods use the enzyme rhodanese. The storage of cyanogenic glucosides and spatially separated hydrolytic enzymes (β-glucosidases and α-hydroxynitrile lyases) are important for an effective hydrogen cyanide release for defensive purposes. Accordingly, such hydrolytic enzymes are also present in many cyanogenic arthropods, and spatial separation has been shown in a few species. Although much knowledge regarding presence, biosynthesis, hydrolysis and detoxification of cyanogenic glucosides in arthropods has emerged in recent years, many exciting unanswered questions remain regarding the distribution, roles apart from defence, and convergent evolution of the metabolic pathways involved.

Taste sensitivity and divergence in host plant acceptance between adult females and larvae of Papilio hospiton.

On the island of Sardinia the lepidopteran Papilio hospiton uses Ferula communis as exclusive host plant. However, on the small island of Tavolara, adult females lay eggs on Seseli tortuosum, a plant confined to the island. When raised in captivity on Seseli only few larvae grew beyond the first-second instar. Host specificity of lepidopterans is determined by female oviposition preferences, but also by larval food acceptance, and adult and larval taste sensitivity may be related to host selection in both cases. Aim of this work was: (i) to study the taste sensitivity of larvae and ovipositing females to saps of Ferula and Seseli; (ii) to cross-compare the spike activity of gustatory receptor neurons (GRNs) to both taste stimuli; (iii) to evaluate the discriminating capability between the two saps and determine which neural code/s is/are used. The results show that: (i) the spike responses of the tarsal GRNs of adult females to both plant saps are not different and therefore they cannot discriminate the two plants; (ii) larval L-lat GRN shows a higher activity in response to Seseli than Ferula, while the opposite occurs for the phagostimulant neurons, and larvae may discriminate between the two saps by means of multiple neural codes; (iii) the number of eggs laid on the two plants is the same, but the larval growth performance is better on Ferula than Seseli. Taste sensitivity differences may explain the absence of a positive relationship between oviposition preferences by adult females and plant acceptance and growth performance by larvae.