Land-use-associated stressors interact to reduce bumblebee health at the individual and colony level

In agricultural landscapes, bees face a variety of stressors, including insecticides and poor-quality food. Although both stressors individually have been shown to affect bumblebee health negatively, few studies have focused on stressor interactions, a scenario expected in intensively used agricultural landscapes. Using the bumblebee Bombus terrestris, a key pollinator in agricultural landscapes, we conducted a fully factorial laboratory experiment starting at nest initiation. We assessed the effects of food quality and insecticides, alone and in interaction, on health traits at various levels, some of which have been rarely studied. Pollen with a diluted nutrient content (low quality) reduced ovary size and delayed colony development. Wing asymmetry, indicating developmental stress, was increased during insecticide exposure and interactions with poor food, whereas both stressors reduced body size. Both stressors and their interaction changed the workers' chemical profile and reduced worker interactions and the immune response. Our findings suggest that insecticides combined with nutritional stress reduce bumblebee health at the individual and colony levels, thus possibly affecting colony performance, such as development and reproduction, and the stability of plant-pollinator networks. The synergistic effects highlight the need of combining stressors in risk assessments and when studying the complex effects of anthropogenic stressors on health outcomes.

A landmark-free analysis of the pelvic girdle in Sulawesi ricefishes (Adrianichthyidae): How 2D and 3D geometric morphometrics can complement each other in the analysis of a complex structure

Geometric morphometrics (GM) enable the quantification of morphological variation on various scales. Recent technical advances allow analyzing complex three-dimensional shapes also in cases where landmark-based approaches are not appropriate. Pelvic girdle bones (basipterygia) of Sulawesi ricefishes are 3D structures that challenge traditional morphometrics. We hypothesize that the pelvic girdle of ricefishes experienced sex-biased selection pressures in species where females provide brood care by carrying fertilized eggs supported by elongated pelvic fins ("pelvic brooding"). We test this by comparing pelvic bone shapes of both sexes in species exhibiting pelvic brooding and the more common reproductive strategy "transfer brooding," by using landmark-free 2D and 3D GM, as well as qualitative shape descriptions. Both landmark-free approaches revealed significant interspecific pelvic bone variation in the lateral process, medial facing side of the pelvic bone, and overall external and internal wing shape. Within pelvic brooders, the three analyzed species are clearly distinct, while pelvic bones of the genus Adrianichthys are more similar to transfer brooding Oryzias. Female pelvic brooding Oryzias exhibit prominent, medially pointing tips extending from the internal wing and basipterygial plate that are reduced or absent in conspecific males, Adrianichthys and transfer brooding Oryzias, supporting our hypothesis that selection pressures affecting pelvic girdle shape are sex-biased in Sulawesi ricefishes. Furthermore, both sexes of pelvic brooding Oryzias have overall larger pelvic bones than other investigated ricefishes. Based on these differences, we characterized two reproductive strategy- and sex-dependent pelvic girdle types for Sulawesi ricefishes. Morphological differences between the investigated pelvic brooding genera Adrianichthys and Oryzias provide additional evidence for two independent origins of pelvic brooding. Overall, our findings add to a better understanding on traits related to pelvic brooding in ricefishes and provide a basis for upcoming studies on pelvic girdle function and morphology.

Signatures of increasing environmental stress in bumblebee wings over the past century: Insights from museum specimens

Determining when animal populations have experienced stress in the past is fundamental to understanding how risk factors drive contemporary and future species' responses to environmental change. For insects, quantifying stress and associating it with environmental factors has been challenging due to a paucity of time-series data and because detectable population-level responses can show varying lag effects. One solution is to leverage historic entomological specimens to detect morphological proxies of stress experienced at the time stressors emerged, allowing us to more accurately determine population responses. Here we studied specimens of four bumblebee species, an invaluable group of insect pollinators, from five museums collected across Britain over the 20th century. We calculated the degree of fluctuating asymmetry (FA; random deviations from bilateral symmetry) between the right and left forewings as a potential proxy of developmental stress. We: (a) investigated whether baseline FA levels vary between species, and how this compares between the first and second half of the century; (b) determined the extent of FA change over the century in the four bumblebee species, and whether this followed a linear or nonlinear trend; (c) tested which annual climatic conditions correlated with increased FA in bumblebees. Species differed in their baseline FA, with FA being higher in the two species that have recently expanded their ranges in Britain. Overall, FA significantly increased over the century but followed a nonlinear trend, with the increase starting c. 1925. We found relatively warm and wet years were associated with higher FA. Collectively our findings show that FA in bumblebees increased over the 20th century and under weather conditions that will likely increase in frequency with climate change. By plotting FA trends and quantifying the contribution of annual climate conditions on past populations, we provide an important step towards improving our understanding of how environmental factors could impact future populations of wild beneficial insects.

Flexible conservatism in the skull modularity of convergently evolved myrmecophagous placental mammals

Background

The skull of placental mammals constitutes one of the best studied systems for phenotypic modularity. Several studies have found strong evidence for the conserved presence of two- and six-module architectures, while the strength of trait correlations (integration) has been associated with major developmental processes such as somatic growth, muscle-bone interactions, and tooth eruption. Among placentals, ant- and termite-eating (myrmecophagy) represents an exemplar case of dietary convergence, accompanied by the selection of several cranial morphofunctional traits such as rostrum elongation, tooth loss, and mastication loss. Despite such drastic functional modifications, the covariance patterns of the skull of convergently evolved myrmecophagous placentals are yet to be studied in order to assess the potential consequences of this dietary shift on cranial modularity.

Results

Here, we performed a landmark-based morphometric analysis of cranial covariance patterns in 13 species of myrmecophagous placentals. Our analyses reveal that most myrmecophagous species present skulls divided into six to seven modules (depending on the confirmatory method used), with architectures similar to those of non-myrmecophagous placentals (therian six modules). Within-module integration is also similar to what was previously described for other placentals, suggesting that most covariance-generating processes are conserved across the clade. Nevertheless, we show that extreme rostrum elongation and tooth loss in myrmecophagid anteaters have resulted in a shift in intermodule correlations in the proximal region of the rostrum. Namely, the naso-frontal and maxillo-palatine regions are strongly correlated with the oro-nasal module, suggesting an integrated rostrum conserved from pre-natal developmental processes. In contrast, the similarly toothless pangolins show a weaker correlation between the anterior rostral modules, resembling the pattern of toothed placentals.

Conclusions

