<title>Background</title> <p>The complex, interactive communities displayed by highly social insects represent an important and highly successful evolutionary innovation [##UREF##0##1##,##UREF##1##2##]. The success of social insects arises largely from their extraordinary cooperative and helping behaviors [##UREF##1##2##, ####UREF##2##3##, ##UREF##3##4##, ##UREF##4##5####4##5##]. However, the intricate social organization displayed by social insects may also be exploited and manipulated. For instance, many social insects are subject to social parasitism [##REF##21227369##6##,##UREF##5##7##]. Social parasites benefit from brood care or other resources at the expense of the society of a social host species [##UREF##6##8##].</p> <p>In hymenopteran social insects (ants, some bees, and some wasps), social parasitism may take several forms [##UREF##2##3##,##UREF##5##7##, ####UREF##6##8##, ##UREF##7##9##, ##UREF##8##10####8##10##]. The most extreme social parasites are completely dependent on their host taxa. Queens of these obligate parasites enter active host colonies, kill the resident queen, and use the remaining worker force of the host to rear their own parasitic offspring. These obligate parasites do not produce their own workers, and are completely reliant on the workers of the host to complete their life cycle. Alternatively, some social parasites display facultatively parasitic behavior. Queens of these species may usurp colonies of their hosts in some cases, but may also reproduce independently under other conditions. These facultative social parasites are of particular interest, because they may represent an intermediate stage in the evolution of socially parasitic behavior [##UREF##9##11##,##UREF##10##12##].</p> <p>The social wasp <italic>Vespula squamosa</italic>, commonly known as the southern yellowjacket, is a facultative social parasite [##UREF##9##11##,##UREF##11##13##, ####UREF##12##14##, ##UREF##13##15####13##15##]. <italic>Vespula squamosa </italic>is found throughout the eastern part of the United States extending south to Honduras (Fig. ##FIG##0##1##) [##UREF##14##16##,##UREF##15##17##]. It is thought that <italic>V. squamosa </italic>queens can found colonies independently under some circumstances [##UREF##9##11##], because no known hosts live in the southern part of the range of <italic>V. squamosa</italic>. However, throughout most of its range, <italic>V. squamosa </italic>is considered to be parasitic [##UREF##11##13##].</p> <p><italic>Vespula squamosa </italic>queens parasitize host taxa by usurping active nests established by other queens [##UREF##9##11##,##UREF##10##12##,##UREF##12##14##,##UREF##13##15##,##UREF##16##18##]. <italic>Vespula squamosa </italic>queens are known to primarily parasitize species in their own genus. However, there are reports of <italic>V. squamosa </italic>queens usurping colonies of distantly related taxa, such as the hornet <italic>Vespa crabro </italic>[##UREF##16##18##]. Where parasitism occurs, <italic>V. squamosa </italic>queens emerge relatively late in the season and seek out already established host colonies to usurp [##UREF##10##12##, ####UREF##11##13##, ##UREF##12##14####12##14##]. Once host colonies are located, the <italic>V. squamosa </italic>queen kills the resident queen and assumes possession of the colony. The remaining host workers help the <italic>V. squamosa </italic>queen rear her own worker and sexual offspring. Eventually the usurped workers die and the entire colony comes to be inhabited by <italic>V. squamosa </italic>individuals.</p> <p>In the southeastern part of the United States, the principal host of <italic>V. squamosa </italic>is the eastern yellowjacket, <italic>V. maculifrons </italic>(Fig. ##FIG##0##1##) [##UREF##12##14##,##UREF##13##15##]. In some areas, 40% of <italic>V. maculifrons </italic>colonies fall victim to <italic>V. squamosa </italic>parasitism [##UREF##17##19##]. Moreover, 80% of <italic>V. squamosa </italic>colonies show clear evidence of having originated by parasitism of <italic>V. maculifrons </italic>colonies in this part of the country. In fact, the frequency of parasitism is likely higher, because usurped nests can only be detected if nest take-over occurs subsequent to significant cell construction by the host [##UREF##13##15##].</p> <p>Social parasitism may play an important role in determining the levels of genetic variation and structure of the two interacting taxa [##UREF##5##7##,##REF##11308089##20##, ####REF##11308088##21##, ##UREF##18##22####18##22##]. For example, <italic>V. squamosa </italic>queens represent a significant mortality factor for <italic>V. maculifrons</italic>, owing to the relatively high rate of nest takeover. Nest parasitism may thereby potentially depress the amount of genetic variation in <italic>V. maculifrons </italic>populations by limiting <italic>V. maculifrons </italic>population size. Moreover, the number of <italic>V. squamosa </italic>nesting sites is potentially constrained by the presence of <italic>V. maculifrons </italic>colonies in the population, because <italic>V. squamosa </italic>queens primarily found nests by usurping those already initiated by <italic>V. maculifrons</italic>. Evolutionary interactions of social parasites are also known to be influenced by rates of migration for a parasite and its host [##UREF##5##7##]. The interaction between migration and local adaptation may lead to a coevolutionary arms race between parasites and their hosts with gene flow as the primary currency [##UREF##19##23##]. In sum, parasite range and effective population size are related to the degree of population structure of its host. However, the degree that this will influence demography and genetic structure is currently unknown in a facultative social parasite. Thus, our aim was to understand the genetic structure and levels of genetic variation in <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>to determine if the parasitic lifestyle differentially affected gene flow and population size in the two taxa.</p> <p>The relationship between a social parasite and its host also likely affects the breeding systems of both taxa. For instance, signaling systems between hosts and parasites provide an interesting paradox associated with genetic diversity generated by the breeding system of the host [##UREF##8##10##]. Increased within-colony genetic diversity of a host species arising as a result of multiple queens within colonies or multiple mating by queens may lead to a superior defense against pathogens or enhanced division-of-labor (e.g., [##REF##15218093##24##, ####UREF##20##25##, ##UREF##21##26##, ##UREF##22##27####22##27##]). However, high levels of within-colony genetic diversity may also enhance the ability of social parasites to invade host colonies, because it leads to a greater diversity of recognition cues [##UREF##23##28##,##UREF##24##29##]. In addition, the breeding system of the parasite may also be affected by the host-parasite interaction. Specifically, Sumner <italic>et al</italic>. [##REF##14999273##30##] recently found that parasitism may affect mate number of the parasitic queen. They discovered that queens of a socially parasitic ant mated singly, whereas queens of the closely-related host species mated multiply. These results suggested that benefits to multiple mating were only accrued in a free-living lifestyle. In contrast, the obligate parasite reverted to single mating, because multiple mating apparently did not provide benefits to the parasitic lifestyle.</p> <p>In <italic>Vespula </italic>species, within colony genetic diversity is directly related to queen mate number because annual <italic>Vespula </italic>colonies are always headed by a single queen. Indeed, queens of all <italic>Vespula </italic>taxa mate with multiple males [##UREF##25##31##, ####UREF##26##32##, ##REF##11412365##33##, ##REF##16050107##34##, ##UREF##27##35##, ##REF##17561915##36##, ##REF##17767595##37####17767595##37##]. We thus aimed to investigate if the host-parasite relationship between <italic>V. squamosa </italic>and <italic>V. maculifrons </italic>altered mate number of the parasite <italic>V. squamosa </italic>relative to its host <italic>V. maculifrons</italic>. In this system, however, because <italic>V. squamosa </italic>colonies exist as free-living entities for most of their life cycle (i.e., after usurpation and colony takeover is complete), we predicted that <italic>V. squamosa </italic>would not have reverted to a mating system with reduced number of mates, as has been found in other obligate social parasites.</p> <p>Overall, the primary objectives of this study were to compare levels of genetic structure and breeding systems of a social wasp and its social parasite. We expected that the parasitic interaction would lead to differences in the genetic structures, but not the breeding systems, of the two species. We conclude this study by discussing how the observed patterns of genetic variation provide insight into how parasitic interactions affect genetic structure within and between species.</p>
<title>Methods</title> <title>Samples</title> <p>We collected 37 <italic>V. maculifrons </italic>and 13 <italic>V. squamosa </italic>colonies in and around the city of Atlanta, GA, USA (Fig. ##FIG##1##2##). Colonies were anesthetized with ether, extracted from the ground, and brought back to the lab. Several workers from each colony were then placed in 95% ethanol for subsequent genetic analysis. Colony collection occurred late in the season (September – November), when usurpation would have been complete, so we did not find mixed colonies containing workers of both species. In addition, most colonies had initiated the reproductive phase of their development and were already producing new queens and males.</p> <title>Microsatellite loci</title> <p>We determined the utility of several previously developed microsatellite loci in providing genetic information in <italic>V. maculifrons </italic>and <italic>V. squamosa</italic>. Specifically, we investigated if the loci cloned by Thorèn <italic>et al. </italic>[##REF##11348506##38##] in <italic>V. rufa</italic>, Hasegawa and Takahashi [##UREF##28##39##] in <italic>Vespa mandarinia</italic>, and Daly <italic>et al. </italic>[##UREF##29##40##] in <italic>V. vulgaris </italic>were genetically informative in our study taxa. For these assays, we determined if each locus would PCR amplify in <italic>V. maculifrons </italic>and <italic>V. squamosa</italic>. If a locus amplified, we then used agarose gel electrophoresis to determine if the locus displayed size variation in a subset of individuals. In the end, we obtained the genotype of all sampled workers of both species at the eight loci, LIST2003, LIST2004, LIST2013, LIST2019, LIST2020, Rufa 5, VMA-3, and VMA-6, which amplified and were polymorphic in both taxa.</p> <title>Estimates of variability</title> <p>Members of social insect colonies do not represent independent genetic samples because they are related. Consequently, it is necessary to remove the effects of relatedness among colony members to effectively analyze population level data. We avoided the problem of genetic nonindependence of colonymates by inferring the genotypes of the queen and male(s) that produced workers from each colony. Paternity analyses in hymenopteran taxa, such as <italic>Vespula</italic>, are particularly straightforward because males are haploid and full siblings always display the same multilocus haplotype derived from their father. Our analyses thus resulted in a set of diploid (queen) and haploid (male) individuals that were genetically independent. We used the genotypes of these inferred parental individuals to obtain estimates of genetic variability.</p> <p>Based on the reconstructed male and queen genotypes, we calculated the allelic richness at the eight loci in <italic>V. maculifrons </italic>and <italic>V. squamosa</italic>. In addition, we used the rarefaction method of Petit <italic>et al. </italic>[##UREF##30##41##] to correct for the difference in the number of genes sampled in the two species. This method estimates the number of alleles that would be observed in an equal sample of genes from multiple groups based on the number of alleles observed in the actual unequal samples obtained. The use of the standardized allele richness allows for comparisons in the diversity of samples of different sizes. We then calculated the effective number of alleles and gene diversity of each locus to obtain general measures of the overall variability for the markers.</p> <title>Population structure</title> <p>To avoid problems caused by nonindependence of workers sampled from the same colony, we used a resampling technique that yielded unbiased measures of population genetic structure [##REF##11412365##33##]. Specifically, we used reduced data sets that included only one individual per colony for population analyses. Briefly, a computer program was written to randomly select a single individual's multilocus genotype from each colony, yielding a new data set with the number of individuals sampled equal to the number of colonies. This procedure was repeated 25 times to produce 25 such data sets. Each of these data sets was then used to calculate the population statistic of interest, and the median of the 25 values was taken as the unbiased estimate. Probability tests implemented by the program GENEPOP 3.4 [##UREF##31##42##] were used in conjunction with this resampling procedure to determine the significance of genotypic disequilibrium between microsatellite markers and deviations from Hardy-Weinberg equilibrium.</p> <p>We next determined if <italic>V. maculifrons </italic>or <italic>V. squamosa </italic>populations showed evidence of isolation by distance over the geographical scale considered. We first calculated estimates of <italic>F</italic><sub>ST </sub>between all pairs of colonies for each species using GENEPOP. Thus, colonies, which consisted of groups of sampled workers, were considered as 'populations' for these analyses. Significance of the correlation between geographical and genetic distance was assessed with a Mantel test. The relationship was deemed to be significant if the observed correlation was greater than 5% of the 10,000 randomly generated values. Spearman's rank order correlation coefficient (<italic>r</italic><sub>S</sub>) was used to quantify the association between the genetic distances between nests and the geographical distances that separated them.</p> <title>Mating system</title> <p>We used the molecular genetic data derived from the worker genotypes to investigate how the breeding systems of <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>differed. Because annual <italic>Vespula </italic>colonies are headed by single queens, within colony genetic diversity provides insight into the queen's mating history. We thus used the genetic data to determine the number of males to which queens were mated. We then compared mate number in the two taxa via a Kruskal-Wallis test. We next calculated the effective paternity (<italic>k</italic><sub><italic>e</italic>3</sub>) for workers using the method of Nielsen <italic>et al. </italic>[##REF##14629394##43##]. The metric <italic>k</italic><sub><italic>e</italic>3 </sub>includes information on the number of times a queen mates and the unequal contributions of a queen's male mates to offspring. We also tested if the estimates of <italic>k</italic><sub><italic>e</italic>3 </sub>obtained for <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>differed using a Kruskal-Wallis test.</p> <p>Additionally, we investigated if the distribution of mate number in the two species matched known distributions, as might be expected if particular biological processes (such as constant mate search time or equal probability of mating with particular numbers of males) affected the mating behavior of queens [##UREF##32##44##]. Specifically, we fitted a Gaussian distribution to queen mate number. We expected that mate number would match this distribution if number of mates were a random function of time and mating opportunity. In contrast, deviations of mate number from the Gaussian distribution might indicate the action of selection operating in these species.</p> <p>We calculated the magnitude of reproductive skew of males mated to queens using the metric <italic>B </italic>of Nonacs [##UREF##33##45##]. The 95% confidence intervals for each estimate of <italic>B </italic>were estimated to determine the significance of skew. Values of <italic>B </italic>were considered to be significant if the 95% confidence intervals failed to overlap 0. Both the skew and significance of skew were calculated using the program Skew Calculator [##UREF##33##45##]. As was the case with our estimates of <italic>k</italic><sub><italic>e</italic>3</sub>, we tested if estimates of <italic>B </italic>differed between species via a Kruskal-Wallis test.</p> <p>We calculated the relatedness of workers belonging to the same colony using the program RELATEDNESS [##UREF##34##46##]. We then estimated the relatedness of the queens to their male mates and the relatedness of males mated to single queens to determine if inbreeding occurred or if related males mated the same queen. Standard errors for all estimates were obtained by jackknifing over loci. Finally, we compared our relatedness estimates with similar data collected over two decades ago to determine how the parasitic relationship may have changed mating structure over time.</p>
<title>Results</title> <title>Samples</title> <p>We sampled 40.56 ± 23.51 (mean ± SD) and 45.62 ± 8.31 workers from each of 37 <italic>V. maculifrons </italic>and 13 <italic>V. squamosa </italic>colonies, respectively (Fig. ##FIG##1##2##). Our total sample size thus consisted of 1501 and 463 workers from the two taxa. The maximum distance between colonies in this study spanned approximately 120 km. Moreover, the sampling locations of <italic>V. maculifrons </italic>showed considerable overlap with those of <italic>V. squamosa </italic>(Fig. ##FIG##1##2##).</p> <title>Microsatellite loci</title> <p>We determined if the 43 microsatellite loci developed in related species within the Vespidae were informative in <italic>V. maculifrons </italic>or <italic>V. squamosa</italic>. Overall, we found that 16 loci amplified and were variable in both species (Appendix). We used eight of these 16 loci for our study (see Methods). We analyzed reduced data sets consisting of only one individual per colony to determine if loci were in Hardy-Weinberg equilibrium in the two species. We found no evidence of disequilibrium among workers in either <italic>Vespula </italic>taxon (<italic>P </italic>&gt; 0.1 for all loci in both species; overall combined <italic>P </italic>= 0.2388 in <italic>V. maculifrons </italic>and <italic>P </italic>= 0.5304 in <italic>V. squamosa</italic>). In addition, there was no evidence of linkage disequilibrium between any pair of loci (<italic>P </italic>&gt; 0.3 for all pairs of loci in both species). The eight loci that we ultimately utilized for this study showed substantial variability in the two taxa, with gene diversities (<italic>h</italic>) often exceeding 0.7 (Table ##TAB##0##1##). The exceptions were the loci Rufa-5 in <italic>V. maculifrons </italic>and LIST2019 in <italic>V. squamosa</italic>, where the variability was somewhat less. Regardless, the variability of the loci was high enough such that the probability of any two males having the same multilocus genotype (nondetection error) was extremely low (&lt;&lt; 0.0001).</p> <title>Estimates of variability</title> <p>The standardized number of alleles (<italic>A</italic><sub>50</sub>, defined as the number of alleles expected to be observed if only 50 genes had been sampled from a given population; Table ##TAB##0##1##; [##UREF##30##41##]) did not differ significantly between the taxa [12.8 ± 6.3 (mean ± SD) for <italic>V. maculifrons </italic>and 11.5 ± 4.45 for <italic>V. squamosa</italic>; paired <italic>t</italic>-test, <italic>t </italic>= 0.74, <italic>P </italic>= 0.48] suggesting that genetic variability was not substantially different between the parasite and its host. In addition, the effective number of alleles did not show strong differences between species (9.2 ± 5.90 for <italic>V. maculifrons </italic>and 7.0 ± 3.31 for <italic>V. squamosa</italic>; paired <italic>t</italic>-test, <italic>t </italic>= 1.27, <italic>P </italic>= 0.24). Moreover, the variability at particular loci was similar in the two species. Specifically, the estimates of <italic>A</italic><sub>50 </sub>and expected heterozygosity (<italic>h</italic>) were significantly correlated across loci in the two taxa (Spearman's correlation coefficient; <italic>A</italic><sub>50</sub>, <italic>r</italic><sub>S </sub>= 0.738, <italic>P </italic>= 0.0366; <italic>h</italic>, <italic>r</italic><sub>S </sub>= 0.762, <italic>P </italic>= 0.0280).</p> <title>Population structure</title> <p>The overall estimates of genetic differentiation between colonies in <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>were <italic>F</italic><sub>ST </sub>= 0.1953 and 0.1914, respectively. We found no evidence of genetic isolation by distance among colonies in either species (Fig. ##FIG##2##3##). The correlation between geographic distance and genetic distance was low and nonsignificant in <italic>V. maculifrons </italic>(<italic>F</italic><sub>ST</sub>, <italic>r</italic><sub>S </sub>= 0.0176, <italic>P </italic>= 0.4694) and <italic>V. squamosa </italic>(<italic>F</italic><sub>ST</sub>, <italic>r</italic><sub>S </sub>= 0.0036, <italic>P </italic>= 0.4748). Thus, both of these <italic>Vespula </italic>species are apparently able to mate at random over the range in which sampling occurred.</p> <title>Mating system</title> <p>As expected, we found that the genotype distributions of workers within colonies always were consistent with the presence of a single queen mated to multiple males. <italic>V. maculifrons </italic>queens mated with 5.64 ± 1.27 males (mean ± SD; range of 3 – 8, Fig. ##FIG##3##4##), whereas <italic>V. squamosa </italic>queens mated with 7.25 ± 1.86 males (range of 5 – 12; Fig. ##FIG##3##4##). The mate number of queens in the two taxa differed significantly from each other (Kruskal Wallis <italic>S </italic>= 412.5, <italic>P </italic>= 0.0037). The estimates of effective mate number (<italic>k</italic><sub><italic>e</italic>3</sub>) in <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>were 4.96 ± 1.40 and 5.58 ± 1.66, respectively (range in <italic>k</italic><sub><italic>e</italic>3 </sub>of 2.61 – 8.82 in <italic>V. maculifrons </italic>and 3.75 – 10.15 in <italic>V. squamosa</italic>; Fig. ##FIG##3##4##). In contrast to the actual mate number, the estimates of effective mate number did not differ significantly between species (<italic>S </italic>= 343, <italic>P </italic>= 0.2482).</p> <p>We then turned our attention to the reproductive skew of males mated to the same queen. The mean magnitude of skew was <italic>B </italic>= 0.0273 ± 0.0409 in <italic>V. maculifrons </italic>and <italic>B </italic>= 0.0433 ± 0.0411 in <italic>V. squamosa</italic>. These values were not significantly different from each other (S = 371, <italic>P </italic>= 0.0685). However, the estimates of skew were significantly greater than zero in only 13 of the 36 <italic>V. maculifrons </italic>colonies, whereas the estimates of skew were significant in 10 of the 12 <italic>V. squamosa </italic>colonies. Thus, the proportions of colonies in which the skew was significant in the two taxa differed (<italic>G</italic><sub>1 </sub>= 8.553, <italic>P </italic>= 0.0034). Nevertheless, although paternity skew was statistically significant in many colonies, the overall magnitude of <italic>B </italic>was low and the actual mate numbers were similar to the effective mate numbers in both species.</p> <p>We investigated if the distribution of queen mate number in the two species fit known distributions to determine if simple biological processes could explain patterns of queen mating behavior [##UREF##32##44##]. We found that the distribution of mate number in <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>differed significantly from that expected under a Gaussian distribution (Shapiro-Wilk test; <italic>V. maculifrons</italic>, <italic>W </italic>= 0.920, <italic>P </italic>= 0.0125; <italic>V. squamosa</italic>, <italic>W </italic>= 0.852, <italic>P </italic>= 0.039). The failure of these distributions to fit the data arose from an observed excess of queens that mated with an intermediate number of males (Fig. ##FIG##3##4##).</p> <p>Nestmate <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>workers had relatedness values significantly greater than zero (Table ##TAB##1##2##). However, the relatedness of queens to their male mates did not differ significantly from zero (Table ##TAB##1##2##), as judged by the fact that the 95% confidence intervals (<italic>r </italic>± 1.96 * SEM) overlapped zero. This indicated that inbreeding does not usually occur in either of these species. Furthermore, we found that relatedness of <italic>V. maculifrons </italic>males mated to single queens was not significantly different from zero (Table ##TAB##1##2##). But we did find weak evidence that relatedness among <italic>V. squamosa </italic>males mated to single females was significantly greater than zero. Much of this signal, however, was due to a single colony in which the relatedness among male mates was surprisingly high (<italic>r </italic>= 0.312).</p>
<title>Discussion</title> <p>The interaction between a parasite and its host may lead to an evolutionary arms race between the taxa. Microparasites often have large population sizes and fast generation times, and may be able to win the evolutionary arms race with their hosts [##REF##15969711##47##]. In contrast, social parasites may have reduced effective population sizes owing to their dependence on host colonies and the 1:1 ratio of parasite to host colony replacement. Moreover, the difference between the level of local adaptation between a host and its parasite is determined by the rates of gene flow within host and parasite populations [##UREF##19##23##]. Limited gene flow may lead to local coadaptation between hosts and parasites [##UREF##19##23##,##UREF##35##48##,##REF##12969483##49##]. In contrast, theoretical models predict that increasing rates of gene flow among host populations can influence parasite population demography such that population densities are reduced, potentially leading to the extinction of the parasite [##UREF##5##7##,##REF##12778544##50##].</p> <p>This study represents one of the first joint analyses of levels of genetic variation and breeding system in a parasitic social wasp and its primary host. The main focus of this study was to understand if the parasite <italic>V. squamosa </italic>displayed significant differences in its mating system, levels of genetic variation, and population genetic structure from its host <italic>V. maculifrons</italic>. Subsequently, our investigation led us to more closely compare and contrast the mating system of the two species in order to understand what factors might affect mating behavior. Overall, we found strong similarities in the genetic structure and mating system of the host, <italic>V. maculifrons</italic>, and its parasite, <italic>V. squamosa</italic>. In addition, our investigation of mating biology uncovered patterns consistent with the effects of selection operating on mate number in both taxa.</p> <title>Population genetic variation and structure</title> <p>One of the more striking results arising from our investigation of levels of genetic variation was that measures of allelic diversity were not appreciably different between <italic>V. squamosa </italic>and <italic>V. maculifrons</italic>. Population genetics theory suggests that the amount of genetic variation maintained within a population is related to the effective population size [##UREF##36##51##]. Therefore, the shared patterns of diversity among these species suggest that long-term effective population size is not substantially different in the two taxa.</p> <p>Furthermore, we failed to find evidence for genetic isolation by distance in either <italic>Vespula </italic>species. The absence of isolation by distance within both taxa suggests that these <italic>Vespula </italic>wasps exhibit substantial levels of gene flow within the sampling range of this study. This finding is consistent with the high flight capacity of <italic>Vespula </italic>wasps, which are naturally capable of dispersing over considerable distances [##UREF##37##52##,##UREF##38##53##]. In addition, hibernating <italic>Vespula </italic>queens may readily be moved via accidental human transportation thereby increasing their effective migratory abilities [##UREF##39##54##,##UREF##40##55##].</p> <p>However, our finding of a lack of geographic structure in <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>differs somewhat from that found in a previous study of the congener <italic>V. germanica </italic>in Australia [##REF##11412365##33##]. In that investigation, <italic>V. germanica </italic>was found to display genetic isolation by distance on approximately the same scale investigated in this study. We note, however, that <italic>V. germanica </italic>is a recent invader to Australia [##UREF##41##56##]. Thus, the population dynamics of <italic>V. germanica </italic>in Australia likely differ from those of <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>in their native range in the United States. In particular, the observed isolation by distance in <italic>V. germanica </italic>populations in Australia may reflect non-equilibrium conditions that will ultimately fade as gene flow swamps local genetic differentiation [##UREF##15##17##].</p> <p>The overall similarity of levels of genetic variation and low levels of structure in <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>stand in contrast to expectations derived from other studies of social hymenopteran parasites. For example, ant social parasites often have low population sizes and dispersal capabilities [##UREF##5##7##]. In accord with these predictions, Trontti et al. [##UREF##42##57##] found that an inquiline parasite species exhibited much greater genetic substructuring than was found in its host. Brandt et al. [##REF##17498232##58##] also found substantial genetic structure in two species of parasitic ants, although the magnitude of structuring did not differ substantially from that of the host species. Thus our findings indicate that social parasites need not show levels of variation distinctly different from their hosts.</p> <title>Queen mate number and mating behavior</title> <p>We compared the mating system of <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>in order to investigate if there were differences between the two species that might be associated with parasitic behaviors. As expected, we found that all colonies of both species were headed by only a single queen, and that all queens of both species were polyandrous [##UREF##27##35##, ####REF##17561915##36##, ##REF##17767595##37####17767595##37##]. Our estimates of worker relatedness for <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>(0.373 and 0.357, respectively; Table ##TAB##1##2##) did not differ significantly from earlier estimates obtained for these species (0.320 and 0.403; [##UREF##27##35##]). Thus, we found no evidence for temporal variation in the mating systems of these two taxa from samples obtained some 23 years apart. Although only tangential to this study, this result is noteworthy. Few studies have investigated temporal stability of mating structure. Here, we found that structure did not change over a 23 generation time-span. This stability in general mating structure suggests that selective factors affecting mating behavior have remained constant over that time.</p> <p>We also did not find any convincing evidence of inbreeding between queens and males of both the parasite and its host. In addition, the relatedness of males mated to the same queen tended to be low in both taxa (Table ##TAB##1##2##). These data are consistent with our understanding of the mating system of <italic>Vespula </italic>wasps. Namely, although inbreeding is possible under laboratory conditions [##UREF##43##59##], matings in natural settings tend to take place away from the nest and result in outbreeding [##UREF##16##18##,##UREF##26##32##,##UREF##44##60##]. Thus, the parasitic behavior of <italic>V. squamosa </italic>has not led to a cycle of inbreeding as occurs in some socially parasitic ants [##UREF##5##7##,##UREF##45##61##].</p> <p>Our primary interest in studying queen mate number was to test hypotheses related to how mate number should evolve in social parasites. We predicted that mate number would be similar in the two taxa, because <italic>V. squamosa </italic>colonies exist as free-living entities for a substantial part of their life cycle. Therefore, they should retain benefits of multiple mating [##REF##17767595##37##], such as superior task performance of parasite defense associated with increased within-colony genetic diversity, unlike obligate social parasites that are more intimately associated with their hosts [##REF##14999273##30##]. Interestingly, we found that queen mate number differed significantly between the species, with parasitic queens tending to have more mates than host queens. This finding may result from a regime of strong selection of parasites on hosts. In this case, hosts of social parasites might be selected to develop less genetically diverse colonies to increase colony uniformity [##UREF##5##7##] and hence increase parasite recognition. However, this is not likely to be the case in this system, as levels of within-colony genetic diversity, the currency of multiple mating, did not differ substantially between the species. Thus, overall, our results suggest that <italic>V. squamosa </italic>does indeed bear the benefits of multiple mating likely as a result of the free-living portion of its life cycle [##REF##14999273##30##].</p> <p>We also found that effective mating frequency of queens was similar between the species. This suggests that unequal use of sperm reduces the effective mate numbers to approximately equal levels in both species, despite the fact that actual queen mate number differed significantly in the two taxa. Although the magnitude of paternity skew was similar between species, such skew was significant more frequently in <italic>V. squamosa</italic>. Male skew has important implications for the reproductive success of males and, potentially, for how colony members allocate resources to the production of sexuals [##UREF##46##62##]. The underlying mechanisms resulting in the variation of male reproductive skew between the two species remain unclear at this time, but may result from female choice of male sperm or male-male competition occurring among sperm within the female reproductive tract.</p> <p>Finally, one striking feature of our data was the similarity in the shape of the distribution of mate number in both species. In particular, the distribution of the number of male mates (Fig. ##FIG##3##4##) indicated that the extremes at both ends of the distribution had been somewhat truncated. This suggests that the number of mates observed does not simply result from random interactions between queens and males. One explanation for the observed pattern is that mate number is under stabilizing selection. Specifically, mating to too few males may be maladaptive [##REF##14575328##63##], whereas mating to too many males may be costly [##UREF##47##64##]. In fact, Goodisman <italic>et al</italic>. [##REF##17767595##37##] identified a positive correlation between a fitness correlate and number of matings by queens in <italic>V. maculifrons</italic>, indicating that there is a selective advantage to polyandry. Regardless, the patterns of mating behavior found here suggest that similar factors may influence queen mate choice decisions in <italic>V. squamosa</italic>. However, queens in each species may achieve the optimal effective mate number via different evolutionary strategies. That is, <italic>V. squamosa </italic>queens mate with more males, but have greater skew in sperm usage, while <italic>V. maculifrons </italic>queens mate with fewer males, but have less skew.</p>
<title>Conclusion</title> <p>The overall similarity in population genetic makeup of <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>may reflect the fact that <italic>V. squamosa </italic>represents an intermediate stage in the evolution of social parasitism. Taylor [##UREF##10##12##] proposed that parasitism in <italic>Vespula </italic>develops through four stages. First, queens display intraspecific, facultative, temporary social parasitism, whereby conspecific queenless nests are taken over by queens searching for a nesting site. Second, a species may evolve interspecific, facultative, temporary social parasitism. In this case, a queen searching for a nest occasionally takes over the queenless nest of a heterospecific host. Third, a species may evolve to interspecific, obligatory, temporary parasitism. In this third evolutionary step, the parasitic queen loses her ability to found new colonies but still produces her own workers once she takes over an established heterospecific nest. And finally, a species may evolve to display interspecific, obligatory, permanent parasitism, where the parasitic worker caste is completely lost.</p> <p>Extreme specialists (obligate, workerless parasites) are usually rare, restricted locally, and highly genetically structured. However, the genetic fingerprint of <italic>V. squamosa </italic>may be somewhat less defined by its parasitic life style because it falls somewhere between stages two and three in Taylor's proposed series of the evolution of parasitism. In support of this hypothesis, Hölldobler and Wilson [##UREF##2##3##] point out that facultatively parasitic ants, or those that are more primitively parasitic, tend to be widely distributed, as seems to be the case with <italic>V. squamosa</italic>.</p> <p>Future research in the study of the social parasite <italic>V. squamosa </italic>should aim to determine the extent that this species is facultatively versus obligately parasitic. In particular, a study of geographic variation in parasitic behavior and frequency would be informative. Additionally, it has been suggested that social parasites that exploit multiple hosts may form host races. A previous study failed to find evidence for such host race formation in <italic>Polistes </italic>[##REF##16135131##65##]. However, closer study of <italic>V. squamosa</italic>, which parasitizes multiple hosts [##UREF##16##18##], would provide another system in which to study how parasitism can lead to increased biodiversity of social insects.</p>
<title>Background</title> <p>Social insects dominate ecological communities because of their sophisticated group behaviors. However, the intricate behaviors of social insects may be exploited by social parasites, which manipulate insect societies for their own benefit. Interactions between social parasites and their hosts lead to unusual coevolutionary dynamics that ultimately affect the breeding systems and population structures of both species. This study represents one of the first attempts to understand the population and colony genetic structure of a parasite and its host in a social wasp system.</p> <title>Results</title> <p>We used DNA microsatellite markers to investigate gene flow, genetic variation, and mating behavior of the facultative social parasite <italic>Vespula squamosa </italic>and its primary host, <italic>V. maculifrons</italic>. Our analyses of genetic variability uncovered that both species possessed similar amounts of genetic variation and failed to show genetic structure over the sampling area. Our analysis of mating system of <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>revealed high levels of polyandry and no evidence for inbreeding in the two species. Moreover, we found no significant differences between estimates of worker relatedness in this study and a previous investigation conducted over two decades ago, suggesting that the selective pressures operating on queen mate number have remained constant. Finally, the distribution of queen mate number in both species deviated from simple expectations suggesting that mate number may be under stabilizing selection.</p> <title>Conclusion</title> <p>The general biology of <italic>V. squamosa </italic>has not changed substantially from that of a typical, nonparasitic <italic>Vespula </italic>wasp. For example, population sizes of the host and its parasite appear to be similar, in contrast to other social parasites, which often display lower population sizes than their hosts. In addition, parasitism has not caused the mating behavior of <italic>V. squamosa </italic>queens to deviate from the high levels of multiple mating that typify <italic>Vespula </italic>wasps. This stands in contrast to some socially parasitic ants, which revert to mating with few males. Overall, the general similarity of the genetic structure of <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>presumably reflects the fact that <italic>V. squamosa </italic>is still capable of independent colony founding and thus reflects an intermediate stage in the evolution of social parasitism.</p>
<title>Authors' contributions</title> <p>EAH, JLK, and MADG conceived the study and participated in its design and coordination. EAH and JLK carried out laboratory analyses. EAH and MADG analyzed the data and drafted the manuscript. All authors read and approved the final manuscript.</p> <title>Appendix</title> <p>Utility of 43 microsatellite loci in V. maculifrons and V. squamosaa (Table ##TAB##2##3##)</p>
<title>Acknowledgements</title> <p>This research was supported by the Georgia Institute of Technology and the United States National Science Foundation (#DEB-0640690). We thank E. A. Matthews for help collecting wasps, D. Bhatka, J. Rekau, and K. A. Sankovich for laboratory support, and T. N. Thirer for assistance with microsatellite analysis.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Distribution of <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>in North America.</bold> The range of the parasite <italic>V. squamosa </italic>largely overlaps with that of the host <italic>V. maculifrons </italic>(adapted from [##UREF##15##17##]).</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Locations of 37 <italic>V. maculifrons </italic>and 13 <italic>V. squamosa </italic>colonies collected in this study.</bold> All samples were obtained from the state of Georgia (inset, see also Fig. 1). Lines denote county boundaries within the state. The greatest pairwise distance between colonies exceeds 120 km and colonies of both species were sampled from overlapping regions.</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Relationship between genetic differentiation and geographic distance of workers sampled from <italic>Vespula </italic>colonies.</bold> There was no evidence for genetic isolation by distance in these species over this range.</p></caption></fig> <fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>Distribution of (A) observed mate number and (B) effective mate number (<italic>k</italic><sub><italic>e</italic>3</sub>) for <italic>V. maculifrons </italic>and <italic>V. squamosa </italic>queens.</bold> Within each species the distribution effective mate number is reduced relative to observed mate number because of unequal sperm use by queens. The distributions of observed mate number differed significantly from a Gaussian distribution in both species due to an excess of queens mated to intermediate numbers of males.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Variability metrics of microsatellite loci in <italic>V. maculifrons </italic>(<italic>Vmac</italic>) and <italic>V. squamosa </italic>(<italic>Vsqu</italic>).</p></caption><table frame="hsides" rules="groups"><thead><tr><th/><th align="center" colspan="2"><italic>N</italic><sup>a</sup></th><th align="center" colspan="2"><italic>A</italic><sub><italic>obs</italic></sub><sup>b</sup></th><th align="center" colspan="2"><italic>A</italic><sub>50</sub><sup>c</sup></th><th align="center" colspan="2"><italic>A</italic><sub><italic>eff</italic></sub><sup>d</sup></th><th align="center" colspan="2"><italic>h</italic><sup>e</sup></th><th align="center" colspan="2">Range (bp)<sup>f</sup></th></tr><tr><th align="left">Locus</th><th align="right"><italic>Vmac</italic></th><th align="right"><italic>Vsqu</italic></th><th align="right"><italic>Vmac</italic></th><th align="right"><italic>Vsqu</italic></th><th align="right"><italic>Vmac</italic></th><th align="right"><italic>Vsqu</italic></th><th align="left"><italic>Vmac</italic></th><th align="left"><italic>Vsqu</italic></th><th align="left"><italic>Vmac</italic></th><th align="left"><italic>Vsqu</italic></th><th align="left"><italic>Vmac</italic></th><th align="left"><italic>Vsqu</italic></th></tr></thead><tbody><tr><td align="left">LIST2003</td><td align="right">259</td><td align="right">107</td><td align="right">31</td><td align="right">17</td><td align="right">20.319</td><td align="right">14.865</td><td align="left">17.642</td><td align="left">8.375</td><td align="left">0.943</td><td align="left">0.889</td><td align="left">185–224</td><td align="left">173–204</td></tr><tr><td align="left">LIST2004</td><td align="right">253</td><td align="right">113</td><td align="right">20</td><td align="right">8</td><td align="right">13.073</td><td align="right">6.623</td><td align="left">8.691</td><td align="left">4.155</td><td align="left">0.885</td><td align="left">0.766</td><td align="left">147–181</td><td align="left">119–149</td></tr><tr><td align="left">LIST2013</td><td align="right">256</td><td align="right">107</td><td align="right">20</td><td align="right">18</td><td align="right">11.839</td><td align="right">13.911</td><td align="left">7.241</td><td align="left">6.357</td><td align="left">0.862</td><td align="left">0.850</td><td align="left">184–213</td><td align="left">173–201</td></tr><tr><td align="left">LIST2019</td><td align="right">224</td><td align="right">108</td><td align="right">11</td><td align="right">4</td><td align="right">8.373</td><td align="right">3.711</td><td align="left">3.907</td><td align="left">1.534</td><td align="left">0.744</td><td align="left">0.352</td><td align="left">126–152</td><td align="left">122–132</td></tr><tr><td align="left">LIST2020</td><td align="right">256</td><td align="right">110</td><td align="right">19</td><td align="right">19</td><td align="right">12.276</td><td align="right">15.061</td><td align="left">7.360</td><td align="left">10.168</td><td align="left">0.864</td><td align="left">0.910</td><td align="left">226–268</td><td align="left">322–356</td></tr><tr><td align="left">Rufa-5</td><td align="right">213</td><td align="right">107</td><td align="right">3</td><td align="right">11</td><td align="right">2.223</td><td align="right">9.066</td><td align="left">1.206</td><td align="left">5.011</td><td align="left">0.171</td><td align="left">0.808</td><td align="left">136–140</td><td align="left">153–175</td></tr><tr><td align="left">VMA-3</td><td align="right">219</td><td align="right">77</td><td align="right">16</td><td align="right">14</td><td align="right">12.543</td><td align="right">13.389</td><td align="left">9.556</td><td align="left">10.420</td><td align="left">0.895</td><td align="left">0.916</td><td align="left">260–286</td><td align="left">254–282</td></tr><tr><td align="left">VMA-6</td><td align="right">256</td><td align="right">109</td><td align="right">39</td><td align="right">20</td><td align="right">22.071</td><td align="right">15.329</td><td align="left">17.722</td><td align="left">10.287</td><td align="left">0.944</td><td align="left">0.911</td><td align="left">262–323</td><td align="left">267–295</td></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>Relatedness (± SEM) estimates for <italic>V. maculifrons </italic>and <italic>V. squamosa</italic>.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Relatedness</th><th align="left">V. maculifrons</th><th align="left">V. squamosa</th></tr></thead><tbody><tr><td align="left">Nestmate workers</td><td align="left">0.373 ± 0.009 ***</td><td align="left">0.357 ± 0.019 ***</td></tr><tr><td align="left">Queens to their male mates</td><td align="left">0.014 ± 0.010</td><td align="left">0.016 ± 0.020</td></tr><tr><td align="left">Among males mated to single queens</td><td align="left">0.023 ± 0.014</td><td align="left">0.053 ± 0.024 *</td></tr></tbody></table></table-wrap> <table-wrap id="T3" position="float"><label>Table 3</label><caption><p>Utility of 43 microsatellite loci in <italic>V. maculifrons </italic>and <italic>V. squamosa</italic><sup>a</sup>.</p></caption><table frame="hsides" rules="groups"><tbody><tr><td align="left">Locus</td><td align="left"><italic>V. maculifrons</italic></td><td align="left"><italic>V. squamosa</italic></td></tr><tr><td/><td/><td/></tr><tr><td align="left">LIST2001</td><td align="left">Var</td><td align="left">-</td></tr><tr><td align="left">LIST2002</td><td align="left">Var</td><td align="left">-</td></tr><tr><td align="left">LIST2003</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">LIST2004</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">LIST2006</td><td align="left">Var</td><td align="left">+++</td></tr><tr><td align="left">LIST2007</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">LIST2008</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">LIST2009</td><td align="left">+++</td><td align="left">+++</td></tr><tr><td align="left">LIST2010</td><td align="left">Var</td><td align="left">+++</td></tr><tr><td align="left">LIST2011</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">LIST2012</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">LIST2013</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">LIST2014</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">LIST2015</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">LIST2016</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">LIST2017</td><td align="left">Var</td><td align="left">-</td></tr><tr><td align="left">LIST2018</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">LIST2019</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">LIST2020</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td/><td/><td/></tr><tr><td align="left">Rufa 1</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">Rufa 2</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">Rufa 3</td><td align="left">+++</td><td align="left">Var</td></tr><tr><td align="left">Rufa 4</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">Rufa 5</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">Rufa 6</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">Rufa 7</td><td align="left">+++</td><td align="left">+++</td></tr><tr><td align="left">Rufa 8</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">Rufa 9</td><td align="left">-</td><td align="left">Var</td></tr><tr><td align="left">Rufa 10</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">Rufa 11</td><td align="left">Var</td><td align="left">+++</td></tr><tr><td align="left">Rufa 12</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">Rufa 13</td><td align="left">-</td><td align="left">Var</td></tr><tr><td align="left">Rufa 14</td><td align="left">-</td><td align="left">+++</td></tr><tr><td align="left">Rufa 15</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">Rufa 16</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">Rufa 17</td><td align="left">+++</td><td align="left">Var</td></tr><tr><td align="left">Rufa 18</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">Rufa 19</td><td align="left">Var</td><td align="left">+++</td></tr><tr><td/><td/><td/></tr><tr><td align="left">VMA-3</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">VMA-4</td><td align="left">-</td><td align="left">-</td></tr><tr><td align="left">VMA-6</td><td align="left">Var</td><td align="left">Var</td></tr><tr><td align="left">VMA 7</td><td align="left">+++</td><td align="left">+++</td></tr><tr><td align="left">VMA-8</td><td align="left">Var</td><td align="left">-</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p><sup>a</sup>Number of genes sampled obtained from the inferred genotypes of females and males. Unequal sample sizes among loci within each species result from missing or incomplete data.</p><p><sup>b</sup>Total number of observed alleles.</p><p><sup>c</sup>Number of alleles that would have been observed if only 50 genes had been sampled.</p><p><sup>d</sup>Effective number of alleles.</p><p><sup>e</sup>Expected heterozygosity.</p><p><sup>f</sup>Range in allele size.</p></table-wrap-foot> <table-wrap-foot><p>Estimate differs significantly from zero at the 0.05 (*) or 0.001 (***) level.</p></table-wrap-foot> <table-wrap-foot><p><sup>a</sup>Designations indicate that a particular locus failed to PCR-amplify (-), PCR-amplified but was not variable (+++), or both PCR-amplified and was variable (VAR). Variability was assessed from visualization of alleles via agarose gel electrophoresis.</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-239-1"/> <graphic xlink:href="1471-2148-8-239-2"/> <graphic xlink:href="1471-2148-8-239-3"/> <graphic xlink:href="1471-2148-8-239-4"/>
[{"surname": ["Maynard Smith", "Szathmary"], "given-names": ["J", "E"], "source": ["The major transitions in evolution"], "year": ["1998"], "publisher-name": ["Oxford , Oxford University Press"]}, {"surname": ["Wilson"], "given-names": ["EO"], "source": ["The insect societies"], "year": ["1971"], "publisher-name": ["Cambridge , Harvard University press"]}, {"surname": ["H\u00f6lldobler", "Wilson"], "given-names": ["B", "EO"], "source": ["The Ants"], "year": ["1990"], "publisher-name": ["Cambridge, Massachusetts , The Belknap Press of Harvard University Press"]}, {"surname": ["Ross", "Matthews"], "given-names": ["KG", "RW"], "source": ["The social biology of wasps"], "year": ["1991"], "publisher-name": ["Ithaca , Comstock Publishing Associates"]}, {"surname": ["Bourke", "Franks"], "given-names": ["AFG", "NR"], "source": ["Social Evolution in Ants"], "year": ["1995"], "publisher-name": ["Princeton, New Jersey , Princeton University Press"]}, {"surname": ["Brandt", "Foitzik", "Fischer-Blass", "Heinze"], "given-names": ["M", "S", "B", "J"], "article-title": ["The coevolutionary dynamics of obligate ant social parasite system-between prudence and antagonism"], "source": ["Biol Rev Biol Rev"], "year": ["2005"], "volume": ["80"], "issue": ["2"], "fpage": ["251"], "lpage": ["267"]}, {"surname": ["Schmid-Hempel"], "given-names": ["P"], "source": ["Parasites in social insects"], "year": ["1998"], "publisher-name": ["Princeton, NJ , Princeton University Press"]}, {"surname": ["Cervo"], "given-names": ["R"], "article-title": ["Polistes wasps and their social parasites: an overview"], "source": ["Annales Zoologici Fennici"], "year": ["2006"], "volume": ["43"], "issue": ["5-6"], "fpage": ["531"], "lpage": ["549"]}, {"surname": ["Lenoir", "D'Ettorre", "Errard", "Hefetz"], "given-names": ["A", "P", "C", "A"], "article-title": ["Chemical ecology and social parasitism in ants"], "source": ["Annual Review of Entomology"], "year": ["2001"], "volume": ["46"], "fpage": ["573"], "lpage": ["599"], "pub-id": ["10.1146/annurev.ento.46.1.573"]}, {"surname": ["Carpenter", "Perera"], "given-names": ["JM", "EP"], "article-title": ["Phylogenetic relationships among yellowjackets and the evolution of social parasitism (Hymenoptera : Vespidae, Vespinae)"], "source": ["American Museum Novitates"], "year": ["2006"], "fpage": ["Cover1"], "lpage": ["19"]}, {"surname": ["Taylor"], "given-names": ["LH"], "article-title": ["Observations on social parasitism in the genus Vespula Thompson"], "source": ["Annals of the Entomological Society of America"], "year": ["1939"], "volume": ["32"], "fpage": ["304"], "lpage": ["315"]}, {"surname": ["Matthews"], "given-names": ["RW"], "person-group": ["Jaisson P"], "article-title": ["Social parasitism in yellowjackets"], "source": ["Social insects in the tropics"], "year": ["1982"], "publisher-name": [" Univ. Paris-Nord"], "fpage": ["193"], "lpage": ["202"]}, {"surname": ["MacDonald", "Matthews"], "given-names": ["JF", "RW"], "article-title": ["Nesting biology of the southern yellowjacket, Vespula squamosa (Hymenoptera: Vespidae): Social parasitism and independent founding"], "source": ["Journal of the Kansas entomological society"], "year": ["1984"], "volume": ["57"], "issue": ["1"], "fpage": ["134"], "lpage": ["151"]}, {"surname": ["MacDonald", "Matthews"], "given-names": ["JF", "RW"], "article-title": ["Vespula squamosa: A yellow jacket wasp evolving towards parasitism"], "source": ["Science"], "year": ["1975"], "volume": ["190"], "fpage": ["1003"], "lpage": ["1004"]}, {"surname": ["Hunt", "Cave", "Borjas"], "given-names": ["JH", "RD", "GR"], "article-title": ["First records from Honduras of a yellowjacket wasp, Vespula squamosa (Drury) (Hymenoptera : Vespidae : Vespinae)"], "source": ["Journal of the Kansas Entomological Society"], "year": ["2001"], "volume": ["74"], "issue": ["2"], "fpage": ["118"], "lpage": ["119"]}, {"surname": ["Akre", "Greene", "MacDonald", "Landolt", "Davis"], "given-names": ["RD", "A", "JF", "PJ", "HG"], "source": ["The yellow jackets of America north of Mexico"], "year": ["1980"], "volume": ["Agriculture handbook 552"], "publisher-name": [" U. S. Department of Agriculture"], "fpage": ["102"]}, {"surname": ["Greene"], "given-names": ["A"], "person-group": ["Ross KG, Matthews RW"], "article-title": ["Dolichovespula and Vespula"], "source": ["The social biology of wasps"], "year": ["1991"], "publisher-name": ["Ithaca , Comstock publishing associates"], "fpage": ["263"], "lpage": ["305"]}, {"surname": ["MacDonald", "Matthews"], "given-names": ["JF", "RW"], "article-title": ["Nesting biology of the Eastern yellowjacket, Vespula maculifrons (Hymenoptera: Vespidae)"], "source": ["Journal of the Kansas entomological society"], "year": ["1981"], "volume": ["54"], "issue": ["3"], "fpage": ["433"], "lpage": ["457"]}, {"surname": ["Foitzik", "Heinze"], "given-names": ["S", "J"], "article-title": ["Microgeographic genetic structure and interspecific parasitism in the ant Leptothorax nylanderi"], "source": ["Ecological entomology"], "year": ["2001"], "volume": ["26"], "fpage": ["449"], "lpage": ["456"], "pub-id": ["10.1046/j.1365-2311.2001.00354.x"]}, {"surname": ["Gandon", "Michalakis"], "given-names": ["S", "Y"], "article-title": ["Local adaptation, evolutionary potential and host-parasite coevolution: interactions between migration, mutation, population size and generation time"], "source": ["Journal of Evolutionary Biology"], "year": ["2002"], "volume": ["15"], "issue": ["3"], "fpage": ["451"], "lpage": ["462"], "pub-id": ["10.1046/j.1420-9101.2002.00402.x"]}, {"surname": ["Seeley", "Tarpy"], "given-names": ["TD", "DR"], "article-title": ["Queen promiscuity lowers disease within honeybee colonies"], "source": ["Proceedings of the Royal Society B-Biological Sciences"], "year": ["2007"], "volume": ["274"], "issue": ["1606"], "fpage": ["67"], "lpage": ["72"], "pub-id": ["10.1098/rspb.2006.3702"]}, {"surname": ["Hughes", "Sumner", "Van Borm", "Boomsma"], "given-names": ["WOH", "S", "S", "JJ"], "article-title": ["Worker caste polymorphism has a genetic basis in Acromyrmex leaf-cutting ants"], "source": ["Proceedings of the national academy of sciences USA"], "year": ["2003"], "volume": ["100"], "issue": ["16"], "fpage": ["9394"], "lpage": ["9397"], "pub-id": ["10.1073/pnas.1633701100"]}, {"surname": ["Crozier", "Fjerdingstad"], "given-names": ["RH", "EJ"], "article-title": ["Polyandry in social Hymenoptera-disunity in diversity?"], "source": ["Annals Zoologici Fennici"], "year": ["2001"], "volume": ["38"], "fpage": ["267"], "lpage": ["285"]}, {"surname": ["Bekkevold", "Boomsma"], "given-names": ["D", "JJ"], "article-title": ["Evolutionary transition to a semelparous life history in the socially parasitic ant Acromyrmex insinuator"], "source": ["Journal of evolutionary biology"], "year": ["2000"], "volume": ["13"], "fpage": ["615"], "lpage": ["623"], "pub-id": ["10.1046/j.1420-9101.2000.00201.x"]}, {"surname": ["Gardner", "Schonrogge", "Elmes", "Thomas"], "given-names": ["MG", "K", "GW", "JA"], "article-title": ["Increased genetic diversity as a defence against parasites is undermined by social parasites: Microdon mutabilis hoverflies infesting Formica lemani ant colonies"], "source": ["Proceedings of the Royal Society B-Biological Sciences"], "year": ["2007"], "volume": ["274"], "issue": ["1606"], "fpage": ["103"], "lpage": ["110"], "pub-id": ["10.1098/rspb.2006.3694"]}, {"surname": ["Foster", "Ratnieks"], "given-names": ["KR", "FLW"], "article-title": ["Paternity, reproduction and conflict in vespine wasps: a model system for testing kin selection predictions"], "source": ["Behavioral ecology and sociobiology"], "year": ["2001"], "volume": ["50"], "fpage": ["1"], "lpage": ["8"], "pub-id": ["10.1007/s002650100336"]}, {"surname": ["Goodisman", "Matthews", "Crozier"], "given-names": ["MAD", "RW", "RH"], "article-title": ["Mating and reproduction in the wasp Vespula germanica"], "source": ["Behavioral ecology and sociobiology"], "year": ["2002"], "volume": ["51"], "fpage": ["497"], "lpage": ["502"], "pub-id": ["10.1007/s00265-002-0470-6"]}, {"surname": ["Ross"], "given-names": ["KG"], "article-title": ["Kin selection and the problem of sperm utilization in social insects"], "source": ["Nature"], "year": ["1986"], "volume": ["323"], "issue": ["6091"], "fpage": ["798"], "lpage": ["800"], "pub-id": ["10.1038/323798a0"]}, {"surname": ["Hasegawa", "Takahashi"], "given-names": ["E", "JI"], "article-title": ["Microsatellite loci for genetic research in the hornet Vespa mandarinia and related species"], "source": ["Molecular ecology notes"], "year": ["2002"], "volume": ["2"], "fpage": ["306"], "lpage": ["308"], "pub-id": ["10.1046/j.1471-8286.2002.00224.x"]}, {"surname": ["Daly", "Archer", "Watts", "Speed", "Hughes", "Barker", "Jones", "Odgaard", "Kemp"], "given-names": ["D", "ME", "PC", "MP", "MR", "FS", "J", "K", "SJ"], "article-title": ["Polymorphic microsatellite loci for eusocial wasps (Hymenoptera: Vespidae)"], "source": ["Molecular ecology notes"], "year": ["2002"], "volume": ["2"], "fpage": ["273"], "lpage": ["275"], "pub-id": ["10.1046/j.1471-8286.2002.00220.x"]}, {"surname": ["Petit", "El Mousadik", "Pons"], "given-names": ["RJ", "A", "O"], "article-title": ["Identifying populations for conservation on the basis of genetic markers"], "source": ["Conservation Biology"], "year": ["1998"], "volume": ["12"], "issue": ["4"], "fpage": ["844"], "lpage": ["855"], "pub-id": ["10.1046/j.1523-1739.1998.96489.x"]}, {"surname": ["Raymond", "Rousset"], "given-names": ["M", "F"], "article-title": ["GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism"], "source": ["Journal of heredity"], "year": ["1995"], "volume": ["86"], "fpage": ["248"], "lpage": ["249"]}, {"surname": ["Goodisman"], "given-names": ["MAD"], "article-title": ["A theoretical analysis of variation in multiple mating in social insects"], "source": ["Sociobiology"], "year": ["2007"], "volume": ["49"], "issue": ["3"], "fpage": ["107"], "lpage": ["119"]}, {"surname": ["Nonacs"], "given-names": ["P"], "article-title": ["Measuring the reliability of skew indices: is there one best index?"], "source": ["Animal Behaviour"], "year": ["2003"], "volume": ["65"], "fpage": ["615"], "lpage": ["627"], "pub-id": ["10.1006/anbe.2003.2096"]}, {"surname": ["Queller", "Goodnight"], "given-names": ["DC", "KF"], "article-title": ["Estimating relatedness using genetic markers"], "source": ["Evolution"], "year": ["1989"], "volume": ["43"], "issue": ["2"], "fpage": ["258"], "lpage": ["275"], "pub-id": ["10.2307/2409206"]}, {"surname": ["Foitzik", "DeHeer", "Hunjan", "Herbers"], "given-names": ["S", "CJ", "DN", "JM"], "article-title": ["Coevolution in host-parasite systems: behavioural strategies of slave-making ants and their hosts"], "source": ["Proceedings of the royal society of London B"], "year": ["2001"], "volume": ["268"], "fpage": ["1139"], "lpage": ["1146"], "pub-id": ["10.1098/rspb.2001.1627"]}, {"surname": ["Nei"], "given-names": ["M"], "source": ["Molecular evolutionary genetics"], "year": ["1987"], "publisher-name": ["New York , Columbia University Press"], "fpage": ["512"]}, {"surname": ["Edwards"], "given-names": ["RE"], "source": ["Social wasps: their biology and control"], "year": ["1980"], "publisher-name": ["East Grinstead , Rentokil"]}, {"surname": ["Moller"], "given-names": ["H"], "article-title": ["Lessons for invasion theory from social insects"], "source": ["Biological conservation"], "year": ["1996"], "volume": ["78"], "fpage": ["125"], "lpage": ["142"], "pub-id": ["10.1016/0006-3207(96)00022-5"]}, {"surname": ["Crosland"], "given-names": ["MW"], "article-title": ["The spread of the social wasp, Vespula germanica, in Australia"], "source": ["New Zealand journal of zoology"], "year": ["1991"], "volume": ["18"], "fpage": ["375"], "lpage": ["388"]}, {"surname": ["Thomas"], "given-names": ["CR"], "article-title": ["The European wasp (Vespula germanica Fab.) in New Zealand"], "source": ["New Zealand department of scientific and industrial research information series"], "year": ["1960"], "volume": ["27"], "publisher-name": ["Auckland "], "fpage": ["1"], "lpage": ["74"]}, {"surname": ["Spradbery", "Maywald"], "given-names": ["JP", "GF"], "article-title": ["The distribution of the European or German wasp, Vespula germanica (F.) (Hymenoptera: Vespidae), in Australia: Past, present and future"], "source": ["Australian journal of zoology"], "year": ["1992"], "volume": ["40"], "fpage": ["495"], "lpage": ["510"], "pub-id": ["10.1071/ZO9920495"]}, {"surname": ["Trontti", "Aron", "Sundstr\u00f6m"], "given-names": ["K", "S", "L"], "article-title": ["The genetic population structure of the ant Plagiolepis xene-implications for genetic vulnerability of obligate social parasites"], "source": ["Conservation Genetics"], "year": ["2006"], "volume": ["7"], "issue": ["2"], "fpage": ["241"], "lpage": ["250"], "pub-id": ["10.1007/s10592-005-9003-y"]}, {"surname": ["Ross"], "given-names": ["KG"], "article-title": ["Laboratory studies of the mating biology of the Eastern yellowjacket, Vespula maculifrons (Hymenoptera: Vespidae)"], "source": ["Journal of the Kansas entomological society"], "year": ["1983"], "volume": ["56"], "issue": ["4"], "fpage": ["523"], "lpage": ["537"]}, {"surname": ["Spradbery"], "given-names": ["JP"], "source": ["Wasps: An account of the biology and natural history of solitary and social wasps"], "year": ["1973"], "publisher-name": ["London , Sidgwick & Jackson"], "fpage": ["408"]}, {"surname": ["Bourke", "Franks"], "given-names": ["AFG", "NR"], "article-title": ["Alternative adaptations, sympatric speciation and the evolution of parasitic, inquiline ants"], "source": ["Biological journal of the Linnean Society"], "year": ["1991"], "volume": ["43"], "fpage": ["157"], "lpage": ["178"], "pub-id": ["10.1111/j.1095-8312.1991.tb00591.x"]}, {"surname": ["Boomsma", "Sundstr\u00f6m"], "given-names": ["JJ", "L"], "article-title": ["Patterns of paternity skew in Formica ants"], "source": ["Behavioral ecology and sociobiology"], "year": ["1998"], "volume": ["42"], "fpage": ["85"], "lpage": ["92"], "pub-id": ["10.1007/s002650050415"]}, {"surname": ["Thornhill", "Alcock"], "given-names": ["R", "J"], "source": ["The evolution of insect mating systems"], "year": ["2001"], "publisher-name": ["Cambridge, MA , Harvard University Press"], "fpage": ["547"]}]
BMC Evol Biol. 2008 Aug 20; 8:239