Starting with “mitochondrial Eve” in 1987, genetics has played an increasingly important role in studies of the last two million years of human evolution. It initially appeared that genetic data resolved the basic models of recent human evolution in favor of the “out-of-Africa replacement” hypothesis in which anatomically modern humans evolved in Africa about 150,000 years ago, started to spread throughout the world about 100,000 years ago, and subsequently drove to complete genetic extinction (replacement) all other human populations in Eurasia. Unfortunately, many of the genetic studies on recent human evolution have suffered from scientific flaws, including misrepresenting the models of recent human evolution, focusing upon hypothesis compatibility rather than hypothesis testing, committing the ecological fallacy, and failing to consider a broader array of alternative hypotheses. Once these flaws are corrected, there is actually little genetic support for the out-of-Africa replacement hypothesis. Indeed, when genetic data are used in a hypothesis-testing framework, the out-of-Africa replacement hypothesis is strongly rejected. The model of recent human evolution that emerges from a statistical hypothesis-testing framework does not correspond to any of the traditional models of human evolution, but it is compatible with fossil and archaeological data. These studies also reveal that any one gene or DNA region captures only a small part of human evolutionary history, so multilocus studies are essential. As more and more loci became available, genetics will undoubtedly offer additional insights and resolutions of human evolution.

Phenotypic plasticity describes an organism's ability to produce multiple phenotypes in direct response to its environmental conditions. Over the past 15 years empiricists have found that this plasticity frequently exhibits geographic variation and often possesses a significant heritable genetic basis. However, few studies have examined both of these aspects of plasticity simultaneously. Here, we examined both the geographic and genetic variations of the plasticity for diapause incidence (the proportion of eggs that enter an arrested state of development capable of surviving over the winter) relative to temperatures and photoperiods associated with long and short season environments across six populations of the striped ground cricket, Allonemobius socius, using a half-sibling split brood quantitative genetic design. We found that plasticity, as measured by the slope of the reaction norm, was greater in the southern-low altitude region (where populations are bivoltine) relative to the southern-high and northern-low altitude regions (where populations are univoltine). However, the heritability of plasticity was only significantly different from zero in univoltine populations that experienced “intermediate” natal season lengths. These patterns suggest that selection may favor the plasticity of diapause incidence in bivoltine regions, but act against plasticity in regions in which populations are univoltine. Furthermore, our data suggest that under “intermediate” natal season length conditions, the interplay between local adaptation and gene flow may keep the plasticity of diapause incidence low (but still significant) while maintaining its genetic variation. As such, this study not only provides a novel observation into the geographic variation of phenotypic plasticity, but also provides much needed groundwork for tests of its adaptive significance.

A population in which there is stabilizing selection acting on quantitative traits toward an intermediate optimum becomes monomorphic in the absence of mutation. Further, genotypes that show least environmental variation are also favored, such that selection is likely to reduce both genetic and environmental components of phenotypic variance. In contrast, intraspecific competition for resources is more severe between phenotypically similar individuals, such that those deviating from prevailing phenotypes have a selective advantage. It has been shown previously that polymorphism and phenotypic variance can be maintained if competition between individuals is “effectively” stronger than stabilizing selection. Environmental variance is generally observed in quantitative traits, so mechanisms to explain its maintenance are sought, but the impact of competition on its magnitude has not previously been studied. Here we assume that a quantitative trait is subject to selection for an optimal value and to selection due to competition. Further, we assume that both the mean and variance of the phenotypic value depend on genotype, such that both may be affected by selection. Theoretical analysis and numerical simulations reveal that environmental variance can be maintained only when the genetic variance (in mean phenotypic value) is constrained to a very low level. Environmental variance will be replaced entirely by genotypic variance if a range of genotypes that vary widely in mean phenotype are present or become so by mutation. The distribution of mean phenotypic values is discrete when competition is strong relative to stabilizing selection; but more genotypes segregate and the distribution can approach continuity as competition becomes extremely strong. If the magnitude of the environmental variance is not under genetic control, there is a complementary relationship between the levels of environmental and genetic variance such that the level of phenotypic variance is little affected.