These results reveal that despite some integration shifts related to extreme functional and morphological features of myrmecophagous skulls, cranial modular architectures have conserved the typical mammalian scheme.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02030-9.

Presence of morphological integration and modularity of the forcipular apparatus in Lithobius melanops (Chilopoda: Lithobiomorpha: Lithobiidae)

The presence of morphological integration and modularity of the forcipular apparatus, despite its evolutionary significance, has not been analyzed in centipedes. This morphological structure has a crucial role in feeding and defense, thanks to its poisonous part (forcipules), which is important for catching the prey. The aims of our study were: i) to test the hypothesis of modularity of the forcipular apparatus in centipede Lithobius melanops; and ii) to investigate the influence of allometry on overall morphological integration in the aforementioned species using a geometric morphometric approach. The presence of fluctuating asymmetry was obtained by Procrustes ANOVA. Allometry was significant only for the symmetric component of the forcipular apparatus. The modularity hypothesis was not accepted, because the covariance coefficients for symmetric and asymmetric components were lower than 89.5% and 72.1% (respectively) of other RV coefficients obtained by a random contiguous partition of the forcipular apparatus. Results of the present study indicate that allometry does increase the level of morphological integration in the forcipular apparatus. According to our results, the forcipular coxosternite and forcipules could not be considered as separate modules; namely, they probably share similar developmental pathways and function in different forms of behavior and survival in L. melanops.

Untangling the relationship between developmental and evolutionary integration

Patterns of integration and modularity among organismal traits are prevalent across the tree of life, and at multiple scales of biological organization. Over the past several decades, researchers have studied these patterns at the developmental, and evolutionary levels. While their work has identified the potential drivers of these patterns at different scales, there appears to be a lack of consensus on the relationship between developmental and evolutionary integration. Here, we review and summarize key studies and build a framework to describe the conceptual relationship between these patterns across organismal scales and illustrate how, and why some of these studies may have yielded seemingly conflicting outcomes. We find that among studies that analyze patterns of integration and modularity using morphological data, the lack of consensus may stem in part from the difficulty of fully disentangling the developmental and functional causes of integration. Nonetheless, in some empirical systems, patterns of evolutionary modularity have been found to coincide with expectations based on developmental processes, suggesting that in some circumstances, developmental modularity may translate to evolutionary modularity. We also advance an extension to Hallgrímsson et al.'s palimpsest model to describe how patterns of trait modularity may shift across different evolutionary scales. Finally, we also propose some directions for future research which will hopefully be useful for investigators interested in these issues.

Effect of urbanization and its environmental stressors on the intraspecific variation of flight functional traits in two bumblebee species

The way urbanization shapes the intraspecific variation of pollinator functional traits is little understood. However, this topic is relevant for investigating ecosystem services and pollinator health. Here, we studied how urbanization affects the functional traits of workers in two bumblebee species (Bombus terrestris and B. pascuorum) sampled in 37 sites along a gradient of urbanization in North Italy (an area of 1800 km² including the metropolitan context of Milan and other surrounding capital districts). Namely, we investigated the effect of land use composition, configuration, air temperature, flower resource abundance, and air pollutants on the variation of traits related to flight performance and of stress during insect development (i.e., wing size, wing shape and size fluctuating asymmetry). The functional traits of the two bumblebees responded idiosyncratically to urbanization. Urban temperatures were associated with smaller wing sizes in B. pascuorum and with more accentuated fluctuating asymmetry of wing size in B. terrestris. Moreover, flower abundance correlated with bigger wings in B. terrestris and with less asymmetric wing size in B. pascuorum. Other traits did not vary significantly, and other urban variables played minor effects. These species-specific variation patterns highlight that environmental stressor linked to urbanization negatively impact the traits related to flight performance and development stability of these syntopic bumblebees, with possible consequences on the pollination service they provide.

Drosophila Wing Integration and Modularity: A Multi-Level Approach to Understand the History of Morphological Structures

Simple Summary

The diverse components of any morphological structure are integrated with respect to each other since they have developed, functioned, and evolved together, a phenomenon known as integration. However, this integration is not absolute but organized in units (i.e., modules) that are relatively independent while participating to generate a structure that acts as a functional whole. Even though most of the studies on modularity and integration have focused on variation among individuals within populations, there are more levels of variation that exhibit modularity and integration, deriving from distinct sources such as genetic variation, phenotypic plasticity, fluctuating asymmetry, evolutionary change, among others. Consequently, the present study focused on analysing the integration and modularity of the wing shape of some of the best-known model organisms, i.e., the genus Drosophila, at the static, developmental, and evolutionary levels to acquire a better insight about how modularity and integration act at different analytical levels. The strong integration and overall similarities observed in the variation pattern at multiple levels suggest a shared mechanism underlying the observed variation in Drosophila’s wing shape and added a new piece of evidence of stasis in the evolutionary history of Drosophila wing.

Abstract

Static, developmental, and evolutionary variation are different sources of morphological variation which can be quantified using morphometrics tools. In the present study we have carried out a comparative multiple level study of integration (i.e., static, developmental, and evolutionary) to acquire insight about the relationships that exist between different integration levels, as well as to better understand their involvement in the evolutionary processes related to the diversification of Drosophila’s wing shape. This approach was applied to analyse wing evolution in 59 species across the whole genus in a large dataset (~10,000 wings were studied). Static integration was analysed using principal component analysis, thus providing an integration measurement for overall wing shape. Developmental integration was studied between wing parts by using a partial least squares method between the anterior and posterior compartments of the wing. Evolutionary integration was analysed using independent contrasts. The present results show that all Drosophila species exhibit strong morphological integration at different levels. The strong integration and overall similarities observed at multiple integration levels suggest a shared mechanism underlying this variation, which could result as consequence of genetic drift acting on the wing shape of Drosophila.

Integration of Vestibular and Auditory Information in Ontogenesis

Background: We analyzed the hypothesis that the vestibular and auditory systems are integrative functions. Methods: The study involved 383 children (5.5 ± 2.4 years old). We assessed the conduct of auditory information by recording the auditory brain stem response (ABR), post-rotational nystagmus (PRN), and cervical vestibular evoked myogenic potentials (cVEMP), and calculated the integration of the parameters. All procedures were carried out using the JACOBI 4 software package. Results: We have found out that PRN, ABR, and cVEMP represent three different groups of integrative functions, each of which is conditioned by its own integrative mechanism. We have proven that PRN and ABR are associated with age, but no relationship was found between cVEMP and age. Conclusion: According to our data, the severity of ABR and PRN depended on age, while cVEMP was not associated with age. The functional immaturity of the child’s vestibular system, which probably arose in utero, often becomes apparent only at school when reading and writing must be mastered. These skills require maturity of both the vestibule ocular and vestibule spinal functions of the vestibular system.

