Genetic structure of native ant supercolonies varies in space and time

Ant supercolonies are the largest cooperative units known in nature. They consist of networks of interconnected nests with hundreds of reproductive queens, where individuals move freely between nests, cooperate across nest boundaries and show little aggression towards non‐nestmates. The combination of high queen numbers and free mixing of workers, queens and brood between nests results in extremely low nestmate relatedness. In such low‐relatedness societies, cooperative worker behaviour appears maladaptive because it may aid random individuals instead of relatives. Here, we provide a comprehensive picture of genetic substructure in supercolonies of the native wood ant Formica aquilonia using traditional population genetic as well as network analysis methods. Specifically, we test for spatial and temporal variation in genetic structure of different classes of individuals within supercolonies and analyse the role of worker movement in determining supercolony genetic networks. We find that relatedness within supercolonies is low but positive when viewed on a population level, which may be due to limited dispersal of individuals and/or ecological factors such as nest site limitation and competition against conspecifics. Genetic structure of supercolonies varied with both sample class and sampling time point, which indicates that mobility of individuals varies according to both caste and season and suggests that generalizing has to be carried out with caution in studies of supercolonial species. Overall, our analysis provides novel evidence that native wood ant supercolonies exhibit fine‐scale genetic substructure, which may explain the maintenance of cooperation in these low‐relatedness societies.     

Limited dispersal and an unexpected aggression pattern in a native supercolonial ant

Understanding how social groups function requires studies on how individuals move across the landscape and interact with each other. Ant supercolonies are extreme cooperative units that may consist of thousands of interconnected nests, and their individuals cooperate over large spatial scales. However, the inner structure of suggested supercolonial (or unicolonial) societies has rarely been extensively studied using both genetic and behavioral analyses. We describe a dense supercolony‐like aggregation of more than 1,300 nests of the ant Formica (Coptoformica) pressilabris. We performed aggression assays and found that, while aggression levels were generally low, there was some aggression within the assumed supercolony. The occurrence of aggression increased with distance from the focal nest, in accordance with the genetically viscous population structure we observe by using 10 DNA microsatellite markers. However, the aggressive interactions do not follow any clear pattern that would allow specifying colony borders within the area. The genetic data indicate limited gene flow within and away from the supercolony. Our results show that a Formica supercolony is not necessarily a single unit but can be a more fluid mosaic of aggressive and amicable interactions instead, highlighting the need to study internest interactions in detail when describing supercolonies.     

Genetic consequences of natal dispersal in the colonial lesser kestrel

Dispersal is a life-history trait that plays a fundamental role in population dynamics, influencing evolution, species distribution, and the genetics and structure of populations. In spite of the fact that dispersal has been hypothesized to be an efficient behavioural mechanism to avoid inbreeding, the expected relationship between dispersal and mate relatedness still remains controversial. Here, we examine the genetic consequences of natal dispersal, namely the higher chance of obtaining genetically less similar mates as a result of moving from natal to breeding sites, in a lesser kestrel (Falco naumanni) population. Relatedness between individuals tended to decrease with distance between their breeding colonies, indicating that the study population follows an 'isolation-by-distance' pattern of spatial genetic structure. Such a fine-scale genetic structure generates a scenario in which individuals can potentially increase the chance of obtaining genetically less similar mates by dispersing over larger distances from their natal colony. Using dispersal information and genotypic data, we showed that mate relatedness decreased with natal dispersal distance, an effect that remained significant both while including and excluding philopatric individuals from the data set. These results, together with the well known detrimental consequences of reduced genetic diversity in the study population, suggest that dispersal may have evolved, at least in part, to avoid the negative fitness consequences of mating with genetically similar individuals.     

Restricted effective queen dispersal at a microgeographic scale in polygynous populations of the ant Formica exsecta

