Host plant utilization, host range oscillations and diversification in nymphalid butterflies: a phylogenetic investigation

Molecular phylogeny of Hesperiidae (Lepidoptera) with an emphasis on Asian and African genera

Despite considerable research efforts in recent years, the deeper phylogenetic relationships among skipper butterflies (Hesperiidae) remain unresolved. This is primarily because of limited sampling, especially within Asian and African lineages. In this study, we consolidated previous data and extensively sampled Asian and African taxa to elucidate the phylogenetic relationships within Hesperiidae. The molecular dataset comprised sequences from two mitochondrial and two nuclear gene regions from 563 species that represented 353 genera. Our analyses revealed seven subfamilies within Hesperiidae: Coeliadinae, Euschemoninae, Eudaminae, Pyrginae, Heteropterinae, Trapezitinae, and Hesperiinae. The systematics of most tribes and genera aligned with those of prior studies. However, notable differences were observed in several tribes and genera. Overall, the position of taxa assigned to incertae sedis in Hesperiinae is largely clarified in this study. Our results strongly support the monophyly of the tribe Tagiadini (Pyrginae), and the systematics of some genera are clarified with comprehensive discussion. We recognize 15 tribes within the subfamily Hesperiinae. Of these, nine tribes are discussed in detail: Aeromachini, Astictopterini, Erionotini, Unkanini (new status), Ancistroidini, Ismini (confirmed status), Plastingini (new status), Gretnini (confirmed status), and Eetionini (confirmed status). We propose four subtribes within Astictopterini: Hypoleucina subtrib.n., Aclerosina, Cupithina, and Astictopterina. Furthermore, we describe a new genus (Hyarotoidesgen.n.) and reinstate two genera (Zeareinst.stat. and Separeinst.stat.) as valid. Additionally, we propose several new combinations: Zea mythecacomb.n.,Sepa bononiacomb.n. & reinst.stat., and Sepa umbrosacomb.n. This study, with extensive sampling of Asian and African taxa, greatly enhances the understanding of the knowledge of the skipper tree of life.

Niche breadth explains the range size of European-centred butterflies, but dispersal ability does not

Aim

The breadth of ecological niches and dispersal abilities have long been discussed as important determinants of species' range sizes. However, studies directly comparing the relative effects of both factors are rare, taxonomically biased and revealed inconsistent results.

Location

Europe.

Time Period

Cenozoic.

Major Taxa

Butterflies, Lepidoptera.

Methods

We relate climate, diet and habitat niche breadth and two indicators of dispersal ability, wingspan and a dispersal tendency index, to the global range size of 369 European-centred butterfly species. The relative effects of these five predictors and their variation across the butterfly phylogeny were assessed by means of phylogenetic generalized least squares models and phylogenetically weighted regressions respectively.

Results

Climate niche breadth was the most important single predictor, followed by habitat and diet niche breadth, while dispersal tendency and wingspan showed no relation to species' range size. All predictors together explained 59% of the variation in butterfly range size. However, the effects of each predictor varied considerably across families and genera.

Main Conclusions

Range sizes of European-centred butterflies are strongly correlated with ecological niche breadth but apparently independent of dispersal ability. The magnitude of range size-niche breadth relationships is not stationary across the phylogeny and is often negatively correlated across the different dimensions of the ecological niche. This variation limits the generalizability of range size-trait relationships across broad taxonomic groups.

Host use by 2 sibling species of bogus yucca moths in relation to plant hardness and saponin content

Plant defenses allow plants to deter or kill their insect herbivores and are considered to be a major driver of host use for herbivorous insects in both ecological and evolutionary time. Many closely related species of insect herbivores differ in their ability to respond to plant defenses and in some cases are specialized to specific plant species. Here we tested whether both mechanical and chemical plant defenses are a major factor in determining the host range of 2 sibling species of Prodoxid bogus yucca moths, Prodoxus decipiens (Riley) and Prodoxus quinquepunctellus (Chambers) that feed within the inflorescence stalk of Yucca species. These 2 moth species have separate suites of host plant species, yet narrowly overlap geographically and share 1 Yucca species, Y. glauca. We surveyed the lignin and cellulose content, the force required to the puncture the stalk tissue, and saponin concentration across 5 Yucca species used as hosts. Lignin, cellulose concentrations, and stalk hardness differed among Yucca species but did not correlate with host use patterns by the moths. Saponin concentrations in the stalk tissue were relatively low for yuccas (<1%) and did not differ among species. The results suggest that these moth species should be able to use each other's hosts for egg deposition. Additional factors such as larval development or competition among larvae for feeding space may serve to keep moth species from expanding onto plants used by its sibling species.

Metabolic consequences of various fruit-based diets in a generalist insect species

Most phytophagous insect species exhibit a limited diet breadth and specialize on a few or a single host plant. In contrast, some species display a remarkably large diet breadth, with host plants spanning several families and many species. It is unclear, however, whether this phylogenetic generalism is supported by a generic metabolic use of common host chemical compounds ('metabolic generalism') or alternatively by distinct uses of diet-specific compounds ('multi-host metabolic specialism')? Here, we simultaneously investigated the metabolomes of fruit diets and of individuals of a generalist phytophagous species, Drosophila suzukii, that developed on them. The direct comparison of metabolomes of diets and consumers enabled us to disentangle the metabolic fate of common and rarer dietary compounds. We showed that the consumption of biochemically dissimilar diets resulted in a canalized, generic response from generalist individuals, consistent with the metabolic generalism hypothesis. We also showed that many diet-specific metabolites, such as those related to the particular color, odor, or taste of diets, were not metabolized, and rather accumulated in consumer individuals, even when probably detrimental to fitness. As a result, while individuals were mostly similar across diets, the detection of their particular diet was straightforward. Our study thus supports the view that dietary generalism may emerge from a passive, opportunistic use of various resources, contrary to more widespread views of an active role of adaptation in this process. Such a passive stance towards dietary chemicals, probably costly in the short term, might favor the later evolution of new diet specializations.

