Patterns of morphological variation highlight the effect of natural selection on eyespots modularity in the butterfly Morpho telemachus.

Morphological correlations can not only stem from developmental constraints but also from selective pressures. Butterfly eyespots are repeated wing color pattern elements, widespread across species. As developmental serial homologs, they are controlled by similar developmental pathways imposing correlations among eyespots: selection on a single eyespot may induce correlated responses in all eyespots. We study the variations in the ventral eyespots of Morpho telemachus, where two different selective regimes are likely to act: while most eyespots are always-visible, two eyespots are conditionally displayed: hidden at rest, they can be exposed when the butterflies are threatened, or during sexual interactions. We investigate how such contrasted selection across eyespots can alter the covariations imposed by their shared developmental origin. We quantified eyespots covariations within a large population of M. telemachus and compared the observed patterns to those found in M. helenor, where all eyespots are always-visible and thus probably affected by a similar selection regime. We found that M. telemachus conditionally displayed eyespots are less variable than always-visible eyespots and that these two eyespots form a separate variational module in this species, in contrast to M. helenor. Our results suggest that eyespots covariations were shaped by selection, highlighting how natural selection may promote the evolution of modularity.

Long-term evolution of quantitative traits in the Drosophila melanogaster species subgroup.

Quantitative genetics aims at untangling the genetic and environmental effects on phenotypic variation. Trait heritability, which summarizes the relative importance of genetic effects, is estimated at the intraspecific level, but theory predicts that heritability could influence long-term evolution of quantitative traits. The phylogenetic signal concept bears resemblance to heritability and it has often been called species-level heritability. Under certain conditions, such as trait neutrality or contribution to phylogenesis, within-species heritability and between-species phylogenetic signal should be correlated. Here, we investigate the potential relationship between these two concepts by examining the evolution of multiple morphological traits for which heritability has been estimated in Drosophila melanogaster. Specifically, we analysed 42 morphological traits in both sexes on a phylogeny inferred from 22 nuclear genes for nine species of the melanogaster subgroup. We used Pagel's λ as a measurement of phylogenetic signal because it is the least influenced by the number of analysed taxa. Pigmentation traits showed the strongest concordance with the phylogeny, but no correlation was found between phylogenetic signal and heritability estimates mined from the literature. We obtained data for multiple climatic variables inferred from the geographical distribution of each species. Phylogenetic regression of quantitative traits on climatic variables showed a significantly positive correlation with heritability. Convergent selection, the response to which depends on the trait heritability, may have led to the null association between phylogenetic signal and heritability for morphological traits in Drosophila. We discuss the possible causes of discrepancy between both statistics and caution against their confusion in evolutionary biology.

Divergence of climbing escape flight performance in Morpho butterflies living in different microhabitats.

Habitat specialization can influence the evolution of animal movement in promoting divergent locomotor abilities adapted to contrasting environmental conditions, differences in vegetation clutter or predatory communities. While the effect of habitat on the evolution of locomotion and particularly escape performance has been well investigated in terrestrial animals, it remains understudied in flying animals. Here, we investigated whether specialization of Morpho butterfly species into different vertical strata of the Amazonian forest affects the performance of upward escape flight manoeuvres. Using stereoscopic high-speed videography, we compared the climbing flight kinematics of seven Morpho species living either in the forest canopy or in the understory. We show that butterflies from canopy species display strikingly higher climbing speed and steeper ascent angle compared with understory species. Although climbing speed increased with wing speed and angle of attack, the higher climb angle observed in canopy species was best explained by their higher body pitch angle, resulting in more upward-directed aerodynamic thrust forces. Climb angle also scales positively with weight-normalized wing area, and this weight-normalized wing area was higher in canopy species. This shows that a combined divergence in flight behaviour and morphology contributes to the evolution of increased climbing flight abilities in canopy species.

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

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

The effect of captivity on craniomandibular and calcaneal ontogenetic trajectories in wild boar.

