The cost of males in Daphnia pulex

Sex is paradoxical, because asexuals should replace their sexual ancestors by avoiding the demographic cost of producing males (hereafter referred to as the cost‐of‐males). Despite the large body of theoretical and empirical work dealing with the paradox of sex, the cost‐of‐males assumption has been rarely tested. In the present study, we tested the cost‐of‐males assumption in the cladoceran Daphnia pulex. Populations of this species consist of both cyclically parthenogenetic (i.e. sexuals) and obligately parthenogenetic (i.e. asexuals) lineages. In addition, some of the asexual lineages produce only female offspring, whereas others produce functional males, which can mate with sexual females. We compared the reproductive investment of sexuals, male‐producing asexuals, and non‐male‐producing asexuals when raised separately under various environmental conditions. We also determined the outcome of competition between pair‐wise combinations of these reproductive modes. The cost of males was evident when sexual and asexual females were raised separately: sexuals produced fewer female offspring. However, there was no cost of males when reproductive modes were raised in pairs, as sexuals won the competition with asexuals. Our results directly relate to the field conditions experienced by D. pulex. Sexuals might suffer the cost of males at the beginning of the season, when resource competition is low; but when conditions deteriorate as the population approaches carrying capacity, sexuals seem to be better competitors in spite of male production.     

Can phosphorus limitation contribute to the maintenance of sex? A test of a key assumption

Why sex is so common remains unclear; what is certain is that the predominance of sex despite its profound costs means that it must confer major advantages. Here, we use elemental and nucleic acid assays to evaluate a key element of a novel, integrative hypothesis considering whether sex might be favoured because of differences in body composition between sexuals and asexuals. We found that asexual Potamopyrgus antipodarum, a New Zealand snail, have markedly higher bodily phosphorus and nucleic acid content per unit mass than sexual counterparts. These differences coincide with and are almost certainly linked to the higher ploidy of the asexuals. Our results are the first documented body composition differences between sexual and asexual organisms, and the first detected phenotypic difference between sexual and asexual P. antipodarum, an important natural model system for the study of the maintenance of sex. These findings also verify a central component of our hypothesis that competition between diploid sexuals and polyploid asexuals could be influenced by phosphorus availability.     

Deleterious mutation accumulation in asexual Timema stick insects

Sexual reproduction is extremely widespread in spite of its presumed costs relative to asexual reproduction, indicating that it must provide significant advantages. One postulated benefit of sex and recombination is that they facilitate the purging of mildly deleterious mutations, which would accumulate in asexual lineages and contribute to their short evolutionary life span. To test this prediction, we estimated the accumulation rate of coding (nonsynonymous) mutations, which are expected to be deleterious, in parts of one mitochondrial (COI) and two nuclear (Actin and Hsp70) genes in six independently derived asexual lineages and related sexual species of Timema stick insects. We found signatures of increased coding mutation accumulation in all six asexual Timema and for each of the three analyzed genes, with 3.6- to 13.4-fold higher rates in the asexuals as compared with the sexuals. In addition, because coding mutations in the asexuals often resulted in considerable hydrophobicity changes at the concerned amino acid positions, coding mutations in the asexuals are likely associated with more strongly deleterious effects than in the sexuals. Our results demonstrate that deleterious mutation accumulation can differentially affect sexual and asexual lineages and support the idea that deleterious mutation accumulation plays an important role in limiting the long-term persistence of all-female lineages.     

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.     

Asexuality Associated with Marked Genomic Expansion of Tandemly Repeated rRNA and Histone Genes

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits—mode of genomic inheritance—influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum, which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S–5.8S–28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA–histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction—with potential functional consequences related to the loss of sexual reproduction.     

Recent and ancient asexuality in Timema walkingsticks

Determining the evolutionary age of asexual lineages should help in inferring the temporal scale under which asexuality and sex evolve and assessing selective factors involved in the evolution of asexuality. We used 416 bp of the mitochondrial COI gene to infer phylogenetic relationships of virtually all known Timema walkingstick species, including extensive intraspecific sampling for all five of the asexuals and their close sexual relatives. The asexuals T. douglasi and T. shepardii were very closely related to each other and evolutionarily young (less than 0.5 million years old). For the asexuals T. monikensis and T. tahoe, evidence for antiquity was weak since only one population of each was sampled, intraspecific divergences were low, and genetic distances to related sexuals were high: maximum-likelihood molecular-clock age estimates ranged from 0.26 to 2.39 million years in T. monikensis and from 0.29-1.06 million years in T. tahoe. By contrast, T. genevieve was inferred to be an ancient asexual, with an age of 0.81 to 1.42 million years. The main correlate of the age of asexual lineages was their geographic position, with younger asexuals being found further north.     