Linkage mapping and genome annotation give novel insights into gene family expansions and regional recombination rate variation in the painted lady (Vanessa cardui) butterfly

Characterization of gene family expansions and crossing over is crucial for understanding how organisms adapt to the environment. Here, we develop a high-density linkage map and detailed genome annotation of the painted lady butterfly (Vanessa cardui) - a non-diapausing, highly polyphagous species famous for its long-distance migratory behavior and almost cosmopolitan distribution. Our results reveal a complex interplay between regional recombination rate variation, gene duplications and transposable element activity shaping the genome structure of the painted lady. We identify several lineage specific gene family expansions. Their functions are mainly associated with protein and fat metabolism, detoxification, and defense against infection - critical processes for the painted lady's unique life-history. Furthermore, the detailed recombination maps allow us to characterize the regional recombination landscape, data that reveal a strong effect of chromosome size on the recombination rate, a limited impact of GC-biased gene conversion and a positive association between recombination and short interspersed elements.

A Revised Phylogenetic Classification for Viola (Violaceae)

The genus Viola (Violaceae) is among the 40–50 largest genera among angiosperms, yet its taxonomy has not been revised for nearly a century. In the most recent revision, by Wilhelm Becker in 1925, the then-known 400 species were distributed among 14 sections and numerous unranked groups. Here, we provide an updated, comprehensive classification of the genus, based on data from phylogeny, morphology, chromosome counts, and ploidy, and based on modern principles of monophyly. The revision is presented as an annotated global checklist of accepted species of Viola, an updated multigene phylogenetic network and an ITS phylogeny with denser taxon sampling, a brief summary of the taxonomic changes from Becker’s classification and their justification, a morphological binary key to the accepted subgenera, sections and subsections, and an account of each infrageneric subdivision with justifications for delimitation and rank including a description, a list of apomorphies, molecular phylogenies where possible or relevant, a distribution map, and a list of included species. We distribute the 664 species accepted by us into 2 subgenera, 31 sections, and 20 subsections. We erect one new subgenus of Viola (subg. Neoandinium, a replacement name for the illegitimate subg. Andinium), six new sections (sect. Abyssinium, sect. Himalayum, sect. Melvio, sect. Nematocaulon, sect. Spathulidium, sect. Xanthidium), and seven new subsections (subsect. Australasiaticae, subsect. Bulbosae, subsect. Clausenianae, subsect. Cleistogamae, subsect. Dispares, subsect. Formosanae, subsect. Pseudorupestres). Evolution within the genus is discussed in light of biogeography, the fossil record, morphology, and particular traits. Viola is among very few temperate and widespread genera that originated in South America. The biggest identified knowledge gaps for Viola concern the South American taxa, for which basic knowledge from phylogeny, chromosome counts, and fossil data is virtually absent. Viola has also never been subject to comprehensive anatomical study. Studies into seed anatomy and morphology are required to understand the fossil record of the genus.

Geographic isolation drives speciation in Nearctic aphids

Across herbivorous insect clades, species richness and host-use diversity tend to positively covary. This could be because host-use divergence drives speciation, or because it raises the ecological limits on species richness. To evaluate these hypotheses, we performed phylogenetic path model analyses of the species diversity of Nearctic aphids. Here, we show that variation in the species richness of aphid clades is caused mainly by host-use divergence, whereas variation in speciation rates is caused more by divergence in non-host-related niche variables. Aphid speciation is affected by both the evolution of host and non-host-related niche components, but the former is largely caused by the latter. Thus, our analyses suggest that host-use divergence can both raise the ecological limits on species richness and drive speciation, although in the latter case, host-use divergence tends to be a step along the causal path leading from non-host-related niche evolution to speciation.

Using phylogenetic path modelling, geographic isolation is revealed to be the main driver of Nearctic aphid speciation, with subsequent host-use divergence raising the limits on clade species richness.

Non-parallel morphological divergence following colonization of a new host plant

Adaptation to new ecological niches is known to spur population diversification and may lead to speciation if gene flow is ceased. While adaptation to the same ecological niche is expected to be parallel, it is more difficult to predict whether selection against maladaptive hybridization in secondary sympatry results in parallel divergence also in traits that are not directly related to the ecological niches. Such parallelisms in response to selection for reproductive isolation can be identified through estimating parallelism in reproductive character displacement across different zones of secondary contact. Here, we use a host shift in the phytophagous peacock fly Tephritis conura, with both host races represented in two geographically separate areas East and West of the Baltic Sea to investigate convergence in morphological adaptations. We asked (i) if there are consistent morphological adaptations to a host plant shift and (ii) if the response to secondary sympatry with the alternate host race is parallel across contact zones. We found surprisingly low and variable, albeit significant, divergence between host races. Only one trait, the length of the female ovipositor, which serves an important function in the interaction with the hosts, was consistently different between host races. Instead, co-existence with the other host race significantly affected the degree of morphological divergence, but the divergence was largely driven by different traits in different contact zones. Thus, local stochastic fixation or reinforcement could generate trait divergence, and additional evidence is needed to conclude whether divergence is locally adaptive.

Expanding the Menu: Are Polyphagy and Gene Family Expansions Linked across Lepidoptera?

