Skip to main content

The emergence of ecotypes in a parasitoid wasp: a case of incipient sympatric speciation in Hymenoptera?

This article has been updated

Abstract

Background

To understand which reproductive barriers initiate speciation is a major question in evolutionary research. Despite their high species numbers and specific biology, there are only few studies on speciation in Hymenoptera. This study aims to identify very early reproductive barriers in a local, sympatric population of Nasonia vitripennis (Walker 1836), a hymenopterous parasitoid of fly pupae. We studied ecological barriers, sexual barriers, and the reduction in F1-female offspring as a postmating barrier, as well as the population structure using microsatellites.

Results

We found considerable inbreeding within female strains and a population structure with either three or five subpopulation clusters defined by microsatellites. In addition, there are two ecotypes, one parasitizing fly pupae in bird nests and the other on carrion. The nest ecotype is mainly formed from one of the microsatellite clusters, the two or four remaining microsatellite clusters form the carrion ecotype. There was slight sexual isolation and a reduction in F1-female offspring between inbreeding strains from the same microsatellite clusters and the same ecotypes. Strains from different microsatellite clusters are separated by a reduction in F1-female offspring. Ecotypes are separated only by ecological barriers.

Conclusions

This is the first demonstration of very early reproductive barriers within a sympatric population of Hymenoptera. It demonstrates that sexual and premating barriers can precede ecological separation. This indicates the complexity of ecotype formation and highlights the general need for more studies within homogenous populations for the identification of the earliest barriers in the speciation process.

Peer Review reports

Background

To understand which barriers initiate speciation is one of the main questions of evolutionary research [1, 2]. For many decades it was assumed that geographic isolation of populations is the essential first step for speciation. Speciation in sympatry, i.e. speciation without geographic isolation, was considered to occur only in special cases like polyploidy [3] or hybrid speciation [4], but was excluded as main speciation mechanism, because gene flow should counteract selection [5]. Only by the end of the last century the concept of sympatric speciation became increasingly accepted [6, 7]. By then, compelling examples of sister taxa with substantial reproductive isolation have accumulated, that must have emerged in overlapping geographic ranges without an allopatric phase, according to the criteria defined by [1]. Now, the debate has shifted to the question, how frequently sympatric speciation occurs [7].

While allopatric speciation is probably driven by incompatibilities according to the Dobzhansky–Muller model [8], and mutation order processes [9], there are numerous theoretic models on sympatric speciation [6]. The dominant idea is that disruptive natural selection leads to assortative mating and reproductive isolation between subpopulations [1, 7]. Prominent examples are ecologically divergent fish [10, 11] or phytophagous insects, which separated via host shift [12, 13]. Specifically in phytophagous insects, sympatric host races, i.e. intermediate stages between polymorphic populations and full species [14], are taken as evidence that genetically differentiated populations can exist despite gene flow which indicates that sympatric speciation is common [12, 14]. As an alternative to natural selection, it is also assumed that disruptive sexual selection, e.g. runaway sexual selection processes [15] or sexual conflict [16] might lead to separation of subpopulations in sympatry. However, because evidence so far indicates that this process is often the consequence of an interaction between sexual and natural selection [17,18,19,20,21] (but also see [22]) it was suggested to drop sexual selection as unique speciation mechanism [21].

The Hymenoptera (bees, wasps and ants) occur in almost all terrestrial ecosystems in large species numbers where they often play a fundamental ecological role as herbivores, pollinators, or natural enemies [23]. They comprise about 153.088 species [24], which is about 8% of all described species [25]. The species richness within the Hymenoptera is mostly due to the monophyletic group of parasitoid wasps [26]. These reproduce by laying their eggs on or in their hosts, mostly other insects, which are consumed and killed by the developing parasitoid larvae [27, 28]. Currently, there are about 90.000 described species of parasitoid wasps (calculated according to data from [24, 29]). However, parasitoid wasps are heavily understudied due to their typical small size, which hampered traditional taxonomic studies based on morphology. Therefore, the actual number is probably much higher. Cryptic species are abundant [30] and new species are continuously being discovered (e.g. [31,32,33,34,35,36]). Based on the observation that most parasitoids are very host specific and almost all holometabolous insect species are attacked by at least one or even several parasitoid wasp species, [29] calculated that parasitoid wasps comprise between 833,000 to 1,107,487 species. This makes Hymenoptera the most species-rich order of all animals.

The reason for the high species diversity in parasitoid wasps is unclear. A comparison of the diversity in carnivorous parasitic insect lineages with their non-parasitic sister groups found no indication that parasitic insects diversify more rapidly than predatory, saprophagous, or phytophagous insects. Therefore the authors concluded that there is no evidence that parasitism itself is the cause of the spectacular diversity of parasitic Hymenoptera [37]. Obviously, other traits must have promoted speciation in this group. Already in 1968, Askew hypothesized for the largest parasitoid superfamily Chalcidoidea with 24.788 described species [24] and an estimated number of 500.000 species [38], that one key factor for their high speciation rate and high species numbers is inbreeding, i.e. mating between siblings [39]. He suggested that it could increase reproductive isolation of lineages by genetic drift [40, 41], and promote speciation similar to geographic barriers in allopatric populations. Interestingly, this hypothesis is analogous to the idea that the high diversity in the ambrosia beetle tribe Xyleborini results from haplodiploidy and reduced gene flow between lineages due to close inbreeding [42]. While inbreeding can cause severe fitness reductions in many hymenopterous species having single-locus complementary sex determination (sl-CSD), which is ancestral in Hymenoptera [43], it is frequent in Chalcidoidea which miss sl-CSD [44] and where deleterious alleles are eliminated in haploid males [45]. In Chalcidoidea, females often place most or all of their eggs in one host or host patch and the emerging offspring readily mate with each other [27, 28, 39]. In addition, many species are monandrous, i.e. females mate only once [46, 47] and a second mating with non-related males after dispersal from the natal patch is rare. As an adaptation to inbreeding, females produce offspring with a strongly female biased sex ratio to avoid competition between their male offspring [48, 49]. The hypothesis that inbreeding is a major cause for the diversity of Chalcidoidea was formulated already in 1968 and is regularly cited in the relevant literature on speciation of Hymenoptera or parasitoid wasps [27, 28, 50,51,52]. However, it has never been explicitly studied. Specifically in Chalcidoidea it is unclear if very early reproductive barriers between sympatric inbreeding female lines emerge before any other barriers, such as ecological separation, as would be predicted by the hypothesis.

Despite their high species numbers and their specific biology, there are only few studies on speciation and the evolution of reproductive barriers in Hymenoptera. So far, sympatric speciation was convincingly demonstrated only for three braconid species [53,54,55,56]. The most detailed studies on the evolution of reproductive barriers exist for the genera Lariophagus (Chalcidoidea: Pteromalidae), parasitoids of beetle larvae [52, 57, 58], and Nasonia (Chalcidoidea: Pteromalidae), parasitoids of fly pupae in bird nests or on carrion [59,60,61,62,63,64,65]. Because the full genome is sequenced for three Nasonia species [66], this genus serves as a model for other biological traits, e.g. wing morphology [67], circadian rhythm and diapause [68, 69], chemical communication [64, 70], sex determination [71], evolution of venom [72], and learning [73, 74]. Reproductive barriers have been studied in Nasonia only at the species level [61, 75,76,77,78]. They are based mainly on genetic incompatibilities in nuclear-nuclear or nuclear-mitochondrial interactions according to the Dobzhansky–Muller model [59, 62, 79], and on endosymbiont induced cytoplasmic incompatibility (CI). CI is assumed to have initiated the separation between N. longicornis and N. giraulti [60] and represents the main reproductive barrier between all Nasonia-species in general [80]. While CI in Nasonia is caused by different strains of Wolbachia, which is the best studied reproductive manipulating endosymbiont [81], there are at least three other endosymbionts which are able to cause the same effect [82,83,84,85]. In crosses between infected males and uninfected females (unidirectional), or between males and females infected with different bacterial strains (bidirectional) all these endosymbionts cause mortality of the fertilized egg (female mortality type) or, in haplodiploid Hymenoptera, the conversion of a diploid female egg into a haploid male egg (male development type) [86]. In addition, CI can also be caused by a higher endosymbiont density in male spermatocytes as compared to the female egg [87].

So far, there are no studies on reproductive barriers within one species of the Nasonia complex, e.g. on the level of ecotypes. For N. vitripennis (Walker 1836), evidence for the existence of two sympatric ecotypes with small molecular differences using RAPD-markers and some reproductive isolation due to a reduction in F1-hybrid females and hybrid breakdown in the F2-generation was provided by [88]. One ecotype was assumed to parasitize fly pupae in bird nests, mostly Protocalliphora flies which are parasitic on nestlings but also carrion flies feeding on dead nestlings [89,90,91,92]. The other ecotype attacks only carrion flies (genera Calliphora, Lucilia etc.) next to carrion [93]. This suggests that ecological separation in sympatry plays a role in divergence in N. vitripennis. However, detailed studies on these ecotypes are missing and it is unclear whether ecological separation is in fact the first barrier to arise. As described above for Chalcidoidea parasitoids in general, inbreeding is common in N. vitripennis. Females mate only once, mostly with their brothers directly after emergence from the host, or even within the host [94]. While males are brachypterous, cannot fly and therefore stay at the natal patch [95], females disperse after mating to locate new host patches, thereby covering distances up to at least 2 km [96]. Most females parasitize hosts only in one patch [97] and lay usually between 15 and 30 eggs per host [98] and up to 800 eggs in total within their lifetime [95] of about 2 to 3 weeks [99].

We studied early barriers within a population of N. vitripennis to elucidate the very early barriers of separation leading to ecotype formation and speciation. To exclude the influence of geography, we studied the population within the area of a local park, i.e. in sympatry according to most definitions [5, 100]. Ecological isolation as a premating barrier was studied with respect to phenology, differences in the reaction to olfactory host habitat cues, and host acceptance. These are important steps for successful reproduction in parasitoids [101, 102]. In addition, we studied premating sexual isolation in mating experiments, and postmating isolation by assessing F1-female offspring within and between inbreeding strains and ecotypes. To address genetic divergence within the population, we also studied the population structure using microsatellites.

Results

Parasitization in the field

Bait bags with fly pupae as hosts for N. vitripennis were exposed in the 35 ha large area of Hohenheim Park (Stuttgart, Germany) in bird nest boxes and next to small carrion hidden on the ground to attract wasps for parasitization. After 1 week, fly pupae were taken to the laboratory where wasps emerged, indicating parasitization in the field. Fly pupae exposed in bird nests, parasitization started in April, increased until the end of May and dropped until the end of July (Fig. 1). On carrion, parasitization occurred during the whole sampling period until September with some breaks in between, probably due to adverse weather conditions. Host pupae from control baits that were attached to twigs without surrounding nest material or carrion never yielded any N. vitripennis offspring.

Fig. 1
figure 1

Percentage of Lucilia sericata baits that had been parasitized by Nasonia vitripennis females in bird nests (yellow) and next to carrion (blue) over a sampling period of 24 weeks (end of April–mid of October 2012) in the Hohenheim Park

Ecological isolation

Wasps that emerged from host pupae, which were exposed in the field in the bait bags (see above), were used to establish laboratory strains for experiments. For each strain, only wasps emerging from one host pupa were taken into culture. Because wasps emerging from one host pupa are the offspring from one or two females [97], these strains simulated the natural situation in which offspring emerging from one host are inbreeding [97], but still maintain at least some genetic diversity. Up to the experiments, all strains were reared under identical conditions in the laboratory on the same host for 20 to 25 generations. To prevent the impact of experience, experimental wasps were dissected as pupa from the host 1 to 2 days prior to hatching to prevent early learning of cues [57].

