Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

The evolution of selfing from outcrossing ancestors in Brassicaceae: what have we learned from variation at the S-locus?

Evolutionary transitions between mating systems have occurred repetitively and independently in flowering plants. One of the most spectacular advances of the recent empirical literature in the field was the discovery of the underlying genetic machinery, which provides the opportunity to retrospectively document the scenario of the outcrossing to selfing transitions in a phylogenetic perspective. In this review, we explore the literature describing patterns of polymorphism and molecular evolution of the locus controlling self-incompatibility (S-locus) in selfing species of the Brassicaceae family in order to document the transition from outcrossing to selfing, a retrospective approach that we describe as the 'mating system genes approach'. The data point to strikingly contrasted scenarios of transition from outcrossing to selfing. We also perform original analyses of the fully sequenced genomes of four species showing self-compatibility, to compare the orthologous S-locus region with that of functional S-locus haplotypes. Phylogenetic analyses suggest that all species we investigated evolved independently towards loss of self-incompatibility, and in most cases almost intact sequences of either of the two S-locus genes suggest that these transitions occurred relatively recently. The S-locus region in Aethionema arabicum, representing the most basal lineage of Brassicaceae, showed unusual patterns so that our analysis could not determine whether self-incompatibility was lost secondarily, or evolved in the core Brassicaceae after the split with this basal lineage. Although the approach we detail can only be used when mating system genes have been identified in a clade, we suggest that its integration with phylogenetic and population genetic approaches should help determine the main routes of this predominant mating system shift in plants.

Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri.

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.

Molecular evolution within and between self-incompatibility specificities.

Genes under multiallelic balancing selection have sharply contrasted evolutionary dynamics across timescales, with much longer coalescence time among functionally distinct allelic lines but much shorter coalescence time among gene copies within allelic lines as compared with the genomic background. In this paper, we combine theoretical and empirical approaches to investigate patterns of molecular evolution within and between self-incompatibility (SI) specificities. We first use numerical simulations to investigate coalescence times within allelic lines in a subdivided population for a sporophytic SI system. We then report on a comprehensive analysis of nucleotide polymorphism among gene copies within five distinct allelic lines in the closely related Arabidopsis halleri and Arabidopsis lyrata. In line with our model predictions, we find that the observed level of polymorphism among gene copies was generally low but differed among allelic lines. The data provide compelling direct evidence for recombination and/or gene conversion not only within the two most recessive allelic lines but also between two closely related but distinct allelic lines, suggesting that recombination at the Arabidopsis SI locus is possible in the absence of large sequence divergence among haplotypes. We observed shared polymorphic sites between the two species in one allelic line and strikingly similar haplotypes in another allelic line. We discuss whether convergent evolution may have led to this pattern and suggest that these observations are consistent with ongoing or very recent introgression, as previously documented.

High paternal diversity in the self-incompatible herb Arabidopsis halleri despite clonal reproduction and spatially restricted pollen dispersal.

The number of sires fertilizing a given dam is a key parameter of the mating system in species with spatially restricted offspring dispersal, since genetic relatedness among maternal sibs determines the intensity of sib competition. In flowering plants, the extent of multiple paternity is determined by factors such as floral biology, properties of the pollen vector, selfing rate, spatial organization of the population, and genetic compatibility between neighbours. To assess the extent of multiple paternity and identify ecological factors involved, we performed a detailed study of mating patterns in a small population of a self-incompatible clonal herb, Arabidopsis halleri. We mapped and genotyped 364 individuals and 256 of their offspring at 12 microsatellite loci and jointly analysed the level of multiple paternity, pollen and seed dispersal, and spatial genetic structure. We found very low levels of correlated paternity among sibs (P(full-sib) = 3.8%) indicating high multiple paternity. Our estimate of the outcrossing rate was 98.7%, suggesting functional self-incompatibility. The pollen dispersal distribution was significantly restricted (mean effective pollen dispersal distance: 4.42 m) but long-distance successful pollination occurred and immigrating pollen was at most 10% of all pollination events. Patterns of genetic structure indicated little extent of clonal reproduction, and a low but significant spatial genetic structure typical for a self-incompatible species. Overall, in spite of restricted pollen dispersal, the multiple paternity in this self-incompatible species was very high, a result that we interpret as a consequence of high plant density and high pollinator service in this population.

