Broad-Scale Analysis Contradicts the Theory That Generation Time Affects Molecular Evolutionary Rates in Plants

Abstract Several studies of plant taxa have concluded that generation time, including annual/perennial life history, may explain molecular evolutionary rate variation in selectively neutral DNA. Unlike in animals, there is little theoretical basis for why generation-time effects would exist in plants. Furthermore, previous reports fail to establish the generality of a generation-time effect in plants because of the small size of the datasets, a large proportion of which compared very widely divergent taxa differing in many characteristics other than generation time. Using 24 phylogenetically independent species pairs, each containing a species with an annual and a species with a perennial life history, and nine species pairs, each containing a tree species with a short and a long minimum generation time, we found no evidence that generation time is related to molecular evolutionary rate variation of the nuclear 18S ITS1 and ITS2 regions. This analysis strongly contradicts the growing belief that evolutionary rates are affected by generation time in plants. Possible reasons for the absence of generation-time effects are discussed, including an evaluation of the cell-division theory.     

A comparison of life‐history and parental care in temperate and tropical wrens

Comparing closely related species that live in different environments is a powerful way to understand selective pressures that influence life‐history evolution. We examined a suite of life‐history traits and parental care in neotropical buff‐breasted wrens Cantorchilus leucotis and north‐temperate Carolina wrens Thryothorus ludovicianus (Family Troglodytidae), to test hypotheses about life‐history evolution. As expected, buff‐breasted wrens exhibited smaller clutch sizes and higher annual adult survival than Carolina wrens. We found minimal support for the nest predation hypothesis, as nest survival and age‐corrected provisioning rates to whole broods were similar between species, and number of breeding attempts and breeding season length were greater in temperate wrens. Critical predictions of the food limitation hypothesis were not supported; in particular age‐corrected provisioning rates per nestling were higher in the tropical than temperate species. The adult survival and offspring quality hypothesis garnered the most support, as buff‐breasted wrens exhibited greater age‐corrected provisioning rates per nestling, a longer nestling period, longer re‐nesting intervals following nest success, and lower annual fecundity than Carolina wrens. Despite similarly prolonged breeding seasons, reproductive strategies differ between species with buff‐breasted wrens investing considerably in single broods to optimize first‐year survival and Carolina wrens investing in multiple small broods to optimize annual fecundity.     

Veronica: Parallel morphological evolution and phylogeography in the Mediterranean

The genus Veronica s. lat. comprises about 450 species (including about 180 species from the southern hemisphere Hebe-complex), many of which grow in the Mediterranean area. Their extreme variability in morphology, life form and habitats has led to many suggestions regarding evolution and biogeography. Difficulties arise from parallel syndromes, widespread among alpine species and lowland perennials, and particularly among annual species of the genus. We have used sequences of the plastid trnL-F region and nuclear ribosomal ITS sequences to differentiate between different clades of Veronica and reveal cases of parallel evolution. Based on this data, cases of parallel evolution have been found in biogeographical patterns among the alpine species of Veroniceae, in which species from European mountains have affinities to those in the Central Asian/Himalayan region whereas alpine species from Turkey are probably more recently derived from lowland southwestern Asian taxa. Different subspecies of Veronica bombycina gained their characteristic morphology independently and parallel in adaptation to their alpine environment. Pinnatifid leaves have been gained parallel in perennial grassland species of Veronica. Finally, parallel evolutionary trends in many characters, not only morphological but also molecular characters, are common among annual species of Veronica.We wish to thank the Studienstiftung des deutschen Volkes for a doctoral scholarship to DCA and the Spanish Ministerio de Ciencia y Tecnologia for a postdoctoral grant (programa FPI) to MMO. This work was partially supported by the Junta de Castilla y Leon and the European Union (FSE) through the research project SA117/01. We also thank Prof. F. Ehrendorfer for critical comments and valuable discussions. Further, DCA thanks Daniela Hanfland and Yoshiki Nakamoto for help in the field. Manfred A. Fischer, Lena Struwe, Bernhard Dickore, Christoph Dobes, Eberhard Fischer, Niels Köster, Gerald Schneeweiss, Tod Stuessy, Tim Utteridge and Matsugo Yokota have provided plant material for this study.     

