Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

Genetic clusters and sex-biased gene flow in a unicolonial Formica ant

Barbara Holzer, Laurent Keller and Michel Chapuisat*

Author Affiliations

Department of Ecology and Evolution, University of Lausanne, Biophore, UNIL-Sorge, 1015 Lausanne, Switzerland

For all author emails, please log on.

BMC Evolutionary Biology 2009, 9:69  doi:10.1186/1471-2148-9-69


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/9/69


Received:12 November 2008
Accepted:31 March 2009
Published:31 March 2009

© 2009 Holzer et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

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.

Results

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.

Conclusion

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.

Background

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 [5].

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 [6]. 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 [18], 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 [9]. 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 [19], whereas L. humile forms supercolonies of much smaller size than in the introduced range [9]. If these small supercolonies are closed breeding units that compete which each other, explaining the evolution and stability of unicoloniality becomes less problematic [9].

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 [31] and Polyrhachis [32]. 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 [1] 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 [34].

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 [40]. 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 [41]. 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 [41]. 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 [41]. 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 [48].

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 [42]. 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 [52]. 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.

Results

Haplotype distribution

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.

thumbnailFigure 1. Genetic clustering of nests. The cryptic genetic structure of three populations of F. paralugubris at Bois de Peney (A), La Dôle (B) and Chalet à Roch (C). Panels 1 (A.1, B.1 and C.1) show synthetic maps of the clusters detected when analysing the microsatellite genotypes with the program Geneland. Large numbers indicate the cluster, grey tones the cluster area, and small numbers the sampled nests. Symbols in color represent sampled nests belonging to the same cluster, and small black dots indicate nests in the vicinity of sampled nests. The asterisk in the population Chalet à Roch (C.1) indicates a nest with a posterior probability close to 0.5 to belong to one or the other cluster. Panels 2 show the distribution and frequency per nest of the mitochondrial haplotypes (violet: H1; orange: H2; white: H3; black: H4). The upper parts of panels 3 (Mst) show the population structure detected when analysing the microsatellite genotypes with the program Structure. Each individual is represented by a vertical line, which is partitioned into k coloured segments that represent the individuals's estimated membership fractions in k clusters. The bottom parts of panels 3 (Haplotypes) show the haplotype for each individual, given as a violet (H1), orange (H2), white (H3) or black (H4) bar. Thin lines in panels 3 separate individuals from different nests. Microsatellite clusters in panels 2 and 3 correspond to the ones detected with Geneland in panels 1.

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).

Table 1. Genetic differentiation among nests and among clusters

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.

Table 2. Estimates of genetic relatedness

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, [44]).

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.

Table 3. Sex specific gene flow

Discussion

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 [53].

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 [15]. The absence of differentiation within native supercolonies of L. humile might be linked to their small size and extensive queen movement [9], whereas the very large supercolonies in the introduced range are probably out of equilibrium and affected by repeated human transport or habitat disturbance [15].

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 [34]. 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 [34]). 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 [64]. 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 [42]. 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 [7]. 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].

Conclusion

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.

Methods

Field collection

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 [33] 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.

Genetic analysis

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 [33]. The microsatellites include four loci developed for F. paralugubris (FL12, FL20, FL21, FL29, [72]) and seven loci originally developed for Formica exsecta (FE7, FE8, FE11, FE19, FE37, FE38, FE42, [73]). See [33] 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: EU600789EU600792]. 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 [74]. 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 [66].

Table 4. Primers to amplify the mitochondrial DNA of F. paralugubris

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 [75] using the EMBL-EBI ClustalW v.1.83 website http://www.ebi.ac.uk/Tools/clustalw webcite. 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. [33]. 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.

