Geographic ranges, population structure, and ages of sexual and parthenogenetic snail lineages

Abstract A sexual reproduction is thought to doom organisms to extinction due to mutation accumulation and parasite exploitation. Theoretical models suggest that parthenogens may escape the negative effects of conspecific and biological enemiecs through escape in space. Through intensive sequencing of a mitochondrial DNA (mtDNA) and a nuclear intron locus in sexual and pathenogenetic freshwater snails (Campelom), I examine three questionss: (1) Are sexual mtDNA lineage more restricted geographically than parthenogenetic mtDNA lineages? (2) Are independent pathenogenetic lineages shorter lived than sexual lineages? (3) Do pathenogens have higher intraindividual nuclear sequence diversity and form well‐differentiated monophyletic groups as expected under the Meselson effect? Geographic ranges of parthenogenetic lineages are significantly larger than geographic ranges of sexual lineages. Based on coalescence times under different deographic assumptions, asexual lineages are short lived, but there is variation in clonal ages. Although alternative explanations exit, these results suggest that asexual lineages may persist in the short term through dispersal, and that various constraints may cause geographic restriction of sexual lineagess. Both allotriploid and diploid Campleloma parthenogens have significantly higher allelic divergence within individuals, but show limited nuclear sequence divergence from sexual ancestors. In contrast to previous allozyme evidence for nonhybrid origins of diploid Campeloma parthenogens, cryptic hybridization may account for elevated heterozygosity.     

Radical amino acid mutations persist longer in the absence of sex

Harmful mutations are ubiquitous and inevitable, and the rate at which these mutations are removed from populations is a critical determinant of evolutionary fate. Closely related sexual and asexual taxa provide a particularly powerful setting to study deleterious mutation elimination because sexual reproduction should facilitate mutational clearance by reducing selective interference between sites and by allowing the production of offspring with different mutational complements than their parents. Here, we compared the rate of removal of conservative (i.e., similar biochemical properties) and radical (i.e., distinct biochemical properties) nonsynonymous mutations from mitochondrial genomes of sexual versus asexual Potamopyrgus antipodarum, a New Zealand freshwater snail characterized by coexisting and ecologically similar sexual and asexual lineages. Our analyses revealed that radical nonsynonymous mutations are cleared at higher rates than conservative changes and that sexual lineages eliminate radical changes more rapidly than asexual counterparts. These results are consistent with reduced efficacy of purifying selection in asexual lineages allowing harmful mutations to remain polymorphic longer than in sexual lineages. Together, these data illuminate some of the population‐level processes contributing to mitochondrial mutation accumulation and suggest that mutation accumulation could influence the outcome of competition between sexual and asexual lineages.     

Mutation rates and intraspecific divergence of the mitochondrial genome of Pristionchus pacificus

Evolutionary reconstruction of the natural history of an organism ultimately requires knowledge about the development, population genetics, ecology, and phylogeny of the species. Such investigations would benefit from studies of mutational processes because mutations are the source of natural variation. The nematode Pristionchus pacificus has been developed as a model organism in evolutionary biology by comparing its development with Caenorhabditis elegans. Pristionchus pacificus and related species are associated with scarab beetles, and their ecology and phylogeny are well known. More than 200 P. pacificus isolates from all over the world are available for this cosmopolitan species. We generated mutation accumulation (MA) lines in P. pacificus to study spontaneous mutation rates in the mitochondrial genome and compared mutation rate estimates with natural variation between nine representative isolates of the species. The P. pacificus mitochondrial genome is 15,955 bp in length and is typical for nematodes. Pristionchus pacificus has all known mitochondrial genes and contains an unusual suppressor transfer RNA (tRNA) for the codon UAA. This has most likely influenced the spectrum of observable mutations because 6 of 12 mutations found in the 82 MA lines analyzed are nonsense mutations that can be suppressed by the suppressor tRNA. The overall mutation rate in P. pacificus is 7.6 × 10?? per site per generation and is less than one order of magnitude different from estimates in C. elegans and Drosophila. Using this mutation rate estimate in a comparison of the mitochondrial genome of nine P. pacificus isolates, we calculate the minimum time to the most recent common ancestor at 10?-10? generations. The combination of mutation rate analysis with intraspecific divergence provides a powerful tool for the reconstruction of the natural history of P. pacificus, and we discuss the ecological implication of these findings.     