Slow algae, fast fungi: exceptionally high nucleotide substitution rate differences between lichenized fungi Omphalina and their symbiotic green algae Coccomyxa

Omphalina basidiolichens are obligate mutualistic associations of a fungus of the genus Omphalina (the exhabitant) and a unicellular green alga of the genus Coccomyxa (the inhabitant). It has been suggested that symbiotic inhabitants have a lower rate of genetic change compared to exhabitants because the latter are more exposed to abiotic environmental variation and competition from other organisms. In order to test this hypothesis we compared substitution rates in the nuclear ribosomal internal transcribed spacer region (ITS1, 5.8S, ITS2) among fungal species with rates among their respective algal symbionts. To ensure valid comparisons, only taxon pairs (12) with a common evolutionary history were used. On average, substitution rates in the ITS1 portion of Omphalina pairs were 27.5 times higher than rates in the corresponding pairs of Coccomyxa since divergence from their respective ancestor at the base of the Omphalina/Coccomyxa lineage. Substitution rates in the 5.8S and the ITS2 portions were 2.4 and 18.0 times higher, respectively. The highest rate difference (43.0) was found in the ITS1 region. These are, to our knowledge, the highest differences of substitution rates reported for symbiotic organisms. We conclude that the Omphalina model system conforms to the proposed hypothesis of lower substitution rates in the inhabitant, but that the mode of transmission of the inhabitant (vertical versus horizontal) could be a prevailing factor in the regulation of unequal rates of nucleotide substitution between co-evolving symbionts. Our phylogenetic study of Coccomyxa revealed three main lineages within this genus, corresponding to free-living Coccomyxa, individuals isolated from basidiolichens Omphalina and Coccomyxa isolated from ascolichens belonging to the Peltigerales.     

The contribution of recombination to heterozygosity differs among plant evolutionary lineages and life-forms

Background  

Despite its role as a generator of haplotypic variation, little is known about how the rates of recombination evolve across taxa. Recombination is a very labile force, susceptible to evolutionary and life trait related processes, which have also been correlated with general levels of genetic diversity. For example, in plants, it has been shown that long-lived outcrossing taxa, such as trees, have higher heterozygosity (H _e) at SSRs and allozymes than selfing or annual species. However, some of these tree taxa have surprisingly low levels of nucleotide diversity at the DNA sequence level, which points to recombination as a potential generator of genetic diversity in these organisms. In this study, we examine how genome-wide and within-gene rates of recombination evolve across plant taxa, determine whether such rates are influenced by the life-form adopted by species, and evaluate if higher genome-wide rates of recombination translate into higher H _e values, especially in trees.     

Genome size variation and evolution in Veronica

BACKGROUND AND AIMS: The amount of DNA per chromosome set is known to be a fairly constant characteristic of a species. Its interspecific variation is enormous, but the biological significance of this variation is little understood. Some of the characters believed to be correlated with DNA amount are alpine habitat, life history and breeding system. In the present study, the aim is to distinguish between direct causal connections and chance correlation of the amount of DNA in the genus Veronica. METHODS: Estimates of DNA amount were analysed for 42 members of Veroniceae in connection with results from a phylogenetic analysis of plastid trnL-F DNA sequences and tested correlations using standard statistical tests, phylogenetically independent contrasts and a model-based generalized least squares method to distinguish the phylogenetic effect on the results. KEY RESULTS: There appears to be a lower upper limit for DNA amount in annuals than in perennials. Most DNAC-values in Veroniceae are below the mean DNA C-value for annuals in angiosperms as a whole. However, the long-debated correlation of low genome size with annual life history is not significant (P = 0.12) using either standard statistical tests or independent contrasts, but it is significant with the generalized least squares method (P < 0.01). CONCLUSIONS: The correlation of annual life history and low genome size found in earlier studies could be due to the association of annual life history and selfing, which is significantly correlated with low genome size using any of the three tests applied. This correlation can be explained by models showing a reduction in transposable elements in selfers. A significant correlation of higher genome sizes with alpine habitats was also detected.     

Is genome downsizing associated with diversification in polyploid lineages of Veronica?

