Email updates

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

Open Access Research article

Geographic parthenogenesis and plant-enemy interactions in the common dandelion

Koen JF Verhoeven* and Arjen Biere

Author Affiliations

Department of Terrestrial Ecology, Netherlands Institute of Ecology (NIOO-KNAW), Droevendaalsesteeg 10, Wageningen, 6708 PB, The Netherlands

For all author emails, please log on.

BMC Evolutionary Biology 2013, 13:23  doi:10.1186/1471-2148-13-23


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


Received:26 June 2012
Accepted:16 January 2013
Published:28 January 2013

© 2013 Verhoeven and Biere; licensee BioMed Central Ltd.

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

Abstract

Background

Many species with sexual and asexual variants show a pattern of geographic parthenogenesis where asexuals have broader and higher-latitude distribution than sexuals. Because sexual reproduction is often considered a costly evolutionary strategy that is advantageous in the face of selection by coevolving pests and pathogens, one possible explanation for geographic parthenogenesis is that populations at higher latitudes are exposed to fewer pests and pathogens. We tested this hypothesis in the common dandelion (Taraxacum officinale), a species with well-established geographic parthenogenesis, by screening prevalence and effects of several specialized pests and pathogens in natural dandelion populations.

Results

We did a population survey of 18 dandelion populations along a geographic transect that ranged from the area where sexual and asexual dandelions co-occur northward into the area where only asexuals occur. In addition we used four southern and four northern populations in a 8x8 cross-inoculation greenhouse experiment in which plants were exposed experimentally to each other’s natural field soil microbial communities. The cross-inoculation experiment indicated a higher pathogenicity of soil microbial communities from the southern, mostly sexual, populations compared to soil microbial communities from the northern asexual populations. Northern dandelion populations also showed reduced infestation by a specialized seed-eating weevil. A similar trend of reduced rust fungus infection in northern populations was observed but this trend was not statistically significant.

Conclusions

The prevalence of pests and pathogens decreased along the south-to-north axis of geographic parthenogenesis. This highlights the potential of biotic interactions in shaping patterns of geographic parthenogenesis.

Keywords:
Red Queen; Taraxacum; Plant-insect interactions; Plant-pathogen interactions; Soil feedback

Background

Geographic parthenogenesis, or the broader and higher-latitude geographic distribution of asexuals compared to their sexual counterparts in taxa where both types coexist, is well-established in both plants and animals [1-5]. Several non-exclusive explanations have been proposed to account for this pattern [6-10]. Popular hypotheses include higher phenotypic plasticity in asexuals due to efficient selection for general-purpose genotypes [11]; better colonizing abilities of asexuals that facilitate range expansion e.g. into previously glaciated areas [1,12]; and reduced pressure from pests and pathogens at higher latitudes that allows asexuals to capitalize on their reproductive advantage compared to sexuals [13,14].

The latter ‘biotic interactions’ hypothesis rests on two assumptions. First, it builds on the idea that sexual reproduction is inherently costly compared to asexual reproduction but is advantageous in the face of specialized coevolving pests and pathogens, presumably because sex generates genetically variable offspring which reduces the risk of infection. This Red Queen (RQ) hypothesis for the evolutionary maintenance of sexual reproduction [13,15,16] has received empirical support in several systems and ecological studies have confirmed key assumptions and predictions of the hypothesis [17-20]. Still, as modeling studies show [21], it remains under debate how general the RQ hypothesis is in explaining host sexual reproduction and many alternative hypotheses for the maintenance of sexual reproduction exist [22].

The second assumption is that herbivory and pathogen pressure is less intense at higher latitudes. This has been a commonly accepted idea with considerable empirical support [23-25]. Recently, however, the existence of a latitudinal trend in plant herbivory and associated plant defenses was challenged by meta-analysis [26] and large-scale empirical screening, which both failed to support the idea [27]. However, as argued by Johnson and Rasmann [28], the patterns predicted by Coley and Aide [25] may be more readily observed by comparing geographical trends in species with similar life history and roughly comparable habitats than in the type of meta-analyses conducted by Moles et al. [26], which include a wide variety of plant growth forms, habitats, and sections of continuous latitudinal ranges. Plant-pathogen interactions have been less studied than plant-herbivore interactions in a latitudinal context [24] and relatively little empirical evidence exists to support the idea that plant pathogen pressure is lower at higher latitudes. These studies have typically been performed either across small latitudinal ranges, spanning less than one [29] up to fifteen degrees latitude [30-32], or on wider geographic scales but focusing on microbial diversity in general rather than on pathogen pressure [33]. If a pattern of decreasing herbivory and pathogen pressure with increasing latitude is in fact absent then the biotic interactions hypothesis cannot be a general explanation for geographic parthenogenesis.