Ecological constraints on effective dispersal have been suggested to be a key factor influencing social evolution in animal societies as well as the shift from single queen colonies (monogyny) to multiple queen colonies (polygyny) in ants. However, little is known about the effective dispersal patterns of ant queens. Here we investigate the microgeographic genetic structure of mitochondrial haplotypes in polygynous populations of the ant Formica exsecta, both between pastures and among nests within pastures. An analysis of molecular variance revealed a very high genetic differentiation (phiST = 0.72) between pastures, indicating that queens rarely disperse successfully between pastures, despite the fact that pastures were sometimes as close as 1 km. Most of the pastures contained only a single haplotype, and haplotypes were frequently distinct between nearby pastures and even between groups of nests within the same pasture. In the three pastures that contained several haplotypes, haplotypes were not randomly distributed, the genetic differentiation between nests being phiST = 0.17, 0.52, and 0.69. This indicates that most queens are recruited within their parental colonies. However, a large proportion of nests contained more than one haplotype, demonstrating that colonies will sometimes accept foreign queens. The relatedness of mitochondrial genes among nestmates varied between 0.62 and 0.75 when relatedness was measured within each pasture and ranged between 0.72 and 1.0 when relatedness was assessed with all pastures as a reference population. Neighboring nests were more genetically similar than distant ones, and there was significant isolation by distance. This pattern may be due to new nests being formed by budding or by limited effective queen dispersal, probably on foot between neighboring nests. These results show that effective queen dispersal is extremely restricted even at a small geographical scale, a pattern consistent with the idea that ecological constraints are an important selective force leading to the evolution and maintenance of polygyny.     

Societies Drifting Apart? Behavioural,Genetic and Chemical Differentiation between Supercolonies in the Yellow Crazy Ant Anoplolepis gracilipes

Background

In populations of most social insects, gene flow is maintained through mating between reproductive individuals from different colonies in periodic nuptial flights followed by dispersal of the fertilized foundresses. Some ant species, however, form large polygynous supercolonies, in which mating takes place within the maternal nest (intranidal mating) and fertilized queens disperse within or along the boundary of the supercolony, leading to supercolony growth (colony budding). As a consequence, gene flow is largely confined within supercolonies. Over time, such supercolonies may diverge genetically and, thus, also in recognition cues (cuticular hydrocarbons, CHC''s) by a combination of genetic drift and accumulation of colony-specific, neutral mutations.

Methodology/Principal Findings

We tested this hypothesis for six supercolonies of the invasive ant Anoplolepis gracilipes in north-east Borneo. Within supercolonies, workers from different nests tolerated each other, were closely related and showed highly similar CHC profiles. Between supercolonies, aggression ranged from tolerance to mortal encounters and was negatively correlated with relatedness and CHC profile similarity. Supercolonies were genetically and chemically distinct, with mutually aggressive supercolony pairs sharing only 33.1%±17.5% (mean ± SD) of their alleles across six microsatellite loci and 73.8%±11.6% of the compounds in their CHC profile. Moreover, the proportion of alleles that differed between supercolony pairs was positively correlated to the proportion of qualitatively different CHC compounds. These qualitatively differing CHC compounds were found across various substance classes including alkanes, alkenes and mono-, di- and trimethyl-branched alkanes.

Conclusions

We conclude that positive feedback between genetic, chemical and behavioural traits may further enhance supercolony differentiation through genetic drift and neutral evolution, and may drive colonies towards different evolutionary pathways, possibly including speciation.     

MICROSATELLITES REVEAL HIGH POPULATION VISCOSITY AND LIMITED DISPERSAL IN THE ANT FORMICA PARALUGUBRIS

We used microsatellites to study the fine-scale genetic structure of a highly polygynous and largely unicolonial population of the ant Formica paralugubris. Genetic data indicate that long-distance gene flow between established nests is limited and new queens are primarily recruited from within their natal nest. Most matings occur between nestmates and are random at this level. In the center of the study area, budding and permanent connections between nests result in strong population viscosity, with close nests being more similar genetically than distant nests. In contrast, nests located outside of this supercolony show no isolation by distance, suggesting that they have been initiated by queens that participated in mating flights rather than by budding from nearby nests in our sample population. Recruitment of nestmates as new reproductive individuals and population viscosity in the supercolony increase genetic differentiation between nests. This in turn inflates relatedness estimates among worker nestmates (r = 0.17) above what is due to close pedigree links. Local spatial genetic differentiation may favor the maintenance of altruism when workers raise queens that will disperse on foot and compete with less related queens from neighboring nests or disperse on the wing and compete with unrelated queens.     

Weak large‐scale population genetic structure in a philopatric seabird,the European Shag Phalacrocorax aristotelis

