Introduction {#Sec1}
============

The genus *Polistes* (Hymenoptera: Vespidae) comprises more than 200 wasp species representing a primitively eusocial group of insects (Carpenter [@CR2]). Several factors have led to the intensive study of these species using either traditional (morphologic, behavioral and/or ecologic data) or molecular approaches. These are, for example, interesting population structure and genetics, a wide range of diversity in nesting habits and life cycle, and their role in ecosystems (Pardi [@CR40]; Ono [@CR38]; Sayama and Takahashi [@CR51]; Hughes et al. [@CR19]; Miyano et al. [@CR34]; Nagamati Jr. et al. [@CR35]; Tsuchida et al. [@CR59]).

There are a limited number of studies that have been devoted to identification and the analysis of microsatellite loci in *Polistes* (Strassmann et al. [@CR56]; Henshaw [@CR15]; Saigo and Tsuchida [@CR50]; Henshaw et al. [@CR16]; Uddin and Tsuchida [@CR61]; Komatsu et al. [@CR23]). Here, for the first time, we performed microsatellite analyses of genetic diversity and structure of the population of Eurasiatic *Polistes nimpha* (Christ 1791) (Polistinae).

*Polistes nimpha* is one of the three species of the genus that can be found in Poland (Oleksa and Wiśniowski [@CR37]). Small nests of this species are located on dead stalks of plants predominantly from the Apiaceae and Hypericaceae families, less often on small trees or on blades of grass (up to 50 cm aboveground surface).

The life cycle of *P. nimpha* is similar to that of the other representatives of the genus *Polistes* (see West-Eberhard [@CR68]; Wilson [@CR69]). Spring foundresses of *Polistes* wasps can make different choices regarding reproductive options (e.g., Seppä et al. [@CR53]). Besides haplometrotic nest foundation (nest is founded by a single foundress), pleometrotic nest foundation (nest is started by several queens) is a widespread phenomenon among many species of *Polistes* (Makino and Sayama [@CR28]). The results obtained so far indicate, that in the majority of cases *Polistes* colonies are founded by a single foundress (Miyano and Hasegawa [@CR33]; Sayama and Takahashi [@CR51]), but in *P. nimpha* pleometrotic nest foundation in sheltered places such as buildings was also noted (Cervo and Turillazzi [@CR4]). Moreover, nest usurpation (queen replacement), either intra- or interspecific, in this species was observed (Cervo et al. [@CR5]; Cervo [@CR3]; Lorenzi et al. [@CR25]). Although the other behavioral choices of foundresses during the pre-emergency period were described in *P. dominula* Christ 1791 (Nonacs and Revee 1995) and *P. carolina* (L. 1767) (Seppä et al. [@CR53]), there is no evidence of visiting, deserting, switching, moving, adopting, or nest joining in *P. nimpha*. There is also no evidence of how reproduction is partitioned in the nest between queens and subordinate foundresses and/or potential joiners, movers, or adopters. Whether queens of *P. nimpha* mate with a single or multiple males also has not been tested (but see Turillazzi and Cervo [@CR60]).

This study represents the first attempt to understand genetic diversity and population structure and relatedness/relationship between individuals of *Polistes nimpha.* The purpose of this research is twofold: (1) to determine the amplification efficiency and polymorphism of *P. dominula* microsatellites (Henshaw [@CR15]) in *P. nimpha* and (2) to explore genetic evidence for possible reproductive options in this species.

Materials and methods {#Sec2}
=====================

Study area {#Sec3}
----------

The study area was located at the edges of a military training field near the city of Poznań, Poland (Fig. [1](#Fig1){ref-type="fig"}). Three area plots (17,400; 5300; and 14,200 m^2^) where nest sites occur were selected in April 2014.Fig. 1Location of *Polistes nimpha* nests (1--15) in Poland

Nest searching and monitoring {#Sec4}
-----------------------------

Established plots were searched thoroughly for nests of *P. nimpha*. The initial nest searching began on May 12, 2014, and each recorded nest was revisited at approximately 7--10 (21)-day intervals. Nests were watched at forenoon hours 8:00 a.m. to 12:00 noon during the active hours of the inhabitants. One week after the nest identification, from May 19 to May 23, 2014, foundresses of *P. nimpha* were marked individually. Each nest was observed for 45 min per week during the following 5 weeks until days after the first workers' appearances. Nest contents, a number of cells, and causes of nest failure were documented. Nests were considered as inactive if foundress was not present on the nest for two consecutive following checks. All nests have been marked with numbered tags. Moreover, the nest locations were measured at their latitude and longitude coordinates.

Sample collection {#Sec5}
-----------------

For genetic analysis, fifty-nine individuals were collected on June 9, 2014 from 15 nests (two pupae and two larvae per nest), except one nest from which one pupa and two larvae were taken. The pupae and larvae were preserved in 95 % alcohol and stored at low temperature (4 °C). Vouchers are deposited in the Department of Systematic Zoology, Faculty of Biology, Adam Mickiewicz University in Poznań, Poland.

DNA extraction and microsatellite amplification {#Sec6}
-----------------------------------------------

Genomic DNA was extracted from the legs of each pupa (right leg of the second pair) or from ca. $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ 1/3{-}1/2 $$\end{document}$ of the larvae body using DNeasy Blood and Tissue Kit (Qiagen GmbH, Hilden, Germany) according to manufacturer's instructions.

The PCR primers used were those developed for *P. dominula* (Henshaw [@CR15]). PCRs were carried out in 5 μl reaction volumes containing 1× Type-it Microsatellite Kit (Qiagen), 0.07 μM forward primer, 0.2 μM reverse primer, 0.14 μM fluorescent-labeled M13 primer, and 1 μl (6--65 ng) of DNA template using a thermocycling profile of one cycle of 5 min at 96 °C followed by 30 steps of 30 s at 95 °C, 90 s at 50 °C, and 30 s at 72 °C, with a final step of 30 min at 65 °C.

Amplicons were fluorescently labeled in multiplexed reactions using a modified M13-tailing method (Oetting et al. [@CR36]). The amplified alleles were separated on an ABI PRISM 3130XL (Applied Biosystems) with Genescan 600LIZ size standard and scored with the program Peak Scanner v 1.0 (Applied Biosystems).

Genetic diversity analysis {#Sec7}
--------------------------

Genetic analyses were completed on a pooled sample of 59 individuals (pupae and larvae) collected from 15 nests. To describe informativeness of genetic markers related to expected heterozygosity, the polymorphism information content (PIC) values were calculated using the program Cervus. The program Micro-Checker v. 2.2.3 (Van Oosterhout et al. [@CR62]) was used to test for the presence of null alleles, stuttering during the PCR amplification, and large allele dropout.

Genotype frequencies at any locus are a function of allele frequencies in the absence of migration, mutation, assortative mating, and natural selection. This situation is termed Hardy--Weinberg equilibrium. Deviations from this equilibrium can imply inbreeding, physical mixing of populations (Wahlund effect), population stratification, and/or null alleles. The deviation from Hardy--Weinberg proportions caused by heterozygote excess was tested with a global test (*U* test), as implemented in the program GenPop v. 4.2 (Raymond and Rousset [@CR45]; Rousset [@CR48], [http://kimura.univ-montp2.fr/\~rousset/Genepop.htm](http://kimura.univ-montp2.fr/%7erousset/Genepop.htm)), using 1000 dememorization steps, 100 batches, and 1000 iterations per batch. To assess the Hardy--Weinberg equilibrium for each locus, the goodness-of-fit Chi-square tests (Weir [@CR66]) with Bonferroni correction, as implemented in the program Cervus v. 3.0.6., were performed.

To assign individuals to groups based on their multilocus genotypes, a Bayesian clustering algorithm was used, as implemented in the BAPS software version 6.0 at <http://www.helsinki.fi/bsg/software/BAPS> (Corander et al. [@CR6]). The data were analyzed using a genetic admixture model for unlinked markers with the a priori upper bound for the number of clusters equal to 15. A principal coordinate analysis (PCoA), available in the program GeneAlEx v. 6.501 (Peakall and Smouse [@CR41], [@CR42], <http://www.anu.edu.auG/BoZo/enAlEx>) was conducted using Nei's unbiased distance pairwise matrix to test whether patterns in microsatellite data support the partitioning of the sample into BAPS clusters.

As suggested by Frantz et al. ([@CR11]), it is important to test the analyzed data for isolation-by-distance before applying the Bayesian method of clustering. The relationship of genetic similarity and geographic distance is expected to be linear under the isolation-by-distance hypothesis (Rousset [@CR47]). To check this hypothesis, a paired Mantel test (Mantel [@CR30]) was performed for analyzed individuals with the program GeneAlEx. The test was based on linear genetic distances and geographic distances calculated on UTM coordinates of Lat./Log.; regression significance was tested with 999 permutations.

Kinship analysis, which relies on a model that assigns dyads to discrete relationship categories (Weir et al. [@CR67]), was performed, based on three estimators: Wang's ([@CR63]), Queller and Goodnight's ([@CR44]) and the maximum-likelihood estimator (Kalinowski and Taper [@CR21]). These estimators are implemented in the programs Coancestry v. 1.0.1.5 (Wang [@CR64]) and ML-Relate (Kalinowski and Taper [@CR21]). A correlation coefficient (*r*) was calculated to measure associations between the relatedness estimators using Statistica v. 10.0 (Statsoft). A *t* test was performed with the same program to determine whether differences between means of relatedness are significant between distinguished groups.

Results {#Sec8}
=======

Systematic survey of nests {#Sec9}
--------------------------

One hundred fifty-seven nests of *Polistes nimpha* were initially found on the selected area plots. A repeat systematic survey of nests has shown that the number of nests recorded decreased substantially by 85.30 % in a short time (the first 3 weeks) (Fig. [2](#Fig2){ref-type="fig"}). At the time when pupae and larvae were collected (June 9, 2014), only 17 nests had survived (Fig. [2](#Fig2){ref-type="fig"}). Materials for genetic analyses were collected from 15 nests because in the two remaining nests the number of individuals was low. More intense sampling during this phase might diminish substantially or destroy the local population.Fig. 2Number of nests of *Polistes nimpha* on the study area (Poland) recorded from May to August 2014

The nest destruction was not witnessed, but trace evidence examinations allowed us to conclude that activities of large animals (e.g., foxes), invertebrates (e.g., ants) and/or humans had destroyed the nests. In over 80 % of the cases, the mammals' activity was the most important predictor of nest mortality. For a large proportion of nests, the reason for their failure is unknown. The degree of destruction varied substantially across sites; only stalks of nests remained (the comb had been damaged), or inhabitants had been eaten, but the nest (comb) had remained untouched (Table [1](#Tab1){ref-type="table"}).Table 1Percentage of destroyed and abandoned nests of *Polistes nimpha* and external sources of nests destruction based on observations made at three plots in the Wielkopolska region (Poland) in 2014Causes of mortalityDate of observationMay 12Jun 5Jun 11Jun 18Jun 25Jul 2Jul 9Jul 18Jul 28Aug 6Aug 13Aug 20Mammals' activity070.700.64000001.9100.640Ants' activity01.910000000000Anthropogenic activity000000.64000000Abandoned nests01.9100000000.6400Undetected causes of nest failure010.833.182.550.641.91000000Total  % of nestsActive10014.6510.838.287.645.095.095.093.181.270.640.64Destroyed083.443.822.550.642.55001.911.910.640

Due to the considerable nest mortality (Fig. [2](#Fig2){ref-type="fig"}), only 42 foundresses (of 157 initially identified nests) were marked. Eleven marked foundresses originated from nests sampled for genetic analyses. Based on the mobility behavior of marked foundresses, we can conclude that all these individuals have been observed only on their natal nests. However, not marked, visiting, or usurping (not aggressively behaved) wasps were observed during weekly inspections of localized colonies of *P. nimpha* on May 8, 2014 and May 14, 2014.

Microsatellite cross-species amplification {#Sec10}
------------------------------------------

High rates of microsatellite amplification were observed. From a total of 13 primer pairs, 12 pairs that amplify without any problems were chosen for genetic analysis. Table S1 shows the characteristics of the pursued and nonamplifiable pair of primers. Eleven microsatellite loci (*Pdom 1, 2, 7, 20, 25, 93, 117, 121, 127, 139, 140*) were polymorphic, and one locus (*Pdom 151*) was monomorphic (Table S1). Only eight polymorphic loci had a size variation consistent with the occurrence of strictly stepwise mutations in the repeat array. The remaining primers (*Pdom 25, 117, 121*) showed more complex patterns of evolution (nonstepwise mutations). These three loci were excluded; thus the truncated set contains data for eight microsatellite loci collected from 59 individuals. Large allele dropout and stuttering were not detected. At the locus *Pdom 7*, missing data for four individuals probably appeared because of mutations in the primer binding sites.

Genetic diversity and structure within the pooled sample {#Sec11}
--------------------------------------------------------

Standard measures of genetic diversity in *P. nimpha* population were calculated. The mean PIC score of the microsatellite loci (*n* = 8) was equal to 0.545, ranging from 0.033 (*Pdom 93*) to above 0.5 (*Pdom 1, 139, 7, 127*). The latter four loci are highly polymorphic. The mean number of different alleles per locus was 6.875, ranging from 2 (*Pdom 93*) to 11 (*Pdom 7*). The mean expected heterozygosity was 0.5765 (unbiased expected heterozygosity UHe = 0.577), meaning that there is a \>57 % chance of being a heterozygote under random mating conditions. The average value of the observed heterozygosity (*H*~o~) was slightly higher than Hardy--Weinberg expectations, *H*~o~ = 0.588. However, no significant deviations from Hardy--Weinberg equilibrium were observed for the alternative hypothesis of heterozygote excess (H1), (*P* = 0.433 ± SE 0.020). Chi-square tests for individual loci showed a significant deviation (*P* \< 0.05) for *Pdom 20* (Chi-square test value using Yates's correction = 7.622), (Table S2).

To relate observable structure (nests) to the genetic structure of analyzed population, a model-based clustering was employed. Taking genetic information into account, the Bayesian clustering method, based on the individual-level clustering model for unlinked markers, gave the probability of \>0.898 of being 10 genetic clusters within the pooled sample of 15 nests (log marginal likelihood of optimal partition = --1025.9824). The inferred six clusters (4, 6--10) were coherent in terms of wasps' nests (i.e., each cluster contains individuals originated from one nest). The remaining clusters (1--3, 5) were composed of all individuals originated from two nests (clusters 1, 2, and 5) or three nests (cluster 3). The principal coordinate analysis (PCoA) results broadly confirmed the grouping based on BAPS analysis. In the PCoA, the individuals were grouped accordingly, but only 40.97 % of the variation in the input data has been captured by the first three axes (PCoA 1--3), (Figs. [3](#Fig3){ref-type="fig"} and [4](#Fig4){ref-type="fig"}). All individuals studied belong to the same generation and based on genetic data all were heterozygous at least at one of analyzed loci. Thus, they represent presumably the same sex (worker-like females), (see Balloux and Lugon-Moulin [@CR1], Lorenzi [@CR24]).Fig. 3Principal coordinate plot of *Polistes nimpha* individuals based on Nei's standard distances. Axis 1 explains 17.54 % of the variation, and axis 2 explains 12.78 % of the variation. Clusters of individuals (clusters 1--10 marked with *different colors*) were obtained in a BAPS analysis (Bayesian analysis of population structure) (color figure online)Fig. 4Principal coordinate plot of *Polistes nimpha* individuals based on Nei's standard distances. Axis 1 explains 17.54 % of the variation, and axis 3 explains 10.65 % of the variation. Clusters of individuals (clusters 1--10 marked with *different* *colors*) were obtained in a BAPS analysis (Bayesian analysis of population structure) (color figure online)

In the case of organism's limited dispersal ability, we expect an isolation-by-distance type of population structure, in which genetic similarity is contrarily related to geographical distance calculations (Meirmans and Van Tienderen, [@CR31]). The comparison of linear genetic distances with geographic distances by the Mantel test revealed a significant correlation between genetic and geographic distances (*R*^2^ = 0.0413, *P* \< 0.001); thus, it can be concluded that genetic differences among individuals do increase linearly with geographic distance.

The three metrics of relatedness employed were highly correlated (*n* = 1711, all *r* \> 0.91). Based on the method-of-moments estimators, the BAPS clusters composed of individuals originated from two or three nests (considered as a group 1, nos. 1--3, 5) expressed significantly lower levels of relatedness coefficients (*P* \< 0.05) than did the second group composed of one-nest clusters (nos. 4, 6--10) (Table S3; Figs. [5](#Fig5){ref-type="fig"} and [6](#Fig6){ref-type="fig"}). The ML estimator did not show significant differences between these groups.Fig. 5Mean and 95 % confidence interval for relatedness coefficients were based on the method-of-moments approach (Queller and Goodnight [@CR44]; Wang [@CR63]) and maximum-likelihood approach (ML estimator) (Kalinowski and Taper [@CR21]) between all possible dyads of *Polistes nimpha* within clusters of individuals obtained in a BAPS analysis (Bayesian analysis of population structure)Fig. 6Percentage of pairwise comparisons between individuals of *Polistes nimpha* divided into relationship categories (full-sibs, half-sibs and unrelated) based on the maximum-likelihood estimator of Kalinowski and Taper ([@CR21]) within clusters of individuals obtained in a BAPS analysis (Bayesian analysis of population structure), and between clusters (these denoted as "11"). Legend: full-siblings (*black circles*), half-siblings (*triangles*), and unrelated individuals (*white circles*)

Discriminations among plausible relationships (full-sibs, half-sibs, and unrelated) were based on the maximal likelihood. All pairwise comparisons between individuals classified to mixed cluster no. 1 (composed of two nests, nos. 2 and 9, Fig. [1](#Fig1){ref-type="fig"}) were classified as full-siblings or parental offspring (here, both categories were considered jointly as full-siblings). However, the mean value of relatedness within this cluster was significantly lower than that within the one-nest clusters (Table S4). In the mixed clusters (nos. 2--3 and 5), the percentage of pairwise comparisons classified as full-siblings ranged from 46.15 to 57.14 %; as half-siblings, from 23.81 to 35.71 %; and as unrelated individuals, from 7.14 to 24.615 % (Fig. [6](#Fig6){ref-type="fig"}). In one-nest clusters (nos. 6--10), only full-siblings were found. In one-nest cluster no. 4, the majority of pairwise comparisons were classified as full-sibs, whereas one-sixth were classified as half-siblings; no unrelated individuals were found.

Discussion {#Sec12}
==========

Eight loci out of 12 analyzed can be effectively used to genotype individuals of *P. nimpha.* These most variable perfect loci have the average expected and observed heterozygosities over 0.57, that are comparable but slightly lower than those detected for *P. dominula* (Henshaw [@CR15]; Henshaw et al. [@CR16]), *P. chinensis antennalis* (Saigo and Tsuchida [@CR50]), and *P. satan* (Komatsu et al. [@CR23]).

In the cases of limited number of microsatellite loci (5--20), the method-of-moments estimators is commonly preferred (Csilléry et al. [@CR7]). Nevertheless, different relatedness estimators employed here yielded highly correlated estimates of relatedness between pairs of *P. nimpha* individuals. Therefore, besides the commonly used, the Queller and Goodnight's ([@CR44]) estimator that was developed for haploid--diploid social Hymenoptera, the maximum-likelihood (ML) estimators (Milligan [@CR32], Kalinowski and Taper [@CR21]) can be successfully employed in the further studies of social insects.

The Bayesian analysis provided support for the existence of genetic structure in the analyzed population of *P. nimpha* by partitioning individuals into 10 clusters. This grouping was also broadly confirmed by principal coordinate analyses. The existence of mixed BAPS clusters (i.e., composed of two or three nests) in which only some individuals that originated from different nests are full-sibs and the others are unrelated may indicate three main reproductive tactics in the population of this species: short visits of foundresses, cooperative (pleometrotic) nest foundations, and/or nest usurpation. These tactics in *Polistes* wasps are hypothesized to have different fitness payoffs over time, depending on environmental shifts (Starks and Fefferman [@CR54]; Seppä et al. [@CR53]).

Pleometrotic nest foundation was described for *Polistes* species, e.g., European *P. dominula* (Pardi [@CR39]; Nonacs and Reeve 1995), American *P.* c*arolina* (L.) (Seppä et al. [@CR53]), neotropical *P. versicolor* (Olivier) (De Oliveira et al. [@CR8] and references therein), and *P. nimpha* (Turillazzi and Cervo 1985). As was hypothesized by Gamboa ([@CR12]) and Tibbetts and Reeve ([@CR58]), multiple foundresses increase colonies' productivity and survivorship and make them less likely to be usurped. Nevertheless, this tactic may negatively affect the nest if foundresses are infected by parasites (Hodges et al. [@CR18] and references therein). To some extent, pleometrotic nest foundation is contradicted by our evidence of marked foundresses of *P. nimpha* that were seen only on natal (single) nests. Nevertheless, this kind of cooperation cannot be fully excluded due to methodological reasons (marking of wasps was performed shortly after eggs were laid but before offspring emerged). So far, in *P. nimpha*, pleometrotic nest foundation was detected only at sites where the aggregation of sibling females was possible during the winter (i.e., "in covered and sheltered" places, such as buildings) (Turillazzi and Cervo 1985: 49). The winter destruction of plants and difficulties in the aggregation of females that build nests on this kind of substratum were hypothesized to be main factors for nest foundation by a single foundress (Turillazzi and Cervo 1985). Based on materials, we collected from the nests built on vegetation; the single-foundress colonies comprise only 50 % of detected BAPS clusters. Thus, the causes of pleometrotic nest foundation have yet to be investigated and determined.

Nevertheless, other possibilities also should be considered to interpret population structure of *P. nimpha* from our genetic data. The possibilities of laying eggs by visiting, nonaggressive wasps in the nests of sibling females or intraspecific usurpation cannot be excluded based on our direct field observations. It should be underlined, however, that in the case of unmarked wasps, there are always difficulties in differentiating a visitor, a joiner, a switcher, an adopter, or a usurper (Nonacs and Reeve 1985). The case of cluster no. 1 being composed of two nests in which full-siblings were found is a prime example of ambiguity because the interpretation could be twofold: eggs were laid by one queen who was fertilized by related males, or eggs were laid by a pair of full sisters. The former phenomenon was discovered for the first time by Seppä et al. ([@CR52]) in *P. biglumis* L.

The results of the comparison of the genetic distances to the geographic distances with the Mantel test suggest that individuals from a particular nest are more genetically related to those inhabiting nearby nests than they are to individuals originating from nests situated in further parts of the study area. Foundresses of the studied nests may be the offspring of one or several closely related queens that returned to their natal sites after hibernation. Cervo and Turillazzi ([@CR4]) observed that new colonies of *P. nimpha* in synanthropic conditions were initiated very closely to the natal nest. It has been shown that the recognition of the natal nest is possible due to the hydrocarbon profile specifics of each *Polistes* colony (Gamboa [@CR13]; Sumana et al. [@CR57]). Philopatric tendencies in some species of *Polistes* wasps were observed previously by West-Eberhard ([@CR68]), Klahn ([@CR22]), Hirose and Yamasaki ([@CR17]), Makino et al. ([@CR29]) and Cervo and Turillazzi ([@CR4]). For example, most colonies of Japanese *Polistes* wasps are initiated within a distance of 50--80 m (164.5--263 feet) from the natal nest (Hirose and Yamasaki [@CR17]; Makino et al. [@CR29]). Nevertheless, no isolation by distance among populations of *P. olivaceus* (De Geer) was detected (Uddin and Tsuchida [@CR61]).

However, the tendency of individuals originating from two or three different nests to cluster together can reflect close affinity of some individuals or be a result of stochastic error (Kalinowski [@CR20]). The number of clusters detected by Bayesian clustering algorithms does not always correspond to the number of biologically meaningful populations (groups) in the studied sample (Francois and Durand [@CR10]). Inference of population structure is sensitive to the number of loci scored, the number of populations, the number of individuals per population, the sampling strategy, and the chosen analysis method (e.g., Evanno et al. [@CR9]; Rodriguez-Ramilo et al. [@CR46]; Francois and Durand [@CR10]; Kalinowski [@CR20]; Putman and Carbone [@CR43]). Moreover, as was shown by Frantz et al. ([@CR11]), the Bayesian methods employed in BAPS overestimate genetic structure when one or more loci significantly deviate from the Hardy--Weinberg equilibrium and when isolation by distance exhibits high levels (*b* ≤ −0.01).

Similar to the majority of wildlife studies, this study suffers from several limitations resulting from the small number of available loci for analysis, the relatively small sample size, and the small size of the area from which studied individuals were collected. Moreover, based on data sets analyzed here, isolation-by-distance and a significant deviation from Hardy--Weinberg proportions for one locus were detected. Thus, the presented results should be treated with caution until more data are collected.

A large number of colonies established by paper wasps are relatively soon completely destroyed or substantially damaged by natural enemies such as predators and parasites and by human activity of various kinds (Strassmann [@CR55]; Makino [@CR26]; Cervo and Turillazzi [@CR4]; Rusina [@CR49]) The highest rates of nest mortality, attributable to an external agent, were observed in *P. nimpha* during the first weeks after nest identification. After nest destruction, foundresses can choose the following options: nesting alone (nest rebuilding), usurping established colonies, or joining other queens (Nonacs and Reeve 1985). Based on our data, no nests rebuilding or joining has been observed (see also Makino [@CR27]). Therefore, dispersing or dying of foundresses has to be considered to explain this observation.

In conclusion, our study demonstrates for the first time a useful panel of microsatellite markers to perform population genetics studies dealing with *P. nimpha.* There are three main possibilities that come into play to explain our genetic results and direct field observations: cooperative nest foundation, visitation, and/or usurpation events. So far, there is no conclusive evidence to exclude or support any of these possibilities. A philopatric behavior of foundresses is suggested based on genetic data. However, further and more extensive studies in terms of the number of individuals and loci are needed.

Electronic supplementary material
=================================

 {#Sec13}

###### 

Supplementary material 1 (PDF 122 kb)

###### 

Supplementary material 2 (PDF 195 kb)

###### 

Supplementary material 3 (PDF 111 kb)

###### 

Supplementary material 4 (PDF 191 kb)

The research was supported by the Dean of Faculty of Biology at Adam Mickiewicz University in Poznań, Poland, Grant No. GDWB-06/2013 to Krzysztof B. Kozyra.