Insect wing 3D printing

Insects have acquired various types of wings over their course of evolution and have become the most successful terrestrial animals. Consequently, the essence of their excellent environmental adaptability and locomotive ability should be clarified; a simple and versatile method to artificially reproduce the complex structure and various functions of these innumerable types of wings is necessary. This study presents a simple integral forming method for an insect-wing-type composite structure by 3D printing wing frames directly onto thin films. The artificial venation generation algorithm based on the centroidal Voronoi diagram, which can be observed in the wings of dragonflies, was used to design the complex mechanical properties of artificial wings. Furthermore, we implemented two representative functions found in actual insect wings: folding and coupling. The proposed crease pattern design software developed based on a beetle hindwing enables the 3D printing of foldable wings of any shape. In coupling-type wings, the forewing and hindwing are connected to form a single large wing during flight; these wings can be stored compactly by disconnecting and stacking them like cicada wings.

Morphological integration of the head capsule in the millipede Megaphyllum unilineatum (C. L. Koch, 1838) (Diplopoda: Julida): can different modules be recognized?

Covariation of multiple morphological traits and modularity have been widely studied in the field of evolutionary developmental biology. Subunits of a morphological structure can evolve separately from each other in a modular fashion. The aims of our study therefore were: i) to test the hypothesis of modularity in the dorsal part of the head capsule and the gnathochilarium separately during late postembryogenesis in the julidan millipede Megaphyllum unilineatum (C. L. Koch, 1838) using geometric morphometrics; and ii) to investigate the influence of allometry on overall morphological integration in the dorsal part of the head capsule and the gnathochilarium in the mentioned species. Individuals from different ontogenetic stadia (stadium VI - stadium XI) were included in the analyses. Significant influence of fluctuating asymmetry on the dorsal part of the head capsule shape was detected by Procrustes ANOVA. Regressions were significant for the symmetric component of both analysed morphological traits, while non-significant regression was detected for the asymmetric component of the head capsule's dorsal part. Hypotheses of modularity for the dorsal part of the head capsule and the gnathochilarium are rejected because our results indicate that a small proportion of alternate partitions has higher covariation between subsets of structure than between the hypothesized modules. Contrary to our expectations, results of the present study show that allometry does not increase the level of morphological integration in the dorsal part of the head capsule and the gnathochilarium in M. unilineatum. Based on the obtained results, we conclude that the dorsal part of the head capsule and the gnathochilarium are not composed of independent modules and that in the case of the capsule's dorsal part, developmental processes affect morphological integration in different ways at different levels of shape variation.

Including Distorted Specimens in Allometric Studies: Linear Mixed Models Account for Deformation

Allometry—patterns of relative change in body parts—is a staple for examining how clades exhibit scaling patterns representative of evolutionary constraint on phenotype, or quantifying patterns of ontogenetic growth within a species. Reconstructing allometries from ontogenetic series is one of the few methods available to reconstruct growth in fossil specimens. However, many fossil specimens are deformed (twisted, flattened, and displaced bones) during fossilization, changing their original morphology in unpredictable and sometimes undecipherable ways. To mitigate against post burial changes, paleontologists typically remove clearly distorted measurements from analyses. However, this can potentially remove evidence of individual variation and limits the number of samples amenable to study, which can negatively impact allometric reconstructions. Ordinary least squares (OLS) regression and major axis regression are common methods for estimating allometry, but they assume constant levels of residual variation across specimens, which is unlikely to be true when including both distorted and undistorted specimens. Alternatively, a generalized linear mixed model (GLMM) can attribute additional variation in a model (e.g., fixed or random effects). We performed a simulation study based on an empirical analysis of the extinct cynodont, Exaeretodon argentinus, to test the efficacy of a GLMM on allometric data. We found that GLMMs estimate the allometry using a full dataset better than simply using only non-distorted data. We apply our approach on two empirical datasets, cranial measurements of actual specimens of E. argentinus (n = 16) and femoral measurements of the dinosaur Tawa hallae (n = 26). Taken together, our study suggests that a GLMM is better able to reconstruct patterns of allometry over an OLS in datasets comprised of extinct forms and should be standard protocol for anyone using distorted specimens.

Exploring directional and fluctuating asymmetry in the human palate during growth

Objectives

Palate morphology is constantly changing throughout an individual's lifespan, yet its asymmetry during growth is still little understood. In this research, we focus on the study of palate morphology by using 3D geometric morphometric approaches to observe changes at different stages of life, and to quantify the impact of directional and fluctuating asymmetry on different areas at different growth stages.

Materials and Methods

The sample consists of 183 individuals (1–72 years) from two identified human skeletal collections of 19th and early 20th Century Italian contexts. A 3D‐template of 41 (semi)landmarks was applied on digital palate models to observe morphological variation during growth.

Results

Asymmetrical components of the morphological structure appears multidirectional on the entire palate surface in individuals <2 years old and become oriented (opposite bilateral direction) between 2 and 6 years of age. Specifically, directional asymmetry differentially impacts palate morphology at different stages of growth. Both the anterior and posterior palate are affected by mild alterations in the first year of life, while between 2 and 6 years asymmetry is segregated in the anterior area, and moderate asymmetry affects the entire palatal surface up to 12 years of age. Our results show that stability of the masticatory system seems to be reached around 13–35 years first by females and then males. From 36 years on both sexes show similar asymmetry on the anterior area. Regarding fluctuating asymmetry, inter‐individual variability is mostly visible up to 12 years of age, after which only directional trends can be clearly observed at a group level.

Discussion

Morphological structure appears instable during the first year of life and acquires an opposite asymmetric bilateral direction between 2 and 6 years of age. This condition has been also documented in adults; when paired with vertical alteration, anterior/posterior asymmetry seems to characterize palate morphology, which is probably due to mechanical factors during the lifespan. Fluctuating asymmetry is predominant in the first period of life due to a plausible relationship with the strength of morphological instability of the masticatory system. Directional asymmetry, on the other hand, shows that the patterning of group‐level morphological change might be explained as a functional response to differential inputs (physiological forces, nutritive and non‐nutritive habits, para‐masticatory activity as well as the development of speech) in different growth stages. This research has implications with respect to medical and evolutionary fields. In medicine, palate morphology should be considered when planning orthodontic and surgical procedures as it could affect the outcome. As far as an evolutionary perspective is concerned the dominance of directional asymmetries in the masticatory system could provide information on dietary and cultural habits as well as pathological conditions in our ancestors.