The study of genome size evolution in a phylogenetic context in related polyploid and diploid lineages can help us to understand the advantages and disadvantages of genome size changes and their effect on diversification. Here, we contribute 199 new DNA sequences and a nearly threefold increase in genome size estimates in polyploid and diploid Veronica (Plantaginaceae) (to 128 species, c. 30% of the genus) to provide a comprehensive baseline to explore the effect of genome size changes. We reconstructed internal transcribed spacer (ITS) and trnL‐trnL‐trnF phylogenetic trees and performed phylogenetic generalized least squares (PGLS), ancestral character state reconstruction, molecular dating and diversification analyses. Veronica 1C‐values range from 0.26 to 3.19 pg. Life history is significantly correlated with 1C‐value, whereas ploidy and chromosome number are strongly correlated with both 1C‐ and 1Cx‐values. The estimated ancestral Veronica 1Cx‐value is 0.65 pg, with significant genome downsizing in the polyploid Southern Hemisphere subgenus Pseudoveronica and two Northern Hemisphere subgenera, and significant genome upsizing in two diploid subgenera. These genomic downsizing events are accompanied by increased diversification rates, but a ‘core shift’ was only detected in the rate of subgenus Pseudoveronica. Polyploidy is important in the evolution of the genus, and a link between genome downsizing and polyploid diversification and species radiations is hypothesized. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 178 , 243–266.     

Disparate Evolution of Paralogous Introns in the Xdh Gene of Drosophila

Drosophila nuclear introns are commonly assumed to change according to a single rate of substitution, yet little is known about the evolution of these non-coding sequences. The hypothesis of a uniform substitution rate for introns seems to be at odds with recent findings that the nucleotide composition of introns varies at a scale unknown before, and that their base content variation is correlated with that of the adjacent exons. However, no direct attempt at comparing substitution rates in introns seems to have been addressed so far. We have studied the rate of nucleotide substitution over a region of the Xdh gene containing two adjacent short, constitutively spliced introns, in several species of Drosophila and related genera. The two introns differ significantly in base composition and substitution rate, with one intron evolving at least twice as fast as the other. In addition, the substitution pattern of the introns is positively associated with that of the surrounding coding regions, evidencing that the molecular evolution of these introns is impacted by the region in which they are embedded. The observed differences cannot be attributed to selection acting differently at the level of the secondary structure of the pre-mRNA. Rather, they are better accounted for by locally heterogeneous patterns of mutation. Received: 26 July 1999 / Accepted: 21 August 1999     

Physical leaf defenses – altered by Zea life‐history evolution,domestication, and breeding – mediate oviposition preference of a specialist leafhopper

Plant anti‐herbivore defenses are known to be affected by life‐history evolution, as well as by domestication and breeding in the case of crop species. A suite of plants from the maize genus Zea (Poaceae) and the specialist herbivore Dalbulus maidis (DeLong & Wolcott) (Hemiptera: Cicadellidae) were used to test the hypothesis that anti‐herbivore defenses are affected by plant life‐history evolution and human intervention through domestication and breeding for high yield. The suite of plants included a maize (Zea mays ssp. mays L.) commercial hybrid, a maize landrace, two populations of the annual Balsas teosinte (Z. mays ssp. parviglumis Iltis & Doebley), and perennial teosinte (Z. diploperennis Iltis, Doebley & Guzman). Leaf toughness, pubescence, and oviposition preference were compared among the suite of host plants looking for effects of transitions in life history (i.e., from perennial to annual life cycle), domestication (i.e., from wild annual to domesticated annual), and breeding (i.e., from landrace to hybrid maize) on defense against D. maidis. Results on leaf toughness suggested that the life‐history and domestication transitions weakened the plant's resistance to penetration by the mouthparts and ovipositor of D. maidis, whereas results on pubescence suggested that this putative defense was strengthened with the breeding transition, contrary to expectations. Results on oviposition preference of D. maidis coincided with the expectation that life‐history and domestication transitions would lead to preference for Balsas teosinte over perennial teosinte, and of landrace maize over Balsas teosinte. Also, a negative correlation suggested that oviposition preference is significantly influenced by leaf toughness. Overall, the results suggested that Zea defenses against the specialist herbivore D. maidis were variably affected by plant life‐history evolution, domestication, and breeding, and that chemical defense may play a role in Zea defense against D. maidis because leaf toughness and pubescence only partially explained its host preferences.     

Five Hundred Sixty-Five Triples of Chicken,Human, and Mouse Candidate Orthologs

