Allelic diversity and gene genealogy at the self-incompatibility locus in the Solanaceae

Balancing selection and wild gene pool contribute to resistance in global rice germplasm against planthopper

Interactions and co-evolution between plants and herbivorous insects are critically important in agriculture. Brown planthopper (BPH) is the most severe insect of rice, and the biotypes adapt to feed on different rice genotypes. Here, we present genomics analyses on 1,520 global rice germplasms for resistance to three BPH biotypes. Genome-wide association studies identified 3,502 single nucleotide polymorphisms (SNPs) and 59 loci associated with BPH resistance in rice. We cloned a previously unidentified gene Bph37 that confers resistance to BPH. The associated loci showed high nucleotide diversity. Genome-wide scans for trans-species polymorphisms revealed ancient balancing selection at the loci. The secondarily evolved insect biotypes II and III exhibited significantly higher virulence and overcame more rice varieties than the primary biotype I. In response, more SNPs and loci evolved in rice for resistance to biotypes II and III. Notably, three exceptional large regions with high SNP density and resistance-associated loci on chromosomes 4 and 6 appear distinct between the resistant and susceptible rice varieties. Surprisingly, these regions in resistant rice might have been retained from wild species Oryza nivara. Our findings expand the understanding of long-term interactions between rice and BPH and provide resistance genes and germplasm resources for breeding durable BPH-resistant rice varieties.

Variation in Genetic Mechanisms for Plumage Polymorphism in Skuas (Stercorarius)

Coloration is evolutionarily labile and so provides an excellent trait for examining the repeatability of evolution. Here, we investigate the repeatability of the evolution of polymorphic variation in ventral plumage coloration in skuas (Stercorarius: Stercorariidae). In 2 species, arctic (S. parasiticus) and pomarine skuas (S. pomarinus), plumage polymorphism was previously shown to be associated with coding changes at the melanocortin-1 receptor (MC1R) locus. Here, we show that polymorphism in a third species, the south polar skua (S. maccormicki), is not associated with coding variation at MC1R or with variation at a Z-linked second candidate locus, tyrosinase-related protein 1 (TYRP1). Hence, convergent evolution of plumage polymorphisms in skuas is only partly repeatable at the level of the genetic locus involved. Interestingly, the pattern of repeatability in skuas is aligned not with phylogeny but with the nature of the phenotypic variation. In particular, south polar skuas show a strong sex bias to coloration that is absent in the other species, and it may be that this has a unique genetic architecture.

A molecularphylogeny offorktail damselflies(genus Ischnura)revealsa dynamic macroevolutionary history of female colour polymorphisms

Colour polymorphisms are popular study systems among biologists interested in evolutionary dynamics, genomics, sexual selection and sexual conflict. In many damselfly groups, such as in the globally distributed genus Ischnura (forktails), sex-limited female colour polymorphisms occur in multiple species. Female-polymorphic species contain two or three female morphs, one of which phenotypically matches the male (androchrome or male mimic) and the other(s) which are phenotypically distinct from the male (heterochrome). These female colour polymorphisms are thought to be maintained by frequency-dependent sexual conflict, but their macroevolutionary histories are unknown, due to the lack of a robust molecular phylogeny. Here, we present the first time-calibrated phylogeny of Ischnura, using a multispecies coalescent approach (StarBEAST2) and incorporating both molecular and fossil data for 41 extant species (55% of the genus). We estimate the age of Ischnura to be between 13.8 and 23.4 millions of years, i.e. Miocene. We infer the ancestral state of this genus as female monomorphism with heterochrome females, with multiple gains and losses of female polymorphisms, evidence of trans-species female polymorphisms and a significant positive relationship between female polymorphism incidence and current geographic range size. Our study provides a robust phylogenetic framework for future research on the dynamic macroevolutionary history of this clade with its extraordinary diversity of sex-limited female polymorphisms.

Population Genomics of the "Arcanum" Species Group in Wild Tomatoes: Evidence for Separate Origins of Two Self-Compatible Lineages

Given their diverse mating systems and recent divergence, wild tomatoes (Solanum section Lycopersicon) have become an attractive model system to study ecological divergence, the build-up of reproductive barriers, and the causes and consequences of the breakdown of self-incompatibility. Here we report on a lesser-studied group of species known as the "Arcanum" group, comprising the nominal species Solanum arcanum, Solanum chmielewskii, and Solanum neorickii. The latter two taxa are self-compatible but are thought to self-fertilize at different rates, given their distinct manifestations of the morphological "selfing syndrome." Based on experimental crossings and transcriptome sequencing of a total of 39 different genotypes from as many accessions representing each species' geographic range, we provide compelling evidence for deep genealogical divisions within S. arcanum; only the self-incompatible lineage known as "var. marañón" has close genealogical ties to the two self-compatible species. Moreover, there is evidence under multiple inference schemes for different geographic subsets of S. arcanum var. marañón being closest to S. chmielewskii and S. neorickii, respectively. To broadly characterize the population-genomic consequences of these recent mating-system transitions and their associated speciation events, we fit demographic models indicating strong reductions in effective population size, congruent with reduced nucleotide and S-locus diversity in the two independently derived self-compatible species.

Inferences on specificity recognition at the Malus×domestica gametophytic self-incompatibility system

In Malus × domestica (Rosaceae) the product of each SFBB gene (the pollen component of the gametophytic self-incompatibility (GSI) system) of a S-haplotype (the combination of pistil and pollen genes that are linked) interacts with a sub-set of non-self S-RNases (the pistil component), but not with the self S-RNase. To understand how the Malus GSI system works, we identified 24 SFBB genes expressed in anthers, and determined their gene sequence in nine M. domestica cultivars. Expression of these SFBBs was not detected in the petal, sepal, filament, receptacle, style, stigma, ovary or young leaf. For all SFBBs (except SFBB15), identical sequences were obtained only in cultivars having the same S-RNase. Linkage with a particular S-RNase was further established using the progeny of three crosses. Such data is needed to understand how other genes not involved in GSI are affected by the S-locus region. To classify SFBBs specificity, the amino acids under positive selection obtained when performing intra-haplotypic analyses were used. Using this information and the previously identified S-RNase positively selected amino acid sites, inferences are made on the S-RNase amino acid properties (hydrophobicity, aromatic, aliphatic, polarity, and size), at these positions, that are critical features for GSI specificity determination.

