Twigs on the tree of life? Neutral and selective models for integrating macroevolutionary patterns with microevolutionary processes in the analysis of asexuality

Neutral models characterize evolutionary or ecological patterns expected in the absence of specific causal processes, such as natural selection or ecological interactions. In this study, we describe and evaluate three neutral models that can, in principle, help to explain the apparent 'twigginess' of asexual lineages on phylogenetic trees without involving the negative consequences predicted for the absence of recombination and genetic exchange between individuals. Previously, such phylogenetic twiggyness of asexual lineages has been uncritically interpreted as evidence that asexuality is associated with elevated extinction rates and thus represents an evolutionary dead end. Our first model uses simple phylogenetic simulations to illustrate that, with sexual reproduction as the ancestral state, low transition rates to stable asexuality, or low rates of ascertained 'speciation' in asexuals, can generate twiggy distributions of asexuality, in the absence of high extinction rates for asexual lineages. The second model, developed by Janko et   al . (2008 ), shows that a dynamic equilibrium between origins and neutral losses of asexuals can, under some conditions, generate a relatively low mean age of asexual lineages. The third model posits that the risk of extinction for asexual lineages may be higher than that of sexuals simply because asexuals inhabit higher latitudes or altitudes, and not due to effects of their reproductive systems. Such neutral models are useful in that they allow quantitative evaluation of whether empirical data, such as phylogenetic and phylogeographic patterns of sex and asexuality, indeed support the idea that asexually reproducing lineages persist over shorter evolutionary periods than sexual lineages, due to such processes as mutation accumulation, slower rates of adaptive evolution, or relatively lower levels of genetic variability.     

Clonal turnover versus clonal decay: a null model for observed patterns of asexual longevity, diversity and distribution

Phylogenetic and phylogeographic studies suggest that a majority of asexual organisms are evolutionarily recent offshoots of extant sexual taxa and that old clonal lineages tend to be isolated from their sexual and younger asexual counterparts. These observations have often been interpreted as support for the long-term disadvantages of asexuality resulting from the mechanisms of clonal decay. Although clonal decay is likely to be an important mechanism that limits the temporal and spatial distribution of asexual lineages, we argue here that contemporary phylogenetic analyses, which are mostly restricted to simple comparisons of "recent" and "ancient" clones, need to be tested against an appropriate null model of neutrality. We use computer simulations to show that many empirical observations of the distribution of asexuality do not in fact reject a null model of the neutral turnover of clones spawned by sexual relatives. In particular, neutral clonal turnover results in qualitatively similar pattern of clonal spatial distribution and age structure, as does a process that includes clonal decay. Although there are important quantitative differences between predictions made by the two models, we show that published empirical data are still inadequate to distinguish between them. Further work on sexual-asexual complexes is therefore required before clonal turnover can be rejected as a parsimonious explanation of the spatial distribution and age structure of asexual lineages.     

The accumulation of deleterious mutations within the frozen niche variation hypothesis

The frozen niche variation hypothesis proposes that asexual clones exploit a fraction of a total resource niche available to the sexual population from which they arise. Differences in niche breadth may allow a period of coexistence between a sexual population and the faster reproducing asexual clones. Here, we model the longer term threat to the persistence of the sexual population from an accumulation of clonal diversity, balanced by the cost to the asexual population resulting from a faster rate of accumulation of deleterious mutations. We use Monte-Carlo simulations to quantify the interaction of niche breadth with accumulating deleterious mutations. These two mechanisms may act synergistically to prevent the extinction of the sexual population, given: (1) sufficient genetic variation, and consequently niche breadth, in the sexual population; (2) a relatively slow rate of accumulation of genetic diversity in the clonal population; (3) synergistic epistasis in the accumulation of deleterious mutations.     

Population genetics of complex life-cycle parasites: an illustration with trematodes