Deciphering the plastic (i.e., nonheritable) changes induced by human control over wild animals in the archeological record is challenging. Previous studies detected morphological markers associated with captivity in the cranium, mandible, and calcaneus of adult wild boar (Sus scrofa) but the developmental trajectories leading up to these changes during ontogeny remain unknown. To assess the impact of growth in a captive environment on morphological structures during postnatal ontogeny, we used an experimental approach focusing on the same three structures and taxon. We investigated the form and size differences of captive-reared and wild-caught wild boar during growth using three-dimensional landmark-based geometric morphometrics. Our results provide evidence of an influence of captivity on the morphology of craniomandibular structures, as wild specimens are smaller than captive individuals at similar ages. The food resources inherent to anthropogenic environments may explain some of the observed differences between captive-reared and wild specimens. The calcaneus presents a different contrasted pattern of plasticity as captive and wild individuals differ in terms of form but not in terms of size. The physically more constrained nature of the calcaneus and the direct influence of mobility reduction on this bone may explain these discrepancies. These results provide new methodological perspectives for bioarchaeological approaches as they imply that the plastic mark of captivity can be observed in juvenile specimens in the same way it has been previously described in adults.

Genome assembly of 3 Amazonian Morpho butterfly species reveals Z-chromosome rearrangements between closely related species living in sympatry.

The genomic processes enabling speciation and species coexistence in sympatry are still largely unknown. Here we describe the whole-genome sequencing and assembly of 3 closely related species from the butterfly genus Morpho: Morpho achilles (Linnaeus, 1758), Morpho helenor (Cramer, 1776), and Morpho deidamia (Höbner, 1819). These large blue butterflies are emblematic species of the Amazonian rainforest. They live in sympatry in a wide range of their geographical distribution and display parallel diversification of dorsal wing color pattern, suggesting local mimicry. By sequencing, assembling, and annotating their genomes, we aim at uncovering prezygotic barriers preventing gene flow between these sympatric species. We found a genome size of  480 Mb for the 3 species and a chromosomal number ranging from 2n = 54 for M. deidamia to 2n = 56 for M. achilles and M. helenor. We also detected inversions on the sex chromosome Z that were differentially fixed between species, suggesting that chromosomal rearrangements may contribute to their reproductive isolation. The annotation of their genomes allowed us to recover in each species at least 12,000 protein-coding genes and to discover duplications of genes potentially involved in prezygotic isolation like genes controlling color discrimination (L-opsin). Altogether, the assembly and the annotation of these 3 new reference genomes open new research avenues into the genomic architecture of speciation and reinforcement in sympatry, establishing Morpho butterflies as a new eco-evolutionary model.

Convergent morphology and divergent phenology promote the coexistence of Morpho butterfly species.

The coexistence of closely-related species in sympatry is puzzling because ecological niche proximity imposes strong competition and reproductive interference. A striking example is the widespread wing pattern convergence of several blue-banded Morpho butterfly species with overlapping ranges of distribution. Here we perform a series of field experiments using flying Morpho dummies placed in a natural habitat. We show that similarity in wing colour pattern indeed leads to interspecific territoriality and courtship among sympatric species. In spite of such behavioural interference, demographic inference from genomic data shows that sympatric closely-related Morpho species are genetically isolated. Mark-recapture experiments in the two most closely-related species unravel a strong temporal segregation in patrolling activity of males. Such divergence in phenology reduces the costs of reproductive interference while simultaneously preserving the benefits of convergence in non-reproductive traits in response to common ecological pressures. Henceforth, the evolution of multiple traits may favour species diversification in sympatry by partitioning niche in different dimensions.

Adaptive evolution of flight in Morpho butterflies.