Environmental gradients structure Daphnia pulex × pulicaria clonal distribution

The rarity of eukaryotic asexual reproduction is frequently attributed to the disadvantage of reduced genetic variation relative to sexual reproduction. However, parthenogenetic lineages that evolved repeatedly from sexual ancestors can generate regional pools of phenotypically diverse clones. Various theories to explain the maintenance of this genetic diversity as a result of environmental and spatial heterogeneity [frozen niche variation (FNV), general-purpose genotype] are conceptually similar to community ecological explanations for the maintenance of regional species diversity. We employed multivariate statistics common in community ecological research to study population genetic structure in the freshwater crustacean, Daphnia pulex × pulicaria. This parthenogenetic hybrid arose repeatedly from sexual ancestors. Daphnia pulex × pulicaria populations harboured substantial genetic variation among populations and the clonal composition at each pond corresponded to nutrient levels and invertebrate predator densities. The interclonal selection process described by the FNV hypothesis likely structured our D. pulex × pulicaria populations.     

Nucleotide Polymorphism and Within-Gene Recombination in Daphnia magna and D. pulex, Two Cyclical Parthenogens

Theory predicts that partially asexual organisms may make the “best of both worlds”: for the most part, they avoid the costs of sexual reproduction, while still benefiting from an enhanced efficiency of selection compared to obligately asexual organisms. There is, however, little empirical data on partially asexual organisms to test this prediction. Here we examine patterns of nucleotide diversity at eight nuclear loci in continentwide samples of two species of cyclically parthenogenetic Daphnia to assess the effect of partial asexual reproduction on effective population size and amount of recombination. Both species have high nucleotide diversities and show abundant evidence for recombination, yielding large estimates of effective population sizes (300,000–600,000). This suggests that selection will act efficiently even on mutations with small selection coefficients. Divergence between the two species is less than one-tenth of previous estimates, which were derived using a mitochondrial molecular clock. As the two species investigated are among the most distantly related species of the genus, this suggests that the genus Daphnia may be considerably younger than previously thought. Daphnia has recently received increased attention because it is being developed as a model organism for ecological and evolutionary genomics. Our results confirm the attractiveness of Daphnia as a model organism, because the high nucleotide diversity and low linkage disequilibrium suggest that fine-scale mapping of genes affecting phenotypes through association studies should be feasible.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

The Australian scincid lizard Menetia greyii: a new instance of widespread vertebrate parthenogenesis

Molecular data derived from allozymes and mitochondrial nucleotide sequences, in combination with karyotypes, sex ratios, and inheritance data, have revealed the widespread Australian lizard Menetia greyii to be a complex of sexual and triploid unisexual taxa. Three sexual species, three presumed parthenogenetic lineages, and one animal of uncertain status were detected amongst 145 animals examined from south-central Australia, an area representing less than one-seventh of the total distribution of the complex. Parthenogenesis appears to have originated via interspecific hybridization, although presumed sexual ancestors could only be identified in two cases. The allozyme and mtDNA data reveal the presence of many distinct clones within the presumed parthenogenetic lineages. This new instance of vertebrate parthenogenesis is a first for the Scincidae and only the second definitive case of unisexuality in an indigenous Australian vertebrate.     

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.     

Mutation accumulation and the effect of copia insertions in Drosophila melanogaster

Repeated efforts to estimate the genomic deleterious mutation rate per generation (U) in Drosophila melanogaster have yielded inconsistent estimates ranging from 0.01 to nearly 1. We carried out a mutation-accumulation experiment with a cryopreserved control population in hopes of resolving some of the uncertainties raised by these estimates. Mutation accumulation (MA) was carried out by brother sister mating of 150 sublines derived from two inbred lines. Fitness was measured under conditions chosen to mimic the ancestral laboratory environment of these genotypes. We monitored the insertions of a transposable element, copia, that proved to accumulate at the unusually high rate of 0.24 per genome per generation in one of our MA lines. Mutational variance in fitness increased at a rate consistent with previous studies, yielding a mutational coefficient of variation greater than 3%. The performance of the cryopreserved control relative to the MA lines was inconsistent, so estimates of mutation rate by the Bateman-Mukai method are suspect. Taken at face value, these data suggest a modest decline in fitness of about 0.3% per generation. The element number of copia was a significant predictor of fitness within generations; on average, insertions caused a 0.76% loss in fitness, although the confidence limits on this estimate are wide.     

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.     

