- Research article
- Open Access
Phylogeography of social polymorphism in a boreo-montane ant
© The Author(s). 2016
- Received: 18 February 2016
- Accepted: 13 June 2016
- Published: 23 June 2016
The disjunct distribution of several Palearctic species has been widely shaped by the changes in climatic conditions during the Quaternary. The observed genetic differentiation or reproductive isolation between extant populations may be the outcome of their contemporary geographic separation or reproductive incompatibility due to differences in phenotypic traits which have evolved in isolated refugia. In the boreal ant Leptothorax acervorum, colonies from central and peripheral populations differ in social structure: colonies from Central and Northern Europe may contain several equally reproductive queens (facultative polygyny), while in colonies from peripheral populations in Spain only one the most dominant of several queens lays eggs (functional monogyny). By reconstructing the specie’s evolutionary and demographic history in Southwestern Europe we examine whether variation in social organization is associated with restricted gene flow between the two social forms.
We show that multi-queen colonies from all so far known inner Iberian populations of L. acervorum are functionally monogynous, whereas multi-queen colonies from all Pyrenean populations are polygynous, like those from other previously studied areas in Europe. Our analyses revealed complex spatial-genetic structure, but no association between spatial-genetic structure and social organization in SW-Europe. The population in the western Pyrenees diverged most strongly from other Iberian populations. Moreover, microsatellite data suggest the occurrence of recent bottlenecks in Pyrenean and inner Iberian populations.
Our study shows a lack of reproductive isolation between the two social forms in SW-Europe. This in turn suggests that demographic and spatial patterns in genetic variation as well as the distribution of social phenotypes are better explained by co-variation with climatic, ecological, and historical factors. Moreover, we for the first time show the existence of substantial spatial-genetic structure in L. acervorum, suggesting the existence of multiple refugia in SW-Europe, including two extra-Mediterranean refugia in France. While gene flow among inner Iberian refugia may have been larger during the late glacial, extra-Mediterranean refugia in southern France may have contributed to the post-glacial recolonization of W-Europe.
- Extra-Mediterranean refugia
- Gene flow
- Leptothorax acervorum
- Plastic behaviour
- Population structure
- Quaternary climate change
- Reproductive skew
Quaternary climate changes and their dramatic effects on demography and distributional ranges have left behind their footprints on the evolution of Palearctic biota. Temperate species escaped the last glacial maximum (LGM 23-18 ka BP) by retreating southwards to suitable habitats in Iberia, Italy and Balkans and recolonized Northern and Central Europe from these refugia after the retreat of the ice sheet [1, 2]. In contrast, boreal species could survive in periglacial areas in Central Europe and expanded their ranges into more southern areas. During postglacial warming, southern populations of boreal species often became “trapped” in mountainous regions, leading to their present arctic-alpine or boreo-montane distribution [3, 4].
Relict populations of boreal species in Mediterranean refugia may be exposed to other environmental conditions than populations from the centre of their range and thus accumulate unique genotypes and specific adaptations [5, 6]. This is particularly evident in the boreo-montane ant Leptothorax acervorum (Fabricius, 1793) (Hymenoptera, Formicidae). Across its range, a large fraction of L. acervorum colonies contain several reproductive queens (the others are queenless or have a single queen). Interestingly, the partitioning of reproduction among nestmate queens in such multi-queen colonies (“reproductive skew”) varies among populations [7, 8]. In colonies from boreal and temperate Eurasia, all queens contribute equally to the brood (low skew, “facultative polygyny”, ). In contrast, queens in colonies from the range margin (e.g. Central Spain and Hokkaido, Japan) form social and reproductive hierarchies and only the top-ranking queen reproduces (high skew, “functional monogyny”, [8, 10, 11]). In addition, mating behaviour differs between populations from Central Spain and Central Europe .
From these pronounced behavioural differences we expected to find genetic differentiation with reproductive isolation and incipient speciation between the two social forms. Alternatively, a lack of congruence between spatial patterns of genetic structure and social behaviour would indicate environment-induced phenotypic plasticity underlying geographic variation in behaviour (e.g. [12–15]).
The present distribution of L. acervorum reaches as far north as the tundra-taiga ecotone at the North Cape and the Lena delta [16, 17], suggesting that it was capable of surviving the ice ages in refugia in Central Europe. This should have left conspicuous traces on the genetic structure and diversity of extant populations, whose reconstruction should result in a better understanding of the respective contributions of potential extra-Mediterranean and “classical” southern-European refugia towards the current distribution of genetic diversity of W-Palaearctic populations of L. acervorum.
The aims of this study, therefore, were (i) to map the distribution of the two social phenotypes of L. acervorum in inner Iberian mountains and Pyrenees, (ii) to analyse the micro-evolutionary history of the two social forms using genetic markers, and (iii) to infer the population history of L. acervorum in SW-Europe in a broader phylogeographic context.
Ant sampling, DNA extraction and queen ovary dissection
Sampling information for the study area of Leptothorax acervorum in SW-Europe, Germany and England
Location (N, E)
No. of samples (STR/mtDNA)
15 / 14
15 / 3
Cantabrian Mts. I
15 / 4
Cantabrian Mts. II
15 / 3
Sierra de Gúdar
15 / 3
Sierra de Albarracin
2008, 2009, 2010
15 / 3
Sierra de la Demanda
15 / 3
Massif Central I
* / 2
Massif Central II
* / 6
Massif Central III
* / 3
* / 5
* / 3
* / 8
* / 2
One worker each from 105 colonies from the Iberian Peninsula and the Pyrenees was genotyped at ten polymorphic microsatellite loci. The mitochondrial DNA dataset contains sequences of 62 workers representing all sampled locations and geographic regions. In addition, we used DNA-samples from England (E, see  and Table 1). All specimens were preserved in 95 % ethanol or frozen at −20 °C until DNA extraction. Genomic DNA was extracted from whole workers using a cetylmethyl ammonium bromide (CTAB)-protocol .
Additionally, for the determination of social structure and reproductive skew we examined how many of the queens co-occurring in single multi-queen nests (a colony) had been reproductively active by dissecting their ovaries (528 queens from 66 colonies from all over Pyrenees (PY) and four inner Iberian populations (SG, SD, SNW I, SNW II). From the 12 Pyrenean colonies one was collected from PY I and four from PY II. The remaining seven colonies were collected from three additional Pyrenean locations (no. of colonies per location; 1, 2, and 4, respectively). Previous analysis had already confirmed functional monogyny for the SA population [8, 20]. Ovary dissections were conducted right after transfer of colonies to the laboratory (August 2009: SG; November 2009: PY, SD; February 2010: SD; July 2010: SNW I, SNW II, SD). For further details on ovary dissections and determination of reproductive skew see Additional file 1.
Four of the ten used microsatellite loci were developed specifically for L. acervorum (GA1, GA2, GT1, GT2; . The remaining loci represent cross-amplifications from other ants (GT218, GT223; , 2MS67; , 2MS46; Suefuji et al., unpublished – 2MS46fwd: 5′- GCTCACTACTATGCTGCCAGC -3′, 2MS46rev: 5′- TCCGTCTATCCCTTCCTGCAA -3′, L18;  and Myrt3; . PCR conditions were chosen as follows: initial denaturation 5 min at 95 °C; 35 cycles of 60 s at 95 °C, 45 s at the locus-specific annealing temperature of 45–60 °C, elongation of 45 s at 72 °C and a locus-specific final extension step of 30–120 s at 72 °C. Total reaction volume was 20 μL with 1 μL DNA template. PCR products were either analysed on an ABI PRISM 310 automated sequencer (GA1, GA2, GT218, GT223) and subsequently genotyped using GENESCAN 3.1 (Applied Biosystems) or sent to GATC Biotech AG (GT1, GT2, L18, Myrt3a, 2MS67, 2MS46 (II)) and subsequently genotyped using Peak Scanner Software v1.0 (Applied Biosystems).
Population structure, gene flow and demographic history
Number of alleles (k), allelic richness (A), number of private alleles (A P ), and expected (HE) and observed (HO) heterozygosities were calculated per population and locus using FSTAT 184.108.40.206  and GENALEX 6.5 . The genealogical relationships among study populations were analysed by neighbour-joining trees in POPULATIONS 1.2.31  using D A distance . Bootstrap values were obtained by 2000 replications over loci. To study population structure in more detail we estimated overall and pairwise F-statistics using GENEPOP as well as GENALEX to calculate allelic-diversity corrected F ST – value analogues (Hedrick’s standardized G ST , G ST corrected for small sample sizes G” ST , and Jost’s D, ). To test the hypothesis that genetic differentiation is equal or higher within than between regions (defined here as geographic regions – Iberian Peninsula (IB) and Pyrenees (PY) or high skew vs. low skew populations), we conducted an analysis of molecular variance (AMOVA) using ARLEQUIN v220.127.116.11 . Significance of results was evaluated over 10,000 replicates.
We applied Bayesian clustering to assess fine-scale population genetic structure (as well as identifying distinct populations) based upon multi-locus genotypic data. Cluster analysis was run in STRUCTURE 2.3.4  allowing individuals to have mixed ancestry (admixture model) without using sampling locations as prior information and without correlated allele frequencies among populations. As recommended by (, and cf. STRUCTURE documentation), we run a first analysis including microsatellite data from all locations and in a second analysis excluded samples from the most divergent population (PY II). Potential population cluster values (K) varied from 1 to 10 with ten runs per value of K, burn-in and sampling period were set to 300,000 (first analysis) and to 200,000 generations (second analysis) accordingly. For each analysis the optimal K value was assessed by following the ΔK-method as described by  and implemented in STRUCTURE HARVESTER web-v0.6.93 . DISTRUCT v1.1  was used to graphically display the output. Additionally, we used a Mantel test  to evaluate genetic isolation by distance in the microsatellite data, with significance of results evaluated by 999 matrix permutations in GENALEX. Genetic distance matrix among populations was calculated as linearized F ST-values [(F ST (1- F ST)−1] in FSTAT whereas geographical distances among populations in kilometres were calculated in GENALEX (based upon the coordinates shown in Table 1).
Because of the patchy distribution of L. acervorum in “mountainous islands” in Central Spain, we tested for the occurrence of bottlenecks in each sampled population. First, we used Wilcoxon’s sign rank test, which tests for an excess of heterozygosity, implemented in BOTTLENECK 1.2.02 . Second, we calculated the M-ratio statistic  to test for a severe reduction in effective population size. Finally, we used the coalescent-based Bayesian MCMC algorithm as implemented in MIGRATE-N v 3.6.10  for the inference of demographic parameters, the analyses of demographic changes over time and estimation of pairwise migration rates. For details on all methods (including parameter settings) see Additional file 1.
Mitochondrial DNA analyses
The primers C1-J-2183 and A8-N-3914  were used to amplify an 1641 bp long mitochondrial DNA fragment, starting from within the COI gene, including the complete COII sequence, and finishing in the very beginning of the ATPase 8 gene. PCR was carried out in a total reaction volume of 15 μL using the BIO-X-ACT Short Mix (Bioline) and 1 μL DNA template. PCR conditions consisted of an initial denaturation 4 min at 94 °C; 38 cycles of 75 s at 94 °C, 75 s at 50 °C (annealing), elongation of 150 s at 72 °C, and a final extension step of 5 min at 72 °C. PCR products were sent to LGC Genomics for purification and Sanger sequencing.
Chromatograms were assessed and edited using CHROMAS LITE 2.1.1 (Technelysium) and subsequently concatenated in BIOEDIT . Sequences were aligned manually and automatically by using the algorithm CLUSTAL W as implemented in BIOEDIT. All sequences could be aligned unambiguously, and no indels were found except in three sequences from the eastern Pyrenees (PY I) that contained a 1 bp-deletion. Since this appeared in a noncoding region only, it was considered valid and was used as a fifth mutational state in network analysis. The absence of unusual stop codons as potential evidence for the presence of pseudogenes (numts, e.g. ) was checked in the ARTEMIS GENOME BROWSER using the invertebrate mitochondrial codon table .
Phylogeography and population structure
For the quantification of genetic polymorphism we used the following standard diversity indices: number of haplotypes (h), number of private haplotypes (h P ), number of segregating sites (S), haplotype diversity (H) and nucleotide diversity (π) for each locality, for each geographic region, and for the whole mtDNA dataset except the extremely divergent haplotypes using DNASP v5.10.01 . We estimated genetic differentiation with and without specifying groups (i.e. geographic regions as defined above), following  by calculating Φ ST – values from mean pairwise differences (N ST ) and haplotype frequencies (G ST ) in ARLEQUIN. By comparing both statistics of genetic differentiation it is possible to test whether the mtDNA-sequence dataset contains a signal of phylogeographic structure (if N ST > > G ST ) beyond that in haplotype frequencies alone (see  for detailed discussion). AMOVA with pairwise comparisons between regions was carried out to test the hypothesis that genetic differentiation is equal or higher within regions (Φ SC ) than between regions (Φ CT , regions as defined in the microsatellite data). When both measures are equal, their ratio is expected to be one. Consequently, if their ratio is larger than one, then genetic differentiation is larger within than between regions (which in case of the PY – IB comparison would also mean: no stronger genetic differentiation between social forms). Significance of AMOVA results was evaluated over 10,000 replicates.
Mismatch distributions and test statistics for selective neutrality (Tajima’s D and Fu’s F S ) were calculated for each region and for the whole dataset to measure whether data deviate from expectations of neutral evolution or equilibrium population size. For further details on methods see Additional file 1.
Ovary dissections and reproductive skew
Reproductive status and skew of Pyrenean (PY) and inner Iberian colonies of L. acervorum
Genetic diversity indices
Genetic diversity indices for 10 microsatellite loci in Iberian populations of L. acervorum
Genetic diversity indices and neutrality tests for mtDNA sequences of L. acervorum in SW- and C-Europe
Fu’s F S
Inner Iberian (IB)
Southern France (FR)
PY + IB + FR
Population structure and phylogeography
Analysis of genealogical relationships among populations (based upon microsatellite data) resulted in an unrooted neighbour-joining tree with nodes supported by bootstrap values between 10 and 96 % (Fig. 2a). The longest internal branch separates the two Pyrenean populations from the remaining Iberian populations. The population from the western Pyrenees (PY II) diverges most among all analysed populations.
Total genetic differentiation (F ST ) was low (0.075) and remained moderate even after correction for multiple alleles (G’ ST : 0.257, G” ST : 0.265, D est: 0.205; all values significant with P < 0.01). The population from the western Pyrenees (PY II) again showed a higher differentiation than all inner Iberian populations (see Additional file 2).
STRUCTURE runs using all locations strongly supported K = 2 genetic clusters (ΔK 2 : 250.6, ΔK 3 : 8.6; Additional file 2). With a probability of ≥80 %, 92 % of the individuals from the inner Iberian populations could be assigned to the first cluster, and 95 % of the individuals from the western Pyrenees (PY II) could be assigned to the second cluster (Fig. 2b). In contrast, genotypes from the eastern Pyrenees (PY I) suggested pronounced admixture (Fig. 2b). An analysis without the western Pyrenean (PY II) samples supported four genetic clusters (K = 4) within the remaining microsatellite dataset (ΔK 4 : 5.9, next smaller ΔK 2 : 5.0, see Additional file 2). Overall, the level of admixture was higher among the remaining populations (Fig. 2c). In particular, PY I, SNW I, and SG showed evidence of high admixture. However, 70 % of the individuals from the remaining Iberian System (SA and SD) could be assigned to a third cluster with a probability ≥80 %, while 67 % of the samples from SNW II could be assigned to a fourth cluster. Tests for isolation by distance (IBD) patterns in the microsatellite data were non-significant for Euclidean (r xy = 0.010, P = 0.5) and log-transformed Euclidean distance (r xy = 0.083, P = 0.4).
The mtDNA-haplotype network has a decentralized structure and shows a short to intermediate expanded genealogy with only limited geographic segregation of haplotypes (i.e. phylogeographic structure, Fig. 3). The only exceptions are two extremely divergent haplotypes from PY I and a nearby population (FR III), which differed from the core network in 31 (23 + 8) and 26 mutations, respectively. In total, the network comprises 28 haplotypes (22 haplotypes without extremely divergent, English, and German haplotypes), of which three are dominant. The first is predominantly distributed in southern France and western Pyrenees (7 + 1 sequences, among them two sequences from FR III), the second from the Pyrenees (6 + 1 sequences from PY I and PY II, respectively) and the third predominantly in inner Iberia (five sequences in total). The remaining French samples formed a second predominantly French haplogroup (Fig. 3).
Results of analysis of molecular variance (AMOVA) by pairwise comparison between regions for L. acervorum
Percentage of variation
IB – PY
G CT (Φ CT )
G SC (Φ SC )
G ST (Φ ST )
IB – PY
PY – FR
FR – IB
IB – PY – FR
Analyses of historical demographic events using summary statistical tests of mtDNA variation provided no significant evidence for demographic changes or selective events within each geographic region (Table 4).
Spatial distribution: behavioural polymorphism and genetic structure
In this study, we provide data on the distribution of the two social forms (facultative polygyny vs. functional monogyny) of the ant Leptothorax acervorum in the Iberian Peninsula and Pyrenees. We show that colonies are functionally monogynous in all studied inner Iberian mountain ranges, as previously shown for the populations of Sierra de Albarracin (SA) and Sierra de Gúdar (SG) [8, 11, 20]. In contrast, colonies from the Pyrenees are facultatively polygynous, like those from other areas in Central and Northern Europe . This supports previous studies according to which high skew may be an adaptation to difficult queen dispersal in patchy habitats, as commonly found at the margin of a species’ range (see [10, 50] for further discussion).
On a first glance, our population genetic data seem to support an association of genetic differentiation and social organization in L. acervorum (e.g. Fig. 2a). However, a more detailed analysis shows stronger genetic differentiation within than between social forms (Fig. 4 & Table 5). In particular, western Pyrenees (PY II) revealed the highest pairwise differentiation to all inner Iberian populations while individuals from eastern Pyrenees (PY I) revealed pronounced admixture from both genetic clusters. Furthermore, haplotype network analysis (mtDNA) failed to find a clear pattern of differentiation between the two social forms.
These results suggest a lack of reproductive isolation between both social forms, contrasting with the stability of social phenotype in the laboratory . Studies with 11 genetic markers might miss genes underlying variation in social structure. For instance, in two ant species queen number is controlled by non-recombining “social chromosomes” [51, 52]. Nevertheless, wherever a genetic basis for social polymorphism has been demonstrated (i.e. via genetic differentiation, phenotypic stability, or genomic rearrangements), different social phenotypes are also associated with mtDNA differentiation (for a literature review see Additional file 2). The lack of mtDNA divergence in L. acervorum might therefore speak against a robust genetic basis of the social polymorphism. Furthermore, experimental manipulation of colonies from a low skew population elicited queen antagonism  and high skew , suggesting that queens can adapt their behaviour to environmental changes. It remains unknown whether queens from functionally monogynous populations behave similarly flexibly, have lost their behavioural plasticity, or have a changed threshold for queen fighting. Finally, phylogenetic analysis showed that functional monogyny evolved convergently in several lineages of the genus Leptothorax and thus corroborates skew as a flexible trait that can evolve with environmental changes .
Phylogeography within Iberian Peninsula
Patterns of genetic diversity and genetic structure in L. acervorum in inner Iberian mountains match recent studies that revealed substantial geographic complexity in genetic diversity, divergence, and extent of admixture among and within Iberian species (“refugia within refugia hypothesis”; see [55, 56] for reviews). The fragmented Iberian landscape with many mountain ranges and river basins allowed the survival of populations in multiple refugia and facilitated their partial differentiation in allopatry, followed by secondary contact between previously differentiated genomes . L. acervorum is a cold-adapted ant species [16, 17] and in the Iberian Peninsula lives mostly in mountainous pinewoods or pine-dominated forests above 1500 m a. s. l. Populations may have co-expanded and -contracted with Pinus sylvestris-dominated forests during glacial and interglacial periods. This suggests larger population sizes and higher levels of gene flow among inner Iberian populations during glacial and transitional phases (e.g. early Holocene). In contrast, during the last interglacial populations would have fragmented and declined.
Our molecular analysis supports such a scenario: on a local scale, microsatellites suggest no current gene flow between inner Iberian populations. In addition, the patchy distributions of private alleles (A P ) and the occurrence of populations with high vs. low admixture (PY I, SG, SNW I vs. PY II, SNW II, SA, SD) are in accordance with the existence of multiple refugia and varying connectivity among them throughout the late Pleistocene. In contrast, mtDNA data did not reveal a pronounced genetic differentiation among Iberian populations.
The westernmost Cantabrian population (SNW II) provides a notable exception to this pattern. It shows a distinct mtDNA haplotype and forms a distinct genetic cluster in the multi-locus microsatellite data, suggesting a longer period of evolution in isolation. This resembles the situation in the capercaillie (Tetrao urogallus), a bird characteristic of Eurasian coniferous forests [57, 58]. Its Cantabrian population inhabits deciduous forests and is the most genetically distinct and depauperate in Europe except for the Pyrenean population, to which it is closely related [59, 60].
Because of the fragmented distribution of L. acervorum in the Iberian mountains, we tested for the occurrence of bottlenecks. The mtDNA demographic analyses (summary statistical tests, mismatch distributions, raggedness) did not provide clear evidence for historical changes in effective population sizes. In contrast, M-ratios in microsatellites were far lower than the critical M value of 0.68 (as derived by ) in all populations, including the Pyrenees. When compared to simulated critical M values, all seven populations seem to have experienced bottlenecks at values of θ (i.e. mutation rate scaled effective population size) that are reasonably large for L. acervorum (θ = 20–40, corresponding to pre-bottleneck N e of ~50,000–100,000). In addition, MIGRATE-N analysis confirmed the occurrence of bottlenecks in all seven populations, while excess heterozygosity confirmed recent bottlenecks in only two populations from the Iberian system (SA and SG). The different results between tests can be explained by differences in test procedures and the low power of the heterozygosity excess test due to low sample size or deviations in mutation models (for a detailed discussion see [40, 61] and references therein). We therefore conclude that by the above mentioned properties the heterozygosity excess test is more likely underestimating the number of bottlenecked populations compared to the other used bottleneck test procedures (in particular when bottlenecks did not occur recently). Furthermore, the consistently lower observed than expected heterozygosities (i.e. lower frequency of heterozygote genotypes than expected under Hardy-Weinberg equilibrium) are in agreement with the occurrence of bottlenecks in all inner Iberian and Pyrenean populations. Finally, the occurrence of bottlenecks in the recent past is in accordance with postglacial range contractions or shifts and the current patchy distribution (above 1500 m. s. l.) of L. acervorum in the Iberian Peninsula, as indicated by both genetic datasets.
Relationships of Iberian populations to the adjacent regions (Pyrenees, W-Alps and Central Europe)
In contrast to the clear signal of population structure (microsatellites) among Iberian populations showed here for the first time, mtDNA sequence data did not reveal any spatial-genetic structure. This matches former studies on L. acervorum across Europe [18, 20, 62]. The data do not support the existence of a single SW-European refugium but two glacial refugia in France alone. One haplotype, found exclusively in the western Pyrenees (PY II), southernmost Massif Central (FR III) and Mont Ventoux (FR IV), suggests a refugium extending from the lower Rhône valley westward to the foothills and lowlands north of the Pyrenees. A second refugium in the upper Rhône valley or adjacent areas in the north is indicated by a haplogroup found in the remaining French localities (FR I, II, V). Glacial refugia north of the Mediterranean peninsulas have recently been recognized for several species [4, 63–67]. In particular, artic-alpine species show patterns of genetic association similar to L. acervorum in southern and Central France (e.g. ).
Subsequent to the late glacial (post-10 ka BP), the first French refugium may have started to fragment and L. acervorum shifted its range towards more suitable habitats in the surrounding mountains and their foothills, while at the same time populations of the second refugium would have expanded their range at least towards Massif Central and western Alps. A closely related haplotype from England suggests a direct connection to the first French refugium (including two sequences that only differ by one mutation from the main haplotype of the latter, Fig. 3). Analogous scenarios of postglacial range dynamics have been described for other European mountain species (e.g. ). For example, patterns of genetic structure within and between W-European high mountains of the mountain butterfly Erebia epiphron suggest (i) recolonization of the central and western Pyrenees and western Alps from refugial lowlands in between and (ii) of the eastern Pyrenees from the unglaciated foothills south-east and east of it .
Similar to the latter scenario, our analysis revealed a differential contribution of the first French refugium and the inner Iberian refugial system to the postglacial recolonization of Pyrenees not just in colonization routes but also in patterns of gene flow by both genetic markers. In contrast to the French refugium, inner Iberian populations share no haplotypes with the Pyrenees (in particular the eastern Pyrenees), while microsatellite data indicate historical gene flow between inner Iberia and eastern Pyrenees, but no or only very restricted gene flow between the western Pyrenees and inner Iberian populations. These patterns suggest that during the early Holocene L. acervorum still could maintain gene flow between unglaciated parts of the eastern Pyrenees (PY I) and mountains of the Iberian System (e.g. via a network of suitable habitat patches in the Ebro basin). Subsequently, gene flow gradually broke down in parallel with the increasing aridification and the decline of pinewoods in the Ebro basin during the second half of Holocene [69–71]. The limited gene flow between inner Iberia and W-Pyrenees (PY II) at the same time, can be explained by the stronger influence of oceanic climate and consequently the replacement of pinewoods by broad-leafed tree taxa in the Cantabrian Mountains and W-Pyrenees [70, 72, 73].
Our study gives important insight into the evolutionary and demographic history of L. acervorum in the Iberian Peninsula, the genetic relationship among its two social phenotypes and the recolonization of its Western-Palearctic range from glacial refugia north of the Pyrenees. First, our analyses show a lack of reproductive isolation between both social forms, which indicates that behavioural flexibility rather than a hard-wired genetic basis is underlying variation in social phenotype. Second, we demonstrate the occurrence of substantial genetic structure, which is congruent with the complex Iberian topography and SW-European climate history. Moreover, our results highlight the vulnerability of L. acervorum in Iberian mid Mountain ranges to climate change and human induced climatic warming in particular, where the species inhabits suitable microhabitat patches over 1500 m a. s. l. only. A complete picture of the postglacial migration of this Holarctic ant with samples from across its whole range in Eurasia and Northern America as well as new marker sets might help to better understand the dynamics of centre-to-periphery pattern of genetic variability and the history of boreal insect species with high dispersal capacity.
We thank I. Rakitchenkova, M. Schneider, I. Miranda da Silva, and A. Trindl for helping during laboratory work; M. Suefuji for kindly providing unpublished microsatellite primers; A. Buschinger, M. Füßl, M. Haubner, S. Hutschenreuther, P. Koller, R. Landstorfer, S. Leopold, D. Rothgänger, G. Schneider, I. Miranda da Silva and C. Wanke for assistance during field sampling, and A. Keikhosravi, C. Löser, C. Oberprieler, C. Schubart and N. Thiercelin for helpful discussions and suggestions. The work was financed by DFG (He 1623/25).
Availability of data and materials
Voucher samples are stored at the Senckenberg Museum für Naturkunde Görlitz. DNA sequences generated in this study are available on GenBank (for accession numbers see Additional file 2). Further data supporting the conclusions of this article is included within Additional file 2 as well or are available upon request to the authors.
Conceived and designed the study: JH, JT. Performed the field work: JT, JH. Analysed the data: JT, SA. Wrote the paper: JT, JH, SA. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Permission for sampling the species in the EU was not necessary. Animals were preserved in ethanol and/or kept under near-natural conditions.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hewitt GM. The structure of biodiversity–insights from molecular phylogeography. Front Zool. 2004;1:1–16.View ArticleGoogle Scholar
- Schmitt T. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Front Zool. 2007;4:11.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmitt T. Biogeographical and evolutionary importance of the European high mountain systems. Front Zool. 2009;6:9.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmitt T, Varga Z. Extra-Mediterranean refugia: the rule and not the exception. Front Zool. 2012;9:22.View ArticlePubMedPubMed CentralGoogle Scholar
- Eckert CG, Samis KE, Lougheed SC. Genetic variation across species’ geographical ranges: the central–marginal hypothesis and beyond. Mol Ecol. 2008;17:1170–88.View ArticlePubMedGoogle Scholar
- Cassel-Lundhagen A. Peripheral relict populations of widespread species; evolutionary hotspots or just more of the same? In: Habel JC, Assmann T, editors. Relict species. Berlin: Springer; 2010. p. 267–75.View ArticleGoogle Scholar
- Heinze J, Lipski N, Schlehmeyer K, Hōlldobler B. Colony structure and reproduction in the ant, Leptothorax acervorum. Behav Ecol. 1995;6:359–67.View ArticleGoogle Scholar
- Felke M, Buschinger A. Social organization, reproductive behavior and ecology of Leptothorax acervorum (Hymenoptera, Formicidae) from the Sierra de Albarracin in central Spain. Insectes Soc. 1999;46:84–91.View ArticleGoogle Scholar
- Buschinger A. Mono- und Polygynie bei Arten der Gattung Leptothorax Mayr (Hymenoptera Formicidae). Insectes Soc. 1968;15:217–25.View ArticleGoogle Scholar
- Ito F. Functional monogyny of Leptothorax acervorum in northern Japan. Psyche. 1990;97:203–11.View ArticleGoogle Scholar
- Trettin J, Haubner M, Buschinger A, Heinze J. Queen dominance and worker policing control reproduction in a threatened ant. BMC Ecol. 2011;11:21.View ArticlePubMedPubMed CentralGoogle Scholar
- Chapuisat M, Bocherens S, Rosset H. Variable queen number in ant colonies: no impact on queen turnover, inbreeding, and population genetic differentiation in the ant Formica selysi. Evolution. 2004;58:1064–72.View ArticlePubMedGoogle Scholar
- Foster SA. The geography of behaviour: an evolutionary perspective. Trends Ecol Evol. 1999;14:190–5.View ArticlePubMedGoogle Scholar
- Shoemaker DD, Ross KG. Effects of social organization on gene flow in the fire ant Solenopsis invicta. Nature. 1996;383:613–6.View ArticleGoogle Scholar
- West-Eberhard MJ. Developmental plasticity and evolution. New York: Oxford University Press; 2003.Google Scholar
- Czechowski W, Radchenko A, Czechowska W, Vepsäläinen K. The ants of Poland with reference to the myrmecofauna of Europe. Warsaw: Natura Optima Dux Foundation; 2012.Google Scholar
- Seifert B. Ameisen Mittel- Und Nordeuropas. Lutra: Tauer; 2007.Google Scholar
- Brandt M, Fischer-Blass B, Heinze J, Foitzik S. Population structure and the co-evolution between social parasites and their hosts. Mol Ecol. 2007;16:2063–78.View ArticlePubMedGoogle Scholar
- Sambrook J, Russell DW. Molecular cloning: a laboratory manual. New York: Cold Spring Harbor Laboratory Press; 2001.Google Scholar
- Gill RJ, Arce A, Keller L, Hammond RL. Polymorphic social organization in an ant. Proc R Soc Lond B Biol Sci. 2009;276:4423–31.View ArticleGoogle Scholar
- Bourke AFG, Green H, Bruford M. Parentage, reproductive skew and queen turnover in a multiple–queen ant analysed with microsatellites. Proc R Soc Lond B Biol Sci. 1997;264:277–83.View ArticleGoogle Scholar
- Hamaguchi K, Ito Y, Takenaka O. GT dinucleotide repeat polymorphisms in a polygynous ant, Leptothorax spinosior and their use for measurement of relatedness. Naturwissenschaften. 1993;80:179–81.View ArticlePubMedGoogle Scholar
- Suefuji M, Trindl A, Heinze J. Characterization of 13 microsatellite markers for the threatened, slave-making ant Myrmoxenus ravouxi (Formicidae: Myrmicinae). Conserv Genet Resour. 2011;3:229–31.View ArticleGoogle Scholar
- Foitzik S, Haberl M, Gadau J, Heinze J. Mating frequency of Leptothorax nylanderi ant queens determined by microsatellite analysis. Insectes Soc. 1997;44:219–27.View ArticleGoogle Scholar
- Evans JD. Parentage analyses in ant colonies using simple sequence repeat loci. Mol Ecol. 1993;2:393–7.View ArticlePubMedGoogle Scholar
- Weir BS. Genetic data analysis II: methods for discrete population genetic data. Sunderland: Sinauer Associates, Inc. Publishers; 1996.Google Scholar
- Rousset F. GENEPOP’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.View ArticlePubMedGoogle Scholar
- Van Oosterhout C, Hutchinson WF, Wills DP, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.View ArticleGoogle Scholar
- Goudet J. FSTAT (Version 1.2): a computer program to calculate F-statistics. J Hered. 1995;86:485–6.Google Scholar
- Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28:2537–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Langella O. Populations version 1.2.31: a population genetic software. Roscoff: CNRS UPR9034; 1999.Google Scholar
- Nei M, Chesser RK. Estimation of fixation indices and gene diversities. Ann Hum Genet. 1983;47:253–9.View ArticlePubMedGoogle Scholar
- Jost L. GST and its relatives do not measure differentiation. Mol Ecol. 2008;17:4015–26.View ArticlePubMedGoogle Scholar
- Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.View ArticlePubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMedPubMed CentralGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.View ArticlePubMedGoogle Scholar
- Earl DA, von Holdt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.View ArticleGoogle Scholar
- Rosenberg NA. DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4:137–8.View ArticleGoogle Scholar
- Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27:209–20.PubMedGoogle Scholar
- Piry S, Luikart G, Cornuet J-M. BOTTLENECK: a program for detecting recent effective population size reductions from allele data frequencies. J Hered. 1999;90:502–3.View ArticleGoogle Scholar
- Garza J, Williamson E. Detection of reduction in population size using data from microsatellite loci. Mol Ecol. 2001;10:305–18.View ArticlePubMedGoogle Scholar
- Beerli P. Comparison of Bayesian and maximum-likelihood inference of population genetic parameters. Bioinformatics. 2006;22:341–5.View ArticlePubMedGoogle Scholar
- Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994;87:651–701.View ArticleGoogle Scholar
- Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.Google Scholar
- Song H, Buhay JE, Whiting MF, Crandall KA. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proc Natl Acad Sci. 2008;105:13486–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream M-A, Barrell B. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16:944–5.View ArticlePubMedGoogle Scholar
- Clement M, Posada D, Crandall KA. TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000;9:1657–9.View ArticlePubMedGoogle Scholar
- Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.View ArticlePubMedGoogle Scholar
- Pons O, Petit RJ. Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics. 1996;144:1237–45.PubMedPubMed CentralGoogle Scholar
- Trettin J, Seyferth T, Heinze J. Behavioral plasticity in Ant queens: environmental manipulation induces aggression among normally peaceful queens in the socially polymorphic ant Leptothorax acervorum. PLoS One. 2014;9:4.View ArticleGoogle Scholar
- Purcell J, Brelsford A, Wurm Y, Perrin N, Chapuisat M. Convergent genetic architecture underlies social organization in ants. Curr Biol. 2014;24:2728–32.View ArticlePubMedGoogle Scholar
- Wang J, Wurm Y, Nipitwattanaphon M, Riba-Grognuz O, Huang Y-C, Shoemaker D, Keller L. A Y-like social chromosome causes alternative colony organization in fire ants. Nature. 2013;493:664–8.View ArticlePubMedGoogle Scholar
- Walter B, Heinze J. Queen–worker ratio affects reproductive skew in a socially polymorphic ant. Ecol Evol. 2015;5:5609–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Heinze J, Gratiashvili N. High skew in the Caucasus: functional monogyny in the ant Leptothorax scamni. Insectes Soc. 2015. doi:10.1007/s00040-015-0415-5.Google Scholar
- Gómez A, Lunt DH. Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. In: Weiss S, Ferrand N, editors. Phylogeography of southern European refugia. Berlin: Springer; 2007. p. 155–88.View ArticleGoogle Scholar
- Nieto FG. Southern European glacial refugia: a tale of tales. Taxon. 2011;60:365–72.Google Scholar
- Duriez O, Sachet J-M, Ménoni E, Pidancier N, Miquel C, Taberlet P. Phylogeography of the capercaillie in Eurasia: what is the conservation status in the Pyrenees and Cantabrian Mounts? Conserv Genet. 2007;8:513–26.View ArticleGoogle Scholar
- Rodríguez-Muñoz R, Mirol P, Segelbacher G, Fernández A, Tregenza T. Genetic differentiation of an endangered capercaillie (Tetrao urogallus) population at the Southern edge of the species range. Conserv Genet. 2007;8:659–70.View ArticleGoogle Scholar
- Blanco-Fontao B, Fernández-Gil A, Obeso JR, Quevedo M. Diet and habitat selection in Cantabrian Capercaillie (Tetrao urogallus cantabricus): ecological differentiation of a rear-edge population. J Ornithol. 2010;151:269–77.View ArticleGoogle Scholar
- Rubiales JM, García-Amorena I, Álvarez SG, Manzaneque FG. The late Holocene extinction of Pinus sylvestris in the western Cantabrian Range (Spain). J Biogeogr. 2008;35:1840–50.View ArticleGoogle Scholar
- Peery MZ, Kirby R, Reid BN, Stoelting R, Doucet-Bëer E, Robinson S, Vásquez-Carrillo C, Pauli JN, Palsbøll PJ. Reliability of genetic bottleneck tests for detecting recent population declines. Mol Ecol. 2012;21:3403–18.View ArticlePubMedGoogle Scholar
- Foitzik S, Bauer S, Laurent S, Pennings P. Genetic diversity, population structure and sex-biased dispersal in three co-evolving species. J Evol Biol. 2009;22:2470–80.View ArticlePubMedGoogle Scholar
- Birks HJB, Willis KJ. Alpines, trees, and refugia in Europe. Plant Ecol Divers. 2008;1:147–60.View ArticleGoogle Scholar
- Kropf M, Comes HP, Kadereit JW. Past, present and future of mountain species of the French Massif Central–the case of Soldanella alpina L. subsp. alpina (Primulaceae) and a review of other plant and animal studies. J Biogeogr. 2012;39:799–812.View ArticleGoogle Scholar
- Leppänen J, Vepsäläinen K, Anthoni H, Savolainen R. Comparative phylogeography of the ants Myrmica ruginodis and Myrmica rubra. J Biogeogr. 2013;40:479–91.View ArticleGoogle Scholar
- Provan J, Bennett KD. Phylogeographic insights into cryptic glacial refugia. Trends Ecol Evol. 2008;23:564–71.View ArticlePubMedGoogle Scholar
- Stewart JR, Lister AM. Cryptic northern refugia and the origins of the modern biota. Trends Ecol Evol. 2001;16:608–13.View ArticleGoogle Scholar
- Schmitt T, Hewitt GM, Müller P. Disjunct distributions during glacial and interglacial periods in mountain butterflies: Erebia epiphron as an example. J Evol Biol. 2006;19:108–13.View ArticlePubMedGoogle Scholar
- Davis BA, Stevenson AC. The 8.2 ka event and Early–Mid Holocene forests, fires and flooding in the Central Ebro Desert, NE Spain. Quat Sci Rev. 2007;26:1695–712.View ArticleGoogle Scholar
- Carrión JS, Fernández S, González-Sampériz P, Gil-Romera G, Badal E, Carrión-Marco Y, López-Merino L, López-Sáez JA, Fierro E, Burjachs F. Expected trends and surprises in the Lateglacial and Holocene vegetation history of the Iberian Peninsula and Balearic Islands. Rev Palaeobot Palynol. 2010;162:458–75.View ArticleGoogle Scholar
- Pérez-Obiol R, Jalut G, Julià R, Pèlachs A, Iriarte MJ, Otto T, Hernández-Beloqui B. Mid-Holocene vegetation and climatic history of the Iberian Peninsula. The Holocene. 2011;21:75–93.View ArticleGoogle Scholar
- Crawford RMM. Trees by the sea: advantages and disadvantages of oceanic climates. Proc R Ir Acad. 2005;105B:129–39.View ArticleGoogle Scholar
- Rubiales JM, García-Amorena I, Hernández L, Génova M, Martínez F, Manzaneque FG, Morla C. Late Quaternary dynamics of pinewoods in the Iberian mountains. Rev Palaeobot Palynol. 2010;162:476–91.View ArticleGoogle Scholar