Describing natural selection on phenotypic traits under varying environmental conditions is essential for a quantitative assessment of the scale at which adaptation might occur and of the impact of environmental variability on evolution. Here we analyzed patterns of multivariate selection via fecundity and viability on three reproductive traits (laying date, clutch size, and egg weight) in a population of great tits (Parus major). We quantified selection under different environmental conditions using (1) local variation in breeding density and (2) distinct areas of the population's habitat. We found that selection gradients were generally stronger for fecundity than for viability selection. We also found correlational selection acting on the combination of laying date and clutch size; this is the first documented evidence of such selection acting on these two traits in a passerine bird. Our analyses showed that both local breeding density and habitat significantly influenced selection patterns, hence favoring different patterns of reproductive investment at a small-scale relative to typical dispersal distances in this species. Canonical rotation of the nonlinear selection matrices yielded similar conclusions as traditional nonlinear selection analyses, and also showed that the main axes of selection and fitness surfaces varied over space within the population. Our results emphasize the importance of quantifying different forms of selection, and of including variation in environmental conditions at small scales to gain a better understanding of potential evolutionary dynamics in wild populations. This study suggests that the fitness landscape for this species is relatively rugged at scales relevant to the life histories of individual birds and their close relatives.

Trait variation and covariation are understood to influence the response of populations to natural selection on generational time scales, but their role, if any, in shaping long-term macroevolutionary divergence is still unclear. The present study uses the rich fossil record of the ostracode genus Poseidonamicus to reconstruct in great detail the evolutionary history of a set of landmark-based morphometric characters. This reconstruction included two kinds of evolutionary inferences: ancestor–descendant transitions among populations repeatedly sampled at the same location and divergence between lineages measured as independent contrasts on a phylogeny. This reconstructed history was then used to test if evolutionary changes were concentrated in directions (traits or combinations of traits) with high phenotypic variance. Two different statistics of association between evolution and variation tested the null hypothesis that evolutionary changes occur in random directions with respect to trait variability. The first of these measured the similarity between the directions of evolutionary change and the axis of maximum variance, and the second measured the degree to which evolutionary changes were concentrated in directions of high phenotypic variation. Randomization tests indicated that both kinds of evolutionary inferences (ancestor–descendant and phylogenetic contrasts) occurred preferentially in directions of high phenotypic variance (and close to the axis of maximal variation), suggesting that within-population variation can structure long-term divergence. This effect decayed after a few million years, but at least for one metric, never disappeared completely. These results are consistent with Schluter's genetic constraints model in which evolutionary trajectories on adaptive landscapes are deflected by variation within and covariation among traits.

In ectotherms, lower temperatures in high-latitude environments would theoretically reduce the annual growth rates of individuals. If slower growth and resultant smaller body size reduce fitness, individuals in higher latitudes may evolve compensatory responses. Two alternative models of such latitudinal compensation are possible: Model I: thermal reaction norms for growth rates of high-latitude individuals may be horizontally shifted to a lower range of temperatures, or Model II: reaction norms may be vertically shifted so that high-latitude individuals can grow faster across all temperatures. Model I is expected when annual growth rates in the wild are only a function of environmental temperatures, whereas Model II is expected when individuals in higher latitudes can only grow during a shorter period of a year. A variety of mixed strategies of these two models are also possible, and the magnitude of horizontal versus vertical variation in reaction norms among latitudinal populations will be indicative of the importance of “temperature” versus “seasonality” in the evolution of latitudinal compensation. However, the form of latitudinal compensation may be affected by possible genetic constraints due to the genetic architecture of reaction norms. In this study, we examine the inter- and intrapopulation variations in thermal reaction norms for growth rate of the medaka fish Oryzias latipes. Common-environment experiments revealed that average reaction norms differed primarily in elevation among latitudinal populations in a manner consistent with Model II (adaptation to “seasonality”), suggesting that natural selection in high latitudes prefers individuals that grow faster even within a shorter growing season to individuals that have longer growing seasons by growing at lower temperatures. However, intrapopulation variation in reaction norms was also vertical: some full-sibling families grew faster than others across all temperatures examined. This tendency in intrapopulation genetic variation for thermal reaction norms may have restricted the evolution of latitudinal compensation, irrespective of the underlying selection pressure.