Allelic diversity of S-RNase alleles in diploid potato species

The S-ribonuclease sequences of 16 S-alleles derived from diploid types of Solanum are presented. A phylogenetic analysis and partial phenotypic analysis support the conclusion that these are functional S-alleles. S-Ribonucleases (S-RNases) control the pistil specificity of the self-incompatibility (SI) response in the genus Solanum and several other members of the Solanaceae. The nucleotide sequences of S-RNases corresponding to a large number of S-alleles or S-haplotypes have been characterised. However, surprisingly, few S-RNase sequences are available for potato species. The identification of new S-alleles in diploid potato species is desirable as these stocks are important sources of traits such as biotic and abiotic resistance. S-RNase sequences are reported here from three distinct diploid types of potato: cultivated Solanum tuberosum Group Phureja, S. tuberosum Group Stenotomum, and the wild species Solanum okadae. Partial S-RNase sequences were obtained from pistil RNA by RT-PCR or 3'RACE (Rapid Amplification of cDNA Ends) using a degenerate primer. Full-length sequences were obtained for two alleles by 5'RACE. Database searches with these sequences identified 16 S-RNases in total, all of which are novel. The sequence analysis revealed all the expected features of functional S-RNases. Phylogenetic analysis with selected published S-RNase and S-like-RNase sequences from the Solanaceae revealed extensive trans-generic evolution of the S-RNases and a clear distinction from S-like-RNases. Pollination tests were used to confirm the self-incompatibility status and cross-compatibility relationships of the S. okadae accessions. All the S. okadae accessions were found to be self-incompatible as expected with crosses amongst them exhibiting both cross-compatibility and semi-compatibility consistent with the S-genotypes determined from the S-RNase sequence data. The progeny analysis of four semi-compatible crosses examined by allele-specific PCR provided further confirmation that these are functional S-RNases.

Coalescence with Background and Balancing Selection in Systems with Bi- and Uniparental Reproduction: Contrasting Partial Asexuality and Selfing

Uniparental reproduction in diploids, via asexual reproduction or selfing, reduces the independence with which separate loci are transmitted across generations. This is expected to increase the extent to which a neutral marker is affected by selection elsewhere in the genome. Such effects have previously been quantified in coalescent models involving selfing. Here we examine the effects of background selection and balancing selection in diploids capable of both sexual and asexual reproduction (i.e., partial asexuality). We find that the effect of background selection on reducing coalescent time (and effective population size) can be orders of magnitude greater when rates of sex are low than when sex is common. This is because asexuality enhances the effects of background selection through both a recombination effect and a segregation effect. We show that there are several reasons that the strength of background selection differs between systems with partial asexuality and those with comparable levels of uniparental reproduction via selfing. Expectations for reductions in Ne via background selection have been verified using stochastic simulations. In contrast to background selection, balancing selection increases the coalescence time for a linked neutral site. With partial asexuality, the effect of balancing selection is somewhat dependent upon the mode of selection (e.g., heterozygote advantage vs. negative frequency-dependent selection) in a manner that does not apply to selfing. This is because the frequency of heterozygotes, which are required for recombination onto alternative genetic backgrounds, is more dependent on the pattern of selection with partial asexuality than with selfing.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

Fixation properties of subdivided populations with balancing selection

In subdivided populations, migration acts together with selection and genetic drift and determines their evolution. Building upon a recently proposed method, which hinges on the emergence of a time scale separation between local and global dynamics, we study the fixation properties of subdivided populations in the presence of balancing selection. The approximation implied by the method is accurate when the effective selection strength is small and the number of subpopulations is large. In particular, it predicts a phase transition between species coexistence and biodiversity loss in the infinite-size limit and, in finite populations, a nonmonotonic dependence of the mean fixation time on the migration rate. In order to investigate the fixation properties of the subdivided population for stronger selection, we introduce an effective coarser description of the dynamics in terms of a voter model with intermediate states, which highlights the basic mechanisms driving the evolutionary process.

The founding of Mauritian endemic coffee trees by a synchronous long-distance dispersal event

The stochastic process of long-distance dispersal is the exclusive means by which plants colonize oceanic islands. Baker's rule posits that self-incompatible plant lineages are unlikely to successfully colonize oceanic islands because they must achieve a coordinated long-distance dispersal of sufficiently numerous individuals to establish an outcrossing founder population. Here, we show for the first time that Mauritian Coffea species are self-incompatible and thus represent an exception to Baker's rule. The genus Coffea (Rubiaceae) is composed of approximately 124 species with a paleotropical distribution. Phylogenetic evidence strongly supports a single colonization of the oceanic island of Mauritius from either Madagascar or Africa. We employ Bayesian divergence time analyses to show that the colonization of Mauritius was not a recent event. We genotype S-RNase alleles from Mauritian endemic Coffea, and using S-allele gene genealogies, we show that the Mauritian allelic diversity is confined to just seven deeply divergent Coffea S-RNase allelic lineages. Based on these data, we developed an individual-based model and performed a simulation study to estimate the most likely number of founding individuals involved in the colonization of Mauritius. Our simulations show that to explain the observed S-RNase allelic diversity, the founding population was likely composed of fewer than 31 seeds that were likely synchronously dispersed from an ancestral mainland species.

Nonmonotonic effects of migration in subdivided populations

The influence of migration on the stochastic dynamics of subdivided populations is still an open issue in various evolutionary models. Here, we develop a self-consistent mean-field-like method in order to determine the effects of migration on relevant nonequilibrium properties, such as the mean fixation time. If evolution strongly favors coexistence of species (e.g., balancing selection), the mean fixation time develops an unexpected minimum as a function of the migration rate. Our analysis hinges only on the presence of a separation of time scales between local and global dynamics, and therefore, it carries over to other nonequilibrium processes in physics, biology, ecology, and social sciences.

Extensive trans-specific polymorphism at the mating type locus of the root decay fungus Heterobasidion