Sexual Dimorphism and Morphological Modularity in Acanthoscelides obtectus (Say, 1831) (Coleoptera: Chrysomelidae): A Geometric Morphometric Approach

Simple Summary

The seed beetle Acanthoscelides obtectus used in this study is a worldwide pest species that inhabits storage facilities and fields of beans. Knowing that sexual dimorphism is very common among insects, we investigated the level of morphological differences between the sexes. Such an approach allowed us to look into the modular organization of this organism. As expected, the females were larger than the males. The observed two modular organization (thorax and abdomen) was sex specific, indicating that reproductive function has the central role in forming the patterns of modularity. It seems that natural selection is driving force for females, while males are influenced more by sexual selection.

Abstract

Sexual dimorphism and specific patterns of development contribute in a great manner to the direction and degree of the sexual differences in body size and shape in many insects. Using a landmark-based geometric morpohometrics approach, we investigated sex-specific morphological size and shape variation in the seed beetle, Acanthoscelides obtectus. We also tested the functional hypothesis of the two morphological modules—thorax and abdomen in both sexes. Female-biased sexual dimorphism in size was shown, while differences in shape were reflected in the wider thorax and abdomen and shorter abdomen in females in comparison to males. The functional hypothesis of a two-module body was confirmed only in females before correction for size, and in both sexes after the allometry correction. Our results indicate that reproductive function has the central role in forming the patterns of modularity. We hypothesize that high morphological integration of the abdomen in females results from intense stabilizing selection, while the more relaxed integration in males is driven by the higher intensity of sexual selection.

Morphological integration of the canine region within the hominine alveolar arch

The early hominin record is characterized by numerous shifts in dental proportions (e.g., canine reduction and megadontia) linked to changes in diet and social behavior. Recent studies suggest that hominins exhibit a reduction in the magnitude of covariation between the anterior and posterior dental components compared with other extant great apes. They point toward, but do not directly test, the relative independence of canine morphology within the hominin alveolar arch. This study focuses specifically on the how the canine region covaries with other regions of the dental arch because the canine region has drastically reduced in size and changed in shape across human evolution. We examine extant primate species most commonly used as a comparative framework for fossil hominin morphology: Gorilla gorilla (n = 27), Pan troglodytes (n = 27), and Homo sapiens (n = 30). We used geometric morphometric methods to test for size and shape covariation between the canine region with other dental regions. We also examined the influence of sexual dimorphism and allometry on intraspecific and interspecific patterns of covariation. The analysis of size and shape covariation between the mandibular canine and other individual tooth regions elucidated complex, species-specific, and sex-specific morphological relationships in the mandibular alveolar arch. There was little evidence to support different patterns of morphological integration between humans on the one hand and nonhuman apes on the other. Canine region morphology was relatively independent from other dental regions across species based on shape and did not significantly covary more with either the incisor or postcanine region in any species. The size correlations between the canine and other dental regions were moderate to high. The species-specific results of this study question the ability to make a priori assumptions about morphological integration in the extant hominin mandibular alveolar arch and its application to the fossil record.

Designing future crops: challenges and strategies for sustainable agriculture

Crop production is facing unprecedented challenges. Despite the fact that the food supply has significantly increased over the past half-century, ~8.9 and 14.3% people are still suffering from hunger and malnutrition, respectively. Agricultural environments are continuously threatened by a booming world population, a shortage of arable land, and rapid changes in climate. To ensure food and ecosystem security, there is a need to design future crops for sustainable agriculture development by maximizing net production and minimalizing undesirable effects on the environment. The future crops design projects, recently launched by the National Natural Science Foundation of China and Chinese Academy of Sciences (CAS), aim to develop a roadmap for rapid design of customized future crops using cutting-edge technologies in the Breeding 4.0 era. In this perspective, we first introduce the background and missions of these projects. We then outline strategies to design future crops, such as improvement of current well-cultivated crops, de novo domestication of wild species and redomestication of current cultivated crops. We further discuss how these ambitious goals can be achieved by the recent development of new integrative omics tools, advanced genome-editing tools and synthetic biology approaches. Finally, we summarize related opportunities and challenges in these projects.

The weakest link: Haploid honey bees are more susceptible to neonicotinoid insecticides

Neonicotinoid insecticides are currently of major concern for the health of wild and managed insects that provide key ecosystem services like pollination. Even though sublethal effects of neonicotinoids are well known, there is surprisingly little information on how they possibly impact developmental stability, and to what extent genetics are involved. This holds especially true for haploid individuals because they are hemizygous at detoxification loci and may be more susceptible. Here we take advantage of haplodiploidy in Western honey bees, Apis mellifera, to show for the first time that neonicotinoids affect developmental stability in diploid females (workers), and that haploid males (drones) are even more susceptible. Phenotypic fore wing venation abnormalities and fluctuating wing asymmetry, as measures of developmental instability, were significantly increased under field-realistic neonicotinoid-exposure of colonies. The higher susceptibility of haploid drones suggests that heterozygosity can play a key role in the ability to buffer the sublethal effects of neonicotinoids. Aiming to improve conservation efforts, our findings highlight the urgent need to better understand the role that genetics plays at enabling non-target organisms to cope with insecticide exposure.

Morphological consequences of artificial cranial deformation: Modularity and integration