The Human Genome Project has provided abundant gene sequence information on human and important model organisms. The chicken is well positioned from an evolutionary standpoint to serve as a link between higher and lower organisms, particularly mammals, and amphibia and fish. In this study we used stringent criteria to select 565 triples of chicken, human, and mouse candidate orthologs. We analyze the sequences with respect to nucleotide and amino acid similarities. This analysis also allows measurement of evolutionary distances of different proteins. We found that chicken-human and chicken-mouse sequence identities are highly correlated; similarly for chicken-human and chicken-mouse evolutionary distances. With chicken as the out-group, we found that mouse has a higher substitution rate than human, supporting the generation-time effect hypothesis. We also described the transversion bias, which is the preference for some transversions than others in nucleotide substitutions. We demonstrated that there are statistically significant properties in the differences of orthologous sequences. The differential patterns, in combination with sequence similarity analysis, may lead to the identification of genes that are very divergent from the mammalian orthologs.     

Paraphyly of Veronica (Veroniceae; Scrophulariaceae): Evidence from the Internal Transcribed Spacer (ITS) Sequences of Nuclear Ribosomal DNA

Veronica (Veroniceae; Scrophulariaceae) and segregated genera, such as Hebe from New Zealand has been debated intensively in the past. We conducted an analysis of sequence data from the internal transcribed spacer region of nuclear ribosomal DNA (ITS) to evaluate the validity of segregate genera and the monophyly of Veronica. According to the results presented here, Veronica is paraphyletic, with the Hebe complex, Synthyris, and Paederota nested within the larger Veronica clade. Pseudolysimachion is in a basal polytomy of the expanded Veronica clade in the strict consensus tree and might be nested within Veronica as well. Clades within Veronica do not correspond to sections traditionally recognized. This study provides a first estimation of the phylogeny of Veroniceae using molecular data and can serve as a starting point for future investigations of Veronica and relatives. Received 24 July 2000/ Accepted in revised form 19 October 2000     

Relative rates of nucleotide substitution at the rbcl locus of monocotyledonous plants

Summary We subjected 35 rbcL nucleotide sequences from monocotyledonous taxa to maximum likelihood relative rate tests and estimated relative differences in rates of nucleotide substitution between groups of sequences without relying on knowledge of divergence times between taxa. Rate tests revealed that there is a hierarchy of substitution rate at the rbcL locus within the monocots. Among the taxa analyzed the grasses have the most rapid substitution rate; they are followed in rate by the Orchidales, the Liliales, the Bromeliales, and the Arecales. The overall substitution rate for the rbcL locus of grasses is over 5 times the substitution rate in the rbcL of the palms. The substitution rate at the third codon positions in the rbcL of the grasses is over 8 times the third position rate in the palms. The pattern of rate variation is consistent with the generation-time-effect hypothesis. Heterogenous rates of substitution have important implications for phylogenetic reconstruction.Offprint requests to: M.T. Clegg     

Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex

Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less‐fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M. decorus, where it is strongly associated with life history differences between co‐occurring M. decorus and annual M. guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre‐date this inversion, we mapped QTLs for life history traits that differ between annual M. guttatus and a more distantly related, collinear perennial species, M. tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M. tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.     

Ascomycota has a faster evolutionary rate and higher species diversity than Basidiomycota

Differences in rates of nucleotide or amino acid substitutions among major groups of organisms are repeatedly found and well documented. A growing body of evidence suggests a link between the rate of neutral molecular change within populations and the evolution of species diversity. More than 98% of terrestrial fungi belong to the phyla Ascomycota or Basidiomycota. The former is considerably richer in number of species than the latter. We obtained DNA sequences of 21 protein-coding genes from the lichenized fungus Rhizoplaca chrysoleuca and used them together with sequences from GenBank for subsequent analyses. Three datasets were used to test rate discrepancies between Ascomycota and Basidiomycota and that within Ascomycota: (i) 13 taxa including 105 protein-coding genes, (ii) nine taxa including 21 protein-coding genes, and (iii) nuclear LSU rDNA of 299 fungal species. Based on analyses of the 105 protein-coding genes and nuclear LSU rDNA datasets, we found that the evolutionary rate was higher in Ascomycota than in Basidiomycota. The differences in substitution rates between Ascomycota and Basidiomycota were significant. Within Ascomycota, the species-rich Sordariomycetes has the fastest evolutionary rate, while Leotiomycetes has the slowest. Our results indicate that the main contribution to the higher substitution rates in Ascomycota does not come from mutualism, ecological conditions, sterility, metabolic rate or shorter generation time, but is possibly caused by the founder effect. This is another example of the correlation between species number and evolutionary rates, which is consistent with the hypothesis that the founder effect is responsible for accelerated substitution rates in diverse clades.     