Accurate inferences on population genetics data require a sound underlying theoretical null model. Organisms alternating sexual and asexual reproduction during their life-cycle have been largely neglected in theoretical population genetic models, thus limiting the biological interpretation of population genetics parameters measured in natural populations. In this article, we derive the expectations of those parameters for the life-cycle of monoecious trematodes, a group comprising several important human and livestock parasites that obligatorily alternate sexual and asexual reproduction during their life-cycle. We model how migration rates between hosts, sexual and asexual mutation rates, adult selfing rate and the variance in reproductive success of parasites during the clonal phase affect the amount of neutral genetic diversity of the parasite (effective population size) and its apportionment within and between definitive hosts (using F-statistics). We demonstrate, in particular, that variance in reproductive success of clones, a parameter that has been completely overlooked in previous population genetics models, is very important in shaping the distribution of the genetic variability both within and among definitive hosts. Within definitive hosts, the parameter F(IS) (a measure of the deviation from random mating) is decreased by high variance in clonal reproductive success of larvae but increased by high adult self-fertilisation rates. Both clonal multiplication and selfing have similar effects on between-host genetic differentiation (F(ST)). Migration occurring before and after asexual reproduction can have different effects on the patterns of F(IS), depending on values of the other parameters such as the mutation rate. While the model applies to any hermaphroditic organism alternating sexual and clonal reproduction (e.g. many plants), the results are specifically discussed in the light of the limited population genetic data on monoecious trematodes available to date and their previous interpretation. We hope that our model will encourage more empirical population genetics studies on monoecious trematodes and other organisms with similar life-cycles.     

On the fate of sexual traits under asexuality

Environmental shifts and life‐history changes may result in formerly adaptive traits becoming non‐functional or maladaptive. In the absence of pleiotropy and other constraints, such traits may decay as a consequence of neutral mutation accumulation or selective processes, highlighting the importance of natural selection for adaptations. A suite of traits are expected to lose their adaptive function in asexual organisms derived from sexual ancestors, and the many independent transitions to asexuality allow for comparative studies of parallel trait maintenance versus decay. In addition, because certain traits, notably male‐specific traits, are usually not exposed to selection under asexuality, their decay would have to occur as a consequence of drift. Selective processes could drive the decay of traits associated with costs, which may be the case for the majority of sexual traits expressed in females. We review the fate of male and female sexual traits in 93 animal lineages characterized by asexual reproduction, covering a broad taxon range including molluscs, arachnids, diplopods, crustaceans and eleven different hexapod orders. Many asexual lineages are still able occasionally to produce males. These asexually produced males are often largely or even fully functional, revealing that major developmental pathways can remain quiescent and functional over extended time periods. By contrast, for asexual females, there is a parallel and rapid decay of sexual traits, especially of traits related to mate attraction and location, as expected given the considerable costs often associated with the expression of these traits. The level of decay of female sexual traits, in addition to asexual females being unable to fertilize their eggs, would severely impede reversals to sexual reproduction, even in recently derived asexual lineages. More generally, the parallel maintenance versus decay of different trait types across diverse asexual lineages suggests that neutral traits display little or no decay even after extended periods under relaxed selection, while extensive decay for selected traits occurs extremely quickly. These patterns also highlight that adaptations can fix rapidly in natural populations of asexual organisms, in spite of their mode of reproduction.     

Parasite resistance predicts fitness better than fecundity in a natural population of the freshwater snail Potamopyrgus antipodarum

The cost of males should give asexual females an advantage when in competition with sexual females. In addition, high‐fecundity asexual genotypes should have an advantage over low‐fecundity clones, leading to reduction in clonal diversity over time. To evaluate fitness components in a natural population, we measured the annual reproductive rate of individual sexual and asexual female Potamopyrgus antipodarum, a New Zealand freshwater snail, in field enclosures that excluded competitors and predators. We used allozyme genotyping to assign the asexual females to particular clonal genotypes. We found that the most fecund asexual clones had similar or higher fecundity as the top 10% of sexual families, suggesting that fecundity selection, even without the cost of males, would lead to replacement of the sexual population by clones. Consequently, we expected that the clones with the highest fecundity would dominate the natural population. Counter to this prediction, we found that high annual reproductive rates did not correlate with the frequency of clones in the natural population. When we exposed the same clones to parasites in the laboratory, we found that resistance to infection was positively correlated with the frequency of clones in the population. The correlation between fecundity and parasite resistance was negative, suggesting a trade‐off between these two traits. Our results thus suggest that parasite resistance is an important short‐term predictor of the success of asexual P. antipodarum in this population.     

Genetic variation in sexual and clonal lineages of a freshwater snail