The cranium is an anatomically complex structure. One source of its complexity is due to its modular organization. Cranial modules are distinct and partially independent units that interact substantially during ontogeny thus generating morphological integration. Artificial Cranial Deformation (ACD) occurs when the human skull is intentionally deformed, through the use of different deforming devices applied to the head while it is developing. Hence, ACD provides an interesting example to assess the degree to which biomechanical perturbations of the developing neurocranium impact on the degree of morphological integration in the skull as a whole. The main objective of this study was to assess how ACD affects the morphological integration of the skull. This was accomplished by comparing a sample of non-deformed crania and two sets of deformed crania (i.e. antero-posterior and oblique). Both developmental and static modularity and integration were assessed through Generalized Procrustes Analysis by considering the symmetric and asymmetric components of variation in adults, using 3D landmark coordinates as raw data. The presence of two developmental modules (i.e. viscero and neurocranium) in the skull was tested. Then, in order to understand how ACD affects morphological integration, the covariation pattern between the neuro and viscerocranium was examined in antero-posterior, oblique and non-deformed cranial categories using Partial Least-Squares. The main objective of this study was to assess how ACD affects the morphological integration of the skull. This was accomplished by comparing a sample of deformed (i.e. antero-posterior and oblique) and non-deformed crania. Hence, differences in integration patterns were compared between groups. The obtained results support the modular organization of the human skull in the two analyzed modules. The integration analyses show that the oblique ACD style differentially affects the static morphological integration of the skull by increasing the covariance between neuro and viscerocranium in a more constrained way than in antero-posterior and non-deformed skulls. In addition, the antero-posterior ACD style seems to affect the developmental integration of the skull by directing the covariation pattern in a more defined manner as compared to the other cranial categories.

Frequent nonrandom shifts in the temporal sequence of developmental landmark events during teleost evolutionary diversification

Morphological transformations can be generated by evolutionary changes in the sequence of developmental events. In this study, we examined the evolutionary dynamics of the developmental sequence on a macroevolutionary scale in teleosts. Using the information from previous reports describing the development of 31 species, we extracted the developmental sequences of 19 landmark events involving the formation of phylogenetically conserved body parts; we then inferred ancestral developmental sequences by two different parsimony-based methods-event-pairing and continuous analysis. The phylogenetic comparisons of these sequences revealed event-dependent heterogeneity in the frequency of sequence changes. Most of the sequence changes occurred as exchanges of temporally neighboring events. These heterochronic changes in developmental sequences accumulated along evolutionary time, but the precise distribution of the changes over the teleostean phylogeny remains unclear due to technical limitations.

Urbanisation and wing asymmetry in the western honey bee (Apis mellifera, Linnaeus 1758) at multiple scales

Changes in the mean and variance of phenotypic traits like wing and head morphology are frequently used as indicators of environmental stress experienced during development and may serve as a convenient index of urbanization exposure. To test this claim, we collected adult western honey bee (Apis mellifera Linnaeus 1758, Hymenoptera, Apidae) workers from colonies located across an urbanization gradient, and quantified associations between the symmetries of both wing size and wing shape, and several landscape traits associated with urbanization. Landscape traits were assessed at two spatial scales (three km and 500 m) and included vegetation and anthropogenic land cover, total road length, road proximity and, population and dwelling density. We then used geometric morphometric techniques to determine two wing asymmetry scores-centroid size, a measure of wing size asymmetry and Procrustes distance, a measure of wing shape asymmetry. We found colony dependent differences in both wing size and shape asymmetry. Additionally, we found a negative association between wing shape asymmetry and road proximity at the three km buffer, and associations between wing shape asymmetry and road proximity, anthropogenic land cover and vegetation cover at the 500 m buffer. Whilst we were unable to account for additional variables that may influence asymmetry including temperature, pesticide presence, and parasitism our results demonstrate the potential usefulness of wing shape asymmetry for assessing the impact of certain landscape traits associated with urbanization. Furthermore, they highlight important spatial scale considerations that warrant investigation in future phenotypic studies assessing urbanization impact.

Stressful conditions reveal decrease in size, modification of shape but relatively stable asymmetry in bumblebee wings

Human activities can generate a wide variety of direct and indirect effects on animals, which can manifest as environmental and genetic stressors. Several phenotypic markers have been proposed as indicators of these stressful conditions but have displayed contrasting results, depending, among others, on the phenotypic trait measured. Knowing the worldwide decline of multiple bumblebee species, it is important to understand these stressors and link them with the drivers of decline. We assessed the impact of several stressors (i.e. natural toxin-, parasite-, thermic- and inbreeding- stress) on both wing shape and size and their variability as well as their directional and fluctuating asymmetries. The total data set includes 650 individuals of Bombus terrestris (Hymenoptera: Apidae). Overall wing size and shape were affected by all the tested stressors. Except for the sinigrin (e.g. glucosinolate) stress, each stress implies a decrease of wing size. Size variance was affected by several stressors, contrary to shape variance that was affected by none of them. Although wing size directional and fluctuating asymmetries were significantly affected by sinigrin, parasites and high temperatures, neither directional nor fluctuating shape asymmetry was significantly affected by any tested stressor. Parasites and high temperatures led to the strongest phenotype modifications. Overall size and shape were the most sensitive morphological traits, which contrasts with the common view that fluctuating asymmetry is the major phenotypic marker of stress.

A simple developmental model recapitulates complex insect wing venation patterns

Insect wings are typically supported by thickened struts called veins. These veins form diverse geometric patterns across insects. For many insect species, even the left and right wings from the same individual have veins with unique topological arrangements, and little is known about how these patterns form. We present a large-scale quantitative study of the fingerprint-like "secondary veins." We compile a dataset of wings from 232 species and 17 families from the order Odonata (dragonflies and damselflies), a group with particularly elaborate vein patterns. We characterize the geometric arrangements of veins and develop a simple model of secondary vein patterning. We show that our model is capable of recapitulating the vein geometries of species from other, distantly related winged insect clades.

Analysis of modularity and integration suggests evolution of dragonfly wing venation mainly in response to functional demands

Insect wings show a high variability in wing venation. Selection for function, developmental pathways and phylogeny likely influenced wing vein diversification, however, quantitative data to estimate these influences and their interplay are missing. Here, it is tested how dragonfly wing vein configuration is influenced by functional demands, development, phylogeny and allometry using the concepts of modularity and integration. In an evolutionary context, modules are sets of characters that evolve in relative independence to other characters, while integration refers to a high degree of association between subparts of a structure. Results show allometric and phylogenetic signal in the wing shape variation, however, patterns of integration and modularity are not influenced by these two factors. Overall, dragonfly wings are highly integrated structures with almost no modular signal. Configuration changes in one wing vein or wing area thus influence wing shape as a whole. Moreover, the fore- and hindwings correlate with each other in their evolutionary shape variation supporting biomechanical data of wing interdependence. Despite the overall high degree of evolutionary integration, functional hypotheses of modularity could be confirmed for two wing areas, the arculus-triangle complex at the base of the wing which is responsible for passive wing folding especially during flapping flight and the location of the pterostigma, a coloured wing cell which is more heavy that other wing cells and passively regulates wing pitch as well as critical flight speeds during gliding. Although evolving as distinct modules, these specific vein regions also show high integration and evolve at the same rates like the whole wing which suggests an influence of these structures on the shape evolution of the rest of the wing. Their biomechanical role as passive regulators of wing corrugation and wing pitch suggests that these structures decisively influenced the evolution of advanced modern flight styles and explains their retention once they had evolved early within the lineage Odonatoptera.