Quantifying population genetic structure is fundamental to testing hypotheses regarding gene flow, population divergence and dynamics across large spatial scales. In species with highly mobile life‐history stages, where it is unclear whether such movements translate into effective dispersal among discrete philopatric breeding populations, this approach can be particularly effective. We used seven nuclear microsatellite loci and mitochondrial DNA (ND2) markers to quantify population genetic structure and variation across 20 populations (447 individuals) of one such species, the European Shag, spanning a large geographical range. Despite high breeding philopatry, rare cross‐sea movements and recognized subspecies, population genetic structure was weak across both microsatellites and mitochondrial markers. Furthermore, although isolation‐by‐distance was detected, microsatellite variation provided no evidence that open sea formed a complete barrier to effective dispersal. These data suggest that occasional long‐distance, cross‐sea movements translate into gene flow across a large spatial scale. Historical factors may also have shaped contemporary genetic structure: cluster analyses of microsatellite data identified three groups, comprising colonies at southern, mid‐ and northern latitudes, and similar structure was observed at mitochondrial loci. Only one private mitochondrial haplotype was found among subspecies, suggesting that this current taxonomic subdivision may not be mirrored by genetic isolation.     

Low relatedness among cooperatively breeding workers of the greenhead ant Rhytidoponera metallica

The greenhead ant Rhytidoponera metallica has long been recognized as posing a potential challenge to kin selection theory, because it has large queenless colonies where apparently many of the morphological workers are mated and reproducing. However, this species has never been studied genetically and important elements of its breeding system and kin structure remain uncertain. We used microsatellite markers to measure the relatedness among nestmates, unravel the fine‐scale population genetic structure, and infer the breeding system of R. metallica. The genetic relatedness among worker nestmates is very low but significantly greater than zero (r=0.082 ± 0.015), which demonstrates that nests contain many distantly related breeders. The inbreeding coefficient is very close to and not significantly different from zero, indicating random mating and lack of microgeographic genetic differentiation. On average, closely located nests are not more similar genetically than distant nests, which is surprising, as new colonies form by budding and female dispersal is restricted. Lack of inbreeding and absence of population viscosity indicates high gene flow mediated by males. Overall, the genetic pattern detected in R. metallica suggests that a high number of moderately related workers mate with unrelated males from distant nests. This breeding system results in the lowest relatedness among nestmates reported for social insect species where breeders and helpers are not morphologically differentiated.     

Small-scale spatial genetic structure in an ant species with sex-biased dispersal

In a population of the monogynous, polyandrous ant Cataglyphis cursor , we analysed the spatial genetic structure of queens, colony fathers and workers at a microgeographical scale to infer the extent of sex-biased dispersal and to assess the impact of limited dispersal on the patterns of relatedness within the colony. To this end, four microsatellite markers were scored for the queen and an average of 26 workers from each of 35 mapped colonies. We used pair-wise kinship coefficients between all pairs of genotypes, including the reconstructed colony father genotypes (1) to test and quantify isolation by distance patterns within each sex or caste through the analysis of kinship–distance curves, and (2) to compute the average relatedness between categories of colony members. The kinship–distance curve was much steeper for colony queens than colony fathers, indicating male-biased dispersal. However, colony fathers also displayed a non-random spatial genetic structure, so that even males show some dispersal limitation at the scale of the population, which extends over less than 250 m. The degree of relatedness between the different sexes and castes of colonies was well predicted from the number of mates per queen and the inbreeding of queens, and the impact of limited dispersal was very weak at this scale of observation. We discuss the interest of kinship–distance curves to assess sex-biased dispersal on a local scale and we compare our results with large-scale analyses of genetic structure in Cataglyphis cursor and other monogynous ant species.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 93 , 465–473.     

Population structure and intraspecific aggression in the invasive ant species Anoplolepis gracilipes in Malaysian Borneo

