<title>Introduction</title> <p>The genomes of incipient species diverge at heterogeneous rates, and recently diverged model species are key systems to investigate the causes of this heterogeneity ##REF##11102384##[1]##–##REF##15677748##[3]##. Hybridization followed by introgression between recently diverged plant and animal species with incomplete reproductive barriers is one of the main processes generating the genomic heterogeneity in species divergence ##UREF##1##[4]##. Indeed, some regions appear to be crossing the species barriers more readily than the genomic background (in Helianthus ##REF##17277373##[5]##, Anopheles ##REF##12947038##[6]##, Quercus ##REF##15579711##[7]##, Mytilus ##REF##12535095##[8]##, Mus ##UREF##2##[9]## and Drosophila ##REF##15965264##[10]##). Although much of this heterogeneity may be accounted for by stochasticity of the genetic drift process, natural selection may also play an important role. In particular, because introgressive hybridization brings genetic material from one species into the co-adapted background of another species, some chromosomal fragments are expected to be selected against and resist introgression ##REF##8582621##[11]##.</p> <p>On the other hand, selection can also promote introgression when a transferred chromosome fragment is advantageous in the recipient species. In such a situation, introgression can potentially mediate the transfer of adaptations. Examples of adaptive introgression involving the transfer of transgenes conferring adaptations such as herbicide or insect resistance <italic>via</italic> hybridization with close relatives of crop species ##UREF##3##[12]## have been documented, but other examples in natural populations are strikingly rare ##REF##15004269##[13]##. In the Louisiana Iris species complex for instance, detailed experimental studies provided support for the transfer of adaptations (flood and shade tolerance) between <italic>Iris fulva</italic> and <italic>I. hexagona</italic>
##REF##16415358##[14]##. In Helianthus, a recent experimental study reported that herbivore resistance traits have introgressed from <italic>Heliantus debilis</italic> to <italic>H. annuus,</italic> thereby increasing adaptation of their naturally occurring hybrid <italic>H. annuus taxanus</italic>
##REF##16649157##[15]##. All these documented examples are thus associated with strong directional selection for adaptive traits recently evolved in one of the species and then transmitted horizontally. Theory predicts that adaptive introgression should also be a general property of alleles at genes evolving under multi-allelic balancing selection, such as the vertebrate MHC system, plant disease resistance or self-incompatibility (SI) genes ##REF##11006634##[16]##. In these systems, rare alleles enjoy a strong selective advantage ##REF##17246937##[17]##. Assuming that a given allele is absent from one of two related species, introgression of this allele would then be as strongly favored as a new allele arising by mutation, unless this is impeded by linked genes that are not well adapted to the recipient species. Thus, in multi-allelic systems evolving under balancing selection, repeated exchanges of alleles promoted by adaptive introgression may be expected between closely related species, as long as fertile hybrids can be formed. Therefore, in the course of evolution of strong reproductive isolation between incipient species, such genomic regions should be among the last to stop introgressing.</p> <p>In this study, we test whether multi-allelic balancing selection mediates introgression between closely related species. We do this by contrasting divergence of a portion of the gene controlling self-incompatibility specificity (<italic>SRK</italic>) with the background level of genomic divergence in two closely related plant species. The study system consists of two closely related <italic>Arabidopsis</italic> species, <italic>A. lyrata</italic> and <italic>A. halleri</italic>, whose genomes diverged approximately 2 million years ago ##REF##17404224##[18]##. The two species have overlapping distributions in Northern Europe ##UREF##4##[19]## and relatively recent introgression has been demonstrated for a small fraction of nuclear genes ##REF##15020431##[20]##. SI prevents self-fertilization and some matings among relatives through recognition and rejection of pollen expressing identical specificity. Molecular and genetic analyses of the self-incompatibility locus (S-locus) in <italic>A. lyrata</italic> and <italic>A. halleri</italic> identified many specificities, and the <italic>SRK</italic> sequences often form monophyletic pairs of high sequence similarity, each of which probably represent the same SI specificity in the two species derived from one specificity in their common ancestor. We refer to these pairs as trans-specifically shared pairs of S-alleles. We use divergence at fourfold degenerate sites between alleles within trans-specifically shared pairs to estimate the divergence corresponding to the time of the last introgression event for S-alleles between the two species, and we find that introgression has occurred at a higher rate or continued over more extended periods of time at the S-locus than at the rest of the nuclear genome.</p>
<title>Material and Methods</title> <title>Extent of trans-Specific Allele Sharing at <italic>SRK</italic>
</title> <p>We surveyed sequence diversity at <italic>SRK</italic> in two species-wide samples in <italic>A. halleri</italic> and <italic>A. lyrata</italic> over a total of about 570 nucleotides from the 3′ end of the S-domain using the strategy detailed in ##UREF##6##[31]##. We identified and sequenced five and eight new putative S-alleles in <italic>A. halleri</italic> and <italic>A. lyrata</italic>, respectively. Overall, we analyzed 30 <italic>SRK</italic> sequences in <italic>A. halleri</italic> and 38 sequences in <italic>A. lyrata</italic>. In each case, the nucleotide sequence was obtained as a consensus over three independently obtained sequence products. All identified sequences in <italic>A. halleri</italic> and <italic>A. lyrata</italic> were amino-acid translated and aligned by ClustalW in BioEdit 7.0.5 ##UREF##14##[52]## and adjusted by eye. On the overall set of sequences at <italic>SRK</italic>, we used MEGA 4 ##REF##17488738##[53]## to reconstruct a phylogeny using the Neighbor-Joining method based on the total number of differences per site or on the number of either synonymous or non-synonymous differences.</p> <title>Divergence within trans-Specifically Shared Pairs of S-Alleles</title> <p>Within each pair of trans-specifically shared sequences at <italic>SRK</italic>, we estimated the number of synonymous nucleotide differences per synonymous site between the <italic>A. halleri</italic> and the <italic>A. lyrata</italic> copy using the method of ##REF##3444411##[54]## with MEGA 4. A homogeneous substitution process across all pairs is expected to result in an accumulation of nucleotide differences according to the Poisson distribution. We used Fisher's dispersion index to test whether the distribution of nucleotide differences across trans-specifically shared sequence pairs could be explained by the stochasticity of the substitution process alone. We used Fisher's exact test of independence to test whether synonymous and non-synonymous differences hit HV regions equally frequently.</p> <title>Inference on Introgression Patterns at the S-Locus</title> <p>Background genomic divergence was estimated by the species-wide average nucleotide divergence at fourfold degenerate sites (<italic>K</italic>
<sub>4fold</sub>) between the two species for 12 reference genes that had been previously sequenced ##REF##15020431##[20]##,##UREF##9##[39]##,##REF##12694288##[40]## and two genes that are members of the S-domain gene family (<italic>Aly10.1, Aly10.2</italic>).</p> <p>To determine whether difference in effective population size and thus coalescence time between S-alleles and genomic background may suffice to explain the low divergence of S-alleles, we applied the isolation with migration model of Nielsen and Wakeley ##REF##11404349##[38]## to polymorphism at fourfold degenerate sites in both species for the eleven reference genes plus Aly9 (12 genes in total, see ##TAB##1##table 2##) as implemented in the IM program ##REF##11404349##[38]##. We chose to focus on fourfold degenerate sites only because differences in substitution rates have been reported among codon positions ##UREF##15##[55]##. The program DNAsp ##REF##14668244##[56]## was used to generate a datafile containing fourfold degenerate sites only. The procedure was run with 10 different random seeds to ensure proper convergence of the six free parameters, <italic>i.e. θ<sub>A</sub></italic>, <italic>θ<sub>lyrata</sub></italic>, <italic>θ<sub>halleri</sub>, t, m</italic>
<sub>hal→lyr</sub>, <italic>m</italic>
<sub>lyr→hal</sub> (polymorphism <italic>θ</italic> = 4<italic>Nμ</italic> in the common ancestor of <italic>A. halleri</italic> and <italic>A. lyrata</italic>, polymorphism in <italic>A. lyrata</italic>, polymorphism in <italic>A. halleri</italic>, splitting time and the rate at which genes introgressed into <italic>A. lyrata</italic> from <italic>A. halleri</italic> and into <italic>A. halleri</italic> from <italic>A. lyrata</italic> as time moves forward, respectively). The HKY mutation model ##REF##3934395##[57]## was used. To single out the <italic>N<sub>A</sub></italic> estimate, we estimated the average per fourfold degenerate site mutation rate (<italic>μ</italic>) as follows. We used <italic>A. thaliana</italic> as outgroup to estimate the average net nucleotide divergence at fourfold degenerate sites between <italic>A. thaliana</italic> and <italic>A. halleri</italic> and between <italic>A. thaliana</italic> and <italic>A. lyrata</italic> for each reference gene. Assuming that the lineages leading to <italic>A. thaliana</italic> and the common ancestor of <italic>A. lyrata</italic> and <italic>A. halleri</italic> separated 5 million years ago ##REF##11018155##[58]##, we obtained a mutation rate estimate per site per year for each reference gene. We computed an average mutation rate per site per year (<italic>μ</italic>) by taking the geometric mean over genes. A mutation rate per generation was computed assuming a mean generation time of two years.</p> <p>The maximum likelihood estimates were then used to simulate divergence between two species isolated since one million generations but still capable of introgression. Ten thousand replicates of pairs of genes with the same number of nucleotides as the real data were performed using SIMCOAL2 ##REF##15117750##[59]##. The genomic background divergence was first used to confirm that the simulations parameters were appropriate. We then determined whether the observed divergence for S-alleles was consistent with the overall genomic rate of introgression by simulating the evolution of S-alleles in this system assuming that 50 S-alleles segregate in the species, and thus that the effective population size of each allelic class is reduced by a factor 50. To remain conservative in this analysis, S-alleles were simulated under the 2.5% low boundary of the 95% credible interval for <italic>N<sub>A</sub></italic> obtained from IM.</p> <p>Using the maximum likelihood estimate for <italic>N<sub>A</sub></italic>, we then aimed to determine the extent to which introgression is increased for S-alleles relative to the genomic background. We did so by gradually increasing <italic>m</italic>
<sub>hal→lyr</sub> and <italic>m</italic>
<sub>lyr→hal</sub> for S-alleles by a multiplicative factor from one to ten until the simulated data came close to the observed divergence.</p> <p>The sequences reported in this paper have been deposited in the GenBank database under accession numbers EU878008- EU878026.</p>
<title>Results</title> <title>Extent of trans-Specific Allele Sharing at <italic>SRK</italic>
</title> <p>Our species-wide survey of sequence diversity reveals that a large fraction of alleles at the pistil self-incompatibility specificity-determining gene <italic>SRK</italic> (S-locus receptor kinase) are trans-specifically shared between the two species (##FIG##0##Figure 1##). Overall, we find 30 sets of <italic>SRK</italic> sequences in <italic>A. halleri</italic> and 38 sets of <italic>SRK</italic> sequences in <italic>A. lyrata.</italic> As is typical for S-alleles ##REF##7982569##[21]##, the sequences fall into sets of nearly identical ones (presumably representing the same specificity, ##REF##7982569##[21]##–##REF##12832651##[23]##) and ones with many differences from all other sequences (presumably representing functionally distinct specificities), with the most similar pairs within <italic>A. halleri</italic> and <italic>A. lyrata</italic> showing 44 and 51 differences, respectively, over a total of about 570 nucleotides. We then compared nucleotide sequences between S-alleles from the two species and find that the mismatch distribution (##FIG##1##Figure 2##) is clearly bimodal. Most comparisons are in line with intraspecific comparisons and range between 45 and 218 differences over a total of about 570 nucleotides (see also ##SUPPL##4##Figure S5##), but the distribution shows a distinct set of 18 highly similar interspecific pairs of sequences (indicated by brackets in ##FIG##0##Figure 1##) with at most 12 nucleotide differences. The numbers of non-synonymous differences within the 18 highly similar pairs of S-alleles ranged from 0 to 9 over a total of 380 non-synonymous sites. These sequences are more similar than pairs of alleles known to have retained the same specificity when comparing the closely related <italic>Brassica oleracea</italic> and <italic>B. rapa</italic>
##UREF##5##[24]##–##REF##16381659##[27]##. Even if these sequences currently occur in two different (but closely related) species, we therefore hypothesize that these pairs have retained identical specificity. We refer to these 18 pairs of S-alleles as “trans-specifically shared” pairs of alleles and note that they represent 60% and 47% of S-alleles found to date in <italic>A. halleri</italic> and <italic>A. lyrata</italic>, respectively. Two of these pairs (<italic>AlSRK37</italic>/<italic>AhSRK04</italic> and <italic>AlSRK16</italic>/<italic>AhSRK10</italic>) were previously identified and shown additionally to be shared trans-specifically with <italic>A. thaliana</italic>
##REF##16782760##[28]##. Phylogenetic reconstructions show that both synonymous (##SUPPL##1##Figure S2##) and non-synonymous (##SUPPL##2##Figure S3##) differences are strikingly low within trans-specifically shared pairs and high among pairs. Note that, by definition, <italic>SRK</italic> alleles we consider as trans-specific pairs are determined based on those S allele pairs that have the fewest differences, so the procedure could potentially lead to ascertainment bias. Yet, close examination of the next best candidates (<italic>AhSRK03/AlSRK28</italic>, <italic>AhSRK28/AlSRK03</italic>, <italic>AhSRK23/AlSRK06</italic> and <italic>AhSRK20/AlSRK04,</italic>
##FIG##0##Figure 1##) suggests that none of these pairs is likely to represent pairs of trans-specifically shared alleles (detailed arguments are presented in T<sc>ext</sc> S1).</p> <title>Divergence within trans-Specifically Shared Pairs of S-Alleles</title> <p>Within <italic>SRK</italic>, several hypervariable (HV) regions have been identified in the domain responsible for binding the pollen protein (S-domain) and shown to be targets of positive selection, suggesting they are involved in determination of specificity ##REF##12399400##[29]##,##REF##17609367##[30]##. Accordingly, HV regions from different specificities within species typically show an excess of non-synonymous substitutions ##REF##12399400##[29]##,##UREF##6##[31]##,##REF##17683611##[32]##. In sharp contrast, we find that as compared to synonymous differences, non-synonymous differences are relatively less frequent in HV regions than in non-HV regions (on average 0.7 and 2.3 differences in HV and non-HV regions respectively for non-synonymous differences, <italic>versus</italic> 1.1 and 1.6 differences respectively for synonymous differences, ##TAB##0##Table 1##). This contrast is significant by Fisher's exact test of independence (odds ratio = 2.5, <italic>p</italic> = 0.029), suggesting that sequence pairs that putatively encode the same specificity tend to have similar HV region sequences for non-synonymous sites, but might differ at synonymous sites in these regions, whereas other regions may differ at both types of sites.</p> <p>If introgression occurs, then divergence might also be affected by the dominance of the S-alleles. Indeed, complex patterns of dominance relationships generally occur among alleles in sporophytic SI systems ##REF##14659710##[33]## and Billiard et al. ##REF##17237502##[34]## reported asymmetric selective pressures for dominant and recessive S-alleles because rare dominant S-alleles will tend to express their specificity more often than rare recessive ones (a process similar to “Haldane's sieve”-the bias against the establishment of recessive beneficial mutations ##UREF##7##[35]##,##UREF##8##[36]##). Hypothesizing that the introgression rate thus differs between dominant and recessive S-alleles, we tested for an effect of dominance on divergence between the two species. The range of variation observed for nucleotide differences across pairs of trans-specifically shared S-alleles cannot be explained fully by the stochasticity of the substitution process (Fisher's dispersion index = 2.03, P = 0.0103), but there was no obvious relationship between number of nucleotide differences and level of dominance of the S-alleles, as inferred from the phylogeny of alleles as suggested by ##REF##15858208##[37]##. Thus, we find no evidence that dominance affects S-allele divergence between the two species.</p> <title>Comparison of Introgression Rate between S-Locus and Genomic Background</title> <p>To test whether balancing selection resulted in adaptive introgression of S-alleles between the two species, we compared levels of divergence at fourfold degenerate sites between trans-specifically shared S-alleles with that of the genomic background, estimated from twelve unlinked control genes and two S-gene family members. These two sets of control genes give similar mean values of divergence (<italic>K</italic>
<sub>4fold</sub> = 0.0743 and <italic>K</italic>
<sub>4fold</sub> = 0.0904, respectively, ##TAB##1##Table 2##), which are about four times higher than the average for trans-specifically shared pairs of S-alleles (<italic>K</italic>
<sub>4fold</sub> = 0.0193, ##TAB##0##Table 1##).</p> <p>Because a large number of S-alleles are actively maintained within species by balancing selection, each S-allele has individually a small effective population size ##REF##7982569##[21]##. Thus, estimates of divergence for S-alleles and reference genes cannot be compared directly because of differences in effective population sizes (##FIG##2##Figure 3##). To take this into account, we used coalescent simulations to test whether our data are compatible with a null model of speciation (the “isolation with migration” model of Nielsen and Wakeley, ##REF##11404349##[38]##) that assumes identical introgression rate for S-alleles and the genomic background. Under this model, we first used previously published species-wide polymorphism data in <italic>A. halleri</italic> and <italic>A. lyrata</italic> from ##REF##15020431##[20]##,##UREF##9##[39]##,##REF##12694288##[40]## to estimate rates of introgression, splitting time <italic>t</italic> as well as <italic>θ<sub>A</sub></italic> = 4<italic>N<sub>A</sub>μ</italic>, where <italic>N<sub>A</sub></italic> is the effective population size in their common ancestor and <italic>μ</italic> the substitution rate. The maximum likelihood estimates for directional rates of introgression are <italic>m<sub>hal→lyr</sub></italic> = 2.775×10<sup>−7</sup>, <italic>m<sub>lyr→hal</sub></italic> = 2.912×10<sup>−7</sup>, and <italic>θ<sub>A</sub></italic> = 1.7975 (##TAB##2##Table 3##). The <italic>t</italic> estimate is 2,533,980 years [1,307,952–5,166,833], which is entirely consistent with the previous 2 Myrs estimate by Koch &amp; Matschinger ##REF##17404224##[18]##. All estimates converge satisfactorily based on 10 replicate runs with different random seeds. To single out the <italic>N<sub>A</sub></italic> estimate, we then used <italic>A. thaliana</italic> as outgroup to obtain a substitution rate at fourfold degenerate sites of <italic>μ</italic> = 1.296×10<sup>−8</sup> substitutions per nucleotide per year ([9.218×10<sup>−9</sup>–1.781×10<sup>−8</sup>] as 95% credible interval). The resulting estimate for <italic>N<sub>A</sub></italic> is 253,892 with [13,772–663,510] as 95% credible interval. Based on these parameters, we then simulated the evolution of two species exchanging migrants at the rate estimated above. The simulations were entirely consistent with the data for the genomic background (<italic>K</italic> = 0.0678 [0.0423–0.0955], ##FIG##3##Figure 4##). In sharp contrast, conservatively assuming a reduction of effective population size for S-alleles by a factor 50 (as expected if 50 different S-alleles segregate in each species) only led to a modest reduction in divergence (<italic>K</italic> = 0.0465, ##FIG##3##Figure 4##), whose 95% credible interval [0.0305–0.0640] did not comprise the observed value for <italic>K</italic> (<italic>K</italic>
<sub>4fold = </sub>0.0193). Hence, the data are not consistent with equal introgression for S-alleles and the genomic background. This result is robust to the conservative use of the lower boundary of the 95%CI for either <italic>N</italic>
<sub>A</sub> or <italic>t</italic>. Increasing the rates of introgression for S-alleles led to a sharp reduction in divergence between <italic>A. halleri</italic> and <italic>A. lyrata.</italic> The simulations best fitted the data when the directional rates of introgression were empirically increased for S-alleles by a factor 5, with divergence value closely approaching the observed data (<italic>K</italic> = 0.0182, ##FIG##3##Figure 4##). A simpler analysis also confirmed that average net interspecific divergence ##UREF##10##[41]## for S-alleles was lower than that at the genomic background (##SUPPL##6##Text S1##).</p> <p>For three pairs of S-alleles (<italic>AlSRK01/AhSRK01</italic>, <italic>AlSRK34/AhSRK05</italic>, <italic>AlSRK37/AhSRK04</italic>) we also surveyed intra-allelic variation in at least 10 copies from each species. We found very little diversity among allelic copies within each surveyed allele in each species (average synonymous diversity = 0.0064, data not shown) in accordance with their low expected effective population sizes. We examined the sequences for shared polymorphisms, and found none in any of these S-allele pairs. This suggests old and infrequent, rather than recent, introgression events since the separation of <italic>A. lyrata</italic> and <italic>A. halleri</italic>. Moreover, the estimated divergence among pairs of S-alleles was more heterogeneous than expected based on the Poisson distribution, suggesting that the last introgression event occurred at different times for different alleles.</p>
<title>Discussion</title> <title>Impact of Founder Events at Speciation</title> <p>The possibility of introgression of S-alleles may have important consequences for the extent of allelic diversity maintained within self-incompatible species. If introgression occurs, hybridizing species effectively share a common pool of S-alleles. If hybridization is restricted, the two species together can maintain more S-alleles than each species individually ##REF##9611223##[42]##. Such a process could be especially important in the first stages of the split because reproductive barriers may then be more leaky, and also because allelic diversity at the S-locus within incipient species may be decreased if founding events were associated with speciation. This process could be responsible for maintaining many highly divergent allelic lineages at the S locus within plant families, where trans-generic sharing of allelic lineages seems to be the rule, and loss of ancestral allelic lineages through strong bottlenecks within particular genera the exception, as has been described in the Solanaceae ##UREF##11##[43]##.</p> <p>It can therefore be misleading to use a species' extant number of lineages at a gene under balancing selection to estimate the minimum population size at speciation. For instance, using polymorphism data for MHC in humans, Takahata ##REF##2320564##[44]## predicted that the number of breeding individuals in the human lineage could not be as small as 50-100 at any time of its evolutionary history, assuming two extant ancestral allelic lineages at HLA-B. According to our hypothesis of adaptive introgression mediated by balancing selection, variation can be efficiently “rescued”, and stronger founder events at speciation would still be compatible with extant variation at HLA-B, if some interbreeding occurred with the chimpanzee lineage after the split. Although identifying the functional types of alleles may not be simple in that case (and recombination may confine the effect of balancing selection to a small region around the selected sites themselves), a detailed analysis of MHC alleles in the great apes would be of great interest to survey whether adaptive introgression mediated by balancing selection has indeed occurred in primates.</p> <title>Shared Chloroplast Haplotypes: Distinguishing between Introgression and Ancestral Polymorphism</title> <p>A recent study by Koch and Matschinger ##REF##17404224##[18]## reported that, whereas <italic>A. lyrata</italic> and <italic>A. halleri</italic> were well separated in phylogenetic trees based on the nuclear encoded ITS region, several cpDNA haplotypes are shared between both species ##REF##17404224##[18]##. This was interpreted as ancestral polymorphism segregating for the chloroplast but not the nucleus. However, this interpretation is at odds with the smaller effective population size expected for the chloroplast (approximately 1/2 for hermaphroditic species, ##REF##6840539##[45]##) and the consequent low expected variability. Indeed most studies in plants have found low sequence diversity for chloroplast genes, taking into account their low mutation rate ##REF##3480529##[46]##, and also stronger differentiation among populations for chloroplast than nuclear markers ##UREF##12##[47]##. In line with our results from S-alleles, we suggest the alternative interpretation that introgression occurred more readily for the chloroplast than nuclear genes, as has been reported in several instances (e.g. ##UREF##13##[48]##,##REF##16780430##[49]##). The haplotype network of chloroplast sequences reported by Koch and Matschinger ##REF##17404224##[18]## also showed greater sharing of more basal haplotypes, suggesting that chloroplast introgression has become less common in recent times.</p> <title>Evolution of New Specificities of Self-Incompatibility Genes</title> <p>Our results also shed light on the evolution of self-incompatibility specificities. Indeed, our data strongly suggest that purifying selection prevents the substitution of non-synonymous differences within HV regions, supporting a role for these regions in determining specificity. More specifically, the strength of purifying selection seems higher on the HV regions than on the rest of the sequence, and this could be related to strong selection against mutations altering specificities. Mechanisms selecting against mutant S-alleles with altered pistil specificities have been discussed by Uyenoyama et al. ##REF##11290732##[50]##.</p> <p>Inter-species exchanges of S-alleles may, however, be important in the evolution of new specificities. Chookajorn et al. ##REF##14694189##[51]## suggested that new specificities could evolve if sufficient variation could be maintained within the pollen (or pistil) S gene for enough time to allow variants of the other gene to co-evolve with them. Due to the small effective population size of individual S-alleles, this hypothesis requires population structure with very limited migration ##REF##11006634##[16]##. Speciation with some introgression of S-alleles leads to precisely the strongly subdivided population needed for this mechanism to work. Under this hypothesis, two alleles could slowly evolve to different specificities in two isolated species and then add to the number of S-alleles in each species after reciprocal introgression. Data testing the specificities of sequence pairs in the two species that differ at few amino acids might help determine whether new specificities have indeed arisen in one species or the other since they split.</p>
<p>Conceived and designed the experiments: VC JB MHS XV. Performed the experiments: VC JB. Analyzed the data: VC JB MHS XV. Contributed reagents/materials/analysis tools: MHS XV. Wrote the paper: VC JB MHS XV.</p> <p>Recently diverged species typically have incomplete reproductive barriers, allowing introgression of genetic material from one species into the genomic background of the other. The role of natural selection in preventing or promoting introgression remains contentious. Because of genomic co-adaptation, some chromosomal fragments are expected to be selected against in the new background and resist introgression. In contrast, natural selection should favor introgression for alleles at genes evolving under multi-allelic balancing selection, such as the MHC in vertebrates, disease resistance, or self-incompatibility genes in plants. Here, we test the prediction that negative, frequency-dependent selection on alleles at the multi-allelic gene controlling pistil self-incompatibility specificity in two closely related species, <italic>Arabidopsis halleri</italic> and <italic>A. lyrata</italic>, caused introgression at this locus at a higher rate than the genomic background. Polymorphism at this gene is largely shared, and we have identified 18 pairs of S-alleles that are only slightly divergent between the two species. For these pairs of S-alleles, divergence at four-fold degenerate sites (<italic>K = </italic>0.0193) is about four times lower than the genomic background (<italic>K = </italic>0.0743). We demonstrate that this difference cannot be explained by differences in effective population size between the two types of loci. Rather, our data are most consistent with a five-fold increase of introgression rates for S-alleles as compared to the genomic background, making this study the first documented example of adaptive introgression facilitated by balancing selection. We suggest that this process plays an important role in the maintenance of high allelic diversity and divergence at the S-locus in flowering plant families. Because genes under balancing selection are expected to be among the last to stop introgressing, their comparison in closely related species provides a lower-bound estimate of the time since the species stopped forming fertile hybrids, thereby complementing the average portrait of divergence between species provided by genomic data.</p> <title>Author Summary</title> <p>The role of natural selection in promoting or preventing genomic divergence between nascent species remains highly debated. As long as reproductive barriers remain incomplete, genetic material from one species is indeed exposed to natural selection into the genomic background of the other species. In some cases, genomic co-adaptations developing independently in each species are believed to select against such transfers. Yet, theory predicts that the transfer of some chromosomal fragments may be favored by natural selection. In particular, this should occur for alleles at genes evolving under a particular form of natural selection, i.e., multi-allelic balancing selection. We test this prediction using two closely related <italic>Arabidopsis</italic> species, and find a four-fold lower divergence at alleles at the gene controlling pistil self-incompatibility specificity than at the genomic background. We conclude that alleles at this gene have been transferred more readily between the two species than the genomic background. We suggest that natural selection may efficiently allow the maintenance of high allelic diversity and divergence across many species at S-loci as well as at all other loci under multi-allelic balancing selection, such as the MHC in vertebrates or disease resistance genes in plants.</p>
<title>Supporting Information</title>
<p>We thank Violaine Llaurens and Maria Valeria Ruggiero for providing sequences from several S-alleles in <italic>A. halleri</italic>; Aude Darracq for a Python script; Jody Hey for advice on IM analyses; Freddy B. Christiansen, Nicolas Bierne and Thomas Bataillon for discussions and Christian R. Landry, Sylvain Billiard and Deborah Charlesworth for detailed comments on the manuscript.</p>
<fig id="pgen-1000168-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000168.g001</object-id><label>Figure 1</label><caption><title>Phylogeny of the 68 <italic>SRK</italic> sequences of <italic>A. lyrata</italic> and <italic>A. halleri.</italic>
</title><p>The phylogeny was obtained by the neighbour-joining method on pairwise proportion of nucleotide divergence after Jukes-Cantor's correction. Brackets indicate interspecific pairs of sequences assumed to represent “trans-specifically shared S-alleles”, <italic>i.e.</italic> alleles assumed to have evolved from a single S-allele in the direct ancestor of <italic>A. lyrata</italic> and <italic>A. halleri</italic>.</p></caption></fig> <fig id="pgen-1000168-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000168.g002</object-id><label>Figure 2</label><caption><title>Distribution of the number of pairwise nucleotide differences for <italic>SRK</italic> sequences in interspecific comparisons between <italic>A. halleri</italic> and <italic>A. lyrata</italic>.</title><p>Note the distinct peak of highly similar sequences observed. The vertical arrow represents the chosen threshold to define “trans-specifically shared” pairs of sequences (≤12 nucleotide differences). Note also that the two pairs of sequences with intermediate nucleotide differences (45 between <italic>AlSRK03</italic> and <italic>AhSRK28</italic>, and 50 between <italic>AlSRK28</italic> and <italic>AhSRK03</italic>) cannot represent trans-specifically shared S-alleles because they are not monophyletic (see ##FIG##0##Figure 1##).</p></caption></fig> <fig id="pgen-1000168-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000168.g003</object-id><label>Figure 3</label><caption><title>Divergence process between <italic>Arabidopsis lyrata</italic>, <italic>A. halleri</italic> and <italic>A. thaliana</italic> at unlinked genes (genomic background) and trans-specifically shared pairs of S-alleles.</title><p>Divergence times were taken from Koch et al. ##REF##17404224##[18]##,##REF##11018155##[58]##. <italic>θ<sub>lyrata</sub></italic>, <italic>θ<sub>halleri</sub></italic> and <italic>θ<sub>A</sub></italic>, refer to polymorphism (<italic>θ = </italic>4<italic>Nμ·</italic> in <italic>A. lyrata</italic>, <italic>A. halleri</italic> and their common ancestor. As compared to unlinked genes, divergence between trans-specifically shared S-alleles is affected by two confounding factors: (1) lower effective population size than the genomic background reducing coalescence time in the ancestral species, and (2) expected higher introgression as represented by thicker dark grey arrows.</p></caption></fig> <fig id="pgen-1000168-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000168.g004</object-id><label>Figure 4</label><caption><title>Predicted nucleotide divergence between <italic>A. halleri</italic> and <italic>A. lyrata</italic> for the genomic background (grey line), S-alleles with the same rate of introgression as the genomic background (dotted line) and S-alleles with 5-fold increased rate of introgression relative to the genomic background (black line).</title><p>10,000 coalescent simulations were performed for each case using maximum likelihood parameter estimates obtained under the “isolation with migration” model, except for the dotted line, where the 2.5% low ancestral population size estimate was used in order to be conservative. Observed nucleotide divergence for the genomic background and S-alleles are represented by grey and black stars on the x-axis, respectively.</p></caption></fig>
<table-wrap id="pgen-1000168-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000168.t001</object-id><label>Table 1</label><caption><title>Divergence between <italic>Arabidopsis halleri</italic> and <italic>A. lyrata</italic> at trans-specifically shared <italic>SRK</italic> alleles at synonymous (<italic>K</italic>
<sub>S</sub>), non synonymous (<italic>K</italic>
<sub>A</sub>) and fourfold degenerate sites (<italic>K</italic>
<sub>4fold</sub>).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Coding sequence length</td><td align="left" rowspan="1" colspan="1">Nucleotides in HV</td><td colspan="2" align="left" rowspan="1">Synonymous differences</td><td colspan="2" align="left" rowspan="1">Non-synonymous differences</td><td colspan="3" align="left" rowspan="1">Nucleotide divergence</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">in HV</td><td align="left" rowspan="1" colspan="1">not in HV</td><td align="left" rowspan="1" colspan="1">in HV</td><td align="left" rowspan="1" colspan="1">not in HV</td><td align="left" rowspan="1" colspan="1">
<italic>K</italic>
<sub>S</sub>
</td><td align="left" rowspan="1" colspan="1">
<italic>K</italic>
<sub>A</sub>
</td><td align="left" rowspan="1" colspan="1">
<italic>K</italic>
<sub>4fold</sub>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">AlSRK01</td><td align="left" rowspan="1" colspan="1">AhSRK01</td><td align="left" rowspan="1" colspan="1">578</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.0080</td><td align="left" rowspan="1" colspan="1">0.0045</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK03</td><td align="left" rowspan="1" colspan="1">AhSRK03</td><td align="left" rowspan="1" colspan="1">598</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.0021</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK08</td><td align="left" rowspan="1" colspan="1">AhSRK19</td><td align="left" rowspan="1" colspan="1">593</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">0.0581</td><td align="left" rowspan="1" colspan="1">0.0086</td><td align="left" rowspan="1" colspan="1">0.0969</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK11</td><td align="left" rowspan="1" colspan="1">AhSRK11</td><td align="left" rowspan="1" colspan="1">565</td><td align="left" rowspan="1" colspan="1">122</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.0241</td><td align="left" rowspan="1" colspan="1">0.0185</td><td align="left" rowspan="1" colspan="1">0.0423</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK13</td><td align="left" rowspan="1" colspan="1">AhSRK29</td><td align="left" rowspan="1" colspan="1">537</td><td align="left" rowspan="1" colspan="1">91</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.0444</td><td align="left" rowspan="1" colspan="1">0.0072</td><td align="left" rowspan="1" colspan="1">0.0330</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK14</td><td align="left" rowspan="1" colspan="1">AhSRK09</td><td align="left" rowspan="1" colspan="1">598</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.0064</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK15</td><td align="left" rowspan="1" colspan="1">AhSRK21</td><td align="left" rowspan="1" colspan="1">590</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.0331</td><td align="left" rowspan="1" colspan="1">0.0022</td><td align="left" rowspan="1" colspan="1">0.0163</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK16</td><td align="left" rowspan="1" colspan="1">AhSRK10</td><td align="left" rowspan="1" colspan="1">592</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.0158</td><td align="left" rowspan="1" colspan="1">0.0022</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK17</td><td align="left" rowspan="1" colspan="1">AhSRK02</td><td align="left" rowspan="1" colspan="1">557</td><td align="left" rowspan="1" colspan="1">118</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">0.0251</td><td align="left" rowspan="1" colspan="1">0.0093</td><td align="left" rowspan="1" colspan="1">0.0160</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK20</td><td align="left" rowspan="1" colspan="1">AhSRK24</td><td align="left" rowspan="1" colspan="1">551</td><td align="left" rowspan="1" colspan="1">106</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">0.0446</td><td align="left" rowspan="1" colspan="1">0.0093</td><td align="left" rowspan="1" colspan="1">0.0347</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK21</td><td align="left" rowspan="1" colspan="1">AhSRK06</td><td align="left" rowspan="1" colspan="1">516</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.0260</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK22</td><td align="left" rowspan="1" colspan="1">AhSRK26</td><td align="left" rowspan="1" colspan="1">568</td><td align="left" rowspan="1" colspan="1">122</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.0081</td><td align="left" rowspan="1" colspan="1">0.0137</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK28</td><td align="left" rowspan="1" colspan="1">AhSRK28</td><td align="left" rowspan="1" colspan="1">548</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.0258</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.0300</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK31</td><td align="left" rowspan="1" colspan="1">AhSRK16</td><td align="left" rowspan="1" colspan="1">509</td><td align="left" rowspan="1" colspan="1">91</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.0194</td><td align="left" rowspan="1" colspan="1">0.0075</td><td align="left" rowspan="1" colspan="1">0.0174</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK34</td><td align="left" rowspan="1" colspan="1">AhSRK05</td><td align="left" rowspan="1" colspan="1">539</td><td align="left" rowspan="1" colspan="1">111</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.0084</td><td align="left" rowspan="1" colspan="1">0.0048</td><td align="left" rowspan="1" colspan="1">0.0163</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK37</td><td align="left" rowspan="1" colspan="1">AhSRK04</td><td align="left" rowspan="1" colspan="1">592</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.0315</td><td align="left" rowspan="1" colspan="1">0.0087</td><td align="left" rowspan="1" colspan="1">0.0287</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK39</td><td align="left" rowspan="1" colspan="1">AhSRK18</td><td align="left" rowspan="1" colspan="1">573</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">AlSRK42</td><td align="left" rowspan="1" colspan="1">AhSRK12</td><td align="left" rowspan="1" colspan="1">552</td><td align="left" rowspan="1" colspan="1">124</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.0261</td><td align="left" rowspan="1" colspan="1">0.0116</td><td align="left" rowspan="1" colspan="1">0.0163</td></tr><tr><td colspan="2" align="left" rowspan="1">
<bold>Average</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>569</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>118</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>1.1</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>1.6</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>0.7</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>2.3</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>0.0221</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>0.0065</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>0.0193</bold>
</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000168-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000168.t002</object-id><label>Table 2</label><caption><title>Divergence between <italic>Arabidopsis halleri</italic> and <italic>A. lyrata</italic> at reference genes and members of the S-gene family at fourfold degenerate sites (<italic>K</italic>
<sub>4fold</sub>).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Coding sequence length</td><td colspan="2" align="left" rowspan="1">Number of sequences analysed</td><td align="left" rowspan="1" colspan="1">
<italic>K</italic>
<sub>4fold</sub>
</td><td align="left" rowspan="1" colspan="1">References</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>A. lyrata</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>A. halleri</italic>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Genomic background</td><td align="left" rowspan="1" colspan="1">CAD</td><td align="left" rowspan="1" colspan="1">956</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">0.0483</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">CHI</td><td align="left" rowspan="1" colspan="1">264</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">0.1468</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">CHS</td><td align="left" rowspan="1" colspan="1">1177</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">0.0881</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">DFR</td><td align="left" rowspan="1" colspan="1">346</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">0.0610</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">F3H</td><td align="left" rowspan="1" colspan="1">450</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">0.1223</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">FAH1</td><td align="left" rowspan="1" colspan="1">1054</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">0.0739</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">GS</td><td align="left" rowspan="1" colspan="1">906</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">0.0861</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">MAML</td><td align="left" rowspan="1" colspan="1">388</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">0.0312</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">CAUL</td><td align="left" rowspan="1" colspan="1">246</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1">0.0184</td><td align="left" rowspan="1" colspan="1">39, 40</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">HAT4</td><td align="left" rowspan="1" colspan="1">340</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">34</td><td align="left" rowspan="1" colspan="1">0.0531</td><td align="left" rowspan="1" colspan="1">39, 40</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">ScADH</td><td align="left" rowspan="1" colspan="1">515</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1">34</td><td align="left" rowspan="1" colspan="1">0.0323</td><td align="left" rowspan="1" colspan="1">39, 40</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Aly9</td><td align="left" rowspan="1" colspan="1">443</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">28</td><td align="left" rowspan="1" colspan="1">0.1296</td><td align="left" rowspan="1" colspan="1">39, 40</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>Average</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>0.0743</bold>
</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">S-gene family</td><td align="left" rowspan="1" colspan="1">Aly10.1</td><td align="left" rowspan="1" colspan="1">936</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0.1043</td><td align="left" rowspan="1" colspan="1">this study</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Aly10.2</td><td align="left" rowspan="1" colspan="1">466</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">0.0765</td><td align="left" rowspan="1" colspan="1">this study</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>Average</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>0.0904</bold>
</td><td align="left" rowspan="1" colspan="1"/></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000168-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000168.t003</object-id><label>Table 3</label><caption><title>Estimates of <italic>θ = </italic>4<italic>Nμ</italic>, effective population sizes, splitting time and rates of introgression using the isolation with migration model ##REF##11404349##[38]##.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Parameter</td><td align="left" rowspan="1" colspan="1">ML estimate</td><td align="left" rowspan="1" colspan="1">95% CI</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Common ancestor</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> θ<sub>A</sub></italic>
</td><td align="left" rowspan="1" colspan="1">1.7975</td><td align="left" rowspan="1" colspan="1">0.0975–4.6975</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> N<sub>A</sub></italic>
</td><td align="left" rowspan="1" colspan="1">253,892</td><td align="left" rowspan="1" colspan="1">13,772–663,510</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>A. lyrata</italic>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> θ<sub>lyrata</sub></italic>
</td><td align="left" rowspan="1" colspan="1">1.2225</td><td align="left" rowspan="1" colspan="1">0.7635–1.8405</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> N<sub>lyrata</sub></italic>
</td><td align="left" rowspan="1" colspan="1">172,675</td><td align="left" rowspan="1" colspan="1">107,842–259,966</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>A. halleri</italic>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> θ<sub>halleri</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.7785</td><td align="left" rowspan="1" colspan="1">0.4635–1.2195</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> N<sub>halleri</sub></italic>
</td><td align="left" rowspan="1" colspan="1">109,961</td><td align="left" rowspan="1" colspan="1">65,468–172,251</td></tr><tr><td align="left" rowspan="1" colspan="1">Splitting time</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> t</italic> (years)</td><td align="left" rowspan="1" colspan="1">2,533,980</td><td align="left" rowspan="1" colspan="1">1,307,952–5,166,833</td></tr><tr><td align="left" rowspan="1" colspan="1">Rates of introgression</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> m<sub>hal→lyr</sub></italic>
</td><td align="left" rowspan="1" colspan="1">2.775×10<sup>−7</sup>
</td><td align="left" rowspan="1" colspan="1">5.186×10<sup>−7</sup>–7.510×10<sup>−7</sup>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> m<sub>lyr→hal</sub></italic>
</td><td align="left" rowspan="1" colspan="1">2.912 ×10<sup>−7</sup>
</td><td align="left" rowspan="1" colspan="1">2.035×10<sup>−7</sup>–1.059×10<sup>−7</sup>
</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pgen.1000168.s001"><label>Figure S1</label><caption><p>Phylogeny of SRK sequences from the species <italic>A. lyrata</italic> (n = 38), <italic>A. halleri</italic> (n = 30) and <italic>Capsella grandiflora</italic> (n = 7, shown in bold). The phylogeny was obtained by the neighbour-joining method on the proportion of amino-acid differences. Brackets indicate the position of two trans-specifically shared pairs of S-alleles between <italic>A. lyrata</italic> and <italic>A. halleri</italic> that are interrupted by the branching of one S-alleles from <italic>C. grandiflora</italic> (thick lines).</p><p>(1.16 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000168.s002"><label>Figure S2</label><caption><p>Phylogenies of 68 SRK sequences from <italic>A. lyrata</italic> and <italic>A. halleri</italic>. The phylogeny was obtained by the neighbour-joining method on synonymous differences. Bootstrap support was obtained by 1,000 independent replicates.</p><p>(1.30 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000168.s003"><label>Figure S3</label><caption><p>Phylogenies of 68 SRK sequences from <italic>A. lyrata</italic> and <italic>A. halleri</italic>. The phylogeny was obtained by the neighbour-joining method on non-synonymous differences. Bootstrap support was obtained by 1,000 independent replicates.</p><p>(1.36 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000168.s004"><label>Figure S4</label><caption><p>Bootstrap distribution (10,000 replicates) of net divergence for SRK alleles (average across 18 S alleles pairs, in black) and the genomic background (average across 12 control genes, in grey).</p><p>(0.86 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000168.s005"><label>Figure S5</label><caption><p>Distribution of the number of pairwise nucleotide differences for SRK sequences in interspecific comparisons between <italic>A. halleri</italic> and <italic>A. lyrata</italic>, excluding the 18 pairs of sequences considered as transspecific pairs.</p><p>(0.43 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000168.s006"><label>Table S1</label><caption><p>Net divergence estimation for the 12 control genes.</p><p>(0.05 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000168.s007"><label>Text S1</label><caption><p>Supplemental material.</p><p>(0.05 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label/><p>All estimates were Jukes &amp; Cantor corrected. HV refers to hypervariable regions as defined by Nishio and Kusaba ##UREF##16##[60]##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt102"><label/><p>All estimates were Jukes &amp; Cantor corrected.</p></fn><fn id="nt103"><p><bold>*:</bold> Note that Aly 9 is a member of the S-domain gene family, but polymorphism data for this gene was used here to increase the genomic background dataset.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt104"><label/><p>
<italic>N<sub>A = </sub></italic>Effective population size in the common ancestor of <italic>A. lyrata</italic> and <italic>A. halleri.</italic>
</p></fn><fn id="nt105"><label/><p>
<italic>N<sub>lyrata = </sub></italic>Effective population size in <italic>A. lyrata</italic>.</p></fn><fn id="nt106"><label/><p>
<italic>N<sub>halleri = </sub></italic>Effective population size in <italic>A. halleri</italic>.</p></fn><fn id="nt107"><label/><p>
<italic>m<sub>hal→lyr</sub> = </italic>Rate at which genes come into <italic>A. lyrata</italic> from <italic>A. halleri</italic> as time moves forward.</p></fn><fn id="nt108"><label/><p>
<italic>m<sub>lyr→hal</sub></italic> = Rate at which genes come into <italic>A. halleri</italic> from <italic>A. lyrata</italic> as time moves forward.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>Funding from the PPF Bioinformatique from Université des Sciences et Technologies de Lille 1 to VC, from the CNRS-Environment and Sustainable Development Department (ATIP-plus grant) and from the French National Research Agency (ANR-06-BIOD grant) to XV, and from the Danish National Research Council to MHS is gratefully acknowledged.</p></fn></fn-group>
<graphic xlink:href="pgen.1000168.g001"/> <graphic xlink:href="pgen.1000168.g002"/> <graphic id="pgen-1000168-t001-1" xlink:href="pgen.1000168.t001"/> <graphic id="pgen-1000168-t002-2" xlink:href="pgen.1000168.t002"/> <graphic xlink:href="pgen.1000168.g003"/> <graphic xlink:href="pgen.1000168.g004"/> <graphic id="pgen-1000168-t003-3" xlink:href="pgen.1000168.t003"/>
<media xlink:href="pgen.1000168.s001.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000168.s002.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000168.s003.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000168.s004.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000168.s005.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000168.s006.docx"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000168.s007.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["2"], "element-citation": ["\n"], "surname": ["Machado", "Kliman", "Markert", "Hey"], "given-names": ["CA", "RM", "JA", "J"], "year": ["2002"], "article-title": ["Inferring the history of speciation from multilocus DNA sequence data: the case of "], "italic": ["Drosophila pseudoobscura"], "source": ["Mol Biol Evol"], "volume": ["19"], "fpage": ["472488"]}, {"label": ["4"], "element-citation": ["\n"], "surname": ["Arnold"], "given-names": ["ML"], "year": ["2006"], "article-title": ["\n"], "italic": ["Evolution Through Genetic Exchange"]}, {"label": ["9"], "element-citation": ["\n"], "surname": ["Payseur", "Nachman"], "given-names": ["BA", "MW"], "year": ["2005"], "article-title": ["The genomics of speciation: investigating the molecular correlates of X chromosome introgression across the hybrid zone between "], "italic": ["Mus domesticus", "Mus musculus"], "source": ["Biol J Linn Soc Lond"], "volume": ["84"], "fpage": ["523"], "lpage": ["534"]}, {"label": ["12"], "element-citation": ["\n"], "surname": ["Stewart", "Halfhill", "Warwick"], "given-names": ["CN", "MD", "SI"], "year": ["2003"], "article-title": ["Transgene introgression from genetically modified crops to their wild relatives."], "source": ["Nature Reviews Genetics"], "volume": ["4"], "fpage": ["806"], "lpage": ["817"]}, {"label": ["19"], "element-citation": ["\n"], "surname": ["Hoffmann"], "given-names": ["MH"], "year": ["2005"], "article-title": ["Evolution of the realized climatic niche in the genus Arabidopsis (Brassicaceae)."], "source": ["Evolution Int J Org Evolution"], "volume": ["59"], "fpage": ["1425"], "lpage": ["1436"]}, {"label": ["24"], "element-citation": ["\n"], "surname": ["Kimura", "Sato", "Fujimoto", "Nishio"], "given-names": ["R", "K", "R", "T"], "year": ["2002"], "article-title": ["Recognition specificity of self-incompatibility maintained after the divergence of "], "italic": ["Brassica oleraceae", "Brassica rapa."], "source": ["Plant J. 29"], "volume": ["(2)"], "fpage": ["215"], "lpage": ["223"]}, {"label": ["31"], "element-citation": ["\n"], "surname": ["Schierup", "Mable", "Awadalla", "Charlesworth"], "given-names": ["MH", "BK", "P", "D"], "year": ["2001"], "article-title": ["Identification and characterization of a polymorphic receptor kinase gene linked to the self-incompatibility locus of "], "italic": ["Arabidopsis lyrata"], "source": ["Genetics"], "volume": ["157"], "fpage": ["387"], "lpage": ["399"]}, {"label": ["35"], "element-citation": ["\n"], "surname": ["Haldane"], "given-names": ["JBS"], "year": ["1924"], "article-title": ["A mathematical theory of natural and artificial selection, Part I."], "source": ["Trans Camb Philos Soc"], "volume": ["23"], "fpage": ["19"], "lpage": ["41"]}, {"label": ["36"], "element-citation": ["\n"], "surname": ["Turner"], "given-names": ["JRG"], "year": ["1977"], "article-title": ["Butterfly mimicry: the genetical evolution of an adaptation,"], "source": ["Evol Biol"], "volume": ["11"], "fpage": ["163"], "lpage": ["206"]}, {"label": ["39"], "element-citation": ["\n"], "surname": ["Ruggiero", "Jacquemin", "Castric", "Vekemans"], "given-names": ["MV", "B", "V", "X"], "year": ["2008"], "article-title": ["Hitch-hiking to a locus under balancing selection: high sequence diversity and low population subdivision at the S-locus genomic region in "], "italic": ["Arabidopsis halleri"], "source": ["Genet. Res."], "volume": ["90"], "fpage": ["37"], "lpage": ["46"]}, {"label": ["41"], "element-citation": ["\n"], "surname": ["Nei"], "given-names": ["M"], "year": ["1987"], "article-title": ["Molecular evolutionary genetics."], "publisher-loc": ["New York"], "publisher-name": ["Columbia University Press"], "fpage": ["1"], "lpage": ["88"]}, {"label": ["43"], "element-citation": ["\n"], "surname": ["Igic", "Bohs", "Kohn"], "given-names": ["B", "L", "JR"], "year": ["2004"], "article-title": ["Historical inferences from the self-incompatibility locus."], "source": ["New Phytologist"], "volume": ["161"], "fpage": ["97"], "lpage": ["105"]}, {"label": ["47"], "element-citation": ["\n"], "surname": ["Newton", "Allnutt", "Gillies", "Lowe", "Ennos"], "given-names": ["AC", "TR", "ACM", "AJ", "RA"], "year": ["1999"], "article-title": ["Molecular phylogeography, intraspecific variation and the conservation of tree species."], "source": ["Tr Ecol Evol"], "volume": ["14"], "fpage": ["140"], "lpage": ["145"]}, {"label": ["48"], "element-citation": ["\n"], "surname": ["Martinsen", "Whitham", "Turek", "Keim"], "given-names": ["GD", "TG", "RJ", "P"], "year": ["2001"], "article-title": ["Hybrid populations selectively filter gene introgression between species."], "source": ["Evolution Int J Org Evolution"], "volume": ["55"], "fpage": ["1325"], "lpage": ["1335"]}, {"label": ["52"], "element-citation": ["\n"], "surname": ["Hall"], "given-names": ["TA"], "year": ["1999"], "article-title": ["BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT."], "source": ["Nucl Acids Symp Ser"], "volume": ["41"], "fpage": ["95"], "lpage": ["98"]}, {"label": ["55"], "element-citation": ["\n"], "surname": ["Nei", "Kumar"], "given-names": ["M", "S"], "year": ["2000"], "article-title": ["Molecular Evolution and Phylogenetics (Oxford Univ. Press, New York)"]}, {"label": ["60"], "element-citation": ["\n"], "surname": ["Nishio", "Kusaba"], "given-names": ["T", "M"], "year": ["2000"], "article-title": ["Sequence diversity of "], "italic": ["SLG", "SRK", "Brassica oleracea"], "source": ["Ann Bot"], "volume": ["85"], "fpage": ["141"], "lpage": ["146"]}]
PLoS Genet. 2008 Aug 29; 4(8):e1000168