MULLER'S RATCHET AND MUTATIONAL MELTDOWNS

We extend our earlier work on the role of deleterious mutations in the extinction of obligately asexual populations. First, we develop analytical models for mutation accumulation that obviate the need for time-consuming computer simulations in certain ranges of the parameter space. When the number of mutations entering the population each generation is fairly high, the number of mutations per individual and the mean time to extinction can be predicted using classical approaches in quantitative genetics. However, when the mutation rate is very low, a fixation-probability approach is quite effective. Second, we show that an intermediate selection coefficient (s) minimizes the time to extinction. The critical value of s can be quite low, and we discuss the evolutionary implications of this, showing that increased sensitivity to mutation and loss of capacity for DNA repair can be selectively advantageous in asexual organisms. Finally, we consider the consequences of the mutational meltdown for the extinction of mitochondrial lineages in sexual species.     

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.     

Incongruence between morphological and mitochondrial-DNA characters suggests hybrid origins of parthenogenetic weevil lineages (Genus aramigus)

An expanded matrix of morphological characters for the genus Aramigus (Coleoptera: Curculionidae), which includes numerous polyploid parthenogenetic lineages, was compared and combined with a published matrix of mitochondrial DNA (mtDNA) characters. The matrix of morphological characters provides little resolution of the A. tessellatus and A. uruguayensis species complexes but does resolve previously unresolved relationships among other morphologically defined species (A. globoculus + A. intermedius, A. curtulus + A. planioculus). The morphological and mtDNA characters are significantly incongruent (0.435 < or = IM < or = 0.463; IMF = 0.0735), according to the tests of Farris et al. (P = 0.010) and Templeton (P < 0.005), probably because of hybrid origins of polyploid parthenogenetic lineages. For the few sexual lineages included in both matrices, morphology and mtDNA provide congruent estimates of phylogeny. In spite of recent injunctions against combining data sets that are incongruent because of differing histories, the results of the combined analyses were used to select one of the most-parsimonious mtDNA trees as the best estimate of maternal-lineage genealogy and to reconstruct the evolution of parthenogenesis under the assumption that transitions from sexuality to parthenogenesis are irreversible. Where cytogenetically justified, as in weevils, the irreversibility assumption is useful for producing conservative estimates of the age of parthenogenetic lineages in spite of potential sampling bias against sexuals.     

Origin and Genetic Diversity of Diploid Parthenogenetic Artemia in Eurasia

There is wide interest in understanding how genetic diversity is generated and maintained in parthenogenetic lineages, as it will help clarify the debate of the evolution and maintenance of sexual reproduction. There are three mechanisms that can be responsible for the generation of genetic diversity of parthenogenetic lineages: contagious parthenogenesis, repeated hybridization and microorganism infections (e.g. Wolbachia). Brine shrimps of the genus Artemia (Crustacea, Branchiopoda, Anostraca) are a good model system to investigate evolutionary transitions between reproductive systems as they include sexual species and lineages of obligate parthenogenetic populations of different ploidy level, which often co-occur. Diploid parthenogenetic lineages produce occasional fully functional rare males, interspecific hybridization is known to occur, but the mechanisms of origin of asexual lineages are not completely understood. Here we sequenced and analysed fragments of one mitochondrial and two nuclear genes from an extensive set of populations of diploid parthenogenetic Artemia and sexual species from Central and East Asia to investigate the evolutionary origin of diploid parthenogenetic Artemia, and geographic origin of the parental taxa. Our results indicate that there are at least two, possibly three independent and recent maternal origins of parthenogenetic lineages, related to A. urmiana and Artemia sp. from Kazakhstan, but that the nuclear genes are very closely related in all the sexual species and parthenogegetic lineages except for A. sinica, who presumable took no part on the origin of diploid parthenogenetic strains. Our data cannot rule out either hybridization between any of the very closely related Asiatic sexual species or rare events of contagious parthenogenesis via rare males as the contributing mechanisms to the generation of genetic diversity in diploid parthenogenetic Artemia lineages.     