Evolutionary expansions and contractions of gene families are often correlated with key innovations and/or ecological characteristics. In butterflies and moths (Lepidoptera), expansions of gene families involved in detoxification of plant specialized metabolites are hypothesized to facilitate a polyphagous feeding style. However, analyses supporting this hypothesis are mostly based on a limited number of lepidopteran species. We applied a phylogenomics approach, using 37 lepidopteran genomes, to analyze if gene family evolution (gene gain and loss) is associated with the evolution of polyphagy. Specifically, we compared gene counts and evolutionary gene gain and loss rates of gene families involved in adaptations with plant feeding. We correlated gene evolution to host plant family range (phylogenetic diversity) and specialized metabolite content of plant families (functional metabolite diversity). We found a higher rate for gene loss than gene gain in Lepidoptera, a potential consequence of genomic rearrangements and deletions after (potentially small-scale) duplication events. Gene family expansions and contractions varied across lepidopteran families, and were associated to host plant use and specialization levels. Within the family Noctuidae, a higher expansion rate for gene families involved in detoxification can be related to the large number of polyphagous species. However, gene family expansions are observed in both polyphagous and monophagous lepidopteran species and thus seem to be species-specific in the taxa sampled. Nevertheless, a significant positive correlation of gene counts of the carboxyl- and choline esterase and glutathione-S-transferase detoxification gene families with the level of polyphagy was identified across Lepidoptera.

Behavioral and Physiological Plasticity Provides Insights into Molecular Based Adaptation Mechanism to Strain Shift in Spodoptera frugiperda

How herbivorous insects adapt to host plants is a key question in ecological and evolutionary biology. The fall armyworm, (FAW) Spodoptera frugiperda (J.E. Smith), although polyphagous and a major pest on various crops, has been reported to have a rice and corn (maize) feeding strain in its native range in the Americas. The species is highly invasive and has recently established in China. We compared behavioral changes in larvae and adults of a corn population (Corn) when selected on rice (Rice) and the molecular basis of these adaptational changes in midgut and antennae based on a comparative transcriptome analysis. Larvae of S. frugiperda reared on rice plants continuously for 20 generations exhibited strong feeding preference for with higher larval performance and pupal weight on rice than on maize plants. Similarly, females from the rice selected population laid significantly more eggs on rice as compared to females from maize population. The most highly expressed DEGs were shown in the midgut of Rice vs. Corn. A total of 6430 DEGs were identified between the populations mostly in genes related to digestion and detoxification. These results suggest that potential adaptations for feeding on rice crops, may contribute to the current rapid spread of fall armyworm on rice crops in China and potentially elsewhere. Consistently, highly expressed DEGs were also shown in antennae; a total of 5125 differentially expressed genes (DEGs) s were identified related to the expansions of major chemosensory genes family in Rice compared to the Corn feeding population. These results not only provide valuable insight into the molecular mechanisms in host plants adaptation of S. frugiperda but may provide new gene targets for the management of this pest.

Phylogenetic reconstruction of ancestral ecological networks through time for pierid butterflies and their host plants

The study of herbivorous insects underpins much of the theory that concerns the evolution of species interactions. In particular, Pieridae butterflies and their host plants have served as a model system for studying evolutionary arms races. To learn more about the coevolution of these two clades, we reconstructed ancestral ecological networks using stochastic mappings that were generated by a phylogenetic model of host-repertoire evolution. We then measured if, when, and how two ecologically important structural features of the ancestral networks (modularity and nestedness) evolved over time. Our study shows that as pierids gained new hosts and formed new modules, a subset of them retained or recolonised the ancestral host(s), preserving connectivity to the original modules. Together, host-range expansions and recolonisations promoted a phase transition in network structure. Our results demonstrate the power of combining network analysis with Bayesian inference of host-repertoire evolution to understand changes in complex species interactions over time.

Phenotypic plasticity in chemical defence of butterflies allows usage of diverse host plants

Host plant specialization is a major force driving ecological niche partitioning and diversification in insect herbivores. The cyanogenic defences of Passiflora plants keep most herbivores at bay, but not the larvae of Heliconius butterflies, which can both sequester and biosynthesize cyanogenic compounds. Here, we demonstrate that both Heliconius cydno chioneus and H. melpomene rosina have remarkable plasticity in their chemical defences. When feeding on Passiflora species with cyanogenic compounds that they can readily sequester, both species downregulate the biosynthesis of these compounds. By contrast, when fed on Passiflora plants that do not contain cyanogenic glucosides that can be sequestered, both species increase biosynthesis. This biochemical plasticity comes at a fitness cost for the more specialist H. m. rosina, as adult size and weight for this species negatively correlate with biosynthesis levels, but not for the more generalist H. c. chioneus. By contrast, H. m rosina has increased performance when sequestration is possible on its specialized host plant. In summary, phenotypic plasticity in biochemical responses to different host plants offers these butterflies the ability to widen their range of potential hosts within the Passiflora genus, while maintaining their chemical defences.

Emerging infectious disease: An underappreciated area of strategic concern for food security

Emerging infectious diseases (EIDs) increasingly threaten global food security and public health. Despite technological breakthroughs, we are losing the battle with (re)emerging diseases as treatment costs and production losses rise. A horizon scan of diseases of crops, livestock, seafood and food-borne illness suggests these costs are unsustainable. The paradigm of coevolution between pathogens and particular hosts teaches that emerging diseases occur only when pathogens evolve specific capacities that allow them to move to new hosts. EIDs ought to be rare and unpredictable, so crisis response is the best we can do. Alternatively, the Stockholm Paradigm suggests that the world is full of susceptible but unexposed hosts that pathogens could infect, given the opportunity. Global climate change, globalized trade and travel, urbanization and land-use changes (often associated with biodiversity loss) increase those opportunities, making EID frequent. We can, however, anticipate their arrival in new locations and their behaviour once they have arrived. We can 'find them before they find us', mitigating their impacts. The DAMA (Document, Assess, Monitor, Act) protocol alters the current reactive stance and embodies proactive solutions to anticipate and mitigate the impacts of EID, extending human and material resources and buying time for development of new vaccinations, medications and control measures.

Plant phylogeny drives arboreal caterpillar assemblages across the Holarctic