Developmental pathways inferred from modularity, morphological integration and fluctuating asymmetry patterns in the human face

Facial asymmetries are usually measured and interpreted as proxies to developmental noise. However, analyses focused on its developmental and genetic architecture are scarce. To advance on this topic, studies based on a comprehensive and simultaneous analysis of modularity, morphological integration and facial asymmetries including both phenotypic and genomic information are needed. Here we explore several modularity hypotheses on a sample of Latin American mestizos, in order to test if modularity and integration patterns differ across several genomic ancestry backgrounds. To do so, 4104 individuals were analyzed using 3D photogrammetry reconstructions and a set of 34 facial landmarks placed on each individual. We found a pattern of modularity and integration that is conserved across sub-samples differing in their genomic ancestry background. Specifically, a signal of modularity based on functional demands and organization of the face is regularly observed across the whole sample. Our results shed more light on previous evidence obtained from Genome Wide Association Studies performed on the same samples, indicating the action of different genomic regions contributing to the expression of the nose and mouth facial phenotypes. Our results also indicate that large samples including phenotypic and genomic metadata enable a better understanding of the developmental and genetic architecture of craniofacial phenotypes.

Wing morphology variations in a natural population of Phlebotomus tobbi Adler and Theodor 1930

Cutaneous leishmaniasis (CL) is highly endemic in the Cukurova region, located on the crossroads of main refugee routes from the Middle East to Europe on the eastern Mediterranean part of Turkey. Our purpose was to investigate the phenotypic variation of Phlebotomus tobbi, the known vector of CL in the region, during one active season. Sand flies and microclimatic data were collected monthly from May to October, 2011, from five locations in six villages in the study area. A geometric morphometric approach was used to investigate wing morphology. Shape analyses revealed that males collected in May and June comprised one group, while specimens collected in August, September, and October formed a second group. Specimens from July were found to be distributed within these two groups. A similar distribution pattern was observed for females, but specimens from October were represented as the third district group. Significant size variation was detected for both sexes between months. Wing size and temperature were negatively correlated for females, but there was no temperature effect for males. Wing size of both sexes was increased in correlation to increasing relative humidity. Males were found to have smaller wings with increasing population density.

Most Colorful Example of Genetic Assimilation? Exploring the Evolutionary Destiny of Recurrent Phenotypic Accommodation

Evolution of adaptation requires both generation of novel phenotypic variation and retention of a locally beneficial subset of this variation. Such retention can be facilitated by genetic assimilation, the accumulation of genetic and molecular mechanisms that stabilize induced phenotypes and assume progressively greater control over their reliable production. A particularly strong inference into genetic assimilation as an evolutionary process requires a system where it is possible to directly evaluate the extent to which an induced phenotype is progressively incorporated into preexisting developmental pathways. Evolution of diet-dependent pigmentation in birds-where external carotenoids are coopted into internal metabolism to a variable degree before being integrated with a feather's developmental processes-provides such an opportunity. Here we combine a metabolic network view of carotenoid evolution with detailed empirical study of feather modifications to show that the effect of physical properties of carotenoids on feather structure depends on their metabolic modification, their environmental recurrence, and biochemical redundancy, as predicted by the genetic assimilation hypothesis. Metabolized carotenoids caused less stochastic variation in feather structure and were more closely integrated with feather growth than were dietary carotenoids of the same molecular weight. These patterns were driven by the recurrence of organism-carotenoid associations: commonly used dietary carotenoids and biochemically redundant derived carotenoids caused less stochastic variation in feather structure than did rarely used or biochemically unique compounds. We discuss implications of genetic assimilation processes for the evolutionary diversification of diet-dependent animal coloration.

Methods for Studying Morphological Integration and Modularity

Morphological integration and modularity are closely related concepts about how different traits of an organism are correlated. Integration is the overall pattern of intercorrelation; modularity is the partitioning of integration into evolutionarily or developmentally independent blocks of traits. Modularity and integration are usually studied using quantitative phenotypic data, which can be obtained either from extant or fossil organisms. Many methods are now available to study integration and modularity, all of which involve the analysis of patterns found in trait correlation or covariance matrices. We review matrix correlation, random skewers, fluctuating asymmetry, cluster analysis, Euclidean distance matrix analysis (EDMA), graphical modelling, two-block partial least squares, RV coefficients, and theoretical matrix modelling and discuss their similarities and differences. We also review different coefficients that are used to measure correlations. We apply all the methods to cranial landmark data from and ontogenetic series of Japanese macaques,Macaca fuscatato illustrate the methods and their individual strengths and weaknesses. We conclude that the exploratory approaches (cluster analyses of various sorts) were less informative and less consistent with one another than were the results of model testing or comparative approaches. Nevertheless, we found that competing models of modularity and integration are often similar enough that they are not statistically distinguishable; we expect, therefore, that several models will often be significantly correlated with observed data.

Network architecture associated with the highly specialized hindlimb of frogs

Network analyses have been increasingly used in the context of comparative vertebrate morphology. The structural units of the vertebrate body are treated as discrete elements (nodes) of a network, whose interactions at their physical contacts (links) determine the phenotypic modules. Here, we use the network approach to study the organization of the locomotor system underlying the hindlimb of frogs. Nodes correspond to fibrous knots, skeletal and muscular units. Edges encode the ligamentous and monoaxial tendinous connections in addition to joints. Our main hypotheses are that: (1) the higher centrality scores (measured as betweenness) are recorded for fibrous elements belonging to the connective system, (2) the organization of the musculoskeletal network belongs to a non-trivial modular architecture and (3) the modules in the hindlimb reflect functional and/or developmental constraints. We confirm all our hypotheses except for the first one, since bones overpass the fibrous knots in terms of centrality. Functionally, there is a correlation between the proximal-to-distal succession of modules and the progressive recruitment of elements involved with the motion of joints during jumping. From a developmental perspective, there is a correspondence between the order of the betweenness scores and the ontogenetic chronology of hindlimbs in tetrapods. Modular architecture seems to be a successful organization, providing of the building blocks on which evolution forges the many different functional specializations that organisms exploit.

Size, shape, and form: concepts of allometry in geometric morphometrics