Invasive species are one of the main sources of the ongoing global loss of biodiversity. Invasive ants are known as particularly damaging invaders and their introductions are often accompanied by population-level behavioural and genetic changes that may contribute to their success. Anoplolepis gracilipes is an invasive ant that has just recently received increased attention due to its negative impact on native ecosystems. We examined the behaviour and population structure of A. gracilipes in Sabah, Malaysia. A total of 475 individuals from 24 colonies were genotyped with eight microsatellite markers. Intracolonial relatedness was high, ranging from 0.37 to 1 (mean +/- SD: 0.82 +/- 0.04), while intercolonial relatedness was low (0.0 +/- 0.02, range -0.5-0.76). We compared five distinct sampling regions in Sabah and Brunei. A three-level hierarchical F-analysis revealed high genetic differentiation among colonies within the same region, but low genetic differentiation within colonies or across regions. Overall levels of heterozygosity were unusually high (mean H(O) = 0.95, mean H(E) = 0.71) with two loci being entirely heterozygous, indicating an unusual reproductive system in this species. Bioassays revealed a negative correlation between relatedness and aggression, suggesting kinship as one factor facilitating supercolony formation in this species. Furthermore, we genotyped one individual per nest from Sabah (22 nests), Sarawak (one nest), Brunei (three nests) and the Philippines (two nests) using two mitochondrial DNA markers. We found six haplotypes, two of which included 82.1% of all sequences. Our study shows that the sampled area in Sabah consists of a mosaic of differently interrelated nests in different stages of colony establishment. While some of the sampled colonies may belong to large supercolonies, others are more likely to represent recently introduced or dispersed propagules that are just beginning to expand.     

The socially parasitic ant Polyergus mexicanus has host‐associated genetic population structure and related neighbouring colonies

The genetic structure of populations can be both a cause and a consequence of ecological interactions. For parasites, genetic structure may be a consequence of preferences for host species or of mating behaviour. Conversely, genetic structure can influence where conspecific interactions among parasites lay on a spectrum from cooperation to conflict. We used microsatellite loci to characterize the genetic structure of a population of the socially parasitic dulotic (aka “slave‐making”) ant (Polyergus mexicanus), which is known for its host‐specificity and conspecific aggression. First, we assessed whether the pattern of host species use by the parasite has influenced parasite population structure. We found that host species use was correlated with subpopulation structure, but this correlation was imperfect: some subpopulations used one host species nearly exclusively, while others used several. Second, we examined the viscosity of the parasite population by measuring the relatedness of pairs of neighbouring parasitic ant colonies at varying distances from each other. Although natural history observations of local dispersal by queens suggested the potential for viscosity, there was no strong correlation between relatedness and distance between colonies. However, 35% of colonies had a closely related neighbouring colony, indicating that kinship could potentially affect the nature of some interactions between colonies of this social parasite. Our findings confirm that ecological forces like host species selection can shape the genetic structure of parasite populations, and that such genetic structure has the potential to influence parasite‐parasite interactions in social parasites via inclusive fitness.     

Variation in the level of aggression, chemical and genetic distance among three supercolonies of the Argentine ant in Europe

In their invasive ranges, Argentine ant populations often form one geographically vast supercolony, genetically and chemically uniform within which there is no intraspecific aggression. Here we present regional patterns of intraspecific aggression, cuticular hydrocarbons (CHCs) and population genetics of 18 nesting sites across Corsica and the French mainland. Aggression tests confirm the presence of a third European supercolony, the Corsican supercolony, which exhibits moderate to high levels of aggression, depending on nesting sites, with the Main supercolony, and invariably high levels of aggression with the Catalonian supercolony. The chemical analyses corroborated the behavioural data, with workers of the Corsican supercolony showing moderate differences in CHCs compared to workers of the European Main supercolony and strong differences compared to workers of the Catalonian supercolony. Interestingly, there were also clear genetic differences between workers of the Catalonian supercolony and the two other supercolonies at both nuclear and mitochondrial markers, but only very weak genetic differentiation between nesting sites of the Corsican and Main supercolonies (F(ST) = 0.06). A detailed comparison of the genetic composition of supercolonies also revealed that, if one of the last two supercolonies derived from the other, it is the Main supercolony that derived from the Corsican supercolony rather than the reverse. Overall, these findings highlight the importance of conducting more qualitative and quantitative analyses of the level of aggression between supercolonies, which has to be correlated with genetic and chemical data.     

Hierarchical analysis of population genetic structure in the monogynous ant Cataglyphis cursor using microsatellite and mitochondrial DNA markers

