A biogeographic genetic approach for testing the role of reinforcement: the case of Drosophila pseudoobscura and D. persimilis

Abstract.— The role of reinforcement in speciation can be explained by two distinct models. In model I, two diverged populations hybridize and produce fertile hybrids that successfully backcross (hybridization with gene flow). In model II, two populations hybridize but succeeding backcrosses are unproductive (hybridization without gene flow). Using Drosophila persimilis and D. pseudoobscura , we have tested model I by comparing the extent of heterospecific introgression in sympatric versus allopatric populations. We show that certain expectations of this particular model of reinforcement, which is based on hybridization and gene flow between divergent populations after secondary contact, are not realized in these two species. The evidence consists of the similarity of genetic distances as well as proportions of unique/rare alleles between sympatric and allopatric heterospecific populations and a negative correlation between genetic distance and geographical distance between heterospecific populations, which suggests ecological differentiation. This approach in quantifying differential gene flow has important consequences to studies that compare sympatric and allopatric isolation using genetic distance. Following model I, one would expect a pattern of higher prezygotic isolation in sympatric species compared to allopatric species of the same genetic distance simply as a result of an underestimation of genetic distance due to introgression between sympatric populations. We suggest more parsimonious explanations such as reinforcement without genetic exchange (model II) and ecological differentiation, which require high levels of preexisting reproductive isolation between populations.     

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

ABSTRACT: BACKGROUND: Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated ('allopatric') populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations. RESULTS: As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation. CONCLUSIONS: Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of 'genomic islands of speciation'.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Influence of voltine ecotype and geographic distance on genetic and haplotype variation in the Asian corn borer

Diapause is an adaptive dormancy strategy by which arthropods endure extended periods of adverse climatic conditions. Seasonal variation in larval diapause initiation and duration in Ostrinia furnacalis may influence adult mating generation number (voltinism) across different local environments. The degree to which voltine ecotype, geographic distance, or other ecological factors influence O. furnacalis population genetic structure remains uncertain. Genetic differentiation was estimated between voltine ecotypes collected from 8 locations. Mitochondrial haplotypes were significantly different between historically allopatric univoltine and bivoltine locations, but confounded by a strong correlation with geographic distance. In contrast, single nucleotide polymorphism (SNP) genotypes show low but significant levels of variation and a lack of influence of geographic distance between allopatric voltine locations. Regardless, 11 of 257 SNP loci were predicted to be under selection, suggesting population genetic homogenization except at loci proximal to factors putatively under selection. These findings provide evidence of haplotype divergent voltine ecotypes that may be maintained in allopatric and sympatric areas despite relatively high rates of nuclear gene flow, yet influence of voltinism on maintenance of observed haplotype divergence remains unresolved.     

Sympatric,parapatric or allopatric: the most important way to classify speciation?

The most common classification of modes of speciation begins with the spatial context in which divergence occurs: sympatric, parapatric or allopatric. This classification is unsatisfactory because it divides a continuum into discrete categories, concentrating attention on the extremes, and it subordinates other dimensions on which speciation processes vary, such as the forces driving differentiation and the genetic basis of reproductive isolation. It also ignores the fact that speciation is a prolonged process that commonly has phases in different spatial contexts. We use the example of local adaptation and partial reproductive isolation in the intertidal gastropod Littorina saxatilis to illustrate the inadequacy of the spatial classification of speciation modes. Parallel divergence in shell form in response to similar environmental gradients in England, Spain and Sweden makes this an excellent model system. However, attempts to demonstrate 'incipient' and 'sympatric' speciation involve speculation about the future and the past. We suggest that it is more productive to study the current balance between local adaptation and gene flow, the interaction between components of reproductive isolation and the genetic basis of differentiation.     

Character displacement in sailfin mollies, <Emphasis Type="Italic">Poecilia latipinna</Emphasis>: allozymes and behavior