Incompatibility systems in which individuals bearing identical alleles reject each other favor the maintenance of a diversity of alleles. Mushroom mating type loci (MAT) encode for dozens or hundreds of incompatibility alleles whose loss from the population is greatly restricted through negative frequency selection, leading to a system of alleles with highly divergent sequences. Here, we use DNA sequences of homeodomain (HD) encoding genes at the MAT locus of five closely related species of the root rot basidiomycete Heterobasidion annosum sensu lato to show that the extended coalescence time of MAT alleles greatly predates speciation in the group, contrasting loci outside of MAT that show allele divergences largely consistent with the species phylogeny with those of MAT, which show rampant trans-species polymorphism. We observe a roughly 6-fold greater genealogical depth and polymorphism of MAT compared with non-MAT that argues for the maintenance of balanced polymorphism for a minimum duration of 24 My based on a molecular-clock calibrated species phylogeny. As with other basidiomycete HD genes, balancing selection appears to be concentrated at the specificity-determining region in the N-terminus of the protein based on identification of codons under selection and the absence of recombination within the region. However, the elevated polymorphism extends into the nonspecificity determining regions as well as a neighboring non-MAT gene, the mitochondrial intermediate peptidase (MIP). In doing so, increased divergence should decrease recombination among alleles and as a by-product create incompatibilities in the functional domains not involved in allele recognition but in regulating sexual development.

Learning from rejection: the evolutionary biology of single-locus incompatibility

The self-incompatibility (S-) locus of flowering plants is among the most polymorphic known. PCR methods can now be used to estimate both the number of alleles in natural populations and their sequence diversity. The number of alleles provides an estimate of recent effective population size, thus the S-locus provides a tool for examining how species characteristics affect population size. Sequence relationships among alleles provide another estimate of population size extending millions of years into the past. Relationships between S-alleles and related genes provide a means of dating the age of origin of incompatibility systems and determining which, if any, angiosperm families share incompatibility by homology.

Testing for "snowballing" hybrid incompatibilities in Solanum: impact of ancestral polymorphism and divergence estimates

Two recent high-profile studies offered empirical evidence for a "snowballing" accumulation of postzygotic incompatibilities in Drosophila and Solanum (tomatoes). Here we present a reanalysis of the Solanum data that is motivated by population genetic principles. Specifically, the high levels of intraspecific nucleotide polymorphism in wild tomato species and presumably large effective population size throughout the divergence history of this clade imply that ancestral polymorphism should be taken into account when evaluating sequence divergence between species. Based on our reanalyses of synonymous-site divergence between the four focal Solanum species and a wide range of ancestral polymorphism, we assessed under which conditions the reported accumulation of seed sterility factors supports the snowball effect. Our results highlight the pivotal impact of levels of ancestral polymorphism and alternate divergence values, and they illustrate that robust tests of the snowball effect in Solanum require genome-wide estimates of divergence.

Evolutionary genetics of an S-like polymorphism in Papaveraceae with putative function in self-incompatibility

BACKGROUND

Papaver rhoeas possesses a gametophytic self-incompatibility (SI) system not homologous to any other SI mechanism characterized at the molecular level. Four previously published full length stigmatic S-alleles from the genus Papaver exhibited remarkable sequence divergence, but these studies failed to amplify additional S-alleles despite crossing evidence for more than 60 S-alleles in Papaver rhoeas alone.

METHODOLOGY/PRINCIPAL FINDINGS

Using RT-PCR we identified 87 unique putative stigmatic S-allele sequences from the Papaveraceae Argemone munita, Papaver mcconnellii, P. nudicuale, Platystemon californicus and Romneya coulteri. Hand pollinations among two full-sib families of both A. munita and P. californicus indicate a strong correlation between the putative S-genotype and observed incompatibility phenotype. However, we also found more than two S-like sequences in some individuals of A. munita and P. californicus, with two products co-segregating in both full-sib families of P. californicus. Pairwise sequence divergence estimates within and among taxa show Papaver stigmatic S-alleles to be the most variable with lower divergence among putative S-alleles from other Papaveraceae. Genealogical analysis indicates little shared ancestral polymorphism among S-like sequences from different genera. Lack of shared ancestral polymorphism could be due to long divergence times among genera studied, reduced levels of balancing selection if some or all S-like sequences do not function in incompatibility, population bottlenecks, or different levels of recombination among taxa. Preliminary estimates of positive selection find many sites under selective constraint with a few undergoing positive selection, suggesting that self-recognition may depend on amino acid substitutions at only a few sites.

CONCLUSIONS/SIGNIFICANCE

Because of the strong correlation between genotype and SI phenotype, sequences reported here represent either functional stylar S-alleles, tightly linked paralogs of the S-locus or a combination of both. The considerable complexity revealed in this study shows we have much to learn about the evolutionary dynamics of self-incompatibility systems.

Differential strengths of selection on S-RNases from Physalis and Solanum (Solanaceae)

Background

The S-RNases of the Solanaceae are highly polymorphic self-incompatibility (S-) alleles subject to strong balancing selection. Relatively recent diversification of S-alleles has occurred in the genus Physalis following a historical restriction of S-allele diversity. In contrast, the genus Solanum did not undergo a restriction of S-locus diversity and its S-alleles are generally much older. Because recovery from reduced S-locus diversity should involve increased selection, we employ a statistical framework to ask whether S-locus selection intensities are higher in Physalis than Solanum. Because different S-RNase lineages diversify in Physalis and Solanum, we also ask whether different sites are under selection in different lineages.

Results

Maximum-likelihood and Bayesian coalescent methods found higher intensities of selection and more sites under significant positive selection in the 48 Physalis S-RNase alleles than the 49 from Solanum. Highest posterior densities of dN/dS (ω) estimates show that the strength of selection is greater for Physalis at 36 codons. A nested maximum likelihood method was more conservative, but still found 16 sites with greater selection in Physalis. Neither method found any codons under significantly greater selection in Solanum. A random effects likelihood method that examines data from both taxa jointly confirmed higher selection intensities in Physalis, but did not find different proportions of sites under selection in the two datasets. The greatest differences in strengths of selection were found in the most variable regions of the S-RNases, as expected if these regions encode self-recognition specificities. Clade-specific likelihood models indicated some codons were under greater selection in background Solanum lineages than in specific lineages of Physalis implying that selection on sites may differ among lineages.

Conclusions