Contrasting patterns of synonymous and nonsynonymous sequence evolution in asexual and sexual freshwater snail lineages

In asexual lineages, both synonymous and nonsynonymous sequence polymorphism may be reduced due to severe founder effects when asexual lineages originate. However, mildly deleterious (nonsynonymous) mutations may accumulate after asexual lineages are formed, because the efficiency of purifying selection is reduced even in the nonrecombining mitochondrial genome. Here we examine patterns of synonymous and nonsynonymous mitochondrial sequence polymorphism in asexual and sexual lineages of the freshwater snail Campeloma. Using clade-specific estimates, we found that synonymous sequence polymorphism was significantly reduced by 75% in asexuals relative to sexuals, whereas nonsynonymous sequence polymorphism did not differ significantly between sexuals and asexuals. Two asexual clades had high negative values for Tajima's D statistic. Coalescent simulations confirmed that various bottleneck scenarios can account for this result. We also used branch-specific estimates of the ratio of amino acid to silent substitutions, K(a)/K(s). Our study revealed that K(a)/K(s) ratios are six times higher in terminal branches of independent asexual lineages compared to sexuals. Coalescent-based reconstruction of gene networks for all sexual and asexual clades indicated that nonsynonymous mutations occurred at a higher frequency in recently derived asexual haplotypes. These findings suggest that patterns of synonymous and nonsynonymous nucleotide polymorphism in asexual snail lineages may be shaped by both severe founder effect and relaxed purifying selection.     

Novel distribution pattern between coexisting sexual and obligate asexual variants of the true estuarine macroalga Ulva prolifera

Where sexual and asexual forms coexist within a species, the asexuals are often found to be prevalent in marginal habitats. This asexual distribution pattern has received evolutionary attention linked to the paradox of sex. In many marine species, there is a distributional trend of asexual modes being more common in lower salinity waters regarded as the ecogeographic marginal, being explained by negative effects of low salinities on sexual reproductive success. However, the distribution pattern of estuarine species recently adapted to low salinity waters has remained unknown. The brackish macroalga Ulva prolifera being a major benthic component of estuarine ecosystems includes a sexual variant and obligate asexual variants by means of motile spores. We examined the sexual–asexual distribution pattern of this alga along a salinity gradient in river estuaries. Surprisingly, opposite to the distributional trend of marine organisms, the results clearly showed the persistent predominance of sexuals in the lower salinity reaches than the asexuals. In addition, a frequent alternating of dioecious gametophytes and sporophytes in the sexual population was observed, suggesting the sexual reproductive process would be robustly performed in the low salinity waters. Considering U. prolifera had evolved from the ancestral marine species to become a true estuarine species of which the core habitat is the low salinity reaches, in a broad sense its sexual–asexual distribution pattern would be involved in asexual marginal occupations of the species range previously reported in other organisms. Based on the frozen niche variation model, we can give a concise explanation for the evolutionary process of this pattern; multiple asexuals with frozen genotypic variation had arisen from sexual ancestors undergoing low salinity adaptation and share the estuarine habitat with the sexuals at present.     

LOWER MITE INFESTATIONS IN AN ASEXUAL GECKO COMPARED WITH ITS SEXUAL ANCESTORS

What advantage do sexually reproducing organisms gain from their mode of reproduction that compensates for their twofold loss in reproductive rate relative to their asexual counterparts? One version of the Red Queen hypothesis suggests that selective pressure from parasites is strongest on the most common genotype in a population, and thus genetically identical clonal lineages are more vulnerable to parasitism over time than genetically diverse sexual lineages. Our surveys of the ectoparasites of an asexual gecko and its two sexual ancestral species show that the sexuals have a higher prevalence, abundance, and mean intensity of mites than asexuals sharing the same habitat. Our experimental data indicate that in one sexual/asexual pair this pattern is at least partly attributable to higher attachment rates of mites to sexuals. Such a difference may occur as a result of exceptionally high susceptibility of the sexuals to mites because of their low genetic diversity (relative to other more-outbred sexual species) and their potentially high stress levels, or as a result of exceptionally low susceptibility of the asexuals to mites because of their high levels of heterozygosity.     