Temperature coupling exists when changes in male signal production with temperature are paralleled by changes in female response. Such thermal effects have been observed in various ectothermic animals producing acoustic, visual, and electric signals in which the signal rate may be subject to stabilizing selection imposed by female preference. Often, coupling was considered as an adaptive function wherein male and female thermal effects coevolved under selection pressure favoring species recognition, although this assumption has not been tested definitively. We investigated thermal effects on pulse-pair rate in male song and female acceptance threshold for male song rate in an acoustic moth, Achroia grisella, in which male song rate is subject to directional selection. Male song rate and female acceptance threshold do exhibit parallel increases as temperature rises from 18°C to 30°C, but female thresholds are much lower than male song rates and the thermal effect on female response cannot augment species recognition. In further investigations using inbred lines of A. grisella we found that the male and female thermal effects are genetically correlated, and we discuss the likely sources of this covariance. We consider several explanations for the occurrence of temperature coupling in this species and suggest that it represents an emergent property arising from the neuromuscular responses to temperature that are common to several physiological systems.

Parallel speciation can occur when traits determining reproductive isolation evolve independently in different populations that experience a similar range of environments. However, a common problem in studies of parallel evolution is to distinguish this hypothesis from an alternative one in which different ecotypes arose only once in allopatry and now share a sympatric scenario with substantial gene flow between them. Here we show that the combination of a phylogenetic approach with life-history data is able to disentangle both hypotheses in the case of the intertidal marine snail Littorina saxatilis on the rocky shores of Galicia in northwestern Spain. In this system, numerous phenotypic and genetic differences have evolved between two sympatric ecotypes spanning a sharp ecological gradient, and as aside effect of the former have produced partial reproductive isolation. A mitochondrial phylogeny of these populations strongly suggests that the two sympatric ecotypes have originated independently several times. Building upon earlier work demonstrating size-based assortative mating as the main contributor to reproductive isolation among ecotypes, our analysis provides strong evidence that divergent selection across a sharp ecological gradient promoted the parallel divergence of body size and shape between two sympatric ecotypes. Thus, divergent selection occurring independently in different populations has produced the marine equivalent of host races, which may represent the first step in speciation.

Coevolutionary outcomes between interacting species are predicted to vary across landscapes, as environmental conditions, gene flow, and the strength of selection vary among populations. Using a combination of molecular, experimental, and field approaches, we describe how broad-scale patterns of environmental heterogeneity, genetic divergence, and regional adaptation have the potential to influence coevolutionary processes in the Linum marginale–Melampsora lini plant–pathogen interaction. We show that two genetically and geographically divergent pathogen lineages dominate interactions with the host across Australia, and demonstrate a hybrid origin for one of the lineages. We further demonstrate that the geographic divergence of the two lineages of M. lini in Australia is related to variation among lineages in virulence, life-history characteristics, and response to environmental conditions. When correlated with data describing regional patterns of variation in host resistance diversity and mating system these observations highlight the potential for gene flow and geographic selection mosaics to generate and maintain coevolutionary diversification in long-standing host–pathogen interactions.

The integrity of species is not fixed and may vary geographically. Here we investigate the geographic distributions and interactions of species in the Tegenaria atrica group (Araneae: Agelenidae). Detailed mapping of T. saeva and T. gigantea in England and Wales shows them to be broadly allopatric in southern England with a tightly defined, and possibly long-standing, narrow zone of parapatry in central southern England. In the north of England (Yorkshire), by contrast, the species are broadly sympatric as a result of recent range expansions. GIS techniques are used to map the species distributions and to quantify, we believe for the first time, the intimacy of interspecific interactions. The extent and nature of hybridization in these two areas is examined through regression and multivariate analyses of morphology. We show that the relative incidence of hybridization is much greater in Yorkshire than within the parapatric zone in the south. Clear patterns of asymmetric introgression are observed in both northern and southern England, with a greater impact of T. gigantea on T. saeva than vice versa. We find no sign of morphological reproductive character displacement at the zone of parapatry that might indicate reinforcement, although we cannot exclude more subtle effects, for example via cuticular pheromones. The integrity of these two species seems to be breaking down in northern England, a process that might gain momentum as the gene pools become more similar.