Sexual reproduction within natural populations of most plants and animals continues to remain an enigma in evolutionary biology. That the enigma persists is not for lack of testable hypotheses but rather because of the lack of suitable study systems in which sexual and asexual females coexist. Here we review our studies on one such organism, the freshwater snail Potamopyrgus antipodarum (Gray). We also present new data that bear on hypotheses for the maintenance of sex and its relationship to clonal diversity. We have found that sexual populations of the snail are composed of diploid females and males, while clonal populations are composed of a high diversity of triploid apomictic females. Sexual and asexual individuals coexist in stable frequencies in many ‘mixed’ populations; genetic data indicate that clones from these mixed populations originated from the local population of sexual individuals without interspecific hybridization. Field data show that clonal and sexual snails have completely overlapping life histories, but individual clonal genotypes are less variable than individuals from the sympatric sexual population. Field data also show segregation of clones among depth‐specific habitat zones within a lake, but clonal diversity remains high even within habitats. A new laboratory experiment revealed extensive clonal variation in reproductive rate, a result which suggests that clonal diversity would be low in nature without some form of frequency‐dependent selection. New results from a long‐term field study of a natural, asexual population reveal that clonal diversity remained nearly constant over a 10‐year period. Nonetheless, clonal turnover occurs, and it occurs in a manner that is consistent with parasite‐mediated, frequency‐dependent selection. Reciprocal cross‐infection experiments have further shown that parasites are more infective to sympatric host snails than to allopatric snails, and that they are also more infective to common clones than rare clones within asexual host populations. Hence we suggest that sexual reproduction in these snails may be maintained, at least in part, by locally adapted parasites. Parasite‐mediated selection possibly also contributes to the maintenance of local clonal diversity within habitats, while clonal selection may be responsible for the distribution of clones among habitats. © 2003 The Linnean Society of London. Biological Journal of the Linnean Society 2003, 79 , 165–181.     

Do clones degenerate over time? Explaining the genetic variability of asexuals through population genetic models

Background  

Quest for understanding the nature of mechanisms governing the life span of clonal organisms lasts for several decades. Phylogenetic evidence for recent origins of most clones is usually interpreted as proof that clones suffer from gradual age-dependent fitness decay (e.g. Muller's ratchet). However, we have shown that a neutral drift can also qualitatively explain the observed distribution of clonal ages. This finding was followed by several attempts to distinguish the effects of neutral and non-neutral processes. Most recently, Neiman et al. 2009 (Ann N Y Acad Sci.:1168:185-200.) reviewed the distribution of asexual lineage ages estimated from a diverse array of taxa and concluded that neutral processes alone may not explain the observed data. Moreover, the authors inferred that similar types of mechanisms determine maximum asexual lineage ages in all asexual taxa. In this paper we review recent methods for distinguishing the effects of neutral and non-neutral processes and point at methodological problems related with them.     

Neutral processes forming large clones during colonization of new areas

In species reproducing both sexually and asexually clones are often more common in recently established populations. Earlier studies have suggested that this pattern arises due to natural selection favouring generally or locally successful genotypes in new environments. Alternatively, as we show here, this pattern may result from neutral processes during species’ range expansions. We model a dioecious species expanding into a new area in which all individuals are capable of both sexual and asexual reproduction, and all individuals have equal survival rates and dispersal distances. Even under conditions that favour sexual recruitment in the long run, colonization starts with an asexual wave. After colonization is completed, a sexual wave erodes clonal dominance. If individuals reproduce more than one season, and with only local dispersal, a few large clones typically dominate for thousands of reproductive seasons. Adding occasional long‐distance dispersal, more dominant clones emerge, but they persist for a shorter period of time. The general mechanism involved is simple: edge effects at the expansion front favour asexual (uniparental) recruitment where potential mates are rare. Specifically, our model shows that neutral processes (with respect to genotype fitness) during the population expansion, such as random dispersal and demographic stochasticity, produce genotype patterns that differ from the patterns arising in a selection model. The comparison with empirical data from a post‐glacially established seaweed species (Fucus radicans) shows that in this case, a neutral mechanism is strongly supported.     

Aging asexual lineages and the evolutionary maintenance of sex

Finite populations of asexual and highly selfing species suffer from a reduced efficacy of selection. Such populations are thought to decline in fitness over time due to accumulating slightly deleterious mutations or failing to adapt to changing conditions. These within‐population processes that lead nonrecombining species to extinction may help maintain sex and outcrossing through species level selection. Although inefficient selection is proposed to elevate extinction rates over time, previous models of species selection for sex assumed constant diversification rates. For sex to persist, classic models require that asexual species diversify at rates lower than sexual species; the validity of this requirement is questionable, both conceptually and empirically. We extend past models by allowing asexual lineages to decline in diversification rates as they age, that is nonrecombining lineages “senesce” in diversification rates. At equilibrium, senescing diversification rates maintain sex even when asexual lineages, at young ages, diversify faster than their sexual progenitors. In such cases, the age distribution of asexual lineages contains a peak at intermediate values rather than showing the exponential decline predicted by the classic model. Coexistence requires only that the average rate of diversification in asexuals be lower than that of sexuals.     

Dynamic Formation of Asexual Diploid and Polyploid Lineages: Multilocus Analysis of Cobitis Reveals the Mechanisms Maintaining the Diversity of Clones

Given the hybrid genomic constitutions and increased ploidy of many asexual animals, the identification of processes governing the origin and maintenance of clonal diversity provides useful information about the evolutionary consequences of interspecific hybridization, asexuality and polyploidy. In order to understand the processes driving observed diversity of biotypes and clones in the Cobitis taenia hybrid complex, we performed fine-scale genetic analysis of Central European hybrid zone between two sexual species using microsatellite genotyping and mtDNA sequencing. We found that the hybrid zone is populated by an assemblage of clonally (gynogenetically) reproducing di-, tri- and tetraploid hybrid lineages and that successful clones, which are able of spatial expansion, recruit from two ploidy levels, i.e. diploid and triploid. We further compared the distribution of observed estimates of clonal ages to theoretical distributions simulated under various assumptions and showed that new clones are most likely continuously recruited from ancestral populations. This suggests that the clonal diversity is maintained by dynamic equilibrium between origination and extinction of clonal lineages. On the other hand, an interclonal selection is implied by nonrandom spatial distribution of individual clones with respect to the coexisting sexual species. Importantly, there was no evidence for sexually reproducing hybrids or clonally reproducing non-hybrid forms. Together with previous successful laboratory synthesis of clonal Cobitis hybrids, our data thus provide the most compelling evidence that 1) the origin of asexuality is causally linked to interspecific hybridization; 2) successful establishment of clones is not restricted to one specific ploidy level and 3) the initiation of clonality and polyploidy may be dynamic and continuous in asexual complexes.     