In plants, little empirical support exists that geographic parthenogenesis is correlated with relaxed pest and pathogen pressure in high-latitude asexual populations compared to low-latitude sexual populations [34-36]. Here, we test for this pattern in the common dandelion (Taraxacum officinale sensu lato, Asteraceae), a species with well-established geographic parthenogenesis [37-39], by linking dandelion sex-asex variation to geographic variation in the occurrence of several specialized dandelion pathogens and herbivores. Sexual T. officinale are diploid self-incompatible plants whereas asexual T. officinale are triploid obligate apomicts; both forms co-occur in south-central Europe but only the asexuals extend much further northward into Scandinavia. First, we experimentally evaluated pathogenic effects of natural soil microbial communities by comparing pathogenicity between soils from southern (predominantly sexual) populations and soils from northern (asexual) populations. Soil pathogens are little explored in a RQ context but they play important roles in shaping plant communities [40]. Second, we surveyed natural dandelion populations along the south–north axis of geographic parthenogenesis to record field prevalence of infection by a specialized rust fungus and a specialized seed-eating weevil. Infection prevalence does not necessarily reflect pathogen pressure and the risk of infection, because prevalence also depends on host susceptibility which may differ between populations. However, prevalence is often used as an initial proxy for population-level differences in infection risk, and this can be a reasonable proxy especially when between-population differences in prevalence are large [41].

The overall hypothesis that we test is that the asexual dandelion populations in the northern part of the transect are exposed to lower herbivore and pathogen pressure compared to the mostly sexual dandelion populations from the southern part of the transect. This tests a key assumption of the hypothesis that antagonists are causally responsible for dandelion geographic parthenogenesis.

Results

Sex-asex variation between dandelion populations

Ploidy analysis confirmed the previously established pattern of geographic parthenogenesis in this species. Populations from the northern part of the transect (Sweden, Denmark, Netherlands, Belgium) were all triploid and thus apomictic. The southern part (Switzerland, eastern France, southern Germany) contained apomictic, sexual and mixed populations (Table 1).

Table 1. Characteristics and infection proportions of theT. officinalepopulations

Soil inoculation effects on plant growth

Effects of the soil inocula, measured across plants from all eight populations and expressed as the log ratio of biomass in inoculated soil compared to control soil (log2[inoculated soil/control soil]), ranged from −0.01 to −0.21 for the southern soils (from predominantly sexual populations) compared to +0.08 to −0.13 for the northern soils (from apomictic populations; Figure 1). Soils from the southern sites tended to cause stronger plant growth suppression than soils from the northern sites (contrast test ‘south versus north soils’, F1,110 = 3.9, p = 0.051, see Table 2). Individual soils showed large variation in effects on different plant populations, with strong growth suppression effects observed more often in soils from southern populations. A significant effect of leaf length indicated, as expected, that plants that are larger at the start of the soil inoculation treatment have higher final biomass; by including this factor as a covariate in the statistical model our test for soil inoculum effects is independent of size differences between individuals at the start of the treatment. There was no evidence that plant populations showed a different growth response to soils from their local home sites compared to soils from other sites (Figure 1, no distinct pattern in the effects observed in diagonal versus off-diagonal cells; contrast test ‘own versus foreign soils’, F8,774 = 1.5, p = 0.15).

thumbnailFigure 1. Population-level effects of field soil inocula on plant growth. Shown are log2 response ratios of shoot biomass obtained in inoculated soils relative to control soils, based on least squares means for each population x soil combination as obtained from the statistical model. Negative values indicate growth suppression compared to control soil and positive values indicate growth stimulation. Highlighted are suppressive effects (grey cells) and strong suppressive effects (dark grey cells, log ratio < −0.3). Highlighted diagonal cells represent plants growing in soil from their own home site.

Table 2. Split-plot analysis of the soil feedback experiment, testing for soil effects on plant shoot biomass

To test whether prevalence of a specialized rust fungus and the specialized Taraxacum weevil Glocianus punctiger differ along the north–south transect we compared the proportions of infected plants per population between the northern part (8 populations) and the southern part of the transect (10 populations, see Table 1). Northern populations showed lower levels of weevil infection (mean proportion of weevil-infected plants in northern populations: 0.13; in southern populations: 0.37; t-test of north–south comparison: t = −2.29, df = 10, p = 0.045). Rust infection also tended to be lower in northern populations (mean proportion of rust-infected plants in northern populations: 0.07; in southern populations: 0.26) but this difference was not statistically significant (t-test of north–south comparison: t = −1.64, df = 16, p = 0.12). Consistent with the observed north–south patterns in rust and weevil incidence, infection levels by rust and weevils tended to be lower in apomictic populations than in sexual populations; however, these relations were not statistically significant (Figure 2; logistic regression parameters and test statistics from PROC GENMOD in SAS 9.1).