Substitution Model of Sequence Evolution for the Human Immunodeficiency Virus Type 1 Subtype B gp120 Gene over the C2-V5 Region

Phylogenetic analyses frequently rely on models of sequence evolution that detail nucleotide substitution rates, nucleotide frequencies, and site-to-site rate heterogeneity. These models can influence hypothesis testing and can affect the accuracy of phylogenetic inferences. Maximum likelihood methods of simultaneously constructing phylogenetic tree topologies and estimating model parameters are computationally intensive, and are not feasible for sample sizes of 25 or greater using personal computers. Techniques that initially construct a tree topology and then use this non-maximized topology to estimate ML substitution rates, however, can quickly arrive at a model of sequence evolution. The accuracy of this two-step estimation technique was tested using simulated data sets with known model parameters. The results showed that for a star-like topology, as is often seen in human immunodeficiency virus type 1 (HIV-1) subtype B sequences, a random starting topology could produce nucleotide substitution rates that were not statistically different than the true rates. Samples were isolated from 100 HIV-1 subtype B infected individuals from the United States and a 620 nt region of the env gene was sequenced for each sample. The sequence data were used to obtain a substitution model of sequence evolution specific for HIV-1 subtype B env by estimating nucleotide substitution rates and the site-to-site heterogeneity in 100 individuals from the United States. The method of estimating the model should provide users of large data sets with a way to quickly compute a model of sequence evolution, while the nucleotide substitution model we identified should prove useful in the phylogenetic analysis of HIV-1 subtype B env sequences. Received: 4 October 2000 / Accepted: 1 March 2001     

Little to no inbreeding depression in a tapeworm with mixed mating

Meta‐studies on hermaphrodites have found a negative relationship between primary selfing rates and levels of inbreeding depression (ID) and, thus, generally support purging in inbred systems. However, in plants, high among‐taxa variance in ID results in no difference in the mean ID between outcrossing and mixed‐mating taxa. Selective interference likely explains high ID among mixed‐mating taxa, whereas low levels of ID among mixed‐mating taxa are not as stressed. Among animal hermaphrodites, primarily molluscs, there are little data on mixed‐mating systems. To fill a taxonomic and mating system gap, we tested for ID in a mixed‐mating tapeworm, Oochoristica javaensis. We provide a direct estimate of ID across infection of an intermediate host by comparing selfing rates at two life history stages. We found little to no evidence for ID, and the level of ID falls in line with what is reported for highly selfing species even though O. javaensis has mixed mating. We discuss this result within the context of kin mating in O. javaensis. Our results emphasize that primary selfing rates alone may be insufficient to classify the inbreeding history in all species when testing for a relationship to ID. Mixed‐mating taxa, and possibly some outcrossing taxa, may exhibit low levels of ID if biparental inbreeding is also driving purging. We advocate that ID studies report estimates of inbreeding history (e.g. F_IS or identity disequilibrium) from nature‐derived adult samples to provide context rather than relying on primary selfing rates alone.     

MULTIPLE ORIGINS OF SELF-COMPATIBILITY IN LINANTHUS SECTION LEPTOSIPHON (POLEMONIACEAE): PHYLOGENETIC EVIDENCE FROM INTERNAL-TRANSCRIBED-SPACER SEQUENCE DATA