The effect of spontaneous mutations on competitive ability

Understanding the impact of spontaneous mutations on fitness has many theoretical and practical applications in biology. Although mutational effects on individual morphological or life‐history characters have been measured in several classic genetic model systems, there are few estimates of the rate of decline due to mutation for complex fitness traits. Here, we estimate the effects of mutation on competitive ability, an important complex fitness trait, in a model system for ecological and evolutionary genomics, Daphnia. Competition assays were performed to compare fitness between mutation‐accumulation (MA) lines and control lines from eight different genotypes from two populations of Daphnia pulicaria after 30 and 65 generations of mutation accumulation. Our results show a fitness decline among MA lines relative to controls as expected, but highlight the influence of genomic background on this effect. In addition, in some assays, MA lines outperform controls providing insight into the frequency of beneficial mutations.     

Deleterious mutation in related species of the plant genus Amsinckia with contrasting mating systems

Theory for the evolution of modifiers of the rate of mutation suggests that a lower rate of mutation may evolve after the breakdown of mechanisms that enforce outcrossing. Mutation accumulation (MA) experiments were conducted to compare deleterious mutation parameters in two closely related species of the plant genus Amsinckia, a group that exhibits wide variation in the mating system. One of the two species studied (A. douglasiana) is predominantly outcrossed in natural populations, where as the other species (A. gloriosa) is predominantly self-pollinated. Progeny assays of flower number per plant from generation 1 lines (control) and generation 11 lines (MA treatment) were conducted in both species. Dry weight measurements of progeny from the control and MA treatment in A. douglasiana also were made. Estimation of mutation parameters was conducted using maximum likelihood under the assumption of a gamma distribution of mutational effects. The two species exhibited similar rates and effects of deleterious mutation affecting flower number. Estimates of mutation rate for dry weight in A. douglasiana are close to those for flower number. Overall, the estimates of mutation parameters observed in these species are intermediate within the range reported for fitness components in other eukaryotes. The results are discussed within the context of evolutionary change in deleterious mutation accompanying mating system evolution and with respect to previous estimates of mutation parameters based on assays of inbreeding depression and the assumption of mutation-selection equilibrium.     

Distribution of nucleotide substitutions in human mitochondrial DNA genes

To analyze the distribution pattern of nucleotide substitutions in human mitochondrial DNA (mtDNA), mutational spectra of the mitochondrial genes were reconstructed. The reconstruction procedure is based on the mutation distribution data for 47 monophyletic mtDNA clusters, to which 794 examined mtDNA sequences encoding for tRNAs, rRNAs, and mitochondrial proteins are attributed. One of specific features of mitochondrial mutational spectra revealed was homoplasy of the mutations (the mean mutation number per variable nucleotide site in the coding region varied from 1.09 to 1.43). It was established that in the mtDNA genes maximum mutational constraint fell onto the guanine bases, albeit the content of these bases in the mtDNA L-chains was minimal. Maximal bias towards parallel G to A transitions was observed for rRNA genes, with the protein- and tRNA-encoding genes ranking next. Despite the fact that the differences in the average G-nucleotides content and variability between the genes of two mtDNA segments located between the OriH and OriL were statistically significant, the results did not provide the conclusion that the G-nucleotide instability observed in the mtDNA L-spectra was determined by the mechanism of asynchronous mtDNA replication, along with the deamination of cytosines in the H-chain regions, which remained single-stranded during replication.     

Genetic variation and de novo mutations in the parthenogenetic Caucasian rock lizard Darevskia unisexualis