Assemblages of insect herbivores are structured by plant traits such as nutrient content, secondary metabolites, physical traits, and phenology. Many of these traits are phylogenetically conserved, implying a decrease in trait similarity with increasing phylogenetic distance of the host plant taxa. Thus, a metric of phylogenetic distances and relationships can be considered a proxy for phylogenetically conserved plant traits and used to predict variation in herbivorous insect assemblages among co-occurring plant species.Using a Holarctic dataset of exposed-feeding and shelter-building caterpillars, we aimed at showing how phylogenetic relationships among host plants explain compositional changes and characteristics of herbivore assemblages.Our plant-caterpillar network data derived from plot-based samplings at three different continents included >28,000 individual caterpillar-plant interactions. We tested whether increasing phylogenetic distance of the host plants leads to a decrease in caterpillar assemblage overlap. We further investigated to what degree phylogenetic isolation of a host tree species within the local community explains abundance, density, richness, and mean specialization of its associated caterpillar assemblage.The overlap of caterpillar assemblages decreased with increasing phylogenetic distance among the host tree species. Phylogenetic isolation of a host plant within the local plant community was correlated with lower richness and mean specialization of the associated caterpillar assemblages. Phylogenetic isolation had no effect on caterpillar abundance or density. The effects of plant phylogeny were consistent across exposed-feeding and shelter-building caterpillars.Our study reveals that distance metrics obtained from host plant phylogeny are useful predictors to explain compositional turnover among hosts and host-specific variations in richness and mean specialization of associated insect herbivore assemblages in temperate broadleaf forests. As phylogenetic information of plant communities is becoming increasingly available, further large-scale studies are needed to investigate to what degree plant phylogeny structures herbivore assemblages in other biomes and ecosystems.

Plant Chemistry Determines Host Preference and Performance of an Invasive Insect

Understanding how host plant chemistry affects invasive insects is crucial for determining the physiological mechanism of host use and predicting invasive insect outbreak and damage on hosts. Here, we examined the effects of plant nutrition and defensive chemicals on host preference and performance of adults and larvae of the invasive potato tuberworm, Phthorimaea operculella (Zeller; Lepidoptera: Gelechiidae), on four native (Solanum tuberosum, Nicotiana tabacum, Datura stramonium, and Solanum lycopersicum) and three new (Solanum melongena, Physalis alkekengi, and Lycium barbarum) host plants. We found that adults preferred to oviposit on S. tuberosum and N. tabacum leaves and the soil around these native host plants over other hosts. Larvae performed well on S. tuberosum and N. tabacum, reaching higher pupa weight and having better survival. Larvae performed poorly on S. melongena, S. lycopersicum, P. alkekengi, D. stramonium, and L. barbarum, with lower pupa weight and lower survival. Solanum tuberosum and N. tabacum had higher leaf soluble proteins than other plants and lower leaf total phenolics than S. lycopersicum, D. stramonium, and L. barbarum. Moreover, carbon content and soluble protein were positively associated with larval survival, while defensive traits (lignin and total phenolics) negatively affected larval survival. These findings provide insights into understanding of biochemical mechanisms of interactions between invasive insects and host plants, indicating the importance of considering plant chemistry when assessing invasive insect host use and damage.

Bayesian Inference of Ancestral Host-Parasite Interactions under a Phylogenetic Model of Host Repertoire Evolution

Intimate ecological interactions, such as those between parasites and their hosts, may persist over long time spans, coupling the evolutionary histories of the lineages involved. Most methods that reconstruct the coevolutionary history of such interactions make the simplifying assumption that parasites have a single host. Many methods also focus on congruence between host and parasite phylogenies, using cospeciation as the null model. However, there is an increasing body of evidence suggesting that the host ranges of parasites are more complex: that host ranges often include more than one host and evolve via gains and losses of hosts rather than through cospeciation alone. Here, we develop a Bayesian approach for inferring coevolutionary history based on a model accommodating these complexities. Specifically, a parasite is assumed to have a host repertoire, which includes both potential hosts and one or more actual hosts. Over time, potential hosts can be added or lost, and potential hosts can develop into actual hosts or vice versa. Thus, host colonization is modeled as a two-step process that may potentially be influenced by host relatedness. We first explore the statistical behavior of our model by simulating evolution of host-parasite interactions under a range of parameter values. We then use our approach, implemented in the program RevBayes, to infer the coevolutionary history between 34 Nymphalini butterfly species and 25 angiosperm families. Our analysis suggests that host relatedness among angiosperm families influences how easily Nymphalini lineages gain new hosts. [Ancestral hosts; coevolution; herbivorous insects; probabilistic modeling.].

Tolerance of Novel Toxins through Generalized Mechanisms: Simulating Gradual Host Shifts of Butterflies

Organisms encounter a wide range of toxic compounds in their environments, from chemicals that serve anticonsumption or anticompetition functions to pollutants and pesticides. Although we understand many detoxification mechanisms that allow organisms to consume toxins typical of their diet, we know little about why organisms vary in their ability to tolerate entirely novel toxins. We tested whether variation in generalized stress responses, such as antioxidant pathways, may underlie variation in reactions to novel toxins and, if so, their associated costs. We used an artificial diet to present cabbage white butterfly caterpillars (Pieris rapae) with plant material containing toxins not experienced in their evolutionary history. Families that maintained high performance (e.g., high survival, fast development time, large body size) on diets containing one novel toxic plant also performed well when exposed to two other novel toxic plants, consistent with a generalized response. Variation in constitutive (but not induced) expression of genes involved in oxidative stress responses was positively related to performance on the novel diets. While we did not detect reproductive trade-offs of this generalized response, there was a tendency to have less melanin investment in the wings, consistent with the role of melanin in oxidative stress responses. Taken together, our results support the hypothesis that variation in generalized stress responses, such as genes involved in oxidative stress responses, may explain the variation in tolerance to entirely novel toxins and may facilitate colonization of novel hosts and environments.

