Testing the role of ecological selection on colour pattern variation in the butterfly Parnassius clodius

No evidence for the melanin desiccation hypothesis in a larval Lepidopteran

Water regulation is an important physiological challenge for insects due to their small body sizes and large surface area to volume ratios. Adaptations for decreasing cuticular water loss, the largest avenue of loss, are especially important. The melanin desiccation hypothesis states that melanin molecules in the cuticle may help prevent water loss, thus offering protection from desiccation. This hypothesis has much empirical support in Drosophila species, but remains mostly untested in other taxa, including Lepidoptera. Because melanin has many other important functions in insects, its potential role in desiccation prevention is not always clear. In this study we investigated the role of melanin in desiccation prevention in the white-lined Sphinx moth, Hyles lineata (Lepidoptera, Sphingidae), which shows high plasticity in the degree of melanin pigmentation during the late larval instars. We took advantage of this plasticity and used density treatments to induce a wide range of cuticular melanization; solitary conditions induced low melanin pigmentation while crowded conditions induced high melanin pigmentation. We tested whether more melanic larvae from the crowded treatment were better protected from desiccation in three relevant responses: i) total water loss over a desiccation period, ii) change in hemolymph osmolality over a desiccation period, and iii) evaporation rate of water through the cuticle. We did not find support for the melanin desiccation hypothesis in this species. Although treatment influenced total water loss, this effect did not occur via degree of melanization. Interestingly, this implies that crowding, which was used to induce high melanin phenotypes, may have other physiological effects that influence water regulation. There were no differences between treatments in cuticular evaporative water loss or change in hemolymph osmolality. However, we conclude that osmolality may not sufficiently reflect water loss in this case. This study emphasizes the context dependency of melanin's role in desiccation prevention and the importance of considering how it may vary across taxa. In lepidopteran larvae that are constantly feeding phytophagous insects with soft cuticles, melanin may not be necessary for preventing cuticular water loss.

Taxonomic Structure and Wing Pattern Evolution in the Parnassius mnemosyne Species Complex (Lepidoptera, Papilionidae)

In our study, using the analysis of DNA barcodes and morphology (wing color, male genitalia, and female sphragis shape), we show that the group of species close to P. mnemosyne comprises the western and eastern phylogenetic lineages. The eastern lineage includes P. stubbendorfii, P. glacialis, and P. hoenei. The western lineage includes three morphologically similar species: P. mnemosyne (Western Eurasia), P. turatii (southwestern Europe), and P. nubilosusstat. nov. (Turkmenistan and NE Iran), as well as the morphologically differentiated P. ariadne (Altai). The latter species differs from the rest of the group in the presence of red spots on the wings. Parnassius mnemosyne s.s. is represented by four differentiated mitochondrial clusters that show clear association with specific geographic regions. We propose to interpret them as subspecies: P. mnemosyne mnemosyne (Central and Eastern Europe, N Caucasus, N Turkey), P. mnemosyne adolphi (the Middle East), P. mnemosyne falsa (Tian Shan), and P. mnemosyne gigantea (Gissar-Alai in Central Asia). We demonstrate that in P. ariadne, the red spots on the wing evolved as a reversion to the ancestral wing pattern. This reversion is observed in Altai, where the distribution areas of the western lineage, represented by P. ariadne, and the eastern lineage, represented by P. stubbendorfii, overlap. These two species hybridize in Altai, and we hypothesize that the color change in P. ariadne is the result of reinforcement of prezygotic isolation in the contact zone. The lectotype of Parnassius mnemosyne var. nubilosus Christoph, 1873, is designated.

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large‐scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD‐sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species‐specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

The impact of life stage and pigment source on the evolution of novel warning signal traits

