Habitat choice precedes host plant choice - niche separation in a species pair of a generalist and a specialist butterfly

Evolution of realized niche breadth diversity driven by community dynamics

Why many herbivorous insects are host plant specialists, with non-negligible exceptions, is a conundrum of evolutionary biology, especially because the host plants are not necessarily optimal larval diets. Here, I present a novel model of host plant preference evolution of two insect species. Because habitat preference evolution is contingent upon demographic dynamics, I integrate the evolutionary framework with the modern coexistence theory. The results show that the two insect species can evolve into a habitat specialist and generalist, when they experience both negative and positive frequency-dependent community dynamics. This happens because the joint action of positive and negative frequency dependence creates multiple (up to nine) eco-evolutionary equilibria. Furthermore, initial condition dependence due to positive frequency dependence allows specialization to poor habitats. Thus, evolved habitat preferences do not necessarily correlate with the performances. The model provides explanations for counterintuitive empirical patterns and mechanistic interpretations for phenomenological models of niche breadth evolution.

Phylogeography of Two Enigmatic Sulphur Butterflies, Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897 (Lepidoptera, Pieridae), with Relations to Wolbachia Infection

The genus Colias Fabricius, 1807 includes numerous taxa and forms with uncertain status and taxonomic position. Among such taxa are Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897, interpreted in the literature either as conspecific forms, as subspecies of different but morphologically somewhat similar Colias species or as distinct species-level taxa. Based on mitochondrial and nuclear DNA markers, we reconstructed a phylogeographic pattern of the taxa in question. We recover and include in our analysis DNA barcodes of the century-old type specimens, the lectotype of C. tamerlana deposited in the Natural History Museum (Museum für Naturkunde), Berlin, Germany (ZMHU) and the paralectotype of C. tamerlana and the lectotype of C. mongola deposited in the Zoological Institute, Russian Academy of Sciences, St. Petersburg, Russia (ZISP). Our analysis grouped all specimens within four (HP_I-HP_IV) deeply divergent but geographically poorly structured clades which did not support nonconspecifity of C. mongola-C. tamerlana. We also show that all studied females of the widely distributed haplogroup HP_II were infected with a single Wolbachia strain belonging to the supergroup B, while the males of this haplogroup, as well as all other investigated specimens of both sexes, were not infected. Our data highlight the relevance of large-scale sampling dataset analysis and the need for testing for Wolbachia infection to avoid erroneous phylogenetic reconstructions and species misidentification.

High-density linkage maps and chromosome level genome assemblies unveil direction and frequency of extensive structural rearrangements in wood white butterflies (Leptidea spp.)

Karyotypes are generally conserved between closely related species and large chromosome rearrangements typically have negative fitness consequences in heterozygotes, potentially driving speciation. In the order Lepidoptera, most investigated species have the ancestral karyotype and gene synteny is often conserved across deep divergence, although examples of extensive genome reshuffling have recently been demonstrated. The genus Leptidea has an unusual level of chromosome variation and rearranged sex chromosomes, but the extent of restructuring across the rest of the genome is so far unknown. To explore the genomes of the wood white (Leptidea) species complex, we generated eight genome assemblies using a combination of 10X linked reads and HiC data, and improved them using linkage maps for two populations of the common wood white (L. sinapis) with distinct karyotypes. Synteny analysis revealed an extensive amount of rearrangements, both compared to the ancestral karyotype and between the Leptidea species, where only one of the three Z chromosomes was conserved across all comparisons. Most restructuring was explained by fissions and fusions, while translocations appear relatively rare. We further detected several examples of segregating rearrangement polymorphisms supporting a highly dynamic genome evolution in this clade. Fusion breakpoints were enriched for LINEs and LTR elements, which suggests that ectopic recombination might be an important driver in the formation of new chromosomes. Our results show that chromosome count alone may conceal the extent of genome restructuring and we propose that the amount of genome evolution in Lepidoptera might still be underestimated due to lack of taxonomic sampling.