Potential for Acanthoscelides obtectus to Adapt to New Hosts Seen in Laboratory Selection Experiments

Effective pest management strategies for a targeted pest species must rely on accurate, reliable and reproducible estimates of population dynamics. Importance of such approaches is even more conspicuous when assessing pest’s potential to utilize other stored products. Using an experimental evolution approach, we have focused our attention on a common bean pest, the seed beetle (Acanthoscelides obtectus). We looked into the potential to invade and sustain population growth on two suboptimal host plants (chickpeas and mung beans). Such an approach simulates steps of the host-shift process in storages. By analyzing population dynamics during initial encountering with a new host plant, we detected a population drop for both novel hosts. However, transgenerational development in a novel environment resulted in a constant population growth in chickpeas, but not in mung bean populations. Reversal of chickpea selected populations to original host plant has led to a severe decrease in population parameters due to low viability of immatures, while the opposite trend was detected in mung bean populations. This paper highlights the importance of good practice in estimating population dynamics for economically important species. With special emphasis on storage pest species, we discuss how this approach can be useful for estimating invading potential of pest insects.

Unifying host-associated diversification processes using butterfly-plant networks

Explaining the exceptional diversity of herbivorous insects is an old problem in evolutionary ecology. Here we focus on the two prominent hypothesised drivers of their diversification, radiations after major host switch or variability in host use due to continuous probing of new hosts. Unfortunately, current methods cannot distinguish between these hypotheses, causing controversy in the literature. Here we present an approach combining network and phylogenetic analyses, which directly quantifies support for these opposing hypotheses. After demonstrating that each hypothesis produces divergent network structures, we then investigate the contribution of each to diversification in two butterfly families: Pieridae and Nymphalidae. Overall, we find that variability in host use is essential for butterfly diversification, while radiations following colonisation of a new host are rare but can produce high diversity. Beyond providing an important reconciliation of alternative hypotheses for butterfly diversification, our approach has potential to test many other hypotheses in evolutionary biology.

Unifying host-associated diversification processes using butterfly-plant networks

Explaining the exceptional diversity of herbivorous insects is an old problem in evolutionary ecology. Here we focus on the two prominent hypothesised drivers of their diversification, radiations after major host switch or variability in host use due to continuous probing of new hosts. Unfortunately, current methods cannot distinguish between these hypotheses, causing controversy in the literature. Here we present an approach combining network and phylogenetic analyses, which directly quantifies support for these opposing hypotheses. After demonstrating that each hypothesis produces divergent network structures, we then investigate the contribution of each to diversification in two butterfly families: Pieridae and Nymphalidae. Overall, we find that variability in host use is essential for butterfly diversification, while radiations following colonisation of a new host are rare but can produce high diversity. Beyond providing an important reconciliation of alternative hypotheses for butterfly diversification, our approach has potential to test many other hypotheses in evolutionary biology.

Herbivorous insects could diversify through radiations after major host switches or through constant variability in new host use. With phylogenetic and network analyses, Braga et al. show that variability in host use supports most butterfly diversification, while rare radiations can further boost diversity.

Pollen metabarcoding as a tool for tracking long-distance insect migrations

Insects account for a large portion of Earth's biodiversity and are key players for ecosystems, notably as pollinators. While insect migration is suspected to represent a natural phenomenon of major importance, remarkably little is known about it, except for a few flagship species. The reason for this situation is mainly due to technical limitations in the study of insect movement. Here, we propose using metabarcoding of pollen carried by insects as a method for tracking their migrations. We developed a flexible and simple protocol allowing efficient multiplexing and not requiring DNA extraction, one of the most time-consuming part of metabarcoding protocols, and apply this method to the study of the long-distance migration of the butterfly Vanessa cardui, an emerging model for insect migration. We collected 47 butterfly samples along the Mediterranean coast of Spain in spring and performed metabarcoding of pollen collected from their bodies to test for potential arrivals from the African continent. In total, we detected 157 plant species from 23 orders, most of which (82.8%) were insect-pollinated. Taxa present in Africa-Arabia represented 73.2% of our data set, and 19.1% were endemic to this region, strongly supporting the hypothesis that migratory butterflies colonize southern Europe from Africa in spring. Moreover, our data suggest that a northwards trans-Saharan migration in spring is plausible for early arrivals (February) into Europe, as shown by the presence of Saharan floristic elements. Our results demonstrate the possibility of regular insect-mediated transcontinental pollination, with potential implications for ecosystem functioning, agriculture and plant phylogeography.

Evolution of Gustatory Receptor Gene Family Provides Insights into Adaptation to Diverse Host Plants in Nymphalid Butterflies

The host plant range of herbivorous insects is a major aspect of insect–plant interaction, but the genetic basis of host range expansion in insects is poorly understood. In butterflies, gustatory receptor genes (GRs) play important roles in host plant selection by ovipositing females. Since several studies have shown associations between the repertoire sizes of chemosensory gene families and the diversity of resource use, we hypothesized that the increase in the number of genes in the GR family is associated with host range expansion in butterflies. Here, we analyzed the evolutionary dynamics of GRs among related species, including the host generalist Vanessa cardui and three specialists. Although the increase of the GR repertoire itself was not observed, we found that the gene birth rate of GRs was the highest in the lineage leading to V. cardui compared with other specialist lineages. We also identified two taxon-specific subfamilies of GRs, characterized by frequent lineage-specific duplications and higher non-synonymous substitution rates. Together, our results suggest that frequent gene duplications in GRs, which might be involved in the detection of plant secondary metabolites, were associated with host range expansion in the V. cardui lineage. These evolutionary patterns imply that the capability to perceive various compounds during host selection was favored during adaptation to diverse host plants.

Host specificity of avian haemosporidian parasites is unrelated among sister lineages but shows phylogenetic signal across larger clades