Our understanding of how novel warning color traits evolve in natural populations is largely based on studies of reproductive stages and organisms with endogenously produced pigmentation. In these systems, genetic drift is often required for novel alleles to overcome strong purifying selection stemming from frequency‐dependent predation and positive assortative mating. Here, we integrate data from field surveys, predation experiments, population genomics, and phenotypic correlations to explain the origin and maintenance of geographic variation in a diet‐based larval pigmentation trait in the redheaded pine sawfly (Neodiprion lecontei), a pine‐feeding hymenopteran. Although our experiments confirm that N. lecontei larvae are indeed aposematic—and therefore likely to experience frequency‐dependent predation—our genomic data do not support a historical demographic scenario that would have facilitated the spread of an initially deleterious allele via drift. Additionally, significantly elevated differentiation at a known color locus suggests that geographic variation in larval color is currently maintained by selection. Together, these data suggest that the novel white morph likely spread via selection. However, white body color does not enhance aposematic displays, nor is it correlated with enhanced chemical defense or immune function. Instead, the derived white‐bodied morph is disproportionately abundant on a pine species with a reduced carotenoid content relative to other pine hosts, suggesting that bottom‐up selection via host plants may have driven divergence among populations. Overall, our results suggest that life stage and pigment source can have a substantial impact on the evolution of novel warning signals, highlighting the need to investigate diverse aposematic taxa to develop a comprehensive understanding of color variation in nature.

Local adaptation and migratory habits balance spatial-genetic structure between continental and insular chestnut tiger butterflies in East Asia

Geographic and climatic differences between islands and continents may affect the evolution of their biota, and promote divergent selection in species distributed in both landscapes. To assess spatial-genetic structure, we genotyped 18 expressed sequence tag-simple sequence repeat (EST-SSR) loci and sequenced two mtDNA markers (ND5 and COI) and one nuclear marker (EF1α) in two subspecies of the butterfly Parantica sita. Compared with nuclear markers, mtDNA had a stronger signal of population structure. Approximate Bayesian computation (ABC) suggested that a continuous-gene-flow model best described the data. According to this model, the two subspecies diverged approximately 23.1 kya, with ten times more introgression from the continental (ssp. sita) to the insular subspecies (ssp. niphonica) than vice versa. Ecological niche modeling was performed to predict the paleo- and current potential distributions and elucidate the geohistorical process, which revealed a northeastern, insular origin. Winter precipitation and annual temperature range were the main determinants of the subspecies distributions. Maximum-likelihood population-effects models showed that the population differentiation of the insular and continental subspecies was primarily affected by environmental resistance and local climate. Sex-biased migration capacity and long-term precipitation-driven divergence between the continental and insular lineages shaped the current genetic structure of P. sita. Evidence from the nuclear markers confirmed inter-subspecific gene flow despite adaptive divergence between the subspecies. These results imply that the continental subspecies is still capable of returning to the island and introgressing with the insular subspecies.

Drainage basins serve as multiple glacial refugia for alpine habitats in the Sierra Nevada Mountains, California

The evolutionary histories of alpine species are often directly associated with responses to glaciation. Deep divergence among populations and complex patterns of genetic variation have been inferred as consequences of persistence within glacier boundaries (i.e., on nunataks), while shallow divergence and limited genetic variation are assumed to result from expansion from large refugia at the edge of ice shields (i.e., massifs de refuge). However, for some species, dependence on specific microhabitats could profoundly influence their spatial and demographic response to glaciation, and such a simple dichotomy may obscure the localization of actual refugia. In this study, we use the Nebria ingens complex (Coleoptera: Carabidae), a water-affiliated ground beetle lineage, to test how drainage basins are linked to their observed population structure. By analysing mitochondrial COI gene sequences and genome-wide single nucleotide polymorphisms, we find that the major drainage systems of the Sierra Nevada Mountains in California best explain the population structure of the N. ingens complex. In addition, we find that an intermediate morphotype within the N. ingens complex is the product of historical hybridization of N. riversi and N. ingens in the San Joaquin basin during glaciation. This study highlights the importance of considering ecological preferences in how species respond to climate fluctuations and provides an explanation for discordances that are often observed in comparative phylogeographical studies.