Synopsis We analyzed variation in allozymes and mating preferences in 12 populations across much of the range of the sailfin molly, Poecilia latipinna. Sailfin mollies can be sympatric with its sexual parasite Amazon mollies, P. formosa. Amazon mollies must co-exist and mate with bisexual males of closely related species (including sailfin mollies) to induce embryogenesis but inheritance is strictly maternal. Where sailfin and Amazon mollies are sympatric there is evidence of reproductive character displacement as males show a significantly stronger mating preference for sailfin molly females over Amazon mollies compared to preferences of males from allopatric populations. From the allozyme data we found a moderate amount of genetic variation across all populations but this variation did not reveal significant partitioning between sympatric and allopatric populations. Additionally, we found no evidence for isolation by distance as genetic distance was not significantly correlated with geographic distance. While allozyme variation also did not significantly correlate with male mating preferences, there was a significant correlation between male mating preferences and geographic distance. This correlation between mating preferences and geographic distance may have arisen from coevolution with Amazon mollies resulting in reproductive character displacement. Taken together, the distribution of genetic and behavioral variation among sympatric and allopatric populations suggests that behavioral evolution has outpaced evolution at the allozyme loci we examined in P. latipinna.     

MULTIPLE MODES OF SPECIATION INVOLVED IN THE PARALLEL EVOLUTION OF SYMPATRIC MORPHOTYPES OF LAKE WHITEFISH (COREGONUS CLUPEAFORMIS,SALMONIDAE)

We performed a phylogenetic analysis of mtDNA variation among seven sympatric pairs of dwarf and normal morphotypes of whitefish from northern Québec and the St. John River drainage to address three questions relevant to understanding their radiation. Are all sympatric pairs reproductively isolated? Do phylogenetic analyses confirm that sympatric whitefish morphotypes found in eastern North America represent the outcome of polyphyletic evolutionary events? If so, did all sympatric pairs from the St. John River drainage originate from the same scenario of allopatric divergence and secondary contact? The hypothesis of genetic differentiation was supported for all sympatric pairs from the St. John River drainage, whereas lack of mtDNA diversity precluded any test of reproductive isolation for northern Québec populations. Patterns of mtDNA variation confirmed that dwarf and normal morphotypes evolved in parallel among independent, yet closely related, lineages, thus providing indirect evidence for the role of natural selection in promoting phenotypic radiation in whitefish. Patterns of mtDNA diversity among sympatric pairs of the St. John River indicated a complex picture of whitefish evolution that implied sympatric divergence and multiple allopatric divergence/secondary contact events on a small geographic scale. These results suggests that ecological opportunities, namely trophic niche availability, may promote population divergence in whitefish.     

REPRODUCTIVE CHARACTER DISPLACEMENT OF MALE STICKLEBACK MATE PREFERENCE: REINFORCEMENT OR DIRECT SELECTION?

Reproductive isolation can evolve between species as a byproduct of adaptation to different niches, through reinforcement, and by direct selection on mating preferences. We investigated the potential role of direct selection in the reproductive isolation between sympatric species of threespine sticklebacks. Each sympatric pair consists of a small-bodied limnetic species and large-bodied benthic species. We compared the mate preferences and courtship behavior of males from one sympatric limnetic population and two allopatric populations. We used limnetic-like allopatric populations to control for the effects of ecological character displacement and adaptation to different niches on mate preferences. The sympatric limnetic males preferred the small limnetic females, whereas the allopatric limnetic-like males preferred the large benthic females, suggesting that adaptation to the limnetic niche does not automatically confer a preference for small limnetic females. This reproductive character displacement of male preference is consistent with the predictions of both reinforcement and direct selection on mate preferences. To test for direct selection, we assessed a prediction of one proposed mechanism: predation by benthic females on eggs guarded by limnetic males. The allopatric males come from populations in which there is no egg predation. Sympatric limnetic males were more aggressive toward benthic females than toward limnetic females, whereas the allopatric limnetic-like males did not treat the two types of females differently. The contrast in male behavior suggests that egg predation has shaped male preferences. Direct selection is potentially more effective than indirect selection via reinforcement, and it is likely that it has been important in building up reproductive isolation between limnetic and benthic sticklebacks.     