Allometry refers to the size-related changes of morphological traits and remains an essential concept for the study of evolution and development. This review is the first systematic comparison of allometric methods in the context of geometric morphometrics that considers the structure of morphological spaces and their implications for characterizing allometry and performing size correction. The distinction of two main schools of thought is useful for understanding the differences and relationships between alternative methods for studying allometry. The Gould-Mosimann school defines allometry as the covariation of shape with size. This concept of allometry is implemented in geometric morphometrics through the multivariate regression of shape variables on a measure of size. In the Huxley-Jolicoeur school, allometry is the covariation among morphological features that all contain size information. In this framework, allometric trajectories are characterized by the first principal component, which is a line of best fit to the data points. In geometric morphometrics, this concept is implemented in analyses using either Procrustes form space or conformation space (the latter also known as size-and-shape space). Whereas these spaces differ substantially in their global structure, there are also close connections in their localized geometry. For the model of small isotropic variation of landmark positions, they are equivalent up to scaling. The methods differ in their emphasis and thus provide investigators with flexible tools to address specific questions concerning evolution and development, but all frameworks are logically compatible with each other and therefore unlikely to yield contradictory results.

Evolutionary modularity and morphological integration in the haptoral anchor structures of Ligophorus spp. (Monogenea: Dactylogyridae)

An important question in the study of phenotypic evolution is whether characters are independent of each other or behave and evolve as integrated modules. Morphological integration and modularity provide a powerful framework for the analysis of the evolution of morphological traits. We used geometric morphometrics and phylogenetically independent contrasts (PIC) to test four different modularity hypotheses in the haptoral anchors of 14 monogenean species of Ligophorus. Integration between the modular units identified was further evaluated with two-block partial least squares analysis. Roots and points represented two modules in the dorsal and ventral anchors, but modularity was not statistically supported when parasite phylogeny was accounted for, which may indicate convergent evolution related to host characteristics and gill morphology. In contrast, PIC revealed medial and lateral modules in ventral anchors only. Moreover, we found evidence for ventral and dorsal anchor pairs forming two modules, supporting the notion that they play different functional roles. Integration between all identified modules was strong. We conclude that there is modular structure in the anchors of Ligophorus spp., accounted by adaptive and phylogenetic factors acting at different levels, and ventral and dorsal anchors evolve as integrated modules with specific roles in attachment.

Size, shape, and form: concepts of allometry in geometric morphometrics

Allometry refers to the size-related changes of morphological traits and remains an essential concept for the study of evolution and development. This review is the first systematic comparison of allometric methods in the context of geometric morphometrics that considers the structure of morphological spaces and their implications for characterizing allometry and performing size correction. The distinction of two main schools of thought is useful for understanding the differences and relationships between alternative methods for studying allometry. The Gould–Mosimann school defines allometry as the covariation of shape with size. This concept of allometry is implemented in geometric morphometrics through the multivariate regression of shape variables on a measure of size. In the Huxley–Jolicoeur school, allometry is the covariation among morphological features that all contain size information. In this framework, allometric trajectories are characterized by the first principal component, which is a line of best fit to the data points. In geometric morphometrics, this concept is implemented in analyses using either Procrustes form space or conformation space (the latter also known as size-and-shape space). Whereas these spaces differ substantially in their global structure, there are also close connections in their localized geometry. For the model of small isotropic variation of landmark positions, they are equivalent up to scaling. The methods differ in their emphasis and thus provide investigators with flexible tools to address specific questions concerning evolution and development, but all frameworks are logically compatible with each other and therefore unlikely to yield contradictory results.

Phenotypic Buffering in a Monogenean: Canalization and Developmental Stability in Shape and Size of the Haptoral Anchors of Ligophorus cephali (Monogenea: Dactylogyridae)

Phenotypic variation results from the balance between sources of variation and counteracting regulatory mechanisms. Canalization and developmental stability are two such mechanisms, acting at two different levels of regulation. The issue of whether or not they act concurrently as a common developmental buffering capacity has been subject to debate. We used geometric morphometrics to quantify the mechanisms that guarantee phenotypic constancy in the haptoral anchors of Ligophorus cephali. Canalization and developmental stability were appraised by estimating inter- and intra-individual variation, respectively, in size and shape of dorsal and ventral anchors. The latter variation was estimated as fluctuating asymmetry (FA) between anchor pairs. The general-buffering-capacity hypothesis was tested by two different methods based on correlations and Principal Components Analyses of the different components of size and shape variation. Evidence for FA in the dorsal and ventral anchors in both shape and size was found. Our analyses supported the hypothesis of a general developmental buffering capacity. The evidence was more compelling for shape than for size and, particularly, for the ventral anchors than for the dorsal ones. These results are in line with previous studies of dactylogyrids suggesting that ventral anchors secure a firmer, more permanent attachment, whereas dorsal anchors are more mobile. Because fixation to the host is crucial for survival in ectoparasites, we suggest that homeostatic development of the ventral anchors has been promoted to ensure the morphological constancy required for efficient attachment. Geometric morphometrics can be readily applied to other host-monogenean models, affording not only to disentangle the effects of canalization and developmental stability, as shown herein, but to further partition the environmental and genetic components of the former.

Morphometrics, 3D Imaging, and Craniofacial Development

Recent studies have shown how volumetric imaging and morphometrics can add significantly to our understanding of morphogenesis, the developmental basis for variation, and the etiology of structural birth defects. On the other hand, the complex questions and diverse imaging data in developmental biology present morphometrics with more complex challenges than applications in virtually any other field. Meeting these challenges is necessary in order to understand the mechanistic basis for variation in complex morphologies. This chapter reviews the methods and theory that enable the application of modern landmark-based morphometrics to developmental biology and craniofacial development, in particular. We discuss the theoretical foundations of morphometrics as applied to development and review the basic approaches to the quantification of morphology. Focusing on geometric morphometrics, we discuss the principal statistical methods for quantifying and comparing morphological variation and covariation structure within and among groups. Finally, we discuss the future directions for morphometrics in developmental biology that will be required for approaches that enable quantitative integration across the genotype-phenotype map.

Environmentally induced changes in correlated responses to selection reveal variable pleiotropy across a complex genetic network

