Translator Disclaimer
1 March 2014 Development of Microsatellite Markers and Detection of Genetic Variation between Goniozus Wasp Populations
Author Affiliations +

Molecular genetic markers reveal differences between genotypes according to the presence of alleles (the same or different) at target loci. Microsatellite markers are especially useful codominant markers that have been used in a wide range of studies to elucidate the population structure and dynamics of a range of organisms, including agriculturally beneficial insects such as parasitic wasps (parasitoids). In the present study, twelve primer pairs were designed for the south Asian , Goniozus nephantidis (Muesebeck) (Hymenoptera: Bethylidae), and 24 for its New World congener, Goniozus legneri Gordh, parasitoids of the larvae of the lepidopteran coconut pest Opisina arenosella Walker (Lepidoptera: Crytophasidae) and other lepidopteran pests, respectively, in order to investigate polymorphism within and between populations. The wasps fingerprinted were a total of 85 G. nephantidis and G. legneri, including individuals belonging to three putatively different strains of G. legneri. Annealing gradient tests (50–65°C) were conducted to study the quality of the PCR amplification across an annealing temperature gradient using a mixed genotype DNA template from each species separately. Seven primer pairs, which amplified clear products of approximately the expected size of G. nephantidis and 18 of G. legneri, were then selected for capillary analysis for fragment size determination on a Beckmann CEQ 8000. Neither G. nephantidis nor G. legneri were polymorphic within populations. However, there were six primer pairs that did show polymorphism between G. legneri populations that originated from different geographical areas within South America (Uruguay and Chile). Furthermore, one primer pair revealed diversity between the two strains collected within Chile. One of the markers was subsequently used to provide unbiased assessment of primary sex ratio in G. legneri.


Among the natural enemies of agricultural insect pests, hymenopteran parasitoids are one of the most important classes of biological control agents and are also widely used in studies of evolutionary ecology and basic population biology (Godfray 1994; Jervis 2005; Hochberg and Ives 2000; Wajnberg et al. 2008; Hardy et al. 2013). The distributions and population structures of parasitoids are influenced by a wide range of factors, such as geological and geographical components, ecological processes, and evolutionary and genetic aspects (Zink 2002; Bond and Stockman 2008). Successful population genetic, ecological, and evolutionary studies can be achieved through the availability of suitable molecular markers, which are important indicators of relationships between both individuals and populations (Carvalho 1998). Such markers can reveal differences between genotypes through the application of a range of random markers not linked a priori to traits.

Among the classes of widely used genetic markers are ‘microsatellite’ markers, also known as simple sequence repeats (SSRs). These consist of short, repeated units of around two to six base pairs in length with an array that can be up to 200 bp long and are found in both coding and non-coding regions in all prokaryotic and eukaryotic genomes (Tautz 1989; Arcot et al. 1995; Beukeboom and Zwaan 2005). To date, they have been developed in a number of parasitoids, especially braconids (Baker et al. 2003; Anton et al. 2006; Lozier et al. 2006).

Microsatellite markers, which are codominant (Loxdale and Lushai 1998), have many advantages over other marker types because not only do they generally have a high number of alleles per locus, which can identify polymorphism, but they can detect high levels of heterozygosity and have high mutation rates (Hancock 1999). As such, microsatellite markers have been used to determine the genetic diversity and differentiation between populations through measuring the degree of heterozygosity in parasitoid species (e.g. 0.04–0.44 in Cotesia melitaearum, 0.171–0.629 in Neotypus melanocephalus, and 0.170–0.367 in Lysiphlebus hirticornis; Kankare et al. 2005; Anton et al. 2007; Nyabuga et al. 2010) as well as the degree of gene flow and dispersal between populations (Avise 1994; McCoy et al. 2001; Molbo et al. 2003; Kankare et al. 2005; Zavodna et al. 2005; Drescher et al. 2010; Nyabuga et al. 2010).

In the present study, we designed a suite of microsatellites for screening two species of bethylid wasps for genetic polymorphisms within and between populations. In principle such markers could also prove useful for pest control applications (Aebi et al. 2008; Ugelvig et al. 2008; Lozier et al. 2009; Zygouridis et al. 2009; Nicholls et al. 2010; Lavandero et al. 2011), evaluating the effect of kinship on social behaviors (Lizé et al. 2012), and for measuring population parameters, such as levels of inbreeding, which have not been directly evaluated but are important in the understanding of reproductive decisions (Hardy and Cook 1995; Hardy et al. 1998, 1999, 2000). The first direct application of these markers has been to provide assessment, as described elsewhere (Khidr et al. 2013), of the sex of individual parasitoid eggs in order to evaluate maternal sex allocation without the biasing influence of developmental mortality.

The Bethylidae is a family of parasitoid hymenopteran wasps that has been thought to comprise four extant subfamilies, Bethylinae, Epyrinae, Pristocerinae, and Mestitiinae (Evans 1964), with over 2000 described species (Gordh and Móczár 1990). Recently, the higher level phylogeny of bethylids has been estimated using molecular data from 33 species, resulting in a split of the sub-family Mestitiinae into two separate sub-families, the Mestitiinae and the Cephalonomiini (Carr et al. 2010). Bethylid wasps attack almost exclusively the immature stages of coleopterans and lepidopterans, many of which are pests of important agricultural commodities such as coffee, coconut, sugarcane, apple, walnut, and almonds (Gordh 1982; Batchelor at al. 2005; Venkatesan et al. 2007; Zaviezo et al. 2007). In this study, we have focused on Goniozus nephantidis (Muesebeck) (Hymenoptera: Bethylidae), a parasitoid of the lepidopteran larvae of the coconut pest Opisina arenosella Walker (Lepidoptera: Crytophasidae) in the Indian sub-continent, and G. legneri Gordh, a parasitoid of several New World lepidopteran pests of walnuts, pistachio nuts, almonds, and apples (Steffan et al. 2001; Garrido et al. 2005; Zaviezo et al. 2007).