In olfactometer experiments we studied 11 strains originating from bird nests and 10 strains from carrion. For each strain, about 20 females were tested for their response to nest and carrion odour. Each female was tested only once. The habitat of origin had a strong effect on the reaction of the wasps. More nest strains than carrion strains significantly reacted to nest odour (Fisher’s Exact Test, p = 0.003), and more carrion strains than nest strains significantly reacted to carrion odour (Fisher’s Exact Test, p = 0.012) (Fig. 2). Statistics for all tested strains are given in Additional file 1: Tables S1 and S2. From the strains which were used in all the further experiments (see below), nest strains N3 and N9 only reacted significantly to odours from bird nests, carrion strains A1, A7, and A19 only reacted to carrion, and nest strain N2 reacted to both (Additional file 1: Table S1).

Fig. 2
figure 2

Number of strains of Nasonia vitripennis with significant reaction to the odour of bird nests and carrion in olfactometer experiments. Wasp strains originate from bird nests (yellow) or next to carrion (blue). For each strain about 20 wasps were tested. Coloured parts of the bars refer to the number of reacting strains

Host acceptance was tested with pupae of L. sericata carrion flies (Fig. 3), but not with pupae from nest flies. Parasitic nest flies of Protocalliphora, the putative main host of N. vitripennis in bird nests, are difficult to rear in the laboratory. Wasps from three strains originating from bird nests (N2, N3, N9) and three strains originating from carrion (A1, A7, A19) were tested. Generally, the number of parasitized fly pupae and offspring emerging from the pupae was higher in the three carrion strains compared to the three nest strains (Additional file 2: Figs. S1, S2). Generalized linear models (family “poisson”) followed by ANOVA revealed no difference between strains within each group (nest or carrion). Therefore, data from all nest strains were pooled and compared to the pooled carrion strains using Whitney–Mann U-test, as data were not normally distributed. This revealed highly significant differences between wasps from nests strains and wasps from carrion strains for the number of parasitized fly pupae (Fig. 3; U-test: W = 3189; p = 0.008) and the number of offspring emerging from the pupae (U-test: W = 3047; p = 0.004).

Fig. 3
figure 3

Number of pupae of Lucilia sericata carrion flies parasitized by Nasonia vitripennis wasps (A) and offspring emerging from these pupae (B). Wasps were collected in bird nests (yellow, n = 90) or next to carrion (blue, n = 90). The box and whisker plots show minimum, maximum, 1st and 3rd quartile, median and mean as asterisk. **: significant difference at p < 0.01 (Mann–Whitney U-test)

Genetic analyses

A population genetic analysis using microsatellites was performed using wasps from six laboratory strains (strains from bird nests: N2, N3, N9; strains from carrion: A1, A7, A19; one wasp per strain) as well as other wasps that had emerged from field-collected host baits in bird nests (n = 71) and next to carrion (n = 26) and were not taken into culture. Each single wasp originated from a different host pupa and a different bait. Thus, each analysed female is representative of one female strain. In total, 103 females were studied. The analysis with STRUCTURE v2.3.4 was performed three times. In all three runs there was strong support for either three (k = 3) or five (k = 5) subpopulation clusters (Additional file 2: Fig. S3). The estimated membership coefficients for each individual in each cluster were almost identical in the three runs. Most individuals of subpopulation 1 originated from bird nests, and most individuals of subpopulations 2, 3, 4, and 5 were found next to carrion (Fig. 4A, B). Laboratory strains N2 and N9 were assigned to subpopulation 1, strain A1 to subpopulation 2, strains N3 and A7 to subpopulation 3, and strain A19 to either subpopulation 2 (k = 3) or 5 (k = 5). In subpopulations 1, 2, and 3, between 12 and 30% of the individuals were collected in the alternative habitat and there is no significant difference between subpopulations in the distribution of individuals between main and alternative microhabitat (Fisher Exact test for k = 3: p = 1; k = 5: p = 0.3298; for data see Additional File 1: Table S3).

Fig. 4
figure 4

Results from a population genetic analysis using microsatellites of a population of the parasitoid wasp Nasonia vitripennis in the park of Hohenheim, Germany, using STRUCTURE v2.3. Wasps were collected in bird nests or next to carrion. In the bar plots each individual is represented by a vertical line coloured according to its probability of assignment to one of either k = 3 (A) or k = 5 (B) subpopulations. Origin of wasps from bird nests or carrion is indicated above each graph. The affiliation of experimental strains N2, N3 N9, A1, A7, and A19 to subpopulations within the dataset is indicated below each graph

Calculations of genetic differentiation revealed considerable inbreeding coefficients (FIS, GIS) within the individuals representing one female strain each, ranging from 0.318 to 0.413 (Table 1). Fst-values between subpopulations defined by their origin from bird nests and carrion are always lower than between subpopulations identified by microsatellites. Between subpopulations identified by microsatellites with k = 3, there is low differentiation between subpopulations 1 and 2, as well as 1 and 3, but higher differentiation between 2 and 3. Higher levels of differentiation also exist between most of the subpopulations identified in the k = 5 model. Hedrick’s standardized G′st values also indicate low differentiation between the two microhabitats, but higher differentiation between the subpopulations identified using microsatellites, with highest values ranging up to 0.809. Subpopulation 2 was excluded from the analysis with k = 5 due to the small number of individuals.

Table 1 Genetic differentiation among putative subpopulations of Nasonia vitripennis in the park of Hohenheim, Stuttgart, Germany

Sexual isolation

Sexual isolation was studied with wasps from strains originating from bird nests and carrion (strains from bird nests: N2, N3, N9; strains from carrion: A1, A7, A19). Single pairs of females and males were tested in all possible combinations. Courtship behaviour by males was observed in all experimental couples, in which the male located the female (n = 1059). Thus, there is no isolation due to male mate selection. There was a small, but highly significant decline in the number of females that accepted males for copulations when males were from different strains than females (Fig. 5A; ANOVA based on generalized linear model, family “binomial”: Chi2 = 25.00, df = 3, p = 0.000). Single comparisons using Tukey-test revealed significant differences between same strain couples and couples from different strains but the same microsatellite cluster, but no differences towards couples from different strains and the same microhabitat. When females and males were from different strains, there were no differences with respect to their origin from the same or from different microsatellite clusters, or from the same or different microhabitats (Fig. 5B; Table 2). Thus, there is a small, significant level of sexual isolation between partners from different strains caused by female mate choice, but no isolation between partners from different microsatellite clusters or different microhabitats.

Fig. 5
figure 5

Occurrence of copulation in N. vitripennis depending on origin of mating partners. A Females and males originate from the same strain (green), or from different strains but the same microsatellite cluster with k = 3 (yellow) or k = 5 (orange) subpopulations, or the same microhabitat (blue). B Females and males originate from different strains but the same microsatellite cluster with k = 3 (yellow) or k = 5 (orange), or from different microsatellite clusters (yellow/orange), or originate from the same (blue) or different microhabitats (white/blue). Coloured parts of the bar indicate the percentage of couples with copulation. Numbers in bars indicate numbers of replicates. Different letters above bars in A indicate significant differences at p < 0.01 (Tukey-test based on a generalized linear model, family “binomial”). In B, ANOVA was based on a generalized linear model, family “binomial”

Table 2 Test statistics of experiments on sexual isolation among females and males from different strains in Nasonia vitripennis

Postmating isolation

Postmating isolation was studied with wasps from the same strains as above (strains from bird nests: N2, N3, N9; strains from carrion: A1, A7, A19). Females and males were mated in all possible combinations. Only three out of the 540 couples did not produce F1-female offspring. The number of F1-female offspring was significantly influenced by the origin of mating partners from the same strain, from different strains but the same microsatellite cluster with k = 3 or k = 5 subpopulations, or the same microhabitat (Fig. 6A; Table 3). The strain and the habitat of the parental female did not influence the number of F1-females (Table 3). Couples from the same strain had significantly more F1-female offspring than couples from different strains but from the same microsatellite cluster with k = 5 subpopulations, and from the same habitat. There was no difference between couples from the same strains and couples from different strains but from the same microsatellite cluster with k = 3 subpopulations.

Fig. 6
figure 6

