- Research article
- Open Access
Genetic clusters and sex-biased gene flow in a unicolonial Formica ant
BMC Evolutionary Biologyvolume 9, Article number: 69 (2009)
Animal societies are diverse, ranging from small family-based groups to extraordinarily large social networks in which many unrelated individuals interact. At the extreme of this continuum, some ant species form unicolonial populations in which workers and queens can move among multiple interconnected nests without eliciting aggression. Although unicoloniality has been mostly studied in invasive ants, it also occurs in some native non-invasive species. Unicoloniality is commonly associated with very high queen number, which may result in levels of relatedness among nestmates being so low as to raise the question of the maintenance of altruism by kin selection in such systems. However, the actual relatedness among cooperating individuals critically depends on effective dispersal and the ensuing pattern of genetic structuring. In order to better understand the evolution of unicoloniality in native non-invasive ants, we investigated the fine-scale population genetic structure and gene flow in three unicolonial populations of the wood ant F. paralugubris.
The analysis of geo-referenced microsatellite genotypes and mitochondrial haplotypes revealed the presence of cryptic clusters of genetically-differentiated nests in the three populations of F. paralugubris. Because of this spatial genetic heterogeneity, members of the same clusters were moderately but significantly related. The comparison of nuclear (microsatellite) and mitochondrial differentiation indicated that effective gene flow was male-biased in all populations.
The three unicolonial populations exhibited male-biased and mostly local gene flow. The high number of queens per nest, exchanges among neighbouring nests and restricted long-distance gene flow resulted in large clusters of genetically similar nests. The positive relatedness among clustermates suggests that kin selection may still contribute to the maintenance of altruism in unicolonial populations if competition occurs among clusters.
Many social animals live in small family groups of closely related individuals, and the role of kinship in promoting the evolution of reproductive altruism and eusociality by kin selection has long been recognized [1–4]. However, some ant species have an extraordinary social organization, called 'unicoloniality', which is characterised by an absence of behavioural boundary among interconnected nests that contain many queens and exchange workers, brood and fertile queens .
The evolution and maintenance of unicolonial populations constitutes one of the enduring major challenges for kin selection, and more generally for modern, gene-centred evolutionary theory [6–9]. This is because high queen number and mixing of individuals among nests can lead to extremely low relatedness among nestmates [6, 8, 10, 11]. In short, if workers help completely unrelated individuals, selection should favour selfish females that develop into queens or produce males, and adaptive altruistic worker behaviour cannot be maintained by kin selection . Hence, the long-term maintenance of altruism critically depends on the actual degree of relatedness among nestmates, which in turn depends on the number of breeders, dispersal of queens and males, and scale at which competition takes place [12–14]. It is therefore of considerable interest to study the breeding system, dispersal pattern and genetic structuring in unicolonial populations.
Studies of unicolonial ants have mostly focussed on invasive species such as the argentine ant Linepithema humile [9, 15–17], the fire ant Solenopsis invicta , or the little fire ant Wasmannia auropunctata [19, 20]. In their invasive range, some of these unicolonial species form very large supercolonies [9, 15, 19]. Supercolonies constitute closed breeding units that are large enough to prevent direct interactions between individuals from distant nests within the same supercolony . Members of supercolonies are aggressive towards members of other supercolonies, but not towards members of their supercolony. In their native range, W. auropunctata seems to be multi-colonial , whereas L. humile forms supercolonies of much smaller size than in the introduced range . If these small supercolonies are closed breeding units that compete which each other, explaining the evolution and stability of unicoloniality becomes less problematic .
Several native, non-invasive ant species have been reported to form unicolonial populations in the genera Formica [21–25], Lasius [7, 26], Myrmica [27–30], Tetramorium  and Polyrhachis . Interestingly, most of these populations of native species do not seem to have clear-cut behavioural boundaries (i.e. aggression) delineating supercolonies, even over large geographical areas in some cases [e.g. [7, 33]]. Moreover, several of the non-invasive unicolonial ant species show significant isolation by distance within populations [i.e. "population viscosity", [1, 22, 29, 32]].
Hamilton  first proposed that population viscosity due to limited dispersal of individuals from their birth place increases the relatedness among interacting individuals, and may thus favour the evolution of altruistic behaviour. However, formal models show that this effect critically depends on the definition of altruism and on details of the life cycle, such as if altruism is provided by workers or focal queens, and occurs before or after dispersal [34–36]. Altruism is more likely to evolve when workers provide it, when there are few breeders, and when the beneficiaries of altruism disperse (e.g. winged queens and males). When altruists directly compete with the individuals they help, the benefits of helping related individuals is reduced by increased competition among relatives [13, 36, 37]. However, if the number of queens is moderate and migration restricted, population viscosity can promote the evolution of altruistic sterile workers under a simple life-cycle in which workers help queens after dispersal .
These models show that the impact of population genetic structure on the evolution of altruism critically depends on specificities of the life-cycle, and in particular on the patterns of dispersal and behavioural interactions of queens, males and workers. It is therefore important to obtain detailed empirical data on sex-specific dispersal and on cryptic genetic discontinuities that may correspond to cooperative units in unicolonial species. Such data will be useful to develop new models and design behavioural studies that will help to assess the role of kin selection in structured population.
Here, we focus on the genetic structure of Formica paralugubris, a wood ant native to the European alpine area. Each nest contains hundreds of reproductive queens [22, 38, 39]. Winged queens and males are produced once per year in spring, and a mating flight occurs in early summer. Formica workers have retained ovaries and are able to lay haploid eggs, but the proportion of egg-laying workers varies across species and worker-laid brood can be policed by other workers . It would be of great interest to study whether F. paralugubris workers gain direct fitness by producing males. However, such studies are hampered by the difficulty to detect worker-produced males with genetic markers when there are many queens per nest.
In F. paralugubris, young queens have alternative reproductive strategies . They can mate and stay within their natal nest [41, 42], or fly to open meadows where they mate before seeking adoption in established nests of F. paralugubris or Serviformica species . Later on, some adult dealate queens may also disperse on foot to join neighbouring nests or establish new nests with the help of workers [colony budding, [11, 22]].
The dispersal behaviour of males has been little studied so far. They may mate within their natal nest or fly away to mate on open meadows, or possibly in other nests . In addition to direct male dispersal, any flight by mated queens will also contribute to increase male gene flow, as sperm is carried along in the spermatheca of queens. A comparison of biparentally-inherited nuclear genetic markers and maternally-inherited mitochondrial markers would permit to evaluate the relative contribution of each sex to gene flow in this ant species with complex and potentially sex-biased dispersal [e.g. [43–45]]. The estimates of numbers of migrants should however be considered with caution and only compared in relative terms, particularly when the markers differ in their levels of genetic variation [46, 47] and the system departs from the assumptions of an island model .
F. paralugubris forms unicolonial populations, as demonstrated by the very high number of queens per nest, movement of queens and workers among neighbouring nests, and lack of aggression between workers from distant nests within populations [22, 33, 38, 39, 49]. Microsatellite analyses in one population revealed that distant nests were genetically differentiated, resulting in low but significant relatedness among nestmate workers [22, 42]. This genetic differentiation indicates that long-distance gene flow between established nests is limited and that most young queens mate and stay within or close to their natal nest . Importantly, such restricted dispersal of queens and males may produce clusters of related nests.
The development of Bayesian clustering methods for population genetics now enables the detection of cryptic genetic population structure using the genotypes of individuals as the sole source of information [50, 51]. In short, these methods attempt to partition individuals into groups at Hardy-Weinberg equilibrium. In a recently developed extension of these methods, geo-referenced individuals are partitioned into sets of panmictic populations, which allows to detect the spatial location of cryptic genetic discontinuities . These new tools are especially useful in species like unicolonial ants where no apparent landscape and behavioural borders are present among nests.
The main aim of our study was to investigate and compare the fine-scale population genetic structure of three unicolonial populations of F. paralugubris. We therefore developed mitochondrial markers that, together with the nuclear microsatellite markers available, allowed us to detect cryptic genetic clusters of nests, obtain more precise information on genetic differentiation within populations, and separate male from female gene flow. These new data should help to evaluate the role of kin selection in native unicolonial populations.
We screened approximately 5165 bp of the mitochondrial DNA of F. paralugubris for sequence polymorphisms in each of 11 individuals coming from 11 nests spread over seven populations and two mountain ranges. Despite the large distance of up to 76 km among sampled individuals, we detected only three polymorphic sites that, in combination, resulted in four haplotypes (H1 to H4). All four haplotypes were present in the population at La Dôle. Three haplotypes (H1, H3, H4) were found at Champs Simon, and two haplotypes (H3, H4) at Bois de Peney and Chalet à Roch. We found only one haplotype (H3) in the population at Bois de Ban, and another single haplotype (H1) in the two alpine populations (Château d'Oex and Lac de l'Hongrin). We detected two haplotypes in many of the nests (40%, 67% and 85% at Bois de Peney, La Dôle and Chalet à Roch, respectively; Fig. 1 panels 2 and 3), which indicates that nests occasionally accept foreign queens.
Genetic clustering of nests within populations
The analyses of the microsatellite genotypes with the program Geneland revealed the presence of cryptic genetic clusters of nests (Fig. 1 panels 1). The number of inferred clusters was the same in all ten runs performed for each population, with posterior probabilities ranging from 53 to 72%. The number and composition of clusters was identical when analysing the data with the program Structure, which confirms that the clustering obtained with Geneland was robust and did not depend on the spatial position of nests. The analysis with Structure also showed some admixture among clusters, as several individuals had a large fraction of their genotype assigned to a cluster other than the one in which they were sampled (Fig. 1 panels 3).
The amount of genetic differentiation among nest clusters varied among populations, but was significant in all cases (Table 1). The limits of the nest clusters based on microsatellite genotypes largely coincided with discontinuities in the frequency distribution of the mitochondrial haplotypes (Fig. 1 panels 2 and 3). As a result, each population exhibited significant mitochondrial DNA haplotype frequency differentiation among the microsatellite-based clusters of nests (Table 1).
The population at Bois de Peney had the highest number of clusters, as four genetically distinct groups of nests emerged when analysing the microsatellite genotypes with Geneland (p = 52.7%, Fig. 1 panel A.1) and Structure (Fig. 1 panel A.3). Three of these nest clusters were almost exclusively composed of individuals with the same mitochondrial haplotype, whereas one cluster had a high frequency of another haplotype (Fig. 1 panels A.2 and A.3). The population at La Dôle had three well-defined nest clusters based on microsatellite differentiation (p = 72.0%, Fig. 1 panels B.1 and B.3). The four mitochondrial haplotypes were unevenly distributed among these clusters: we detected only two haplotypes per cluster, and two haplotypes were restricted to single clusters (Fig. 1 panels B.2 and B.3). Finally, the population at Chalet à Roch had two rather ill-defined nest clusters based on microsatellite differentiation (p = 60.0%, Fig. 1 panels C.1 and C.3), with many of the individuals having a large fraction of their genotype assigned to the other cluster (Fig. 1 panel C.3). Two mitochondrial haplotypes intermingled throughout the transect, but had different frequencies in each cluster (Fig. 1 panels C.2 and C.3).
The hierarchical F-analyses of variance of the microsatellite genotypes revealed significant genetic differentiation among nests within clusters and among clusters within populations at our three study sites (Table 1). The lowest genetic differentiation among clusters was found at Chalet à Roch and the highest at La Dôle (Fcluster-pop = 3.7% and 10%, respectively). The significant Find-pop coefficient indicates non-random mating at the population level. This deficit of heterozygotes is consistent with a Wahlund effect over small spatial scales due to genetic differentiation among nests and among clusters. In contrast, Find-cluster was not significantly different from zero, which is in agreement with the hypothesis of random mating within clusters that underlies the way clusters are identified. This result suggests that clusters roughly correspond to breeding units in Hardy-Weinberg equilibrium, despite the significant genetic differentiation among nests. Finally, slightly negative and close to zero Find-nest coefficients are consistent with the expected excess of heterozygous within brood coupled with high number of queens that mate at random within nests.
In line with the significant genetic differentiation among nests and clusters, the relatedness among nestmate workers and the relatedness among clustermates measured with respect to the allele frequencies in the total population were significantly larger than zero in all three populations (Table 2). Interestingly, the relatedness among clustermates was fairly high, and accounted for most of the relatedness among nestmates measured with respect to the total population. Conversely, the relatedness among nestmates measured with respect to the cluster was moderate, but still significantly greater than zero in each population (Table 2). Together, these results reveal a multi-level pattern of genetic structure, with fairly high gene flow among neighbouring nests within clusters and restricted gene flow among clusters within unicolonial populations. This additional level of structure above the nest level results in a high relatedness among members of the same cluster.
The mitochondrial haplotypes showed strong and significant genetic differentiation among nests within clusters, and among clusters within populations (Table 1). The lowest genetic differentiation among clusters was also observed in the Chalet à Roch population (Table 1). Overall, the F-estimates for the mitochondrial haplotypes were four to six times greater than the ones for the nuclear makers, which indicates male-biased gene flow (see below; when both sexes disperse equally, a ratio Fst mitochondrial/Fst nuclear of approximately three is expected because of the smaller effective population size for the mitochondrial genome, ).
Gene flow through males and females
The contrasting pattern of genetic differentiation between nuclear and mitochondrial markers revealed that the effective gene flow among nests was significantly male-biased in all three populations (Table 3). This result is conservative, as the gene diversity was slightly higher for mitochondrial markers than for nuclear ones (Table 3), which, other things being equal, should tend to decrease the Fst for mitochondrial markers [46, 47]. Rough estimates of relative gene flow through each sex suggest that the gene flow through males was on average three to nine times higher than the one through females, depending on the population (Table 3). The overall number of genetically effective migrants suggests that long-distance gene flow is restricted, given the very high number of queens per nest and presumably high level of genetic exchanges between neighbouring nests.
The evolution and maintenance of unicoloniality critically depends on patterns of dispersal and population genetic structuring, which in turn affect the degree of relatedness among nestmates. Here, we provide genetic evidence for cryptic genetic clustering in three unicolonial populations of the native wood ant F. paralugubris. This clustering is associated with male-biased gene flow, local exchanges between neighbouring nests and restricted long-distance gene flow.
The Bayesian analyses of the microsatellite genotypes revealed the presence of genetic clusters of nests within each of the three populations. Moreover, the use of a new model that incorporates geographic information permitted us to visualize the spatial location of cryptic genetic discontinuities within populations. In line with these findings, the F-statistics showed a Wahlund effect at the population level, as expected from the microgeographical genetic differentiation among nest clusters. In contrast, the genotypes of workers within clusters were in Hardy-Weinberg equilibrium, as expected from the model assumptions of the clustering methods. This analysis suggests that the clusters correspond to randomly mating breeding populations, with many queens and fairly high levels of local gene flow through queens and males, whereas gene flow between clusters is restricted.
The structure in clusters of nests has important consequences for the functioning and evolution of unicolonial populations. The crucial point is that clusters of nests are genetically differentiated, which results in moderate but significant levels of relatedness among members of the same cluster. In addition, nests within clusters are also slightly genetically differentiated, which further increases the relatedness among nestmates. The differentiation among nests within clusters may be explained by a finite number of queens per nest coupled with restricted queen and worker dispersal [22, 42]. It is also possible that workers and queens tend to segregate along kin lines during colony budding .
The significant genetic structuring among nests within clusters and among clusters within populations in F. paralugubris contrasts with the absence of genetic differentiation within supercolonies of the extensively studied invasive and unicolonial argentine ant L. humile . The absence of differentiation within native supercolonies of L. humile might be linked to their small size and extensive queen movement , whereas the very large supercolonies in the introduced range are probably out of equilibrium and affected by repeated human transport or habitat disturbance .
Beside F. paralugubris, several non-invasive ant species had significant genetic structuring within unicolonial populations [25, 27–30, 54]. The finding of genetic structuring within native unicolonial populations and the resulting elevated relatedness among nestmates despite high queen numbers per nest beg the question of the potential impact of kin selection on the maintenance of altruism in these systems.
The restricted effective gene flow and ensuing genetic differentiation among distant nests within populations may promote the maintenance of altruism in two ways. First, elevated relatedness due to population viscosity may contribute to the maintenance of altruistic sterile workers . However, given the high number of queens and low relatedness among nestmates in our study populations, the long-term maintenance of altruism would require very large benefits to costs ratios, implying that the synergetic benefits of cooperation are huge (see equation 2.4 in ). Second, the maintenance of reproductive altruism could be explained by kin selection if competition occurs among distant, genetically differentiated nests, or among clusters [42, 55, 56]. Whether clusters actually compete for resources remains to be investigated. Competition among clusters does not seem to involve aggressive worker behaviour [33, 49, 57]. However, given the high number of queens and males that fly away from the nests in early summer, it seems likely that competition among clusters will occur after the dispersal flight, when males try to mate and queens seek to establish new colonies or to be adopted into nests [41, 58, 59].
The use of biparentally-inherited nuclear microsatellites and maternally-inherited mitochondrial markers allowed us to roughly estimate the amount of gene flow through males and females. Previous studies based on microsatellites had revealed limited long-distance gene flow and isolation by distance in the population at Chalet à Roch [22, 42, 60]. Our new analyses of three populations confirmed that long-distance gene flow is restricted. First, there was significant genetic structuring for both types of markers at all levels. Second, the number of genetically effective migrants among all nests within each population was moderate, given the high numbers of queens per nest and high rate of exchanges between neighbouring nests.
The comparatively much higher degree of differentiation at mitochondrial than nuclear markers further revealed that effective gene flow is male-biased in all three populations. These data indicate that queens are more philopatric and disperse less frequently or less far than males do. It is likely that male gene flow occurs mostly within clusters, which would explain both the Hardy-Weinberg equilibrium within clusters and significant genetic differentiation among clusters. Male-biased gene flow seems to be common in other polygnous ants in which queens and workers disperse by budding [18, 44, 61–63].
The presence of more than one mitochondrial haplotype in many nests suggests that queens are occasionally accepted into foreign nests. This result contrasts with kin selection models suggesting that colonies should only accept daughter queens . However, it does not imply that gene flow through females is high, because, once established, matriline associations can persist indefinitely and spread in the population if queen number is high, queens are adopted into their natal nests and nest budding occurs, as is the case in our study species [22, 42]. Multiple mitochondrial haplotypes per nest have been found in many ant species, which indicates that unrelated queens can associate in ants, be they unicolonial or not [44, 65, 66].
Overall, the data suggest a social organisation in large breeding units containing moderately related individuals. At the same time, the movement of workers among neighbouring nests and absence of aggression within populations [33, 49] are clear signs of unicoloniality. The evolutionary pathway to this type of social organisation with no aggression, mixing of workers among neighbouring nests and restricted long-distance dispersal of queens and males remains unclear, and has been suggested to depend on both ecological and genetical factors [7, 42, 67, 68].
The life-history of F. paralugubris is at least partly consistent with an 'ecological pathway' hypothesis by which unicoloniality is ultimately mediated by habitat saturation [7, 42, 67, 69]. F. paralugubris has a sessile life-history in a stable habitat. Workers build large, long lasting nests. They also tend aphids, and these socially-expandable resources may in turn lead to increased nest densities. This progressive habitat saturation makes independent colony foundation increasingly difficult, promoting the reacceptance of queens by the maternal colony and establishment of new nests by budding, which in turn leads to the formation of large networks of nests . In addition, the absence of aggression reduces the costs of territoriality, allowing higher worker densities and effective habitat monopolisation [16, 70, 71]. It has even been suggested that under certain ecological conditions, abandoning aggression while retaining some discrimination ability might be a first stage in the formation of supercolonies or unicolonial populations . In line with this argument, the absence of aggression was not due to a complete lack of discrimination in F. paralugubris and L. austriacus, as shown by the significantly longer antennation bouts among non-nestmates [7, 33].
This study revealed restricted, male-biased gene flow and cryptic genetic clustering within three native unicolonial populations of wood ants. This peculiar form of social organisation might be driven by progressive habitat saturation that favours intranidal mating, acceptance of daughter queens and colony budding. The positive relatedness among clustermates suggests that kin selection may still contribute to the maintenance of altruism in unicolonial populations as long as competition occurs among clusters, which seems likely during dispersal flights. Whether kin selection will be strong enough to maintain altruistic worker behaviour over evolutionary times remains an open question.
We focussed on three large populations of F. paralugubris in the Swiss Jura Mountains. We sampled workers from 20 nests at Bois de Peney, 20 nests at Chalet à Roch and 18 nests at La Dôle in July 2003. In addition, for the analysis of mitochondrial variability we sampled workers from two other populations in the Jura and two populations in the Swiss Alps (ten nests at Bois de Ban, seven nests at Champs Simon, nine nests at Château d'Oex, and two nests near Lac de l'Hongrin). All sites were situated between 950 and 1310 meters above sea level (see  for a map). We recorded the coordinates of each nest with a global positioning system (D-GPS; Garmin Ltd. Romsey, UK). The distances between nests ranged from 15 m to 76 km.
We used the microsatellite genotypes of eight workers per nest that we had already obtained for a previous analysis of the large-scale population genetic structure of F. paralugubris . The microsatellites include four loci developed for F. paralugubris (FL12, FL20, FL21, FL29, ) and seven loci originally developed for Formica exsecta (FE7, FE8, FE11, FE19, FE37, FE38, FE42, ). See  for details on microsatellite genotyping.
We screened for mitochondrial DNA variation by sequencing the cytochrome oxidase subunit I and II (COI-COII; 1368 bp) [GenBank: EU600788], the region of NADH dehydrogenase from subunit 4 to subunit 6 (ND4–ND6; 1497 bp with a gap of 100 bp) [GenBank: EU600793], and the region from ND6 to ND1, which includes the cytochrome b (2427 bp and 2430 bp) [GenBank: EU600789–EU600792]. The primers used to amplify the DNA are listed in Table 4. We designed primers for the COI-COII region on conserved regions in the alignment of sequences from Apis mellifera, Formica fusca and F. truncorum, except for the primer COI-RLR-F which came from Roehrdanz & Degrugillier . We designed primers for the ND6 to ND1 and ND4 to ND6 regions on conserved regions in the alignment of sequences from A. mellifera and F. exsecta, except for the primers Cytb-Fe which came from Liautard & Keller .
In the initial screen for polymorphism, we sequenced 11 workers sampled from 11 independent nests from the seven study populations. We extracted DNA from head and thorax with a standard phenol chloroform protocol. PCR products were purified using the QIAquick purification kit (Qiagen) and sequenced using the BigDye terminator ready-reaction kit and an ABI 3100 automated sequencer (Applied Biosystems). Sequence data were edited and compiled with the software Sequencher 3.0 (Gene Codes Corporation). To confirm the amplification of the target region of the mtDNA, the final sequences were aligned with the honeybee sequence  using the EMBL-EBI ClustalW v.1.83 website http://www.ebi.ac.uk/Tools/clustalw. We detected one deletion/insertion of three base pairs and two single nucleotide polymorphisms (SNPs) in the region between ND6 and ND1. We found no polymorphism in the other regions sequenced.
For the routine typing of the three polymorphic sites in the ND6 to ND1 region, we designed the specific primers CbM3, PB300 T/C and PB30 A/G (Table 4). The DNA of the workers was extracted with a Chelex protocol as described in Holzer et al. . To analyze the three base pair deletion/insertion, we designed new primers labelled with fluorescent dyes (CbM3), which amplified a fragment that was either 163 or 166 bp long. The fragment was scored on an ABI 377XL sequencer. To genotype the two SNPs PB300 T/C and PB30 A/G, we used specific locked nucleic acids (LNA) primers (Sigma-Proligo). The PCR products were electrophoresed on a standard ethidium bromide 1.5% agarose gel. We routinely scored the mitochondrial haplotypes in the three largest populations (Bois de Peney, Chalet à Roch and La Dôle), using the same eight individuals as for the microsatellites. In the other populations, we genotyped two individuals per nest.
To search for cryptic genetic structure among groups of nests within populations, we analysed the microsatellite data set with the methods implemented in the computer programs Geneland [52, 76] and Structure v. 2 . These individual-centred methods aim to detect cryptic genetic structure, admixed individuals and migrants on the basis of the fact that migrants will have multilocus genotypes that differ from the ones expected for native individuals .
Geneland divides up geo-referenced individuals into clusters corresponding to panmictic units without any a priori information on the shape and limits of these clusters [52, 76]. The set of individuals are attributed to "virtual" clusters represented by random unions of convex polygons of a Voronoi tessellation. Allele frequencies are estimated iteratively and the spatial organization that best corresponds to panmictic units is identified with a Markov chain Monte Carlo procedure. Simulations showed that this method correctly identifies unconnected spatial groups that belong to the same panmictic cluster, as well as large and uniform clusters corresponding to homogenous populations . We first determined the number of clusters (k) within the populations with a maximum number of clusters of ten. Each run was repeated ten times with 100000 iterations. To allow the assignment of individuals sharing the same nest to a different cluster, the 'uncertainty of coordinates' was set to 20, as neighbouring nests are often within a distance range of 20 m. Once k was determined, we re-run the program ten times with the fixed number of k inferred from the previous runs.
The computer program Structure detects the proportion of the genome of an individual originating from each inferred cluster without taking any spatial information into account [admixture analysis, ]. Individuals are assigned sequentially to clusters, or to two or more clusters if their genotypes indicate that they are admixed. We used the admixture model with a burn-in of 50000 iterations and 1000000 iterations to estimate the parameters. The procedure was repeated ten times for each k.
The relatedness among worker nestmates was estimated with the computer program Relatedness 5.0.8 . Nests where weighted equally. Standard errors and 95% confidence intervals were obtained by jackknifing over nests for populations, and over loci for clusters. Relatedness was measured with respect to allele frequencies in clusters and in populations. This double way to measure relatedness is justified by the fact that part of the competition may be local and occur only within clusters, whereas some degree of competition may occur at the population scale. Specifically, the relatedness rnestmates (population) and rclustermates (population) estimate the probability by which an allele is more likely to be shared among nest- and clustermates than expected from the allele frequency in the population to which the nest or cluster belongs. The relatedness rnestmates (cluster) is the probability by which an allele is more likely to be shared among nestmates than expected from the allele frequency in the cluster.
The genetic differentiation at the microsatellite markers was further analyzed in a three-level hierarchical F-analysis of variance, as implemented in the computer program Hierfstat [79, 80]. The three hierarchical levels were the individuals, nests, and clusters previously identified in the analysis of microsatellites with Geneland. We analysed each population separately. Find-nest and Find-cluster are the correlations of genes within individuals, as compared to random genes within the nest and cluster, respectively. Find-pop is the correlation of genes within individuals, as compared to random genes within the population. Fnest-cluster estimates the genetic differentiation among nests within clusters, and Fcluster-pop the genetic differentiation among clusters within populations. The genetic differentiation at the mitochondrial DNA was investigated in a two-level hierarchical analysis, as implemented in the computer program Arlequin . The two levels were the nests and clusters previously identified in the analysis of the microsatellites with Geneland.
The use of both nuclear and mitochondrial markers allowed us to roughly estimate the relative contribution of each sex to gene flow. We estimated the gene flow among nests within populations from the genetic differentiation among all nests within each population at the biparentally-inherited nuclear microsatellites (Fst nucl) and maternally-inherited mitochondrial (Fst mit) markers. Because more diverse genetic markers have lower Fst (the maximal possible value of Fst being given by 1-Hs) [46, 47], we compared the gene (or haplotype) diversity Hs of microsatellite and mitochondrial markers. The relative strength of gene flow through females or males was calculated according to Seppä et al. , with the effective gene flow of females estimated as N f m f ≈ (1 - Fst mit)/(2 Fst mit) , where N f is the number of females breeding in a population and m f the immigration rate. The gene flow of males was calculated as N m m m ≈ [(1 - Fst nucl)/(2 Fst nucl)] - [(1 - Fst mit)/Fst mit)] by applying Wright's island model  modified for haplodiploid species . Because our study system does not fulfil many of the assumptions of the island model, the estimates N f m f and N m m m should not be considered as actual number of migrants , but as relative estimates of gene flow through each sex.
Hamilton WD: The genetical evolution of social behaviour. II. J Theor Biol. 1964, 7: 17-52. 10.1016/0022-5193(64)90039-6.
Foster KR, Wenseleers T, Ratnieks FLW: Kin selection is the key to altruism. Trends Ecol Evol. 2006, 21: 57-60. 10.1016/j.tree.2005.11.020.
Hughes WOH, Oldroyd BP, Beekman M, Ratnieks FLW: Ancestral monogamy shows kin selection is key to the evolution of eusociality. Science. 2008, 320: 1213-1216. 10.1126/science.1156108.
Crozier RH: Advanced eusociality, kin selection and male haploidy. Aust J Entomol. 2008, 47: 2-8. 10.1111/j.1440-6055.2007.00621.x.
Hölldobler B, Wilson EO: The number of queens: an important trait in ant evolution. Naturwissenschaften. 1977, 64: 8-15. 10.1007/BF00439886.
Queller DC, Strassmann JE: Kin selection and social insects. Bioscience. 1998, 48: 165-175. 10.2307/1313262.
Steiner FM, Schlick-Steiner BC, Moder K, Stauffer C, Arthofer W, Buschinger A, Espadaler X, Christian E, Einfinger K, Lorbeer E, et al: Abandoning aggression but maintaining self-nonself discrimination as a first stage in ant supercolony formation. Curr Biol. 2007, 17: 1903-1907. 10.1016/j.cub.2007.09.061.
Crozier RH, Pamilo P: Evolution of social insect colonies: sex allocation and kin selection. 1996, Oxford: Oxford University Press
Pedersen JS, Krieger MJB, Vogel V, Giraud T, Keller L: Native supercolonies of unrelated individuals in the invasive Argentine ant. Evolution. 2006, 60: 782-791.
Keller L: Social life: the paradox of multiple-queen colonies. Trends Ecol Evol. 1995, 10: 355-360. 10.1016/S0169-5347(00)89133-8.
Bourke AFG, Franks NR: Social evolution in ants. 1995, Princeton: Princeton University Press
Ross KG: Molecular ecology of social behaviour: analyses of breeding systems and genetic structure. Mol Ecol. 2001, 10: 265-284. 10.1046/j.1365-294x.2001.01191.x.
Queller DC: Genetic relatedness in viscous populations. Evol Ecol. 1994, 8: 70-73. 10.1007/BF01237667.
West SA, Griffin AS, Gardner A: Evolutionary explanations for cooperation. Curr Biol. 2007, 17: R661-R672. 10.1016/j.cub.2007.06.004.
Giraud T, Pedersen JS, Keller L: Evolution of supercolonies: the Argentine ants of southern Europe. Proc Natl Acad Sci USA. 2002, 99: 6075-6079. 10.1073/pnas.092694199.
Holway DA, Suarez AV, Case TJ: Loss of intraspecific aggression in the success of a widespread invasive social insect. Science. 1998, 282: 949-952. 10.1126/science.282.5390.949.
Tsutsui ND, Suarez AV, Holway DA, Case TJ: Reduced genetic variation and the success of an invasive species. Proc Natl Acad Sci USA. 2000, 97: 5948-5953. 10.1073/pnas.100110397.
Ross KG, Shoemaker DD, Krieger MJB, DeHeer CJ, Keller L: Assessing genetic structure with multiple classes of molecular markers: A case study involving the introduced fire ant Solenopsis invicta. Mol Biol Evol. 1999, 16: 525-543.
Le Breton J, Delabie JHC, Chazeau J, Dejean A, Jourdan H: Experimental evidence of large-scale unicoloniality in the tramp ant Wasmannia auropunctata (Roger). J Insect Behav. 2004, 17: 263-271. 10.1023/B:JOIR.0000028575.28700.71.
Errard C, Delabie J, Jourdan H, Hefetz A: Intercontinental chemical variation in the invasive ant Wasmannia auropunctata (Roger) (Hymenoptera Formicidae): a key to the invasive success of a tramp species. Naturwissenschaften. 2005, 92: 319-323. 10.1007/s00114-005-0628-y.
Elias M, Rosengren R, Sundström L: Seasonal polydomy and unicoloniality in a polygynous population of the red wood ant Formica truncorum. Behav Ecol Sociobiol. 2005, 57: 339-349. 10.1007/s00265-004-0864-8.
Chapuisat M, Goudet J, Keller L: Microsatellites reveal high population viscosity and limited dispersal in the ant Formica paralugubris. Evolution. 1997, 51: 475-482. 10.2307/2411120.
Higashi S, Yamauchi K: Influence of a supercolonial ant Formica (Formica) yessensis Forel on the distribution of other ants in Ishikari Coast. Jap J Ecol. 1979, 29: 257-264.
Rosengren R, Pamilo P: The evolution of polygyny and polydomy in mound-building Formica ants. Acta Entomol Fennica. 1983, 42: 65-77.
Kümmerli R, Keller L: Contrasting population genetic structure for workers and queens in the putatively unicolonial ant Formica exsecta. Mol Ecol. 2007, 16: 4493-4503. 10.1111/j.1365-294X.2007.03514.x.
Yamauchi K, Kinomura K, Miyake S: Sociobiological studies of the polygynic ant Lasius sakagami. 1 General features of its polydomous system. Ins Soc. 1981, 28: 279-296. 10.1007/BF02223629.
Seppä P, Pamilo P: Gene flow and population viscosity in Myrmica ants. Heredity. 1995, 74: 200-209. 10.1038/hdy.1995.28.
Seppä P, Walin L: Sociogenetic organization of the red ant Myrmica rubra. Behav Ecol Sociobiol. 1996, 38: 207-217. 10.1007/s002650050234.
Hammen van der T, Pedersen JS, Boomsma JJ: Convergent development of low-relatedness supercolonies in Myrmica ants. Heredity. 2002, 89: 83-89. 10.1038/sj.hdy.6800098.
Pedersen JS, Boomsma JJ: Genetic analysis of colony structure in polydomous and polygynous ant populations. Biol J Linn Soc. 1999, 66: 115-144. 10.1111/j.1095-8312.1999.tb01919.x.
Steiner FM, Schlick-Steiner BC, Buschinger A: First record of unicolonial polygyny in Tetramorium cf. caespitum (Hymenoptera, Formicidae). Ins Soc. 2003, 50: 98-99. 10.1007/s000400300015.
van Zweden JS, Carew ME, Henshaw MT, Robson SKA, Crozier RH: Social and genetic structure of a supercolonial weaver ant, Polyrhachis robsoni, with dimorphic queens. Ins Soc. 2007, 54: 34-41. 10.1007/s00040-007-0909-x.
Holzer B, Chapuisat M, Kremer N, Finet C, Keller L: Unicoloniality, recognition and genetic differentiation in a native Formica ant. J Evol Biol. 2006, 19: 2031-2039. 10.1111/j.1420-9101.2006.01133.x.
Lehmann L, Ravigné V, Keller L: Population viscosity can promote the evolution of altruistic sterile helpers and eusociality. Proc R Soc B. 2008, 275: 1887-1895. 10.1098/rspb.2008.0276.
Rousset F: Genetic structure and selection in subdivided populations. 2004, Princeton: Princeton University Press
Taylor PD: Altruism in viscous populations – an inclusive fitness model. Evol Ecol. 1992, 6: 352-356. 10.1007/BF02270971.
Griffin AS, West SA: Kin selection: fact and fiction. Trends Ecol Evol. 2002, 17: 15-21. 10.1016/S0169-5347(01)02355-2.
Cherix D: Note préliminaire sur la structure, la phénologie et le régime alimentaire d'une super-colonie de Formica lugubris Zett. Ins Soc. 1980, 27: 226-236. 10.1007/BF02223666.
Cherix D: Pseudogynes (= sécrétergates) et répartition des individus à l'intérieur d'une fourmilière de Formica lugubris Zett. (Hymenoptera: Formicidae). Ins Soc. 1983, 30: 184-192. 10.1007/BF02223868.
Helanterä H, Sundström L: Worker reproduction in Formica ants. Am Nat. 2007, 170: E14-E25. 10.1086/518185.
Cherix D, Chautems D, Fletcher DJC, Fortelius W, Gris G, Keller L, Passera L, Rosengren R, Vargo EL, Walter F: Alternative reproductive strategies in Formica lugubris Zett. (Hymenoptera Formicidae). Ethol Ecol Evol. 1991, 61-66. Special Issue 1
Chapuisat M, Keller L: Extended family structure in the ant Formica paralugubris : the role of the breeding system. Behav Ecol Sociobiol. 1999, 46: 405-412. 10.1007/s002650050636.
Goodisman MAD, Ross KG: A test of queen recruitment models using nuclear and mitochondrial markers in the fire ant Solenopsis invicta. Evolution. 1998, 52: 1416-1422. 10.2307/2411311.
Seppä P, Gyllenstrand M, Corander J, Pamilo P: Coexistence of the social types: Genetic population structure in the ant Formica exsecta. Evolution. 2004, 58: 2462-2471.
Prugnolle F, de Meeus T: Inferring sex-biased dispersal from population genetic tools: a review. Heredity. 2002, 88: 161-165. 10.1038/sj.hdy.6800060.
Hedrick PW: Perspective: Highly variable loci and their interpretation in evolution and conservation. Evolution. 1999, 53: 313-318. 10.2307/2640768.
Hedrick PW: A standardized genetic differentiation measure. Evolution. 2005, 59: 1633-1638.
Whitlock MC, McCauley DE: Indirect measures of gene flow and migration: FST ≠ 1/(4Nm+1). Heredity. 1999, 82: 117-125. 10.1038/sj.hdy.6884960.
Chapuisat M, Bernasconi C, Hoehn S, Reuter M: Nestmate recognition in the unicolonial ant Formica paralugubris. Behav Ecol. 2005, 16: 15-19. 10.1093/beheco/arh128.
Excoffier L, Heckel G: Computer programs for population genetics data analysis: a survival guide. Nature Reviews Genetics. 2006, 7: 745-758. 10.1038/nrg1904.
Beaumont MA, Rannala B: The Bayesian revolution in genetics. Nature Reviews Genetics. 2004, 5: 251-261. 10.1038/nrg1318.
Guillot G, Estoup A, Mortier F, Cosson JF: A spatial statistical model for landscape genetics. Genetics. 2005, 170: 1261-1280. 10.1534/genetics.104.033803.
Seppä P, Fernandez-Escudero I, Gyllenstrand N, Pamilo P: Colony fission affects kinship in a social insect. Behav Ecol Sociobiol. 2008, 62: 589-597. 10.1007/s00265-007-0483-2.
Beye M, Neumann P, Chapuisat M, Pamilo P, Moritz RFA: Nestmate recognition and the genetic relatedness of nests in the ant Formica pratensis. Behav Ecol Sociobiol. 1998, 43: 67-72. 10.1007/s002650050467.
Queller DC: Does population viscosity promote kin selection?. Trends Ecol Evol. 1992, 7: 322-324. 10.1016/0169-5347(92)90120-Z.
Griffin AS, West SA, Buckling A: Cooperation and competition in pathogenic bacteria. Nature. 2004, 430: 1024-1027. 10.1038/nature02744.
Holzer B, Meunier J, Keller L, Chapuisat M: Stay or drift? Queen acceptance in the ant Formica paralugubris. Ins Soc. 2008, 55: 392-396. 10.1007/s00040-008-1017-2.
Chapuisat M: Mating frequency of ant queens with alternative dispersal strategies, as revealed by microsatellite analysis of sperm. Mol Ecol. 1998, 7: 1097-1105. 10.1046/j.1365-294x.1998.00422.x.
Cherix D, Fletcher DJC, Chautems D, Fortelius W, Gris G, Keller L, Rosengren R, Vargo E, Walter F: Attraction of the sexes in Formica lugubris Zett. (Hymenoptera: Formicidae). Ins Soc. 1993, 40: 319-324. 10.1007/BF01242368.
Pamilo P, Chautems D, Cherix D: Genetic differentiation of disjunct populations of the ants Formica aquilonia and Formica lugubris in Europe. Ins Soc. 1992, 39: 15-29. 10.1007/BF01240528.
Sundström L, Seppä P, Pamilo P: Genetic population structure and dispersal patterns in Formica ants – a review. Ann Zool Fenn. 2005, 42: 163-177.
Rüppell O, Stratz M, Baier B, Heinze J: Mitochondrial markers in the ant Leptothorax rugatulus reveal the population genetic consequences of female philopatry at different hierarchical levels. Mol Ecol. 2003, 12: 795-801. 10.1046/j.1365-294X.2003.01769.x.
Seppä P, Fernandez-Escudero I, Gyllenstrand N, Pamilo P: Obligatory female philopatry affects genetic population structure in the ant Proformica longiseta. Ins Soc. 2006, 53: 362-368. 10.1007/s00040-006-0882-9.
Nonacs P: Queen number in colonies of social Hymenoptera as a kin-selected adaptation. Evolution. 1988, 42: 566-580. 10.2307/2409040.
Heinze J, Keller L: Alternative reproductive strategies: a queen perspective in ants. Trends Ecol Evol. 2000, 15: 508-512. 10.1016/S0169-5347(00)01995-9.
Liautard C, Keller L: Restricted effective queen dispersal at a microgeographic scale in polygynous populations of the ant Formica exsecta. Evolution. 2001, 55: 2484-2492.
Chapman RE, Bourke AFG: The influence of sociality on the conservation biology of social insects. Ecol Lett. 2001, 4: 650-662. 10.1046/j.1461-0248.2001.00253.x.
Jackson DE: Social evolution: pathways to ant unicoloniality. Curr Biol. 2007, 17: R1063-R1064. 10.1016/j.cub.2007.10.031.
Ross KG, Keller L: Ecology and evolution of social organization: insights from fire ants and other highly eusocial insects. Annu Rev Ecol Syst. 1995, 26: 631-656. 10.1146/annurev.es.26.110195.003215.
Debout G, Schatz B, Elias M, McKey D: Polydomy in ants: what we know, what we think we know, and what remains to be done. Biol J Linn Soc. 2007, 90: 319-348. 10.1111/j.1095-8312.2007.00728.x.
Holway DA, Suarez AV: Colony-structure variation and interspecific competitive ability in the invasive Argentine ant. Oecologia. 2004, 138: 216-222. 10.1007/s00442-003-1414-1.
Chapuisat M: Characterization of microsatellite loci in Formica lugubris B and their variability in other ant species. Mol Ecol. 1996, 5: 599-601. 10.1111/j.1365-294X.1996.tb00354.x.
Gyllenstrand N, Gertsch PJ, Pamilo P: Polymorphic microsatellite DNA markers in the ant Formica exsecta. Mol Ecol Notes. 2002, 2: 67-69. 10.1046/j.1471-8286.2002.00152.x.
Roehrdanz RL, Degrugillier ME: Long sections of mitochondrial DNA amplified from fourteen orders of insects using conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1998, 91: 771-778.
Crozier RH, Crozier YC: The mitochondrial genome of the honeybee Apis mellifera : complete sequence and genome organization. Genetics. 1993, 133: 97-117.
Guillot G, Mortier F, Estoup A: GENELAND: a computer package for landscape genetics. Mol Ecol Notes. 2005, 5: 712-715. 10.1111/j.1471-8286.2005.01031.x.
Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.
Queller DC, Goodnight KF: Estimating relatedness using genetic markers. Evolution. 1989, 242: 258-275. 10.2307/2409206.
de Meeus T, Goudet J: A step-by-step tutorial to use HierFstat to analyse populations hierarchically structured at multiple levels. Infection Genetics and Evolution. 2007, 7: 731-735. 10.1016/j.meegid.2007.07.005.
Goudet J: HIERFSTAT, a package for R to compute and test hierarchical F -statistics. Mol Ecol Notes. 2005, 5: 184-186. 10.1111/j.1471-8286.2004.00828.x.
Schneider S, Roessli D, Excoffier L: Arlequin: a software for population genetics data analysis, version 2.000. 2000, Genetics and Biometry Lab, Dept. of Anthropology, University of Geneva, 2.000
Ennos RA: Estimating the relative rates of pollen and seed migration among plant populations. Heredity. 1994, 72: 250-259. 10.1038/hdy.1994.35.
Wright S: The genetical structure of populations. Ann Eug. 1951, 15: 323-354.
Berg LM, Lascoux M, Pamilo P: The infinite island model with sex-differentiated gene flow. Heredity. 1998, 81: 63-68. 10.1046/j.1365-2540.1998.00358.x.
We thank Corinne Peter for help in the laboratory, and Claudie Doums, Thomas Broquet, Valérie Vogel and three anonymous reviewers for comments on the manuscript. This study was supported by several grants from the Swiss National Science Foundation to LK and MC.
BH, LK and MC jointly conceived, designed and coordinated the study. BH carried out the sampling and laboratory analyses, analyzed the data and wrote a first draft of the manuscript. LK and MC supervised the analyses and MC finalized the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.