Both G. nephantidis and G. legneri have each been used in biocontrol programs (Dharmaraju 1963; Legner and Silveira-Guido 1983; Gothilf and Mazor 1987; Lyla et al. 2006) and in a range of behavioral ecological studies (e.g., Hardy and Cook 1995; Goubault et al. 2006, 2007; Humphries et al. 2006; Bentley et al. 2009; Lizé et al. 2012). Their basic life histories are similar; both are gregarious idiobiont ectoparasitoids exhibiting sub-social behavior, such as maternal care and defense of the developing brood (Hardy and Blackburn 1991; Bentley et al. 2009), and appear to conform closely, but probably not exactly, to single foundress local mate competition (Hamilton 1967; Hardy and Cook 1995; Hardy et al. 1998, 1999, 2000). Males usually emerge before females and have sufficient capacity to inseminate their sisters (Hardy et al. 1999, 2000), and, as with many other bethylids, the sex ratios of these species are generally female biased with a low degree of sex ratio variation between broods (subbinomial variance: Green et al. 1982; Hardy and Mayhew 1998; Hardy et al. 1998; Khidr et al. 2013).

Materials and methods

Parasitoid origins and cultures

The culture of G. nephantidis used in this study had been maintained in the laboratory for more than 20 years on the facultative host Corcyra cephalonica Stainton (Lepidoptera: Pyralidae), as described in Lizé et al. (2012). Corcyra cephalonica was also used as the facultative host for G. legneri. Three strains of G. legneri were used. One, termed strain ‘U’, was obtained from a commercial insectary in the USA and kept in our laboratory for more than eight years. The original material is believed to have been collected from southern Uruguay in 1978 (Gordh 1982; Gordh et al. 1983; Legner and Silveira-Guido 1983). Two further strains of G. legneri were brought to our laboratory in May 2009 from Santiago, Chile. One strain was collected directly from walnut trees and was termed ‘C-field’, and the other strain was termed ‘C-lab’, as it had been maintained in a Chilean insectary for several years following collection from a field site near Santiago (Zaviezo et al. 2007). All cultures were maintained in a constant environment room at 25–27°C, 12:12 L:D, and with high relative humidity maintained by an ambient temperature water bath.

Design and preparation of the primers

