The transition to self-compatibility in Arabidopsis thaliana and evolution within S-haplotypes over 10 Myr

A recent investigation found evidence that the transition of Arabidopsis thaliana from ancestral self-incompatibility (SI) to full self-compatibility occurred very recently and suggested that this occurred through a selective fixation of a nonfunctional allele (PsiSCR1) at the SCR gene, which determines pollen specificity in the incompatibility response. The main evidence is the lack of polymorphism at the SCR locus in A. thaliana. However, the nearby SRK gene, which determines stigma specificity in self-incompatible Brassicaceae species, has extremely high sequence diversity, with 3 very divergent SRK haplotypes, 2 of them present in multiple strains. Such high diversity is extremely unusual in this species, and it suggests the possibility that multiple, different SRK haplotypes may have been preserved from A. thaliana's self-incompatible ancestor. To study the evolution of S-haplotypes in the A. thaliana lineage, we searched the 2 most closely related Arabidopsis species Arabidopsis lyrata and Arabidopsis halleri, in which most populations have retained SI, and found SRK sequences corresponding to all 3 A. thaliana haplogroup sequences. Our molecular evolutionary analyses of these 3 S-haplotypes provide an independent estimate of the timing of the breakdown of SI and again exclude an ancient transition to selfing in A. thaliana. Comparing sequences of each of the 3 haplogroups between species, we find that 2 of the 3 SRK sequences (haplogroups A and B) are similar throughout their length, suggesting that little or no recombination with other SRK alleles has occurred since these species diverged. The diversity difference between the SCR and SRK loci in A. thaliana, however, suggests crossing-over, either within SRK or between the SCR and SRK loci. If the loss of SI involved fixation of the PsiSCR1 sequence, the exchange must have occurred during its fixation. Divergence between the species is much lower at the S-locus, compared with reference loci, and we discuss two contributory possibilities. Introgression may have occurred between A. lyrata and A. halleri and between their ancestral lineage and A. thaliana, at least for some period after their split. In addition, the coalescence times of sequences of individual S-haplogroups are expected to be less than those of alleles at non-S-loci.