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.

Disentangling the effects of mating systems and mutation rates on cytoplasmic [correction of cytoplamic] diversity in gynodioecious Silene nutans and dioecious Silene otites.

Many flowering plant species exhibit a variety of distinct sexual morphs, the two most common cases being the co-occurrence of females and males (dioecy) or the co-occurrence of hermaphrodites and females (gynodioecy). In this study, we compared DNA sequence variability of the three genomes (nuclear, mitochondrial and chloroplastic) of a gynodioecious species, Silene nutans, with that of a closely related dioecious species, Silene otites. In the light of theoretical models, we expect cytoplasmic diversity to differ between the two species due to the selective dynamics that acts on cytoplasmic genomes in gynodioecious species: under an epidemic scenario, the gynodioecious species is expected to exhibit lower cytoplasmic diversity than the dioecious species, while the opposite is expected in the case of balancing selection maintaining sterility cytoplasms in the gynodioecious species. We found no difference between the species for nuclear gene diversity, but, for the cytoplasmic loci, the gynodioecious S. nutans had more haplotypes, and higher nucleotide diversity, than the dioecious relative, S. otites, even though the latter has a relatively high rate of mitochondrial synonymous substitutions, and therefore presumably a higher mutation rate. Therefore, as the mitochondrial mutation rate cannot account for the higher cytoplasmic diversity found in S. nutans, our findings support the hypothesis that gynodioecy in S. nutans has been maintained by balancing selection rather than by epidemic-like dynamics.

Structural and content diversity of mitochondrial genome in beet: a comparative genomic analysis.

Despite their monophyletic origin, mitochondrial (mt) genomes of plants and animals have developed contrasted evolutionary paths over time. Animal mt genomes are generally small, compact, and exhibit high mutation rates, whereas plant mt genomes exhibit low mutation rates, little compactness, larger sizes, and highly rearranged structures. We present the (nearly) whole sequences of five new mt genomes in the Beta genus: four from Beta vulgaris and one from B. macrocarpa, a sister species belonging to the same Beta section. We pooled our results with two previously sequenced genomes of B. vulgaris and studied genome diversity at the species level with an emphasis on cytoplasmic male-sterilizing (CMS) genomes. We showed that, contrary to what was previously assumed, all three CMS genomes belong to a single sterile lineage. In addition, the CMSs seem to have undergone an acceleration of the rates of substitution and rearrangement. This study suggests that male sterility emergence might have been favored by faster rates of evolution, unless CMS itself caused faster evolution.

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.

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.

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.

Heterozygote deficiencies in small lacustrine populations of brook charr Salvelinus Fontinalis Mitchill (Pisces, Salmonidae): a test of alternative hypotheses.

Empirical studies of natural populations have commonly reported departures from Hardy-Weinberg expected proportions of heterozygote individuals. Recent advances in statistical population genetics now offer the potential to exploit individual multilocus genotypic information to test more rigorously for possible sources of heterozygote deficiencies. In a previous study in lacustrine brook charr (Salvelinus fontinalis), we reported stronger deficits in small than in large lakes. In the present paper, we propose a methodology for empirically testing alternative hypotheses to identify the cause of the deficits observed in three of the smallest lakes (85, 109 and 182 ha) analysed. First, as in several salmonid species, brook charr may exhibit a trophic polymorphism in north temperate lakes. If morphs are genetically divergent, indiscriminate sampling of both forms would result in less heterozygote individuals than expected in a randomly mating population (Wahlund effect). Using an individual-based method aiming at detecting cryptic population structure, we can reject this explanation as the sole source of deficits for all three lakes. Secondly, mating among relatives could also be frequent in small lakes and lead to heterozygote deficiencies. Significantly more fish than expected at random had low individual multilocus heterozygosity in two of the lakes, suggesting that inbred fish may have been present. Thirdly, sampling of genetically related fish would also lead to departures from Hardy-Weinberg proportions. In the same two lakes, the distribution of pairwise relatedness coefficients departed from its random expectation, suggesting that non-random sampling of kin may have occurred.

Landscape structure and hierarchical genetic diversity in the brook charr, Salvelinus fontinalis.

Explaining the extent, causes, and consequences of biotic distributions in space is fundamental to our understanding of how species evolve and cope with particular environments. Yet, identifying extrinsic barriers to migration imposed by landscape structure and predicting their impacts on intraspecific genetic diversity remains a major challenge in population biology. In this study, 30 populations (771 individuals) of brook charr (Salvelinus fontinalis, Salmonidae) representing six major river drainages from Maine, USA, were characterized at six microsatellite loci to quantify the role of landscape features, such as habitat size, altitude, contemporary and historical connectivity. in shaping genetic diversity at three spatial scales: within lakes, within river drainages, and among river drainages. Within-population expected heterozygosity was negatively correlated with altitude, whereas no significant correlation was observed with lake size. Conversely, the extent of heterozygote deficiency within lakes was negatively associated with habitat size. The hierarchical analysis of genetic variance revealed that the extent of among-drainage differentiation was unexpectedly low relative to the pronounced population structuring within drainage. Geographically proximate St. John and Penobscot River drainages were characterized by opposite effects of altitude and geographic distance in shaping the pattern of population differentiation within drainages. The geographic pattern of differentiation among drainages could not be accounted for either by an isolation by distance or by a stepwise range expansion model. Overall, this study provided evidence for the role of contemporary landscape features in shaping the observed pattern of genetic diversity at smaller geographic scales (within and among populations within river drainage). On a broader geographic scale, contemporary landscape structure appeared to be only a minor factor determining the observed pattern of genetic structuring among drainages. These results add to the increasing evidence for nonequilibrium conditions between drift and migration in a wide array of animal taxa. The development of more realistic theoretical descriptions of nonequilibrium population structure thus appears to be important to better understand the relative influence of historical and ecological factors in shaping genetic variation in young habitats, such as recently deglaciated areas.