Despite having winged queens, female dispersal in the monogynous ant Cataglyphis cursor is likely to be restricted because colonies reproduce by fission. We investigated the pattern of population genetic structure of this species using eight microsatellite markers and a mitochondrial DNA (mtDNA) sequence, in order to examine the extent of female and nuclear gene flow in two types of habitat. Sampling was carried out at a large spatial scale (16 sites from 2.5 to 120 km apart) as well as at a fine spatial scale (two 4.5-km transects, one in each habitat type). The strong spatial clustering of mtDNA observed at the fine spatial scale strongly supported a restricted effective female dispersal. In agreement, patterns of the mtDNA haplotypes observed at large and fine spatial scales suggested that new sites are colonized by nearby sites. Isolation by distance and significant nuclear genetic structure have been detected at all the spatial scales investigated. The level of local genetic differentiation for mitochondrial marker was 15 times higher than for the nuclear markers, suggesting differences in dispersal pattern between the two sexes. However, male gene flow was not sufficient to prevent significant nuclear genetic differentiation even at short distances (500 m). Isolation-by-distance patterns differed between the two habitat types, with a linear decrease of genetic similarities with distance observed only in the more continuous of the two habitats. Finally, despite these low dispersal capacities and the potential use of parthenogenesis to produce new queens, no signs of reduction of nuclear genetic diversity was detected in C. cursor populations.     

Extreme sex-biased dispersal in the communally breeding,nonmigratory Bechstein's bat (Myotis bechsteinii)

Maternity colonies of the communally breeding, nonmigratory Bechstein's bat consist of closely related females that live together with unrelated females, and average colony relatedness is low despite the absence of immigration. We compared the genetic structure of both nuclear and mitochondrial microsatellites in order to quantify sex-specific dispersal rates. More specifically, we aimed at testing whether male dispersal is able to balance the genetic effect of strong (absolute) female philopatry. Absolute female philopatry, indicated by an extreme mitochondrial DNA population differentiation of 96%, was indeed opposed by strong (possibly complete) male dispersal. Based on our knowledge of the biology of Myotis bechsteinii, we conclude that inbreeding avoidance is likely to be the crucial factor driving male dispersal in this species.     

Identification of polymorphic microsatellite loci in the queenless,ponerine ant Diacamma ceylonense

Diacamma ceylonense is a queenless, ponerine ant whose colonies are headed by a single, mated, egg‐laying worker referred to as the gamergate. Thus, new colonies are a result of dispersal by wingless gamergates. This is expected to influence patterns of colony dispersal and spatial distribution of genetic variablity. In order to facilitate the study of population genetic structure we have identified six unique, polymorphic, microsatellite loci. We have used fluorescence tagged primers to detect polymorphism at these loci.     

Fine‐scale genetic structure reflects sex‐specific dispersal strategies in a population of sociable weavers (Philetairus socius)

Dispersal is a critical driver of gene flow, with important consequences for population genetic structure, social interactions and other biological processes. Limited dispersal may result in kin‐structured populations in which kin selection may operate, but it may also increase the risk of kin competition and inbreeding. Here, we use a combination of long‐term field data and molecular genetics to examine dispersal patterns and their consequences for the population genetics of a highly social bird, the sociable weaver (Philetairus socius), which exhibits cooperation at various levels of sociality from nuclear family groups to its unique communal nests. Using 20 years of data, involving capture of 6508 birds and 3151 recaptures at 48 colonies, we found that both sexes exhibit philopatry and that any dispersal occurs over relatively short distances. Dispersal is female‐biased, with females dispersing earlier, further, and to less closely related destination colonies than males. Genotyping data from 30 colonies showed that this pattern of dispersal is reflected by fine‐scale genetic structure for both sexes, revealed by isolation by distance in terms of genetic relatedness and significant genetic variance among colonies. Both relationships were stronger among males than females. Crucially, significant relatedness extended beyond the level of the colony for both sexes. Such fine‐scale population genetic structure may have played an important role in the evolution of cooperative behaviour in this species, but it may also result in a significant inbreeding risk, against which female‐biased dispersal alone is unlikely to be an effective strategy.     

Genetic evidence for landscape effects on dispersal in the army ant Eciton burchellii

Inhibited dispersal, leading to reduced gene flow, threatens populations with inbreeding depression and local extinction. Fragmentation may be especially detrimental to social insects because inhibited gene flow has important consequences for cooperation and competition within and among colonies. Army ants have winged males and permanently wingless queens; these traits imply male‐biased dispersal. However, army ant colonies are obligately nomadic and have the potential to traverse landscapes. Eciton burchellii, the most regularly nomadic army ant, is a forest interior species: colony raiding activities are limited in the absence of forest cover. To examine whether nomadism and landscape (forest clearing and elevation) affect population genetic structure in a montane E. burchellii population, we reconstructed queen and male genotypes from 25 colonies at seven polymorphic microsatellite loci. Pairwise genetic distances among individuals were compared to pairwise geographical and resistance distances using regressions with permutations, partial Mantel tests and random forests analyses. Although there was no significant spatial genetic structure in queens or males in montane forest, dispersal may be male‐biased. We found significant isolation by landscape resistance for queens based on land cover (forest clearing), but not on elevation. Summed colony emigrations over the lifetime of the queen may contribute to gene flow in this species and forest clearing impedes these movements and subsequent gene dispersal. Further forest cover removal may increasingly inhibit Eciton burchellii colony dispersal. We recommend maintaining habitat connectivity in tropical forests to promote population persistence for this keystone species.     