Likelihood and Bayesian methods provide a statistical approach to testing differential selection across populations or species. These tests appear robust to the levels of polymorphism found in diverse S-allele collections subject to strong balancing selection. As predicted, the intensity of selection at the S-locus was higher in the taxon with more recent S-locus diversification. This is the first confirmation by statistical test of differing selection intensities among self-incompatibility alleles from different populations or species.

Darwin and genetics

Darwin's theory of natural selection lacked an adequate account of inheritance, making it logically incomplete. We review the interaction between evolution and genetics, showing how, unlike Mendel, Darwin's lack of a model of the mechanism of inheritance left him unable to interpret his own data that showed Mendelian ratios, even though he shared with Mendel a more mathematical and probabilistic outlook than most biologists of his time. Darwin's own "pangenesis" model provided a mechanism for generating ample variability on which selection could act. It involved, however, the inheritance of characters acquired during an organism's life, which Darwin himself knew could not explain some evolutionary situations. Once the particulate basis of genetics was understood, it was seen to allow variation to be passed intact to new generations, and evolution could then be understood as a process of changes in the frequencies of stable variants. Evolutionary genetics subsequently developed as a central part of biology. Darwinian principles now play a greater role in biology than ever before, which we illustrate with some examples of studies of natural selection that use DNA sequence data and with some recent advances in answering questions first asked by Darwin.

Mating-system evolution: genies from a bottleneck

Evolutionary shifts from outcrossing to selfing have been frequent in plants, but little is known about how this occurs. Two new studies of the same species point to a recent bottleneck that coincided with one such shift.

Fundamental concepts in genetics: effective population size and patterns of molecular evolution and variation

The effective size of a population, N(e), determines the rate of change in the composition of a population caused by genetic drift, which is the random sampling of genetic variants in a finite population. N(e) is crucial in determining the level of variability in a population, and the effectiveness of selection relative to drift. This article reviews the properties of N(e) in a variety of different situations of biological interest, and the factors that influence it. In particular, the action of selection means that N(e) varies across the genome, and advances in genomic techniques are giving new insights into how selection shapes N(e).

Speciation and the evolution of gamete recognition genes: pattern and process

Proteins on gamete surfaces are major determinants of fertilization success, particularly in free-spawning animals. Molecular analyses of these simple genetic systems show rapid evolution, positive selection, accelerated coalescence and, sometimes, extensive polymorphism. Careful analysis of the behavior of sperm produced by males with different gamete alleles shows that these alleles can deliver significant functional differences. Three forms of allele-specific fertilization advantage have been shown: assortative mating based on gamete type, rare allele advantage and heterozygote superiority. Models suggest that sperm and egg proteins may be coevolutionary partners that can alternate between directional selection for high fertilization ability and cyclic adaptation of eggs and sperm driven by sexual conflict. These processes act within allopatric populations and may accelerate their divergence if gamete adaptations in separate demes reduce cross-fertilization. Reproductive character displacement by reinforcement may play a diversifying role when previously allopatric populations rejoin. In circumstance that might prove to be common, divergence in sympatry can be driven by sexual conflict or by association of mating types with ecological differences. The ecology of fertilization, especially the degree of sperm competition and egg death via polyspermy, are important determinants of the strength and direction of selection on gametes. Free-spawning animals allow careful analysis of gamete recognition -from the behavior of adults and interactions of gametes, to molecular patterns of allele divergence. Future research efforts on the evolutionary consequences of fertilization ecology, and the interaction between extensive variation in egg surface proteins and sperm fertilization ability, are particularly needed.

Evolutionary genetics: changed sex determination in honeybees

In several insects and fish, and probably some mammals, the gene controlling the male-female switch has changed during evolution. It now seems that this has also happened in honeybees, where the sex-determining gene has now been shown to be a duplicate of another Hymenopteran sex-determining gene.

Ecological and evolutionary genomics in the wild tomatoes (Solanum sect. Lycopersicon)

The plant group Solanum section Lycopersicon (the clade containing the domesticated tomato and its wild relatives) is ideal for integrating genomic tools and approaches into ecological and evolutionary research. Wild species within Lycopersicon span broad morphological, physiological, life history, mating system, and biochemical variation, and are separated by substantial, but incomplete postmating reproductive barriers, making this an ideal system for genetic analyses of these traits. This ecological and evolutionary diversity is matched by many logistical advantages, including extensive historical occurrence records for all species in the group, publicly available germplasm for hundreds of known wild accessions, demonstrated experimental tractability, and extensive genetic, genomic, and functional tools and information from the tomato research community. Here I introduce the numerous advantages of this system for Ecological and Evolutionary Functional Genomics (EEFG), and outline several ecological and evolutionary phenotypes and questions that can be fruitfully tackled in this system. These include biotic and abiotic adaptation, reproductive trait evolution, and the genetic basis of speciation. With the modest enhancement of some research strengths, this system is poised to join the best of our currently available model EEFG systems.

Molecular characterization of Lal2, an SRK-like gene linked to the S-locus in the wild mustard Leavenworthia alabamica

Single-locus sporophytic self-incompatibility inhibits inbreeding in many members of the mustard family (Brassicaceae). To investigate the genetics of self-incompatibility in the wild mustard Leavenworthia alabamica, diallel crosses were conducted between full siblings. Patterns of incompatibility were consistent with the action of single-locus sporophytic self-incompatibility. DNA sequences related to S-locus receptor kinase (SRK), the gene involved in self-pollen recognition in mustards, were cloned and sequenced. A single sequence with high identity to SRK and several other groups of sequences (Lal1, Lal2, Lal3, Lal8, and Lal14) were isolated from L. alabamica. We propose that either Lal2 sequences are divergent alleles of SRK or Lal2 is in tight linkage with SRK because (1) Lal2 alleles cosegregate with S-alleles inferred from dialleles in all 97 cases tested in five families; (2) Lal2 sequences are highly diverse at both synonymous and nonsynonymous sites and exhibit patterns of selective constraint similar to those observed at SRK in Brassica and Arabidopsis; and (3) transcripts of one Lal2 allele were detected in leaves and the styles of open flowers, but were most abundant in the stigmas of maturing buds. We discuss the utility of the S-linked polymorphism at Lal2 for studying the evolutionary forces acting on self-incompatibility in Leavenworthia.

The number, age, sharing and relatedness of S-locus specificities in Prunus