The diversity of flying animals suggests that countless combinations of flight morphologies and behaviors have evolved with specific lifestyles, thereby exploiting diverse aerodynamic mechanisms. How morphology, flight behavior, and aerodynamic properties together diversify with contrasting ecology remains to be elucidated. We studied the adaptive codivergence in wing shape, flight behavior, and aerodynamic efficiency among Morpho butterflies living in different forest strata by combining high-speed videography in the field with morphometric analyses and aerodynamic modeling. By comparing canopy and understory species, we show that adaptation to an open canopy environment resulted in increased glide efficiency. Moreover, this enhanced glide efficiency was achieved by different canopy species through distinct combinations of flight behavior, wing shape, and aerodynamic mechanisms, highlighting the multiple pathways of adaptive evolution.

Punctuational ecological changes rather than global factors drive species diversification and the evolution of wing phenotypes in Morpho butterflies.

Assessing the relative importance of geographical and ecological drivers of evolution is paramount to understand the diversification of species and traits at the macroevolutionary scale. Here, we use an integrative approach, combining phylogenetics, biogeography, ecology and quantified phenotypes to investigate the drivers of both species and phenotypic diversification of the iconic Neotropical butterfly genus Morpho. We generated a time-calibrated phylogeny for all known species and inferred historical biogeography. We fitted models of time-dependent (accounting for rate heterogeneity across the phylogeny) and paleoenvironment-dependent diversification (accounting for global effect on the phylogeny). We used geometric morphometrics to assess variation of wing size and shape across the tree and investigated their dynamics of evolution. We found that the diversification of Morpho is best explained when considering variable diversification rates across the tree, possibly associated with lineages occupying different microhabitat conditions. First, a shift from understory to canopy was characterized by an increased speciation rate partially coupled with an increasing rate of wing shape evolution. Second, the occupation of dense bamboo thickets accompanying a major host-plant shift from dicotyledons towards monocotyledons was associated with a simultaneous diversification rate shift and an evolutionary 'jump' of wing size. Our study points to a diversification pattern driven by punctuational ecological changes instead of a global driver or biogeographic history.

Phenotypic plasticity, canalisation and developmental stability of Triatoma infestans wings: effects of a sublethal application of a pyrethroid insecticide.

BACKGROUND: Triatomine control campaigns have traditionally consisted of spraying the inside of houses with pyrethroid insecticides. However, exposure to sublethal insecticide doses after the initial application is a common occurrence and may have phenotypic consequences for survivors. Here, using Triatoma infestans (the main vector of Chagas disease in the Southern Cone of South America) as a model species, we quantified the effects of exposure to a sublethal dose of pyrethroid insecticide on wing morphology. We tested if the treatment (i) induced a plastic effect (change in the character mean); (ii) altered environmental canalisation (higher individual variation within genotypes); (iii) altered genetic canalisation (higher variation among genotypes); and (iv) altered developmental stability (higher fluctuating asymmetry [FA]). METHODS: Each of 25 full-sib families known to be susceptible to pyrethroid insecticides were split in two groups: one to be treated with a sublethal dose of deltamethrin (insecticide-treated group) and the other to be treated with pure acetone (control group). Wings of the emerging adults were used in a landmark-based geometric morphometry analysis to extract size and shape measurements. Average differences among treatments were measured. Levels of variation among families, among individuals within families and among sides within individuals were computed and compared among treatments. RESULTS: Wing size and shape were affected by a sublethal dose of deltamethrin. The treated insects had larger wings and a more variable wing size and shape than control insects. For both wing size and shape, genetic variation was higher in treated individuals. Individual variations and variations in FA were also greater in deltamethrin-treated insects than in control ones for all full-sib families; however, the patterns of shape variation associated with genetic variation, individual variation and FA were different. CONCLUSIONS: Insects exposed to a sublethal dose of deltamethrin presented larger, less symmetrical and less canalised wings. The insecticide treatment jointly impaired developmental stability and genetic and environmental canalisation. The divergent patterns of shape variation suggest that the related developmental buffering processes differed at least partially. The morphological modifications induced by a single sublethal exposure to pyrethroids early in life may impinge on subsequent flight performance and consequently affect the dynamics of house invasion and reinfestation, and the effectiveness of triatomine control operations.