Evolutionary insights into Rhinolophus episcopus (Chiroptera,Rhinolophidae) in China: Isolation by distance,environment, or sensory system?

Evolutionary processes can be influenced by several factors, such as geographic isolation, environmental selection, and sensory variation. For most nocturnal bats, echolocation is the primary sensory system used to prey and communicate, and plays important roles in chiropteran diversification and evolution. Understanding the relative contribution of geography, the environment, and this sensory system to population genetic divergence can elucidate the processes involved in bat incipient speciation and evolution. In this study, we collected spatial and environmental information, echolocation calls, as well as the previously published genetic data (six microsatellite loci and the mitochondrial cytochrome b gene) of widely distributed Rhinolophus episcopus populations to test three hypotheses for nuclear and mitochondrial divergence (isolation by distance, isolation by environment, and isolation by sensory variation) and unveil the factors that drive intraspecific genetic differentiation. The moderate level of nuclear differentiation was correlated with geographic/spatial distance and acoustic variation, whereas the relatively high level of mitochondrial differentiation was mainly associated with acoustic divergence. No significant correlation was observed between genetic divergence and environmental variables. Among the three factors, acoustic divergence explained the highest percentage of both nuclear and mitochondrial divergence. Thus, our results indicate that sensory variation may have played important roles in driving population isolation early in bat speciation, which is consistent with the hypothesis of isolation by sensory variation. Our study emphasizes the need to consider more factors, especially sensory traits, and combine multiple statistical methods in landscape genetic studies to test their potential contributions to driving population divergence.     

A test of the chromosomal rearrangement model of speciation in Drosophila pseudoobscura

Recent studies suggest that chromosomal rearrangements play a significant role in speciation by preventing recombination and maintaining species persistence despite interspecies gene flow. Factors conferring adaptation or reproductive isolation are maintained in rearranged regions in the face of hybridization, while such factors are eliminated from collinear regions. As a direct test of this rearrangement model, we evaluated the genetic basis of hybrid male sterility in a sympatric species pair, Drosophila pseudoobscura pseudoobscura and D. persimilis, and an allopatric species pair, D. pseudoobscura bogotana and D. persimilis. Our results are consistent with the proposed model: virtually all of the sterility factors in the former pair are associated with three inverted regions, whereas sterility factors are present in the collinear regions in the latter pair. These findings indicate recombination and selection may have eliminated sterility factors outside the inverted regions between D. p. pseudoobscura and D. persimilis, suggesting chromosomal rearrangements may facilitate species persistence despite hybridization.     

DOES SELECTION ON FLORAL ODOR PROMOTE DIFFERENTIATION AMONG POPULATIONS AND SPECIES OF THE SEXUALLY DECEPTIVE ORCHID GENUS OPHRYS?

Sexually deceptive orchids from the genus Ophrys attract their pollinators primarily through the chemical mimicry of female hymenopteran sex pheromones, thereby deceiving males into attempted matings with the orchid labellum. Floral odor traits are crucial for the reproductive success of these pollinator-limited orchids, as well as for maintaining reproductive isolation through the attraction of specific pollinators. We tested for the signature of pollinator-mediated selection on floral odor by comparing intra and interspecific differentiation in odor compounds with that found at microsatellite markers among natural populations. Three regions from southern Italy were sampled. We found strong floral odor differentiation among allopatric populations within species, among allopatric species and among sympatric species. Population differences in odor were also reflected in significant variation in the attractivity of floral extracts to the pollinator, Colletes cunicularius. Odor compounds that are electrophysiologically active in C. cunicularius males, especially alkenes, were more strongly differentiated among conspecific populations than nonactive compounds in the floral odor. In marked contrast to these odor patterns, there was limited population or species level differentiation in microsatellites (FST range 0.005 to 0.127, mean FST 0.075). We propose that the strong odor differentiation and lack of genetic differentiation among sympatric taxa indicates selection imposed by the distinct odor preferences of different pollinating species. Within species, low FST values are suggestive of large effective population sizes and indicate that divergent selection rather than genetic drift accounts for the strong population differentiation in odor. The higher differentiation in active versus non-active odor compounds suggests that divergent selection among orchid populations may be driven by local pollinator preferences for those particular compounds critical for pollinator attraction.     