Neutral and selection-driven decay of sexual traits in asexual stick insects

Environmental shifts and lifestyle changes may result in formerly adaptive traits becoming non-functional or maladaptive. The subsequent decay of such traits highlights the importance of natural selection for adaptations, yet its causes have rarely been investigated. To study the fate of formerly adaptive traits after lifestyle changes, we evaluated sexual traits in five independently derived asexual lineages, including traits that are specific to males and therefore not exposed to selection. At least four of the asexual lineages retained the capacity to produce males that display normal courtship behaviours and are able to fertilize eggs of females from related sexual species. The maintenance of male traits may stem from pleiotropy, or from these traits only regressing via drift, which may require millions of years to generate phenotypic effects. By contrast, we found parallel decay of sexual traits in females. Asexual females produced altered airborne and contact signals, had modified sperm storage organs, and lost the ability to fertilize their eggs, impeding reversals to sexual reproduction. Female sexual traits were decayed even in recently derived asexuals, suggesting that trait changes following the evolution of asexuality, when they occur, proceed rapidly and are driven by selective processes rather than drift.     

INBREEDING DEPRESSION UNDER JOINT SELFING,OUTCROSSING, AND ASEXUALITY

Partial asexual reproduction was introduced into a model of inbreeding depression due to nearly recessive lethal mutations in a partially selfing population. The frequencies of asexuality, selfing, and outcrossing were either constant or occurred in cycles of a single sexual generation followed by one or more asexual generations. We found that increasing the degree of asexuality generally increases the inbreeding depression maintained in an equilibrium population with a given selfing rate. This is due to the increase in the number of mutations relative to sexual generations during which selfing-induced purging of mutations may take place. For very high genomic mutation rates, sufficient to produce a threshold rate of self-fertilization for purging recessive lethal mutations, asexuality can have the opposite effect, decreasing equilibrium inbreeding depression, because of an increase in the efficiency of selection against mutations in heterozygotes with asexuality.     

Male offspring production by asexual Potamopyrgus antipodarum, a New Zealand snail