Constraints associated with captivity alter craniomandibular integration in wild boar.

The domestication process is associated with substantial phenotypic changes through time. However, although morphological integration between biological structures is purported to have a major influence on the evolution of new morphologies, little attention has been paid to the influence of domestication on the magnitude of integration. Here, we assessed the influence of constraints associated with captivity, considered as one of the crucial first steps in the domestication process, on the integration of cranial and mandibular structures. We investigated the craniomandibular integration in Western European Sus scrofa using three-dimensional (3D) landmark-based geometric morphometrics. Our results suggest that captivity is associated with a lower level of integration between the cranium and the mandible. Plastic responses to captivity can thus affect the magnitude of integration of key functional structures. These findings underline the critical need to develop integration studies in the context of animal domestication to better understand the processes accountable for the set-up of domestic phenotypes through time.

Convergence in sympatry: Evolution of blue-banded wing pattern in Morpho butterflies.

Species interactions such as mimicry can promote trait convergence but disentangling this effect from those of shared ecology, evolutionary history, and niche conservatism is often challenging. Here by focusing on wing colour pattern variation within and between three butterfly species living in sympatry in a large proportion of their range, we tested the effect of species interactions on trait diversification. These butterflies display a conspicuous iridescent blue coloration on the dorsal side of their wings and a cryptic brownish colour on the ventral side. Combined with an erratic and fast flight, these colour patterns increase the difficulty of capture by predators and contribute to the high escape abilities of these butterflies. We hypothesize that, beyond their direct contribution to predator escape, these wing patterns can be used as signals of escape abilities by predators, resulting in positive frequency-dependent selection favouring convergence in wing pattern in sympatry. To test this hypothesis, we quantified dorsal wing pattern variations of 723 butterflies from the three species sampled throughout their distribution, including sympatric and allopatric situations and compared the phenotypic distances between species, sex and localities. We detected a significant effect of localities on colour pattern, and higher inter-specific resemblance in sympatry as compared to allopatry, consistent with the hypothesis of local convergence of wing patterns. Our results provide support to the existence of escape mimicry in the wild and stress the importance of estimating trait variation within species to understand trait variation between species, and to a larger extent, trait diversification at the macro-evolutionary scale.

Correction to 'Limited thermal plasticity and geographical divergence in the ovipositor of Drosophila suzukii'.

[This corrects the article DOI: 10.1098/rsos.191577.].

Drosophila suzukii wing spot size is robust to developmental temperature.

Phenotypic plasticity is an important mechanism allowing adaptation to new environments and as such it has been suggested to facilitate biological invasions. Under this assumption, invasive populations are predicted to exhibit stronger plastic responses than native populations. Drosophila suzukii is an invasive species whose males harbor a spot on the wing tip. In this study, by manipulating developmental temperature, we compare the phenotypic plasticity of wing spot size of two invasive populations with that of a native population. We then compare the results with data obtained from wild-caught flies from different natural populations. While both wing size and spot size are plastic to temperature, no difference in plasticity was detected between native and invasive populations, rejecting the hypothesis of a role of the wing-spot plasticity in the invasion success. In contrast, we observed a remarkable stability in the spot-to-wing ratio across temperatures, as well as among geographic populations. This stability suggests either that the spot relative size is under stabilizing selection, or that its variation might be constrained by a tight developmental correlation between spot size and wing size. Our data show that this correlation was lost at high temperature, leading to an increased variation in the relative spot size, particularly marked in the two invasive populations. This suggests: (a) that D. suzukii's development is impaired by hot temperatures, in agreement with the cold-adapted status of this species; (b) that the spot size can be decoupled from wing size, rejecting the hypothesis of an absolute constraint and suggesting that the wing color pattern might be under stabilizing (sexual) selection; and (c) that such sexual selection might be relaxed in the invasive populations. Finally, a subtle but consistent directional asymmetry in spot size was detected in favor of the right side in all populations and temperatures, possibly indicative of a lateralized sexual behavior.