The geographical pattern of speciation and the relationship between floral variation and species ranges were investigated in the tribe Sinningieae (Gesneriaceae), which is found mainly in the Atlantic forests of Brazil. Geographical distribution data recorded on a grid system of 0.5 × 0.5 degree intervals and a near-complete species-level phylogenetic tree of Sinningieae inferred from a simultaneous analysis of seven DNA regions were used to address the role of geographical isolation in speciation. Geographical range overlaps between sister lineages were measured across all nodes in the phylogenetic tree and analyzed in relation to relative ages estimated from branch lengths. Although there are several cases of species sympatry in Sinningieae, patterns of sympatry between sister taxa support the predominance of allopatric speciation. The pattern of sympatry between sister taxa is consistent with range shifts following allopatric speciation, except in one clade, in which the overlapping distribution of recent sister species indicates speciation within a restricted geographical area and involving changes in pollinators and habitats. The relationship between floral divergence and regional sympatry was also examined by analyzing floral contrasts, phenological overlap, and the degree of sympatry between sister clades. Morphological contrast between flowers is not increased in sympatry and phenological divergence is more apparent between allopatric clades than between sympatric clades. Therefore, our results failed to indicate a tendency for sympatric taxa to minimize morphological and phenological overlap (geographic exclusion and/or character displacement hypotheses). Instead, they point toward adaptation in phenology to local conditions and buildup of sympatries at random with respect to flower morphology. Additional studies at a lower geographical scale are needed to identify truely coexisting species and the components of their reproductive isolation.

The mating system of flowering plant populations evolves through selection on genetically based phenotypic variation in floral traits. The physical separation of anthers and stigmas within flowers (herkogamy) is expected to be an important target of selection to limit self-fertilization. We investigated the pattern of phenotypic and genetic variation in herkogamy and its effect of self-fertilization in a broad sample of natural populations of Aquilegia canadensis, a species that is highly selfing despite strong inbreeding depression. Within natural populations, plants exhibit substantial phenotypic variation in herkogamy caused primarily by variation in pistil length rather than stamen length. Compared to other floral traits, herkogamy is much more variable and a greater proportion of variation is distributed among rather than within individuals. We tested for a genetic component of this marked phenotypic variation by growing naturally pollinated seed families from five populations in a common greenhouse environment. For three populations, we detected a significant variation in herkogamy among families, and a positive regression between parental herkogamy measured in the field and progeny herkogamy in the greenhouse, suggesting that there is often genetic variation in herkogamy within natural populations. We estimated levels of self-fertilization for groups of flowers that differed in herkogamy and show that, as expected, herkogamy was associated with reduced selfing in 13 of 19 populations. In six of these populations, we performed floral emasculations to show that this decrease in selfing is due to decreased autogamy (within-flower selfing), the mode of selfing that herkogamy should most directly influence. Taken together, these results suggest that increased herkogamy should be selected to reduce the production of low-quality selfed seed. The combination of high selfing and substantial genetic variation for herkogamy in A. canadensis is enigmatic, and reconciling this observation will require a more integrated analysis of how herkogamy influences not only self-fertilization, but also patterns of outcross pollen import and export.

A phylogenetic approach was taken to investigate the evolutionary history of seed appendages in the plant family Polygalaceae (Fabales) and determine which factors might be associated with evolution of elaiosomes through comparisons to abiotic (climate) and biotic (ant species number and abundance) timelines. Molecular datasets from three plastid regions representing 160 species were used to reconstruct a phylogenetic tree of the order Fabales, focusing on Polygalaceae. Bayesian dating methods were used to estimate the age of the appearance of ant-dispersed elaiosomes in Polygalaceae, shown by likelihood optimizations to have a single origin in the family. Topology-based tests indicated a diversification rate shift associated with appearance of caruncular elaiosomes. We show that evolution of the caruncular elaiosome type currently associated with ant dispersal occurred 54.0–50.5 million year ago. This is long after an estimated increase in ant lineages in the Late Cretaceous based on molecular studies, but broadly concomitant with increasing global temperatures culminating in the Late Paleocene–Early Eocene thermal maxima. These results suggest that although most major ant clades were present when elaiosomes appeared, the environmental significance of elaiosomes may have been an important factor in success of elaiosome-bearing lineages. Ecological abundance of ants is perhaps more important than lineage numbers in determining significance of ant dispersal. Thus, our observation that elaiosomes predate increased ecological abundance of ants inferred from amber deposits could be indicative of an initial abiotic environmental function.

Optimality models collapse the vagaries of genetics into simple trade-offs to calculate phenotypes expected to evolve by natural selection. Optimality approaches are commonly criticized for this neglect of genetic details, but resolution of this disagreement has been difficult. The importance of genetic details may be tested by experimental evolution of a trait for which an optimality model exists and in which genetic details can be studied. Here we evolved lysis time in bacteriophage T7, a virus of Escherichia coli. Lysis time is equivalent to the age of reproduction in an organism that reproduces once and then dies. Delaying lysis increases the number of offspring but slows generation time, and this trade-off renders the optimum sensitive to environmental conditions: earlier lysis is favored when bacterial hosts are dense, later lysis is favored when hosts are sparse. In experimental adaptations, T7 evolved close to the optimum in conditions favoring early lysis but not in conditions favoring late lysis. One of the late lysis adaptations exhibited no detectable phenotypic evolution despite genetic evolution; the other evolved only partly toward the expected optimum. Overall, the lysis time of the adapted phages remained closer to their starting values than predicted by the model. From the perspective of the optimality model, the experimental conditions were expected to select changes only along the postulated trade-off, but a trait outside the trade-off evolved as well. Evidence suggests that the model's failure ultimately stems from a violation of the trade-off, rather than a paucity of mutations.