Higher rate of tissue regeneration in polyploid asexual versus diploid sexual freshwater snails

Characterizing phenotypic differences between sexual and asexual organisms is a critical step towards understanding why sexual reproduction is so common. Because asexuals are often polyploid, understanding how ploidy influences phenotype is directly relevant to the study of sex and will provide key insights into the evolution of ploidy-level variation. The well-established association between genome size and cell cycle duration, evidence for a link between genome size and tissue regeneration rate and the growing body of research showing that ploidy influences growth rate and gene expression led us to hypothesize that healing and tissue regeneration might be affected by ploidy-level variation. We evaluated this hypothesis by measuring the rate of regeneration of antenna tissue of Potamopyrgus antipodarum, a New Zealand snail characterized by frequent coexistence between diploid sexuals and polyploid asexuals. Antennae of triploid and presumptive tetraploid asexuals regenerated more rapidly than the antennae of diploid sexuals, but regeneration rate did not differ between triploids and tetraploids. These results suggest either that ploidy elevation has nonlinear positive effects on tissue regeneration and/or that factors associated directly with reproductive mode affect regeneration rate more than ploidy level. The results of this study also indicate that the lower ploidy of sexual P. antipodarum is unlikely to confer advantages associated with more rapid regeneration.     

Complex mitonuclear interactions and metabolic costs of mating in male seed beetles

The lack of evolutionary response to selection on mitochondrial genes through males predicts the evolution of nuclear genetic influence on male‐specific mitochondrial function, for example by gene duplication and evolution of sex‐specific expression of paralogs involved in metabolic pathways. Intergenomic epistasis may therefore be a prevalent feature of the genetic architecture of male‐specific organismal function. Here, we assess the role of mitonuclear genetic variation for male metabolic phenotypes [metabolic rate and respiratory quotient (RQ)] associated with ejaculate renewal, in the seed beetle Callosobruchus maculatus, by assaying lines with crossed combinations of distinct mitochondrial haplotypes and nuclear lineages. We found a significant increase in metabolic rate following mating relative to virgin males. Moreover, processes associated with ejaculate renewal showed variation in metabolic rate that was affected by mitonuclear interactions. Mitochondrial haplotype influenced mating‐related changes in RQ, but this pattern varied over time. Mitonuclear genotype and the energy spent during ejaculate production affected the weight of the ejaculate, but the strength of this effect varied across mitochondrial haplotypes showing that the genetic architecture of male‐specific reproductive function is complex. Our findings unveil hitherto underappreciated metabolic costs of mating and ejaculate renewal, and provide the first empirical demonstration of mitonuclear epistasis on male reproductive metabolic processes.     

Elevated genetic diversity of mitochondrial genes in asexual populations of Bark Lice ('Psocoptera': Echmepteryx hageni)

Asexual reproduction is commonly thought to be associated with low genetic diversity in animals. Echmepteryx hageni (Insecta: 'Psocoptera') is one of several psocopteran species that are primarily parthenogenetic, but also exists in small, isolated sexual populations. We used mitochondrial DNA sequences to investigate the population history and genealogical relationships between the sexual and asexual forms of this species. The asexual population of E. hageni exhibits extremely high mitochondrial haplotype diversity (H=0.98), whereas the sexual forms had significantly lower haplotypic diversity (H=0.25, after correcting for sample size). This diversity in asexuals represents one the greatest genetic diversities reported for asexual animals in the literature. Nucleotide diversities were also higher in asexual compared to sexual populations (π=0.0071 vs. 0.00027). Compared to other reported estimates of π in insects, asexual nucleotide diversity is high, but not remarkably elevated. Three hypotheses might explain the elevated genetic diversity of asexual populations: (i) larger effective population size, (ii) greater mutation rate or (iii) possible recent origin of sexuals. In addition, phylogeographic analysis revealed little geographic structure among asexual E. hageni, although specimens from the upper Midwest form a single clade and are genetically differentiated. The mismatch distribution and neutrality tests indicate a historical population size increase, possibly associated with expansion from glacial refugia.     

Poor male function favours the coexistence of sexual and asexual relatives