thumbnailFigure 2. Population-level association between the proportion of asexual plants and infection prevalence. a. Association with incidence of rust infection (18 populations). b. Association with weevil larvae infection (12 populations). Dashed lines are logistic regressions based on population proportions. Overlapping data points at x,y = 0,0 or 0,1 are jittered for the purpose of visualization.

Discussion

The biotic interactions (RQ) hypothesis for geographic parthenogenesis assumes, first, that herbivory and pathogen pressure is less intense at higher latitudes; and second that pressure from coevolving pests and pathogens selects for host sexual reproduction while absence of this pressure selects for asexual reproduction. Our study tested and tentatively confirmed the first assumption and is further consistent with the second assumption. For all three types of antagonists that we tested, northern dandelion populations on average showed reduced antagonist prevalence or impact compared to southern dandelion populations, and this was generally correlated to dandelion sex-asex variation and geographic parthenogenesis. The statistical support for this pattern varied among antagonists, being significant for weevils, very close to significant for soil pathogens (p = 0.051) and non-significant for rust incidence. However, the consistency of the relationship in all antagonist types lends some confidence to the observed pattern. The lack of statistical significance may be a statistical power issue due to a comparatively limited number of populations included in our study, relative to large inter-population variation. For instance, similar surveys in other systems demonstrated the same pattern using larger numbers of populations [42,43]. A positive correlation between pest/pathogen prevalence and host sexual reproduction has been found in some but not in all systems that have been comprehensively studied [44,45].

Recent studies have challenged the existence of a latitudinal trend in herbivory [26,27]. Large scale studies of latitudinal trends in herbivory often compare different species and systems. It has been pointed out that better control for phylogenetic differences should provide better tests of the latitudinal trends [27,46]. Studies such as ours, which compare populations of the same species along a latitudinal transect, may provide the best possible test cases for the existence of a latitudinal trend in plant-antagonist interactions because the patterns are not confounded by species or community-specific differences along the latitudinal transect.

Our results are consistent with, but do not directly confirm, the assumption that antagonists select for host sexual reproduction. Better support for this hypothesis could come from observing the same correlation between antagonists and host sex-asex variation in randomly selected populations, not in populations sampled along the north–south transect of geographic parthenogenesis (because of possible confounding effects of other factors that might cause geographic parthenogenesis). Population surveys in different species have previously explored a correlation between host sex-asex variation and antagonist prevalence, confirming RQ predictions in some systems [42,43,47] but not in others [44,45,48]. If we look at the populations only in the southern part of the dandelion transect, where sexuals and asexuals coexist, the association between infection rate and the proportion of asexuals also showed a negative trend across populations (β = −0.16 for the asexual-rust association, n = 10 population; β = −0.26 for the asexual-weevil association, n = 7 populations). This conforms to the RQ predictions, however, the trends are not significant and the low number of populations makes it difficult to confidently interpret this result. In general, testing RQ predictions in areas of mixed sex-asex dandelion populations may be complicated by potentially high levels of clonal diversity within apomictic dandelion populations [49]. Because pests and pathogens may select for host diversity rather than host sexual reproduction per se, sexual reproduction is not necessarily the winning strategy under high pest and pathogen pressure when competing with a diverse population of asexual genotypes [50].

Our study provides a first evaluation of the patterns of incidence for some of the antagonists that are relevant in natural dandelion populations. Of the tested antagonists, the rust fungus may be the most relevant in a RQ context. Our survey shows that it is widespread and often infects a large proportion of plants within a population. Moreover, rusts can have severe fitness impacts in natural plant populations [51,52] and can track common host genotypes and potentially drive negative frequency dependent selection within populations [53,54]. Soil pathogens have not previously been considered in plant RQ studies. We show that pathogenic effects of soil communities are common and can result in dandelion growth reduction of 20% or more (log2 response ratio ≤ 0.3). This is the net effect of all soil microbes present in the soil communities. Identifying fitness effects of individual soil pathogen species may not be straightforward, but the widespread occurrence of soil pathogenic effects highlights the impact of soil antagonists for dandelion performance. We detected the weevil G. punctiger in all populations along the transect that we screened for weevil larvae presence. This widespread occurrence, combined with generally high infection rates within populations (Table 1) and high levels of seed predation within infected capitula [55], suggest that G. punctiger is a dominant herbivore in dandelion. However, it is unknown whether the dandelion-weevil interaction shows enough genotype-specificity to play a role in RQ dynamics.

In species with geographic parthenogenesis, including T. officinale, the asexual variant often has an increased ploidy level compared to the sexual variant. It has been argued that ploidy differences, not reproductive mode, can be responsible for patterns of geographic parthenogenesis [1,6]. However, where sex-asex variation exists in plants without ploidy level change, a pattern of geographic parthenogenesis can still be observed [34]. In Boechera holboellii, where polyploidy and apomixis can be uncoupled and diploid and polyploid apomicts exist, diploids extend more extensively into previously glaciated areas than polyploids [56]. Although ploidy and not reproductive mode may contribute to geographic parthenogenesis in some systems [6], such observations suggest that polyploidy per se is not a general and sufficient explanation for geographic parthenogenesis.