Influence of origin of the mating partners on number and sex ratio of F1-offspring in the parasitoid wasp Nasonia vitripennis. A F1-female offspring; B F1-male offspring; C F1-total offspring; D sex ratio. Females and males originate from the same strain (green, n = 90), or from different strains but the same microsatellite cluster with k = 3 (yellow, n = 90) or k = 5 (orange, n = 60) subpopulations, or originate from the same microhabitat (blue, n = 180). Box and whisker plots show minimum, maximum, 1st and 3rd quartile, and median. Different lower case letters indicate significant differences (p < 0.05, generalized linear models with family “quasipoisson” followed by ANOVA and Tukey-test. For test statistics see Tables 3 and 4

Table 3 ANOVA test statistics of experiments on the influence of different factors on the number of F1-female offspring as postmating barrier in Nasonia vitripennis

To study if the reduction in female offspring is due to the number of eggs laid during the experimental period, due to mortality of the F1-female offspring, or due to partial CI, we also analysed the number of F1-male offspring, the total number of F1-offspring (males and females), and the sex ratio. Like with F1-females, couples from the same strain had significantly more F1-male and F1-total offspring than couples from different strains but from the same microsatellite cluster subpopulations, and from the same habitat (Fig. 6; Table 4). There were no differences in sex ratio between the treatments.

Table 4 Test statistics of experiments on postmating isolation between strains in Nasonia vitripennis

To analyse data of couples with females and males from different strains, we used generalized linear mixed model (“negative binomial”) with affiliation to the same or different microsatellite clusters (either based on k = 3 or k = 5 populations), origin of the wasps (same microhabitat, different microhabitats), strain of parental female, and habitat of parental female as factors. This revealed a highly significant impact of affiliation to the same or different microsatellite clusters on the number of F1-female and F1-total offspring, but no significant effect of microhabitat, female strain and female habitat (Table 5; for F1-offspring numbers from each strain combination see Additional file 1: Table S4). In addition, subsequent Tukey-test showed that numbers of F1-female and F1-total offspring in couples from the same microsatellite cluster are significantly higher as compared to couples from different clusters, but that there is no difference in offspring numbers between couples from the same vs. couples from different microhabitats (Fig. 7A, C). The number of F1-male offspring in the treatments with partners from different microsatellite clusters or different microhabitats were always slightly lower (Fig. 7B), and there were significant effects of affiliation to the microsatellite cluster k = 3, of microhabitat, female strain, and female habitat (Fig. 7B, Table 4). The overall mean sex ratio of F1-offspring was 0.23 ± 0.14 (mean ± sd; n = 540) and not influenced by either factor (Fig. 7D, Table 5).

Table 5 Test statistics of experiments on postmating isolation between microsatellites clusters and microhabitats in Nasonia vitripennis
Fig. 7
figure 7

The influence of origin of mating partners on number and sex ratio of F1-offspring in the parasitoid wasp Nasonia vitripennis. A F1-female offspring; B F1-male offspring; C F1-total offspring; D sex ratio. Mating partners originated from the same or from different microsatellite clusters (with k = 3 or k = 5 subpopulations), or from the same or different microhabitats. Box and whisker plots show minimum, maximum, 1st and 3rd quartile, and median. Levels of significance are based on generalized linear models, family “negative binomial”, followed by Tukey-test. For test statistics see Table 5

Isolation indices

Isolation indices were calculated using the olfactometer data for ecological isolation, data on the occurrence of copulations for sexual isolation, and data on the number of F1-female offspring for postmating isolation (Table 6). Indices show considerable total isolation between couples from the same strain and couples from different strains but from the same microsatellite cluster or the same microhabitat. This is either based on sexual barriers, postmating barriers, or both. While different microsatellite clusters are separated by postmating barriers, strains from the two different microhabitats are only separated by ecological barriers.

Table 6 Isolation indices for different barriers between Nasonia vitripennis wasps from different strains but same microhabitat, and wasps from different strains and different microhabitats

Discussion

To answer the question which reproductive barriers initiate speciation is difficult for two reasons. First, because several barriers can arise during the speciation process [1] study systems must consist of populations that are still exchanging genes at the very early phase of separation [20]. However, most studies on early barriers have used young species pairs (e.g. [103]), allopatric populations [104], or ecologically diverged subpopulations (e.g. [105]). In these systems separation has already proceeded which prevents the identification of the earliest barrier. Second, initial barriers are expected to be only weak at the beginning. Therefore, many experimental replicates are needed for their identification, which requires study organisms that are available in large numbers. Thus, there are only very few studies that allow conclusions on barriers within populations which have not yet diverged [106]. To identify very early barriers of separation in Hymenoptera, we studied ecological isolation, sexual isolation, postmating isolation and population structure within a sympatric population of the parasitoid wasp N. vitripennis. Our study revealed three levels of isolation within the population: Between female strains, between subpopulation clusters based on microsatellites, and between ecotypes.

Ecological isolation

Our study confirms the existence of two different ecotypes in N. vitripennis, one occurring in bird nests, and the other in carrion. This has already been suggested previously [107, 108]. However, we could not find any premating or postmating reproductive barriers between these ecotypes except of ecological traits, and the carrion ecotype does not seem to represent a distinct genetic cluster (see below). In contrast, as shown in our microsatellite study, it consists of several, genetically distinct subpopulations.

The ecological differences between the ecotypes refer to their phenology (1), their reaction to host habitat cues (2) and their host acceptance (3). In bird nests, wasps were present from May until end of July after nestlings had left, while they occurred on carrion from May until September. We assume that this difference in phenology is due to the fact that nest wasps are adapted to the presence of pupae of Protocalliphora flies which are parasitic on nestlings that occur in nests only until the end of June [109]. In contrast, carrion wasps might be adapted to the longer availability of carrion and pupae from carrion flies. However, we cannot exclude that this finding is an artefact caused by the loss of attractiveness of nest material due to the complete evaporation of the involved volatiles after the second brood has left the nest.

In olfactometer experiments on host habitat location, the reaction of wasps depended on the microhabitat of origin. All strains reacted to the odour of the microhabitat in which they were collected, but only some strains reacted to the odour of the foreign microhabitat. Because wasps were reared for 20 to 25 generations on the same hosts and were dissected before emergence from the pupa to avoid early adult experience [57] this reaction can not have been learned, but must be based either on genetic differences between the strains, or differences in the microbiome [110]. The importance of odours for host habitat location is highlighted by the absence of parasitization in control fly pupae, which were offered without habitat odours in the field. This demonstrates that hosts are only found in those microhabitats, for which an olfactory reaction is present.

In experiments on host acceptance with pupae of the carrion associated fly species L. sericata [111], the number of parasitized pupae and the number of emerging total offspring was significantly higher in wasps from carrion as compared to wasps from bird nests. Again, the higher acceptance of carrion fly pupae by carrion wasps must be inherited, because all experimental wasps were reared on the same hosts for many generations and dissected in the pupal stage from the host. It remains to be examined if the reduced offspring numbers in nest wasps is due to reduced parasitization (host acceptance), or because of a higher mortality during development, i.e. an inferior host suitability of carrion flies. Although technical reasons prevented us from conducting the reciprocal experiment with Protocalliphora, which is the main host of wasps from bird nests, these data clearly demonstrate a better adaptation of wasps collected from carrion to pupae of carrion flies. However, it might well be that the level of ecological segregation is underestimated due to the lack of experiments with Protocalliphora.

Based on these and other traits, the ecotypes meet many criteria for host races of plant feeding insects [14, 112, 113]: They are associated with different habitats and hosts which is demonstrated by their specific reaction to habitat odours, their host acceptance and possibly their phenology. They display a correlation between host choice and mate choice due to their inbreeding behaviour between siblings, at least one (the carrion type) has a higher fitness on its specific hosts, and they occur in sympatry. In contrast to host races, however, only one of the ecotypes is genetically differentiated.

Sexual isolation

While males always performed courtship behaviour regardless of the origin of the females, females were less likely to accept males for copulation when they originated from different strains but the same microsatellite cluster. In contrast, there were no significant differences in the number of copulations between strains from the same or from different microsatellite clusters, and between the same or different microhabitats. Thus, there is slight sexual isolation between inbreeding strains caused by mate choice decisions from the females, but no isolation between microsatellite clusters or microhabitats. This mate choice decision by the female is most likely based on the male mandibular pheromone which is applied on the female antennae by the males during courtship [114]. It must have a genetic basis and can not have been induced by the developmental experience of the females as described for the related L. distinguendus [58], because all wasps were reared under identical conditions and dissected as pupa from the host prior to hatching.

Postmating isolation

Couples consisting of females and males from the same strain had more F1-female offspring as compared to couples from different strains but from at least one of the identified microsatellite cluster or the same microhabitat. Likewise, females produced more F1-female offspring after copulations with males from the same microsatellite cluster as compared to copulations with males from a foreign cluster. The origin of mating partners from the same or from a different microhabitat had no effect. Obviously, postmating reproductive barriers exist between inbreeding female strains, and between microsatellite clusters, but not between ecotypes.

Generally, the reduction in F1-female offspring in couples with partners from different female strains and microsatellite clusters can have different reasons: (1) Genetic incompatibilities according to the Dobzhansky–Muller (D–M) model due to nuclear–nuclear or nuclear–mitochondrial interactions, which are both described for Nasonia [59, 62, 79]. This would result in a reduction of the total number of F1-offspring due to higher mortality of hybrid F1-females during development, and consequently a change in the normal, female biased sex ratio [49] towards more male offspring. (2) Total or partial cytoplasmic incompatibility (CI) induced by Wolbachia or other endosymbionts causing mortality of the hybrid female zygote (female mortality type), or its conversion to a male (male development type) [86]. All types of CI would lead to changes in sex ratio. Female mortality CI would also result in a reduced total number of F1-offspring similar to D–M incompatibilities. Male development CI would result in an increased number of male offspring and no change in the total number of F1-offspring. (3) A lower reproductive investment by females after matings with males from a different strain. (4) A noncompetitive postmating, prezygotic barrier (PMPZ) acting between copulation and fertilization [1]. Because we observed no differences in sex ratio in all experiments on postmating isolation, D–M incompatibilities and CI of any type is unlikely as explanation. Male development CI can also be excluded because the number of F1-male offspring did not increase and the total number of F1-offspring declined. A lower investment by females is also unlikely as reason because females mate only once and should optimize their own offspring numbers, regardless of the quality of their mates. Therefore we favour PMPZ as explanation. PMPZ involves the interaction of gametes and/or reproductive tissue or reproductive proteins and can result in several barriers like reduced sperm transfer, inefficient sperm storage, or the failure of sperm to fertilize eggs [115]. We assume that PMPZ has reduced the amount of sperm available to females for fertilization of eggs. Therefore female laid fewer F1-female eggs and also less male eggs to maintain an optimal sex ratio to avoid local-mate competition [48, 49]. Interestingly, PMPZ has been shown to evolve very fast [116], and to occur as first or only barrier between closely related taxa or populations of Drosophila [115, 117]. We are currently studying this hypothesis as well as further barriers like sexual isolation, and intrinsic postzygotic barriers like inviability and fecundity of hybrids that could also contribute to isolation between strains.

While our data demonstrate a reduction of F1-female offspring after outbreeding, we could not find an indication for inbreeding depression, i.e. a decrease in fitness after mating with siblings. This is remarkable. Generally, inbreeding is assumed to result in higher levels of homozygosity and consequently inbreeding depression [118]. This effect seems to be lower in haplodiploid taxa like Hymenoptera, mites [119], and ambrosia beetles [42] than in diploid species which is explained by prolonged inbreeding and purging of the genetic load by the regular exposition of deleterious mutations to selection in haploid males [42]. Outbreeding, on the other hand, is generally associated with elevated fitness compared to inbreeding. However, above certain levels of parental dissimilarity, outbreeding depression can occur in hybrids, by genetic intrinsic incompatibility due to underdominance, chromosomal rearrangements or deleterious epistatic interactions, or extrinsic postzygotic barriers like maladaptation of hybrids to paternal habitats [120, 121]. While outbreeding depression is well studied in plants, there are much less studies on animals so far. Interestingly, most animal species with outbreeding depression have a similar biology to N. vitripennis. They either perform selfing as equivalent to inbreeding like freshwater snails [122], or nematodes [123], or they are haplodiploid like mites [124], or do both like in ambrosia beetles [42]. Possibly, inbreeding or selfing as well as haplodiploidy are a prerequisite to outbreeding depression.

Genetic population structure based on microsatellites

The analyses of microsatellite markers indicate a remarkable level of inbreeding within the female lines, and the existence of one subpopulation cluster, which is largely restricted to bird nests, and two or four clusters largely restricted to carrion as microhabitat. While a high level of within patch inbreeding was also found in a Swedish population of N. vitripennis [125], the existence of subpopulation clusters within a local population has not been reported before for N. vitripennis. In support of the identified population structure, Fst and Gst-values showed stronger differentiation between the microsatellite subpopulations but only low genetic differentiation between wasps from the two microhabitats. Remarkably, Fst and Gst-values between microsatellite subpopulations are within the range or even higher than values reported by [96] for geographically distant populations of N. vitripennis within North America or Europe, and even between North America and Europe. In addition, our mating experiments revealed significantly lower F1-offspring numbers in couples from different microsatellite subpopulations as compared to couples from the same subpopulation, while this was not the case for couples within and between different microhabitats. Obviously, microsatellite subpopulations are not only genetically different, but also separated by the same postmating barrier that is present between the different female strains.

The population structure of N. vitripennis and its emergence

Taken together, the data demonstrate that our local population of N. vitripennis is structured on three levels: Between inbreeding strains, between genetic clusters based on microsatellite analysis, and between ecotypes. Strains are characterized by considerable inbreeding coefficients, and are separated by slight sexual isolation and postmating isolation from other strains of the same microsatellite cluster and from other strains of the same ecotype. Strains from different microsatellite clusters are also genetically distinct and separated by postmating isolation. In total, barriers between inbreeding strains lead to a total isolation ranging from 0.128 and 0.328, and between microsatellite clusters from almost zero (0.039) to 0.258 (Table 6). In contrast, the ecotypes of N. vitripennis are not separated by sexual isolation or postmating barriers, but only by considerable ecological barriers leading to a total isolation of 0.525 (Table 6).

Remarkably, these ecotypes consist of different microsatellite subpopulations. Subpopulation 1, which is the largest, was mostly collected in bird nests and probably represents the ecotype that is adapted to bird nests as host microhabitat. Accordingly, laboratory wasp strains N2 and N9 from subpopulation 1 showed a significant reaction to odours from bird nests in the olfactometer experiments, and a lower parasitization of fly pupae from carrion flies in experiments on host acceptance. In contrast, subpopulations 2, 3, 4, and 5 were mostly collected next to carrion and are probably adapted to this microhabitat, forming the carrion ecotypes. Therefore, A1, A7, and A19 belonging to subpopulations 2 and 3 (k = 3), or 2, 3, and 5 (k = 5), only reacted to carrion and had a higher acceptance of carrion fly pupae. However, the separation between the two ecotypes and microhabitats is not strict and between 12 and 30% of the individuals of strains 1, 2, and 3 were also collected in the alternative habitat. For instance, N3 was genetically assigned to the carrion subpopulation 3, but was collected in bird nests and only reacted to nest odours. Possibly, this strain is intermediate between the two ecotypes.

The emergence of premating and postmating reproductive barriers between inbreeding strains as first barrier in our population is consistent with Askew’s inbreeding hypothesis [39], stating that inbreeding could lead to reproductive isolation by genetic drift [40, 41]. Alternatively, postmating barriers could also be explained by CI due to qualitative or quantitative differences in the presence of CI-inducing bacteria like Wolbachia, and premating barriers could result from differences in developmental experience influencing mate choice decisions of the females [58]. However, the influence of CI can be excluded by our data on F1-offspring and F1-sex ration, and the role of experience due to our experimental procedure (see above). Therefore, at present we believe that Askew’s inbreeding hypothesis is the most likely explanation. As mentioned above, mating in N. vitripennis occurs mostly between female and male offspring from the same maternal female directly after emergence from the host [97], or sometimes even within the host [94]. This could lead to reproductive isolation and promote separation. Remarkably, other animal species with outbreeding depression also perform inbreeding [42] or selfing [122].

With respect to the evolution of the microsatellite subpopulations and the resulting ecotypes three scenarios are conceivable. Thereby it is unclear if the initial population consisted of bird nest specialists, from which carrion specialists evolved by host switching, or vice-versa, or if it consisted of generalists, that evolved into microhabitat specialists for bird nests and carrion. The fact that all other species of the genus Nasonia occur in bird nests [126] could indicate that N. vitripennis has originated from a bird nest specialist. However, its larger host range comprising pupae of blowflies, fleshflies, houseflies, and others [126] suggests that it is a generalist. In any case, microsatellite subpopulations could have evolved in sympatry from the inbreeding lines and became adapted to their microhabitat, either once for the bird nest ecotype and several times independently for the carrion ecotype (Fig. 8, scenario A), or only once for each ecotype (scenario B). Thereby the adaptation to the microhabitats could be based on new mutations, or on ancestral alleles present in low frequencies in the original population [127]. However, the observed pattern could also be the result of immigration of ecotypes that evolved their adaptations in allopatry (Fig. 8, scenario C). This agrees with the idea that ecological speciation starts in allopatry and continues in sympatry [128]. Note that scenarios 1 and 3 assume that phenotypes evolved repeatedly in genetically separated clusters [129], i.e. the reaction of wasps from strains A1, A7, and A19 to odours from carrion, and their higher acceptance of carrion fly pupae would have evolved independently. To identify the most likely of these hypothetical scenarios, an analysis of the ecology and the population genetics of N. vitripennis populations from a larger geographic area, e.g. Central Europe, is required.

Fig. 8
figure 8

Hypothetical scenarios for the emergence of the observed pattern of subpopulations and ecotypes of Nasonia vitripennis. Scenario A Inbreeding lines evolve into subpopulations in sympatry and become independently adapted to their respective microhabitats. Scenario B Inbreeding lines evolve into two subpopulations in sympatry. One subpopulation becomes adapted to bird nests as microhabitat and one subpopulation becomes adapted to carrion. Subsequently, the carrion subpopulations splits up into 2 or 4 separate subpopulations. Scenario C Subpopulations became adapted to their respective microhabitats in allopatry and become sympatric by immigration into the study area

While it is unclear if the two different ecotypes evolved in sympatry, our results clearly demonstrate that they are able to co-exist in sympatry. This supports the idea of sympatric speciation also for parasitoid wasps, which was only shown for the braconid species Diachasma alloeum, Utetes canaliculatus, and Diachasmimorpha mellea that have separated in sympatry as a consequence of the divergence of their Rhagoletis hosts [53, 55]. If true, the putative separation in N. vitripennis was not induced by host radiation as in these braconids, but occurred independently due to the presence of two different host-microhabitats. It would therefore represent the only parasitoid example demonstrating that incipient sympatric host race or ecotype formation is possible even in the absence of sympatric host radiation.

Conclusions

Our study sheds light on the earliest steps of ecological separation in a hymenopterous parasitoid. In contrast to the assumption that separation by natural divergent selection generally starts with ecological barriers [130, 131], we found sexual isolation and postmating isolation as first barrier between inbreeding strains that are not separated ecologically. This points to inbreeding as a hitherto underestimated mechanism in evolution and speciation [39, 42, 132,133,134]. In addition, the two ecotypes of N. vitripennis, one occurring in bird nests and the other on carrion, do not form two, but several genetically distinct, sympatric subpopulations. This indicates the complexity of ecotype formation. To further substantiate the hypothesis that this separation occurred in sympatry, the ecology and the population structure of N. vitripennis populations from Central Europe is required. Finally, the study highlights the general need for more studies within seemingly homogenous, sympatric populations for the identification of the earliest barriers in the speciation process [20]. Because barriers at the onset of separation are very weak per definition, a large number of replicates are required to obtain sufficient statistical power for the identification of these barriers. Unfortunately, this is very difficult to realize in most study systems.

Methods

Field experiments

To study phenology and to establish laboratory cultures of N. vitripennis, wasps were collected using bait bags with fly pupae in the park of the University of Hohenheim (Stuttgart, Germany). The park has 35 ha of meadows, trees, understory, and ponds and is very species-rich. A 24-h survey in 2013 (“Geo-Tag der Artenvielfalt”) revealed the presence of at least 365 plant species and 846 animal species in the area. N. vitripennis was collected from the end of March to the end of October 2012 in bird nests and next to carrion. Bait bags were made of aluminium gauze (4.5 cm × 1 cm, 1.5 mm mesh width) and contained five fly pupae of the green blowfly Lucilia sericata (Meigen 1826) (Diptera: Calliphoridae). Fly pupae were obtained from a pet store as larvae, allowed to pupate, and stored at 4 °C until being used. Single bait bags were placed below the nesting material in 29 bird boxes (Schwegler, Germany) situated up to 3.5 m high on trees. The boxes contained nests of tits (Parus major L., Cyanistes caeruleus L., Parus palustris L.) from the current breeding season. One bait bag was also placed next to each of 15 steel grid cages (215 × 115 × 55 mm) containing fresh carrion of small rodents (Mus musculus L., Microtus arvalis (Pallas 1778), Meriones unguiculatus (Milne-Edwards 1867)). Eight control baits without nest boxes or carrion were attached to twigs at different locations in the park. After 1 week, fly pupae from the baits were replaced by new fly pupae, taken to the laboratory and incubated at 25 ± 1 °C and 60% r.h. in a thermostatic cabinet (Aqualytic® AL-185-4). These conditions were used for all cultures. Emerging wasps were determined under a microscope using [135].

Experimental wasps

To establish inbreeding strains for laboratory cultures, wasps emerging from single host pupae from the bait bags were placed in Petri dishes with L. sericata pupae (approx. 6 g) that had been killed by freezing at − 20 °C for a minimum of 1–2 days. By using only wasps from one single host pupa for each strain, we simulated the natural situation in which offspring emerging from one host are often inbreeding [97], but still maintains at least some genetic diversity. Note that this approach is different from generating artificial isofemale lines, which are used to produce genetically homogenous individuals to study the genetics of certain traits (e.g. [136]). Strains from nest boxes and from carrion were kept in different rooms at constant light. Cultures were provided with dead host pupae three times a week. For most experiments, except for the reaction to host habitat odours, we selected three nest strains (N2, N3, N9) and three carrion strains (A1, A7, A19) from the laboratory cultures. These strains were collected in different collection weeks and with different baits. To prevent the impact of experience, all strains had been reared under identical conditions on dead L. sericata pupae for many generations (see below). In addition, experimental wasps were dissected as pupa from the host pupae 1 to 2 days prior to hatching. Therefore, all identified differences between strains must be genetically fixed and can not be attributed to learning behaviour [57].

Reaction to host habitat odours in the olfactometer

In a four-chamber olfactometer we tested for the presence of a positive reaction to the host habitat odour of bird nests in nine nest and 10 carrion strains, and for the presence of a positive reaction to carrion odour in 11 nest strains and 10 carrion strains (Additional file 1: Table S1). Up to the experiments, the strains had been reared for 20 to 25 generations after being collected. Experimental wasps were dissected from the host to exclude experience as described above and kept in groups of 10–20 individuals in perforated Eppendorf tubes prior to the experiments. Each wasp was used only once. We tested the reaction of about 20 female wasps for each strain. Odour samples of nest material had been taken from a nest box in which parasitization of fly pupae by N. vitripennis was recorded. Odour samples from carrion originated from a dead mouse that was exposed in the field for 2 weeks in June/July and was infested by microorganisms, feeding maggots of L. sericata, and other arthropods. Until being used, the nest material and the mouse cadaver were stored in odour-impermeable roasting bags at 4 °C (nest material), or at − 23 °C (carrion). To obtain standardized odour samples, 3.5 g nest and carrion material was tested for its attractiveness to females from nest strain N29 and carrion strain A7, respectively. Odour extracts were produced by extracting 3.5 g attractive material in a 100 ml bottle (Duran®) together with 20 ml dichloromethane at room temperature. After 24 h the extract was filtered and stored in in glass vials (22 ml, Supelco®) at − 23 °C.

The static 4-chamber olfactometer according to [137] consists of a polyvinyl chloride cylinder, which is divided into four chambers (Fig. 9). The cylinder is covered with metal gauze (mesh size 0.2 mm) as walking arena on which the wasps can move freely above the four chambers. Every chamber contains a glass dish (Ø 40 mm, height 7 mm). In one chamber the dish contained a filter paper treated with host habitat odour extract (test), the dish in the opposite chamber contained a filter paper treated with dichloromethane (control). The intermediate chambers served as transition zones and remained empty.

Fig. 9
figure 9

Static 4-chamber olfactometer to test the reaction of female N. vitripennis to odours of bird nest and carrion. 1: glass plate; 2: walking arena; 3: test field; 4: control field; 5: transition zone

Individual wasps were separated into Eppendorf tubes 15 min prior to the experiments. For the experiment, the Eppendorf tube with the wasp was opened and placed in the centre of the arena, which was covered with a glass plate (thickness 2 mm). As soon as the wasp entered the arena, the behaviour (walking, sitting, other) and the location of the wasp (above the test, control, or the two transition chambers) was recorded for 5 min using the software “The Observer® 5.0” (Noldus, 2003, Wageningen, NL). A positive reaction to a host habitat odour was recorded for those strains where allocation times above the test chamber were significantly longer than times above the opposite control chamber. Females which spent < 150 s stationary were considered inactive and excluded from the analysis. After each experiment, the top and bottom of the arena, as well as the bottom of the glass plate were cleaned with 70% ethanol and the olfactometer was rotated by 90°. A maximum of five females from each strain were tested per day. To account for daily response variation due to weather [138], positive controls were performed on each experimental day with females from strain N29 for nest odour and A7 for carrion odour. When these wasps failed to respond, all data from that day were discarded.

Host acceptance of pupae of carrion flies

Experiments were performed with six strains (N2, N3, N9, A1, A7, A19) from the laboratory culture. After being collected, wasp strains were reared on identical host species (dead L. sericata pupae) for 80 to 100 generations. As described above, wasp pupae were dissected from the host pupae and placed into Eppendorf tubes with perforated lids in groups of 10 female and four male pupae for eclosion and subsequent mating. Adult wasps were transferred into Petri dishes (ø 55 mm) and allowed to feed for 24 h on sucrose solution (10%) on filter paper. Thus females had the opportunity to mate and to take up energy for egg maturation. For experiments, single females were placed in Petri dishes (ø 55 mm) with five pupae of L. sericata each. After 6 h, wasps were removed and pupae were kept until emergence. The numbers of fly pupae that had been parasitized by each female and the offspring of each female were counted after 3 weeks. From each strain, 30 wasps were tested.

Sexual isolation

To study sexual isolation, mating experiments were performed with females from three strains originating from bird nests (N2, N3, N9) and three strains originating from carrion (A1, A7, A19) in all possible combinations. After being collected wasp strains were reared on identical host species (dead L. sericata pupae) for 60 to 75 generations. To obtain virgin insects for the experiments, wasp pupae where dissected from host pupae as described above, sexed and kept separated by sex in groups of up to 10 individuals in Eppendorf tubes. Experiments were performed in arenas consisting of the bottom of a Petri dish (ø 55 mm) covered with a glass plate. The arena was placed under a digital camcorder (TK-1480BE Color Video Camera, JVC, Yokohama, Japan). Single virgin females (max. age 3 d) were placed in the arena and allowed to calm down for several minutes before one virgin male was added. The observation was finished when copulation was observed or after 10 min. Per combination between 20 and 52 couples were studied.

Postmating isolation

Experiments were performed with females from three strains from bird nests (N2, N3, N9) and three strains from carrion (A1, A7, A19), which were reared on dead L. sericata pupae for 20 to 100 generations. For copulation, single female and male wasps from all strains in all combinations were placed into small chambers (40 × 20 × 15 mm) made of acrylic glass closed with a transparent lid. After copulation was observed the wasps were left in the chamber until the next day. Then, single females were transferred into Petri dishes (60 × 15 mm) with 10 pupae of L. sericata for parasitization and kept there until their death. No additional food was provided. Host pupae had been killed by freezing at − 20 °C and subsequently thawed at room temperature (see above) to rule out potential differences between the strains based on any ability to suppress the host immune defence. Care was taken to make sure that the host puparia were all of the same size and quality. After 4 weeks the offspring was sexed and counted.

Calculation of isolation indices

Reproductive isolation (RI) is the strength of any reproductive barrier and represents an estimate how much gene flow is reduced by this barrier [1]. Based on Sobel and Chen [139] RI varies from − 1 (the barrier allows for gene flow only between heterospecifics) over 0 (gene flow is random) to 1 (only gene flow between conspecifics). RI was calculated according to the following formula:

$$1-2*\frac{(H)}{\left(H\right)+(C)}$$

H refers to events, which enable gene flow between different inbreeding lines, microsatellite clusters or ecotypes (i.e. frequency of mating or the number of offspring after mating), while C refers to events, which enable gene flow between different inbreeding lines, microsatellite clusters or ecotypes. RI for ecological isolation was calculated based on the data from olfactometer experiments. To calculate RI for sexual isolation, we used the data of the mating experiments presented here. To calculate RI for postmating isolation we used the mean number of female F1 offspring.

Total isolation (T) varies between 0 (no isolation) and 1 (full isolation between populations). It was calculated as the sum of reproductive isolation caused by all different barriers during the life history of the wasps, starting with the premating barriers ecological isolation, and sexual isolation, and ending up with postmating isolation, i.e. the number of F1 ♀ offspring. For each barrier, we used the absolute contribution (AC). It considers this part of gene flow that has not already been prevented by previous stages of reproductive isolation [140]. It was calculated as follows:

AC for the first barrier: AC1 = RI1.

AC for the second barrier: AC2 = RI2 * (1 − AC1).

AC for the third barrier: AC3 = RI3 * [1 − (AC1 + AC2)].

Total isolation (T) was calculated using the following formula:

$$T=\sum_{i=1}^{m}ACi$$

Population genetic analysis

A population genetic analysis was performed with 103 females to study the population substructure. We used one individual wasp from each of the six laboratory strains N2, N3, N9, A1, A7, and A19 as well as 96 wasps (71 from bird nests and 26 from carrion) that had emerged from field-exposed host baits with alive fly pupae (see above) that were not used for further cultures. From these, one female per bait was analysed, each representing one female strain. Females were extracted using the DNeasy Blood & Tissue Kit (Qiagen®, Venlo, NL) following the manufacturer’s instructions. The tissue was mechanically ground and proteinase K enzyme was added. After incubation at 56 °C, the solution was washed in repeated steps with buffer solutions and subsequent centrifuging. All PCRs were carried out in TGradient® 96 (Biometra, Göttingen, GER) and Techne PrimeG® (Cole-Parmer, Stone, UK) thermocycling machines.

Eight microsatellite markers were selected based on [141]. Loci on all five chromosomes of N. vitripennis were chosen to achieve a better coverage of the genome and to account for different selection strengths between chromosomes (Table 7). The Qiagen® Multiplex PCR Kit (Qiagen, Venlo, NL) was used for amplification of microsatellite loci in 96-well plates with a volume of 12.5 µl per reaction (6.25 µl Qiagen® Multiplex PCR Master Mix 2×, 2.5 µl Qiagen® Q-Solution, 1.25 µl primer mix (10 mM) and 1.5 µl RNAse-free water). The primer mix contained two primer pairs with different fluorescent markers (HEX: Hexachloro-Fluorescein, 6FAM: 6-Carboxyfluoresceine, Applied Biosystems, Waltham, MA, US). All sequences were amplified using the same PCR profile. The reaction was initiated through incubation at 95 °C for 15 min, followed by 25 cycles at 94 °C for 30 s, 60 °C for 90 s and extension at 72 °C for 60 s, with a final elongation at 72 °C for 45 min and cool-down to 2 °C until further usage. Fragment length analysis was carried out by Eurofins Genomics (Ebersberg, Germany) using the ROX500 size standard and an ABI-D filter system on a ABI 3130 XL sequencing machine (Applied Biosystems, Waltham, MA, US). Raw fragment length data were examined and the quality of individual signals/peaks checked with PeakScanner v1.0 (Applied Biosystems®, Waltham, MA, US) according to [142]. The software MICROCHECKER v2.2.3 [143] was used to check the data set for false allele sizes as well as existence of null-alleles or allelic dropout.

Table 7 Markers used in the microsatellite analysis

Analyses were performed using the software STRUCTURE v2.3.4 [144] and the online tool CLUMPAK [145]. We tested a k-range of 1–12, each k with 30 iterations. Running length was 1.000.000 MCMC replicates with a burn-in of 100.000 using the admixture model with correlated allele frequencies [146]. As sampling location information [147] we included the origin of the single wasps from bird nests or carrion. The structure analysis was repeated three times to check for consistency of results. Genetic differentiation among populations was quantified using GenAlEx 6.502 [148, 149]. Following earlier microsatellite studies with N. vitripennis [96, 97] we calculated pairwise FST [150,151,152], Hedrick’s GʹST [153, 154], and p-values following G-statistics (as implemented in GenAlEx 6.502) with 999 permutations.

Statistical analyses

Statistical analyses were carried out using R [155]. Numeric data from olfactometer experiments (Additional file 1: Tables S1, S2), i.e. the allocation time above the test, control, or the two transition chambers were analysed using linear mixed models [156], or generalized linear mixed models [157] depending on distribution and homogeneity of variances. Multiple comparisons were performed with post-hoc Tukey’s test [158]. The number of strains from each of the habitats that significantly reacted to one of the two odours in the olfactometer experiments, and the distribution of individuals from the subpopulations between main and alternative microhabitat were analysed using Fisher Exact Test (Fig. 2). As numeric data from experiments onhost acceptance, i.e. numbers of parasitized pupae and offspring did not have a normal distribution, the differences between strains were analysed using generalized linear models (family “poisson”) [157]. Because this revealed no difference between strains within each group (nest or carrion), data between nest and carrion strains were compared using non-parametric Mann–Whitney U-test. The occurrence of copulations were compared using generalized linear models, family “binomial”, [157], followed by ANOVA and Tukey-tests for single comparisons [158]. F1-female offspring numbers from experiments on postmating barriers were compared by using generalized linear models, “quasipoisson” and “negative binomial”, followed by ANOVA and post-hoc Tukey’s tests [158].

Availability of data and materials

The datasets supporting the conclusions of this article will be made available in https://doi.org/10.5061/dryad.8kprr4xn3.

Change history

  • 14 May 2022

    The Funding section has been updated.

References

  1. Coyne JA, Orr HA. Speciation. Sunderland: Sinauer Associates; 2004.

    Google Scholar 

  2. Butlin R, Debelle A, Kerth C, Snook RR, Beukeboom LW, Castillo RC, et al. What do we need to know about speciation? Trends Ecol Evol. 2012;27:27–39.

    Article  PubMed  Google Scholar 

  3. Otto SP, Whitton J. Polyploid incidence and evolution. Annu Rev Genet. 2000;34:401–37.

    Article  CAS  PubMed  Google Scholar 

  4. Mallet J. Hybrid speciation. Nature. 2007;446:279–83.

    Article  CAS  PubMed  Google Scholar 

  5. Bird CE, Fernandez-Silva I, Skillings DJ, Toonen RJ. Sympatric speciation in the post “Modern Synthesis” era of evolutionary biology. Evol Biol. 2012;39:158–80.

    Article  Google Scholar 

  6. Via S. Sympatric speciation in animals: the ugly duckling grows up. Trends Ecol Evol. 2001;16:381–90.

    Article  CAS  PubMed  Google Scholar 

  7. Bolnick DI, Fitzpatrick BM. Sympatric speciation: models and empirical evidence. Annu Rev Ecol Evol Syst. 2007;38:459–87.

    Article  Google Scholar 

  8. Orr HA, Turelli M. The evolution of postzygotic isolation: accumulating Dobzhansky–Muller incompatibilities. Evolution. 2001;55:1085–94.

    CAS  PubMed  Google Scholar 

  9. Nosil P, Flaxman SM. Conditions for mutation-order speciation. Proc R Soc B Biol Sci. 2011;278:399–407.

    Article  Google Scholar 

  10. Malinsky M, Challis RJ, Tyers AM, Schiffels S, Terai Y, Ngatunga BP, et al. Genomic islands of speciation separate cichlid ecomorphs in an East African crater lake. Science. 2015;350:1493–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Kautt AF, Kratochwil CF, Nater A, Machado-Schiaffino G, Olave M, Henning F, et al. Contrasting signatures of genomic divergence during sympatric speciation. Nature. 2020;588:106–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Berlocher SH, Feder JL. Sympatric speciation in phytophagous insects: moving beyond controversy? Annu Rev Entomol. 2002;47:773–815.

    Article  CAS  PubMed  Google Scholar 

  13. Forbes AA, Devine SN, Hippee AC, Tvedte ES, Ward AKG, Widmayer HA, et al. Revisiting the particular role of host shifts in initiating insect speciation. Evolution. 2017;71:1126–37.

    Article  PubMed  Google Scholar 

  14. Drès M, Mallet J. Host races in plant–feeding insects and their importance in sympatric speciation. Philos Trans R Soc Lond B Biol Sci. 2002;357:471–92.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Higashi M, Takimoto G, Yamamura N. Sympatric speciation by sexual selection. Nature. 1999;402:523–6.

    Article  CAS  PubMed  Google Scholar 

  16. Gavrilets S, Waxman D. Sympatric speciation by sexual conflict. Proc Natl Acad Sci. 2002;99:10533–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Ritchie MG. Sexual selection and speciation. Annu Rev Ecol Evol Syst. 2007;38:79–102.

    Article  Google Scholar 

  18. Weissing FJ, Edelaar P, Van Doorn GS. Adaptive speciation theory: a conceptual review. Behav Ecol Sociobiol. 2011;65:461–80.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Safran RJ, Scordato ES, Symes LB, Rodríguez RL, Mendelson TC. Contributions of natural and sexual selection to the evolution of premating reproductive isolation: a research agenda. Trends Ecol Evol. 2013;28:643–50.

    Article  PubMed  Google Scholar 

  20. Scordato ESC, Symes LB, Mendelson TC, Safran RJ. The role of ecology in speciation by sexual selection: a systematic empirical review. J Hered. 2014;105:782–94.

    Article  PubMed  Google Scholar 

  21. Rundle HD, Rowe L. The contribution of sexual selection to ecological and mutation-order speciation. Evolution. 2018;72:2571–5.

    Article  PubMed  Google Scholar 

  22. Maan ME, Seehausen O. Ecology, sexual selection and speciation. Ecol Lett. 2011;14:591–602.

    Article  PubMed  Google Scholar 

  23. Austin A, Dowton M. The Hymenoptera: an introduction. In: Hymenoptera: evolution, biodiversity and biological control. Collingwood: CSIRO Publishing; 2000. p. 3–10.

  24. Aguiar AP, Deans AR, Engel MS, Forshage M, Huber JT, Jennings JT, et al. Order Hymenoptera. In: Zhang Z-Q, editor. Animal biodiversity: an outline of higher-level classification and survey of taxonomic richness (Addenda 2013). Zootaxa. 2013;3703:51–62.

  25. Davis RB, Baldauf SL, Mayhew PJ. The origins of species richness in the Hymenoptera: insights from a family-level supertree. BMC Evol Biol. 2010;10:109.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Peters RS, Krogmann L, Mayer C, Donath A, Gunkel S, Meusemann K, et al. Evolutionary history of the Hymenoptera. Curr Biol. 2017;27:1013–8.

    Article  CAS  PubMed  Google Scholar 

  27. Godfray HCJ. Parasitoids: behavioral and evolutionary ecology. Princeton: Princeton University Press; 1994.

    Book  Google Scholar 

  28. Quicke DLJ. Parasitic wasps. London: Chapman Hall; 1997.

    Google Scholar 

  29. Forbes AA, Bagley RK, Beer MA, Hippee AC, Widmayer HA. Quantifying the unquantifiable: why Hymenoptera, not Coleoptera, is the most speciose animal order. BMC Ecol. 2018;18:1–11.

    Article  Google Scholar 

  30. Gokhman V. Dimensions and borderlines of parasitoid Hymenoptera species: a paradigm shift? Biol Bull Rev. 2018;8:227–33.

    Article  Google Scholar 

  31. Gauld ID, Janzen DH. The systematics and biology of the Costa Rican species of parasitic wasps in the Thyreodon genus-group (Hymenoptera: Ichneumonidae). Zool J Linn Soc. 2004;141:297–351.

    Article  Google Scholar 

  32. Smith MA, Rodriguez JJ, Whitfield JB, Deans AR, Janzen DH, Hallwachs W, et al. Extreme diversity of tropical parasitoid wasps exposed by iterative integration of natural history, DNA barcoding, morphology, and collections. Proc Natl Acad Sci. 2008;105:12359–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Stigenberg J. Revision of the Western Palearctic Meteorini (Hymenoptera, Braconidae), with a molecular characterization of hidden Fennoscandian species diversity. Zootaxa. 2011;3084:1–95.

    Article  Google Scholar 

  34. Butcher BA, Smith MA, Sharkey MJ, Quicke DL. A turbo-taxonomic study of Thai Aleiodes (Aleiodes) and Aleiodes (Arcaleiodes) (Hymenoptera: Braconidae: Rogadinae) based largely on COI barcoded specimens, with rapid descriptions of 179 new species. Zootaxa. 2012;3457:1–232.

    Article  Google Scholar 

  35. Chen H, Talamas EJ, Valerio AA, Masner L, Johnson NF. Revision of the world species of the genus Chromoteleia Ashmead (Hymenoptera, Platygastridae, Scelioninae). ZooKeys. 2018;778:1–95.

    Article  Google Scholar 

  36. van Achterberg K, Schilthuizen M, van der Meer M, Delval R, Dias C, Hoynck M, et al. A new parasitoid wasp, Aphaereta vondelparkensis sp. n. (Braconidae, Alysiinae), from a city park in the centre of Amsterdam. Biodivers Data J. 2020;8:e49017.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Wiegmann BM, Mitter C, Farrell B. Diversification of carnivorous parasitic insects: extraordinary radiation or specialized dead end? Am Nat. 1993;142:737–54.

    Article  Google Scholar 

  38. Noyes JS. Universal chalcidoidea database. 2019. http://www.nhm.ac.uk/chalcidoids.

  39. Askew RR. Considerations on speciation in Chalcidoidea (Hymenoptera). Evolution. 1968;22:642–5.

    Article  CAS  PubMed  Google Scholar 

  40. Unruh T, Messing R. Intraspecific biodiversity in Hymenoptera: implications for conservation and biological control. In: Hymenoptera and biodiversity. Wallingford: CAB International; 1993. p. 27–52.

  41. Conner JK, Hartl DL. A primer of ecological genetics. Sunderland: Sinauer Associates Incorporated; 2004.

    Google Scholar 

  42. Peer K, Taborsky M. Outbreeding depression, but no inbreeding depression in haplodiploid ambrosia beetles with regular sibling mating. Evolution. 2005;59:317–23.

    PubMed  Google Scholar 

  43. Zayed A, Packer L. Complementary sex determination substantially increases extinction proneness of haplodiploid populations. Proc Natl Acad Sci. 2005;102:10742–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Asplen M, Whitfield J, De Boer J, Heimpel GE. Ancestral state reconstruction analysis of hymenopteran sex determination mechanisms. J Evol Biol. 2009;22:1762–9.

    Article  CAS  PubMed  Google Scholar 

  45. Werren JH. The evolution of inbreeding in haplodiploid organisms. In: The natural history of inbreeding and outbreeding. Chicago: The University of Chicago Press; 1993. p. 42–59.

  46. Gordh G, DeBach P. Courtship behavior in the Aphytis lingnanensis group, its potential usefulness in taxonomy, and a review of sexual behavior in the parasitic Hymenoptera (Chalcidoidea: Aphelinidae). Hilgardia. 1978;46:37–75.

    Article  Google Scholar 

  47. Ridley M. Clutch size and mating frequency in parasitic Hymenoptera. Am Nat. 1993;142:893–910.

    Article  Google Scholar 

  48. Charnov EL. Sex ratio evolution in a variable environment. Nature. 1981;289:27–33.

    Article  CAS  PubMed  Google Scholar 

  49. Burton-Chellew MN, Koevoets T, Grillenberger BK, Sykes EM, Underwood SL, Bijlsma K, et al. Facultative sex ratio adjustment in natural populations of wasps: cues of local mate competition and the precision of adaptation. Am Nat. 2008;172:393–404.

    Article  PubMed  Google Scholar 

  50. Machado CA, Robbins N, Gilbert MTP, Herre EA. Critical review of host specificity and its coevolutionary implications in the fig/fig-wasp mutualism. Proc Natl Acad Sci. 2005;102(suppl 1):6558–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Sha Z, Zhu C, Murphy R, Huang D. Diglyphus isaea (Hymenoptera: Eulophidae): a probable complex of cryptic species that forms an important biological control agent of agromyzid leaf miners. J Zool Syst Evol Res. 2007;45:128–35.

    Article  Google Scholar 

  52. König K, Zundel P, Krimmer E, König C, Pollmann M, Gottlieb Y, et al. Reproductive isolation due to prezygotic isolation and postzygotic cytoplasmic incompatibility in parasitoid wasps. Ecol Evol. 2019. https://doi.org/10.1002/ece3.5588.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Forbes AA, Powell THQ, Stelinski LL, Smith JJ, Feder JL. Sequential sympatric speciation across trophic levels. Science. 2009;323:776–9.

    Article  CAS  PubMed  Google Scholar 

  54. Feder JL, Forbes AA. Sequential speciation and the diversity of parasitic insects. Ecol Entomol. 2010;35:67–76.

    Article  Google Scholar 

  55. Hood GR, Forbes AA, Powell THQ, Egan SP, Hamerlinck G, Smith JJ, et al. Sequential divergence and the multiplicative origin of community diversity. Proc Natl Acad Sci. 2015;112:E5980–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Hamerlinck G, Hulbert D, Hood GR, Smith JJ, Forbes AA. Histories of host shifts and cospeciation among free-living parasitoids of Rhagoletis flies. J Evol Biol. 2016;29:1766–79.

    Article  CAS  PubMed  Google Scholar 

  57. König K, Krimmer E, Brose S, Gantert C, Buschlüter I, König C, et al. Does early learning drive ecological divergence during speciation processes in parasitoid wasps? Proc R Soc Lond B Biol Sci. 2015;282:20141850.

    Google Scholar 

  58. König K, Seeger L, Steidle JL. Sexy mouth odour? Male oral gland pheromone in the grain beetle parasitoid Lariophagus distinguendus (Förster) (Hymenoptera: Pteromalidae). BioMed Res Int. 2015. https://doi.org/10.1155/2015/216952.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Breeuwer JAJ, Werrren JH. Hybrid breakdown between two haplodiploid species: the role of nuclear and cytoplasmic genes. Evolution. 1995;49:705–17.

    Article  PubMed  Google Scholar 

  60. Bordenstein SR, O’Hara FP, Werren JH. Wolbachia-induced incompatibility precedes other hybrid incompatibilities in Nasonia. Nature. 2001;409:707–10.

    Article  CAS  PubMed  Google Scholar 

  61. Clark ME, O’Hara FP, Chawla A, Werren JH. Behavioral and spermatogenic hybrid male breakdown in Nasonia. Heredity. 2010;104:289–301.

    Article  CAS  PubMed  Google Scholar 

  62. Koevoets T, Niehuis O, Van De Zande L, Beukeboom LW. Hybrid incompatibilities in the parasitic wasp genus Nasonia: negative effects of hemizygosity and the identification of transmission ratio distortion loci. Heredity. 2012;108:302–11. https://doi.org/10.1038/hdy.2011.75.

    Article  CAS  PubMed  Google Scholar 

  63. Giesbers MCWG, Gerritsma S, Buellesbach J, Diao W, Pannebakker BA, van de Zande L, et al. Prezygotic isolation in the parasitoid wasp genus Nasonia. In: Speciation: natural processes, genetics and biodiversity. Nova Science Publishers; 2013. p. 165–91. https://www.scopus.com/inward/record.uri?eid=2-s2.0-84892127414&partnerID=40&md5=bc825cd613098b46e716e6eb6f9018d0.

  64. Niehuis O, Buellesbach J, Gibson JD, Pothmann D, Hanner C, Mutti NS, et al. Behavioural and genetic analyses of Nasonia shed light on the evolution of sex pheromones. Nature. 2013;494:345–8.

    Article  CAS  PubMed  Google Scholar 

  65. Beukeboom LW, Koevoets T, Morales HE, Ferber S, van de Zande L. Hybrid incompatibilities are affected by dominance and dosage in the haplodiploid wasp Nasonia. Front Genet. 2015. https://doi.org/10.3389/fgene.2015.00140.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Werren JH, Richards S, Desjardins CA, Niehuis O, Gadau J, Colbourne JK, et al. Functional and evolutionary insights from the genomes of three parasitoid Nasonia species. Science. 2010;327:343–8. https://doi.org/10.1126/science.1178028.

    Article  CAS  PubMed  Google Scholar 

  67. Loehlin DW, Oliveira DCSG, Edwards R, Giebel JD, Clark ME, Cattani MV, et al. Non-coding changes cause sex-specific wing size differences between closely related species of Nasonia. PLoS Genet. 2010;6:e1000821.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  68. Bertossa RC, Van De Zande L, Beukeboom LW, Beersma DGM. Phylogeny and oscillating expression of period and cryptochrome in short and long photoperiods suggest a conserved function in Nasonia vitripennis. Chronobiol Int. 2014;31:749–60.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Pegoraro M, Bafna A, Davies NJ, Shuker DM, Tauber E. DNA methylation changes induced by long and short photoperiods in Nasonia. Genome Res. 2016;26:203–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Mair MM, Kmezic V, Huber S, Pannebakker BA, Ruther J. The chemical basis of mate recognition in two parasitoid wasp species of the genus Nasonia. Entomol Exp Appl. 2017;164:1–15.

    Article  CAS  Google Scholar 

  71. Geuverink E, Rensink AH, Rondeel I, Beukeboom LW, van de Zande L, Verhulst EC. Maternal provision of transformer-2 is required for female development and embryo viability in the wasp Nasonia vitripennis. Insect Biochem Mol Biol. 2017;90 Supplement C:23–33.

    Article  CAS  Google Scholar 

  72. Martinson EO, Mrinalini, Kelkar YD, Chang C-H, Werren JH. The evolution of venom by co-option of single-copy genes. Curr Biol. 2017;27:2007-2013.e8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Schurmann D, Sommer C, Schinko AP, Greschista M, Smid H, Steidle JL. Demonstration of long-term memory in the parasitic wasp Nasonia vitripennis. Entomol Exp Appl. 2012;143:199–206.

    Article  Google Scholar 

  74. Schurmann D, Kugel D, Steidle JL. Early memory in the parasitoid wasp Nasonia vitripennis. J Comp Physiol A. 2015;201:375–83.

    Article  Google Scholar 

  75. Breeuwer JAJ, Werren JH. Microorganisms associated with chromosome destruction and reproductive isolation between two insect species. Nat Lond. 1990;346:558–60.

    Article  CAS  Google Scholar 

  76. Bordenstein SR, Werren JH. Bidirectional incompatibility among divergent Wolbachia and incompatibility level differences among closely related Wolbachia in Nasonia. Heredity. 2007;99:278–87.

    Article  CAS  PubMed  Google Scholar 

  77. Koevoets T, Beukeboom LW. Genetics of postzygotic isolation and Haldane’s rule in haplodiploids. Heredity. 2009;102:16–23.

    Article  CAS  PubMed  Google Scholar 

  78. Raychoudhury R, Desjardins CA, Buellesbach J, Loehlin DW, Grillenberger BK, Beukeboom L, et al. Behavioral and genetic characteristics of a new species of Nasonia. Heredity. 2010;104:278–88.

    Article  CAS  PubMed  Google Scholar 

  79. Ellison CK, Niehuis O, Gadau J. Hybrid breakdown and mitochondrial dysfunction in hybrids of Nasonia parasitoid wasps. J Evol Biol. 2008;21:1844–51.

    Article  CAS  PubMed  Google Scholar 

  80. Dittmer J, van Opstal EJ, Shropshire JD, Bordenstein SR, Hurst GD, Brucker RM. Disentangling a holobiont–recent advances and perspectives in Nasonia wasps. Front Microbiol. 2016;7:1478.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Werren JH, Baldo L, Clark ME. Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol. 2008;6:741–51.

    Article  CAS  PubMed  Google Scholar 

  82. Werren JH. Biology of Wolbachia. Annu Rev Entomol. 1997;42:587–609.

    Article  CAS  PubMed  Google Scholar 

  83. Hunter MS, Perlman SJ, Kelly SE. A bacterial symbiont in the Bacteroidetes induces cytoplasmic incompatibility in the parasitoid wasp Encarsia pergandiella. Proc R Soc Lond B Biol Sci. 2003;270:2185–90.

    Article  Google Scholar 

  84. Takano S, Tuda M, Takasu K, Furuya N, Imamura Y, Kim S, et al. Unique clade of alphaproteobacterial endosymbionts induces complete cytoplasmic incompatibility in the coconut beetle. Proc Natl Acad Sci. 2017;114:6110–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Rosenwald LC, Sitvarin MI, White JA. Endosymbiotic Rickettsiella causes cytoplasmic incompatibility in a spider host. Proc R Soc B. 2020;287:20201107.

    Article  PubMed  PubMed Central  Google Scholar 

  86. Vavre F, Fleury F, Varaldi J, Fouillet P, Bouletreau M. Infection polymorphism and cytoplasmic incompatibility in Hymenoptera-Wolbachia associations. Heredity. 2002;88:361–5.

    Article  CAS  PubMed  Google Scholar 

  87. Breeuwer JA, Werren JH. Cytoplasmic incompatibility and bacterial density in Nasonia vitripennis. Genetics. 1993;135:565–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Schröder H. Crossing experiments and molecularbiological researches at two ecotypes of Nasonia vitripennis (Hymenoptera, Pteromalidae). Mitteilungen Dtsch Ges Für Allg Angew Entomol. 2000;12:165–8.

    Google Scholar 

  89. Abraham R, Peters RS. Nistkästen als Lebensraum für Insekten, besonders Fliegen und ihre Schlupfwespen, Nestboxes as habitat for insects, especially for flies and their parasitoids. Vogelwarte. 2008;46:195–205.

    Google Scholar 

  90. Schurmann D, Collatz J, Hagenbucher S, Ruther J, Steidle JLM. Olfactory host finding, intermediate memory and its potential ecological adaptation in Nasonia vitripennis. Naturwissenschaften. 2009;96:383–91.

    Article  CAS  PubMed  Google Scholar 

  91. Peters RS, Abraham R. The food web of parasitoid wasps and their non-phytophagous fly hosts in birds’ nests (Hymenoptera: Chalcidoidea, and Diptera: Cyclorrhapha). J Nat Hist. 2010;44:625–38.

    Article  Google Scholar 

  92. Garrido-Bautista J, Moreno-Rueda G, Baz A, Canal D, Camacho C, Cifrián B, et al. Variation in parasitoidism of Protocalliphora azurea (Diptera: Calliphoridae) by Nasonia vitripennis (Hymenoptera: Pteromalidae) in Spain. Parasitol Res. 2020;119:559–66.

    Article  PubMed  Google Scholar 

  93. Grassberger M, Frank C. Initial study of arthropod succession on pig carrion in a Central European urban habitat. J Med Entomol. 2004;41:511–23.

    Article  CAS  PubMed  Google Scholar 

  94. Ruther J, McCaw J, Böcher L, Pothmann D, Putz I. Pheromone diversification and age-dependent behavioural plasticity decrease interspecific mating costs in Nasonia. PLoS ONE. 2014;9:e89214–e89214.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  95. Whiting AR. The biology of the parasitic wasp Mormoniella vitripennis [= Nasonia brevicornis](Walker). Q Rev Biol. 1967;42:333–406.

    Article  Google Scholar 

  96. Grillenberger BK, Gadau J, Bijlsma R, Van De Zande L, Beukeboom LW. Female dispersal and isolation-by-distance of Nasonia vitripennis populations in a local mate competition context. Entomol Exp Appl. 2009;132:147–54.

    Article  Google Scholar 

  97. Grillenberger B, Koevoets T, Burton-Chellew M, Sykes E, Shuker D, Van de Zande L, et al. Genetic structure of natural Nasonia vitripennis populations: validating assumptions of sex-ratio theory. Mol Ecol. 2008;17:2854–64.

    Article  CAS  PubMed  Google Scholar 

  98. Van der Merwe J. Investigations on the biology and ecology of Mormoniella vitripennis Walk. (Pteromalidae, Hym.). J Entomol Soc S Afr. 1943;6:48–64.

    Google Scholar 

  99. Smith G, Pimentel D. The effect of two host species on the longevity and fertility of Nasonia vitripennis. Ann Entomol Soc Am. 1969;62:305–8.

    Article  Google Scholar 

  100. Fitzpatrick BM, Fordyce JA, Gavrilets S. What, if anything, is sympatric speciation? J Evol Biol. 2008;21:1452–9.

    Article  CAS  PubMed  Google Scholar 

  101. Vinson SB. The general host selection behavior of parasitoid Hymenoptera and a comparison of initial strategies utilized by larvaphagous and oophagous Species. Biol Control. 1998;11:79–96.

    Article  Google Scholar 

  102. Steidle JL, van Loon JJ. Chemoecology of parasitoid and predator oviposition behaviour. In: Chemoecology of insect eggs and egg deposition. Berlin: Blackwell; 2002. p. 291–317.

  103. Turissini D, Liu G, David J, Matute D. The evolution of reproductive isolation in the Drosophila yakuba complex of species. J Evol Biol. 2015;28:557–75.

    Article  CAS  PubMed  Google Scholar 

  104. Tadeo E, Feder JL, Egan SP, Schuler H, Aluja M, Rull J. Divergence and evolution of reproductive barriers among three allopatric populations of Rhagoletis cingulata across eastern North America and Mexico. Entomol Exp Appl. 2015;156:301–11.

    Article  Google Scholar 

  105. Kozak GM, Rudolph AB, Colon BL, Fuller RC. Postzygotic isolation evolves before prezygotic isolation between fresh and saltwater populations of the Rainwater Killifish, Lucania parva. Int J Evol Biol. 2012. https://doi.org/10.1155/2012/523967.

    Article  PubMed  PubMed Central  Google Scholar 

  106. Chin TA, Cáceres CE, Cristescu ME. The evolution of reproductive isolation in Daphnia. BMC Evol Biol. 2019;19:216.

    Article  PubMed  PubMed Central  Google Scholar 

  107. Schröder H. Okologie und intraspezifische Variation bei Nasonia vitripennis Walker (Hymenoptera, Pteromalidae). Mitteilungen Dtsch Ges Fuer Allg Angew Entomol. 1995;10:521–4.

    Google Scholar 

  108. Schröder H, Abraham R. Ecotypes or two species within Nasonia vitripennis (Chalcidoidea, Pteromalidae)? Mitteilungen Dtsch Ges Fuer Allg Angew Entomol. 1997;11:789–92.

    Google Scholar 

  109. Møller A, Flensted-Jensen E, Klarborg K, Mardal W, Nielsen J. Climate change affects the duration of the reproductive season in birds. J Anim Ecol. 2010;79:777–84.

    PubMed  Google Scholar 

  110. Lewis Z, Lizé A. Insect behaviour and the microbiome. Curr Opin Insect Sci. 2015;9:86–90.

    Article  PubMed  Google Scholar 

  111. Holdaway FG. Field populations and natural control of Lucilia sericata. Nature. 1930;126:648–9.

    Article  Google Scholar 

  112. Bush GL. Host race formation and sympatric speciation in Rhagoletis fruit flies (Diptera: Tephritidae). Psyche (Stuttg). 1992;99:335–57.

    Article  Google Scholar 

  113. Jaenike J. Criteria for ascertaining the existence of host races. Am Nat. 1981;117:830–4.

    Article  Google Scholar 

  114. Ruther J, Hammerl T. An oral male courtship pheromone terminates the response of Nasonia vitripennis females to the male-produced sex attractant. J Chem Ecol. 2014;40:56–62.

    Article  CAS  PubMed  Google Scholar 

  115. Jennings JH, Snook RR, Hoikkala A. Reproductive isolation among allopatric Drosophila montana populations. Evolution. 2014;68:3095–108.

    Article  PubMed  Google Scholar 

  116. Turissini DA, McGirr JA, Patel SS, David JR, Matute DR. The Rate of evolution of postmating-prezygotic reproductive isolation in Drosophila. Mol Biol Evol. 2017;35:312–34.

    Article  PubMed Central  CAS  Google Scholar 

  117. Garlovsky MD, Snook RR. Persistent postmating, prezygotic reproductive isolation between populations. Ecol Evol. 2018;8:9062–73.

    Article  PubMed  PubMed Central  Google Scholar 

  118. Keller LF, Waller DM. Inbreeding effects in wild populations. Trends Ecol Evol. 2002;17:230–41.

    Article  Google Scholar 

  119. Antolin MF. A genetic perspective on mating systems and sex ratios of parasitoid wasps. Popul Ecol. 1999;41:29–37.

    Article  Google Scholar 

  120. Waser NM, Williams CF. Inbreeding and outbreeding. In: Evolutionary ecology: concepts and case studies. 2001; p. 84–96.

  121. Edmands S. Does parental divergence predict reproductive compatibility? Trends Ecol Evol. 2002;17:520–7.

    Article  Google Scholar 

  122. Escobar JS, Nicot A, David P. The different sources of variation in inbreeding depression, heterosis and outbreeding depression in a metapopulation of Physa acuta. Genetics. 2008;180:1593–608.

    Article  PubMed  PubMed Central  Google Scholar 

  123. Gimond C, Jovelin R, Han S, Ferrari C, Cutter AD, Braendle C. Outbreeding depression with low genetic variation in selfing Caenorhabditis nematodes. Evolution. 2013;67:3087–101.

    Article  PubMed  Google Scholar 

  124. Atalay D, Schausberger P. Balancing in-and out-breeding by the predatory mite Phytoseiulus persimilis. Exp Appl Acarol. 2018;74:159–69.

    Article  PubMed  PubMed Central  Google Scholar 

  125. Molbo D, Parker E. Mating structure and sex ratio variation in a natural population of Nasonia vitripennis. Proc R Soc Lond B Biol Sci. 1996;263:1703–9.

    Article  Google Scholar 

  126. Werren JH, Loehlin DW. The parasitoid wasp Nasonia: an emerging model system with haploid male genetics. Cold Spring Harb Protoc. 2009;2009:pdb.emo134.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  127. Johannesson K, Panova M, Kemppainen P, André C, Rolan-Alvarez E, Butlin RK. Repeated evolution of reproductive isolation in a marine snail: unveiling mechanisms of speciation. Philos Trans R Soc B Biol Sci. 2010;365:1735–47.

    Article  Google Scholar 

  128. Rundle HD, Nosil P. Ecological speciation. Ecol Lett. 2005;8:336–52.

    Article  Google Scholar 

  129. Bolnick DI, Barrett RD, Oke KB, Rennison DJ, Stuart YE. (Non) parallel evolution. Annu Rev Ecol Evol Syst. 2018;49:303–30.

    Article  Google Scholar 

  130. Nosil P. Ecological speciation. Oxford: Oxford University Press; 2012.

    Book  Google Scholar 

  131. Seehausen O, Butlin RK, Keller I, Wagner CE, Boughman JW, Hohenlohe PA, et al. Genomics and the origin of species. Nat Rev Genet. 2014;15:176–92.

    Article  CAS  PubMed  Google Scholar 

  132. Buschinger A. Evolution, speciation, and inbreeding in the parasitic ant genus Epimyrma (Hymenoptera, Formicidae). J Evol Biol. 1989;2:265–83.

    Article  Google Scholar 

  133. Stegniy V. Hard inbreeding under extreme environmental conditions is the most important factor of microevolution and speciation. Russ J Genet. 2017;53:757–65.

    Article  CAS  Google Scholar 

  134. Olsen KC, Ryan WH, Winn AA, Kosman ET, Moscoso JA, Krueger-Hadfield SA, et al. Inbreeding shapes the evolution of marine invertebrates. Evolution. 2020;74:871–82.

    Article  PubMed  PubMed Central  Google Scholar 

  135. Bouček Z, Rasplus JY. Illustrated key to west-Palearctic genera of Pteromalidae (Hymenoptera: Chalcidoidea). Paris: INRA Editions; 1991.

    Google Scholar 

  136. Pannebakker BA, Cook N, Van Den Heuvel J, Van De Zande L, Shuker DM. Genomics of sex allocation in the parasitoid wasp Nasonia vitripennis. BMC Genomics. 2020;21:1–14.

    Article  CAS  Google Scholar 

  137. Steidle JLM, Schoeller M. Olfactory host location and learning in the granary weevil parasitoid Lariophagus distinguendus (Hymenoptera: Pteromalidae). J Insect Behav. 1997;10:331–42.

    Article  Google Scholar 

  138. Steinberg S, Dicke M, Vet LEM, Wanningen R. Response of the braconid parasitoid Cotesia (=Apanteles) glomerata to volatile infochemicals: effects of bioassay set-up, parasitoid age and experience and barometric flux. Entomol Exp Appl. 1992;63:163–75.

    Article  CAS  Google Scholar 

  139. Sobel JM, Chen GF. Unification of methods for estimating the strength of reproductive isolation. Evolution. 2014;68:1511–22.

    Article  PubMed  Google Scholar 

  140. Ramsey J, Bradshaw HD, Schemske DW. Components of reproductive isolation between the monkeyflowers Mimulus lewisii and M. cardinalis (Phyrmaceae). Evolution. 2003;57:1520–34.

    PubMed  Google Scholar 

  141. van de Zande L, Ferber S, de Haan A, Beukeboom LW, van Heerwaarden J, Pannebakker BA. Development of a Nasonia vitripennis outbred laboratory population for genetic analysis. Mol Ecol Resour. 2014;14:578–87.

    Article  PubMed  Google Scholar 

  142. Flores-Rentería L, Krohn A. Scoring microsatellite loci. In: Microsatellites. Springer; 2013. p. 319–36.

  143. Van Oosterhout C, Hutchinson WF, Wills DP, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.

    Article  CAS  Google Scholar 

  144. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  145. Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015;15:1179–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  146. Falush D, Stephens M, Pritchard J. Inference of population genetic structure: extensions to linked loci and correlated alíele frequencies. Genetics. 2003;164:1567–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  147. Hubisz MJ, Falush D, Stephens M, Pritchard JK. Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009;9:1322–32.

    Article  PubMed  PubMed Central  Google Scholar 

  148. Peakall R, Smouse PE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.

    Article  Google Scholar 

  149. Smouse PE, Peakall R. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28:2537–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  150. Hartl DL, Clark AG, Clark AG. Principles of population genetics. Sunderland: Sinauer Associates; 1997.

    Google Scholar 

  151. Frankham R, Briscoe DA, Ballou JD. Introduction to conservation genetics. Cambridge: Cambridge University Press; 2002.

    Book  Google Scholar 

  152. Frankham R, Ballou JD, Briscoe DA. A primer of conservation genetics. Cambridge: Cambridge University Press; 2004.

    Book  Google Scholar 

  153. Hedrick PW. A standardized genetic differentiation measure. Evolution. 2005;59:1633–8.

    CAS  PubMed  Google Scholar 

  154. Meirmans PG, Hedrick PW. Assessing population structure: FST and related measures. Mol Ecol Resour. 2011;11:5–18.

    Article  PubMed  Google Scholar 

  155. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016. https://www.R-project.org/.

  156. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest Package: tests in linear mixed effects models. J Stat Softw. 2017;82:1–26.

    Article  Google Scholar 

  157. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2014;67:1–48.

    Google Scholar 

  158. Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biom J. 2008;50:346–63.

    Article  PubMed  Google Scholar 

  159. Pritchard JK, Wena X, Falushb D. Documentation for structure software: Version 2.3. 2009. http://pritch.bsd.uchicago.edu/structure.html.

Download references

Acknowledgements

The authors would like to thank Dick Askew and two unknown reviewers for very helpful comments and improving the English of the manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL. PM was financially supported by the Landesgraduiertenförderung Baden-Württemberg.

Author information

Authors and Affiliations

Authors

Contributions

PM contributed to the design of the study, conducted experiments, performed data analyses and drafted the manuscript. MP and VP performed the olfactometer experiments. CR und RR performed host acceptance studies. JW, RB and DM conducted the experiments on postmating barriers. MF contributed to the analysis of the microsatellite data. MH collected the field data. CK contributed to the molecular analyses. LvB analysed olfactometer data. JLMS conceived the concept of the study together with PM, performed statistical analyses, supervised the experiments and finalised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Johannes L. M. Steidle.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Test statistics from olfactometer experiments with strains of Lariophagus distinguendus from bird nests and carrion tested on the odour of bird nests or carrion. The odour field contained samples of bird nests or carrion. Strain abbreviations with “N” refer to strains collected in bird nests, strains with “A” were collected next to carrions. Different letters within the same row indicate significant differences in allocation time between fields (p < 0.005). Single comparisons were made using Tukey-test based on linear mixed models or generalized mixed models with field as factor and observation as random factor (Table S2, separate file). Lines with significant differences between odour field and control field 2 are shaded. Table S2. Test statistics from olfactometer experiments with strains of Nasonia vitripennis from bird nests and carrion tested on the odour of bird nests or carrion. Strain abbreviations with “N” refer to strains collected in bird nests, strains with “A” were collected next to carrions. Table S3. Distribution of individuals between main and alternative microhabitat in the different subpopulations identified by microsatellites for k = 3 and k = 5 subpopulations. The main habitats are bird nests for subpopulation 1 and carrion for subpopulations 2–5. Table S4. Mean number (± S.D.) of F1-female offspring from couples consisting of females and males from different strains. Data from intra-strain couples are shaded.

Additional file 2: Fig. S1.

Box and whisker plot of pupae of Lucilia sericata carrion flies parasitized by Nasonia vitripennis wasps. N2, N3, N9: Wasp strains that were collected in bird nests (yellow, n = 30 per strain). A1, A7, A19 Wasp strains that were collected in next to carrion (blue, n = 30 per strain). The plots show minimum, maximum, 1st and 3rd quartile, median, outliers as circles and mean as asterisk. Fig. S2. Box and whisker plot of offspring emerging from pupae of Lucilia sericata carrion flies parasitized by Nasonia vitripennis wasps. N2, N3, N9: Wasp strains that were collected in bird nests (yellow, n = 30 per strain). A1, A7, A19 Wasp strains that were collected in next to carrion (blue, n = 30 per strain). The plots show minimum, maximum, 1st and 3rd quartile, median, outliers as circles and mean as asterisk. Fig. S3. Putative population structure of the parasitoid N. vitripennis in the Hohenheim Park (Germany) based on Delta K. Results from three microsatellite analyses (Run 1-Run 3) using STRUCTURE v2.3.4 [159], and the online tool CLUMPAK.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Malec, P., Weber, J., Böhmer, R. et al. The emergence of ecotypes in a parasitoid wasp: a case of incipient sympatric speciation in Hymenoptera?. BMC Ecol Evo 21, 204 (2021). https://doi.org/10.1186/s12862-021-01938-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-021-01938-y

Keywords