Unequal allelic frequencies at the self-incompatibility locus within local populations of Prunus avium L.: an effect of population structure?

In this paper, we investigated the genetic structure and distribution of allelic frequencies at the gametophytic self-incompatibility locus in three populations of Prunus avium L. In line with theoretical predictions under balancing selection, genetic structure at the self-incompatibility locus was almost three times lower than at seven unlinked microsatellites. Furthermore, we found that S-allele frequencies in wild cherry populations departed significantly from the expected isoplethic distribution towards which balancing selection is expected to drive allelic frequencies (i.e. identical frequency equal to the inverse of the number of alleles in the population). To assess whether this departure could be caused either by drift alone or by population structure, we used numerical simulations to compare our observations with allelic frequency distributions expected : (1) within a single deme from a subdivided population with various levels of differentiation; and (2) within a finite panmictic population with identical allelic diversity. We also investigated the effects of sample size and degree of population structure on tests of departure from isoplethic equilibrium. Overall, our results showed that the observed allele frequency distributions were consistent with a model of subdivided population with demes linked by moderate migration rate.

Chloroplast DNA variation and postglacial recolonization of common ash (Fraxinus excelsior L.) in Europe.

We used chloroplast polymerase chain reaction-restriction-fragment length polymorphism (PCR-RFLP) and chloroplast microsatellites to assess the structure of genetic variation and postglacial history across the entire natural range of the common ash (Fraxinus excelsior L.), a broad-leaved wind-pollinated and wind-dispersed European forest tree. A low level of polymorphism was observed, with only 12 haplotypes at four polymorphic microsatellites in 201 populations, and two PCR-RFLP haplotypes in a subset of 62 populations. The clear geographical pattern displayed by the five most common haplotypes was in agreement with glacial refugia for ash being located in Iberia, Italy, the eastern Alps and the Balkan Peninsula, as had been suggested from fossil pollen data. A low chloroplast DNA mutation rate, a low effective population size in glacial refugia related to ash's life history traits, as well as features of postglacial expansion were put forward to explain the low level of polymorphism. Differentiation among populations was high (GST= 0.89), reflecting poor mixing among recolonizing lineages. Therefore, the responsible factor for the highly homogeneous genetic pattern previously identified at nuclear microsatellites throughout western and central Europe (Heuertz et al. 2004) must have been efficient postglacial pollen flow. Further comparison of variation patterns at both marker systems revealed that nuclear microsatellites identified complex differentiation patterns in south-eastern Europe which remained undetected with chloroplast microsatellites. The results suggest that data from different markers should be combined in order to capture the most important genetic patterns in a species.

New insights from fine-scale spatial genetic structure analyses in plant populations.

Many empirical studies have assessed fine-scale spatial genetic structure (SGS), i.e. the nonrandom spatial distribution of genotypes, within plant populations using genetic markers and spatial autocorrelation techniques. These studies mostly provided qualitative descriptions of SGS, rendering quantitative comparisons among studies difficult. The theory of isolation by distance can predict the pattern of SGS under limited gene dispersal, suggesting new approaches, based on the relationship between pairwise relatedness coefficients and the spatial distance between individuals, to quantify SGS and infer gene dispersal parameters. Here we review the theory underlying such methods and discuss issues about their application to plant populations, such as the choice of the relatedness statistics, the sampling scheme to adopt, the procedure to test SGS, and the interpretation of spatial autocorrelograms. We propose to quantify SGS by an 'Sp' statistic primarily dependent upon the rate of decrease of pairwise kinship coefficients between individuals with the logarithm of the distance in two dimensions. Under certain conditions, this statistic estimates the reciprocal of the neighbourhood size. Reanalysing data from, mostly, published studies, the Sp statistic was assessed for 47 plant species. It was found to be significantly related to the mating system (higher in selfing species) and to the life form (higher in herbs than trees), as well as to the population density (higher under low density). We discuss the necessity for comparing SGS with direct estimates of gene dispersal distances, and show how the approach presented can be extended to assess (i) the level of biparental inbreeding, and (ii) the kurtosis of the gene dispersal distribution.

Estimating seed vs. pollen dispersal from spatial genetic structure in the common ash.