Epistatic interactions between mutations are widespread. Theoretical investigations have shown that variability in epistatic effects influences fundamental evolutionary processes, yet few empirical studies have identified causes or the extent of this variation. We examined variation in epistatic effects of mutations at two host recognition sites in φX174 bacteriophage. We calculated epistatic effects from the sum of fitness effects (log scale) of two single mutants and their corresponding double mutant for five combinations of mutations in six conditions. We found that epistatic effects differed in sign, degree, and variability across conditions. The data highlight that even between single mutations at the same two sites the sign and variability of epistatic effects are affected by environment. We discuss these findings in the context of studying the role of epistasis in evolution.

This study uses the enteric bacterium Escherichia coli as an experimental system to examine evolutionary responses of bacteria to an environmental acidic-alkaline range between pH 5.3 and 7.8 (15–5000 nM [H ]). Our goal was both to test general hypotheses about adaptation to abiotic variables and to provide insights into how coliform organisms might respond to changing conditions inside and outside of hosts. Six replicate lines of E. coli evolved for 2000 generations at one of four different constant pH conditions: pH 5.3, 6.3, 7.0, or 7.8. Direct adaptation to the evolutionary environment, as well as correlated changes in other environments, was measured as a change in fitness relative to the ancestor in direct competition experiments. The pH 5.3 group had the highest fitness gains, with a highly significant increase of 20%. The pH 7.8 group had far less significant gains and much higher variance among its lines. Analysis of individual lines within these two groups revealed complex patterns of adaptation: all of the pH 5.3 lines exhibited trade-offs (reduced fitness in another environment), but only 33% of the pH 7.8 lines showed such trade-offs and one of the pH 7.8 lines demonstrated exaptation by improving fitness in the pH 5.3 environment. Although there was also prevalent exaptation in other groups to the acidic environment, there were no such cases of exaptation to alkalinity. Comparison across the entire experimental pH range revealed that the most acidic lines, the pH 5.3 group, were all specialists, in contrast to the pH 6.3 lines, which were almost all generalists. That is, although none of the pH 5.3 lines showed any correlated fitness gains, all of the pH 6.3 lines did.

Recent studies have used a variety of theoretical arguments to show that mitochondrial (mt) DNA rarely evolves as a strictly neutral marker and that selection operates on the mtDNA of many species. However, the vast majority of researchers are not convinced by these arguments because data linking mtDNA variation with phenotypic differences are limited. We investigated sequence variation in the three mtDNA and nine nuclear genes (including all isoforms) that encode the 12 subunits of cytochrome c oxidase of the electron transport chain in Drosophila. We then studied cytochrome c oxidase activity as a key aspect of mitochondrial bioenergetics and four life-history traits. In Drosophila simulans, sequence data from the three mtDNA encoded cytochrome c oxidase genes show that there are 76 synonymous and two nonsynonymous fixed differences among flies harboring siII compared with siIII mtDNA. In contrast, 13 nuclear encoded genes show no evidence of genetic subdivision associated with the mtDNA. Flies with siIII mtDNA had higher cytochrome c oxidase activity and were more starvation resistant. Flies harboring siII mtDNA had greater egg size and fecundity, and recovered faster from cold coma. These data are consistent with a causative role for mtDNA variation in these phenotypic differences, but we cannot completely rule out the involvement of nuclear genes. The results of this study have significant implications for the use of mtDNA as an assumed neutral marker and show that evolutionary shifts can involve changes in mtDNA despite the small number of genes encoded in the organelle genome.