Parasites can vary in the number of host species they infect, a trait known as "host specificity". Here we quantify phylogenetic signal-the tendency for closely related species to resemble each other more than distantly related species-in host specificity of avian haemosporidian parasites (genera Plasmodium, Haemoproteus and Leucocytozoon) using data from MalAvi, the global avian haemosporidian database. We used the genetic data (479 base pairs of cytochrome b) that define parasite lineages to produce genus level phylogenies. Combining host specificity data with those phylogenies revealed significant levels of phylogenetic signal while controlling for sampling effects; phylogenetic signal was higher when the phylogenetic diversity of hosts was taken into account. We then tested for correlations in the host specificity of pairs of sister lineages. Correlations were generally close to zero for all three parasite genera. These results suggest that while the host specificity of parasite sister lineages differ, larger clades may be relatively specialised or generalised.

The complete mitochondrial genome of Vanessa indica and phylogenetic analyses of the family Nymphalidae

Vanessa indica is a small butterfly lacking historical molecular and biological research. Vanessa indica belongs to the family Nymphalidae (Lepidoptera: Papilionoidea), which is the largest group of butterflies and are nearly ubiquitous. However, after more than a century of taxonomic and molecular studies, there is no consensus for family classification, and the phylogenetic relationships within Nymphalidae are controversial. The first objective was to sequence and characterize the complete mitochondrial genome of V. indica. The most important objective was to completely reconstruct the phylogenetic relationships for family members within Nymphalidae. The mitochondrial genomic DNA (mtDNA) of V. indica was extracted and amplified by polymerase chain reaction. The complete mitochondrial sequence was annotated and characterized by analyzing sequences with SeqMan program. The phylogenetic analyses were conducted on thirteen protein coding genes (PCGs) in 95 mtDNA of Nymphalidae downloaded from GenBank for reference using the maximum likelihood method and Bayesian inference to ensure the validity of the results. The complete mitogenome was a circular molecule with 15,191 bp consisting of 13 protein coding genes, two ribosomal RNA genes (16S rRNA and 12S rRNA), 22 transfer RNA (tRNA) genes, and an A + T-rich region (D-loop). The nucleotide composition of the genome was highly biased for A + T content, which accounts for 80.0% of the nucleotides. All the tRNAs have putative secondary structures that are characteristic of mitochondrial tRNAs, except tRNA^Ser(AGN). All the PCGs started with ATN codons, except cytochrome c oxidase subunit 1 (COX1), which was found to start with an unusual CGA codon. Four genes were observed to have unusual codons: COX1 terminated with atypical TT and the other three genes terminated with a single T. The A + T rich region of 327 bp consisted of repetitive sequences, including a ATAGA motif, a 19-bp poly-T stretch, and two microsatellite-like regions (TA)₈. The phylogenetic analyses consistently placed Biblidinae as a sister cluster to Heliconiinae and Calinaginae as a sister clade to Satyrinae. Moreover, the phylogenetic tree identified Libytheinae as a monophyletic group within Nymphalidae. The complete mitogenome of V. indica was 15,191 bp with mitochondrial characterizations common for lepidopteran species, which enriched the mitochondria data of Nymphalid species. And the phylogenetic analysis revealed different classifications and relationships than those previously described. Our results are significant because they would be useful in further understanding of the evolutionary biology of Nymphalidae.

Mysteries of host switching: Diversification and host specificity in rodent-coccidia associations

Recent studies show that host switching is much more frequent than originally believed and constitutes an important driver in evolution of host-parasite associations. However, its frequency and ecological mechanisms at the population level have been rarely investigated. We address this issue by analyzing phylogeny and population genetics of an extensive sample, from a broad geographic area, for commonly occurring parasites of the genus Eimeria within the abundant rodent genera Apodemus, Microtus and Myodes, using two molecular markers. At the most basal level, we demonstrate polyphyletic arrangement, i.e. multiple origin, of the rodent-specific clusters within the Eimeria phylogeny, and strong genetic/phylogenetic structure within these lineages determined at least partially by specificities to different host groups. However, a novel and the most important observation is a repeated occurrence of host switches among closely related genetic lineages which may become rapidly fixed. Within the studied model, this phenomenon applies particularly to the switches between the eimerians from Apodemus flavicollis/Apodemus sylvaticus and Apodemus agrarius groups. We show that genetic differentiation and isolation between A. flavicollis/A. sylvaticus and A. agrarius faunas is a secondary recent event and does not reflect host-parasite coevolutionary history. Rather, it provides an example of rapid ecology-based differentiation in the parasite population.

Variably hungry caterpillars: predictive models and foliar chemistry suggest how to eat a rainforest

A long-term goal in evolutionary ecology is to explain the incredible diversity of insect herbivores and patterns of host plant use in speciose groups like tropical Lepidoptera. Here, we used standardized food-web data, multigene phylogenies of both trophic levels and plant chemistry data to model interactions between Lepidoptera larvae (caterpillars) from two lineages (Geometridae and Pyraloidea) and plants in a species-rich lowland rainforest in New Guinea. Model parameters were used to make and test blind predictions for two hectares of an exhaustively sampled forest. For pyraloids, we relied on phylogeny alone and predicted 54% of species-level interactions, translating to 79% of all trophic links for individual insects, by sampling insects from only 15% of local woody plant diversity. The phylogenetic distribution of host-plant associations in polyphagous geometrids was less conserved, reducing accuracy. In a truly quantitative food web, only 40% of pair-wise interactions were described correctly in geometrids. Polyphenol oxidative activity (but not protein precipitation capacity) was important for understanding the occurrence of geometrids (but not pyraloids) across their hosts. When both foliar chemistry and plant phylogeny were included, we predicted geometrid-plant occurrence with 89% concordance. Such models help to test macroevolutionary hypotheses at the community level.