Abandoning aggression but maintaining self-nonself discrimination as a first stage in ant supercolony formation

An ant supercolony is a very large entity with very many queens. Although normal colonies of small extent and few queens remain distinct, a supercolony is integrated harmoniously over a large area [1, 2]. The lack of aggression is advantageous: Aggression is costly, involving direct and indirect losses and recognition errors [3, 4]. Indeed, supercolonial ants are among the ecologically most successful organisms [5-7]. But how supercolonies arise remains mysterious [1, 2, 8]. Suggestions include that reduced within-colony relatedness or reduced self-nonself discrimination would foster supercolony formation [1, 2, 5, 7, 9-12]. However, one risks confusing correlation and causality in deducing the evolution from distinct colonies to supercolonies when observing established supercolonies. It might help to follow up observations of another lack of aggression, that between single-queened colonies in some ant species. We show that the single-queened Lasius austriacus lacks aggression between colonies and occasionally integrates workers across colonies but maintains high within-colony relatedness and self-nonself discrimination. Provided that the ecological framework permits, reduced aggression might prove adaptive for any ant colony irrespective of within-colony relatedness. Abandoning aggression while maintaining discrimination might be a first stage in supercolony formation. This adds to the emphasis of ecology as central to the evolution of cooperation in general [13].     

Coexistence of the social types: genetic population structure in the ant Formica exsecta

The ant Formica exsecta has two types of colonies that exist in sympatry but usually as separate subpopulations: colonies with simple social organization and single queens (M type) or colonial networks with multiple queens (P type). We used both nuclear (DNA microsatellites) and mitochondrial markers to study the transition between the social types, and the contribution of males and females in gene flow within and between the types. Our results showed that the social types had different spatial genetic structures. The M subpopulations formed a fairly uniform population, whereas the P subpopulations were, on average, more differentiated from each other than from the nearby M subpopulations and could have been locally established from the M-type colonies, followed by philopatric behavior and restricted emigration of females. Thus, the relationship between the two social types resembles that of source (M type) and sink (P type) populations. The comparison of mitochondrial (phiST) and nuclear (FST) differentiation indicates that the dispersal rate of males is four to five times larger than that of females both among the P-type subpopulations and between the social types. Our results suggest that evolution toward complex social organization can have an important effect on genetic population structure through changes in dispersal behavior associated with different sociogenetic organizations.     

Queen execution increases relatedness among workers of the invasive Argentine ant,Linepithema humile

Polygyny in social insects can greatly reduce within‐nest genetic relatedness. In polygynous ant species, potential rival queens in colonies with multiple queens are often executed by other queens, workers, or both. The Argentine ant, Linepithema humile, native to South America, forms a “supercolony” that is composed of a large number of nests and is considered to contribute to the ant's invasion success. Currently, four mutually antagonistic supercolonies are contiguously distributed within a small area of Japan. Here, we analyzed the genetic structure and relatedness within and among the four supercolonies using microsatellite markers to clarify how L. humile maintains its supercoloniality. The results of AMOVA and BASP, the F_ST values, and the existence of several private alleles indicated that the L. humile population in the Kobe area had a characteristic genetic structure. Within a given supercolony, there was significant genetic differentiation (F_ST) among workers collected in May and those collected in September. The significant deviation from Hardy–Weinberg equilibrium increased, and the relatedness among workers significantly increased from May to September in all supercolonies. This result suggested that the supercolonies replaced old queens with new ones during the reproductive season, thus supporting the plausibility of queen execution. From the perspective of kin selection, workers collectively eliminate queens, thereby increasing their own inclusive fitness. Restricted gene flow among supercolonies, together with mating with sib and queen execution, could help to maintain the unique social structure of L. humile, the distribution of which is expanding worldwide.