Sexual selection based on signaling requires that signals used by females in mate choice are reliable indicators of a male's heritable total fitness. A signal and the preference for it are expected to be heritable, resulting in the maintenance of genetic covariance between these two traits. However, a recent article has proposed that signals may quickly become unreliable in the presence of both environmental variation and genotype-by-environment interaction (G × E) with crossing reaction norms, potentially compromising the mechanisms of sexual selection. Here we examine the heritability and plasticity of a male dominance advertisement in the bank vole, Clethrionomys glareolus, in stable and changing rearing environments from father to son. The bank vole is naturally exposed to considerable sources of spatial and temporal environmental variation and male reproductive success is determined by both intra- (male–male competition) and inter- (females prefer to mate with dominant males) sexual selection. Significant G × E for male dominance was found with crossing reaction norms. Plasma testosterone level (T), rather than condition, determined a male's dominance and T also showed a significant G × E. Dominance showed a considerable plasticity across environments, but was only heritable under stable conditions. We document a negative between-environments correlation of male dominance, suggesting that when the environment changes between father and son, the dominance signal is unreliable to females and sexual selection may be compromised. We discuss how G × E and environmental variation interacting with other mechanisms may preserve the reliability of signals and thus the mechanism of sexual selection itself.

A weighted-average model, which reliably estimates endocranial volume from three external measurements of the neurocranium of extant taxa in the mammalian order Carnivora, was tested for its applicability to fossil taxa by comparing model-estimated endocranial volumes to known endocast volumes. The model accurately reproduces endocast volumes for a wide array of fossil taxa across the crown radiation of the Carnivora, three stem carnivoramorphan taxa, and Pleistocene fossils of two extant species. Applying this model to fossil taxa without known endocast volumes expanded the sample of fossil taxa with estimated brain volumes in the carnivoran suborder Caniformia from 11 to 60 taxa. This then allowed a comprehensive assessment of the evolution of relative brain size across this clade. An allometry of brain volume to body mass was calculated on phylogenetically independent contrasts for the set of extant taxa, and from this, log-transformed encephalization quotients (logEQs) were calculated for all taxa, extant, and fossil. A series of Mann–Whitney tests demonstrated that the distributions of logEQs for taxa early in caniform evolutionary history possessed significantly lower median logEQs than extant taxa. Median logEQ showed a pronounced shift around the Miocene–Pliocene transition. Support tests, based on likelihood ratios, demonstrated that the variances of these distributions also were significantly lower than among modern taxa, but logEQ variance increased gradually through the history of the clade, not abruptly. Reconstructions of ancestral logEQs using weighted squared-change parsimony demonstrate that increased encephalization is observed across all major caniform clades (with the possible exception of skunks) and that these increases were achieved in parallel, although an “ancestor-descendant differencing” method could not rule out drift as a hypothesis. Peculiarities in the estimated logEQs for the extinct caniform family Amphicyonidae were also investigated; these unusual patterns are likely due to a unique allometry in scaling brain to body size in this single clade.

Although the evolutionary importance of natural hybridization has been debated for decades, it has become increasingly clear that hybridization plays a fundamental role in the evolution of many plant and animal taxa, sometimes resulting in the formation of entirely new species. Although some hybrid species retain the base chromosome number of their parents, others combine the full chromosomal complements of their progenitors. Hybrid speciation can thus produce two fundamentally different types of evolutionary lineages, yet relatively little is known about the factors influencing ploidy level in hybrid neospecies. We estimated genetic divergence between species pairs that have given rise to homoploid and polyploid hybrid species and found that divergence is significantly greater for the parents of polyploids, even after controlling for potentially confounding factors. Our data thus provide the first direct evidence in support of the notion that the extent of genomic divergence between hybridizing species influences the likelihood of diploid versus polyploid hybrid speciation.

Despite important advances in the last few years, the evolution of reproductive isolation (RI) remains an unresolved and critical gap in our understanding of speciation processes. In this study, we investigated the evolution of RI among species of the parasitic fungal species complex Microbotryum violaceum, which is responsible for anther smut disease of the Caryophyllaceae. We found no evidence for significant positive assortative mating by M. violaceum even over substantial degrees of genetic divergence, suggesting a lack of prezygotic isolation. In contrast, postzygotic isolation increased with the genetic distance between mating partners when measured as hyphal growth. Total RI, measured as the ability of the pathogen to infect and produce a diploid progeny in the host plant, was significantly and positively correlated with genetic distance, remaining below complete isolation for most of the species pairs. The results of this study, the first one on the time course of speciation in a fungus, are therefore consistent with previous works showing that RI generally evolves gradually with genetic distance, and thus presumably with time. Interestingly, prezygotic RI due to gamete recognition did not increase with genetic distance, in contrast to the pattern found in plants and animals.