Microsatellite-enriched genomic libraries were created essentially according to Kloda et al. (2004), with the final sequencing step performed using barcoded adaptors and a 1/16th run of non-titanium reagents Roche 454 pyrosequencing (as part of a mixture of 9 different libraries) ( The generated Fasta files were separated in silico to identify the individual libraries, and those for G. legneri and G. nephantidis were searched for microsatellite motifs using the script ( Primer pairs flanking the simple sequence repeats were designed either by Primer 3 (Rozen and Skaletsky 2000) and/or WebSat (Martins et al. 2009) for G. legneri (Table 1) and for G. nephantidis (Table 2). Primers were synthesized by Eurofins MWG Operon ( with a forward primer 5′ extension consisting of the M13 sequence to allow fluorescent labeling of the final product through a three primer reaction (Schuelke 2000), and prepared to 1000× concentration using Sigma-Aldrich ( www.sigmaaldrich) molecular biology grade water (to create primer stocks of 200 pmol/µL). After vortexing and spinning, tubes were placed on ice for 30 min. Primers were kept in a freezer at -20°C, and to produce a 10× primer stock, 5µL of the 1000× stock was mixed with 495 µL sterile distilled water for both forward and reverse primers in separate tubes on ice. The third primer (M13;TGTAAAACGACGGCCAGT-3′) was ordered from Sigma-Aldrich and labeled with dye D4 (blue; WellRed dyes).

DNA extraction

A sample size of 85 adult individuals was used for capillary testing. Seventeen individuals of G. nephantidis were examined, plus a total of 68 individuals for the different strains of G. legneri (25 of ‘U’ strain, 22 of ‘C-lab’ strain, and 21 individuals of the ‘C-field’ strain). In addition, five pooled samples of 20 individuals were used for the annealing gradient test. Individual adult females or pooled adult samples were placed in 1.5 mL Eppendorf tubes (Sarstedt, and then immersed into liquid nitrogen and crushed using a mini pestle to start the extraction. Genomic DNA was then extracted either by using a GenElute plant Genomic kit (Sigma-Aldrich) or by following, with some modifications, methods given in Sambrook et al. (1989) and Vogler and Desalle (1993) before elution/re-suspension into 50–60 µL of sterile distilled water and storage at -20°C.

Polymerase chain reaction (PCR)

PCR reactions were performed in either a Thermo Hybaid Express PCR machine ( or in an ABI PCR 9700 Thermocycler machine (Applied Biosystems, Life Technologies, The Thermo Hybaid Express PCR was used for annealing gradient tests and run with a 15°C gradient by using a total volume of 20 µL for one reaction through mixing different components consisting of 2 µL of 10× forward primer and 2 µL of 10× reverse primer (2 pmol/µL final); 2 µL of 10× PCR buffer; 0.16 µL of dNTP's (mixed dNTP 25 mM final concentration per nucleotide); 2 µL of DNA template (mixture of many individuals); 0.10 µL of Taq DNA polymerase (5 units/µL), and 11.7 µL of sterile distilled water. Master mixes were prepared, where possible, to decrease the effects of pipetting error.

Thus, the optimal annealing temperatures for each of the primer pairs used were determined according to following program: an initial 3 minutes denaturation at 94°C, followed by 35 cycles of 1 minute at 94°C (denaturation), 1 minute at 50–65°C (annealing), and 72°C for 2 minutes (extension), with a final extension of 72°C for 10 minutes at the end of the program.

For genotyping of individual samples amplified in the ABI PCR 9700 Thermocycler, the aforementioned program was used, but the determined optimum temperature for the annealing step was 60°C for the majority of the primers unless otherwise stated. PCR reactions consisted of 0.2 µL of 10× forward primer and 2 µL of 10× reverse primer (2 pmol/µL final), 2 µL (10×) PCR buffer, 0.16 µL dNTP's (each in a 25 mM final concentration), ∼ 0.05 µL M13 Blue Taq of 1000× (53.8 nM concentration), 2 µL of individual genomic DNA (∼ 5 ng/µL), 0.10 Taq DNA polymerase (5 units/µL), and 13.5 µL sterile distilled water. Thus, a fluorescently-labeled M13 tail sequence was added to the 5′-end of the forward primer (Schuelke 2000) to be used for capillary sequencing.

Agarose gel electrophoresis

Samples to be loaded onto the gel were mixed with 6× gel loading blue buffer (Promega, in the ratio of one part sample to one part loading buffer. Loading buffer was added to each well of the plates from the PCR machine and was spun briefly. Thereafter 10 µL from each well was loaded onto a submerged gel that consisted of 2% agarose (molecular grade, Bioline, prepared in 0.5× TBE (Tris-Borate-EDTA, pH 8.0) buffer, followed by addition of 2 µL of ethidium bromide stock before pouring (10 mg/mL stock; Promega). Each primer pair reaction was loaded onto one row of the gel (each primer pair having 12 reactions across a 15°C annealing gradient). Alongside, an appropriate size marker (5 µL of 2-log DNA ladder; New England Biolabs, was loaded in the first lane of each primer pair, and the gel was run at 90 V for approximately 1 hour. Following electrophoresis, bands were visualized and photographed under UV-light in a Bio-Rad Gel Doc 2000 gel box ( ).

Quantification test of DNA templates and Goniozus individuals’ DNA extractions were carried out by comparison with known uncut lambda DNA (BioLabs) (50 ng/µL) loaded in the amounts of 10 µL, 5 µL, 2.5 µL, and 1.5 µL to provide a fluorescence comparison with the unknown samples. The 1% agarose gel was run at 90 V for 75 min, then different individuals of both species were quantified and tested for DNA integrity by ensuring that genomic samples largely ran at limiting mobility. Good quality samples were used as DNA templates.

Capillary sequencing: Preparing fragment samples for analysis

The CEQ 8000 Fragment Analysis Software Version 8 (Beckman Coulter, was used to measure and analyze the PCR product fragment sizes. For the preparation of the sample in half-reactions, for each row of eight samples, 215 µL of SLS (sample loading solution) was added to 2 µL of SS 400 (standard size) mixed by vortexing and spun briefly. 27 µL of this mixture was then added to each well in the row, and 2 µL of multiplexed PCR product was later added. The mixture in each well was overlaid immediately with a drop of mineral oil and placed in the CEQ machine. Later, cluster analysis between different populations of G. legneri was performed using the Multi-Variate Statistical Package version 3.2 (MVSP; Kovach Computing Services,


Primer design in the Microsatellite library

A genomic library consisting of 273 sequences containing microsatellite motifs was screened to design primers for G. legneri. Twenty-four of these sequences were consid ered to be clearly unique sequences with adequate flanking sequence length to design primer pairs flanking the repeat unit. The dinucleotide (GA)n repeat was the predominant marker, followed by (AG)n and (TC)n, while the tri- and tetra-nucleotide microsatellites were frequently present as compound microsatellites. In G. nephantidis, there were 3356 microsatellites, of which 12 were chosen to design primers for the investigation of polymorphism within the population.

Figure 1.

Annealing gradient for primers 7 and 8 of Goniozus legneri (primer labels correspond to those in Table 1). High quality figures are available online.


Figure 2.

Gel plage of PCR products for primers 7 and 8 of Goniozus legneri. High quality figures are available online.


Annealing gradient tests

PCR analysis was performed to optimize annealing gradients for the 12 new G. nephantidis primer pairs and 24 primer pairs for G. legneri strains. Occasionally a number of primer pairs were designed to the same microsatellite repeat sequence to increase the probability of successful amplification. The process was repeated several times to test the reliability of the new primers. Representative results of the annealing electrophoresis gels for G. legneri primers are shown in Figure 1, with the left hand column showing the 2-log DNA ladder (New England Biolabs).

Primers displaying clear bands in annealing tests were selected for PCR amplification and polymorphism testing. In addition, the best annealing temperature for each primer was recorded in order for it to be used for the PCR. According to the results of the annealing tests, the best temperature for all G. nephantidis primers was 60°C, except for primer GnSSR11 at 56°C. No amplification was observed for primers GnSSR1, 2, 3, 4, and 10 in the test. In G. legneri, the optimum temperatures were 54°C, 56°C, and 58°C for primers GlSSR 14, 5, and 22, respectively. The remainder of the primer pairs had optimal annealing temperatures close to 60°C, allowing simultaneous amplification in the same thermoblock. The following primers were excluded from further work: GlSSR1, 4, 6, 10, 19*b, 20*, 21*, 22*, and 23.

DNA quality test

DNA quality tests were conducted for Goniozus that were to be used as a template, and samples were diluted for annealing tests. However, DNA preparations from individual wasps for use in PCR did not usually need dilution because concentrations were generally at or below 10 ng/µL.

Second round of PCR runs

Primers chosen in the annealing test were amplified on the ABI PCR machine at different temperatures according to their annealing test optima. The gel electrophoresis results were visualized using UV light. Some primers were rejected before capillary testing due to lack of amplification in an annealing test or, if they amplified, not displaying clear/discrete single bands on the gel. PCR products of primers 7 and 8 amplified from G. legneri are shown in Figure 2.

Capillary sequencing

The results of the capillary fragment analysis were processed using the CEQ 8000 DNA sequencer to determine amplified product size. Neither G. nephantidis nor G. legneri were polymorphic within strains. Nonetheless, there were six primers that showed clear inter-strain polymorphism in G. legneri (Table 3). For instance, primer GlSSR7 showed a large size difference between strains, while U-strain was 137 bp and both Chilean populations were 153 bp. Furthermore, primer GlSSR5 showed a different size allele at 228 bp for both U and C-lab strains, with 224 bp for the C-field strain (Figure 3). Dendrogram analysis based on the overall microsatellite alleles patterns showed clear differences between populations from the two different geographical locations, Uruguay and Chile, as well as differentiation within the two strains collected in Chile (Clab and C-field). However, both strains were more closely related to each other than either was to the Uruguay strain (Figure 4).

Figure 3.

Examples of capillary tube analysis of primer 7 (A) and primer 8 (B) for Goniozus legneri strains. High quality figures are available online.


Figure 4.

Dendogram showing relationships between and within strains of Goniozus legneri. High quality figures are available online.



The main aim of this study was to develop molecular markers for bethylid wasps so that these markers could subsequently be utilized in evolutionary ecology as well as in applied agricultural research. In general, the percentage of amplified loci decreases with increasing genetic distance between the species tested, making such markers most suitable for closely-related species, such as congeners (Hancock and Simon 2005; Barbara et al. 2007). Since there are well over 100 described species of Goniozus (Gordh and Móczár 1990), such markers have considerable potential in many diverse studies, both pure and applied, on this genus of parasitoids.

The molecular data from this study reveal a lack of variation within strains of both G. nephantidis and G. legneri. Potential reasons for this include a loss of genetic diversity during laboratory maintenance (e.g., due to relatively small populations in culture, occasional population crashes, and the exposure of deleterious recessive genes in haploid males) (Unruh et al. 1984; Shields 1993; Cook 1993; Henter 2003). However, this explanation would most likely apply to G. nephantidis and the U-strain of G. legneri, as these have been maintained in culture for many more years than the strains of G. legneri from Chile (Clab and C-field). The fact that polymorphisms were not found within the much more recently collected Chilean strains suggests that lack of variation is not a result of genetic drift due to small population and time-in-culture effects. A more probable explanation is that there is limited genetic variation within each of the populations from which field collections were made. Genetic homozygosity can result from inbreeding because close relatives mate more frequently than expected by chance given the overall size of the population (Henter 2003; Elias et al. 2010; Mazzi et al. 2011), and both species of Goniozus are known to exhibit high levels of pre-dispersal sibling mating in the laboratory (Hardy et al. 1999, 2000) and have sex ratios that largely conform to theoretical expectations under such ‘local mate competition’ (Gordh et al. 1983; Hardy and Cook 1995; Hardy et al. 1998; Khidr et al. 2013). It is further known that G. nephantidis does not exhibit inbreeding depression in terms of effects on developmental mortality or sex ratio control (Cook 1993). Nonetheless, to date the post-dispersal mating behavior of these two species has not been directly evaluated, although the natural mating systems of these wasps are likely to have a large effect on the evolution of their sex ratios (Hardy 1994; Hardy and Cook 1995; Hardy and Mayhew 1998; Hardy et al. 1998). The lack of withinstrain genetic polymorphism that was observed provides a degree of evidence that sibling-mating is the predominant feature of the natural mating system of both Goniozus species.

For G. legneri strains collected within the same geographical region of Chile, C-lab and C-field, an allelic difference was observed with just one primer of the six amplified. This could relate to the fact that the field strain was collected from carob moth larvae (Ectomyelois ceratoniae Zeller (Lepidoptera: Pyralidae)) feeding on walnuts, while the laboratory strain was derived from a mixture of individuals collected from both carob moth on walnuts and codling moth larvae (Cydia pomonella L. (Lepidoptera: Tortricidae)) feeding on apples (Zaviezo et al. 2007; T. Zaviezo personal observation, I. Hardy personal observation). Genetic diversity in several other parasitoid species has been found to be associated with host and host plant species (e.g., Kavallieratos et al. 2004; Stireman et al. 2006).

Six of the primer pairs tested showed clear microsatellite polymorphisms between G. legneri strains (U and C-field). Genetic differences between these strains could arise due to differences in the species of insect host or host plant they were collected from (see above) and also to geographical differences (Menken 1981; Ruiz-Montoya et al. 2003; Stireman et al. 2006; Pannebakker et al. 2008; Phillips et al. 2008; Lozier et al. 2009; Lavandero et al. 2011). Various biological traits and genetic diversity might associate with populations from different geographic localities due to geographic isolation, influenced by different climatic effects and thus different selection pressures (Diehl and Bush 1984; Hopper et al. 1993; Thompson 1994; Goodisman et al. 2001; Hufbauer et al. 2004).

Establishing that there is genetic polymorphism within parasitoid species opens up a number of possibilities for the use of these markers. For instance, molecular markers have been used to show that the host searching behavior of different Agathis sp. n. populations was not affected by geographical structure and they have the ability to disperse for long distances (Althoff and Thompson 2001), and reciprocal crossing of two geographically and host-species distinct strains of Aphelinus albipodus showed reproductive compatibility and no reduction in fecundity (Wu et al. 2004). Furthermore, genetic relatedness is a crucial factor in the evolution of social behaviors between individuals (Hamilton 1964a, b; Mateo 2004; Lizé et al. 2006; Gardner and West 2007), and relatedness between insects can be usefully assessed using microsatellite markers (Buczkowski et al. 2004; Trindl et al. 2004; Jaquiery et al. 2005; Drescher et al. 2010), leading to key insights into behaviors such as kin-based altruism and aggression assays (Giraud et al. 2002; Tsutsui et al. 2003; Drescher et al. 2010; El-Showk et al. 2010). In our own study system, the development of microsatellite markers provides useful support for empirical work on kin recognition mechanisms (Lizé et al. 2012), as they confirm the assumption that females from different strains of G. legneri derive from populations with different genetic backgrounds, and thus are not as closely related as are females from within the same strain. In general, genetic recognition cues will usually associate with the level of polymorphism (Ratnieks 1991; Buczkowski et al. 2004), and the degree of aggressive behavior between encountered individuals is attuned to the level of genetic diversity recognition loci (Giraud et al. 2002; Drescher et al. 2010).

Microsatellite markers have also been used in studies of hymenopteran sex ratios to identify the sex of eggs (Ratnieks and Keller 1998; Abe et al. 2009). Assessment of the primary sex ratios of the parasitoid Melittobia australica showed that sex allocation is under precise control with the sexes produced in a regular sequence throughout the period of oviposition (Abe et al. 2009). The microsatellites developed in the present study have also been directly applied to the molecular-genetic detection of haploid (male) and diploid (female) eggs in G. legneri (Khidr et al. 2013). This provides an evaluation of primary sex ratios that is unbiased by developmental mortality (a longstanding obstacle in sex allocation research on many species, e.g., Fiala 1980; Hardy and Cook 1995; Hardy et al. 1998; Krackow and Neuhäuser 2008; Abe et al. 2009). The consistent between-strain polymorphisms (U and C-field) and cross-mated mothers were utilized, such that haploid and diploid eggs had different marker compositions. This work showed, for instance, that relationships between sex ratio and group size can be obscured by developmental mortality when the sex of eggs is not assessed directly, and also that male and female eggs may tend to be laid in spatial separation (Khidr et al. 2013).

In conclusion, the G. nephantidis laboratory culture evaluated was found not to be polymorphic in terms of the 12 microsatellite markers presently developed. This likely reflects the limited genetic variability within this population but may be due to a prolonged period in laboratory culture. For G. legneri, no polymorphisms were found within strains using the 23 designed markers. As some strains were recently collected from the field, this finding suggests natural genetic variation is locally limited. However, there were six primers that showed clear between-strain marker polymorphism in G. legneri. Six markers differed between strains collected recently in Chile and strains believed to originate from Uruguay several decades ago, while the two Chilean strains differed in only one microsatellite marker.

These markers have already proved useful for experimental work on kin recognition mecha nisms, as they show that females from different strains genuinely derive from populations with a different genetic background, and also for studies on sex allocation strategies, as consistent between strain polymorphisms allow the molecular-genetic detection of haploid (male) and diploid (female) eggs of crossmated mothers.

Table 1.

Primers designed for Goniozus legneri.


Table 2.

Primers designed for Goniozus nephantidis.


Table 3.

The six polymorphic primers for Goniozus legneri strains.



We thank Yoanna Nabalón for help with field collections, Nariman Ahmad for molecular help, and Julietta Marquez and Alda Romero for assistance with insect rearing. This work was supported by Santander Bank, who provided funds for I. C. W. Hardy to travel to Chile for field collections, and by an Iraqi government Ph.D. studentship awarded to S. K. Khidr.




simple sequence repeat



J Abe , Y Kamimura , M Shimada , SA. West 2009. Extremely female-biased primary sex ratio and precisely constant male production in a parasitoid wasp Melittobia. Animal Behaviour 78: 515–523. Google Scholar


A Aebi , T Shani , C Hansson , J Contreras-Garduno , G Mansion , B. Benrey 2008. The potential of native parasitoids for the control of Mexican bean beetles: A genetic and ecological approach. Biological Control 47: 289–297. Google Scholar


DM Althoff , JN. Thompson 2001. Geographic structure in the searching behavior of a specialist parasitoid: combining molecular and behavioral approaches. Journal of Evolutionary Biology 14: 406–417. Google Scholar


C Anton , J Settele , W. Durka 2006. Nine polymorphic microsatellite loci for the parasitic wasp Neotypus melanocephalus (Hymenoptera: Ichneumonidae). Molecular Ecology Notes 6: 399–401. Google Scholar


C Anton , I Zeisset , M Musche , W Durka , JJ Boomsma , J. Settele 2007. Population structure of a large blue butterfly and its specialist parasitoid in a fragmented landscape. Molecular Ecology 16: 3828–3838. Google Scholar


SS Arcot , Z Wang , JL Weber , PL Deininger , MA. Batzer 1995. Alu repeats: a source for the genesis of primate microsatellites. Genomics 29: 136–144. Google Scholar


JC. Avise 1994. Molecular Markers, Natural History and Evolution. Chapman & Hall. Google Scholar


DA Baker , HD Loxdale , OR. Edwards 2003. Genetic variation and founder effects in the parasitoid wasp, Diaeretiella rapae (M'intosh) (Hymenoptera: Braconidae: Aphidiidae), affecting its potential as a biological control agent. Molecular Ecology 12: 3303–3311. Google Scholar


T Barbara , C Palma-Silva , GM Paggi , F Bered , MF Fay , C. Lexer 2007. Cross-species transfer of nuclear microsatellite markers: potential and limitations. Molecular Ecology 16: 3759–3767. Google Scholar


TP Batchelor , ICW Hardy , JF Barrera , G. Pérez-Lachaud 2005. Insect gladiators II: Competitive interactions within and between bethylid parasitoid species of the coffee berry borer, Hypothenemus hampei (Coleoptera: Scolytidae). Biological Control 33: 194–202. Google Scholar


T Bentley , TT Hull , ICW Hardy , M. Goubault 2009. The elusive paradox: owner-intruder roles, strategies and outcomes in parasitoid contests. Behavioral Ecology 20: 296–304. Google Scholar


LW Beukeboom , BJ. Zwaan 2005. Genetics. In: MA Jervis , Editor. Insects as Natural Enemies: a Practical Perspective. pp 167–218. Springer. Google Scholar


J Bond , A. Stockman 2008. An integrative method for delimiting cohesion species: finding the population-species interface in a group of Californian trapdoor spiders with extreme genetic divergence and geographic structuring. Systematic Biology 57: 628–646. Google Scholar


G Buczkowski , EL Vargo , J. Silverman 2004. The diminutive supercolony: the Argentine ants of the southeastern United States. Molecular Ecology 13: 2235–2242. Google Scholar


M Carr , J Peter , W Young , PJ. Mayhew 2010. Phylogeny of bethylid wasps (Hymenoptera: Bethylidae) inferred from 28S and 16S rRNA genes. Insect Systematics and Evolution 41:55–73. Google Scholar


GR. Carvalho 1998. Advances in Molecular Ecology. IOS Press.. Google Scholar


JM Cook 1993. Experimental test of sex determination in Goniozus nephantidis. Heredity 71: 130–137. Google Scholar


E. Dharmaraju 1963. Biological control of coconut leaf caterpillar (Nephantis serinopa Meyrick) in Ceylon. Bulletin of the Coconut Research Institute Ceylon 21: 1–46. Google Scholar


SR Diehl , GL. Bush 1984. An evolutionary and applied perspective of insect biotypes. Annual Review of Entomology 29: 471–483. Google Scholar


J Drescher , N Büthgen , T Schmitt , J Bühler , H. Feldhaar 2010. Societies drifting apart? Behavioural, genetic and chemical differentiation between supercolonies in the yellow crazy ant Anoplolepis gracilipes. PLoS One 5(10): e13581. doi: 10.1371/journal.pone.0013581 Google Scholar


J Elias , S Dorn , D. Mazzi 2010. No evidence for increased extinction proneness with decreasing effective population size in a parasitoid with complementary sex determination and fertile diploid males. BMC Evolutionary Biology 10: 366. Google Scholar


S EL-Showk , JS van Zweden , P d'Ettorre , L. Sundström 2010. Are you my mother? Kin recognition in the ant Formica fusca. Journal of Evolutionary Biology 23: 397–406. Google Scholar


HE. Evans 1964. A synopsis of the American Bethylidae (Hymenoptera: Aculeata). Bulletin of the Museum of Comparative Zoology. 132: 1–122. Google Scholar


KL. Fiala 1980. On estimating the primary sex ratio from incomplete data. American Naturalist 115: 442–444. Google Scholar


A Gardner , SA. West 2007. Social evolution: the decline and fall of genetic kin recognition. Current Biology 17: R810–R812. Google Scholar


S Garrido , L Cichón , D Fernádez , C. Azevedo 2005. Primera cita de la especie Goniozus legneri (Hymenoptera: Bethylidae) en el Alto Valle de Río Negro, Patagonia Argentina. Revista del la Societad de Entomología de Argentina 64: 14–16. Google Scholar


T Giraud , JS Pedersen , L. Keller 2002. Evolution of supercolonies: the Argentine ants of southern Europe. Proceedings of the National Academy of Sciences USA 99: 6075– 6079. Google Scholar


HCJ. Godfray 1994. Parasitoids: Behavioral and Evolutionary Ecology. Princeton University Press. Google Scholar


MA Goodisman , RW Matthews , RH. Crozier 2001. Hierarchical genetic structure of the introduced wasp Vespula germanica in Australia. Molecular Ecology 10: 1423–1432. Google Scholar


G Gordh 1982. A new species of Goniozus (Hymenoptera: Bethylidae) imported into California for the biological control of the navel orangeworm (Lepidoptera: Pyralidae). Entomological News 93: 136–138. Google Scholar


G Gordh , L. Móczár 1990. A Catalog of the world Bethylidae (Hymenoptera: Aculeata). Memoirs of the American Entomological Institute 46: 1–364. Google Scholar


G Gordh , JB Woolley , RA. Medeved 1983. Biological studies on Goniozus legneri Gordh (Hymenoptera: Bethylidae), primary external parasite of the navel orangeworm Amyelois transitella and pink bollworm Pectinophora gossypiella (Lepidoptera: Pyralidae, Gelechiidae). Contributions of the American Entomological Institute 20: 433–468. Google Scholar


S Gothilf , M. Mazor 1987. Release and recovery of imported parasites of the carob moth, Spectorbates ceratonia (Lepidoptera: Pyralidae) in Israel. Israel Journal of Entomology 21: 19–23. Google Scholar


M Goubault , TP Batchelor , RST Linforth , AJ Taylor , ICW. Hardy 2006. Volatile emission by contest losers revealed by realtime chemical analysis. Proceedings of the Royal of Society of London B 273: 2853–2859. Google Scholar


M Goubault , AFS Mack , ICW. Hardy 2007. Encountering competitors reduces clutch size and increases offspring size in a parasitoid with female-female fighting. Proceedings of the Royal of Society of London B 274: 2571– 2577. Google Scholar


RF Green , G Gordh , BA. Hawkins 1982. Precise sex ratios in highly inbred parasitic wasps. American Naturalist 120: 653–665. Google Scholar


WD. Hamilton 1964a. The genetic evolution of social behaviour. II. Journal of Theoretical Biology 7: 17–52. Google Scholar


WD. Hamilton 1964b. The genetic evolution of social behaviour. I. Journal of Theoretical Biology 7: 1–16. Google Scholar


WD. Hamilton 1967. Extraordinary sex ratios. Science 156: 477–488. Google Scholar


JM Hancock 1999. Microsatellites and other simple sequences: genomic context and mutational mechanisms. In: DB Goldstein , C Schlotterer , Editors. Microsatellites Evolution and Applications. pp. 1–9. Oxford University Press. Google Scholar


JM Hancock , M. Simon 2005. Simple sequence repeats in proteins and their significance for network evolution. Genetics 345: 113–118. Google Scholar


ICW. Hardy 1994. Sex ratio and mating structure in the parasitoid Hymenoptera. Oikos 69: 3–20. Google Scholar


ICW Hardy , TM. Blackburn 1991. Brood guarding in a bethylid wasp. Ecological Entomology 16: 55–62. Google Scholar


ICW Hardy , JM. Cook 1995. Brood sex ratio variance, developmental mortality and virginity in gregarious parasitoid wasp. Oecologia 103: 162–169. Google Scholar


ICW Hardy , LJ Dijkstra , JEM Gillis , PA. Luft 1998. Patterns of sex ratio, virginity and developmental mortality in gregarious parasitoids. Biological Journal of the Linnean Society 64: 239–270. Google Scholar


ICW Hardy , M Goubault , TP. Batchelor 2013. Hymenopteran contests and agonistic behaviour. In: ICW Hardy , M Briffa , Editors. Animal Contests. pp 147–177. Cambridge University Press. Google Scholar


ICW Hardy , PJ. Mayhew 1998. Sex ratio, sexual dimorphism and mating structure in bethylid wasps. Behavioral Ecology and Sociobiology 42: 383–395. Google Scholar


ICW Hardy , JB Pedersen , MK Sejr , UH. Linderoth 1999. Local mating, dispersal and sex ratio in a gregarious parasitoid wasp. Ethology 105: 57–72. Google Scholar


ICW Hardy , S Stokkebo , J Bønløkke-Pedersen , MK. Sejr 2000. Insemination capacity and dispersal in relation to sex allocation decisions in Goniozus legneri (Hymenoptera: Bethylidae): why are there more males in larger broods? Ethology 106: 1021–1032. Google Scholar


HJ. Henter 2003. Inbreeding depression and haplodiploidy: experimental measures in a parasitoid and comparisons across diploid and haplodiploid insect taxa. Evolution 57: 1793– 1803. Google Scholar


ME Hochberg , AR. Ives 2000. Parasitoid Population Biology. Princeton University Press. Google Scholar


KR Hopper , RT Roush , W. Powell 1993. Management of genetics of biological control introductions. Annual Review of Entomology 38: 27–51. Google Scholar


RA Hufbauer , SM Bogdanowicz , RG. Harrison 2004. The population genetics of a biological control introduction: mitochondrial DNA and microsatellite variation in native and introduced populations of Aphidus ervi, a parasitoid wasp. Molecular Ecology 13: 337– 348. Google Scholar


EL Humphries , AJ Hebblethwaite , TP Batchelor , ICW. Hardy 2006. The importance of valuing resources: host weight and contender age as determinants of parasitoid wasp contest outcomes. Animal Behaviour 72: 891–898. Google Scholar


J Jaquiery , V Vogel , L. Keller 2005. Multilevel genetic analyses of two European supercolonies of the Argentine ant, Linepithema humile. Molecular Ecology 14: 589–598. Google Scholar


MA. Jervis 2005. Insects as Natural Enemies: a Practical Perspective. Springer. Google Scholar


M Kankare , S van Nouhuys , I. Hanski 2005. Genetic divergence among host-specific cryptic species in Cotesaia melitaearum aggregate (Hymenoptera: Braconidae), parasitoids of checkerspot butterflies. Annals of the Entomological Society of America 98: 382–394. Google Scholar


NG Kavallieratos , Ž Tomanović , P Starý , CG Athanassiou , GP Sarlis , O Petrović , M Niketić , M. Anagnou-Veroniki 2004. A survey of aphid parasitoids (Hymenoptera: Braconidae: Aphidiinae) of Southeastern Europe and their aphid-plant associations. Applied Entomology and Zoology 39: 527– 563. Google Scholar


SK Khidr , S Mayes , ICW. Hardy 2013. Primary and secondary sex ratios in a gregarious parasitoid with local mate competition. Behavioral Ecology 24(2): 435– 443. Google Scholar


JM Kloda , PDG Dean , DM MacDonald , S. Mayes 2004. Isolation and characterisation of microsatellite loci in Ononis repens, Leguminosae. Molecular Ecology Notes 4: 596–598. Google Scholar


S Krackow , M. Neuhäuser 2008. Insights from complete-incomplete brood sex-ratio disparity. Behavioral Ecology and Sociobiology 62: 469–477. Google Scholar


B Lavandero , CC Figueroa , P Franck , A. Mendez 2011. Estimating gene flow between refuges and crops: a case study of the bological control of Eriosoma lanigerum by Aphelinus mali in apple orchards. PLoS One 6(11): e26694. doi: 10.1371/journal.pone.0026694Google Scholar


EF Legner , A. Silveira-Guido 1983. Establishment of Goniozus emigratus and Goniozus legneri (Hym: Bethylidae) on navel orangeworm, Amyelois transitella (Lep: Phycitidae) in California and biological control potential. Entomophaga 28: 97–106. Google Scholar


A Lizé , D Carval , AM Cortesero , S Fournet , D. Poinsot 2006. Kin discrimination and altruism in the larvae of a parasitoid insect. Proceedings the Royal Society of London B 273(1599): 2381–2386. Google Scholar


A Lizé , SK Khidr , ICW. Hardy 2012. Two components of kin recognition influence parasitoid aggression in resource competition. Animal Behaviour 83: 793–799. Google Scholar


HD Loxdale , G. Lushai 1998. Molecular markers in entomology. Bulletin of Entomological Research 88: 577–600. Google Scholar


JD Lozier , NJ Mills , GK. Roderick 2006. Diand trinucleotide repeat microsatellites for the parasitoid wasp, Aphidius transcaspicus. Molecular Ecology Notes 6: 27–29. Google Scholar


JD Lozier , GK Roderick , NJ. Mills 2009. Molecular markers reveal strong geographic, but not host associated, genetic differentiation in Aphidius transcaspicus, a parasitoid of the aphid genus Hyalopterus. Bulletin of Entomological Research 99: 83–96. Google Scholar


KR Lyla , SP Beevi , T. Venkatesan 2006. Field evaluation of Goniozus nephantidis (Muesebeck) against coconut black-headed caterpillar in Kerala using different release techniques. Biological Control 20: 33–36. Google Scholar


WS Martins , DCS Lucas , KN Fabricio de Souza , DJ. Bertioli 2009. WebSat. A web software for microsatellite marker development. Bioinformation 3: 282–283. Available online: Google Scholar


JM. Mateo 2004. Recognition systems and biological organization: the perception component of social recognition. Annales Zoologici Fennici 41: 729–745. Google Scholar


D Mazzi , F Hatt , S Hein , S. Dorn 2011. Ladies last: diel rhythmicity of adult emergence in a parasitoid with complementary sex determination. Physiological Entomology 36: 47– 53. Google Scholar


KD McCoy , T Boulinier , C Tirard , Y. Michalakis 2001. Host specificity of a generalist parasite: genetic evidence of sympatric host races in seabird tick Ixodes uriae. Journal of Evolutionary Biology 14: 395–405. Google Scholar


SBJ. Menken 1981. Host races and sympatric speciation in small ermine moths, Yponomeutidae. Entomologia Experimentalis et Applicata 30: 280–292. Google Scholar


D Molbo , CA Machado , JG Sevenster , L Keller , EA. Herre 2003. Cryptic species of figpollinating wasps: implications for the evolution of the fig-wasp mutualism, sex allocation, and precision of adaptation. Proceedings of the National Academy of Sciences USA 100: 5867–5872. Google Scholar


JA Nicholls , P Fuentes-Utrilla , A Hayward , G Melika , G Csóka , J-L Nieves-Aldrey , J Pujade-Villar , M Tavakoli , K Schönrogge , GN. Stone 2010. Community impacts of anthropogenic disturbance: natural enemies exploit multiple routes in pursuit of invading herbivore hosts. BMC Evolutionary Biology 10: 322. Google Scholar


F Nyabuga , H Loxdale , DG Heckel , WW. Weisser 2010. Spatial population dynamics of a specialist aphid parasitoid, Lysiphlebus hirticornis Mackauer (Hymenoptera: Braconidae: Aphidiinae): evidence for philopatry and restricted dispersal. Heredity 105: 433–442. Google Scholar


BA Pannebakker , NRT Garrido , BJ Zwaan , JJM. van Alphen 2008. Geographic variation in host-selection behaviour in the Drosophila parasitoid Leptopilina clavipes. Entomologia Experimentalis et Applicata 127: 48–54. Google Scholar


CB Phillips , CJ Vink , A Blanchet , KA. Hoelmer 2008. Hosts are more important than destinations: what genetic variation in Microctonus aethiopoides (Hymenoptera: Braconidae) means for foreign exploration for natural enemies. Molecular Phylogenetics and Evolution 49: 467–476. Google Scholar


FLW. Ratnieks 1991. The evolution of genetic cue diversity in social Hymenoptera. American Naturalist 137: 202–226. Google Scholar


FLW Ratnieks , L. Keller 1998. Queen control of egg fertilization in the honey bee. Behavioral Ecology and Sociobiology 44: 57–61. Google Scholar


S Rozen , H. Skaletsky 2000. Primer3 on the www for general users and for biologist programmers. Methods in Molecular Biology 132: 365–386. Available online: Google Scholar


L Ruiz-Montoya , J Nunez-Farfan , J. Vargas 2003. Host-associated genetic structure of Mexican populations of the cabbage aphid Brevicoryne brassicae L. (Homoptera: Aphididae). Heredity 91: 415–421. Google Scholar


J Sambrook , EF Fritsch , T. Maniatis 1989. Molecular Cloning: A Laboratory Manual , 2nd Edition. Cold Spring Harbor Laboratory Press. Google Scholar


M. Schuelke 2000. An economic method for the fluorescent labelling of PCR fragments. Nature Biotechnology 18(2): 233–234. Google Scholar


WM. Shields 1993. The natural and unnatural history of inbreeding and outbreeding. In: NW Thornhill , Editor. The Natural History of Inbreeding and Outbreeding Theoretical and Empirical Perspectives. pp. 143–169. University of Chicago Press. Google Scholar


SA Steffan , KM Daane , DL. Mahr 2001. 15N-enrichment of plant tissue to mark phytophagous insects, associated parasitoids, and flower-visiting entomophaga. Entomologia Experimentalis et Applicata 98: 173–180. Google Scholar


JO Stireman , JD Nason , S Heard , JM. Seehawer 2006. Cascading host-associated genetic differentiation in parasitoids of phytophagous insects. Proceedings of the Royal Society of London B 273: 523–530. Google Scholar


D Tautz 1999. Hypervariability of simple sequences as a general source for polymorphic DNA markers. Nucleic Acid Research 17: 6463–6471. Google Scholar


JN. Thompson 1994. The Coevolutionary Process. University of Chicago Press. Google Scholar


A Trindl , J Heinze , P. D'Ettorre 2004. Isolation and characterization of five microsatellite loci in the ponerine ant Pachycondyla inversa (Hymenoptera, Formicidae). Molecular Ecology Notes 4: 583–585. Google Scholar


ND Tsutsui , AV Suarez , RK. Grosberg 2003. Genetic diversity, asymmetrical aggression, and cooperation in a widespread invasive species. Proceedings of the National Academy of Sciences USA 100: 1078–1083. Google Scholar


LV Ugelvig , FP Drijfhout , DJC Kronauer , JJ Boomsma , JS Pedersen , S. Cremer 2008. The introduction history of invasive garden ants in Europe: integrating genetic, chemical and behavioural approaches. BMC Biology 6: 11. Google Scholar


TR Unruh , G Gordh , D. Gonzalez 1984. Electrophoretic studies on parasitic Hymenoptera and implications for biological control. International Congress of Entomology Proceedings 17: 705. Google Scholar


T Venkatesan , SK Jalali , K Srinivasmurthy , RJ Rabindra , CB. Dasan 2007. Economics of production of Goniozus nephantidis (Muesebeck), an important parasitoid of coconut black-headed caterpillar, Opisina arenosella (Walker) for bio-factories. Biological Control 21: 53–58. Google Scholar


AP Vogler , R. Desalle 1993. Phylogeographic patterns in coastal North American tiger beetles (Cicindela dorsalis Say) inferred from mitochondrial DNA sequences. Evolution 47: 1192–1202. Google Scholar


E Wajnberg , C Bernstein , JJM. van Alphen 2008. Behavioral Ecology of Insect Parasitoids: from Theoretical Approaches to Field Applications. Blackwell Publishing. Google Scholar


Z Wu , KR Hopper , RJ O'Neil , DJ Voegtlin , DR Prokrym , GE. Heimpel 2004. Reproductive compatibility and genetic variation between two strains of Aphelinus albipodus (Hymenoptera: Aphelinidae), a parasitoid of the soybean aphid, Aphis glycines (Homoptera: Aphididae). Biological Control 31: 311–319. Google Scholar


T Zaviezo , A Romero , D Castro , A. Wagner 2007. First record of Goniozus legneri (Hymenoptera: Bethylidae) in Chile. Ciencia e Investigación Agraria 34: 49–52. Google Scholar


M Zavodna , P Arens , PJ van Dijk , T Partomihardjo , B Vosman , JM. van Damme 2005. Pollinating fig wasps: genetic consequences of island recolonization. Journal of Evolutionary Biology 18: 1234– 1243. Google Scholar


RM. Zink 2002. Methods in comparative phylogeography, and their application to studying evolution in the North American Aridlands. Integrative and Comparative Biology 42: 953–959. Google Scholar


NE Zygouridis , AA Augustinos , FG Zalom , KD. Mathiopoulos 2009. Analysis of olive fly invasion in California based on microsatellite markers. Heredity 102: 402–412. Google Scholar
This is an open access paper. We use the Creative Commons Attribution 3.0 license that permits unrestricted use, provided that the paper is properly attributed.
Sahand K. Khidr, Ian C.W. Hardy, Tania Zaviezo, and Sean Mayes "Development of Microsatellite Markers and Detection of Genetic Variation between Goniozus Wasp Populations," Journal of Insect Science 14(43), 1-17, (1 March 2014).
Received: 23 June 2012; Accepted: 14 March 2013; Published: 1 March 2014

Back to Top