Spatial genetic structure was analysed with five highly polymorphic microsatellite loci in a Romanian population of common ash (Fraxinus excelsior L.), a wind-pollinated and wind-dispersed tree species occurring in mixed deciduous forests over almost all of Europe. Contributions of seed and pollen dispersal to total gene flow were investigated by analysing the pattern of decrease in kinship coefficients among pairs of individuals with geographical distance and comparing it with simulation results. Plots of kinship against the logarithm of distance were decomposed into a slope and a shape component. Simulations showed that the slope is informative about the global level of gene flow, in agreement with theoretical expectations, whereas the shape component was correlated with the relative importance of seed vs. pollen dispersal. Hence, our results indicate that insights into the relative contributions of seed and pollen dispersal to overall gene flow can be gained from details of the pattern of spatial genetic structure at biparentally inherited loci. In common ash, the slope provided an estimate of total gene dispersal in terms of Wright's neighbourhood size of Nb = 519 individuals. No precise estimate of seed vs. pollen flow could be obtained from the shape because of the stochasticity inherent to the data, but the parameter combinations that best fitted the data indicated restricted seed flow, sigmas pound 14 m, and moderate pollen flow, 70 m pound sigmap pound 140 m.

Distinction between cultivated and wild chicory gene pools using AFLP markers.

The cultivation area of industrial chicory, Cichorium intybus L. cv Sativum, coincides with the natural distribution area of its wild relative, C. intybus L., which could lead to gene flow between wild and cultivated types. The genetic diversity within and between the two types has therefore been studied using AFLP genotyping of samples from 12 wild populations collected in Belgium and ten commercial varieties. The genotyping of 233 individuals allowed the identification of 254 AFLP markers. Similar levels of genetic diversity were observed within wild populations and cultivated varieties, suggesting the absence of any strong bottleneck in the history of the cultivated types. The phylogenetic analysis pointed to a monophyletic origin of cultivated varieties as compared to the local wild populations studied, hence the two types of chicory form two separate gene pools. The genotyping of some individuals sampled in ruderal sites clearly showed that they belong to the cultivated gene pool, which suggests the existence of feral or weedy types. The low differentiation observed among wild populations indicates that gene flow might be important in this species.

Genetic structure and mating systems of metallicolous and nonmetallicolous populations of Thlaspi caerulescens.

• Genetic variation structure and breeding system were investigated in metallicolous (MET) and nonmetallicolous (NONMET) populations of the heavy-metal hyperaccumulator Thlaspi caerulescens from Belgium, Luxembourg and the French Mediterranean region. • Allozyme variation showed a clear differentiation between the two ecotypes in Belgium and Luxembourg but not in southern France, in line with the lower degree of geographical isolation between the two ecotypes in this region. • In both regions inbreeding coefficient (FIS ) of NONMET populations was significantly higher compared to MET populations. Pollen/Ovule ratios were significantly higher in MET compared with NONMET populations. • These results suggest that NONMET populations of T. caerulescens are more self-fertile than their MET counterparts. This contrasts with earlier studies on other metal-tolerant species in which selfing rates were higher in MET populations. This pattern may be explained by the fact that both ecotypes are not in sympatry in the populations studied, and therefore reproductive isolation has not been selected to maintain the adaptations to metal-contaminated soils. In addition, higher selfing rates in NONMET populations may have evolved as a mechanism of reproductive assurance, because these populations generally are at low densities.

Data from amplified fragment length polymorphism (AFLP) markers show indication of size homoplasy and of a relationship between degree of homoplasy and fragment size.

We investigate the distribution of sizes of fragments obtained from the amplified fragment length polymorphism (AFLP) marker technique. We find that empirical distributions obtained in two plant species, Phaseolus lunatus and Lolium perenne, are consistent with the expected distributions obtained from analytical theory and from numerical simulations. Our results indicate that the size distribution is strongly asymmetrical, with a much higher proportion of small than large fragments, that it is not influenced by the number of selective nucleotides nor by genome size but that it may vary with genome-wide GC-content, with a higher proportion of small fragments in cases of lower GC-content when considering the standard AFLP protocol with the enzyme MseI. Results from population samples of the two plant species show that there is a negative relationship between AFLP fragment size and fragment population frequency. Monte Carlo simulations reveal that size homoplasy, arising from pulling together nonhomologous fragments of the same size, generates patterns similar to those observed in P. lunatus and L. perenne because of the asymmetry of the size distribution. We discuss the implications of these results in the context of estimating genetic diversity with AFLP markers.