Classical models of the evolution of sex typically assume that an asexual lineage, once derived, is reproductively separate from the sexual lineage from which it was derived. However, many asexuals, including hermaphrodite plants, produce male gametes capable of fertilising the eggs of co-existing sexuals, giving rise to sexual and asexual progeny. This male function of asexuals may be poor, and it has been proposed that this could favour sexuality and adversely affect the successful establishment of asexual lineages. We show that things are more complicated than this; the effect is frequency-dependent and poor male function may sometimes favour asexuality. In a spatially distributed population of flowering plants, it can prevent the successful invasion of either reproductive mode by the other via long-range dispersal. Consequently invasions must be driven by short-range dispersal, and are therefore extremely slow. Thus poor male function favours long-term co-existence of sexuals and asexuals. When coupled with the superior ability of asexuals to colonise virgin territory after an Ice Age, it may explain current ecological distribution patterns.     

Clonal diversity driven by parasitism in a freshwater snail

One explanation for the widespread abundance of sexual reproduction is the advantage that genetically diverse sexual lineages have under strong pressure from virulent coevolving parasites. Such parasites are believed to track common asexual host genotypes, resulting in negative frequency‐dependent selection that counterbalances the population growth‐rate advantage of asexuals in comparison with sexuals. In the face of genetically diverse asexual lineages, this advantage of sexual reproduction might be eroded, and instead sexual populations would be replaced by diverse assemblages of clonal lineages. We investigated whether parasite‐mediated selection promotes clonal diversity in 22 natural populations of the freshwater snail Melanoides tuberculata. We found that infection prevalence explains the observed variation in the clonal diversity of M. tuberculata populations, whereas no such relationship was found between infection prevalence and male frequency. Clonal diversity and male frequency were independent of snail population density. Incorporating ecological factors such as presence/absence of fish, habitat geography and habitat type did not improve the predictive power of regression models. Approximately 11% of the clonal snail genotypes were shared among 2–4 populations, creating a web of 17 interconnected populations. Taken together, our study suggests that parasite‐mediated selection coupled with host dispersal ecology promotes clonal diversity. This, in return, may erode the advantage of sexual reproduction in M. tuberculata populations.     

EVALUATION OF ELEVATED PLOIDY AND ASEXUAL REPRODUCTION AS ALTERNATIVE EXPLANATIONS FOR GEOGRAPHIC PARTHENOGENESIS IN EUCYPRIS VIRENS OSTRACODS

Transitions from sexual to asexual reproduction are often coupled with elevations in ploidy. As a consequence, the importance of ploidy per se for the maintenance and spread of asexual populations is unclear. To examine the effects of ploidy and asexual reproduction as independent determinants of the success of asexual lineages, we sampled diploid sexual, diploid asexual, and triploid asexual Eucypris virens ostracods across a European wide range. Applying nuclear and mitochondrial markers, we found that E. virens consists of genetically highly differentiated diploid sexual populations, to the extent that these sexual clades could be considered as cryptic species. All sexual populations were found in southern Europe and North Africa and we found that both diploid asexual and triploid asexual lineages have originated multiple times from several sexual lineages. Therefore, the asexual lineages show a wide variety of genetic backgrounds and very strong population genetic structure across the wide geographic range. Finally, we found that triploid, but not diploid, asexual clones dominate habitats in northern Europe. The limited distribution of diploid asexual lineages, despite their shared ancestry with triploid asexual lineages, strongly suggests that the wider geographic distribution of triploids is due to elevated ploidy rather than to asexuality.     

Existence of two sexual races in the planarian species switching between asexual and sexual reproduction

In certain planarian species that are able to switch between asexual and sexual reproduction, determining whether a sexual has the ability to switch to the asexual state is problematic, which renders the definition of sexuals controversial. We experimentally show the existence of two sexual races, acquired and innate, in the planarian Dugesia ryukyuensis. Acquired sexuals used in this study were experimentally switched from asexuals. Inbreeding of acquired sexuals produced both innate sexuals and asexuals, but inbreeding of innate sexuals produced innate sexuals only and no asexuals. Acquired sexuals, but not innate sexuals, were forced to become asexuals by ablation and regeneration (asexual induction). This suggests that acquired sexuals somehow retain asexual potential, while innate sexuals do not. We also found that acquired sexuals have the potential to develop hyperplastic and supernumerary ovaries, while innate sexuals do not. In this regard, acquired sexuals were more prolific than innate sexuals. The differences between acquired and innate sexuals will provide a structure for examining the mechanism underlying asexual and sexual reproduction in planarians.     