Polyphagy and diversification in tussock moths: Support for the oscillation hypothesis from extreme generalists

Theory on plasticity driving speciation, as applied to insect-plant interactions (the oscillation hypothesis), predicts more species in clades with higher diversity of host use, all else being equal. Previous support comes mainly from specialized herbivores such as butterflies, and plasticity theory suggests that there may be an upper host range limit where host diversity no longer promotes diversification. The tussock moths (Erebidae: Lymantriinae) are known for extreme levels of polyphagy. We demonstrate that this system is also very different from butterflies in terms of phylogenetic signal for polyphagy and for use of specific host orders. Yet we found support for the generality of the oscillation hypothesis, in that clades with higher diversity of host use were found to contain more species. These clades also consistently contained the most polyphagous single species. Comparing host use in Lymantriinae with related taxa shows that the taxon indeed stands out in terms of the frequency of polyphagous species. Comparative evidence suggests that this is most probably due to its nonfeeding adults, with polyphagy being part of a resulting life history syndrome. Our results indicate that even high levels of plasticity can drive diversification, at least when the levels oscillate over time.

A Neutral Model for the Evolution of Diet Breadth

Variation in diet breadth among organisms is a pervasive feature of the natural world that has resisted general explanation. In particular, trade-offs in the ability to use one resource at the expense of another have been expected but rarely detected. We explore a spatial model for the evolution of specialization, motivated by studies of plant-feeding insects. The model is neutral with respect to the causes and consequences of diet breadth: the number of hosts utilized is not constrained by trade-offs, and specialization or generalization does not confer a direct advantage with respect to the persistence of populations or the probability of diversification. We find that diet breadth evolves in ways that resemble reports from natural communities. Simulated communities are dominated by specialized species, with a predictable but less species-rich component of generalized taxa. These results raise the possibility that specialization might be a consequence of stochastic diversification dynamics acting on spatially segregated consumer-resource associations rather than a trait either favored or constrained directly by natural selection. Finally, our model generates hypotheses for global patterns of herbivore diet breadth, including a positive effect of host richness and a negative effect of evenness in host plant abundance on the number of specialized taxa.

Phylogeny and Evolution of Lepidoptera

Until recently, deep-level phylogeny in Lepidoptera, the largest single radiation of plant-feeding insects, was very poorly understood. Over the past two decades, building on a preceding era of morphological cladistic studies, molecular data have yielded robust initial estimates of relationships both within and among the ∼43 superfamilies, with unsolved problems now yielding to much larger data sets from high-throughput sequencing. Here we summarize progress on lepidopteran phylogeny since 1975, emphasizing the superfamily level, and discuss some resulting advances in our understanding of lepidopteran evolution.

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.

Nutrition shapes life-history evolution across species

Nutrition is a key component of life-history theory, yet we know little about how diet quality shapes life-history evolution across species. Here, we test whether quantitative measures of nutrition are linked to life-history evolution across 96 species of butterflies representing over 50 independent diet shifts. We find that butterflies feeding on high nitrogen host plants as larvae are more fecund, but their eggs are smaller relative to their body size. Nitrogen and sodium content of host plants are also both positively related to eye size. Some of these relationships show pronounced lineage-specific effects. Testis size is not related to nutrition. Additionally, the evolutionary timing of diet shifts is not important, suggesting that nutrition affects life histories regardless of the length of time a species has been adapting to its diet. Our results suggest that, at least for some lineages, species with higher nutrient diets can invest in a range of fitness-related traits like fecundity and eye size while allocating less to each egg as offspring have access to a richer diet. These results have important implications for the evolution of life histories in the face of anthropogenic changes in nutrient availability.

Why stay in a bad relationship? The effect of local host phenology on a generalist butterfly feeding on a low-ranked host

BACKGROUND

In plant-feeding insects, the evolutionary retention of polyphagy remains puzzling. A better understanding of the relationship between these organisms and changes in the metabolome of their host plants is likely to suggest functional links between them, and may provide insights into how polyphagy is maintained.

RESULTS

We investigated the phenological change of Cynoglossum officinale, and how a generalist butterfly species, Vanessa cardui, responded to this change. We used untargeted metabolite profiling to map plant seasonal changes in both primary and secondary metabolites. We compared these data to differences in larval performance on vegetative plants early and late in the season. We also performed two oviposition preference experiments to test females' ability to choose between plant developmental stages (vegetative and reproductive) early and late in the season. We found clear seasonal changes in plant primary and secondary metabolites that correlated with larval performance. The seasonal change in plant metabolome reflected changes in both nutrition and toxicity and resulted in zero survival in the late period. However, large differences among families in larval ability to feed on C. officinale suggest that there is genetic variation for performance on this host. Moreover, females accepted all plants for oviposition, and were not able to discriminate between plant developmental stages, in spite of the observed overall differences in metabolite profile potentially associated with differences in suitability as larval food.

CONCLUSIONS

In V. cardui, migratory behavior, and thus larval feeding times, are not synchronized with plant phenology at the reproductive site. This lack of synchronization, coupled with the observed lack of discriminatory oviposition, obviously has potential fitness costs. However, this "opportunistic" behavior may as well function as a source of potential host plant evolution, promoting for example the acceptance of new plants.

Unravelling the role of host plant expansion in the diversification of a Neotropical butterfly genus

BACKGROUND

Understanding the processes underlying diversification is a central question in evolutionary biology. For butterflies, access to new host plants provides opportunities for adaptive speciation. On the one hand, locally abundant host species can generate ecologically significant selection pressure. But a diversity of host plant species within the geographic range of each population and/or species might also eliminate any advantage conferred by specialization. This paper focuses on four Melinaea species, which are oligophagous on the family Solanaceae: M. menophilus, M. satevis, M. marsaeus, and finally, M. mothone. We examined both female preference and larval performance on two host plant species that commonly occur in this butterfly's native range, Juanulloa parasitica and Trianaea speciosa, to determine whether the different Melinaea species show evidence of local adaptation.