Hybridization and transgressive exploration of colour pattern and wing morphology in Heliconius butterflies.

Hybridization can generate novel phenotypes distinct from those of parental lineages, a phenomenon known as transgressive trait variation. Transgressive phenotypes might negatively or positively affect hybrid fitness, and increase available variation. Closely related species of Heliconius butterflies regularly produce hybrids in nature, and hybridization is thought to play a role in the diversification of novel wing colour patterns despite strong stabilizing selection due to interspecific mimicry. Here, we studied wing phenotypes in first- and second-generation hybrids produced by controlled crosses between either two co-mimetic species of Heliconius or between two nonmimetic species. We quantified wing size, shape and colour pattern variation and asked whether hybrids displayed transgressive wing phenotypes. Discrete traits underlain by major-effect loci, such as the presence or absence of colour patches, generate novel phenotypes. For quantitative traits, such as wing shape or subtle colour pattern characters, hybrids only exceed the parental range in specific dimensions of the morphological space. Overall, our study addresses some of the challenges in defining and measuring phenotypic transgression for multivariate traits and our data suggest that the extent to which transgressive trait variation in hybrids contributes to phenotypic diversity depends on the complexity and the genetic architecture of the traits.

Limited thermal plasticity and geographical divergence in the ovipositor of Drosophila suzukii.

Phenotypic plasticity has been repeatedly suggested to facilitate adaptation to new environmental conditions, as in invasions. Here, we investigate this possibility by focusing on the worldwide invasion of Drosophila suzukii: an invasive species that has rapidly colonized all continents over the last decade. This species is characterized by a highly developed ovipositor, allowing females to lay eggs through the skin of ripe fruits. Using a novel approach based on the combined use of scanning electron microscopy and photogrammetry, we quantified the ovipositor size and three-dimensional shape, contrasting invasive and native populations raised at three different developmental temperatures. We found a small but significant effect of temperature and geographical origin on the ovipositor shape, showing the occurrence of both geographical differentiation and plasticity to temperature. The shape reaction norms are in turn strikingly similar among populations, suggesting very little difference in shape plasticity among invasive and native populations, and therefore rejecting the hypothesis of a particular role for the plasticity of the ovipositor in the invasion success. Overall, the ovipositor shape seems to be a fairly robust trait, indicative of stabilizing selection. The large performance spectrum rather than the flexibility of the ovipositor would thus contribute to the success of D. suzukii worldwide invasion.

Effects of natural wing damage on flight performance in Morpho butterflies: what can it tell us about wing shape evolution?

Flying insects frequently experience wing damage during their life. Such irreversible alterations of wing shape affect flight performance and ultimately fitness. Insects have been shown to compensate for wing damage through various behavioural adjustments, but the importance of damage location over the wings has scarcely been studied. Using natural variation in wing damage, we tested how the loss of different wing parts affects flight performance. We quantified flight performance in two species of large butterflies, Morpho helenor and Morpho achilles, caught in the wild and displaying large variation in the extent and location of wing damage. We artificially generated more severe wing damage in our sample to contrast natural versus higher magnitude wing loss. Wing shape alteration across our sample was quantified using geometric morphometrics to test the effect of different damage distributions on flight performance. Our results show that impaired flight performance clearly depends on damage location over the wings, pointing to a relative importance of different wing parts for flight. A deteriorated forewing leading edge most critically affected flight performance, specifically decreasing flight speed and the proportion of gliding flight. In contrast, the most frequent natural damage, deteriorated wing margin, had no detectable effect on flight behaviour. Damage located on the hindwings - although having a limited effect on flight - was associated with reduced flight height, suggesting that the forewings and hindwings play different roles in butterfly flight. By contrasting harmless and deleterious consequences of various types of wing damage, our study highlights different selective regimes acting on morphological variations of butterfly wings.