When a clonal genome finds its way back to a sexual species: evidence from ongoing but rare introgression in the hybridogenetic water frog complex

Besides several exceptions, asexual metazoans are usually viewed as ephemeral sinks for genomes, which become ‘frozen’ in clonal lineages after their emergence from ancestral sexual species. Here, we investigated whether and at what rate the asexuals are able to introgress their genomes back into the parental sexual population, thus more or less importantly affecting the gene pools of sexual species. We focused on hybridogenetic hybrids of western Palaearctic water frogs (Pelophylax esculentus), which originate through hybridization between P. ridibundus and P. lessonae, but transmit only clonal ridibundus genome into their gametes. Although usually mating with P. lessonae, P. esculentus may upon mating with P. ridibundus or another hybrid produce sexually reproducing P. ridibundus offspring with the introgressed ex‐clonal genome. We compared the rate of nuclear amplified fragment length polymorphism (AFLP) and mitochondrial introgression in two types of populations, that is, those where P. ridibundus occurs in isolation and those where it lives with the hybridogens. Although significant differentiation (Φpt) between sexual and clonal ridibundus genomes suggested limited gene flow between sexuals and hybridogens, a non‐negligible (~5%) proportion of P. ridibundus bore introgressed mtDNA and AFLP markers. Whereas transfer of mtDNA was exclusively unidirectional, introgression of nuclear markers was bidirectional. The proportion of introgressed P. ridibundus was highest in syntopic populations with P. esculentus, proving an ongoing and site‐specific interspecific genetic transfer mediated by hybridogenetic hybrids. It turns out that asexual hybrids are not just a sink for genes of sexual species, but may significantly influence the genetic architecture of their sexual counterparts.     

Phylogenetic evidence for hybrid origins of asexual lineages in an aphid species

Understanding the mode of origin of asexuality is central to ongoing debates concerning the evolution and maintenance of sexual reproduction in eukaryotes. This is because it has profound consequences for patterns of genetic diversity and ecological adaptability of asexual lineages, hence on the outcome of competition with sexual relatives both in short and longer terms. Among the possible routes to asexuality, hybridization is a very common mechanism in animals and plants. Aphids present frequent transitions from their ancestral reproductive mode (cyclical parthenogenesis) to permanent asexuality, but the mode of origin of asexual lineages is generally not known because it has never been thoroughly investigated with appropriate molecular tools. Rhopalosiphum padi is an aphid species with coexisting sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages that are genetically distinct. Previous studies have shown that asexual lineages of R. padi are heterozygous at most nuclear loci, suggesting either that they have undergone long-term asexuality (under which heterozygosity tends to increase) or that they have hybrid origins. To discriminate between these alternatives, we conducted an extensive molecular survey combining the sequence analysis of alleles of two nuclear DNA markers and mitochondrial DNA haplotypes in sexual and asexual lineages of R. padi. Both nuclear and cytoplasmic markers clearly showed that many asexual lineages have hybrid origins, the first such demonstration in aphids. Our results also indicated that asexuals result from multiple events of hybridization between R. padi and an unknown sibling species, and are of recent origin (contradicting previous estimates that asexual R. padi lineages were of moderate longevity). This study constitutes another example that putatively ancient asexual lineages are actually of much more recent origin than previously thought. It also presents a robust approach for testing whether hybrid origin of asexuality is also a common phenomenon in aphids.     

Migration load and the coexistence of ecologically similar sexuals and asexuals

The frozen niche variation hypothesis suggests that sexuals can coexist with closely related, ecologically similar asexuals because sexuals and narrowly adapted asexual clones use different resources. However, because a collection of clones can potentially dominate the entire resource axis, such coexistence is not stable. We show that if the sexual population inhabits multiple selection regimes and asexuals are intrinsically slightly less fit than sexuals, migration load in the sexual population allows sexuals and asexuals to coexist stably at the regional level. By decreasing sexuals' fitness, migration load allows asexuals to invade the sexual population. However, as the sexuals' range contracts, migration load decreases, preventing asexuals from driving sexuals to extinction. This "buffering" effect of migration load is even more relevant in models that include more realistic conditions, such as demographic asymmetries or explicit spatial structure.