As only females contribute directly to population growth, sexual females investing equally in sons and daughters experience a two-fold cost relative to asexuals producing only daughters. Typically, researchers have focused on benefits of sex that can counter this 'cost of males' and thus explain its predominance. Here, we instead ask whether asexuals might also pay a cost of males by quantifying the rate of son production in 45 experimental populations ('lineages') founded by obligately asexual female Potamopyrgus antipodarum. This New Zealand snail is a powerful model for studying sex because phenotypically similar sexual and asexual forms often coexist, allowing direct comparisons between sexuals and asexuals. After 2 years of culture, 23 of the 45 lineages had produced males, demonstrating that asexual P. antipodarum can make sons. We used maximum-likelihood analysis of a model of male production in which only some lineages can produce males to estimate that ~50% of lineages have the ability to produce males and that ~5% of the offspring of male-producing lineages are male. Lineages producing males in the first year of the experiment were more likely to make males in the second, suggesting that some asexual lineages might pay a cost of males relative to other asexual lineages. Finally, we used a simple deterministic model of population dynamics to evaluate how male production affects the rate of invasion of an asexual lineage into a sexual population, and found that the estimated rate of male production by asexual P. antipodarum is too low to influence invasion dynamics.     

No evidence for accumulation of deleterious mutations and fitness degradation in clonal fish hybrids: Abandoning sex without regrets

Despite its inherent costs, sexual reproduction is ubiquitous in nature, and the mechanisms to protect it from a competitive displacement by asexuality remain unclear. Popular mutation‐based explanations, like the Muller's ratchet and the Kondrashov's hatchet, assume that purifying selection may not halt the accumulation of deleterious mutations in the nonrecombining genomes, ultimately leading to their degeneration. However, empirical evidence is scarce and it remains particularly unclear whether mutational degradation proceeds fast enough to ensure the decay of clonal organisms and to prevent them from outcompeting their sexual counterparts. To test this hypothesis, we jointly analysed the exome sequences and the fitness‐related phenotypic traits of the sexually reproducing fish species and their clonal hybrids, whose evolutionary ages ranged from F1 generations to 300 ky. As expected, mutations tended to accumulate in the clonal genomes in a time‐dependent manner. However, contrary to the predictions, we found no trend towards increased nonsynonymity of mutations acquired by clones, nor higher radicality of their amino acid substitutions. Moreover, there was no evidence for fitness degeneration in the old clones compared with that in the younger ones. In summary, although an efficacy of purifying selection may still be reduced in the asexual genomes, our data indicate that its efficiency is not drastically decreased. Even the oldest investigated clone was found to be too young to suffer fitness consequences from a mutation accumulation. This suggests that mechanisms other than mutation accumulation may be needed to explain the competitive advantage of sex in the short term.     

Faster clonal turnover in high‐infection habitats provides evidence for parasite‐mediated selection

According to the Red Queen hypothesis for sex, parasite‐mediated selection against common clones counterbalances the reproductive advantage of asexual lineages, which would otherwise outcompete sexual conspecifics. Such selection on the clonal population is expected to lead to a faster clonal turnover in habitats where selection by parasites is stronger. We tested this prediction by comparing the genetic structure of clonal and sexual populations of freshwater snail Potamopyrgus antipodarum between years 2003 and 2007 in three depth‐specific habitats in Lake Alexandrina (South Island, New Zealand). These habitats differ in the risk of infection by castrating trematodes and in the relative proportion of sexual individuals. As predicted, we found that the clonal structure changed significantly in shallow and mid‐water habitats, where prevalence of infection was high, but not in the deep habitat, where parasite prevalence was low. Additionally, we found that both clonal diversity and evenness of the asexual population declined in the shallow habitat. In contrast, the genetic structure (based on F–statistics) of the coexisting sexual population did not change, which suggests that the change in the clonal structure cannot be related to genetic changes in the sexual population. Finally, the frequency of sexuals had no effect on the diversity of the sympatric clonal population. Taken together, our results show a more rapid clonal turnover in high‐infection habitats, which gives support for the Red Queen hypothesis for sex.     

EVOLUTIONARY FEEDBACKS BETWEEN REPRODUCTIVE MODE AND MUTATION RATE EXACERBATE THE PARADOX OF SEX

