At the brink of supercoloniality: genetic, behavioral, and chemical assessments of population structure of the desert ant Cataglyphis niger

The nesting habits of ants play an important role in structuring ant populations. They vary from monodomy, a colony occupies a single nest, via polydomy, a colony occupies multiple adjacent nests, to supercoloniality, a colony spans over large territories comprising dozen to thousands nests without having any boundaries. The population structure of the desert ant Cataglyphis niger, previously considered to form supercolonies, was studied using genetic, chemical and behavioral tools in plots of 50x50 meters at two distinct populations. At the Palmahim site, the plot comprised 15 nests that according to the genetic analysis constituted three colonies. Likewise at the Rishon Leziyyon site 14 nests constituted 5 genetic colonies. In both sites, both chemical analysis and the behavioral (aggression) tests confirmed the colony genetic architecture. The behavioral tests also revealed that aggression between colonies within a population was higher than that exhibited between colonies of different populations, suggesting the occurrence of the “nasty neighbor” phenomenon. In contrast to supercolony structure previously reported in another population of this species, the presently studied populations were composed of polydomous colonies. However, both the genetic and chemical data revealed that the inter-colonial differences between sites were larger than those within site, suggesting some within-site population viscosity. Thus, C. niger exhibits flexible nesting characteristics, from polydomy to supercoloniality, and can be considered at the brink of supercoloniality. We attribute the differences in population structure among sites to the intensity of intraspecific competition.


