Mutualistic interactions are wide-spread but the mechanisms underlying their evolutionary stability and ecological dynamics remain poorly understood. Cultivation mutualisms in which hosts consume symbionts occur in phylogenetically diverse groups, but often have symbiont monocultures for each host. This is consistent with the prediction that symbionts should avoid coexistence with other strains so that host services continue to benefit relatives, but it is less clear whether hosts should always favor monocultures and what mechanisms they might have to manipulate symbiont diversity. Few mutualisms have been studied in sufficient genetic detail to address these issues, so we decided to characterize symbiont diversity in the complex mutualism between multiple root aphid species and Lasius flavus ants. After showing elsewhere that three of these aphid species have low dispersal and mostly if not exclusively asexual reproduction, we here investigate aphid diversity within and between ant nest mounds.
The three focal species (Geoica utricularia, Forda marginata and Tetraneura ulmi) had considerable clonal diversity at the population level. Yet more than half of the ant mounds contained just a single aphid species, a significantly higher percentage than expected from a random distribution. Over 60% of these single-species mounds had a single aphid clone, and clones tended to persist across subsequent years. Whenever multiple species/clones co-occurred in the same mound, they were spatially separated with more than 95% of the aphid chambers containing individuals of a single clone.
L. flavus “husbandry” is characterized by low aphid “livestock” diversity per colony, especially at the nest-chamber level, but it lacks the exclusive monocultures known from other cultivation mutualisms. The ants appear to eat most of the early instar aphids, so that adult aphids are unlikely to face limited phloem resources and scramble competition with other aphids. We suggest that such culling of carbohydrate-providing symbionts for protein ingestion may maintain maximal host yield per aphid while also benefitting the domesticated aphids as long as their clone-mates reproduce successfully. The cost-benefit logic of this type of polyculture husbandry has striking analogies with human farming practices based on slaughtering young animals for meat to maximize milk-production by a carefully regulated adult livestock population.
Mutualistic symbioses are widespread and of crucial importance in many ecosystems
. Although evolutionary theory to explain the stability of mutualistic interactions has progressed considerably (see
 for a review), consensus on the general underlying mechanisms that keep these interactions stable and cooperative has not been achieved