Evolutionary theory suggests that low mutation rates should favor the persistence of asexuals. Additionally, given the observation that most nonneutral mutations are deleterious, asexuality may strengthen selection for reduced mutation rates. This reciprocal relationship raises the possibility of a positive feedback loop between sex and mutation rate. We explored the consequences of this evolutionary feedback with an individual‐based model in which a sexual population is continually challenged by the introduction of asexual clones. We found that asexuals were more likely to spread in a population when mutation rates were able to evolve relative to a model in which mutation rates were held constant. In fact, under evolving mutation rates, asexuals were able to spread to fixation even when sexuals faced no cost of sex whatsoever. The added success of asexuals was the result of their ability to evolve lower mutation rates and thereby slow the process of mutation accumulation that otherwise limited their spread. Given the existence of ample mutation rate variation in natural populations, our findings show that the evolutionary feedback between sex and mutation rate may intensify the “paradox of sex,” supporting the argument that deleterious mutation accumulation alone is likely insufficient to overcome the reproductive advantage of asexual competitors in the short term.     

Population genetic consequences of the reproductive system in the liverwort <Emphasis Type="Italic">Mannia fragrans</Emphasis>

Ecological factors affecting reproduction and dispersal are particularly important in determining genetic structure of plant populations. Polyoicous reproductive system is not rare in bryophytes; however, to date, nothing is known about its functioning and possible population genetic effects. Using the liverwort Mannia fragrans as a model species, the main aims of this study were to separate the relative importance of the components of the polyoicous reproductive system and to assess its consequences on the genetic structure of populations. High sex expression rates increasing with patch size and strongly female-biased sex ratios were detected. Additional input into clonal growth after production of sex organs was found in males compared to females. Similar clonal traits of the rare bisexual and asexual plants and preference toward newly colonized patches suggest that selection prefers colonizers that first develop organs of both sexes, hence ensuring sexual reproduction when no partner is present. Despite frequent spore production, ISSR markers revealed low genetic diversity, probably resulting from the effective clonal propagation of the species and frequent crossing between genetically identical plants. The presence of numerous rare alleles and unique recombinant haplotypes indicates occasional recombination and mutation. Effective spreading of new haplotypes is probably hampered however by large spore size. Since populations are small and isolated, such haplotypes are probably continuously eliminated by genetic drift. These results suggest that although both sexual and asexual reproductions seem to be effective, asexual components of the reproductive system play a greater role in shaping the genetic composition of the populations.     

A population of sexual Daphnia pulex resists invasion by asexual clones

Asexual reproduction avoids the costs associated with sex, predicting that invading asexual clones can quickly replace sexual populations. Daphnia pulex populations in the Great Lakes area are predominately asexual, but the elimination of sexual populations by invading clones is poorly understood. Asexual clones were detected at low frequency in one rare sexual population in 1995, with some increase in frequency during 2003 and 2004. However, these clones remained at low frequency during further yearly sampling (2005–2013) with no evidence that the resident sexual population was in danger of elimination. There was evidence for hybridization between rare males produced by asexual clones and sexual females with the potential to produce new asexual genotypes and spread the genetic factors for asexuality. In a short-term laboratory competition experiment, the two most common asexual clones did not increase in frequency relative to a genetically diverse sexual population due in part to a greater investment in diapausing eggs that trades-off current population growth for increased contribution to the egg bank. Our results suggest that a successful invasion can be prolonged, requiring a combination of clonal genotypes with high fitness, persistence of clones in the egg bank and negative factors affecting the sexual population such as inbreeding depression resulting from population bottlenecks.     

Sperm-dependent parthenogens delay the spatial expansion of their sexual hosts

It has been known for long time that asexual organisms may affect the distribution of sexual taxa. In fact, such phenomenon is inherent in the concept of geographical parthenogenesis. On the other hand, it was generally hypothesized that sperm-dependent asexuals may not exercise the same effect on related sexual population, due to their dependence upon them as sperm-donors. Recently, however, it became clear that sperm-dependent asexuals may directly or indirectly affect the distribution of their sperm-hosts, but rather in a small scale. No study addressed the large-scale biogeographic effect of the coexistence of such asexuals with the sexual species. In our study we were interested in the effect of sexual–asexual coexistence on the speed of spatial expansion of the whole complex. We expand previously published Lotka–Volterra model of the coexistence of sexual and gynogenetic forms of spined loach (Cobitis; Teleostei) hybrid complex by diffusion. We show that presence of sperm-dependent parthenogens is likely to negatively affect the spatial expansion of sexuals, and hence the whole complex, compared to pure sexual population. Given that most of the known sperm-dependent asexual complexes are distributed in areas prone to climate-induced colonization/extinction events, we conclude that such mechanism may be an important agent in determining the biogeography of sexual taxa and therefore requires further attention including empirical tests.