Conclusions

Our study suggests that dandelion populations at higher latitudes are exposed to reduced herbivores and pathogens and this pattern generally correlates with dandelion geographic parthenogenesis. This highlights the possibility that antagonists drive patterns of geographic parthenogenesis. Whether or not antagonists are in fact causal to dandelion geographic parthenogenesis depends on the validity of the RQ hypothesis for sexual reproduction, which remains to be further demonstrated in this system.

Methods

Dandelion antagonists

In order to drive RQ dynamics, it is thought that the antagonist’s effect on host plant fitness needs to be severe as well as genotype-specific, such that some host genotypes in the population are less susceptible than others. Because little a priori knowledge is available on which antagonists are most relevant to dandelion fitness, our approach was to monitor a range of different antagonists:

Soil microbial pathogens

Effects of soil microbial communities on plants can be tested experimentally using soil inoculums and soil feedback experiments. In such experiments the net effect of a ‘living’ soil (potentially containing both mutualist and antagonist microbial species) on plant growth is evaluated [57]. Many studies have indicated that soil inoculum effects, when compared to appropriate controls, are often negative, implying overall pathogenic effects of soil microbial communities [40]. Pathogen build-up in soils is plant species-specific [58-60] and plants often suffer more from the soil microbial communities that they accumulate themselves than from soil microbial communities associated with other plant species [61]. Intraspecific variation in soil feedbacks has also been demonstrated [62] which suggests that plant genotypes within a population might differ in the pathogen communities that they accumulate and/or in their vulnerability to the same pathogen communities.

Rust fungus

Rust fungi and their host plants show highly coevolved interactions based on gene-for-gene recognition mechanisms [63]. Such mechanisms are characterized by genotypic variation in host sensitivity and pathogen virulence. Although it is debated whether gene-for-gene systems lend themselves easily to RQ dynamics [64-66], rust fungi are considered as potential drivers of negative frequency-dependent selection and RQ dynamics in host plant populations [53,54]. A specialized rust fungus, Puccinia hieracii infects Taraxacum in a genotype-dependent way [67] and several other rust fungus species can also infect T. officinale[68]. Fitness impact of rust infection is unknown in T. officinale but is reported to be significant in related Asteraceae species [51,52]. In the present field study no attempt was made to identify the rust species.

Seed-eating weevil

Glocianus punctiger is a seed-eating weevil that specializes on Taraxacum[69]. It is considered a dominant pre-dispersal predator of Taraxacum seeds [55]. Up to 30% of seeds per capitulum can be lost to weevil predation [55] suggesting considerable fitness consequences of infestation.

Population survey

Field sampling

Taraxacum officinale populations were sampled in hay meadows and grazed pastures from eastern France to southern Sweden during spring 2005 (18 populations in total, Figure 3). In each population 100 plants were sampled along 5 transects of 20 m, selecting the closest individual with ripe seeds at 1-meter intervals. From each plant, one capitulum with ripe seeds was collected for subsequent greenhouse propagation and ploidy analysis (see below). Each plant was checked for visible signs of leaf rust infection, usually orange-brown uredinia pustules but also larger aecia in some populations. In 12 of the populations we additionally screened for Glocianus punctiger infestation by randomly selecting 100 developing capitula and inspecting them for the presence of weevil larvae inside. At the growing sites of four apomictic populations from the northern part of the transect (B1, N1, D2, Sw1) and four predominantly sexual populations from the southern part of the transect (F1, F3, S1, G1; see Figure 3) field soil samples were collected from the upper 15 cm soil layer, taken at several locations per field site and subsequently mixed to yield one bulk sample per site. Soil samples were stored at 4°C until experimentation.

thumbnailFigure 3. Locations of T. officinale sampling sites. Labels correspond to Table 1.

Ploidy analysis

Per population, seeds from 30 collected capitula were germinated and seedlings were grown under standard greenhouse conditions. Ploidy level was determined in leaf tissue from two seedlings from each collected capitulum, using a Partec Ploidy Analyzer ultraviolet flow cytometer to measure nuclear DNA content relative to a reference sample of known diploidy [70]. In virtually all cases the two seedlings were either both diploid or both triploid, indicating respectively sexual or apomictic status of the field-sampled mother plant [70]. In mixed populations with both diploids and triploids we screened an additional 30–70 plants to improve the estimation of the proportion of apomicts in the population.

Soil inoculation experiment

Experimental design