Statistical analysis

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 [77]. 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 [50].

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 [52]. 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, [77]]. 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 [78]. 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 [81]. 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. [44], with the effective gene flow of females estimated as Nf mf ≈ (1 - Fst mit)/(2 Fst mit) [82], where Nf is the number of females breeding in a population and mf the immigration rate. The gene flow of males was calculated as Nm mm ≈ [(1 - Fst nucl)/(2 Fst nucl)] - [(1 - Fst mit)/Fst mit)] by applying Wright's island model [83] modified for haplodiploid species [84]. Because our study system does not fulfil many of the assumptions of the island model, the estimates Nf mf and Nm mm should not be considered as actual number of migrants [48], but as relative estimates of gene flow through each sex.

Authors' contributions

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.

Acknowledgements

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.

References

  1. Hamilton WD: The genetical evolution of social behaviour. II.

    J Theor Biol 1964, 7:17-52. PubMed Abstract | Publisher Full Text OpenURL

  2. Foster KR, Wenseleers T, Ratnieks FLW: Kin selection is the key to altruism.

    Trends Ecol Evol 2006, 21:57-60. PubMed Abstract | Publisher Full Text OpenURL

  3. 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. PubMed Abstract | Publisher Full Text OpenURL

  4. Crozier RH: Advanced eusociality, kin selection and male haploidy.

    Aust J Entomol 2008, 47:2-8. Publisher Full Text OpenURL

  5. Hölldobler B, Wilson EO: The number of queens: an important trait in ant evolution.

    Naturwissenschaften 1977, 64:8-15. Publisher Full Text OpenURL

  6. Queller DC, Strassmann JE: Kin selection and social insects.

    Bioscience 1998, 48:165-175. Publisher Full Text OpenURL

  7. 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. PubMed Abstract | Publisher Full Text OpenURL

  8. Crozier RH, Pamilo P: Evolution of social insect colonies: sex allocation and kin selection. Oxford: Oxford University Press; 1996.

  9. 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. PubMed Abstract OpenURL

  10. Keller L: Social life: the paradox of multiple-queen colonies.

    Trends Ecol Evol 1995, 10:355-360. Publisher Full Text OpenURL

  11. Bourke AFG, Franks NR: Social evolution in ants. Princeton: Princeton University Press; 1995. OpenURL

  12. Ross KG: Molecular ecology of social behaviour: analyses of breeding systems and genetic structure.

    Mol Ecol 2001, 10:265-284. PubMed Abstract | Publisher Full Text OpenURL

  13. Queller DC: Genetic relatedness in viscous populations.

    Evol Ecol 1994, 8:70-73. Publisher Full Text OpenURL

  14. West SA, Griffin AS, Gardner A: Evolutionary explanations for cooperation.

    Curr Biol 2007, 17:R661-R672. PubMed Abstract | Publisher Full Text OpenURL

  15. Giraud T, Pedersen JS, Keller L: Evolution of supercolonies: the Argentine ants of southern Europe.

    Proc Natl Acad Sci USA 2002, 99:6075-6079. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Holway DA, Suarez AV, Case TJ: Loss of intraspecific aggression in the success of a widespread invasive social insect.

    Science 1998, 282:949-952. PubMed Abstract | Publisher Full Text OpenURL

  17. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. 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. PubMed Abstract | Publisher Full Text OpenURL

  19. 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. Publisher Full Text OpenURL

  20. 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. PubMed Abstract | Publisher Full Text OpenURL

  21. 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. Publisher Full Text OpenURL

  22. Chapuisat M, Goudet J, Keller L: Microsatellites reveal high population viscosity and limited dispersal in the ant Formica paralugubris.

    Evolution 1997, 51:475-482. Publisher Full Text OpenURL

  23. 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. OpenURL

  24. Rosengren R, Pamilo P: The evolution of polygyny and polydomy in mound-building Formica ants.

    Acta Entomol Fennica 1983, 42:65-77. OpenURL

  25. 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. PubMed Abstract | Publisher Full Text OpenURL

  26. 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. Publisher Full Text OpenURL

  27. Seppä P, Pamilo P: Gene flow and population viscosity in Myrmica ants.

    Heredity 1995, 74:200-209. Publisher Full Text OpenURL

  28. Seppä P, Walin L: Sociogenetic organization of the red ant Myrmica rubra.

    Behav Ecol Sociobiol 1996, 38:207-217. Publisher Full Text OpenURL

  29. Hammen T, Pedersen JS, Boomsma JJ: Convergent development of low-relatedness supercolonies in Myrmica ants.

    Heredity 2002, 89:83-89. PubMed Abstract | Publisher Full Text OpenURL

  30. Pedersen JS, Boomsma JJ: Genetic analysis of colony structure in polydomous and polygynous ant populations.

    Biol J Linn Soc 1999, 66:115-144. Publisher Full Text OpenURL

  31. Steiner FM, Schlick-Steiner BC, Buschinger A: First record of unicolonial polygyny in Tetramorium cf. caespitum (Hymenoptera, Formicidae).

    Ins Soc 2003, 50:98-99. Publisher Full Text OpenURL

  32. 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. Publisher Full Text OpenURL

  33. 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. PubMed Abstract | Publisher Full Text OpenURL

  34. 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. PubMed Abstract | Publisher Full Text OpenURL

  35. Rousset F: Genetic structure and selection in subdivided populations. Princeton: Princeton University Press; 2004. OpenURL

  36. Taylor PD: Altruism in viscous populations – an inclusive fitness model.

    Evol Ecol 1992, 6:352-356. Publisher Full Text OpenURL

  37. Griffin AS, West SA: Kin selection: fact and fiction.

    Trends Ecol Evol 2002, 17:15-21. Publisher Full Text OpenURL

  38. 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. Publisher Full Text OpenURL

  39. 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. Publisher Full Text OpenURL

  40. Helanterä H, Sundström L: Worker reproduction in Formica ants.

    Am Nat 2007, 170:E14-E25. PubMed Abstract | Publisher Full Text OpenURL

  41. 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, (Special Issue 1):61-66. OpenURL

  42. 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. Publisher Full Text OpenURL

  43. 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. Publisher Full Text OpenURL

  44. 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. PubMed Abstract OpenURL

  45. Prugnolle F, de Meeus T: Inferring sex-biased dispersal from population genetic tools: a review.

    Heredity 2002, 88:161-165. PubMed Abstract | Publisher Full Text OpenURL

  46. Hedrick PW: Perspective: Highly variable loci and their interpretation in evolution and conservation.

    Evolution 1999, 53:313-318. Publisher Full Text OpenURL

  47. Hedrick PW: A standardized genetic differentiation measure.

    Evolution 2005, 59:1633-1638. PubMed Abstract OpenURL

  48. Whitlock MC, McCauley DE: Indirect measures of gene flow and migration: FST ≠ 1/(4Nm+1).

    Heredity 1999, 82:117-125. PubMed Abstract | Publisher Full Text OpenURL

  49. Chapuisat M, Bernasconi C, Hoehn S, Reuter M: Nestmate recognition in the unicolonial ant Formica paralugubris.

    Behav Ecol 2005, 16:15-19. Publisher Full Text OpenURL

  50. Excoffier L, Heckel G: Computer programs for population genetics data analysis: a survival guide.

    Nature Reviews Genetics 2006, 7:745-758. PubMed Abstract | Publisher Full Text OpenURL

  51. Beaumont MA, Rannala B: The Bayesian revolution in genetics.

    Nature Reviews Genetics 2004, 5:251-261. PubMed Abstract | Publisher Full Text OpenURL

  52. Guillot G, Estoup A, Mortier F, Cosson JF: A spatial statistical model for landscape genetics.

    Genetics 2005, 170:1261-1280. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Seppä P, Fernandez-Escudero I, Gyllenstrand N, Pamilo P: Colony fission affects kinship in a social insect.

    Behav Ecol Sociobiol 2008, 62:589-597. Publisher Full Text OpenURL

  54. 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. Publisher Full Text OpenURL

  55. Queller DC: Does population viscosity promote kin selection?

    Trends Ecol Evol 1992, 7:322-324. Publisher Full Text OpenURL

  56. Griffin AS, West SA, Buckling A: Cooperation and competition in pathogenic bacteria.

    Nature 2004, 430:1024-1027. PubMed Abstract | Publisher Full Text OpenURL

  57. Holzer B, Meunier J, Keller L, Chapuisat M: Stay or drift? Queen acceptance in the ant Formica paralugubris.

    Ins Soc 2008, 55:392-396. Publisher Full Text OpenURL

  58. Chapuisat M: Mating frequency of ant queens with alternative dispersal strategies, as revealed by microsatellite analysis of sperm.

    Mol Ecol 1998, 7:1097-1105. Publisher Full Text OpenURL

  59. 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. Publisher Full Text OpenURL

  60. 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. Publisher Full Text OpenURL

  61. 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. OpenURL

  62. 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. PubMed Abstract | Publisher Full Text OpenURL

  63. 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. Publisher Full Text OpenURL

  64. Nonacs P: Queen number in colonies of social Hymenoptera as a kin-selected adaptation.

    Evolution 1988, 42:566-580. Publisher Full Text OpenURL

  65. Heinze J, Keller L: Alternative reproductive strategies: a queen perspective in ants.

    Trends Ecol Evol 2000, 15:508-512. PubMed Abstract | Publisher Full Text OpenURL

  66. 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. PubMed Abstract OpenURL

  67. Chapman RE, Bourke AFG: The influence of sociality on the conservation biology of social insects.

    Ecol Lett 2001, 4:650-662. Publisher Full Text OpenURL

  68. Jackson DE: Social evolution: pathways to ant unicoloniality.

    Curr Biol 2007, 17:R1063-R1064. PubMed Abstract | Publisher Full Text OpenURL

  69. 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. Publisher Full Text OpenURL

  70. 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. Publisher Full Text OpenURL

  71. Holway DA, Suarez AV: Colony-structure variation and interspecific competitive ability in the invasive Argentine ant.

    Oecologia 2004, 138:216-222. PubMed Abstract | Publisher Full Text OpenURL

  72. Chapuisat M: Characterization of microsatellite loci in Formica lugubris B and their variability in other ant species.

    Mol Ecol 1996, 5:599-601. PubMed Abstract | Publisher Full Text OpenURL

  73. Gyllenstrand N, Gertsch PJ, Pamilo P: Polymorphic microsatellite DNA markers in the ant Formica exsecta.

    Mol Ecol Notes 2002, 2:67-69. Publisher Full Text OpenURL

  74. 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. OpenURL

  75. Crozier RH, Crozier YC: The mitochondrial genome of the honeybee Apis mellifera : complete sequence and genome organization.

    Genetics 1993, 133:97-117. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  76. Guillot G, Mortier F, Estoup A: GENELAND: a computer package for landscape genetics.

    Mol Ecol Notes 2005, 5:712-715. Publisher Full Text OpenURL

  77. Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data.

    Genetics 2000, 155:945-959. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Queller DC, Goodnight KF: Estimating relatedness using genetic markers.

    Evolution 1989, 242:258-275. Publisher Full Text OpenURL

  79. 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. Publisher Full Text OpenURL

  80. Goudet J: HIERFSTAT, a package for R to compute and test hierarchical F -statistics.

    Mol Ecol Notes 2005, 5:184-186. Publisher Full Text OpenURL

  81. Schneider S, Roessli D, Excoffier L: Arlequin: a software for population genetics data analysis, version 2.000. 2.000th edition. Genetics and Biometry Lab, Dept. of Anthropology, University of Geneva; 2000.

  82. Ennos RA: Estimating the relative rates of pollen and seed migration among plant populations.

    Heredity 1994, 72:250-259. Publisher Full Text OpenURL

  83. Wright S: The genetical structure of populations.

    Ann Eug 1951, 15:323-354. OpenURL

  84. Berg LM, Lascoux M, Pamilo P: The infinite island model with sex-differentiated gene flow.

    Heredity 1998, 81:63-68. Publisher Full Text OpenURL