In gametophytic self-incompatibility systems, many specificities (different 'lock-and-key' combinations) are maintained by frequency-dependent selection for very long evolutionary times. In Solanaceae, trans-specific evolution (the observation that an allele from one species may be more closely related to an allele from another species than to others from the same species) has been taken as an argument for the very old age of specificities. In this work, by determining, for the first time, the age of extant Prunus species, we show that this reasoning cannot be applied to Prunoideae. Furthermore, since our sample size is large (all S-RNase encoding the female component and SFB encoding the male component GenBank sequences), we were able to estimate the age of the oldest Prunus specificities. By doing so, we show that the lower variability levels at the Prunus S-locus, in comparison with Solanaceae, is due to the younger age of Prunus alleles, and not to a difference in silent mutation rates. We show that the ancestor to extant Prunus species harboured at least 102 specificities, in contrast to the maximum of 33 observed in extant Prunus species. Since the number of specificities that can be maintained in a population depends on the effective population size, this observation suggests a bottleneck in Prunus evolutionary history. Loss of specificities may have occurred during this event. Using only information on amino acid sites that determine specificity differences, and a simulation approach, we show that a model that assumes closely related specificities are not preferentially lost during evolution, fails to predict the observed degree of specificity relatedness.

The evolution of self-incompatibility when mates are limiting

Self-incompatibility (SI) is a genetic barrier to inbreeding that is broadly distributed in angiosperms. In finite populations of SI plants, the loss of S-allele diversity can limit plant reproduction by reducing the availability of compatible mates. Many studies have shown that small or fragmented plant populations suffer from mate limitation. The advent of molecular typing of S-alleles in many species has paved the way to address quantitatively the importance of mate limitation, and to provide greater insight into why and how SI systems breakdown frequently in nature. In this review, we highlight the ecological factors that contribute to mate limitation in SI taxa, discuss their consequences for the evolution and functioning of SI, and propose new empirical research directions.

A tale of two continents: Baker's rule and the maintenance of self-incompatibility in Lycium (Solanaceae)

Over 50 years ago, Baker (1955, 1967) suggested that self-compatible species were more likely than self-incompatible species to establish new populations on oceanic islands. His logic was straightforward and rested on the assumption that colonization was infrequent; thus, mate limitation favored the establishment of self-fertilizing individuals. In support of Baker's rule, many authors have documented high frequencies of self-compatibility on islands, and recent work has solidified the generality of Baker's ideas. The genus Lycium (Solanaceae) has ca. 80 species distributed worldwide, and phylogenetic studies suggest that Lycium originated in South America and dispersed to the Old World a single time. Previous analyses of the S-RNase gene, which controls the stylar component of self-incompatibility, have shown that gametophytically controlled self-incompatibility is ancestral within the genus, making Lycium a good model for investigating Baker's assertions concerning reproductive assurance following oceanic dispersal. Lycium is also useful for investigations of reproductive evolution, given that species vary both in sexual expression and the presence of self-incompatibility. A model for the evolution of gender dimorphism suggests that polyploidy breaks down self-incompatibility, leading to the evolution of gender dimorphism, which arises as an alternative outcrossing mechanism. There is a perfect association of dimorphic gender expression, polyploidy, and self-compatibility (vs. cosexuality, diploidy, and self-incompatibility) among North American Lycium. Although the association between ploidy level and gender expression also holds for African Lycium, to date no studies of mating systems have been initiated in Old World species. Here, using controlled pollinations, we document strong self-incompatibility in two cosexual, diploid species of African Lycium. Further, we sequence the S-RNase gene in 15 individuals from five cosexual, diploid species of African Lycium and recover 24 putative alleles. Genealogical analyses indicate reduced trans-generic diversity of S-RNases in the Old World compared to the New World. We suggest that genetic diversity at this locus was reduced as a result of a founder event, but, despite the bottleneck, self-incompatibility was maintained in the Old World. Maximum-likelihood analyses of codon substitution patterns indicate that positive Darwinian selection has been relatively strong in the Old World, suggesting the rediversification of S-RNases following a bottleneck. The present data thus provide a dramatic exception to Baker's rule, in addition to supporting a key assumption of the Miller and Venable (2000) model, namely that self-incompatibility is associated with diploidy and cosexuality.

Genetic and molecular characterization of three novel S-haplotypes in sour cherry (Prunus cerasus L.)

Tetraploid sour cherry (Prunus cerasus L.) exhibits gametophytic self-incompatibility (GSI) whereby the specificity of self-pollen rejection is controlled by alleles of the stylar and pollen specificity genes, S-RNase and SFB (S haplotype-specific F-box protein gene), respectively. As sour cherry selections can be either self-compatible (SC) or self-incompatible (SI), polyploidy per se does not result in SC. Instead the genotype-dependent loss of SI in sour cherry is due to the accumulation of non-functional S-haplotypes. The presence of two or more non-functional S-haplotypes within sour cherry 2x pollen renders that pollen SC. Two new S-haplotypes from sour cherry, S(33) and S(34), that are presumed to be contributed by the P. fruticosa species parent, the complete S-RNase and SFB sequences of a third S-haplotype, S(35), plus the presence of two previously identified sweet cherry S-haplotypes, S(14) and S(16) are described here. Genetic segregation data demonstrated that the S(16)-, S(33)-, S(34)-, and S(35)-haplotypes present in sour cherry are fully functional. This result is consistent with our previous finding that 'hetero-allelic' pollen is incompatible in sour cherry. Phylogenetic analyses of the SFB and S-RNase sequences from available Prunus species reveal that the relationships among S-haplotypes show no correspondence to known organismal relationships at any taxonomic level within Prunus, indicating that polymorphisms at the S-locus have been maintained throughout the evolution of the genus. Furthermore, the phylogenetic relationships among SFB sequences are generally incongruent with those among S-RNase sequences for the same S-haplotypes. Hypotheses compatible with these results are discussed.

Expression of 10 S-class SLF-like genes in Nicotiana alata pollen and its implications for understanding the pollen factor of the S locus