Plants from the four northern populations (B1, N1, D2, Sw1) and the four southern populations (F1, F3, S1, G1) were exposed to field-collected soils from all eight populations in a cross-inoculation greenhouse experiment. Seeds from 15 plants per population were germinated in petri dishes on moist filter paper and seedlings were transplanted individually into 38 ml seedling tray cells filled with sterilized mineral sand (sterilized by gamma radiation ≥ 25 kGray). Eight days after transplanting, 4 ml of field soil suspension was added to each individual plant. Soil suspensions were made by adding demi water to field-collected soil samples to a final mass ratio of ¼ dry soil and ¾ water (adjusted for soil moisture content as determined by evaporative water loss in oven-dried test samples), thorough mixing (1h at 170 rpm) and passing the mixture through a 200-μm mesh sieve. 15 plants from each of the eight populations were grown with each of the eight soils, resulting in 960 plants in the experiment. An additional 30 plants per population were grown in sterilized sand without added soil suspension (control soils), in order to evaluate the overall effect of soil addition on plant growth. Seedling tray cells were placed in plastic containers in groups of eight (one plant per population) with all plants per container receiving the same soil suspension. Containers were placed on a greenhouse bench in a randomized block design, with three blocks that each contained five replicates of each population x soil treatment (blocks corresponded to three consecutive days for setting up the entire experiment). Plants received 0.5 strength Hoagland nutrient solution 1x per week and demi water 2x per week during the experiment. Plants were harvested four weeks after transplanting. Shoot biomass was determined after oven drying at 70°C for 48 hours.

Statistical analysis

We used a split-plot design to test for soil suspension effects on shoot biomass, with ‘soil origin’ as main-plot factor (8 levels, administered at the level of containers) and ‘plant population’ as split-plot level (8 levels, at the level of individual cells within containers). Of specific interest is the a priori contrast between soils from southern (predominantly sexual) versus northern (apomictic) populations. This explicitly tests the hypothesis that geographic parthenogenesis is correlated with differences in soil pathogenicity, specifically with reduced pathogenicity in northern apomictic populations. Blocking in the experiment (3 levels) was included in the model as a main effect. The interaction between ‘soil origin’ and ‘replicate’, where ‘replicate’ is nested within blocks, was used as error term for testing the soil origin effect. This interaction captures all variability between replicated main-plots (soil origin), which in this experiment is contained both in the differences between main-plot replicates within each block and in the soil x block interaction. At the time of soil suspension administration we measured the length of the first true leaf in each plant and included this as a covariate in the model to account for size effects prior to the soil treatment. Analysis was carried out in SAS version 9.1 for Windows (The SAS Institute, Cary, NC).

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KJFV and AB conceived of the study, KJFV designed and carried out the population survey and the greenhouse soil experiment, performed the statistical analysis and drafted the manuscript. AB participated in drafting the manuscript. Both authors read and approved the final manuscript.

Acknowledgements

We thank Thure Hauser and Francois Felber for their help during the population survey, Gideon van den Ende for assistance with the soil inoculum experiment, and Jac Thissen for statistical advice. Funding for this study was provided by the Netherlands Organisation for Scientific Research (NWO-ALW). This is publication 5368 of the Netherlands Institute of Ecology (NIOO-KNAW).