INTRODUCTION
The diversity of spatial organization of ant populations is associated greatly with differences among species in life history, reproductive strategies and dispersion (Bourke and Franks, 1995;Crozier and Pamilo, 1996;Ross, 2001). The classical, possibly ancestral ant colony is composed of a single, monandrous (single mating) queen that is the sole reproducer and her sterile daughters that act as workers and tend to the brood and all other colony tasks. This basic colony kin structure, leading to high levels of relatedness among colony members, is believed to have evolved via kin selection (Hölldobler and Wilson, 1990;Hughes et al., 2008a;Boomsma, 2009). Typically, each colony-family occupies a nest around which it maintains a territory that it defends against intruders, con-and heterospecific alike.
Genetic studies, however, revealed that breeding systems that lower intracolony relatedness have evolved multiple times in ants (Hughes et al., 2008b). For example, both queen number and queen mating frequency may vary greatly between species, and even within species in a population. Colonies may be headed by a single queen (monogynous) or by two or more queens (polygynous) (Bourke and Franks, 1995;Crozier and Pamilo, 1996;Purcell and Chapuisat, 2013). Species may also vary in the number of males with which the queen(s) mate, from strict monandry to high levels of polyandry . This social polymorphism is the result of multiple selective pressures such as ecological constraints on nesting success, predation risks, or increased genetic diversity to better defend against pathogens or to enhance division of labor and colony task efficiency (Bourke and Franks, 1995;Keller, 1995;Crozier and Pamilo, 1996). Colony queen number is usually associated with profound changes in life history strategies and dispersal behavior (Bourke and Franks, 1995;Keller, 1995;Crozier and Pamilo, 1996;Ross, 2001;Chapuisat et al., 2004). Monogyny is generally (but see exceptions in Peeters and Molet, 2010;Cronin et al., 2013) associated with long-range nuptial flights and independent colony foundation. In contrast, under polygyny, mating occurs inside or in the vicinity of the natal nest, and colony reproduction proceeds by budding whereby the mated queen(s) leave their natal nest with a worker force to found a new colony nearby. At the population level, dispersion by budding usually results in a pattern of genetic isolation-by-distance (Sundstrom et al., 2005;Leniaud et al., 2009). In some ant species, the daughter colonies maintain exchanges of workers, brood, and food with the mother nest creating polydomous colonies, each composed of several nest units (Hölldobler and Wilson, 1990;Debout et al., 2007). In its most extreme form, polydomy can lead to the creation of supercolonies composed of several hundreds of nests, with no aggression between workers that belong to different units over a vast geographic scale (Pedersen et al., 2006;Thomas et al., 2006). Such unicoloniality is typical of many invasive ants, and possibly also the key to their success (Holway et al., 2002).
Supercoloniality poses an evolutionary challenge of understanding the selective forces that result in their prosperity, and whether or not kin selection operates despite the multiple genetic lineages that compose such supercolonies (Helantera et al., 2009;Moffett, 2012). Unicoloniality is also typified by the lack of aggression, even between very distant nests e.g., Linepithema humilis (Giraud et al., 2002;Brandt et al., 2009b;Blight et al., 2012), Wasmannia auropunctata (Errard et al., 2005;Vonshak et al., 2009;Foucaud et al., 2010) or Pheidole megacephala (Fournier et al., 2009). This is attributed to the relaxation or complete loss of nestmate recognition in these population, the mechanism of which is postulated to be either "genetic cleansing" following introduction (Giraud et al., 2002) or a "bottleneck" phenomenon (Suarez et al., 2008). Ample evidence indicate that hydrocarbons constitute at least part of the recognition system (reviewed in Hefetz, 2007), and that there is high degree of congruency in the hydrocarbon recognition cues within a supercolony and distinct compositions between supercolonies (Errard et al., 2005;Brandt et al., 2009a). Recently, it was shown in the unicolonial ant Formica yessensis that profile similarity/disparity as well as incomplete chemical discrimination via the hydrocarbon sensilla (Kidokoro-Kobayashi et al., 2012) are important factors in the establishment of unicoloniality.
Cataglyphis ants have evolved a remarkable number of life histories and reproduction strategies (Lenoir et al., 2009). Some species are strictly monogynous (Pearcy et al., 2004a;Timmermans et al., 2008) while others are facultatively or obligatory polygynous (Leniaud et al., 2011;Eyer et al., 2012). Monogynous species may be organized as independent colonies (Leniaud et al., 2011) while others are organized in polydomous colonies, with one central nest holding the queen and the rest of the nests function as queenless satellites (Dahbi et al., 1997). In contrast to the usual association between monogyny and longrange nuptial flights, young queens of the monogyne C. cursor and C. floricola mate close to their natal nest and colony reproduction proceeds by fission (Lenoir et al., 1988;Hardy et al., 2008;Amor et al., 2011;Chéron et al., 2011). Queens of some species are singly mated (Eyer et al., 2012;Leniaud et al., 2012), whereas in other species they are obligatory multiply mated (Timmermans et al., 2008(Timmermans et al., , 2010Pearcy et al., 2009). Finally, queens of several species use alternative reproductive pathways for the production of the queen and worker castes. New queens are produced asexually by thelytokous parthenogenesis, while workers are produced by normal sexual reproduction (Pearcy et al., 2004a;Leniaud et al., 2012;Eyer et al., 2013).
A recent population genetic study of Cataglyphis niger revealed that it is polygynous, that queens mate multiply and that the population sampled was supercolonial (Leniaud et al., 2011). However, this study encompassed only a single population and a limited sample size of 12 nests in an area of 4 square kilometers. Moreover, neither the chemistry or aggressive behavior of its nest constituents were studied to support the genetic data.
In the present study, we investigated the association between genetic architecture, chemical signature, and nestmate recognition in two, newly sampled populations of C. niger. First, we determined the structure of populations, using microsatellites markers. Second, we measured intra-and inter-population tolerance on the basis of aggression tests between workers. Finally, we examined the pattern of variation of worker cuticular hydrocarbon profiles in relation to population genetic structure and the spatial distribution of colonies within the population and nests within a colony.

METHODS AND MATERIALS FIELD COLLECTION AND SAMPLING
Ant samples were collected at two study sites, representing two disparate populations (∼15 km apart), on the Israeli coastal plain in Rishon Leziyyon (31.98056, 34.742863) and Palmahim (31.935046, 34.705571). All nests at each site (14 and 15 nests for Rishon Leziyyon and Palmahim, respectively) were partially excavated and approximately 50 workers sampled. The spatial distribution and distances between nests and the distances between physical nests in each plot were carefully measured on a defined and marked 50 × 50 m area. A portion of the workers was immediately frozen for further molecular and chemical analysis and the remaining workers were brought back to the laboratory for behavioral bioassays. In the laboratory they were placed in artificial nests composed of a round plastic box 10 cm in diameter and 8.5 cm in height, with a plaster laid floor to obtain sufficient moister. The ants were kept in a rearing room under a controlled temperature of 29 ± 1 • C and a photoperiod of 12L:12D, and were provided with an identical diet of sugar water and minced insects three times a week. The colonies from each site were kept in the laboratory for not longer than 1 week, during which the behavioral assays were conducted.

GENETIC ANALYSES
DNA extraction-Individual ants were ground in digestion solution (100 mM NaCl, 50 mM Tris, 1 mM EDTA, 0.5% SDS, and 200 μg/ml proteinase K) and incubated for 12 h at 55 • C in 100 μl of solution. Genomic DNA was purified through phenol/chloroform extractions and ethanol precipitation following standard protocols (Sambrook and Russell, 2001). The DNA pellet was resuspended in 100 μl of sterile H 2 O.
Genomic DNA genotyping-Genotypes of ants were determined at 9 microsatellite loci: four loci previously described for C. cursor (Cc11, Cc51, Cc54, Cc99; Pearcy et al., 2004b), three for C. hispanica (Ch01, Ch11, Ch12; Darras et al., 2014), and two new markers developed for C. niger (Cn04, Cn08). PCR were performed with a MJ Research PTC-200 thermocycler. Annealing temperature and magnesium concentration were optimized for each locus individually (available upon request). PCR products were genotyped on a ABI 3100 automated sequencer (Applied Biosystems). The size of the different alleles was determined using the Peak Scanner version 1.0 analysis software (Applied Biosystems) and used to construct a multilocus genotype for each individual. Control for genotyping errors due to null alleles and allele drop-outs was performed with Micro-Checker (Van Oosterhout et al., 2004). Possible linkage disequilibrium were tested with GENEPOP ON THE WEB (Rousset, 2008).
Genetic analyses were performed on a sample of 168 workers (X ±SE = 11.2 ± 1.6, N = 14 nests) from Rishon Leziyyon, and 168 workers (12 ± 0, N = 15 nests) from Palmahim. The number of alleles, allele frequencies, observed (H O ) and expected (H E ) heterozygosities were calculated using FSTAT 2.9.3.2 (Goudet, 1995) and averaged across loci. F-statistics were estimated according to the method of Weir and Cockerham (1984) using GENEPOP ON THE WEB (Rousset, 2008). To determine if different nests belonged to the same genetic entity (i.e., polydomy), genotypic frequencies at all collection points were compared using a log-likelihood (G) based test of differentiation, from GENEPOP ON THE WEB (Rousset, 2008). Overall significance was determined using a Fisher's combined probability test; a Bonferroni correction was applied to account for multiple comparisons. Nests were considered as belonging to different colonies if genotypic differentiation was statistically significant (α < 0.0005 and α < 0.0006 after Bonferroni correction for Palmachim and Rishon Leziyyon, respectively). Population structure was investigated by plotting ] coefficients between pairs of nests against the ln of the geographical distance (Slatkin, 1993). The correlation coefficients between genetic differentiation [Fst/(1-Fst)], ln of geographical distances and chemical distances were assessed with a Mantel test as implemented in GENEPOP ON THE WEB (Rousset, 2008). Geneticdistance based PCoA was performed with GENALEX v. 6.41. Relatedness coefficients among workers were estimated using the program FSTAT 2.9.3.2 (Goudet, 1995). The mean relatedness was estimated between workers of the same colony relative to the collection site.

BEHAVIORAL BIOASSAYS
Encounters were conducted in dyads using randomly selected workers from two randomly selected physical nests either within a site or between sites. Different workers were used for each trial, and the order of encounter types was random. Before each test, each of the workers was gently placed in an open glass tube (3 cm in diameter and 8 cm in height) within the neutral arena (10 cm in diameter and 8.5 cm in height) for 90 s for acclimatization. The test started by lifting the enclosing glass tubes and lasted for 3 min. The behavior of the ants was recorded in real time using the program JWatcher (version 0.9; http://www.jwatcher.ucla.edu), as well as videotaped, as backup, using a Logitech HD webcam pro C910 video camera. Encounters between individuals from the same nests were considered as nestmate control assays; encounters between individuals belonging to different physical nests, whether from the same or different sites, were considered as non-nestmate experimental assays. The arena was paved with a filter paper that was replaced every trial, to avoid residual odors.
Aggression was scored using a graded index from 0-4 as follows: 0, antennation; 1, mandibular threat; 2, short biting with jumping; 3, biting and /or seizing the opponent with mandibles; 4, curling the abdomen and spraying formic acid. The aggression index was calculated as follows: where Ai and ti are the aggression score and duration of each act, respectively, and T is the total interaction time defined as the sum of times in which the ants were in physical contact (not necessarily equal to assay duration). For each encounter type (that is, involving workers from two specific nests), we performed 4-6 replicate tests and the average of which was defined as the aggression index. A total of 20 (13 involving different physical nests and 6 involving nestmates) and 17 (12 involving different physical nests and 5 involving nestmates) encounter types were done at the Palmahim and Rishon Leziyyon sites, respectively (For details of the physical nests assayed see Table 1). In addition we performed 9 encounters (37 replicates in total) between nests from the two sites.

CHEMICAL ANALYSES
Previous findings have shown that the cuticular hydrocarbon (CHCs) serve as the recognition cue in C. niger (Lahav et al., 1999)  and that there is complete congruency in hydrocarbon composition between the epicuticle and those of the postpharyngeal glands (PPG) (Soroker et al., 1994(Soroker et al., , 1995. Since the PPG secretion is by far more copious than extractable cuticular hydrocarbons, it enables a much more accurate and fine tuned analysis, and can serve as representative of the CHC composition. Fifteen workers from each nest were used for the chemical analysis. The PPG of individual workers was cleanly dissected and extracted in 100 μl pentane, containing 625 ng docosane and 6250 ng of tetratriacontane as internal standards. Compound identification was done by combined GC/MS (see also Soroker et al., 1995) and was based on their fragmentation patterns, as compared to authentic standards. Sample quantification was done by gas chromatography, using a VF-5ms capillary column that was temperature-programmed from 150 • C (with 1 min initial hold) to 300 • C at a rate of 5 • C per min, with a final hold of 15 min. Peak integration was performed using the program Galaxie Varian 1.9. Compounds that constituted less than 1% across all samples were omitted from the analyses, because their quantification is unreliable. For assessing hydrocarbon profilespecificity we used multivariate statistics, i.e., discriminant analysis using the stepwise forward mode. This module of the program works initially as a principle component analysis for reducing the number of variables to fit the number of cases, and then perform the discriminant function.

RESULTS
Cartography of the nests revealed the occurrence of 14 and 15 nests, with a mean ±SD distance between nests of 7.52 ± 6.22 m and 12.56 ± 10.83 m in Palmahim and Rishon Leziyyon, respectively. In both sites, nests were not evenly distributed within the plot, but were clustered (Figures 1A,B). The index of aggregation R for the Plamahim site was 0.701 and for Rishon Leziyyon was 0.854. However, if we Exclude nest # 291 who was the only one that belonged to the genetic colony C5, the R value at Rishon Leziyyon was reduced to 0.506. Both indicate that nests were not randomly dispersed, but are aggregated (Clark and Evans, 1954).

FIGURE 1 | Spatial distribution of nests at the Palmahim (A) and Rishon Leziyyon (B) sites.
Each nest is represented by a number (e.g., 285; 286 etc.) and each genetic colony is represented by a different symbol, color and number (e.g., C1; C2 etc.).

BEHAVIORAL ASSAYS
Effect of the workers' origin (nest, colony origin, site) on aggression levels was tested using dyadic encounters. The initial selection of nests for the dyadic encounters was done randomly, without knowledge of their genetic compositions. However, since we hypothesize that nests that are genetically indistinguishable show also high degree of tolerance, the analyses of the aggression tests proceeded according to colony classification using the molecular data as follows: nestmates (i.e., workers from the same physical nest); workers of different nests but within the same genetic colony; workers from nests belonging to different genetic colonies within a site; and workers from nests of the two sites. Figure 2A presents the results according to colony and Figure 2B present the pooled results as well as the result of the encounters between colonies of the different populations. The results of both types of analyses were in accordance with colony assignment using the molecular results. Both at the Palmahim and Rishon Leziyyon sites, level of aggression between workers from different nests belonging to the same genetic colony was low and not significantly different from that between nestmates. In contrast, the aggression between workers from nests belonging to different genetic colonies, within a site, was significantly higher than that between workers of different nests within a genetic colony or nestmates (for particular encounters, Figure 2A: Palmahim site One Way ANOVA, F = 19.2, df = 6, p < 0.001, followed by Tukey post hoc test p < 0.05; Rishon Leziyyon site (One Way ANOVA, F = 20.26, df = 9, p < 0.001, followed by Tukey post hoc test p < 0.05) (For pooled encounters, Figure 2B: One Way ANOVA, F = 54.75, df = 6, p < 0.01, followed by Tukey post hoc test p < 0.05). Interestingly, the aggression between workers belonging to colonies from the different sites was significantly lower than that between colonies of the same site.

CHEMICAL ANALYSES
The chemical composition and compound identification of the PPG secretion is presented in Figure 3. The secretion is composed of hydrocarbons ranging from tetracosane to tritriacontane comprising linear, monomethyl, diemethyl, and trimethyl branched components. The secretion was dominated by 2 peaks comprising a mixture of 11-+ 13-+ 15-methylnoacosane and 11,15-+ 13,15-dimethylnonacosane, similar to previous analyses (Soroker et al., 1995).
To examine whether differences in hydrocarbon profiles of workers match the behavioral and molecular data, we performed a disriminant analysis of the genetically defined colonies in Palmahim and Rishon Leziyyon (Figures 4A,B). For that purpose we used all the individuals that belong to a genetic colony irrespective of their nest origin. The analysis performed used the stepwise forwards module, which ranks the variables according to their impact and accordingly reduce their number to fit the appropriate number of cases. All three genetic colonies in Palmahim were completely chemically distinct (Figure 4A; Wilks' Lambda: 0.0090835, F (38, 96) = 23.98068 p < 0.005). Likewise, all five genetic colonies in Rishon Leziyyon were chemically distinct ( Figure 4B; Wilks' Lambda: 0.0000376 F (100, 161) = 19.44184 p < 0.0000). Table 2 compares the Mahalanobis distances between nests, measured irrespective of their grouping, within a genetic colony to that between nests belonging to different genetic colonies. For both populations, the cuticular hydrocarbon profiles of nests within a genetic colony were significantly more similar to each other than to those of nests belonging to a different genetic colony (t-test Palmahim: n 1 = 30, n 2 = 41, t = −7.5779, df = 68, p < 0.0001; t-test Rishon Leziyyon: n 1 = 16, n 2 = 61, t = −6.8134, df = 35, p < 0.0001).
For details of the physical nest assayed see Table 1

DISCUSSION
Polydomy is a nesting habit whereby a colony is composed of several nests that extend over a large territory and is characterized by constant worker movement between nests. It enables the colony a larger foraging area, the extent of which (i.e., the number of nests per colony and their spatial distribution) is limited by competition with other colonies in the population (reviewed by Hölldobler and Wilson, 1990;Debout et al., 2007). Highly extended polydomy can comprise many nests and when it extends over very large geographic areas the species is considered supercolonial. Supercolonial population-structure typifies many invasive ant populations (Holway et al., 2002), but it may also occurs in native populations (Higashi and Yamauchi, 1979;Holzer et al., 2006;Buczkowski and Bennett, 2008;Buczkowski, 2010). Polydomy and supercoloniality can be detected genetically as well as behaviorally and chemically. Members of a supercolony are tolerant to each other even if they originate from different nests at the time of their testing, and exhibit chemical uniformity of cuticular hydrocarbon compositions (Errard et al., 2005;Brandt et al., 2009b;Vonshak et al., 2009;Blight et al., 2012;Kidokoro-Kobayashi et al., 2012). Therefore, assessing the population structure of an ant species is best done by encompassing all three parameters, genetic, behavior and recognition chemistry. Recently we have reported a genetic analysis of colonies of C. niger from a population (Ashqelon), located in the southern coastal plain of Israel. We found a genetic uniformity of 12 nests that span over 4 km 2 , suggesting that the species forms a supercolony (Leniaud et al., 2011). However, this study was based on a single population with limited sample size.
In contrast to our previous finding with the Ashqelon population (Leniaud et al., 2011), our present study combining genetic, behavioral assays and chemical analyses shows that both the Palmahim and Rishon Leziyyon populations are polydomous rather than supercolonial. Our genetic analysis revealed that average relatedness values among workers from the same colony are very similar in both Palmahim and Rishon Leziyyon (r w-w ≈ 0.30). Quite surprisingly, these values are 10 folds of that registered for the Ashqelon population (0.02). Despite the fact that both sites in the present study were somewhat smaller (50 × 50 m each) than that in Ashqelon we could detect 3 genetically different colonies in Palmahim, comprising 4,−7 nests, and 5 colonies in Rishon Leziyyon, comprising 1-5 nests. Since we limited our analyses to a 50 × 50 m plots it is hard to generalize on a mean number of nests per colony. Yet, together these results do not support supercoloniality in the two presently studied populations of C. niger. The marked difference in within-colony genetic relatedness between this work and the population of Ashqelon remains enigmatic. A research constraint we had in the current study was our inability to completely excavate the nests. As a consequence, no queen could be sampled; hence, the number of queens per nest/colony and their mating frequency, two primary factors on relatedness values within colonies, remain unknown. Therefore, our genetic data are limited to the analyses of workers only. Nevertheless comparing the within and between colony genetic similarity of the two populations indicate some genetic viscosity within each, facilitating high polydomy and perhaps bringing the population to the verge of supercoloniality. Thus the differences in population structure between Ashqelon and Palmahim and Rishon Leziyyon lend credence to the hypothesis that C. niger shows large plasticity in social structure, and what we may observe is a transition state from high polydomy to supercoloniality. Similar plasticity was also recorded for Formica truncorum where in certain sites the population is largely polydomus while in other sites it is unicolonial (Rosengren et al., 1985;Elias et al., 2005).
There is a good congruency between the genetic data and both the behavior and the chemistry. Workers from different nests but from the same genetic colony were as tolerant to each other as are nestmates. This was true for both populations and is depicted more pronouncedly when we compared nestmate vs colony-mate vs. between colonies, irrespective of colony identity ( Figure 2B). Interestingly, aggression between colonies of the two populations was significantly lower than that between colonies within population. We consider this phenomenon as "nasty neighbor" (as opposed to the "dear enemy" phenomenon). Nasty neighbor phenomenon is common in social insects, in particular ants where the threat from neighboring nests is far greater than a wondering individual from a colony that is further away (Knaden and Wehner, 2003;Sanada-Morimura et al., 2003;Boulay et al., 2007;Thomas et al., 2007;Wilgenburg, 2007;Newey et al., 2010;Roux et al., 2013). Competition between colonies for space in Significance was tested using t-test for unpaired data with unequal variance. polydomous species is predicted to be even higher compared to monodomous species due to their ability to create new satellite nests as strongholds in the new territory and by that facilitating the expansion of their foraging area.
Nestmate recognition in C. niger, and in all likelihood in other ant species, is mediated by cuticular hydrocarbons (Lahav et al., 1999). We therefore estimated their chemical similarity/disparity according to the genetic colonies. For both sites there was a clear segregation of profiles according to colony identity, with almost no overlap. Comparing the chemical (Mahalanobis) distances confirmed that the distance between workers of a colony, irrespective of their nest origin, was significantly lower than that between workers of different colonies. Given the moderate genetic viscosity within site, we were interested in comparing the sites also with respect to their CHC profiles. It is clear from Figure 5A that the two populations are markedly different, much more than that based on genetics. This suggests that both genetics and the environment, i.e., difference in locality, may affect the composition of the CHCs.
In summary, all three parameters investigated, genetics, behavior and chemistry point to the fact that colonies of C. niger in both Palmahim and Rishon Leziyyon are polydomous in nature. This is in contrast to a previous study in which the population of Ashqelon was indicated by genetics to be unicolonial. We tend to exclude the possibility that in Ashqelon the population is also polydomous rather than unicolonial. The total area that the 12 colonies occupied was 4 km 2 and the mean distance between physical nests was over 50 m, in line with the phenomenon of unicoloniality. The present site span over a smaller area of 2.5 km 2 , and the mean distance between nests ranged between 7 and 12 m, in line with polydomy. These differences may be attributed to several factors, among which we consider intraspecific competition to be a major one. Adephic and climatic conditions among the 3 sites do not differ by much, but the populations in Palmahim and Rishon Leziyyon are situated sites that are highly visited by human, compared to the site in Ashqelon, which seemingly generate greater resources. If the existence and survival of a colony is a balance between foraging and fighting efforts, these conditions may be supportive for the establishment of new colonies and thus mutually limiting the expansion of each. If, on the other hand resources are more limited the establishment of new colonies become more difficult and may be prohibited by the already existing colonies. Since all population show isolation by distance, indicating reproduction by colony splitting, supercolonies may be formed in uncompetitive areas, but limited to polydomy in highly competitive area. The ant Cataglyphis niger thus merits the statement of "being at the brink of supercoloniality."

AUTHORS AND CONTRIBUTORS
Maya Saar conducted the field studies, including the behavioral assays and the chemical analyses. Laurianne Leniaud conducted the molecular population genetic studies. Abraham Hefetz and Serge Aron conceived the idea and supervised the work. All authors contributed to the writing of the manuscript.