[3–7]. While further theoretical work might alleviate this problem, these difficulties also illustrate that mutualistic interactions are highly variable in their ecological contexts
[8–10] and degrees of commitment
[11–13], and that very few of them have been studied in considerable depth (reviewed in
). Two aspects are thought to have important implications for the interaction stability of host-symbiont mutualisms:
The level of sexual reproduction and the degree of independent dispersal of the symbionts, and
Genetic diversity among symbionts of a single host
. In a previous study we investigated the first aspect in the hitherto poorly studied mutualism of Lasius flavus ants farming root-aphids
. The present study focuses on the second aspect.
In cultivation (farming) mutualisms, the host partner promotes the growth of a symbiont that it consumes, either individually or as somatic modules
. While scenarios of ‘enslavement domestication’ have been suggested for the early evolution of such mutualisms
[16, 17], it remains difficult to understand how symbionts would be actively selected to make the transition from free-living to being domesticated. The latter state would imply becoming reproductively isolated from free-living relatives which would require consistent direct benefits to be sustainable. Domestication often also implies losing options for horizontal transmission, having many offspring consumed by the host, and potentially being mixed with other symbiont lineages, consequences that could all discourage life as a symbiont. Domestication mutualisms would thus seem most likely to evolve if symbiont services ultimately benefit the reproduction of close symbiont relatives and if the productivity of domesticated reproduction consistently exceeds the fitness that can be obtained from a free-living life-style. When symbionts are already clonal before domestication, one would therefore expect symbioses to elaborate this form of propagation when making symbionts commit irreversibly to a dependent life-style, which requires new host-serving adaptations that impede survival and reproduction without the host. The ‘trophobiotic organs’ evolved in the aphids of our present study
[18, 19] are examples of such adaptations.
While symbiont interests in being cultivated would be expected to benefit from monopolizing host attention to a group of close relatives, hosts should not necessarily favor the same tendencies towards rearing monocultures, as a more variable community of symbionts might offer a broader spectrum of services or be less vulnerable to parasites (e.g.
). As outlined by in earlier studies
[21, 22], hosts would be selected to enforce monocultures only if scramble competition between multiple symbiont strains would decrease the overall productivity of the symbiotic interaction, i.e. if different symbiont strains would compete for the same limited resource provided by the host. Similar selection pressure towards monoculture farming would apply if coexistence of multiple strains within the same host would allow free-riding by underperforming strains, leading to a direct reduction in overall productivity (e.g.
Incentives for competition or cheating would destabilize mutualistic interactions between symbionts and hosts, unless specific mechanisms of symbiont screening upon admission
 or symbiont rewarding/sanctioning in proportion to performance
[5, 10] can evolve. The relative importance of these mechanisms is controversial, but available data suggest that monocultures are commonly found in the cultivation mutualisms that have been studied, from the gardens of algae-growing damselfish
 to those of fungus-growing termites and ants
[11, 26–28]. In fungus-farming leaf-cutting ants, monocultures appear to be enforced by a combination of incompatibility between genetically different symbiont strains and active symbiont policing by the hosts
[11, 28, 29], whereas a simple mechanism of positive frequency-dependent propagation within established colonies appears sufficient to enforce life-time commitment between a termite host colony and a single symbiont clone
. However, more studies are needed to establish the generality of this principle, particularly for cultivation mutualisms where hosts are able to segregate symbionts in space or time to avoid competition
, so that the benefits of polyculture might surpass the costs.
In the present study we focus on a farming symbiosis that has been known for decades but has rarely been studied: the root aphid husbandry for sugar (honeydew, “milk”) and nitrogen (“meat”) of the Yellow meadow ant Lasius flavus, which is likely to be essential for ant colony growth and reproduction, and involves an entire array of root aphid species
[18, 30–35]. These root aphid species have a number of distinct traits that improve performance as ant symbionts but are never found in free-living aphids, such as the ‘trophobiotic organ’ to hold honeydew for the ants
. The most common species have further lost most if not all sexual reproduction in Northwest Europe, but have maintained low frequencies of winged morphs that may disperse between colonies
. In the present study we use a newly developed set of DNA microsatellite markers
 to assess aphid species number and clonal diversity at the level of single ant nest mounds.
The objectives of our study were to use hierarchical sampling (Figure
1) and DNA microsatellite analysis to: 1. Estimate species- and clone diversity for three focal species of root aphids (Geoica utricularia, Tetraneura ulmi, Forda marginata) within L. flavus nests, soil samples within nests, and single aphid chambers (Figure
1a) within these soil samples, 2. Evaluate whether the observed distributions are consistent with the expectation that symbiont diversity within nests is low, 3. Analyze the extent to which diversity patterns change across sampling levels and years, and 4. Infer which potential mechanisms can lead to the observed diversity patterns.
Aphid diversity and abundance
As shown in Figure
2, considerable aphid diversity existed along the sampled 7 km transect, but the distribution of this diversity across ant mounds deviated significantly from random. At all sampling levels (ant mound, soil sample and chamber) monocultures containing only a single species occurred much more often than expected from a random distribution (Figure
1), with 52% of the sampled mounds and 99% of the aphid chambers containing only a single species. Also genetic diversity within species was always non-randomly distributed over the mounds, as there were more mounds that contained a single multilocus lineage (MLL) than expected based on the distribution of MLLs over transect locations (Figure
1). The same was true for the distribution of multilocus genotypes (MLGs) over mounds, with G. utricularia MLGs occurring significantly more often in monocultures than expected. In the other two species the frequency of MLG-monocultures across mounds was not significantly different from random expectation (Figure
Results of the monoculture analyses
Within species MLL
Between species MLG
For each organization level (between species and between MLLs and MLGs within species) the probability (P) that the observed number of monocultures at a given sampling level (ant mounds, soil samples and aphid chambers) could have resulted from a random distribution of aphids was estimated using a bootstrap approach with 1000 iterations. P values below 0.05 (bold Figures) indicate deviations from a random distribution.
At lower sampling levels within mounds (soil samples, chambers) high percentages of monocultures were also found, both between and within species (Table
1). However, these monoculture percentages did mostly not significantly deviate from randomness, because low aphid diversity at the species, MLL or MLG level across mounds or soil samples will automatically lead to low aphid diversity at the next level below. Figure
3 illustrates this for the spatial distribution of G. utricularia MLGs in one of the nests of Figure
2, showing that most MLGs occurred spatially separated already at the soil sample level, so that aphid chambers could only contain monocultures (Table
4 shows the distribution of aphid numbers per chamber, with most chambers containing only one aphid, but some chambers having as many as 13 aphids (means per chamber ± SE G. utricularia 1.61 ± 0.13, T. ulmi 1.84 ± 0.18, F. marginata 2.39 ± 0.28). Even aphid chambers with rather many aphids often contained monocultures in terms of MLLs (Figure
4), and chambers that did contain polycultures never had more than 2 MLLs.
Annual turnover of aphid clonal lineages
Ten of the ant mounds sampled in 2008 were resampled in 2009 and 2010. In seven of these we found one or more of the focal species in the subsequent years (Figure
5). Most MLGs that we found in later years had already been found in the same mounds in 2008. There were only two exceptions to this apparent continuity over time: in the first mound resampled for T. ulmi (Figure
5b) we found an additional MLG in 2009 that had not been observed in that mound in the previous year, and in the second mound resampled for F. marginata (Figure
5c), we found a MLG that had not been identified before, but which belonged to one of the MLLs that had been observed in 2008 in other nest mounds nearby (colored in green shade, Figure
2). These apparent exceptions might either reflect recent colonization events or might be due to under-sampling in 2008. For example, the overall composition of the nest mound in which the F. marginata MLG was newly observed did not change significantly between 2008 and 2009 (Fisher Exact Test, P = 0.111), likely because the newly observed MLG belonged to a MLL that had a population-wide frequency of 0.044 in 2008. With such a low frequency, it is quite likely that this MLG was missed in an earlier year. In contrast, the overall composition of aphid MLGs in the nest where we found a new MLG for T. ulmi did significantly change between 2008 and 2009 (Fisher Exact Test, P = 0.024). It thus appears less likely that the new MLG was due to under-sampling in 2008, since this MLG occurred at a high frequency overall (0.52). Overall, we infer that clonal lineage composition of aphid livestock in L. flavus ant mounds changes relatively little from year to year. We would have liked to test this with a formal heterogeneity analysis across years, but too low numbers in several cells precluded this.
Aphid distribution and abundance
In half of the ant mounds sampled in 2008 only one of the three focal root aphid species, G. utricularia, T. ulmi or F. marginata was found, despite the other aphid species being present within a radius of 50 m. This level of aphid specificity among ant nests matches earlier findings by
 in a British field survey of the same ant species and its underground aphids. If there were multiple aphid species per mound, we found that they tended to be clustered in separate soil samples and hardly ever occurred in the same aphid chamber. This not only applied for the three most abundant species that we focused on, but also for other rarer species of root aphids. We are confident, therefore, that inclusion of these other aphids would not have changed our overall conclusions. Unfortunately, aphid sample sizes remained low for most mounds and for two of the three focal species, which seems unavoidable as earlier non-destructive large-scale surveys obtained similar numbers of adult root aphids for these species (
Additional file 1). Comparing frequencies and absolute densities across studies is not easy as authors have used different sampling methods in the past
[33, 37]. While these have given very different estimates of root aphid density, we show in the
Additional file 1 that this is almost certainly due to these sampling differences, and that our estimates of adult aphid densities are in line with previous studies. In spite of these sample size limitations, we are confident, therefore, that our results would be repeatable with larger sample sizes at: (1) the mound level, because aphids were generally found scattered throughout the entire mound, so that systematic bias in our non-destructive sampling appears unlikely, and (2) at the chamber level, because Figure
4 shows that within-chamber aphid diversity does not increase with increasing numbers of aphids per chamber (ca. equivalent to chamber size). Moreover, in our statistical analysis we control for any effect of the low sample sizes, by simulating the exact same sample sizes as achieved in the field.
The considerable interaction-specificity, often between single ant colonies and single aphid lineages was also encountered at the genetic level within species. Mounds often only harbored one clonal lineage of a single aphid species and if mounds had multiple aphid clones they were almost always compartmentalized in different chambers. A similar degree of host specificity has also been shown for above-ground aphids tended by ants as opposed to non-tended aphids
. However, complete spatial separation of aphid clones is less frequently observed above ground
[38, 39], probably because these aphids can more easily move around. Our limited sampling across years further indicated a high degree of constancy of distributions of aphid clones over time. It would have been interesting to compare our results with similar studies on other myrmecophilous and non-myrmecophilous root aphids, but to our knowledge such studies have not yet been done.
Within- and between ant mound aphid distribution patterns
The between- and within-mound distribution patterns were very similar for the three aphid species under study (Figure
2), suggesting that similar dispersal and recruitment mechanisms apply. Most aphid chambers contained only a single aphid (Figure
4) and chamber sizes seemed proportional to the number of aphids housed in them (
Additional file 1). It thus appears unlikely that these aphids competed for limiting phloem resources, even in the few cases where different aphid species or MLLs shared a chamber. Rather, the husbanding ants seem to optimize the feeding conditions for each aphid adult, because aphid densities (on average 1.00 per liter soil,
Additional file 1) remained well below densities that would occupy all available root phloem resources. These relatively low numbers of adult aphids can be explained by the ants eating the vast majority of aphid nymphs and only keeping a small number of adults for honeydew production as inferred previously by
The low aphid diversity per mound, the apparent invariance of clonal distributions per mound among years, and the high degree of population viscosity
 are consistent with horizontal transmission of aphids between mounds being infrequent. After successful dispersal and adoption, aphid propagation within mounds would then mostly be in the form of clonal copies of fundatrices (aphid ‘foundresses’) replacing their ancestors. We would thus expect that the genetic diversity of aphid livestock within a given ant mound would slowly increase over the years. The densest L. flavus populations in Northwest Europe are normally found in extensively grazed old pastures that have been stable for centuries and where nest-mounds are large because many generations of L. flavus colonies have contributed to building them. Compared to such populations, the coastal transect that we studied is more variable in age and stability, which appeared to be reflected in the younger parts of the salt-marsh harboring less aphid diversity, at least for G. utricularia (
Additional file 1). Patterns like this would be reminiscent of older trees having richer communities of underground mycorrhiza and leaf-endophytes
[40, 41], but also of above-ground aphid colonies becoming more genetically diverse over the season due to the immigration of new aphid clones
Inferring the evolutionary logic of aphid husbandry in Lasius flavus colonies
Genetic diversity of symbionts has been a central issue in mutualism theory
[3, 21, 22] as diversity levels that simultaneously maximize the fitness of both hosts and symbionts are often expected to be low
[21, 22]. This is indeed what we found throughout our data set (i.e. at the species, MLL and MLG level). Compartmentalization of symbionts is known to promote mutualism stability in other systems (e.g. mycorrhizal mutualisms
), because benefits can be preferentially allocated toward cooperative symbionts. However, many of these conceptual arguments are based on the assumption that symbiont lineages compete and that the collateral damage of such interactions for host fitness maintains selection to suppress symbiont diversity
. While the high root aphid densities per L. flavus mound reported in the literature (
Additional file 1) that inspired this study suggested that such competition might also apply in this system, our results prompt us to reappraise this assumption, because: 1. Aphid husbandry is special, relative to other resource enhancing mutualisms, in that the ant hosts can exploit their aphid symbionts both for sugars (“milking” adults in their prime age) and for proteins (eating young instars and old adults) and 2. Our data suggest that consumption of most of the aphid offspring by the ants reduces total aphid numbers per mound (
Additional file 1) to such extent that the grass-root phloem resource constraints that might have induced aphid competition are unlikely to apply.
Many details of the interaction between L. flavus ants and their communities of mutualistic root aphids remain unknown and deserve further study. However, our present results indicate that the biological details and specific resource constraints of an obligate mutualism may be decisive for the selection factors that determine evolutionary stability over time. Our present data indicate that prevailing paradigms of partner choice and sanctions
[5, 10, 24, 43] may not apply in the ant-aphid mutualism that we studied, because fundamental assumptions of scramble competition between unrelated symbionts
 are not fulfilled (
Additional file 1).
After initial domestication, the aphid clones would have continued to benefit from the symbiosis, because the premature death of most early instar nymphs (which individually are of low value as sugar providers for the ants) reduces competition over resources and reproduction, and extensive clonality ensures that vertical transmission will maintain clonal tenure within nests. This interpretation might explain why L. flavus is reputedly obligately dependent on root aphids
[33, 35], but without having specialized on any of the large number of aphid species that are available, despite the aphids having evolved specialized traits that enhance productivity as ant symbionts but preclude independent life (see
 for details). Testing the validity of our interpretation that early instar aphids are worth more as direct sources of protein than as later sources of carbohydrates will require controlled lab experiments, which might be feasible in spite of the challenges of keeping these ants and aphids in artificial nests
Analogies with human subsistence farming
The results of our study suggest that polyculture aphid husbandry in L. flavus follows similar efficiency principles as modern cattle husbandry practices in humans, where adult cows are kept in numbers that secure maximal milk-productivity in a competition-free environment and where surplus reproduction is slaughtered for meat-consumption soon after birth. How this analogy could come about is interesting to evaluate.
The English name for L. flavus, Yellow meadow ant, indicates lack of pigmentation because the ants are almost never exposed to direct sunlight. This exclusively underground life style, shared by many but far from all Lasius ants
, must have implied that foraging territories became limited to the direct nest environment, so that access to prey was reduced but protection and monopolization of domesticated aphids became easier. Intensification of aphid husbandry thus seems a logical consequence of becoming subterranean and a prudent way of harvesting a small local resource-base that ultimately depends on primary production (grass roots) rather than secondary production (free-living prey capture). Extensive culling of immature aphids for meat not only allowed polyculture practices (by eliminating competition), but may have actively encouraged it when different aphid species would exploit somewhat different plant root niches, when their availability would be unpredictable, or when they would produce honeydew with slightly different chemical composition
The analogies between aphid husbandry in L. flavus and human cultural practices are quite striking as farming husbandry allowed human populations to sustain themselves at much higher densities than hunter-gatherer populations
. Likewise, the density of L. flavus ants in mature grasslands is among the highest known for ants
[35, 47, 48] and appears to be sustainable with only a modest ecological footprint. As in humans, the secret of success appears to be a unique combination of traits, such as the ability to actively engineer nest mound habitat (a form of niche construction
) rather than living in fixed plant structures as other obligately aphid-dependent ants do (e.g.
), and the availability of multiple aphids that could be domesticated without the need to specialize on any one of them. This suggests that ant farming practices for meat
[50, 51] deserve more explicit study, as they may provide remarkable insights into sustainable farming practice.
Farming mutualisms are highly diverse. Some have a long history of coadaptation, specificity and vertical symbiont transmission, whereas others have evolved interdependences based on horizontal symbiont acquisition and low specificity. Many ant species obtain facultative benefits from tending aphids. Some of these interactions have evolved to be highly specific, but the Lasius flavus husbandry system that we studied is unusual in that both ants and root aphids appear to be obligately interdependent and adapted to their respective life styles as farmers and livestock, but without obvious signs of species-by-species interaction specificity.
Our genetic explorations of a large island population with dense populations of L. flavus suggest that the combination of permanently underground nesting, aphid clonality, and very low gene flow between aphid populations of neighboring mounds has allowed these ants to evolve an unusual form of polyculture symbiosis. Species and clonal lineages of aphids appear to be kept apart, which likely gives colonies the possibility to actively manage the diversity and abundance of their livestock. We hypothesize that this allows the ants to secure maximal yield from a subset of mature aphids that are kept for carbohydrates under optimal conditions of phloem feeding and ant care. These selected aphids may then also reproduce at the highest possible rate, so that the ants both secure maximal protein intake by eating the excess of early instar aphids, and replacement of their honeydew-producing livestock when adult aphids age and become less productive.
Many mechanistic details that govern the dynamics of this mutualism await further research. However, we feel that analogies with human husbandry practices based on similar cost-benefit considerations lend sufficient credibility to our interpretations to generate novel interest into natural selection processes that have produced ant farming practices for both meat and carbohydrates.
Natural history of the model system
The subterranean Yellow meadow ant Lasius flavus constructs conspicuous nest mounds (↕ ca. 30 cm, ø ca. 80 cm) in grassland habitats to house both its own colonies and the root aphids on which it depends for honeydew as a source of carbohydrates
[18, 30, 32, 33] and which they eat for protein
[30, 32, 33]. The ants actively protect the aphids (
, A.B.F. Ivens, personal observation) and keep them in specially constructed ‘aphid-chambers’, cavities around grass-roots with one or several aphids (Figure
1a). Thirteen species of root aphids are known to be tended by L. flavus, often with multiple species in the same nest mound (
[18, 31, 33, 34], A.B.F. Ivens, personal observation). Among these, Tetraneura ulmi, Geoica utricularia and Forda marginata are often the most dominant species
[14, 18, 31, 33, 34]. This was also the case at our study site, so we focused our study on these three species. These aphids can also be found in nests of other ants, such as Myrmica sp. and other Lasius species
, albeit in lower numbers than in the typical L. flavus mounds.
Aphid reproductive cycles can be fully asexual (anholocyclic) or include a single sexual phase at the end of the season (holocyclic). In another study
 we showed that the three focal aphid species are predominantly if not completely asexual in our study population (see also below), consistent with all three species having been shown to feed year-round on roots of the grasses Festuca rubra, Agrostis spp. and Elytrigia maritima without requiring a host shift during winter
[14, 18, 31, 33]. The possible winter host shift to Ulmus trees that has previously been described for Tetraneura ulmi (
, O.E. Heie, personal communication) thus appears to be absent in our NW European study population. However, several other mechanisms can account for more limited horizontal aphid dispersal in salt march habitats such as our study site: walking, floating on tidal floods and wind dispersal of winged individuals (alates) that are produced at very low frequencies in all three species. Considerable genetic population viscosity confirms that horizontal dispersal between mounds is generally limited
. However, this appears to be the only dispersal mode available as neither vertical nor horizontal transmission by the tending workers has ever been observed for Lasius ants.
Ant mounds were sampled for aphids on the island of Schiermonnikoog, the Netherlands (53°28′ N, 6°09′ E) in July 2008, 2009 and 2010 along a 7 km transect across most of the salt-marsh on the island (Figure
1b). The westernmost first kilometer close to the inhabited part of the island was grazed by cattle, whereas the remaining transect crossed ungrazed salt-marsh. The transect was subdivided into eight locations (one every km). At every location we sampled five same-sized ant mounds (Ø ca. 60 cm), taking 21 cylindrical soil samples (10 cm deep and Ø 8 cm, volume 0.64 l), according to a fixed sampling scheme (Figure
1b). The average volume of the part of the mounds that was suitable for aphids (i.e. had roots of the appropriate grasses) was 66.7 l (
Additional file 1). We obtained this estimate by adding the volume of the aboveground part of an average mound and the volume of a ring directly surrounding the mound (10 cm wide, 8 cm deep) which is known to often contain root aphids as well
Every soil sample was hand-sorted for ‘aphid chambers’, cavities containing one or more root aphid individuals in spatial isolation from any other aphids (Figure
1a). This sampling scheme resulted in a four-level nested design: transect location, ant nest mound, soil sample and aphid chamber (Figure
1b). In July 2009 and 2010 we resampled 10 of the 40 previously sampled nest mounds, for which we had obtained sufficiently detailed aphid distributions in 2008 to be able to detect changes in later years.
Molecular methods and data analysis
A detailed description of the molecular analysis of the aphids and properties of the genetic markers is provided by
. In short, all collected aphids were genotyped for an array of polymorphic microsatellite markers (Geoica utricularia, eight markers: Gu2, Gu3, Gu5, Gu6, Gu8, Gu9, Gu11, Gu13; Forda marginata, seven markers: Fm1, Fm3, Fm4, Fm6, Gu6, Gu11, Gu13; Tetraneura ulmi, six markers: Tu1, Tu2, Tu3, Tu4, Tu10, Tu11) after DNA extraction from entire aphids using 200 μl 2yChelex® 100 resin (Fluka)
. Following PCR-amplification, products were analyzed on an ABI-PRISM 3130XL (Applied Biosystems) sequencer and chromatograms were analyzed in Genemapper (Applied Biosystems).
When amplification failed, samples were re-run at least two more times. When amplification remained unsatisfactory, the specific microsatellite locus was scored as ‘missing data’. When data were missing for more than half of the loci the individual was omitted from further analysis. In total, we included 239 individuals of Geoica utricularia (2008: 201, 2009: 23, 2010: 15, after omitting a total of 28 individuals), 191 of Forda marginata (2008: 158, 2009: 4, 2010: 29, 11 omitted) and 105 of Tetraneura ulmi (2008: 92, 2009: 7, 2010: 6, 4 omitted).
Diploid clonal multilocus genotypes (MLGs) consist of a unique combination of alleles across all genotyped loci. The genotypic data allowed us to assign every aphid to a MLG using the software MLGSIM2.0
, an updated version of MLGSIM
. A multilocus lineage (MLL) is a group of closely related MLGs that differ by only one or two alleles
. All MLGs could be grouped into MLLs. The complete analysis is detailed in
When a sample only contained aphids from a single species, MLL or MLG, we classified that sample as a ‘monoculture’ at the species, MLL or MLG level. Samples were taken at three ‘sampling levels’: ant mound, soil sample or aphid chamber. To test whether the observed monocultures occurred more frequently than expected under a random distribution, we wrote a bootstrap routine in R 2.13.0
 (routine available upon request). For a given level of sampling, the routine distributed the species, MLLs, or MLGs randomly over samples in 1000 iterations with simulation sample sizes corresponding to the observed sample sizes. The routine thus used the observed frequency distributions of species, MLLs, or MLGs at the sampling level above the focal level (Figure
1b) to estimate the probability (P) that the same or a higher number of monocultures than the observed number would be obtained by chance (one-tailed test). When P was < 0.05, the null hypothesis that the observed number of monocultures resulted from a random distribution of aphids over samples was rejected.
Any behavioral experiments referred to were conducted with insects. In this context, insects are not considered animals, so these experiments do not, as far as we are aware, require an formal approval of an ethics committee.
All authors have long standing interest in the evolutionary biology of cooperation and conflict within and between species. JJB, DJCK and ABFI mainly focus on different aspects of social evolution of ants and their mutualistic partners. FJW and IP are theoretical biologists, focusing on the theory of conflict and cooperation in social systems.
We thank Thijs Janzen and Sylvia Mathiasen for help in the field and in the laboratory, Søren Feodor Nielsen for statistical advice, Ruth Howison for help with the analysis of the elevation data (Additional file
1) and David Nash and two anonymous reviewers for comments on earlier versions of the manuscript. A PhD-Grant of the Centre for Ecological and Evolutionary Studies, University of Groningen and grants by the Langerhuizen Fund and the Nicolaas Mulerius Fund, awarded to ABFI supported this study, as well as a grant by the Danish National Research Foundation to JJB. We are grateful to ‘Vereniging Natuurmonumenten’ for permission to do field work in the national nature reserve Schiermonnikoog.
1Theoretical Biology, Centre for Ecological and Evolutionary Studies, University of Groningen, Nijenborgh 7, Groningen 9747 AG, The Netherlands.2Centre for Social Evolution, Department of Biology, University of Copenhagen, Universitetsparken 15, Copenhagen DK-2100, Denmark.3Laboratory of Insect Social Evolution, The Rockefeller University, 1230 York Avenue, New York, NY 10065, USA.
Theoretical Biology, Centre for Ecological and Evolutionary Studies, University of Groningen
Centre for Social Evolution, Department of Biology, University of Copenhagen
Laboratory of Insect Social Evolution, The Rockefeller University
Stachowicz JJ: Mutualism, facilitation, and the structure of ecological communities.Bioscience 2001, 51:235–246.View Article
Herre EA, Knowlton N, Mueller UG, Rehner SA: The evolution of mutualisms: exploring the paths between conflict and cooperation.Trends Ecol Evol 1999, 14:49–53.PubMedView Article
Edwards DP, Hassall M, Sutherland WJ, Yu DW: Selection for protection in an ant-plant mutualism: host sanctions, host modularity and the principal-agent game.Proc R Soc Lond Ser B 2006, 273:595–602.View Article
Weyl EG, Frederickson ME, Yu DW, Pierce NE: Economic contract theory tests models of mutualism.Proc Natl Acad Sci U S A 2010, 107:15712–15716.PubMedView Article
Kiers ET, Denison RF, Kawakita A, Herre EA: The biological reality of host sanctions and partner fidelity.Proc Natl Acad Sci U S A 2011, 108:E7-E7.PubMedView Article
Weyl EG, Frederickson ME, Yu DW, Pierce NE: Reply to Kiers et al.: Economic and biological clarity in the theory of mutualism.Proc Natl Acad Sci U S A 2011, 108:E8.View Article
Bshary R, Grutter AS: Image scoring and cooperation in a cleaner fish mutualism.Nature 2006, 441:975–978.PubMedView Article
Palmer TM, Doak DF, Stanton ML, Bronstein JL, Kiers ET, Young TP, Goheen JR, Pringle RM: Synergy of multiple partners, including freeloaders, increases host fitness in a multispecies mutualism.Proc Natl Acad Sci U S A 2010, 107:17234–17239.PubMedView Article
Kiers ET, Denison RF: Sanctions, cooperation, and the stability of plant-rhizosphere mutualisms.Annu Rev Ecol Evol Syst 2008, 39:215–236.View Article
Poulsen M, Boomsma JJ: Mutualistic fungi control crop diversity in fungus-growing ants.Science 2005, 307:741–744.PubMedView Article
Thompson J, Fernandez C: Temporal dynamics of antagonism and mutualism in a geographically variable plant-insect interaction.Ecology 2006, 87:103–112.PubMedView Article
Moran NA, McCutcheon JP, Nakabachi A: Genomics and evolution of heritable bacterial symbionts.Annu Rev Genet 2008, 42:165–190.PubMedView Article
Ivens ABF, Kronauer DJC, Pen I, Weissing FJ, Boomsma JJ: Reproduction and dispersal in an ant-associated root aphid community.Molecular Ecology 2012, 21:4257–4269.PubMedView Article
Hata H, Kato M: A novel obligate cultivation mutualism between damselfish and Polysiphonia algae.Biol Lett 2006, 2:593–596.PubMedView Article
West SA, Griffin AS, Gardner A: Evolutionary explanations for cooperation.Curr Biol 2007, 17:R661-R672.PubMedView Article
Login FH, Balmand S, Vallier A, Vincent-Monégat C, Vigneron A, Weiss-Gayet M, Rochat D, Heddi A: Antimicrobial peptides keep insect endosymbionts under control.Science 2011, 334:362–365.PubMedView Article
Heie OE: The Aphidoidea (Hemiptera) of Fennoscandia and Denmark. I.Fauna Entomol Scand 1980, 9:1–236.
Hölldobler B, Wilson EO: The Ants. Harvard University Press, Cambridge, MA; 1990:522–529.
van Borm S, Buschinger A, Boomsma J, Billen J: Tetraponera ants have gut symbionts related to nitrogen-fixing root-nodule bacteria.Proc R Soc Lond Ser B Biol Sci 2002, 269:2023–2027.View Article
Wulff JL: Clonal organisms and the evolution of mutualism. In Population biology and evolution of clonal organisms. Edited by: Jackson JBC, Buss LW, Cook RE. Yale University Press, New Haven, CT; 1985:437–466.
Frank SA: Host-symbiont conflict over the mixing of symbiotic lineages.Proc R Soc London, Ser B 1996, 263:339–344.View Article
Bronstein JL: The exploitation of mutualisms.Ecol Lett 2001, 4:277–287.View Article
Archetti M, Úbeda F, Fudenberg D, Green J, Pierce NE, Yu DW: Let the right one in: a microeconomic approach to partner choice in mutualisms.Am Nat 2011, 177:75–85.PubMedView Article
Hata H, Kato M: Weeding by the herbivorous damselfish Stegastes nigricans in nearly monocultural algae farms.Mar Ecol Prog Ser 2002, 237:227–231.View Article
Bot ANM, Rehner SA, Boomsma JJ: Partial incompatibility between ants and symbiotic fungi in two sympatric species of Acromyrmex leaf-cutting ants.Evolution 2001, 55:1980–1991.PubMed
Mueller UG, Scott JJ, Ishak HD, Cooper M, Rodrigues A: Monoculture of leafcutter ant gardens.PLoS One 2010,5(9):e12668.PubMedView Article
Ivens ABF, Nash DR, Poulsen M, Boomsma JJ: Caste-specific symbiont policing by workers of Acromyrmex fungus-growing ants.Behav Ecol 2009, 20:378–384.View Article
Pontin AJ: A preliminary note on the eating of aphids by ants of the genus Lasius (Hym., Formicidae).Entomol Mon Mag 1958, 94:9–11.
Muir DA: The ant-aphid-plant relationship in West Dunbartonshire.J Anim Ecol 1959, 28:133–140.View Article
Pontin AJ: The prey of Lasius niger (L.) and L. flavus (F.) (Hym., Formicidae).Entomol Mon Mag 1961, 97:135–137.
Pontin AJ: Numbers and distribution of subterranean aphids and their exploitation by the ant Lasius flavus (Fabr.).Ecol Entomol 1978, 3:203–207.View Article
Godske L: Aphids in nests of Lasius flavus F. in Denmark I: faunistic description.Entomol Medd 1991, 59:85–89.
Seifert B: Die Ameisen Mittel- und Nordeuropas. Lutra, Görlitz; 2007:280.
Ivens ABF, Kronauer DJC, Boomsma JJ: Characterisation and cross-amplification of polymorphic microsatellite loci in ant-associated root-aphids.Conserv Genet Resour 2011, 3:73–77.View Article
Godske L: Aphids in nests of Lasius flavus F. in Denmark II: population dynamics.Entomol Medd 1992, 60:21–26.
Yao I: Contrasting patterns of genetic structure and dispersal ability in ant-attended and non-attended Tuberculatus aphids.Biol Lett 2010, 6:282–286.PubMedView Article
Vantaux A, Billen J, Wenseleers T: Levels of clonal mixing in the black bean aphid Aphis fabae, a facultative ant mutualist.Mol Ecol 2011, 20:4772–4785.PubMedView Article
Palfner G, Casanova-Katny M, Read D: The mycorrhizal community in a forest chronosequence of Sitka spruce [Picea sitchensis (Bong.) Carr.] in Northern England.Mycorrhiza 2005, 15:571–579.PubMedView Article
Arnold A, Herre E: Canopy cover and leaf age affect colonization by tropical fungal endophytes: Ecological pattern and process in Theobroma cacao (Malvaceae).Mycologia 2003, 95:388–398.PubMedView Article
Bever JD, Richardson SC, Lawrence BM, Holmes J, Watson M: Preferential allocation to beneficial symbiont with spatial structure maintains mycorrhizal mutualism.Ecol Lett 2009, 12:13–21.PubMedView Article
Noë R, Hammerstein P: Biological markets - supply and demand determine the effect of partner choice in cooperation, mutualism and mating.Behav Ecol Sociobiol 1994, 35:1–11.View Article
Smart L: Asymmetric interactions between ants, aphids and plants. PhD thesis. University of Bath, Bath; 1991:206–280.
Fischer MK, Shingleton AW: Host plant and ants influence the honeydew sugar composition of aphids.Funct Ecol 2001, 15:544–550.View Article
Larsen CS: Biological changes in human populations with agriculture.Annu Rev Anthropol 1995, 24:185–213.View Article
Odum EP, Pontin AJ: Population density of underground ant, Lasius Flavus, as determined by tagging with P-32.Ecology 1961, 42:186–188.View Article
Boomsma JJ, Van Loon AJ: Structure and diversity of ant communities in successive coastal dune valleys.J Anim Ecol 1982, 51:957–974.View Article
Stenberg P, Lundmark M, Saura A: MLGsim: a program for detecting clones using a simulation approach.Mol Ecol Notes 2003, 3:329–331.View Article
R Development Core Team: R: A language and environment for statistical computing Vienna. R Foundation for Statistical Computing, Austria; 2011.
Pontin AJ: Population stabilization and competition between the ants Lasius flavus (F) and L. niger (L).J Anim Ecol 1961, 30:47–54.View Article
Pontin AJ: Further considerations of competition and the ecology of the ants Lasius flavus (F) and L. niger (L).J Anim Ecol 1963, 32:565–574.View Article
Olff H, De Leeuw J, Bakker J, Platerink R, Van Wijnen H, De Munck W: Vegetation succession and herbivory in a salt marsh: changes induced by sea level rise and silt deposition along an elevational gradient.J Ecol 1997, 85:799–814.View Article
Boomsma JJ, Isaaks JA: Effects of inundation and salt on the survival of ants in a sandy coastal plain.Ecol Entomol 1982, 7:121–130.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (
http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.