Allozyme segregation and inter-cytotype reproductive barriers in the polyploid complex Centaurea jacea.

In eastern Belgium, diploid and tetraploid knapweeds (Centaurea jacea L. s. l.) show a parapatric distribution with a contact zone. To compare the success of inter- and intra-cytotype crosses, and to investigate the pattern of allozyme segregation, controlled crosses involving plants from this contact zone were performed. Polysomic inheritance at two loci was observed in tetraploids, suggesting an autopolyploid origin. Two crosses allowed the detection of double reduction events in tetraploids at one locus, but no such event was detected among the 217 progenies. Null alleles were detected in tetraploids at two loci. Both cytotypes were highly self-incompatible. Inter-cytotype crosses were much less successful than intra-cytotype crosses in terms of seed set (five-fold reduction) and germination rates (three-fold reduction), suggesting a 'triploid block' effect. Progenies from inter-cytotype crosses most often had the maternal cytotype and resulted from selfing or pollen contamination. However, six triploids were detected. The significance of these results for the understanding of the geographical distribution of the cytotypes and inter-cytotype gene flow is discussed.

Assessment of genetic structure within and among Bulgarian populations of the common ash (Fraxinus excelsior L.).

We analysed genetic variation within and between populations of the common ash from Bulgaria in order to extract biological information useful in the context of conservation management of eastern European genetic resources of noble hardwood species. A total of 321 trees from three regions of Bulgaria were typed at six highly polymorphic microsatellite loci. Analysis of within-population inbreeding suggests an upper boundary value of 2.7% for the selfing rate. Significant spatial genetic structure consistent with models of isolation by distance was detected within four out of 10 populations as well as among populations. Estimates of neighbourhood size in the range 38-126 individual trees were obtained based on spatial genetic structure analyses at either the intrapopulation or interpopulation level. Differentiation among populations explained only about 8.7% of total genetic diversity. These results are discussed in comparison with data from social broad-leaved species such as oak and beech.

Patterns of allozyme variation in diploid and tetraploid Centaurea jacea at different spatial scales.

The extent and spatial patterns of genetic variation at allozyme markers were investigated within and between diploid and autotetraploid knapweeds (Centaurea jacea L. sensu lato, Asteraceae) at contrasted geographic scales: (1) among populations sampled from a diploid-tetraploid contact zone in the northeastern part of the Belgian Ardennes, and (2) within mixed populations from that zone where diploids and tetraploids coexist. Our data were also compared with a published dataset by Sommer (1990) describing allozyme variation in separate diploid and tetraploid knapweeds populations collected throughout Europe. Genetic diversity was higher in tetraploids. In the Belgian Ardennes and within the mixed populations, both cytotypes had similar levels of spatial genetic structure, they were genetically differentiated, and their distributions of allele frequencies were not spatially correlated. In contrast, at the European scale, diploids and tetraploids did not show differentiated gene pools and presented a strong correlation between their patterns of spatial genetic variation. Numerical simulations showed that the striking difference in patterns observed at small and large geographic scales could be accounted for by a combination of (1) isolation by distance within cytotypes; and (2) partial reproductive barriers between cytotypes and/or recurrent formation of tetraploids. We suggest that this may explain the difficulty of the taxonomic treatment of knapweeds and of polyploid complexes in general.

The effect of subdivision on variation at multi-allelic loci under balancing selection.

Simulations are used to investigate the expected pattern of variation at loci under different forms of multi-allelic balancing selection in a finite island model of a subdivided population. The objective is to evaluate the effect of restricted migration among demes on the distribution of polymorphism at the selected loci at equilibrium, and to compare the results with those expected for a neutral locus. The results show that the expected number of alleles maintained, and numbers of nucleotide differences between alleles, are relatively insensitive to the migration rate, and differentiation remains low even under very restricted migration. However, nucleotide divergence between copies of functionally identical alleles increases sharply when migration decreases. These results are discussed in relation to published surveys of allelic diversity in MHC and plant self-incompatibility systems, and to the possibility of inferring ancient population genetic events and processes. In addition, it is shown that, for sporophytic self-incompatibility systems, it is not necessarily true in a subdivided population that recessive alleles are more frequent than dominant ones.