The S locus of Nicotiana alata encodes a polymorphic series of ribonucleases (S-RNases) that determine the self-incompatibility (SI) phenotype of the style. The pollen product of the S locus (pollen S) in N. alata is unknown, but in species from the related genus Petunia and in self-incompatible members of the Plantaginaceae and Rosaceae, this function has been assigned to an F-box protein known as SLF or SFB. Here we describe the identification of 10 genes (designated DD1-10) encoding SLF-related proteins that are expressed in N. alata pollen. Because our approach to cloning the DD genes was based on sequences of SLFs from other species, we presume that one of the DD genes encodes the N. alata SLF ortholog. Seven of the DD genes were exclusively expressed in pollen and a low level of sequence variation was found in alleles of each DD gene. Mapping studies confirmed that all 10 DD genes were linked to the S locus and that at least three were located in the same chromosomal segment as pollen S. Finally, the different topologies of the phylogenetic trees produced using available SLF-related sequences and those produced using S-RNase sequences suggests that pollen S and the S-RNase have different evolutionary histories.

Studies of self-incompatibility in wild tomatoes: I. S-allele diversity in Solanum chilense (Dun.) Reiche [corrected] (Solanaceae)

We characterized the molecular allelic variation of RNases at the self-incompatibility (SI) locus of Solanum chilense Dun. We recovered 30 S-RNase allele sequences from 34 plants representing a broad geographic sample. This yielded a species-wide estimate of 35 (95% likelihood interval 31-40) S-alleles. We performed crosses to confirm the association with SI function of 10 of the putative S-RNase allele sequences. Results in all cases were consistent with the expectation that these sequences represent functional alleles under single-locus gametophytic SI. We used the allele sequences to conduct an analysis of selection, as measured by the excess of nonsynonymous changes per site, and found evidence for adaptive changes both within the traditionally defined hypervariable regions and downstream, near the 3'-end of the molecule.

Evolution of the complementary sex-determination gene of honey bees: balancing selection and trans-species polymorphisms

The mechanism of sex determination varies substantively among evolutionary lineages. One important mode of genetic sex determination is haplodiploidy, which is used by approximately 20% of all animal species, including >200,000 species of the entire insect order Hymenoptera. In the honey bee Apis mellifera, a hymenopteran model organism, females are heterozygous at the csd (complementary sex determination) locus, whereas males are hemizygous (from unfertilized eggs). Fertilized homozygotes develop into sterile males that are eaten before maturity. Because homozygotes have zero fitness and because common alleles are more likely than rare ones to form homozygotes, csd should be subject to strong overdominant selection and negative frequency-dependent selection. Under these selective forces, together known as balancing selection, csd is expected to exhibit a high degree of intraspecific polymorphism, with long-lived alleles that may be even older than the species. Here we sequence the csd genes as well as randomly selected neutral genomic regions from individuals of three closely related species, A. mellifera, Apis cerana, and Apis dorsata. The polymorphic level is approximately seven times higher in csd than in the neutral regions. Gene genealogies reveal trans-species polymorphisms at csd but not at any neutral regions. Consistent with the prediction of rare-allele advantage, nonsynonymous mutations are found to be positively selected in csd only in early stages after their appearances. Surprisingly, three different hypervariable repetitive regions in csd are present in the three species, suggesting variable mechanisms underlying allelic specificities. Our results provide a definitive demonstration of balancing selection acting at the honey bee csd gene, offer insights into the molecular determinants of csd allelic specificities, and help avoid homozygosity in bee breeding.

Fine scale spatial structuring of sex and mitochondria in Silene vulgaris

Fine scale spatial structure (FSSS) of cytoplasmic genes in plants is thought to be generated via founder events and can be amplified when seeds germinate close to their mother. In gynodioecious species these processes are expected to generate FSSS in sex ratio because maternally inherited cytoplasmic male sterility genes partially influence sex expression. Here we document a striking example of FSSS in both mitochondrial genetic markers and sex in roadside populations of Silene vulgaris. We show that in one population FSSS of sexes influences relative fruit production of females compared to hermaphrodites. Furthermore, FSSS in sex ratio is expected to persist into future generations because offspring sex ratios from females are female-biased whereas offspring sex ratios from hermaphrodites are hermaphrodite-biased. Earlier studies indicated that pollen limitation is the most likely mechanism underlying negative frequency dependent fitness of females. Our results support the theoretical predictions that FSSS in sex ratio can reduce female fitness by decreasing the frequency at which females experience hermaphrodites. We argue that the influence of FSSS on female fitness is complementary to the influence of larger scale population structure on female fitness, and that population structure at both scales will act to decrease female frequencies in gynodioecious species. Better comprehension of the spatial structure of genders and genes controlling sex expression at a local scale is required for future progress toward understanding sex ratio evolution in gynodioecious plants.

Comparative analysis of the within-population genetic structure in wild cherry (Prunus avium L.) at the self-incompatibility locus and nuclear microsatellites