Selection in novel environments can lead to a coordinated evolutionary response across a suite of characters. Environmental conditions can also potentially induce changes in the genetic architecture of complex traits, which in turn could alter the pattern of the multivariate response to selection. We describe a factorial selection experiment using the nematode Caenorhabditis remanei in which two different stress-related phenotypes (heat and oxidative stress resistance) were selected under three different environmental conditions. The pattern of covariation in the evolutionary response between phenotypes or across environments differed depending on the environment in which selection occurred, including asymmetrical responses to selection in some cases. These results indicate that variation in pleiotropy across the stress response network is highly sensitive to the external environment. Our findings highlight the complexity of the interaction between genes and environment that influences the ability of organisms to acclimate to novel environments. They also make clear the need to identify the underlying genetic basis of genetic correlations in order understand how patterns of pleiotropy are distributed across complex genetic networks.

Studying morphological integration and modularity at multiple levels: concepts and analysis

Although most studies on integration and modularity have focused on variation among individuals within populations or species, this is not the only level of variation for which integration and modularity exist. Multiple levels of biological variation originate from distinct sources: genetic variation, phenotypic plasticity resulting from environmental heterogeneity, fluctuating asymmetry from random developmental variation and, at the interpopulation or interspecific levels, evolutionary change. The processes that produce variation at all these levels can impart integration or modularity on the covariance structure among morphological traits. In turn, studies of the patterns of integration and modularity can inform about the underlying processes. In particular, the methods of geometric morphometrics offer many advantages for such studies because they can characterize the patterns of morphological variation in great detail and maintain the anatomical context of the structures under study. This paper reviews biological concepts and analytical methods for characterizing patterns of variation and for comparing across levels. Because research comparing patterns across level has only just begun, there are relatively few results, generalizations are difficult and many biological and statistical questions remain unanswered. Nevertheless, it is clear that research using this approach can take advantage of an abundance of new possibilities that are so far largely unexplored.

The macroevolutionary consequences of phenotypic integration: from development to deep time

Phenotypic integration is a pervasive characteristic of organisms. Numerous analyses have demonstrated that patterns of phenotypic integration are conserved across large clades, but that significant variation also exists. For example, heterochronic shifts related to different mammalian reproductive strategies are reflected in postcranial skeletal integration and in coordination of bone ossification. Phenotypic integration and modularity have been hypothesized to shape morphological evolution, and we extended simulations to confirm that trait integration can influence both the trajectory and magnitude of response to selection. We further demonstrate that phenotypic integration can produce both more and less disparate organisms than would be expected under random walk models by repartitioning variance in preferred directions. This effect can also be expected to favour homoplasy and convergent evolution. New empirical analyses of the carnivoran cranium show that rates of evolution, in contrast, are not strongly influenced by phenotypic integration and show little relationship to morphological disparity, suggesting that phenotypic integration may shape the direction of evolutionary change, but not necessarily the speed of it. Nonetheless, phenotypic integration is problematic for morphological clocks and should be incorporated more widely into models that seek to accurately reconstruct both trait and organismal evolution.

Phenotypic plasticity in haptoral structures of Ligophorus cephali (Monogenea: Dactylogyridae) on the flathead mullet (Mugil cephalus): a geometric morphometric approach

Evaluating phenotypic plasticity in attachment organs of parasites can provide information on the capacity to colonise new hosts and illuminate evolutionary processes driving host specificity. We analysed the variability in shape and size of the dorsal and ventral anchors of Ligophorus cephali from Mugil cephalus by means of geometric morphometrics and multivariate statistics. We also assessed the morphological integration between anchors and between the roots and points in order to gain insight into their functional morphology. Dorsal and ventral anchors showed a similar gradient of overall shape variation, but the amount of localised changes was much higher in the former. Statistical models describing variations in shape and size revealed clear differences between anchors. The dorsal anchor/bar complex seems more mobile than the ventral one in Ligophorus, and these differences may reflect different functional roles in attachment to the gills. The lower residual variation associated with the ventral anchor models suggests a tighter control of their shape and size, perhaps because these anchors seem to be responsible for firmer attachment and their size and shape would allow more effective responses to characteristics of the microenvironment within the individual host. Despite these putative functional differences, the high level of morphological integration indicates a concerted action between anchors. In addition, we found a slight, although significant, morphological integration between roots and points in both anchors, which suggests that a large fraction of the observed phenotypic variation does not compromise the functional role of anchors as levers. Given the low level of genetic variation in our sample, it is likely that much of the morphological variation reflects host-driven plastic responses. This supports the hypothesis of monogenean specificity through host-switching and rapid speciation. The present study demonstrates the potential of geometric morphometrics to provide new and previously unexplored insights into the functional morphology of attachment and evolutionary processes of host-parasite coevolution.

Genetic and morphological variation in sexual and asexual parasitoids of the genus Lysiphlebus - an apparent link between wing shape and reproductive mode

Background

Morphological divergence often increases with phylogenetic distance, thus making morphology taxonomically informative. However, transitions to asexual reproduction may complicate this relationship because asexual lineages capture and freeze parts of the phenotypic variation of the sexual populations from which they derive. Parasitoid wasps belonging to the genus Lysiphlebus Foerster (Hymenoptera: Braconidae: Aphidiinae) are composed of over 20 species that exploit over a hundred species of aphid hosts, including many important agricultural pests. Within Lysiphlebus, two genetically and morphologically well-defined species groups are recognised: the “fabarum” and the “testaceipes” groups. Yet within each group, sexual as well as asexual lineages occur, and in L. fabarum different morphs of unknown origin and status have been recognised. In this study, we selected a broad sample of specimens from the genus Lysiphlebus to explore the relationship between genetic divergence, reproductive mode and morphological variation in wing size and shape (quantified by geometric morphometrics).

Results

The analyses of mitochondrial and nuclear gene sequences revealed a clear separation between the “testaceipes” and “fabarum” groups of Lysiphlebus, as well as three well-defined phylogenetic lineages within the “fabarum” species group and two lineages within the “testaceipes” group. Divergence in wing shape was concordant with the deep split between the “testaceipes” and “fabarum” species groups, but within groups no clear association between genetic divergence and wing shape variation was observed. On the other hand, we found significant and consistent differences in the shape of the wing between sexual and asexual lineages, even when they were closely related.

Conclusions

Mapping wing shape data onto an independently derived molecular phylogeny of Lysiphlebus revealed an association between genetic and morphological divergence only for the deepest phylogenetic split. In more recently diverged taxa, much of the variation in wing shape was explained by differences between sexual and asexual lineages, suggesting a mechanistic link between wing shape and reproductive mode in these parasitoid wasps.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0293-5) contains supplementary material, which is available to authorized users.