Genome‐wide differentiation in closely related populations: the roles of selection and geographic isolation

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as ‘isolation by adaptation’ (IBA). IBA can be complementary to the well‐known pattern of ‘isolation by distance’ (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome‐wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole‐genome sequence for H. rustica, to which we aligned a set of 9493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behaviour and sexual signalling, as well as geographic distance, together explain over 70% of genome‐wide divergence among populations. Controlling for IBD, we find 42% of genomewide divergence is attributable to IBA through pairwise differences in traits related to migratory behaviour and sexual signalling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population‐level differentiation in this group of closely related populations.     

Genome‐wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.     

Sexual conflict does not drive reproductive isolation in experimental populations of Drosophila pseudoobscura

Sexual conflict has been predicted to drive reproductive isolation by generating arbitrary but rapid coevolutionary changes in reproductive traits among allopatric populations. A testable prediction of this proposal is that allopatric populations experiencing different levels of sexual conflict should exhibit different levels of reproductive isolation. We tested this prediction using experimentally evolved populations of the promiscuous Drosophila pseudoobscura. We manipulated sexual conflict by enforcing either monogamy, maintaining natural levels of promiscuity, or elevating promiscuity. Within each treatment, we carried out sympatric and allopatric crosses using replicated populations and examined pre-zygotic (number of mating pairs, mating speed and copulation duration) and post-zygotic (hybrid inviability and sterility) indicators of reproductive isolation. After 50 generations of selection, none of the measures conformed to predictions of sexual conflict driving reproductive isolation. Our results cannot be explained by lack of genetic variation or weak selection and suggest that sexual conflict may not be a widespread engine of speciation.     

Reinforcement shapes clines in female mate discrimination in Drosophila subquinaria

Reinforcement of species boundaries may alter mate recognition in a way that also affects patterns of mate preference among conspecific populations. In the fly Drosophila subquinaria, females sympatric with the closely related species D. recens reject mating with heterospecific males as well as with conspecific males from allopatric populations. Here, we assess geographic variation in behavioral isolation within and among populations of D. subquinaria and use cline theory to understand patterns of selection on reinforced discrimination and its consequences for sexual isolation within species. We find that selection has fixed rejection of D. recens males in sympatry, while significant genetic variation in this behavior occurs within allopatric populations. In conspecific matings sexual isolation is also asymmetric and stronger in populations that are sympatric with D. recens. The clines in behavioral discrimination within and between species are similar in shape and are maintained by strong selection in the face of gene flow, and we show that some of their genetic basis may be either shared or linked. Thus, while reinforcement can drive extremely strong phenotypic divergence, the long‐term consequences for incipient speciation depend on gene flow, genetic linkage of discrimination traits, and the cost of these behaviors in allopatry.     

Patterns of genetic,phenotypic, and acoustic variation across a chiffchaff (Phylloscopus collybita abietinus/tristis) hybrid zone

Characterizing patterns of evolution of genetic and phenotypic divergence between incipient species is essential to understand how evolution of reproductive isolation proceeds. Hybrid zones are excellent for studying such processes, as they provide opportunities to assess trait variation in individuals with mixed genetic background and to quantify gene flow across different genomic regions. Here, we combine plumage, song, mtDNA and whole‐genome sequence data and analyze variation across a sympatric zone between the European and the Siberian chiffchaff (Phylloscopus collybita abietinus/tristis) to study how gene exchange between the lineages affects trait variation. Our results show that chiffchaff within the sympatric region show more extensive trait variation than allopatric birds, with a large proportion of individuals exhibiting intermediate phenotypic characters. The genomic differentiation between the subspecies is lower in sympatry than in allopatry and sympatric birds have a mix of genetic ancestry indicating extensive ongoing and past gene flow. Patterns of phenotypic and genetic variation also vary between regions within the hybrid zone, potentially reflecting differences in population densities, age of secondary contact, or differences in mate recognition or mate preference. The genomic data support the presence of two distinct genetic clades corresponding to allopatric abietinus and tristis and that genetic admixture is the force underlying trait variation in the sympatric region—the previously described subspecies (“fulvescens”) from the region is therefore not likely a distinct taxon. In addition, we conclude that subspecies identification based on appearance is uncertain as an individual with an apparently distinct phenotype can have a considerable proportion of the genome composed of mixed alleles, or even a major part of the genome introgressed from the other subspecies. Our results provide insights into the dynamics of admixture across subspecies boundaries and have implications for understanding speciation processes and for the identification of specific chiffchaff individuals based on phenotypic characters.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.     

Location-specific sympatric morphological divergence as a possible response to species interactions in West Virginia Plethodon salamander communities

1. The competitive interactions of closely related species have long been considered important determinants of community composition and a major cause of phenotypic diversification. However, while patterns such as character displacement are well documented, less is known about how local adaptation influences diversifying selection from interspecific competition. 2. We examined body size and head shape variation among allopatric and sympatric populations of two salamander species, the widespread Plethodon cinereus and the geographically restricted P. nettingi. We quantified morphology from 724 individuals from 20 geographical localities throughout the range of P. nettingi. 3. Plethodon nettingi was more robust in cranial morphology relative to P. cinereus, and sympatric localities were more robust relative to allopatric localities. Additionally, there was significantly greater sympatric head shape divergence between species relative to allopatric communities, and sympatric localities of P. cinereus exhibited greater morphological variation than sympatric P. nettingi. 4. The sympatric morphological divergence and increase in cranial robustness of one species (P. nettingi) were similar to observations in other Plethodon communities, and were consistent with the hypothesis of interspecific competition. These findings suggest that interspecific competition in Plethodon may play an important role in phenotypic diversification in this group. 5. The increase in among-population variance in sympatric P. cinereus suggests a species-specific response to divergent natural selection that is influenced in part by other factors. We hypothesize that enhanced morphological flexibility and ecological tolerance allow P. cinereus to more rapidly adapt to local environmental conditions, and initial differences among populations have allowed the evolutionary response of P. cinereus to vary across replicate sympatric locations, resulting in distinct evolutionary trajectories of morphological change.     

Host-related genetic differentiation in the anther smut fungus Microbotryum violaceum in sympatric, parapatric and allopatric populations of two host species Silene latifolia and S. dioica

We investigated genetic diversity in West European populations of the fungal pathogen Microbotryum violaceum in sympatric, parapatric and allopatric populations of the host species Silene latifolia and S. dioica, using four polymorphic microsatellite loci. In allopatric host populations, the fungus was highly differentiated by host species, exhibiting high values of F(ST) and R(ST), and revealed clear and distinct host races. In sympatric and parapatric populations we found significant population differentiation as well, except for one sympatric population in which the two host species grew truly intermingled. The mean number of alleles per locus for isolates from each of the host species was significantly higher in sympatric/parapatric than in allopatric populations. This suggests that either gene flow between host races in sympatry, or in case of less neutral loci, selection in a more heterogeneous host environment can increase the level of genetic variation in each of the demes. The observed pattern of host-related genetic differentiation among these geographically spread populations suggest a long-term divergence between these host races. In sympatric host populations, both host races presumably come in secondary contact, and host-specific alleles are exchanged depending on the amount of fungal gene flow.