Preference Provides a Plethora of Problems (Don't Panic)

This review was solicited as an autobiography. The "problems" in my title have two meanings. First, they were professional difficulties caused by my decision to study oviposition preferences of butterflies that were not susceptible to traditional preference-testing designs. Until I provided video, my claim that the butterflies duplicate natural post-alighting host-assessment behavior when placed on hosts by hand was not credible, and the preference-testing technique that I had developed elicited skepticism, anger, and derision. The second meaning of "problems" is scientific. Insect preference comes with complex dimensionality that interacts with host acceptability. Part Two of this review describes how my group's work in this area has revealed unexpected axes of variation in plant-insect interactions-axes capable of frustrating attempts to derive unequivocal conclusions from apparently sensible experimental designs. The possibility that these complexities are lurking should be kept in mind as preference and performance experiments are devised.

Multidimensional plasticity in the Glanville fritillary butterfly: larval performance is temperature, host and family specific

Variation in environmental conditions during development can lead to changes in life-history traits with long-lasting effects. Here, we study how variation in temperature and host plant (i.e. the consequences of potential maternal oviposition choices) affects a suite of life-history traits in pre-diapause larvae of the Glanville fritillary butterfly. We focus on offspring survival, larval growth rates and relative fat reserves, and pay specific attention to intraspecific variation in the responses (G × E × E). Globally, thermal performance and survival curves varied between diets of two host plants, suggesting that host modifies the temperature impact, or vice versa. Additionally, we show that the relative fat content has a host-dependent, discontinuous response to developmental temperature. This implies that a potential switch in resource allocation, from more investment in growth at lower temperatures to storage at higher temperatures, is dependent on the larval diet. Interestingly, a large proportion of the variance in larval performance is explained by differences among families, or interactions with this variable. Finally, we demonstrate that these family-specific responses to the host plant remain largely consistent across thermal environments. Together, the results of our study underscore the importance of paying attention to intraspecific trait variation in the field of evolutionary ecology.

Dissecting the Effects of Selection and Mutation on Genetic Diversity in Three Wood White (Leptidea) Butterfly Species

The relative role of natural selection and genetic drift in evolution is a major topic of debate in evolutionary biology. Most knowledge spring from a small group of organisms and originate from before it was possible to generate genome-wide data on genetic variation. Hence, it is necessary to extend to a larger number of taxonomic groups, descriptive and hypothesis-based research aiming at understanding the proximate and ultimate mechanisms underlying both levels of genetic polymorphism and the efficiency of natural selection. In this study, we used data from 60 whole-genome resequenced individuals of three cryptic butterfly species (Leptidea sp.), together with novel gene annotation information and population recombination data. We characterized the overall prevalence of natural selection and investigated the effects of mutation and linked selection on regional variation in nucleotide diversity. Our analyses showed that genome-wide diversity and rate of adaptive substitutions were comparatively low, whereas nonsynonymous to synonymous polymorphism and substitution levels were comparatively high in Leptidea, suggesting small long-term effective population sizes. Still, negative selection on linked sites (background selection) has resulted in reduced nucleotide diversity in regions with relatively high gene density and low recombination rate. We also found a significant effect of mutation rate variation on levels of polymorphism. Finally, there were considerable population differences in levels of genetic diversity and pervasiveness of selection against slightly deleterious alleles, in line with expectations from differences in estimated effective population sizes.

Butterflies embrace maladaptation and raise fitness in colonizing novel host

We illustrate an evolutionary host shift driven by increased fitness on a novel host, despite maladaptation to it in six separate host-adaptive traits. Here, local adaptation is defined as possession of traits that provide advantage in specific environmental contexts; thus individuals can have higher fitness in benign environments to which they are maladapted than in demanding environments to which they are well adapted. A population of the butterfly Euphydryas editha adapted to a long-lived, chemically well-defended host, Pedicularis, had traditionally been under natural selection to avoid the ephemeral, less-defended Collinsia. The lifespan of Collinsia was so short that it senesced before larvae entered diapause. After logging killed Pedicularis in clear-cut patches and controlled burning simultaneously extended Collinsia lifespan, insect fitness on Collinsia in clearings suddenly became higher than on Pedicularis in adjacent unlogged patches. Collinsia was rapidly colonized and preference for it evolved, but insects feeding on it retained adaptations to Pedicularis in alighting bias, two aspects of postalighting oviposition preference, dispersal bias, geotaxis, and clutch size, all acting as maladaptations to Collinsia. Nonetheless, populations boomed on Collinsia in clearings, creating sources that fed pseudosinks in unlogged patches where Pedicularis was still used. After c. 20 years, butterfly populations in clearings disappeared and the metapopulation reverted to Pedicularis-feeding. Here we show, via experimental manipulation of oviposition by local Pedicularis-adapted and imported Collinsia-adapted butterflies, that the highest survival at that time would have been from eggs laid in clearings by butterflies adapted to Collinsia. Second highest were locals on Pedicularis. In third place would have been locals on Collinsia in clearings, because local females maladaptively preferred senescent plants. Collinsia had been colonized despite maladaptation and, after successional changes, abandoned because of it. However, the abandoned Collinsia could still have provided the highest fitness, given appropriate adaptation. The butterflies had tumbled down an adaptive peak.

Female fecundity variation affects reproducibility of experiments on host plant preference and acceptance in a phytophagous insect

Reproducibility is a scientific cornerstone. Many recent studies, however, describe a reproducibility crisis and call for assessments of reproducibility across scientific domains. Here, we explore the reproducibility of a classic ecological experiment-that of assessing female host plant preference and acceptance in phytophagous insects, a group in which host specialization is a key driver of diversification. We exposed multiple cohorts of Pieris napi butterflies from the same population to traditional host acceptance and preference tests on three Brassicaceae host species. Whereas the host plant rank order was highly reproducible, the propensity to oviposit on low-ranked hosts varied significantly even among cohorts exposed to similar conditions. Much variation could be attributed to among-cohort variation in female fecundity, a trait strongly correlated both to female size and to the size of the nuptial gift a female receives during mating. Small males provide small spermatophores, and in our experiment small females that mated with small males had a disproportionally low propensity to oviposit on low-ranked hosts. Hence, our results provide empirical support to the theoretical prediction that female host utilization is strongly affected by non-genetic, environmental variation, and that such variation can affect the reproducibility of ecological experiments even under seemingly identical conditions.

Preference for C4 shade grasses increases hatchling performance in the butterfly, Bicyclus safitza

The Miocene radiation of C4 grasses under high-temperature and low ambient CO 2 levels occurred alongside the transformation of a largely forested landscape into savanna. This inevitably changed the host plant regime of herbivores, and the simultaneous diversification of many consumer lineages, including Bicyclus butterflies in Africa, suggests that the radiations of grasses and grazers may be evolutionary linked. We examined mechanisms for this plant-herbivore interaction with the grass-feeding Bicyclus safitza in South Africa. In a controlled environment, we tested oviposition preference and hatchling performance on local grasses with C3 or C4 photosynthetic pathways that grow either in open or shaded habitats. We predicted preference for C3 plants due to a hypothesized lower processing cost and higher palatability to herbivores. In contrast, we found that females preferred C4 shade grasses rather than either C4 grasses from open habitats or C3 grasses. The oviposition preference broadly followed hatchling performance, although hatchling survival was equally good on C4 or C3 shade grasses. This finding was explained by leaf toughness; shade grasses were softer than grasses from open habitats. Field monitoring revealed a preference of adults for shaded habitats, and stable isotope analysis of field-sampled individuals confirmed their preference for C4 grasses as host plants. Our findings suggest that plant-herbivore interactions can influence the direction of selection in a grass-feeding butterfly. Based on this work, we postulate future research to test whether these interactions more generally contribute to radiations in herbivorous insects via expansions into new, unexploited ecological niches.

Decoupling of female host plant preference and offspring performance in relative specialist and generalist butterflies

The preference-performance hypothesis posits that the host plant range of plant-feeding insects is ultimately limited by larval costs associated with feeding on multiple resources, and that female egg-laying preferences evolve in response to these costs. The trade-off of either using few host plant species and being a strong competitor on them due to effective utilization or using a wide host plant range but being a poor competitor is further predicted to result in host plant specialization. This follows under the hypothesis that both females and offspring are ultimately favoured by utilizing only the most suitable host(s). We develop an experimental approach to identify such trade-offs, i.e. larval costs associated with being a host generalist, and apply a suite of experiments to two sympatric and syntopic populations of the closely related butterflies Pieris napi and Pieris rapae. These butterflies show variation in their level of host specialization, which allowed comparisons between more and less specialized species and between families within species. Our results show that, first, the link between female host preference and offspring performance was not significantly stronger in the specialist compared to the generalist species. Second, the offspring of the host plant specialist did not outperform the offspring of the generalist on the former's most preferred host plant species. Finally, the more generalized species, or families within species, did not show higher survival or consistently higher growth rates than the specialists on the less preferred plants. Thus, the preference and performance traits appear to evolve as largely separated units.

Reproductive isolation and patterns of genetic differentiation in a cryptic butterfly species complex

Molecular studies of natural populations are often designed to detect and categorize hidden layers of cryptic diversity, and an emerging pattern suggests that cryptic species are more common and more widely distributed than previously thought. However, these studies are often decoupled from ecological and behavioural studies of species divergence. Thus, the mechanisms by which the cryptic diversity is distributed and maintained across large spatial scales are often unknown. In 1988, it was discovered that the common Eurasian Wood White butterfly consisted of two species (Leptidea sinapis and Leptidea reali), and the pair became an emerging model for the study of speciation and chromosomal evolution. In 2011, the existence of a third cryptic species (Leptidea juvernica) was proposed. This unexpected discovery raises questions about the mechanisms preventing gene flow and about the potential existence of additional species hidden in the complex. Here, we compare patterns of genetic divergence across western Eurasia in an extensive data set of mitochondrial and nuclear DNA sequences with behavioural data on inter- and intraspecific reproductive isolation in courtship experiments. We show that three species exist in accordance with both the phylogenetic and biological species concepts and that additional hidden diversity is unlikely to occur in Europe. The Leptidea species are now the best studied cryptic complex of butterflies in Europe and a promising model system for understanding the formation of cryptic species and the roles of local processes, colonization patterns and heterospecific interactions for ecological and evolutionary divergence.

Heterospecific courtship, minority effects and niche separation between cryptic butterfly species

Species interacting in varied ecological conditions often evolve in different directions in different local populations. The butterflies of the cryptic Leptidea complex are sympatrically distributed in different combinations across their Eurasian range. Interestingly, the same species is a habitat generalist in some regions and a habitat specialist in others, where a sibling species has the habitat generalist role. Previous studies suggest that this geographically variable niche divergence is generated by local processes in different contact zones. By varying the absolute and relative densities of Leptidea sinapis and Leptidea juvernica in large outdoor cages, we show that female mating success is unaffected by conspecific density, but strongly negatively affected by the density of the other species. Whereas 80% of the females mated when a conspecific couple was alone in a cage, less than 10% mated when the single couple shared the cage with five pairs of the other species. The heterospecific courtships can thus affect the population fitness, and for the species in the local minority, the suitability of a habitat is likely to depend on the presence or absence of the locally interacting species. If the local relative abundance of the different species depends on the colonization order, priority effects might determine the ecological roles of interacting species in this system.

Ecological constraints on female fitness in a phytophagous insect

Although understanding female reproduction is crucial for population demography, determining how and to what relative extent it is constrained by different ecological factors is complicated by difficulties in studying the links between individual behavior, life history, and fitness in nature. We present data on females in a natural population of the butterfly Leptidea sinapis. These data were combined with climate records and laboratory estimates of life-history parameters to predict the relative impact of different ecological constraints on female fitness in the wild. Using simulation models, we partitioned effects of male courtship, host plant availability, and temperature on female fitness. Results of these models indicate that temperature is the most constraining factor on female fitness, followed by host plant availability; the short-term negative effects of male courtship that were detected in the field study were less important in models predicting female reproductive success over the entire life span. In the simulations, females with more reproductive reserves were more limited by the ecological variables. Reproductive physiology and egg-laying behavior were therefore predicted to be co-optimized but reach different optima for females of different body sizes; this prediction is supported by the empirical data. This study thus highlights the need for studying behavioral and life-history variation in orchestration to achieve a more complete picture of both demographic and evolutionary processes in naturally variable and unpredictable environments.

Host-plant-induced larval decision-making in a habitat/host-plant generalist butterfly

Phenotypic plasticity can be a passive response to fluctuating environmental conditions or an active and presumably adaptive (evolved) response selected for in different environments. Here we ask if the larval decision to enter diapause when reared on a host plant associated with a colder habitat is an active or a passive response to host-plant quality or suitability. We compare plasticity in larval propensity to enter diapause of the habitat generalist butterfly Leptidea sinapis and the meadow specialist Leptidea reali in a range of temperatures and long day length on a forest plant, Lathyrus linifolius, and a meadow-associated plant, Lathyrus pratensis. The warmer meadow habitat promotes direct development whereas the colder forest habitat is conducive to diapause. Larvae of L. sinapis had a higher propensity to enter diapause when reared on the forest plant L. linifolius across all temperatures. Conversely, the propensity of L. reali to enter diapause was consistently lower and did not differ between host plants. Larval growth rates were similar between and within butterfly species and between host plants. Hence, we conclude that larval pathway decision-making in L. sinapis is an active response mediated by information from their host plants.

The evolutionary ecology of generalization: among-year variation in host plant use and offspring survival in a butterfly

The majority of phytophagous insects are relatively specialized in their food habits, and specialization in resource use is expected to be favored by selection in most scenarios. Ecological generalization is less common and less well understood, but it should be selected for by (1) rarity of resources, (2) resource inconstancy, or (3) unreliability of resource quality. Here, we test these predictions by studying egg distribution and offspring survival in the orange tip butterfly, Anthocharis cardamines, on different host plants in Sweden over a five-year period. A total of 3800 eggs were laid on 16 of the 18 crucifers available at the field site during the five years. Three main factors explained host plant generalization: (1) a rarity of food resources in which the female encounter rate of individual crucifer plants was low and within-year phenological succession of flowering periods of the different crucifers meant that individual species were suitable for oviposition only within a short time window, which translates to a low effective abundance of individual crucifer species as experienced by females searching for host plants, making specialization on a single crucifer species unprofitable; (2) variation in food resources in which among-year variation in availability of any one host plant species was high; and (3) larval survivorship varied unpredictably among years on all host plants, thereby necessitating a bet-hedging strategy and use of several different host plants. Unpredictable larval survival was caused by variation in plant stand habitat characteristics, which meant that drowning and death from starvation affected different crucifers differently, and by parasitism, which varied by host plant and year. Hence, our findings are in agreement with the theoretical explanation of ecological generalization above, helping to explain why A. cardamines is a generalist throughout its range with respect to genera within the Cruciferae.