RESULT

In choice experiments, M. mothone females used both host plants for oviposition, whereas all other species used J. parasitica almost exclusively. In no choice experiment, M. mothone was the only species that readily accepted T. speciosa as a larval host plant. Larval survival was highest on J. parasitica (82.0 % vs. 60.9 %) and development took longer on T. speciosa (14.12 days vs. 13.35 days), except for M. mothone, which did equally well on both host plants. For all species, average pupal weight was highest on J. parasitica (450.66 mg vs. 420.01 mg), although this difference was least apparent in M. mothone.

CONCLUSION

We did not find that coexisting species of Melinaea partition host plant resources as expected if speciation is primarily driven by host plant divergence. Although M. mothone shows evidence of local adaptation to a novel host plant, T. speciosa, which co-occurs, it does not preferentially lay more eggs on or perform better on this host plant than on host plants used by other Melinaea species and not present in its distributional range. It is likely that diversification in this genus is driven by co-occurring Müllerian mimics and the resulting predation pressure, although this is also likely made possible by greater niche diversity as a consequence of plasticity for potential hosts.

Diet dependent metabolic responses in three generalist insect herbivores Spodoptera spp

Adaption to dietary changes is critical in the evolution of host plant ranges in polyphagous insects. We compared three taxa of lepidopteran herbivores from the predominantly generalist genus Spodoptera showing different degrees of polyphagy: Spodoptera littoralis, with a broad host range including both mono- and dicotyledonous plants, and two Spodoptera frugiperda strains [Corn (i.e. maize) (C) and Rice (R)] adapted primarily to different grass species. When feeding on maize we show a lower performance in the broad generalist taxon compared to the grass adapted taxa. Among these taxa, the maize adapted S. frugiperda C-strain generally performed better than the R-strain on maize leaves. On artificial pinto diet, all taxa performed well. Our RNA-Seq analysis of midgut transcriptomes from 3rd instar larvae feeding on maize showed broader transcriptional readjustments in the generalist S. littoralis compared to grass adapted S. frugiperda strains. Substantial alteration in the expression levels of midgut physiological function related transcripts, such as digestive and detoxifying enzymes, transporters, immunity, and peritrophic membrane associated transcripts, existed in all taxa. We found high background expression of UDP-glucosyl transferases, which are known to neutralize maize leaf toxins, in the maize adapted S. frugiperda C-strain, contributing to its fitness on maize compared to the R-strain. Our findings provide evidence for divergent diet specific response of digestive physiology within these Spodoptera taxa. Unexpectedly, the C- and R-strains of S. frugiperda fed on the same diet showed large differences in expression patterns between these two closely related taxa.

Evolutionary history of host use, rather than plant phylogeny, determines gene expression in a generalist butterfly

BACKGROUND

Although most insect species are specialized on one or few groups of plants, there are phytophagous insects that seem to use virtually any kind of plant as food. Understanding the nature of this ability to feed on a wide repertoire of plants is crucial for the control of pest species and for the elucidation of the macroevolutionary mechanisms of speciation and diversification of insect herbivores. Here we studied Vanessa cardui, the species with the widest diet breadth among butterflies and a potential insect pest, by comparing tissue-specific transcriptomes from caterpillars that were reared on different host plants. We tested whether the similarities of gene-expression response reflect the evolutionary history of adaptation to these plants in the Vanessa and related genera, against the null hypothesis of transcriptional profiles reflecting plant phylogenetic relatedness.

RESULT

Using both unsupervised and supervised methods of data analysis, we found that the tissue-specific patterns of caterpillar gene expression are better explained by the evolutionary history of adaptation of the insects to the plants than by plant phylogeny.

CONCLUSION

Our findings suggest that V. cardui may use two sets of expressed genes to achieve polyphagy, one associated with the ancestral capability to consume Rosids and Asterids, and another allowing the caterpillar to incorporate a wide range of novel host-plants.

Understanding Host-Switching by Ecological Fitting

Despite the fact that parasites are highly specialized with respect to their hosts, empirical evidence demonstrates that host switching rather than co-speciation is the dominant factor influencing the diversification of host-parasite associations. Ecological fitting in sloppy fitness space has been proposed as a mechanism allowing ecological specialists to host-switch readily. That proposal is tested herein using an individual-based model of host switching. The model considers a parasite species exposed to multiple host resources. Through time host range expansion can occur readily without the prior evolution of novel genetic capacities. It also produces non-linear variation in the size of the fitness space. The capacity for host colonization is strongly influenced by propagule pressure early in the process and by the size of the fitness space later. The simulations suggest that co-adaptation may be initiated by the temporary loss of less fit phenotypes. Further, parasites can persist for extended periods in sub-optimal hosts, and thus may colonize distantly related hosts by a "stepping-stone" process.

The diversification of Heliconius butterflies: what have we learned in 150 years?

Research into Heliconius butterflies has made a significant contribution to evolutionary biology. Here, we review our understanding of the diversification of these butterflies, covering recent advances and a vast foundation of earlier work. Whereas no single group of organisms can be sufficient for understanding life's diversity, after years of intensive study, research into Heliconius has addressed a wide variety of evolutionary questions. We first discuss evidence for widespread gene flow between Heliconius species and what this reveals about the nature of species. We then address the evolution and diversity of warning patterns, both as the target of selection and with respect to their underlying genetic basis. The identification of major genes involved in mimetic shifts, and homology at these loci between distantly related taxa, has revealed a surprising predictability in the genetic basis of evolution. In the final sections, we consider the evolution of warning patterns, and Heliconius diversity more generally, within a broader context of ecological and sexual selection. We consider how different traits and modes of selection can interact and influence the evolution of reproductive isolation.