References

  1. Bierzychudek P: Patterns in plant parthenogenesis.

    Experientia 1985, 41(10):1255-1264. Publisher Full Text OpenURL

  2. Buckley TR, Marske K, Attanayake D: Phylogeography and ecological niche modelling of the New Zealand stick insect Clitarchus hookeri (White) support survival in multiple coastal refugia.

    J Biogeogr 2010, 37(4):682-695. Publisher Full Text OpenURL

  3. Law JH, Crespi BJ: The evolution of geographic parthenogenesis in Timema walking-sticks.

    Mol Ecol 2002, 11(8):1471-1489. PubMed Abstract | Publisher Full Text OpenURL

  4. Martins MJF, Vandekerkhove J, Namiotko T: Environmental stability and the distribution of the sexes: insights from life history experiments with the geographic parthenogen Eucypris virens (Crustacea: Ostracoda).

    Oikos 2008, 117(6):829-836. Publisher Full Text OpenURL

  5. Thompson SL, Whitton J: Patterns of recurrent evolution and geographic parthenogenesis within apomictic polyploid Easter daises (Townsendia hookeri).

    Mol Ecol 2006, 15(11):3389-3400. PubMed Abstract | Publisher Full Text OpenURL

  6. Adolfsson S, Michalakis Y, Paczesniak D, Bode SNS, Butlin RK, Lamatsch DK, Martins MJF, Schmit O, Vandekerkhove J, Jokela J: Evaluation of elevated ploidy and asexual reproduction as alternative explanations for geographic parthenogenesis in Eucypris virens ostracods.

    Evolution 2010, 64(4):986-997. PubMed Abstract OpenURL

  7. Britton NF, Mogie M: Poor male function favours the coexistence of sexual and asexual relatives.

    Ecol Lett 2001, 4(2):116-121. Publisher Full Text OpenURL

  8. Haag CR, Ebert D: A new hypothesis to explain geographic parthenogenesis.

    Annales Zoologici Fennici 2004, 41(4):539-544. OpenURL

  9. Kearney M: Hybridization, glaciation and geographical parthenogenesis.

    Trends Ecol Evol 2005, 20(9):495-502. PubMed Abstract | Publisher Full Text OpenURL

  10. Horandl E: The complex causality of geographical parthenogenesis.

    New Phytol 2006, 171(3):525-538. PubMed Abstract OpenURL

  11. Lynch M: Destabilizing hybridization, general purpose genotypes and geographic parthenogenesis.

    Q Rev Biol 1984, 59(3):257-290. Publisher Full Text OpenURL

  12. Stebbins GL: Variation and evolution in plants. New York: Columbia University Press; 1950. OpenURL

  13. Bell G: The masterpiece of nature: the evolution and genetics of sexuality. Berkeley: California Press; 1982. OpenURL

  14. Glesener RR, Tilman D: Sexuality and components of environmental uncertainty: clues from geographic parthenogenesis in terrestrial animals.

    Amer Nat 1978, 112(986):659-673. Publisher Full Text OpenURL

  15. Jeanike J: An hypothesis to account for the maintenance of the sex within populations.

    Evol Theory 1978, 3:191-194. OpenURL

  16. Hamilton WD: Sex versus non-sex versus parasite.

    Oikos 1980, 35:282-290. Publisher Full Text OpenURL

  17. Decaestecker E, Gaba S, Raeymaekers JAM, Stoks R, Van Kerckhoven L, Ebert D, De Meester L: Host-parasite ‘Red Queen’ dynamics archived in pond sediment.

    Nature 2007, 450(7171):870-874. PubMed Abstract | Publisher Full Text OpenURL

  18. Dybdahl MF, Lively CM: Host-parasite coevolution: Evidence for rare advantage and time-lagged selection in a natural population.

    Evolution 1998, 52(4):1057-1066. Publisher Full Text OpenURL

  19. Kerstes NAG, Berenos C, Schmid-Hempel P, Wegner KM: Antagonistic experimental coevolution with a parasite increases host recombination frequency.

    Bmc Evol Bio 2012, 12:18. BioMed Central Full Text OpenURL

  20. Morran LT, Schmidt OG, Gelarden IA, Parrish RC, Lively CM: Running with the Red Queen: Host-parasite coevolution selects for biparental sex.

    Science 2011, 333(6039):216-218. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Lively CM: A review of Red Queen models for the persistence of obligate sexual reproduction.

    J Hered 2010, 101:S13-S20. PubMed Abstract | Publisher Full Text OpenURL

  22. Hartfield M, Keightley PD: Current hypotheses for the evolution of sex and recombination.

    Integrative Zoology 2012, 7(2):192-209. PubMed Abstract | Publisher Full Text OpenURL

  23. Coley PD, Barone JA: Herbivory and plant defenses in tropical forests.

    Ann Rev Ecol Evol Syst 1996, 27:305-335. Publisher Full Text OpenURL

  24. Schemske DW, Mittelbach GG, Cornell HV, Sobel JM, Roy K: Is there a latitudinal gradient in the importance of biotic interactions?

    Ann Rev Ecol Evol Syst 2009, 40:245-269. Publisher Full Text OpenURL

  25. Coley PD, Aide TM: Comparison of herbivory and plant defenses in temperate and tropical broad-leaved forests. In Plant-animal interactions: Evolutionary ecology in tropical and temperate regions. Edited by Price PW, Lewinsohn TM, Fernandes GW, Benson WW. New York: Wiley; 1991:25-49. OpenURL

  26. Moles AT, Bonser SP, Poore AGB, Wallis IR, Foley WJ: Assessing the evidence for latitudinal gradients in plant defence and herbivory.

    Funct Ecol 2011, 25(2):380-388. Publisher Full Text OpenURL

  27. Moles AT, Wallis IR, Foley WJ, Warton DI, Stegen JC, Bisigato AJ, Cella-Pizarro L, Clark CJ, Cohen PS, Cornwell WK, et al.: Putting plant resistance traits on the map: a test of the idea that plants are better defended at lower latitudes.

    New Phytol 2011, 191(3):777-788. PubMed Abstract | Publisher Full Text OpenURL

  28. Johnson MTJ, Rasmann S: The latitudinal herbivory-defence hypothesis takes a detour on the map.

    New Phytol 2011, 191(3):589-592. PubMed Abstract | Publisher Full Text OpenURL

  29. Springer YP: Edaphic quality and plant-pathogen interactions: effects of soil calcium on fungal infection of a serpentine flax.

    Ecol 2009, 90(7):1852-1862. Publisher Full Text OpenURL

  30. Albrectsen BR, Witzell J, Robinson KM, Wulff S, Luquez VMC, Agren R, Jansson S: Large scale geographic clines of parasite damage to Populus tremula L.

    Ecography 2010, 33(3):483-493. OpenURL

  31. Burdon JJ, Oates JD, Marshall DR: Interactions between Avena and Puccinia species. 1. The wild hosts: Avena barbata Pott Ex Link, A. fatua L., A. ludoviciana Durieu.

    J Appl Ecol 1983, 20(2):571-584. Publisher Full Text OpenURL

  32. Seabloom EW, Borer ET, Mitchell CE, Power AG: Viral diversity and prevalence gradients in North American Pacific Coast grasslands.

    Ecology 2010, 91(3):721-732. PubMed Abstract | Publisher Full Text OpenURL

  33. Chu HY, Fierer N, Lauber CL, Caporaso JG, Knight R, Grogan P: Soil bacterial diversity in the Arctic is not fundamentally different from that found in other biomes.

    Environ Microbiol 2010, 12(11):2998-3006. PubMed Abstract | Publisher Full Text OpenURL

  34. Johnson MTJ, Smith SD, Rausher MD: Effects of plant sex on range distributions and allocation to reproduction.

    New Phytol 2010, 186(3):769-779. PubMed Abstract | Publisher Full Text OpenURL

  35. Silvertown J: The evolutionary maintenance of sexual reproduction: Evidence from the ecological distribution of asexual reproduction in clonal plants.

    Int J Plant Sci 2008, 169(1):157-168. Publisher Full Text OpenURL

  36. Johnson MTJ, Smith SD, Rausher MD: Plant sex and the evolution of plant defenses against herbivores.

    Proc Natl Acad Sci USA 2009, 106(43):18079-18084. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Menken SBJ, Smit E, DenNijs HCM: Genetical population structure in plants: Gene flow between diploid sexual and triploid asexual dandelions (Taraxacum section Ruderalia).

    Evolution 1995, 49(6):1108-1118. Publisher Full Text OpenURL

  38. Mogie M, Ford H: Sexual and asexual Taraxacum species.

    Biol J Linn Soc 1988, 35(2):155-168. Publisher Full Text OpenURL

  39. van Dijk PJ: Ecological and evolutionary opportunities of apomixis: insights from Taraxacum and Chondrilla.

    Philos Trans R Soc Lond B Biol Sci 2003, 358(1434):1113-1120. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Kulmatiski A, Beard KH, Stevens JR, Cobbold SM: Plant-soil feedbacks: a meta-analytical review.

    Ecol Lett 2008, 11(9):980-992. PubMed Abstract | Publisher Full Text OpenURL

  41. Lively CM: Trematode infection and the distribution and dynamics of parthenogenetic snail populations.

    Parasitology 2001, 123:S19-S26. PubMed Abstract OpenURL

  42. Lively CM: Parthenogenesis in a fresh-water snail: reproductive assurance versus parasitic release.

    Evolution 1992, 46(4):907-913. Publisher Full Text OpenURL

  43. Lively CM: Evidence from a New Zealand snail for the maintenance of sex by parasitism.

    Nature 1987, 328(6130):519-521. Publisher Full Text OpenURL

  44. Ben-Ami F, Heller J: Spatial and temporal patterns of parthenogenesis and parasitism in the freshwater snail Melanoides tuberculata.

    J Evol Biol 2005, 18(1):138-146. PubMed Abstract | Publisher Full Text OpenURL

  45. Killick SC, Obbard DJ, West SA, Little TJ: Parasitism and breeding system variation in North American populations of Daphnia pulex.

    Ecol Res 2008, 23(1):235-240. Publisher Full Text OpenURL

  46. Rasmann S, Agrawal AA: Latitudinal patterns in plant defense: evolution of cardenolides, their toxicity and induction following herbivory.

    Ecol Lett 2011, 14(5):476-483. PubMed Abstract | Publisher Full Text OpenURL

  47. Kumpulainen T, Grapputo A, Mappes J: Parasites and sexual reproduction in psychid moths.

    Evolution 2004, 58(7):1511-1520. PubMed Abstract OpenURL

  48. Meirmans S, Skorping A, Loyning MK, Kirkendall LR: On the track of the Red Queen: bark beetles, their nematodes, local climate and geographic parthenogenesis.

    J Evol Biol 2006, 19(6):1939-1947. PubMed Abstract | Publisher Full Text OpenURL

  49. van der Hulst RGM, Mes THM, den Nijs JCM, Bachmann K: Amplified fragment length polymorphism (AFLP) markers reveal that population structure of triploid dandelions (Taraxacum officinale) exhibits both clonality and recombination.

    Mol Ecol 2000, 9(1):1-8. PubMed Abstract | Publisher Full Text OpenURL

  50. Lively CM, Howard RS: Selection by parasites for clonal diversity and mixed mating.

    Philos Trans R Soc Lond B Biol Sci 1994, 346(1317):271-280. PubMed Abstract | Publisher Full Text OpenURL

  51. Burdon JJ, Groves RH, Cullen JM: The impact of biological control on the distribution and abundance of Chondrilla juncea in Western Australia.

    J Appl Ecol 1981, 18(3):957-966. Publisher Full Text OpenURL

  52. Roy BA, Bierzychudek P: The potential for rust infection to cause natural selection in apomictic Arabis holboelli (Brassicaceae).

    Oecologia 1993, 95(4):533-541. OpenURL

  53. Chaboudez P, Burdon JJ: Frequency-dependent selection in a wild plant-pathogen system.

    Oecologia 1995, 102(4):490-493. Publisher Full Text OpenURL

  54. Siemens DH, Roy BA: Tests for parasite-mediated frequency-dependent selection in natural populations of an asexual plant species.

    Evol Ecol 2005, 19(4):321-338. Publisher Full Text OpenURL

  55. Honek A, Martinkova Z: Pre-dispersal predation of Taraxacum officinale (dandelion) seed.

    J Ecol 2005, 93(2):335-344. Publisher Full Text OpenURL

  56. Dobes C, Mitchell-Olds T, Koch MA: Intraspecific diversification in North American Boechera stricta (=Arabis drummondii), Boechera x divaricarpa, and Boechera holboellii (Brassicaceae) inferred from nuclear and chloroplast molecular markers - an integrative approach.

    Am J Bot 2004, 91(12):2087-2101. PubMed Abstract | Publisher Full Text OpenURL

  57. Brinkman EP, Van der Putten WH, Bakker EJ, Verhoeven KJF: Plant-soil feedback: experimental approaches, statistical analyses and ecological interpretations.

    J Ecol 2010, 98(5):1063-1073. Publisher Full Text OpenURL

  58. Bever JD: Feedback between plants and their soil communities in an old field community.

    Ecology 1994, 75(7):1965-1977. Publisher Full Text OpenURL

  59. Kardol P, Bezemer TM, van der Putten WH: Temporal variation in plant-soil feedback controls succession.

    Ecol Lett 2006, 9(9):1080-1088. PubMed Abstract | Publisher Full Text OpenURL

  60. Van der Putten WH, Van Dijk C, Peters BAM: Plant-specific soil-borne diseases contribute to succession in foredune vegetation.

    Nature 1993, 362(6415):53-56. Publisher Full Text OpenURL

  61. Klironomos JN: Feedback with soil biota contributes to plant rarity and invasiveness in communities.

    Nature 2002, 417(6884):67-70. PubMed Abstract | Publisher Full Text OpenURL

  62. Felker-Quinn E, Bailey JK, Schweitzer JA: Soil biota drive expression of genetic variation and development of population-specific feedbacks in an invasive plant.

    Ecology 2011, 92(6):1208-1214. PubMed Abstract | Publisher Full Text OpenURL

  63. Dodds P, Thrall P: Recognition events and host-pathogen co-evolution in gene-for-gene resistance to flax rust.

    Funct Plant Biol 2009, 36(5):395-408. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Clay K, Kover PX: The Red Queen Hypothesis and plant/pathogen interactions.

    Annu Rev Phytopathol 1996, 34:29-50. PubMed Abstract | Publisher Full Text OpenURL

  65. Otto SP, Nuismer SL: Species interactions and the evolution of sex.

    Science 2004, 304(5673):1018-1020. PubMed Abstract | Publisher Full Text OpenURL

  66. Agrawal A, Lively CM: Infection genetics: gene-for-gene versus matching-alleles models and all points in between.

    Evol Ecol Res 2002, 4(1):79-90. OpenURL

  67. Waters CW: The control of teliospore and urediniospore formation by experimental methods.

    Phytopathol 1928, 18:157-213. OpenURL

  68. Termorshuizen AJ, Swertz CA: Dutch rust fungi. Doorwerth: Termorshuizen; 2011. OpenURL

  69. Gumovsky AV: Parasitism of Entedon costalis (Hymenoptera: Eulophidae) in Glocianus punctiger (Coleoptera: Curculionidae): an example of intentional discovery of the parasitoid-host association.

    Zootaxa 1964, 2008:40-68. OpenURL

  70. Tas ICQ, Van Dijk PJ: Crosses between sexual and apomictic dandelions (Taraxacum). I. The inheritance of apomixis.

    Heredity 1999, 83:707-714. PubMed Abstract | Publisher Full Text OpenURL