Gametophytic self-incompatibility (SI) systems in plants exhibit high polymorphism at the SI controlling S-locus because individuals with rare alleles have a higher probability to successfully pollinate other plants than individuals with more frequent alleles. This process, referred to as frequency-dependent selection, is expected to shape number, frequency distribution, and spatial distribution of self-incompatibility alleles in natural populations. We investigated the genetic diversity and the spatial genetic structure within a Prunus avium population at two contrasting gene loci: nuclear microsatellites and the S-locus. The S-locus revealed a higher diversity (15 alleles) than the eight microsatellites (4-12 alleles). Although the frequency distribution of S-alleles differed significantly from the expected equal distribution, the S-locus showed a higher evenness than the microsatellites (Shannon's evenness index for the S-locus: E = 0.91; for the microsatellites: E = 0.48-0.83). Also, highly significant deviations from neutrality were found for the S-locus whereas only minor deviations were found for two of eight microsatellites. A comparison of the frequency distribution of S-alleles in three age-cohorts revealed no significant differences, suggesting that different levels of selection acting on the S-locus or on S-linked sites might also affect the distribution and dynamics of S-alleles. Autocorrelation analysis revealed a weak but significant spatial genetic structure for the multilocus average of the microsatellites and for the S-locus, but could not ascertain differences in the extent of spatial genetic structure between these locus types. An indirect estimate of gene dispersal, which was obtained to explain this spatial genetic pattern, indicated high levels of gene dispersal within our population (sigma(g) = 106 m). This high gene dispersal, which may be partly due to the self-incompatibility system itself, aids the effective gene flow of the microsatellites, thereby decreasing the contrast between the neutral microsatellites and the S-locus.

Balancing selection and its effects on sequences in nearby genome regions

Our understanding of balancing selection is currently becoming greatly clarified by new sequence data being gathered from genes in which polymorphisms are known to be maintained by selection. The data can be interpreted in conjunction with results from population genetics models that include recombination between selected sites and nearby neutral marker variants. This understanding is making possible tests for balancing selection using molecular evolutionary approaches. Such tests do not necessarily require knowledge of the functional types of the different alleles at a locus, but such information, as well as information about the geographic distribution of alleles and markers near the genes, can potentially help towards understanding what form of balancing selection is acting, and how long alleles have been maintained.

Trans-specificity at loci near the self-incompatibility loci in Arabidopsis

We compared allele sequences of two loci near the Arabidopsis lyrata self-incompatibility (S) loci with sequences of A. thaliana orthologs and found high numbers of shared polymorphisms, even excluding singletons and sites likely to be highly mutable. This suggests maintenance of entire S-haplotypes for long evolutionary times and extreme recombination suppression in the region.

Ancient polymorphism reveals unidirectional breeding system shifts

Loss of complex characters is thought to be irreversible (Dollo's law). However, hypotheses of irreversible evolution are remarkably difficult to test, especially when character transitions are frequent. In such cases, inference of ancestral states, in the absence of fossil evidence, is uncertain and represents the single greatest constraint for reconstructing the evolutionary history of characters. Breeding system character transitions are of particular interest because they affect the amount and distribution of genetic variation within species. Transitions from obligate outcrossing to partial or predominant self-fertilization are thought to represent one of the most common trends in flowering plants. We use the unique molecular genetic properties (manifested as deep persistent polymorphisms) of the locus that enforces outcrossing to demonstrate that its loss is irreversible in the plant family Solanaceae. We argue that current phylogenetic methods of reconstruction are potentially inadequate in cases where ancestral state information is inferred by using only the phylogeny and the distribution of character states in extant taxa. This study shows in a statistical framework that a particular character transition is irreversible, consistent with Dollo's law.

Plant self-incompatibility systems: a molecular evolutionary perspective

Incompatibility recognition systems preventing self-fertilization have evolved several times in independent lineages of Angiosperm plants, and three main model systems are well characterized at the molecular level [the gametophytic self-incompatibility (SI) systems of Solanaceae, Rosaceae and Anthirrhinum, the very different system of poppy, and the system in Brassicaceae with sporophytic control of pollen SI reactions]. In two of these systems, the genes encoding both components of pollen-pistil recognition are now known, showing clearly that these two proteins are distinct, that is, SI is a lock-and-key mechanism. Here, we review recent findings in the three well-studied systems in the light of these results and analyse their implications for understanding polymorphism and coevolution of the two SI genes, in the context of a tightly linked genome region.

Historical events and allelic polymorphism at the gametophytic self-incompatibility locus in Solanaceae

The historical migration rate of a species is often difficult to estimate with neutral markers, because the relationship between the turnover time of the markers and the age of the species commonly remains unknown. Compared with neutral markers, the plant self-incompatibility locus (S) provides a much better source of data for migration-rate estimation due to its high allelic polymorphism and antiquity. Here, the results from extensive surveys of S alleles in two wild solanaceous species, Solanum carolinense and Physalis longifolia, indicate that historical migration rates have differed significantly between the species; the higher migration rate found in S. carolinense appears to have interacted with the balancing selection at the S locus to result in fewer S alleles being maintained in the species. Historical population growth rates estimated via a modified coalescent approach also suggest a faster growing population for S. carolinense than for P. longifolia, which would have further widened their interspecific difference in S-allelle polymorphism. These historical factors may have reduced the probability of new S alleles to prevailing in S. carolinense, leaving old ones segregating at the S locus with little signature of positive selection being currently detectable.

Genetic diversity and the evolutionary history of plant immunity genes in two species of Zea

Plant pathogenesis-related genes (PR genes) code for enzymes, enzyme inhibitors, and other peptides that confer resistance to pathogens and herbivores. Although several PR genes have been the subject of molecular population genetic analyses, a general understanding of their long-term evolutionary dynamics remains incomplete. Here we analyze sequence data from 17 PR genes from two closely related teosinte species of central Mexico. In addition to testing whether patterns of diversity at individual loci depart from expectations under a neutral model, we compared patterns of diversity at defense genes, as a class, to nondefense genes. In Zea diploperennis, the majority of defense genes have patterns of diversity consistent with neutral expectations while at least two genes showed evidence of recent positive selection consistent with arms-race models of antagonistic coevolution. In Zea mays ssp. parviglumis, by contrast, analyses of both defense and nondefense genes revealed strong and consistent departures from the neutral model, suggestive of nonequilibrium population dynamics or population structure. Nevertheless, we found a significant excess of replacement polymorphism in defense genes compared to nondefense genes. Although we cannot exclude relaxed selective constraint as an explanation, our results are consistent with temporally variable (transient or episodic) selection or geographically variable selection acting on parviglumis defense genes. The different patterns of diversity found in the two Zea species may be explained by parviglumis' greater distribution and population structure together with geographic variation in selection.

Molecular analysis of the conserved C4 region of the S11-RNase of Solanum chacoense

The stylar component to gametophytic self-incompatibility in Solanaceae is an S-RNase. Its primary structure has a characteristic pattern of two hypervariable regions, involved in pollen recognition, and five constant regions. Two of the latter (C2 and C3) constitute the active site, while the highly hydrophobic C1 and C5 are believed to be involved in protein stability. We analyzed the role of the C4 region by site-directed mutagenesis. A GGGG mutant, in which the four charged residues in the C4 region were replaced with glycine, did not accumulate the protein to detectable levels in styles, suggestive of a role in protein stability. A R115G mutant, in which a charged amino acid was eliminated to reduce the potential binding affinity, had no effect on the pollen rejection phenotype. This suggests the C4 does not interact with partners such as potential pollen tube receptors facilitating S-RNase uptake. Finally, a K113R mutant replaced a potential ubiquitination target with arginine. However, this RNase acted as the wild type in both incompatible and compatible crosses. The latter crosses rule out the role of the conserved C4 lysine in ubiquitination.

Plant self-incompatibility in natural populations: a critical assessment of recent theoretical and empirical advances

Self-incompatibility systems in plants are genetic systems that prevent self-fertilization in hermaphrodites through recognition and rejection of pollen expressing the same allelic specificity as that expressed in the pistils. The evolutionary properties of these self-recognition systems have been revealed through a fascinating interplay between empirical advances and theoretical developments. In 1939, Wright suggested that the main evolutionary force driving the genetic and molecular properties of these systems was strong negative frequency-dependent selection acting on pollination success. The empirical observation of high allelic diversity at the self-incompatibility locus in several species, followed by the discovery of very high molecular divergence among alleles in all plant families where the locus has been identified, supported Wright's initial theoretical predictions as well as many of its later developments. In the last decade, however, advances in the molecular characterization of the incompatibility reaction and in the analysis of allelic frequencies and allelic divergence from natural populations have stimulated new theoretical investigations that challenged some important assumptions of Wright's model of gametophytic self-incompatibility. We here review some of these recent empirical and theoretical advances that investigated: (i) the hypothesis that S-alleles are selectively equivalent, and the evolutionary consequences of genetic interactions between alleles; (ii) the occurrence of frequency-dependent selection in female fertility; (iii) the evolutionary genetics of self-incompatibility systems in subdivided populations; (iv) the evolutionary implications of the self-incompatibility locus's genetic architecture; and (v) of its interactions with the genomic environment.

Sheltered load associated with S-alleles in Solanum carolinense

Bud pollinations allowed me to examine the effects of homozygosity at loci in the area of suppressed recombination around the S-locus in Solanum carolinense, whose S-alleles show surprisingly low diversification rates. The total number of seeds produced was lower for incompatible than compatible pollinations, revealing that self-incompatibility was only somewhat overcome by bud pollination. However, low seed set in incompatible crosses was not due solely to the incompatibility response; crosses between distinct plants sharing the same alleles at the S-locus led to dramatically high seed abortion, nearly equal to that found upon selfing. An excess of heterozygotes in the surviving progeny supports the supposition that these high abortion rates are due to sheltered load, that is, previously unexpressed load accumulated due to enforced heterozygosity and recombination suppression around the S-locus. Of the seven alleles examined in total, two showed evidence of severe load and five did not. The magnitude of load was consistent with terminal branch length in some, but not all, cases.

Allelic diversity of S-RNase at the self-incompatibility locus in natural flowering cherry populations (Prunus lannesiana var. speciosa)

In the Rosaceae family, which includes Prunus, gametophytic self-incompatibility (GSI) is controlled by a single multiallelic locus (S-locus), and the S-locus product expressed in the pistils is a glycoprotein with ribonuclease activity (S-RNase). Two populations of flowering cherry (Prunus lannesiana var. speciosa), located on Hachijo Island in Japan's Izu Islands, were sampled, and S-allele diversity was surveyed based on the sequence polymorphism of S-RNase. A total of seven S-alleles were cloned and sequenced. The S-RNases of flowering cherry showed high homology to those of Prunus cultivars (P. avium and P. dulcis). In the phylogenetic tree, the S-RNases of flowering cherry and other Prunus cultivars formed a distinct group, but they did not form species-specific subgroups. The nucleotide substitution pattern in S-RNases of flowering cherry showed no excess of nonsynonymous substitutions relative to synonymous substitutions. However, the S-RNases of flowering cherry had a higher Ka/Ks ratio than those of other Prunus cultivars, and a subtle heterogeneity in the nucleotide substitution rates was observed among the Prunus species. The S-genotype of each individual was determined by Southern blotting of restriction enzyme-digested genomic DNA, using cDNA for S-RNase as a probe. A total of 22 S-alleles were identified. All individuals examined were heterozygous, as expected under GSI. The allele frequencies were, contrary to the expectation under GSI, significantly unequal. The two populations studied showed a high degree of overlap, with 18 shared alleles. However, the allele frequencies differed considerably between the two populations.

Population genetics of sporophytic self-incompatibility in Senecio squalidus L. (Asteraceae) II: a spatial autocorrelation approach to determining mating behaviour in the presence of low S allele diversity

We recently estimated that as few as six S alleles represent the extent of S locus diversity in a British population of the self-incompatible (SI) coloniser Senecio squalidus (Oxford Ragwort). Despite the predicted constraints to mating imposed by such a low number of S alleles, S. squalidus maintains a strong sporophytic self-incompatibility (SSI) system and there is no evidence for a breakdown of SSI or any obvious negative reproductive consequences for this highly successful coloniser. The present paper assesses mating behaviour in an Oxford S. squalidus population through observations of its effect on spatial patterns of genetic diversity and thus the extent to which it is responsible for ameliorating the potentially detrimental reproductive consequences of low S allele diversity in British S. squalidus. A spatial autocorrelation (SA) treatment of S locus and allozyme polymorphism data for four loci indicates that mating events regularly occur at all the distance classes examined from 60 to 480 m throughout the entire sample population. Less SA is observed for S locus data than for allozyme data in accordance with the hypothesis that SSI and low diversity at the S locus are driving these large-scale mating events. The limited population structure at small distances of 60 m and less observed for SA analysis of the Me-2 locus and by F-statistics for all the allozyme data, is evidence of some local relatedness due to limited seed and pollen dispersal in S. squalidus. However, the overall impression of mating dynamics in this S. squalidus population is that of ample potential mating opportunities with many individuals at large population scales, indicating that reproductive success is not seriously affected by few S alleles available for mating interactions.

Evidence for rare recombination at the gametophytic self-incompatibility locus

The gametophytic self-incompatibility locus has been thought to be a nonrecombining genomic region. Inferences have been made, however, about the functional importance of different parts of the S-locus, based on differences in the levels of variability along the gene, and this is valid only if recombination occurs. It is thus important to test whether recombination occurs within and near the S-locus. Several recent attempts to test this have reached conflicting conclusions. In this study, we examine a large data set on sequence variation at the S-locus in several species with gametophytic self-incompatibility systems, in the Solanaceae, Rosaceae and Scrophulariaceae. We use the longest sequences available to test for recombination based on linkage disequilibrium between polymorphic sites in the S-locus. The relationship between linkage disequilibrium and physical distance between the sites suggests rare intragenic exchange in the evolutionary history of four species of Solanaceae and two species of Rosaceae.