Phylogenetic reconstruction based on sequence variation in the internal transcribed spacer (ITS) region of rDNA was used to investigate the evolutionary dynamics of homomorphic self-incompatibility in Linanthus section Leptosiphon (Polemoniaceae), a group of annual plant species. Hand-pollination experiments revealed that five species were self-incompatible and four were self-compatible. Optimization of breeding systems onto the tree resulting from maximum-likelihood analysis, with no assumptions made about the ancestral condition, indicated that self-incompatibility has been lost four times in this section. An alternative tree rearrangement conforming to the hypothesis of three losses of self-incompatibility did not have a significantly lower likelihood than the maximum-likelihood tree as determined by a paired-sites test, but all rearrangements resulting in fewer than three losses were statistically rejected. Linanthus bicolor, a selfing species, was found to be polyphyletic, with populations from different geographic regions occurring in three well-supported clades. Morphological similarity in these distinct lineages is likely to have resulted from convergent evolution of traits associated with self-fertilization. Selection for reproductive assurance is hypothesized to have played an important role in the recurrent transformations from self-incompatibility to selfing in this group of annual species.     

EVOLUTION OF THE MAGNITUDE AND TIMING OF INBREEDING DEPRESSION IN PLANTS

Estimates of inbreeding depression obtained from the literature were used to evaluate the association between inbreeding depression and the degree of self-fertilization in natural plant populations. Theoretical models predict that the magnitude of inbreeding depression will decrease with inbreeding as deleterious recessive alleles are expressed and purged through selection. If selection acts differentially among life history stages and deleterious effects are uncorrelated among stages, then the timing of inbreeding depression may also evolve with inbreeding. Estimates of cumulative inbreeding depression and stage-specific inbreeding depression (four stages: seed production of parent, germination, juvenile survival, and growth/reproduction) were compiled for 79 populations (using means of replicates, N = 62) comprising 54 species from 23 families of vascular plants. Where available, data on the mating system also were collected and used as a measure of inbreeding history. A significant negative correlation was found between cumulative inbreeding depression and the primary selfing rate for the combined sample of angiosperms (N = 35) and gymnosperms (N = 9); the correlation was significant for angiosperms but not gymnosperms examined separately. The average inbreeding depression in predominantly selfing species (δ = 0.23) was significantly less (43%) than that in predominantly outcrossing species (δ = 0.53). These results support the theoretical prediction that selfing reduces the magnitude of inbreeding depression. Most self-fertilizing species expressed the majority of their inbreeding depression late in the life cycle, at the stage of growth/reproduction (14 of 18 species), whereas outcrossing species expressed much of their inbreeding depression either early, at seed production (17 of 40 species), or late (19 species). For species with four life stages examined, selfing and outcrossing species differed in the magnitude of inbreeding depression at the stage of seed production (selfing δ = 0.05, N = 11; outcrossing δ = 0.32, N = 31), germination (selfing δ = 0.02, outcrossing δ = 0.12), and survival to reproduction (selfing δ = 0.04, outcrossing δ = 0.15), but not at growth and reproduction (selfing δ = 0.21, outcrossing δ = 0.27); inbreeding depression in selfers relative to outcrossers increased from early to late life stages. These results support the hypothesis that most early acting inbreeding depression is due to recessive lethals and can be purged through inbreeding, whereas much of the late-acting inbreeding depression is due to weakly deleterious mutations and is very difficult to purge, even under extreme inbreeding.     

Rates of deleterious mutation and the evolution of sex in Caenorhabditis

A variety of models propose that the accumulation of deleterious mutations plays an important role in the evolution of breeding systems. These models make predictions regarding the relative rates of protein evolution and deleterious mutation in taxa with contrasting modes of reproduction. Here we compare available coding sequences from one obligately outcrossing and two primarily selfing species of Caenorhabditis to explore the potential for mutational models to explain the evolution of breeding system in this clade. If deleterious mutations interact synergistically, the mutational deterministic hypothesis predicts that a high genomic deleterious mutation rate (U) will offset the reproductive disadvantage of outcrossing relative to asexual or selfing reproduction. Therefore, C. elegans and C. briggsae (both largely selfing) should both exhibit lower rates of deleterious mutation than the obligately outcrossing relative C. remanei. Using a comparative approach, we estimate U to be equivalent (and < 1) among all three related species. Stochastic mutational models, Muller's ratchet and Hill-Robertson interference, are expected to cause reductions in the effective population size in species that rarely outcross, thereby allowing deleterious mutations to accumulate at an elevated rate. We find only limited support for more rapid molecular evolution in selfing lineages. Overall, our analyses indicate that the evolution of breeding system in this group is unlikely to be explained solely by available mutational models.