The effect of hitch-hiking on genes linked to a balanced polymorphism in a subdivided population.

The effect of multi-allelic balancing selection on nucleotide diversity at linked neutral sites was investigated by simulations of subdivided populations. The motivation is to understand the behaviour of self-recognition systems such as the MHC and plant self-incompatibility. For neutral sites, two types of subdivision are present: (1) into demes (connected by migration), and (2) into classes defined by different functional alleles at the selected locus (connected by recombination). Previous theoretical studies of each type of subdivision separately have shown that each increases diversity, and decreases the relative frequencies of low-frequency variants, at neutral sites or loci. We show here that the two types of subdivision act non-additively when sampling is at the whole population level, and that subdivision produces some non-intuitive results. For instance, in highly subdivided populations, genetic diversity at neutral sites may decrease with tighter linkage to a selected locus or site. Another conclusion is that, if there is population subdivision, balancing selection leads to decreased expected FST values for neutral sites linked to the selected locus. Finally, we show that the ability to detect balancing selection by its effects on linked variation, using tests such as Tajima's D, is reduced when genes in a subdivided population are sampled from the total population, rather than within demes.

Higher impact of female than male migration on population structure in large mammals.

We simulated large mammal populations using an individual-based stochastic model under various sex-specific migration schemes and life history parameters from the blue whale and the Asian elephant. Our model predicts that genetic structure at nuclear loci is significantly more influenced by female than by male migration. We identified requisite comigration of mother and offspring during gravidity and lactation as the primary cause of this phenomenon. In addition, our model predicts that the common assumption that geographical patterns of mitochondrial DNA (mtDNA) could be translated into female migration rates (Nmf) will cause biased estimates of maternal gene flow when extensive male migration occurs and male mtDNA haplotypes are included in the analysis.

Isolation by distance in a continuous population: reconciliation between spatial autocorrelation analysis and population genetics models.

Analysis of the spatial genetic structure within continuous populations in their natural habitat can reveal acting evolutionary processes. Spatial autocorrelation statistics are often used for this purpose, but their relationships with population genetics models have not been thoroughly established. Moreover, it has been argued that the dependency of these statistics on variation in mutation rates among loci strongly limits their interest for inferential purposes. In the context of an isolation by distance process, we describe relationships between a descriptor of the spatial genetic structure used in empirical studies, Moran's I statistic and population genetics parameters. In particular, we point out that, when Moran's I statistic is used to describe correlation in allele frequencies at the individual level, it provides an estimator of Wright's coefficient of relationship. We also show that the latter parameter, as a descriptor of genetic structure, is not influenced by selfing rate or ploidy level. Under specific finite population models, numerical simulations show that values of Moran's I statistic can be predicted from analytical theory. These simulations are also used to estimate the time taken to approach a structure at equilibrium. Finally, we discuss the conditions under which spatial autocorrelation statistics are little influenced by variation in mutation rates, so that they could be used to estimate gene dispersal parameters.

The signature of balancing selection: fungal mating compatibility gene evolution.

A key problem in evolutionary biology has been distinguishing the contributions of current and historical processes to the maintenance of genetic variation. Because alleles at self-recognition genes are under balancing selection, they exhibit extended residence times in populations and thus may provide unique insight into population demographic history. However, evidence for balancing selection and extended residence times has almost exclusively depended on identification of transspecific polymorphisms; polymorphisms retained in populations through speciation events. We present a broadly applicable approach for detecting balancing selection and apply it to the b1 mating type gene in the mushroom fungus Coprinus cinereus. The comparison of neutral molecular variation within and between allelic classes was used to directly estimate the strength of balancing selection. Different allelic classes are defined as encoding different mating compatibility types and are thus potentially subject to balancing selection. Variation within an allelic class, where all alleles have the same mating compatibility type, provided an internal standard of neutral evolution. Mating compatibility in this organism is determined by the complex A mating type locus, and b1 is one of several redundantly functioning genes. Consequently, we conducted numerical simulations of a model with two subloci and varying levels of recombination to show that balancing selection should operate at each sublocus. Empirical data show that strong balancing selection has indeed occurred at the b1 locus. The widespread geographic distribution of identical b1 alleles suggests that their association with differing A mating types is the result of recent recombination events.