Unisexual all-female lizards of the genus Darevskia that are well adapted to various habitats are known to reproduce normally by true parthenogenesis. Although they consist of unisexual lineages and lack effective genetic recombination, they are characterized by some level of genetic polymorphism. To reveal the mutational contribution to overall genetic variability, the most straightforward and conclusive way is the direct detection of mutation events in pedigree genotyping. Earlier we selected from genomic library of D. unisexualis two polymorphic microsatellite containing loci Du281 and Du215. In this study, these two loci were analyzed to detect possible de novo mutations in 168 parthenogenetic offspring of 49 D. unisexualis mothers and in 147 offspring of 50 D. armeniaca mothers. No mutant alleles were detected in D. armeniaca offspring at both loci, and in D. unisexualis offspring at the Du215 locus. There were a total of seven mutational events in the germ lines of four of the 49 D. unisexualis mothers at the Du281 locus, yielding the mutation rate of 0.1428 events per germ line tissue. Sequencing of the mutant alleles has shown that most mutations occur via deletion or insertion of single microsatellite repeat being identical in all offspring of the family. This indicates that such mutations emerge at the early stages of embryogenesis. In this study we characterized single highly unstable (GATA)(n) containing locus in parthenogenetic lizard species D. unisexualis. Besides, we characterized various types of mutant alleles of this locus found in the D. unisexualis offspring of the first generation. Our data has shown that microsatellite mutations at highly unstable loci can make a significant contribution to population variability of parthenogenetic lizards.     

Analysis and implications of mutational variation

Variation from new mutations is important for several questions in quantitative genetics. Key parameters are the genomic mutation rate and the distribution of effects of mutations (DEM), which determine the amount of new quantitative variation that arises per generation from mutation (V _M ). Here, we review methods and empirical results concerning mutation accumulation (MA) experiments that have shed light on properties of mutations affecting quantitative traits. Surprisingly, most data on fitness traits from laboratory assays of MA lines indicate that the DEM is platykurtic in form (i.e., substantially less leptokurtic than an exponential distribution), and imply that most variation is produced by mutations of moderate to large effect. This finding contrasts with results from MA or mutagenesis experiments in which mutational changes to the DNA can be assayed directly, which imply that the vast majority of mutations have very small phenotypic effects, and that the distribution has a leptokurtic form. We compare these findings with recent approaches that attempt to infer the DEM for fitness based on comparing the frequency spectra of segregating nucleotide polymorphisms at putatively neutral and selected sites in population samples. When applied to data for humans and Drosophila, these analyses also indicate that the DEM is strongly leptokurtic. However, by combining the resultant estimates of parameters of the DEM with estimates of the mutation rate per nucleotide, the predicted V _M for fitness is only a tiny fraction of V _M observed in MA experiments. This discrepancy can be explained if we postulate that a few deleterious mutations of large effect contribute most of the mutational variation observed in MA experiments and that such mutations segregate at very low frequencies in natural populations, and effectively are never seen in population samples.     

POLYPHYLETIC ORIGINS OF ASEXUALITY IN DAPHNIA PULEX. I. BREEDING-SYSTEM VARIATION AND LEVELS OF CLONAL DIVERSITY

There is growing evidence that transitions from sexual to asexual reproduction are often provoked by internal genetic factors rather than extrinsic selection pressures. In the cladoceran crustacean Daphnia pulex, the shift to asexuality has been linked to sex-limited meiosis suppression. Most populations of this species reproduce by obligate parthenogenesis, but cyclically parthenogenetic populations persist in the southern portion of its range. The meiosis-suppressor model predicts that asexuality in D. pulex has polyphyletic origins and that the coexistence of cyclically parthenogenetic lines with male-producing obligately asexual clones should be unstable. For the present study, we examined the genotypic structure of D. pulex populations from a region in which there is an abrupt microgeographical shift in breeding system. Populations in Michigan largely reproduce by cyclic parthenogenesis, while those in Ontario are obligately asexual. Allozyme studies on 77 populations from this area revealed 50 obligately asexual clones, divisible into two groups: one derived from a single parent species and the other derived via interspecific hybridization. Although nearly 50% of the clones retained male production, there was, as predicted, no evidence of coexistence between cyclically parthenogenetic populations and male-producing obligately asexual clones. The survey did, however, reveal a low incidence of cyclically parthenogenetic populations in Ontario. The high genotypic diversity of these populations suggests that they are not only resistant to meiosis suppression, but able to rework genetic variation gained from asexual clones into a sexual breeding system.