Genetic and morphological variation in sexual and asexual parasitoids of the genus Lysiphlebus – an apparent link between wing shape and reproductive mode
© Petrovic et al.; licensee BioMed Central. 2015
Received: 1 September 2014
Accepted: 22 January 2015
Published: 4 February 2015
Morphological divergence often increases with phylogenetic distance, thus making morphology taxonomically informative. However, transitions to asexual reproduction may complicate this relationship because asexual lineages capture and freeze parts of the phenotypic variation of the sexual populations from which they derive. Parasitoid wasps belonging to the genus Lysiphlebus Foerster (Hymenoptera: Braconidae: Aphidiinae) are composed of over 20 species that exploit over a hundred species of aphid hosts, including many important agricultural pests. Within Lysiphlebus, two genetically and morphologically well-defined species groups are recognised: the “fabarum” and the “testaceipes” groups. Yet within each group, sexual as well as asexual lineages occur, and in L. fabarum different morphs of unknown origin and status have been recognised. In this study, we selected a broad sample of specimens from the genus Lysiphlebus to explore the relationship between genetic divergence, reproductive mode and morphological variation in wing size and shape (quantified by geometric morphometrics).
The analyses of mitochondrial and nuclear gene sequences revealed a clear separation between the “testaceipes” and “fabarum” groups of Lysiphlebus, as well as three well-defined phylogenetic lineages within the “fabarum” species group and two lineages within the “testaceipes” group. Divergence in wing shape was concordant with the deep split between the “testaceipes” and “fabarum” species groups, but within groups no clear association between genetic divergence and wing shape variation was observed. On the other hand, we found significant and consistent differences in the shape of the wing between sexual and asexual lineages, even when they were closely related.
Mapping wing shape data onto an independently derived molecular phylogeny of Lysiphlebus revealed an association between genetic and morphological divergence only for the deepest phylogenetic split. In more recently diverged taxa, much of the variation in wing shape was explained by differences between sexual and asexual lineages, suggesting a mechanistic link between wing shape and reproductive mode in these parasitoid wasps.
The morphological diversity of the living world, including the variation in size and shape, is generally assumed to be adaptive and shaped by natural selection. However, numerous internal factors, such as shared developmental systems inherited from common ancestors as well as structural and functional constraints, can restrict the response to selection [1-3]. An independently derived phylogeny is required to disentangle the effects of these influences on morphological evolution [4,5]. Although numerous studies have investigated morphological evolution in a phylogenetic context in various groups [6-9], this approach is still relatively rare in the studies of morphological variation in parasitoid wasps, even though they represent a very interesting biological model. Parasitoids are hyperdiverse (together with their hosts and host plants representing half of the world’s biodiversity [10,11]) and they possess complex life cycles and morphologies shaped by intimate host-parasitoid interactions. The frequent occurrence of both sexual and asexual reproduction poses a particular challenge because asexual lineages can freeze and amplify particular portions of morphological variation.
We addressed these issues in aphid parasitoids belonging to the genus Lysiphlebus Foerster (Braconidae: Aphidiinae). They are solitary koinobiont parasitoids, meaning that a single parasitoid larva develops per aphid host and allows the host to continue its growth and development while feeding upon it. Lysiphlebus exploits over a hundred species of mostly small aphids, including many important pests (e.g., soybean aphid, Aphis glycines Matsumura; black bean aphid, A. fabae Scopoli; cotton aphid, A. gossypii Glover) [12,13].
The polygenic basis of wing shape and its adaptive significance (flight, host finding), combined with the relative simplicity of a two-dimensional structure that allows a precise recording and analysis of the shape using geometric morphometrics, makes wing shape a very useful trait in the study of morphological evolution [7,26-28]. Here, we explored the molecular variation (based on nuclear and mitochondrial genes) and morphological variation (based on wing shape) to analyse possible factors, including reproductive mode, that are associated with the evolution of wing shape in aphid parasitoids from the genus Lysiphlebus. Specifically, we addressed the following questions: i) Is genetic differentiation accompanied by parallel changes in the shape of the wing? ii) Are parasitoid wasps assigned to different morphs and/or with different reproductive modes also divergent in their wing shape?
Analyses of mitochondrial sequences supported the separation of the genus into the “fabarum” (clades A, B, C; Figure 2) and the “testaceipes” (clades D, E; Figure 2) species groups, with a mean genetic distance of 7.6% between the groups.
Within the “testaceipes” group, L. orientalis (clade E) and L. testaceipes (clade D) split as separate taxa, with a mean genetic distance of 3.7%. L. testaceipes consisted of two distinct haplotypes with 1% genetic divergence, whereas four haplotypes were recorded in L. orientalis (0.3% mean genetic distance).
Based on the COI gene analysis, the “fabarum” group was composed of 16 distinct haplotypes that are clustered into three well-supported phylogenetic lineages (A, B, C in Figure 2). Overall, the mean genetic distance among specimens within the “fabarum” group was 2.2%.
Most of the haplotypes (12) were grouped into phylogenetic clade A. Within this clade, the tree topology had poor statistical support, and the highest genetic distance between two haplotypes was 1.1% (mean distance was 0.6%). Two haplotypes were found in both asexual and sexual parasitoids (L. fabarum 7 and L. fabarum 8), two were found only in sexual parasitoids, and eight were found only in parasitoids with asexual reproduction (Figure 2, Additional file 1).
Associations of COI haplotypes with the three morphs of L. fabarum were inconsistent. Most haplotypes were associated with one morph, but L. fabarum 1 was associated with the fabarum and cardui morphs, whereas the most common haplotype, L. fabarum 7, was associated with the fabarum and confusus morphs, but also with L. melandriicola.
Clade B was represented only with one haplotype, which is associated with all L. hirticornis specimens. Lysiphlebus hirticornis is clearly separated from clade A with a mean genetic distance of 4.1%.
All cardui morphs with sexual reproduction are grouped in phylogenetic clade C and separated from all other species/morphs. There is a genetic distance of 4.7% between these morphs and the same morphs with an asexual mode of reproduction that are all grouped in clade A. This provides strong evidence that these sexual parasitoids represent a yet undescribed species. The phylogenetic trees showed that the closest relatives of sexual cardui morphs in our material are the specimens belonging to L. hirticornis (2% divergence from L. hirticornis 1 haplotype).
The analysis of nuclear sequences (28S D2 gene) showed a well-supported split (100/99/99, Additional file 2) between outgroups and Lysiphlebus parasitoids, with a sequence divergence ranging from 11.2 to 26.7%.
All Lysiphlebus specimens clustered into two groups corresponding to the “fabarum” group (L. fabarum 1–4, L. cardui 1–2, L. confusus, and L. hirticornis 1 haplotypes) and the “testaceipes” group (L. orientalis 1, 2 and L. testaceipes 1 haplotype) (Additional file 2). The mean genetic divergence of the nuclear 28S D2 sequences between the “fabarum” and “testaceipes” species group was 3.1%.
Within the two species groups, the phylogenetic tree based on 28S D2 was poorly resolved due to insufficient sequence variation (average difference within species groups was only 0.1%).
Mean and standard deviation of wing size (CS) per phylogenetic lineage
mean size ± StDev
1336.36 ± 102.35
1073.80 ± 61.51
1252.34 ± 118.24
1171.66 ± 127.08
1276.62 ± 64.68
The permutation test against the null hypothesis of no phylogenetic signal revealed no evidence of a phylogenetic signal in wing size (tree length = 0.0231, P = 0.674) but did reveal a significant phylogenetic signal in wing shape (tree length = 0.0094, P = 0.037). The superimposition of the phylogenetic tree in the morphospace of the first two PC axes shows that concordance in divergence between genetic and morphological variation is due to the divergence between the two main clades (Figure 3).
Separate PCAs of wing shape variation within the “fabarum” (ABC) and “testaceipes” (DE) groups reveal that the divergence in wing shape is largely associated with reproductive mode. Within clade A, four subgroups can be distinguished a priori, based on the combination of morph (see Figure 1) and reproductive mode: A1 – L. cardui/asexual, A2 – L. confusus/asexual, A3 – L. fabarum/asexual and A4 – L. fabarum/sexual). Clades B (L. hirticornis/sexual) and C (L. cardui/sexual) are uniform with respect to these criteria.
Differences in wing size (below diagonal) and wing shape (above diagonal) between a priori groups based on morph and/or reproductive mode
Mean size ± StDev
A1 L. cardui/asexual
1360.1 ± 110.6
A2 L. confusus/asexual
1396.3 ± 32.0
A3 L. fabarum/asexual
1347.1 ± 97.2
A4 L. fabarum/sexual
1273.8 ± 86.1
B L. hirticornis/sexual
1073.8 ± 61.5
C L. cardui/sexual
1252.3 ± 118.2
Correct classification of individuals into the a priori groups was (values after cross-validation in parentheses): A1 – asexual L. cardui 97% (86%); A2 – asexual L. confusus 100% (93%); A3 – asexual L. fabarum 100% (89%); A4 – sexual L. fabarum 100 (62%); B – sexual L. hirticornis 100 (85%); and C – sexual L. cardui 100 (79%).
L. testaceipes (D) and L. orientalis (E) also differ clearly in wing size and wing shape (Procrustes distance is 0.051, P < 0.001). Correct classification of individuals was D – L. testaceipes 100% (98%) and E – L. orientalis 100% (100%).
The shape changes between pairs of a priori groups are illustrated in Figure 5. As previously described, the sexual and asexual lineages diverge in the shape of the distal part of the wing (described by landmarks 6, 7, 8, 9, 10 and 11), whereas different morphs often diverge in the shape of proximal area of the wing (described by landmarks 1, 2, 3, 4, 5 and 12).
By analysing the sequences of nuclear and mitochondrial genes we confirmed the existence of two previously acknowledged species groups (“testaceipes” and “fabarum”) within the genus Lysiphlebus . In terms of taxonomic characterisation, the barcoding region of the COI gene has proven to be a suitable marker for species identification within the genus Lysiphlebus [15,24,25,29], whereas the more conservative nuclear 28S D2 gene appears to be informative only at the generic or species group level. An unexpected discovery was a distinct group of COI haplotypes (clade C, Figure 2) that comprised all cardui morphs (Figure 1) with a sexual mode of reproduction. It represents a yet undescribed species that deserves further taxonomic treatment.
Marked divergence in wing shape was found between the “testaceipes” and “fabarum” species groups. However, within these groups, especially within the “fabarum” group, there is no clear correspondence between the variation in wing shape and genetic divergence. Additionally, we did not find any concordance between wing shape and the three described morphs in L. fabarum.
A significant finding in our study is that reproductive mode (namely sexual vs. asexual reproduction) is associated with wing shape. Similar shape changes between the taxa with asexual and sexual modes of reproduction were recorded in the “testaceipes” as well as in the “fabarum” group. The relationship between wing shape and reproductive mode is particularly notable in the “fabarum” species group, in which most of the morphological variation could be related to the differences between the sexual and asexual taxa and morphs.
Taxa with sexual reproduction had smaller wings than the taxa with asexual reproduction, and there were significant, size-related allometric changes in the shape of the wing. However, divergence in the wing shape was not solely the result of allometry, because the sexual and asexual lineages differed in wing shape even when the size effect was removed (allometry-free data).
Sexual reproduction (arrhenotoky) is the ancestral and dominant mode of reproduction in Hymenoptera. Asexual reproduction (thelytokous parthenogenesis) is less frequent and often induced by heritable endosymbiotic bacteria [18,19], but apparently not in Lysiphlebus. Sandrock and Vorburger  showed that asexual reproduction in L. fabarum is genetically determined and inherited as a single-locus recessive trait. Asexual reproduction has the potential to spread in parasitoid populations because, very rarely, asexual (thelytokous) females produce fertile haploid males that may carry the thelytoky-inducing allele into sexual populations, a process referred to as ‘contagious parthenogenesis’ . There are currently no data about the underlying processes of asexual reproduction in L. orientalis, but the simplest assumption would be that they are similar to those in L. fabarum, considering that thelytoky is otherwise very rare in the Aphidiinae and that rare males are also observed in natural and laboratory populations of L. orientalis . In this study, the association of mitochondrial haplotypes with the mode of reproduction in the “fabarum” group was inconsistent, and two haplotypes were shared between the asexual and sexual wasps. These results are in concordance with the statement that the “fabarum” group is an evolutionary young sexual-asexual complex with incomplete genetic isolation between the reproductive modes .
The mechanistic basis of the relationship between wing shape and reproductive mode is currently unknown. There is a lack of information on the genetic basis of wing morphogenesis in parasitic wasps, but studies on Drosophila melanogaster show that numerous genes of small effect determine wing shape in flies . Presuming that parasitoid and fly wing morphogenesis is not dramatically different, and considering that reproductive mode in L. fabarum is inherited as a single-locus trait , pleiotropy or co-inherence of the thelytoky-inducing allele and alleles affecting wing shape in the same linkage group may explain the observed link between reproductive mode and wing shape. This remains to be investigated.
Deep genetic divergence of aphid parasitoids from the genus Lysiphlebus (the “fabarum” vs. “testaceipes” group) is accompanied by changes in wing shape. At a shorter timescale within the two species groups, there is no clear correspondence between morphological evolution and genetic divergence. However, we observed a clear association between reproductive mode and wing shape; in both species groups, similar differences exist between sexual and asexual wasps, explaining much of the variation in wing shape. Further studies are necessary to resolve the mechanisms underlying the apparent relationship between reproductive mode and wing shape in Lysiphlebus wasps.
Field sampling and determination of reproductive modes
Lysiphlebus parasitoids were sampled between 2006 and 2011 in the surroundings of Belgrade, Serbia, except for L. testaceipes, which was collected along the Mediterranean coast of Montenegro. In addition to those specimens, for molecular analyses we used samples from other geographical areas (Additional file 1). Plant samples infested with live and mummified aphids were collected in the field and placed into plastic containers covered with nylon mesh . Aphid hosts were mainly from the genera Aphis and Brachycaudus, except for Metopeurum fuscoviridae Stroyan, which is parasitised by the monophagous parasitoid L. hirticornis (Additional file 1). Caged samples were held at 22.5°C until parasitoid emergence. Four collected species have known modes of reproduction: L. hirticornis – sexual; L. melandriicola – sexual; L. testaceipes – sexual; and L. orientalis – asexual. For the fabarum, cardui and confusus morphs of L. fabarum (see Figure 1), in which sexual and asexual lineages may occur, reproductive modes were mostly inferred from the field sex ratios. The complete absence of males was taken as an indicator of asexual reproduction, whereas samples containing males and females were treated as lineages with a sexual mode of reproduction. This approach is not applicable to samples with only a few individuals who are all female. To determine the reproductive modes of females from small samples (<5 individuals), we genotyped them at microsatellite locus Lysi07 , which happens to be linked to reproductive mode in parasitoids of the L. fabarum group . All females that were homozygous for allele 183 at microsatellite locus Lysi07 were treated as asexual, whereas all others were treated as sexual .
In this study we analysed the following Lysiphlebus species/morphs: L. fabarum – sexual and asexual; L. cardui – asexual; L. confusus – asexual; L. hirticornis – sexual. We also analysed sexual lineages of L. cf. cardui, which is here recorded for the first time. Sexual specimens of L. melandriicola were used only for molecular analyses due to insufficient number of individuals available for analyses of the wing shape. All of the above belong to the “fabarum” group. From the “testaceipes” group, we used L. testaceipes - sexual and L. orientalis – asexual (Additional file 1). Only female specimens were included in all molecular and morphometric analyses.
DNA extraction, PCR amplification and sequencing
To determine the genetic variation and phylogenetic relationships among the Lysiphlebus species/morphs, two molecular markers were chosen. The second expansion segment of the nuclear 28S rRNA gene (28S D2) was amplified using the primer pair 28SD2f (5’-AGAGAGAGTTCAAGAGTACGTG-3’)  and 28SD2r (5’-TTGGTCCGTGTTTCAAGACGGG-3’) . Additionally, the barcoding region of the mitochondrial cytochrome oxidase subunit I gene (COI) was amplified using the primers LCO1490 and HCO2198 . Total nucleic acids were extracted from 89 Lysiphlebus specimens (Additional file 1) and three outgroup taxa (Diaeretus essigellae Starý & Zuparko 2002, Areopraon chaitophori Tomanović & Petrović 2009 and Toxares deltiger Haliday 1883) using a nondestructive TES method . PCR reactions were performed in an Eppendorf Mastercycler® (Hamburg, Germany). Fragments of 28S D2 were amplified in a final volume of 20 μl, containing 1 μl of extracted DNA, 14.35 μl of H20, 2 μl of High Yield Reaction Buffer A, 1.5 μl of MgCl2 (2.25 mM), 0.5 μl of dNTP (0.25 mM), 1 μl of each primer (0.5 μM) and 0.15 μl of KAPATaq DNA polymerase (0.0375 U/μl) (Kapa Biosystems Inc., Boston, USA). The PCR protocol included initial denaturation at 95°C for 3 min, 30 cycles consisting of 30 s at 95°C, 30 s at 48°C, 2 min at 72°C, and the final extension at 72°C for 10 min. All COI products were amplified according to the protocol and cycling conditions described by Petrović et al. . The PCR products were checked on 1% agarose gels and purified using the QIAquick PCR Purification Kit (QIAGEN, Valencia, USA) according to the manufacturer’s instructions. DNA sequencing was performed by Macrogen Inc. (Seoul, Korea).
Sequences of both genes were manually edited in FinchTV v.1.4.0 (Geospiza, Inc., Seattle, USA; http://www.geospiza.com) and aligned using the ClustalW program implemented in the MEGA 5.1 software package . The same software was used to confirm the continuous open reading frame in the protein-coding COI gene to exclude the possibility of nuclear copies [40,41]. Kimura’s two-parameter method (K2P) of base substitution  was used to calculate an average genetic distance between the sequences.
Phylogenetic reconstruction was carried out using the methods of Maximum Likelihood (ML), Maximum Parsimony (MP) and Neighbour-Joining (NJ), all computed using the MEGA 5.1 software package. The MP tree was obtained using the Subtree-Pruning-Regrafting (SPR) algorithm . The NJ phylogenetic tree was inferred using the K2P evolutionary distances . For the reconstruction of the ML tree for 28S D2 sequences, we have identified the GTR + G model as the best-fitting model of sequence evolution based on the Bayesian Information Criterion (BIC) and the Akaike Information Criterion corrected (AICc)  as implemented in the software Modeltest 3.5 . Based on the same criteria, the HKY + G model was identified as the best-fitting model for ML reconstruction based on the COI sequences. There were a total of 618 positions in the final dataset for 28S D2 sequences and 628 positions for COI sequences in all phylogenetic reconstructions. Sequences of all haplotypes are deposited in GenBank (accession numbers are given in Additional file 1).
Geometric morphometric analyses
We applied the Generalised Procrustes Analysis (GPA) to obtain the matrix of the shape coordinates (also known as Procrustes coordinates) from which the differences due to position, scale and orientation have been removed [48,49]. The general size was computed as the centroid size (CS), which reflects the amount of dispersion around the centroid of the landmark configuration. The shape variables (Procrustes coordinates) were obtained using the MorphoJ software .
Divergence in wing size and shape
The divergence in wing size among a priori-defined groups based on phylogenetic clades, morphs or reproductive mode was analysed with Analyses of Variance (ANOVA). Post-hoc tests (Tukey’s Studentised Range) were used to test for differences between specific groups. The divergence of wing shape among groups was analysed with a Multivariate Analysis of Variance (MANOVA) using the PROC GLM procedure in SAS . To visualise the patterns of variation in wing shape, a Principle Components Analysis (PCA) on the covariance matrix of wing shape variables was carried out. We also mapped geometric morphometric data onto the molecular phylogeny and tested for a phylogenetic signal in wing shape. The generalised method of least squares [52,53] was used to find the values for the internal nodes of the phylogeny from the shape averages of the terminal taxa [53-55]. The phylogenetic signal in the shape data was tested by a permutation approach using the MorphoJ software , which simulated the null hypothesis of a complete absence of phylogenetic structure by randomly reassigning the phenotypic data to the terminal nodes .
To quantify the shape differences, Procrustes distances were calculated between each pair of analysed groups. The statistical significance of differences in the mean shape between groups was estimated using a permu tation test based on 10,000 iterations. The statistically significant differences after Bonferroni correction are presented.
Differences among the a priori-defined groups (clades, morphs, and reproductive modes) were explored with a Discriminant Function (DF) analysis. We report both original and cross validation percentages to better estimate the uncertainty in assigning individuals to groups based on wing shape . To further explore variation among the sexual and asexual lineages, we performed separate PCA analysis of the covariance matrices of wing shape data for the species groups “fabarum” and “testaceipes”. To determine the degree to which wing shape variation was associated with size variation (allometry), we performed a multivariate regression of wing shape variables on log CS. The multivariate regression was performed for the “fabarum” and “testaceipes” species groups separately. To explore allometry-free shape data, we performed PCA of covariance matrices using residuals obtained from multivariate regression as size-free shape variables.
Availability of supporting data
The data sets supporting the results of this article are available in the Dryad.org repository, doi:10.5061/dryad.9t0g2, http://datadryad.org/review?doi = doi:10.5061/dryad.9t0g2.
Nucleotide sequences are submitted to GenBank database under the accession numbers [GenBank: KP663427- KP663464].
We wish to thank Dr. Olivera Petrović-Obradović (Faculty of Agriculture, University of Belgrade) for the identification of aphids. This study was supported by the SCOPES program of the Swiss National Science Foundation (grant IZ73Z0_128174) and partly supported by the grant III43001 (The Ministry of Education, Science and Technological Development of the Republic of Serbia). The participation by P. Starý was partially supported from the Entomology Institute project RVO: 60077344.
- Lande R, Arnold SJ. The measurement of selection on correlated characters. Evolution. 1983;37:1210–26.View ArticleGoogle Scholar
- Richardson MK, Chipman AD. Developmental constraints in a comparative framework: a test case using variations in phalanx number during amniote evolution. J Exp Zool B Mol Dev Evol. 2003;296B:8–22.View ArticleGoogle Scholar
- Galis F, Metz JAJ. Evolutionary novelties: the making and breaking of pleiotropic constraints. Integr Comp Biol. 2007;47:409–19.View ArticlePubMedGoogle Scholar
- Blomberg SP, Garland T, Ives AR. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution. 2003;57:717–45.View ArticlePubMedGoogle Scholar
- Losos JB. Seeing the forest for the trees: the limitations of phylogenies in comparative biology. Am Nat. 2011;177:709–27.View ArticlePubMedGoogle Scholar
- Losos JB. Ecomorphology, performance capability, and scaling of West Indian Anolis lizards: an evolutionary analysis. Ecol Monograph. 1990;60:369–88.View ArticleGoogle Scholar
- Gidaszewski NA, Baylac M, Klingenberg CP. Evolution of sexual dimorphism of wing shape in the Drosophila melanogaster subgroup. BMC Evol Biol. 2009;9:110.View ArticlePubMed CentralPubMedGoogle Scholar
- Drake AG, Klingenberg CP. Large-scale diversification of skull shape in domestic dogs: disparity and modularity. Am Nat. 2010;175:289–301.View ArticlePubMedGoogle Scholar
- Sherratt E, Gower DJ, Klingenberg CP, Wilkinson M. Evolution of cranial shape in caecilians (Amphibia: Gymnophiona). Evol Biol. 2014;41:528–45.View ArticleGoogle Scholar
- Price PW. Evolutionary biology of parasites. Princeton: Princeton University Press; 1980.Google Scholar
- Strong DR, Lawton JH, Southwood R. Insects on plants. Community patterns and mechanisms. Oxford: Blackwell Scientific; 1984.Google Scholar
- Kavallieratos NG, Tomanović Ž, Starý P, Athanassiou CG, Sarlis GP, Petrović O, et al. A survey of aphid parasitoids (Hymenoptera: Braconidae: Aphidiinae) of southeastern Europe and their aphid-plant associations. App Entomol Zool. 2004;39:527–63.View ArticleGoogle Scholar
- Starý P. Aphid parasitoids of the Czech Republic. Prague: Academia; 2006.Google Scholar
- Starý P, Rakhshani E, Tomanović Ž, Hoelmer K, Kavallieratos NG, Yu J, et al. A new species of Lysiphlebus Förster 1862 (Hymenoptera: Braconidae: Aphidiinae) attacking soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae) from China. J Hymenoptera Res. 2010;19:179–86.Google Scholar
- Petrović A, Mitrović M, Starý P, Petrović-Obradović O, Tomanović Ž, Žikić V, et al. Lysiphlebus orientalis, a new invasive parasitoid in Europe – evidence from molecular markers. Bull Entomol Res. 2013;103:451–7.View ArticlePubMedGoogle Scholar
- Ameri M, Rasekh A, Mohammadi Z. A comparison of life history traits of sexual and asexual strains of the parasitoid wasp, Lysiphlebus fabarum (Braconidae: Aphidiinae). Ecol Entomol. 2014. doi:10.1111/een.12155.Google Scholar
- Žikić V, Petrović A, Ivanović A. Allometric shape changes indicate significant divergence in the wing shape between asexual and sexual lineages of Lysiphlebus fabarum (Marshall) (Hymenoptera: Braconidae). Acta Entomol Serbica. 2014;19:53–62.Google Scholar
- Cook JM. Sex determination in the Hymenoptera-a review of models and evidence. Heredity. 1993;71:421–35.View ArticleGoogle Scholar
- Heimpel GE, de Boer JG. Sex determination in the Hymenoptera. Annu Rev Entomol. 2008;53:209–30.View ArticlePubMedGoogle Scholar
- Starý P, Tomanović Ž, Petrović O. A new parasitoid of root-feeding aphids from the Balkan mountains (Hymenoptera, Braconidae, Aphidiinae). Deut Entomol Z. 1998;45:175–9.View ArticleGoogle Scholar
- Kavallieratos NG, Tomanović Ž. Some rare and endemic aphid parasitoid species (Hymenoptera: Braconidae: Aphidiinae) from the Balkan Peninsula. Acta Entomol Serbica. 2001;6:121–9.Google Scholar
- Belshaw R, Quicke DLJ, Völkl W, Godfray HCJ. Molecular markers indicate rare sex in a predominantly asexual parasitoid wasp. Evolution. 1999;53:1189–99.View ArticleGoogle Scholar
- Starý P. Biology and distribution of microbe-associated thelytokous populations of aphid parasitoids (Hymenoptera, Braconidae, Aphidiinae). J Appl Entomol. 1999;123:231–5.View ArticleGoogle Scholar
- Sandrock C, Schirrmeister B, Vorburger C. Evolution of reproductive mode variation and host associations in a sexual-asexual complex of aphid parasitoids. BMC Evol Biol. 2011;11:348. doi:10.1186/1471-2148-11-348.View ArticlePubMed CentralPubMedGoogle Scholar
- Derocles SAP, Le Ralec A, Plantegenest M, Chaubet B, Cruaud C, Cruaud A, et al. Identification of molecular markers for DNA barcoding in the Aphidiinae (Hym. Braconidae). Mol Ecol Resour. 2012;12:197–208.View ArticlePubMedGoogle Scholar
- Gilchrist AS, Azevedo RB, Partridge L, O’Higgins P. Adaptation and constraint in the evolution of Drosophila melanogaster wing shape. Evol Dev. 2000;2:114–24.View ArticlePubMedGoogle Scholar
- Klingenberg CP, Zaklan SD. Morphological integration between developmental compartments in the Drosophila wing. Evolution. 2000;54:1273–85.View ArticlePubMedGoogle Scholar
- Klingenberg CP, Badyaev AV, Sowry SM, Beckwith NJ. Inferring developmental modularity from morphological integration: analysis of individual variation and asymmetry in bumblebee wings. Amer Nat. 2001;157:11–23.View ArticleGoogle Scholar
- Mitrović M, Petrović A, Kavallieratos NG, Starý P, Petrović-Obradović O, Tomanović Ž, et al. Geographic structure with no evidence for host-associated lineages in European populations of Lysiphlebus testaceipes, an introduced biological control agent. Biol Control. 2013;66:150–8.View ArticleGoogle Scholar
- Sandrock C, Vorburger C. Single-locus recessive inheritance of asexual reproduction in a parasitoid wasp. Curr Biol. 2011;21:433–7.View ArticlePubMedGoogle Scholar
- Engelstädter J, Sandrock C, Vorburger C. Contagious parthenogenesis, automixis, and a sex determination meltdown. Evolution. 2011;65:501–11.View ArticlePubMedGoogle Scholar
- Carreira VP, Soto IM, Mensch J, Fanara JJ. Genetic basis of wing morphogenesis in Drosophila: sexual dimorphism and non-allometric effects of shape variation. BMC Dev Biol. 2011;11:32.View ArticlePubMed CentralPubMedGoogle Scholar
- Kavallieratos NG, Lykouressis DP, Sarlis GP, Stathas GJ, Sanchis-Segovia A, Athanassiou CG. The Aphidiinae (Hymenoptera: Ichneumonoidea: Braconidae) of Greece. Phytoparasitica. 2001;29:306–40.View ArticleGoogle Scholar
- Sandrock C, Frauenfelder N, von Burg S, Vorburger C. Microsatellite DNA markers for the aphid parasitoid Lysiphlebus fabarum and their applicability to related species. Mol Ecol Notes. 2007;7:1080–3.View ArticleGoogle Scholar
- Belshaw R, Quicke DLJ. A molecular phylogeny of the Aphidiinae (Hymenoptera: Braconidae). Mol Phylogenet Evol. 1997;7:281–93.View ArticlePubMedGoogle Scholar
- Campbell BC, Steffen–Campbell JD, Werren JH. Phylogeny of the Nasonia species complex (Hymenoptera: Pteromalidae) inferred froman internal transcribed spacer (ITS2) and 28S rDNA sequences. Insect Mol Biol. 1993;2:225–37.View ArticlePubMedGoogle Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3:294–9.PubMedGoogle Scholar
- Mahuku GS. A simple extraction method suitable for PCR based analysis of plant, fungal, and bacterial DNA. Plant Mol Biol Rep. 2004;22:71–81.View ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang DX, Hewitt GM. Nuclear integrations-challenges for mitochondrial- DNA markers. Trends Ecol Evol. 1996;11:247–51.View ArticlePubMedGoogle Scholar
- Song H, Buhay JE, Whiting MF, Crandall K. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proc Natl Acad Sci U S A. 2008;105:13486–91.View ArticlePubMed CentralPubMedGoogle Scholar
- Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.View ArticlePubMedGoogle Scholar
- Nei M, Kumar S. Molecular evolution and phylogenetics. New York: Oxford University Press; 2000.Google Scholar
- Posada D, Crandall KA. MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998;14:817–8.View ArticlePubMedGoogle Scholar
- Tomanović Ž, Kos K, Petrović A, Starý P, Kavallieratos NG, Žikić V, et al. The relationship between molecular variation and variation in the wing shape of three aphid parasitoid species: Aphidius uzbekistanicus Luzhetzki, Aphidius rhopalosiphi De Stefani Perez and Aphidius avenaphis (Fitch) (Hymenoptera: Braconidae: Aphidiinae). Zool Anz. 2013;252:41–7.View ArticleGoogle Scholar
- Rohlf FJ. TpsDig program, version 2.04, ecology and evolution, SUNY at stony brook. New York: State University of New York; 2005. http://life.bio.sunysb.edu/morph/.Google Scholar
- Sharkey MJ, Wharton RA. Morphology and terminology. In: Wharton RA, Marsh PM, Sharkey MJ, editors. Manual of the new world genera of the family braconidae (hymenoptera). special publication 1. Washington, DC: International Society of Hymenopterists; 1997. p. 19–37.Google Scholar
- Rohlf FJ, Slice E. Extensions of the Procrustes method for the optimal superimposition of landmarks. Syst Zool. 1990;39:40–59.View ArticleGoogle Scholar
- Dryden IL, Mardia KV. Statistical shape analysis. New York: Wiley; 1998.Google Scholar
- Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Resour. 2011;11:353–7.View ArticlePubMedGoogle Scholar
- SAS Institute Inc: Base SAS® 9.1.3. Base SAS® 9.1.3. Procedures Guide. Cary, NC: SAS Institute In; 2009.Google Scholar
- Martins EP, Hansen TF. Phylogenies and the comparative method: a general approach to incorporating phylogenetic information into the analysis of interspecific data. Am Nat. 1997;149:646–67.View ArticleGoogle Scholar
- Rohlf FJ. Comparative methods for the analysis of continuous variables: geometric interpretations. Evolution. 2001;55:2143–60.View ArticlePubMedGoogle Scholar
- Maddison W. Squared-change parsimony reconstructions of ancestral states for continuous-valued characters on a phylogenetic tree. Syst Zool. 1991;40:304–14.View ArticleGoogle Scholar
- McArdle B, Rodrigo AG. Estimating the ancestral states of a continuous-valued character using squared-change parsimony: An analytical solution. Syst Biol. 1994;43:573–8.View ArticleGoogle Scholar
- Klingenberg CP, Gidaszewski NA. Testing and quantifying phylogenetic signals and homoplasy in morphometric data. Syst Biol. 2010;59:245–61.View ArticlePubMedGoogle Scholar
- Manly BFJ. Randomization, bootstrap and Monte Carlo methods in biology. London: Chapman & Hall; 1997.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.