Phenotypic plasticity, global change, and the speed of adaptive evolution.

The role phenotypic plasticity might play in adaptation to the ongoing climate changes is unclear. Plasticity allows for the production of a diversity of intra-generational responses, whose inter-generational evolutionary consequences are difficult to predict. In this article, we review theory and empirical studies addressing this question in insects by considering three scenarios. The first scenario corresponds to adaptive plasticity that should lead to slow or no evolution. The second scenario is the case of non-adaptive phenotypic plasticity to new environmental conditions that should lead either to extinction or, on the contrary, to rapid evolutionary change. The third scenario deals with how plasticity alters the variance selection acts upon. These scenarios are then discussed by highlighting examples of empirical studies on insects. We conclude that more studies are needed to better understand the relationship between phenotypic plasticity and evolutionary processes in insects.

Adaptive evolution of butterfly wing shape: from morphology to behaviour.

Butterflies display extreme variation in wing shape associated with tremendous ecological diversity. Disentangling the role of neutral versus adaptive processes in wing shape diversification remains a challenge for evolutionary biologists. Ascertaining how natural selection influences wing shape evolution requires both functional studies linking morphology to flight performance, and ecological investigations linking performance in the wild with fitness. However, direct links between morphological variation and fitness have rarely been established. The functional morphology of butterfly flight has been investigated but selective forces acting on flight behaviour and associated wing shape have received less attention. Here, we attempt to estimate the ecological relevance of morpho-functional links established through biomechanical studies in order to understand the evolution of butterfly wing morphology. We survey the evidence for natural and sexual selection driving wing shape evolution in butterflies, and discuss how our functional knowledge may allow identification of the selective forces involved, at both the macro- and micro-evolutionary scales. Our review shows that although correlations between wing shape variation and ecological factors have been established at the macro-evolutionary level, the underlying selective pressures often remain unclear. We identify the need to investigate flight behaviour in relevant ecological contexts to detect variation in fitness-related traits. Identifying the selective regime then should guide experimental studies towards the relevant estimates of flight performance. Habitat, predators and sex-specific behaviours are likely to be major selective forces acting on wing shape evolution in butterflies. Some striking cases of morphological divergence driven by contrasting ecology involve both wing and body morphology, indicating that their interactions should be included in future studies investigating co-evolution between morphology and flight behaviour.

Mouse Skull Mean Shape and Shape Robustness Rely on Different Genetic Architectures and Different Loci.

The genetic architecture of skull shape has been extensively studied in mice and the results suggest a highly polygenic and additive basis. In contrast few studies have explored the genetic basis of the skull variability. Canalization and developmental stability are the two components of phenotypic robustness. They have been proposed to be emergent properties of the genetic networks underlying the development of the trait itself, but this hypothesis has been rarely tested empirically. Here we use outbred mice to investigate the genetic architecture of canalization of the skull shape by implementing a genome-wide marginal epistatic test on 3D geometric morphometric data. The same data set had been used previously to explore the genetic architecture of the skull mean shape and its developmental stability. Here, we address two questions: (1) Are changes in mean shape and changes in shape variance associated with the same genomic regions? and (2) Do canalization and developmental stability rely on the same loci and genetic architecture and do they involve the same patterns of shape variation? We found that unlike skull mean shape, among-individual shape variance and fluctuating asymmetry (FA) show a total lack of additive effects. They are both associated with complex networks of epistatic interactions involving many genes (protein-coding and regulatory elements). Remarkably, none of the genomic loci affecting mean shape contribute these networks despite their enrichment for genes involved in craniofacial variation and diseases. We also found that the patterns of shape FA and individual variation are largely similar and rely on similar multilocus epistatic genetic networks, suggesting that the processes channeling variation within and among individuals are largely common. However, the loci involved in these two networks are completely different. This in turn underlines the difference in the origin of the variation at these two levels, and points at buffering processes that may be specific to each level.