The number of self-incompatibility alleles in a finite, subdivided population

The distribution of self-incompatibility systems in angiosperms: the relationship between mating system diversity, life span, growth habit and latitude in a changing global environment

BACKGROUND AND AIMS

There is ample theoretical and experimental evidence that angiosperms harbouring self-incompatibility (SI) systems are likely to respond to global changes in unique ways relative to taxa with other mating systems. In this paper, we present an updated database on the prevalence of SI systems across angiosperms and examine the relationship between the presence of SI and latitude, biomes, life-history traits and management conditions to evaluate the potential vulnerability of SI taxa to climate change and habitat disturbance.

METHODS

We performed literature searches to identify studies that employed controlled crosses, microscopic analyses and/or genetic data to classify taxa as having SI, self-compatibility (SC), partial self-compatibility (PSC) or self-sterility (SS). Where described, the site of the SI reaction and the presence of dimorphic versus monomorphic flowers were also recorded. We then combined this database on the distribution of mating systems with information about the life span, growth habit, management conditions and geographic distribution of taxa. Information about the geographic distribution of taxa was obtained from a manually curated version of the Global Biodiversity Information Facility database, and from vegetation surveys encompassing nine biomes. We employed multinomial logit regression to assess the relationship between mating system and life-history traits, management condition, latitude and latitude-squared using self-compatible taxa as the baseline. Additionally, we employed LOESS regression to examine the relationship between the probability of SI and latitude. Finally, by summarizing information at the family level, we plotted the distribution of SI systems across angiosperms, including information about the presence of SI or dioecy and the inferred reaction site of the SI system when known, as well as the proportion of taxa in a family for which information is available.

KEY RESULTS

We obtained information about the SI status of 5686 hermaphroditic taxa, of which 55% exhibit SC and the remaining 45% harbour SI, SS or PSC. Highlights of the multinomial logit regression include that taxa with PSC have a greater odds of being short-lived (OR = 1.3) or long-lived (OR = 1.57) perennials relative to SC ones, and that SS/SI taxa (pooled) are less likely to be annuals (OR = 0.64) and more likely to be long-lived perennials (OR = 1.32). SS/SI taxa had a greater odds of being succulent (OR = 2.4) or a tree (OR = 2.05), and were less likely to be weeds (OR = 0.34). Further, we find a quadratic relationship between the probability of being self-incompatible with latitude: SI taxa were more common in the tropics, a finding that was further supported by the vegetation surveys, which showed fewer species with SS/SI in temperate and northern latitudes compared with Mediterranean and tropical biomes.

CONCLUSIONS

We conclude that in the short-term habitat fragmentation, pollinator loss and temperature increases may negatively impact plants with SI systems, particularly long-lived perennial and woody species dominant in tropical forests. In the longer term, these and other global changes are likely to select for self-compatible or partially self-compatible taxa, which, due to the apparent importance of SI as a driver of plant diversification across the angiosperm tree of life, may globally influence plant species richness.

Genetic and Habitat Rescue Improve Population Viability in Self-Incompatible Plants

Habitat fragmentation and the acceleration of environmental change threaten the survival of many plant species. The problem is especially pronounced for plant species with self-incompatibility mating systems, which are obligate outcrossers, thus requiring high mate availability to persist. In such situations, plant populations suffering decreased fitness could be rescued by: (a) improving local habitat conditions (habitat rescue), (b) increasing the number of individuals (demographic rescue), or (c) introducing new genetic variation (genetic rescue). In this study, we used a spatially and genetically explicit individual-based model to approximate the demography of a small (N = 250) isolated self-incompatible population using a timescale of 500 years. Using this model, we quantified the effectiveness of the different types of rescues described above, singly and in combination. Our results show that individual genetic rescue is the most effective type of rescue with respect to improving fitness and population viability. However, we found that introducing a high number of individuals (N > 30) to a small population (N = 50) at the brink of extinction through demographic rescue can also have a positive effect on viability, improving average fitness by 55% compared to introducing a low number of individuals (N = 10) over a long timescale (> 500 years). By itself, habitat rescue showed the lowest effects on viability. However, combining genetic and habitat rescue provided the best results overall, increasing both persistence (> 30%) and mate availability (> 50%). Interestingly, we found that the addition of even a small number of new S alleles (20%) can be highly beneficial to increase mate availability and persistence. We conclude that genetic rescue through the introduction of new S alleles and an increase in habitat suitability is the best management strategy to improve mate availability and population viability of small isolated SI plant populations to overcome the effects of demographic stochasticity and positive density dependence.

Analysis of randomly mutated AlSRKb genes reveals that most loss-of-function mutations cause defects in plasma membrane localization

Only very limited information is available on why some nonsynonymous variants severely alter gene function while others have no effect. To identify the characteristic features of mutations that strongly influence gene function, this study focused on SRK which encodes a highly polymorphic receptor kinase expressed in stigma papillary cells that underlies a female determinant of self-incompatibility in Brassicaceae. A set of 300 Arabidopsis thaliana transformants expressing mutated SRKb from A. lyrata was constructed using error-prone PCR and the genotype and self-incompatibility phenotype of each transformant were determined. Almost all the transformants showing the self-incompatibility defect contained mutations in AlSRKb that altered localization to the plasma membrane. The observed mutations occurred in amino acid residues that were highly conserved across S haplotypes and whose predicted locations were in the interior of the protein. Our findings suggested that mutations causing the self-incompatibility defect were more likely to result from changes to AlSRKb biosynthesis than from loss of AlSRKb function. In addition, we examined whether the RandomForest and Extreme Gradient Boosting methods could predict the self-incompatibility phenotypes of SRK mutants.

Self-incompatibility phenotypes of SRK mutants can be predicted with high accuracy

Only very limited information is available on why some non-synonymous variants severely alter gene function while others have no effect. To identify the characteristic features of mutations that strongly influence gene function, this study focused on S-locus receptor kinase, SRK, which encodes a highly polymorphic receptor kinase expressed in stigma papillary cells that underlies a female determinant of self-incompatibility in Brassicaceae. A set of 299 Arabidopsis thaliana transformants expressing mutated SRKb from A. lyrata was constructed and analyzed to determine the genotype and self-incompatibility phenotype of each transformant. Almost all the transformants showing the self-incompatibility defect contained mutations in AlSRKb that altered localization to the plasma membrane. The observed mutations occurred in amino acid residues that were highly conserved across S haplotypes and whose predicted locations were in the interior of the protein. These mutations were likely to underlie the self-incompatibility defect as they caused significant changes to amino acid properties. Such findings suggested that mutations causing the self-incompatibility defect were more likely to result from changes to AlSRKb biosynthesis than from loss of function. In addition, this study showed the RandomForest and Extreme Gradient Boosting methods could predict self-incompatibility phenotypes of SRK mutants with high accuracy.

Revisiting the number of self-incompatibility alleles in finite populations: From old models to new results

Under gametophytic self-incompatibility (GSI), plants are heterozygous at the self-incompatibility locus (S-locus) and can only be fertilized by pollen with a different allele at that locus. The last century has seen a heated debate about the correct way of modelling the allele diversity in a GSI population that was never formally resolved. Starting from an individual-based model, we derive the deterministic dynamics as proposed by Fisher (The genetical theory of natural selection - A complete, Variorum edition, Oxford University Press, 1958) and compute the stationary S-allele frequency distribution. We find that the stationary distribution proposed by Wright (Evolution, 18, 609, 1964) is close to our theoretical prediction, in line with earlier numerical confirmation. Additionally, we approximate the invasion probability of a new S-allele, which scales inversely with the number of resident S-alleles. Lastly, we use the stationary allele frequency distribution to estimate the population size of a plant population from an empirically obtained allele frequency spectrum, which complements the existing estimator of the number of S-alleles. Our expression of the stationary distribution resolves the long-standing debate about the correct approximation of the number of S-alleles and paves the way for new statistical developments for the estimation of the plant population size based on S-allele frequencies.

Non-self recognition-based self-incompatibility can alternatively promote or prevent introgression

Self-incompatibility alleles (S-alleles), which prevent self-fertilisation in plants, have historically been expected to benefit from negative frequency-dependent selection and invade when introduced to a new population through gene flow. However, the most taxonomically widespread form of self-incompatibility, the ribonuclease-based system ancestral to the core eudicots, functions through collaborative non-self recognition, which can affect both short-term patterns of gene flow and the long-term process of S-allele diversification. We analysed a model of S-allele evolution in two populations connected by migration, focussing on comparisons among the fates of S-alleles initially unique to each population and those shared among populations. We found that both shared and unique S-alleles from the population with more unique S-alleles were usually fitter compared with S-alleles from the population with fewer S-alleles. Resident S-alleles often became extinct and were replaced by migrant S-alleles, although this outcome could be averted by pollen limitation or biased migration. Collaborative non-self recognition will usually either result in the whole-sale replacement of S-alleles from one population with those from another or else disfavour introgression of S-alleles altogether.

Massive introgression of major histocompatibility complex (MHC) genes in newt hybrid zones

Variation in the vertebrate major histocompatibility complex (MHC) genes is crucial for fighting pathogen assault. Because new alleles confer a selective advantage, MHC should readily introgress between species, even under limited hybridization. Using replicated transects through two hybrid zones between strongly reproductively isolated European newts, Lissotriton montandoni and L. vulgaris, we demonstrated recent and ongoing MHC class I and II introgression in the Carpathian region. The extent of introgression correlated with the age of contact. In the older zone, MHC similarity between species within transects exceeded similarity between transects within species, implying pervasive introgression - a massive exchange of MHC genes, not limited to specific variants. In simulations, the observed pattern emerged under the combined action of balancing selection and hybridization, but not when these processes acted separately. Thus, massive introgression at advanced stages of divergence can introduce novel and restore previously lost MHC variation, boosting the adaptive potential of hybridizing taxa. In consequence, MHC genes may be the last to stop introgressing between incipient species.

Positive selection and climatic effects on MHC class II gene diversity in hares (Lepus capensis) from a steep ecological gradient

In natural populations, allelic diversity of the major histocompatibility complex (MHC) is commonly interpreted as resulting from positive selection in varying spatiotemporal pathogenic landscapes. Composite pathogenic landscape data are, however, rarely available. We studied the spatial distribution of allelic diversity at two MHC class II loci (DQA, DQB) in hares, Lepus capensis, along a steep ecological gradient in North Africa and tested the role of climatic parameters for the spatial distribution of DQA and DQB proteins. Climatic parameters were considered to reflect to some extent pathogenic landscape variation. We investigated historical and contemporary forces that have shaped the variability at both genes, and tested for differential selective pressure across the ecological gradient by comparing allelic variation at MHC and neutral loci. We found positive selection on both MHC loci and significantly decreasing diversity from North to South Tunisia. Our multinomial linear models revealed significant effects of geographical positions that were correlated with mean annual temperature and precipitation on the occurrence of protein variants, but no effects of co-occurring DQA or DQB proteins, respectively. Diversifying selection, recombination, adaptation to local pathogenic landscapes (supposedly reflected by climate parameters) and neutral demographic processes have shaped the observed MHC diversity and differentiation patterns.

Forward-in-Time, Spatially Explicit Modeling Software to Simulate Genetic Lineages Under Selection

SELECTOR is a software package for studying the evolution of multiallelic genes under balancing or positive selection while simulating complex evolutionary scenarios that integrate demographic growth and migration in a spatially explicit population framework. Parameters can be varied both in space and time to account for geographical, environmental, and cultural heterogeneity. SELECTOR can be used within an approximate Bayesian computation estimation framework. We first describe the principles of SELECTOR and validate the algorithms by comparing its outputs for simple models with theoretical expectations. Then, we show how it can be used to investigate genetic differentiation of loci under balancing selection in interconnected demes with spatially heterogeneous gene flow. We identify situations in which balancing selection reduces genetic differentiation between population groups compared with neutrality and explain conflicting outcomes observed for human leukocyte antigen loci. These results and three previously published applications demonstrate that SELECTOR is efficient and robust for building insight into human settlement history and evolution.

Paternal-specific S-allele transmission in sweet cherry (Prunus avium L.): the potential for sexual selection

Homomorphic self-incompatibility is a well-studied example of a physiological process that is thought to increase population diversity and reduce the expression of inbreeding depression. Whereas theoretical models predict the presence of a large number of S-haplotypes with equal frequencies at equilibrium, unequal allele frequencies have been repeatedly reported and attributed to sampling effects, population structure, demographic perturbation, sheltered deleterious mutations or selection pressure on linked genes. However, it is unclear to what extent unequal segregations are the results of gametophytic or sexual selection. Although these two forces are difficult to disentangle, testing S-alleles in the offspring of controlled crosses provides an opportunity to separate these two phenomena. In this work, segregation and transmission of S-alleles have been characterized in progenies of mixed donors and fully compatible pollinations under field conditions in Prunus avium. Seed set patterns and pollen performance have also been characterized. The results reveal paternal-specific distorted transmission of S-alleles in most of the crosses. Interestingly, S-allele segregation within any given paternal or maternal S-locus was random. Observations on pollen germination, pollen tube growth rate, pollen tube cohort size, seed set dynamics and transmission patterns strongly suggest post-pollination, prezygotic sexual selection, with male-male competition as the most likely mechanism. According to these results, post-pollination sexual selection takes precedence over frequency-dependent selection in explaining unequal S-haplotype frequencies.

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.

A review of molecular approaches for investigating patterns of coevolution in marine host-parasite relationships

Parasites and their relationships with hosts play a crucial role in the evolutionary pathways of every living organism. One method of investigating host-parasite systems is using a molecular approach. This is particularly important as analyses based solely on morphology or laboratory studies of parasites and their hosts do not take into account historical evolutionary interactions that can shape the distribution, abundance and population structure of parasites and their hosts. However, the predominant host-parasite coevolution literature has focused on terrestrial hosts and their parasites, and there still is a lack of studies in marine environments. Given that marine systems are generally more open than terrestrial ones, they provide fascinating opportunities for large-scale (as well as small-scale) geographic studies. Further, patterns and processes of genetic structuring and systematics are becoming more available across many different taxa (but especially fishes) in many marine systems, providing an excellent basis for examining whether parasites follow host population/species structure. In this chapter, we first highlight the factors and processes that challenge our ability to interpret evolutionary patterns of coevolution of hosts and their parasites in marine systems at different spatial, temporal and taxonomic scales. We then review the use of the most commonly utilized genetic markers in studying marine host-parasite systems. We give an overview and discuss which molecular methodologies resolve evolutionary relationships best and also discuss the applicability of new approaches, such as next-generation sequencing and studies utilizing functional markers to gain insights into more contemporary processes shaping host-parasite relationships.

Life history mediates mate limitation and population viability in self-incompatible plant species

Genetically controlled self-incompatibility systems represent links between genetic diversity and plant demography with the potential to directly impact on population dynamics. We use an individual-based spatial simulation to investigate the demographic and genetic consequences of different self-incompatibility systems for plants that vary in reproductive capacity and lifespan. The results support the idea that, in the absence of inbreeding effects, populations of self-incompatible species will often be smaller and less viable than self-compatible species, particularly for shorter-lived organisms or where potential fecundity is low. At high ovule production and low mortality, self-incompatible and self-compatible species are demographically similar, thus self-incompatibility does not automatically lead to reduced mate availability or population viability. Overall, sporophytic codominant self-incompatibility was more limiting than gametophytic or sporophytic dominant systems, which generally behaved in a similar fashion. Under a narrow range of conditions, the sporophytic dominant system maintained marginally greater mate availability owing to the production of S locus homozygotes. While self-incompatibility reduces population size and persistence for a broad range of conditions, the actual number of S alleles, beyond that required for reproduction, is important for only a subset of life histories. For these situations, results suggest that addition of new S alleles may result in significant demographic rescue.

Pollen limitation and reduced reproductive success are associated with local genetic effects in Prunus virginiana, a widely distributed self-incompatible shrub

BACKGROUND AND AIMS

A vast quantity of empirical evidence suggests that insufficient quantity or quality of pollen may lead to a reduction in fruit set, in particular for self-incompatible species. This study uses an integrative approach that combines field research with marker gene analysis to understand the factors affecting reproductive success in a widely distributed self-incompatible species, Prunus virginiana (Rosaceae).

METHODS

Twelve patches of P. virginiana distributed within three populations that differed in degree of disturbance were examined. Two of the sites were small (7-35 km(2)) remnants of forest in an intensively used agricultural landscape, while the third was continuous (350 km(2)) and less disturbed. Field studies (natural and hand cross-pollinations) were combined with marker gene analyses (microsatellites and S-locus) in order to explore potential factors affecting pollen delivery and consequently reproductive success at landscape (between populations) and fine scales (within populations).

KEY RESULTS

Reductions in reproductive output were found in the two fragments compared with the continuous population, and suggest that pollen is an important factor limiting fruit production. Genetic analyses carried out in one of the fragments and in the continuous site suggest that even though S-allele diversity is high in both populations, the fragment exhibits an increase in biparental inbreeding and correlated paternity. The increase in biparental inbreeding in the fragment is potentially attributable to variation in the density of individuals and/or the spatial distribution of genotypes among populations, both of which could alter mating dynamics.

CONCLUSIONS

By using a novel integrative approach, this study shows that even though P. virginiana is a widespread species, fragmented populations can experience significant reductions in fruit set and pollen limitation in the field. Deatiled examination of one fragmented population suggests that these linitations may be explained by an increase in biparental inbreeding, correlated paternity and fine-scale genetic structure. The consistency of the field and fine-scale genetic analyses, and the consistency of the results within patches and across years, suggest that these are important processes driving pollen limitation in the fragment.

Inference of selection based on temporal genetic differentiation in the study of highly polymorphic multigene families

The co-evolutionary arms race between host immune genes and parasite virulence genes is known as Red Queen dynamics. Temporal fluctuations in allele frequencies, or the ‘turnover’ of alleles at immune genes, are concordant with predictions of the Red Queen hypothesis. Such observations are often taken as evidence of host-parasite co-evolution. Here, we use computer simulations of the Major Histocompatibility Complex (MHC) of guppies (Poecilia reticulata) to study the turnover rate of alleles (temporal genetic differentiation, G'_ST). Temporal fluctuations in MHC allele frequencies can be ≥≤order of magnitude larger than changes observed at neutral loci. Although such large fluctuations in the MHC are consistent with Red Queen dynamics, simulations show that other demographic and population genetic processes can account for this observation, these include: (1) overdominant selection, (2) fluctuating population size within a metapopulation, and (3) the number of novel MHC alleles introduced by immigrants when there are multiple duplicated genes. Synergy between these forces combined with migration rate and the effective population size can drive the rapid turnover in MHC alleles. We posit that rapid allelic turnover is an inherent property of highly polymorphic multigene families and that it cannot be taken as evidence of Red Queen dynamics. Furthermore, combining temporal samples in spatial F_ST outlier analysis may obscure the signal of selection.

Impact of negative frequency-dependent selection on mating pattern and genetic structure: a comparative analysis of the S-locus and nuclear SSR loci in Prunus lannesiana var. speciosa

Mating processes of local demes and spatial genetic structure of island populations at the self-incompatibility (S-) locus under negative frequency-dependent selection (NFDS) were evaluated in Prunus lannesiana var. speciosa in comparison with nuclear simple sequence repeat (SSR) loci that seemed to be evolutionarily neutral. Our observations of local mating patterns indicated that male-female pair fecundity was influenced by not only self-incompatibility, but also various factors, such as kinship, pollen production and flowering synchrony. In spite of the mating bias caused by these factors, the NFDS effect on changes in allele frequencies from potential mates to mating pollen was detected at the S-locus but not at the SSR loci, although the changes from adult to juvenile cohorts were not apparent at any loci. Genetic differentiation and isolation-by-distance over various spatial scales were smaller at the S-locus than at the SSR loci, as expected under the NFDS. Allele-sharing distributions among the populations also had a unimodal pattern at the S-locus, indicating the NFDS effect except for alleles unique to individual populations probably due to isolation among islands, although this pattern was not exhibited by the SSR loci. Our results suggest that the NFDS at the S-locus has an impact on both the mating patterns and the genetic structure in the P. lannesiana populations studied.

Microevolution of S-allele frequencies in wild cherry populations: respective impacts of negative frequency dependent selection and genetic drift

Negative frequency dependent selection (NFDS) is supposed to be the main force controlling allele evolution at the gametophytic self-incompatibility locus (S-locus) in strictly outcrossing species. Genetic drift also influences S-allele evolution. In perennial sessile organisms, evolution of allelic frequencies over two generations is mainly shaped by individual fecundities and spatial processes. Using wild cherry populations between two successive generations, we tested whether S-alleles evolved following NFDS qualitative and quantitative predictions. We showed that allelic variation was negatively correlated with parental allelic frequency as expected under NFDS. However, NFDS predictions in finite population failed to predict more than half S-allele quantitative evolution. We developed a spatially explicit mating model that included the S-locus. We studied the effects of self-incompatibility and local drift within populations due to pollen dispersal in spatially distributed individuals, and variation in male fecundity on male mating success and allelic frequency evolution. Male mating success was negatively related to male allelic frequency as expected under NFDS. Spatial genetic structure combined with self-incompatibility resulted in higher effective pollen dispersal. Limited pollen dispersal in structured distributions of individuals and genotypes and unequal pollen production significantly contributed to S-allele frequency evolution by creating local drift effects strong enough to counteract the NFDS effect on some alleles.

Sporophytic self-incompatibility in Senecio squalidus (Asteraceae): S allele dominance interactions and modifiers of cross-compatibility and selfing rates

Understanding genetic mechanisms of self-incompatibility (SI) and how they evolve is central to understanding the mating behaviour of most outbreeding angiosperms. Sporophytic SI (SSI) is controlled by a single multi-allelic locus, S, which is expressed in the diploid (sporophyte) plant to determine the SI phenotype of its haploid (gametophyte) pollen. This allows complex patterns of independent S allele dominance interactions in male (pollen) and female (pistil) reproductive tissues. Senecio squalidus is a useful model for studying the genetic regulation and evolution of SSI because of its population history as an alien invasive species in the UK. S. squalidus maintains a small number of S alleles (7-11) with a high frequency of dominance interactions. Some S. squalidus individuals also show partial selfing and/or greater levels of cross-compatibility than expected under SSI. We previously speculated that these might be adaptations to invasiveness. Here we describe a detailed characterization of the regulation of SSI in S. squalidus. Controlled crosses were used to determine the S allele dominance hierarchy of six S alleles and effects of modifiers on cross-compatibility and partial selfing. Complex dominance interactions among S alleles were found with at least three levels of dominance and tissue-specific codominance. Evidence for S gene modifiers that increase selfing and/or cross-compatibility was also found. These empirical findings are discussed in the context of theoretical predictions for maintenance of S allele dominance interactions, and the role of modifier loci in the evolution of SI.

Expression and inheritance of sporophytic self-incompatibility in synthetic allohexaploid Senecio cambrensis (Asteraceae)

Allopolyploid speciation is common in plants and is frequently associated with shifts from outcrossing, for example self-incompatibility, to inbreeding (i.e. selfing). Senecio cambrensis is a recently evolved allohexaploid species that formed following hybridization between diploid self-incompatible S. squalidus and tetraploid self-compatible S. vulgaris. Studies of reproduction in wild populations of S. cambrensis have concluded that it is self-compatible. Here, we investigated self-compatibility in synthetic lines of S. cambrensis generated via hybridization and colchicine-induced polyploidization and wild S. cambrensis using controlled crossing experiments. Synthetic F(1)S. cambrensis individuals were all self-compatible but, in F(2) and later generations, self-incompatible individuals were identified at frequencies of 6.7-9.2%. Self-incompatibility was also detected in wild sampled individuals at a frequency of 12.2%. The mechanism and genetics of self-incompatibility were tested in synthetic S. cambrensis and found to be similar to those of its paternal parent S. squalidus (i.e. sporophytic). These results show, for the first time, that functional sporophytic self-incompatibility can be inherited and expressed in allopolyploids as early as the second (F(2)) generation. Wild S. cambrensis should therefore be considered as possessing a mixed mating system with the potential for evolution towards either inbreeding or outcrossing.

Modeling multiallelic selection using a Moran model

We present a Moran-model approach to modeling general multiallelic selection in a finite population and show how it may be used to develop theoretical models of biological systems of balancing selection such as plant gametophytic self-incompatibility loci. We propose new expressions for the stationary distribution of allele frequencies under selection and use them to show that the continuous-time Markov chain describing allele frequency change with exchangeable selection and Moran-model reproduction is reversible. We then use the reversibility property to derive the expected allele frequency spectrum in a finite population for several general models of multiallelic selection. Using simulations, we show that our approach is valid over a broader range of parameters than previous analyses of balancing selection based on diffusion approximations to the Wright-Fisher model of reproduction. Our results can be applied to any model of multiallelic selection in which fitness is solely a function of allele frequency.

Effect of variation in self-incompatibility on pollen limitation and inbreeding depression in Flourensia cernua (Asteraceae) scrubs of contrasting density

BACKGROUND AND AIMS

Selection may favour a partial or complete loss of self-incompatibility (SI) if it increases the reproductive output of individuals in the presence of low mate availability. The reproductive output of individuals varying in their strength of SI may also be affected by population density via its affect on the spatial structuring and number of S-alleles in populations. Modifiers increasing levels of self-compatibility can be selected when self-compatible individuals receive reproductive compensation by, for example, increasing seed set and/or when they become associated with high fitness genotypes.

METHODS

The effect of variation in the strength of SI and scrub density (low versus high) on seed set, seed germination and inbreeding depression in seed germination (delta(germ)) was investigated in the partially self-incompatible species Flourensia cernua by analysing data from self-, cross- and open-pollinated florets.

KEY RESULTS

Examination of 100 plants in both high and low scrub densities revealed that 51% of plants were strongly self-incompatible and 49 % varied from being self-incompatible to self-compatible. Seed set after hand cross-pollination was higher than after open-pollination for self-incompatible, partially self-incompatible and self-compatible plants but was uniformly low for strongly self-incompatible plants. Strongly self-incompatible and self-incompatible plants exhibited lower seed set, seed germination and multiplicative female fitness (floral display x seed set x seed germination) in open-pollinated florets compared with partially self-incompatible and self-compatible plants. Scrub density also had an effect on seed set and inbreeding depression: in low-density scrubs seed set was higher after open-pollination and delta(germ) was lower.

CONCLUSIONS

These data suggest that (a) plants suffered outcross pollen limitation, (b) female fitness in partially self-incompatible and self-compatible plants is enhanced by increased mate-compatibility and (c) plants in low-density scrubs received higher quality pollen via open-pollination than plants in high-density scrubs.

Antagonism between local dispersal and self-incompatibility systems in a continuous plant population

Many self-incompatible plant species exist in continuous populations in which individuals disperse locally. Local dispersal of pollen and seeds facilitates inbreeding because pollen pools are likely to contain relatives. Self-incompatibility promotes outbreeding because relatives are likely to carry incompatible alleles. Therefore, populations can experience an antagonism between these forces. In this study, a novel computational model is used to explore the effects of this antagonism on gene flow, allelic diversity, neighbourhood sizes, and identity by descent. I confirm that this antagonism is sensitive to dispersal levels and linkage. However, the results suggest that there is little to no difference between the effects of gametophytic and sporophytic self-incompatibility systems (GSI and SSI) on unlinked loci. More importantly, both GSI and SSI affect unlinked loci in a manner similar to obligate outcrossing without mating types. This suggests that the primary evolutionary impact of self-incompatibility systems may be to prevent selfing, and prevention of biparental inbreeding might be a beneficial side-effect.

Molecular population genetics of the SRK and SCR self-incompatibility genes in the wild plant species Brassica cretica (Brassicaceae)

Self-incompatibility (SI) in plants is a classic example of a trait evolving under strong frequency-dependent selection. As a consequence, population genetic theory predicts that the S locus, which controls SI, should maintain numerous alleles, display a high level of nucleotide diversity, and, in structured populations, show a lower level of among-population differentiation compared to neutral loci. Population-level investigations of DNA sequence variation at the S locus have recently been carried out in the genus Arabidopsis, largely confirming results from theoretical models of S-locus evolutionary dynamics, but no comparable studies have been done in wild Brassica species. In this study, we sequenced parts of the S-locus genes SRK and SCR, two tightly linked genes that are directly involved in the determination of SI specificity in samples from four natural populations of the wild species Brassica cretica. The amount and distribution of nucleotide diversity, as well as the frequency spectrum of putative functional haplotypes, observed at the S locus in B. cretica fit very well with expectations from theoretical models, providing strong evidence for frequency-dependent selection acting on the S locus in a wild Brassica species.

Selection at work in self-incompatible Arabidopsis lyrata. II. Spatial distribution of S haplotypes in Iceland

We survey the distribution of haplotypes at the self-incompatibility (SI) locus of Arabidopsis lyrata (Brassicaceae) at 12 locations spread over the species' natural distribution in Iceland. Previous investigations of the system have identified 34 functionally different S haplotypes maintained by frequency-dependent selection and arranged them into four classes of dominance in their phenotypic expression. On the basis of this model of dominance and the island model of population subdivision, we compare the distribution of S haplotypes with that expected from population genetic theory. We observe 18 different S haplotypes, recessive haplotypes being more common than dominant ones, and dominant ones being shared by fewer populations. As expected, differentiation, although significant, is very low at the S locus even over distances of up to 300 km. The frequency of the most recessive haplotype is slightly larger than expected for a panmictic population, but consistent with a subdivided population with the observed differentiation. Frequencies in nature reflect effects of segregation distortion previously observed in controlled crosses. The dynamics of the S-locus variation are, however, well represented by a 12-island model and our simplified model of dominance interactions.

Repeated adaptive introgression at a gene under multiallelic balancing selection

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, Arabidopsis halleri and A. lyrata, 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 (K = 0.0193) is about four times lower than the genomic background (K = 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.

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 Arabidopsis 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.

Unequal segregation of SRK alleles at the S locus in Brassica cretica

In the Brassicaceae plant family, which includes the Arabidopsis and Brassica genera, self-incompatibility (SI) is controlled by genes at the S locus. Using experimental crosses, we studied the pattern of inheritance of S-locus alleles in the wild species Brassica cretica. Four full-sib families were established and unequal segregation of alleles at the SRK SI gene was found in one family. The segregation distortion acted in favour of a recessive (class II) allele and was best explained by some form of gametic-level selection. Our findings are discussed in the light of theoretical predictions of differential accumulation of deleterious mutations among S-locus alleles.

Genomic consequences of selection on self-incompatibility genes

Frequency-dependent selection at plant self-incompatibility systems is inherent and well understood theoretically. A self-incompatibility locus leads to a strong peak of diversity in the genome, to a unique distribution of diversity across the species and possibly to increased introgression between closely related species. We review recent empirical studies demonstrating these features and relate the empirical findings to theoretical predictions. We show how these features are being exploited in searches for other genes under multi-allelic balancing selection and for inference on recent breakdown of self-incompatibility.

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.

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

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

Patterns of genetic diversity in outcrossing and selfing populations ofArabidopsis lyrata

Arabidopsis lyrata is normally considered an obligately outcrossing species with a strong self-incompatibility system, but a shift in mating system towards inbreeding has been found in some North American populations (subspecies A. lyrata ssp. lyrata). This study provides a survey of the Great Lakes region of Canada to determine the extent of this mating system variation and how outcrossing rates are related to current population density, geographical distribution, and genetic diversity. Based on variation at microsatellite markers (progeny arrays to estimate multilocus outcrossing rates and population samples to estimate diversity measures) and controlled greenhouse pollinations, populations can be divided into two groups: (i) group A, consisting of individuals capable of setting selfed seed (including autogamous fruit set in the absence of pollinators), showing depressed outcrossing rates (T(m) = 0.2-0.6), heterozygosity (H(O) = 0.02-0.06) and genetic diversity (H(E) = 0.08-0.10); and (ii) group B, consisting of individuals that are predominantly self-incompatible (T(m) > 0.8), require pollinators for seeds set, and showing higher levels of heterozygosity (H(O) = 0.13-0.31) and diversity (H(E) = 0.19-0.410). Current population density is not related to the shift in mating system but does vary with latitude. Restricted gene flow among populations was evident among all but two populations (F(ST) = 0.11-0.8). Group A populations were more differentiated from one another (F(ST) = 0.78) than they were from group B populations (F(ST) = 0.59), with 41% of the variation partitioned within populations, 47% between populations, and 12% between groups. No significant relationship was found between genetic and geographical distance. Results are discussed in the context of possible postglacial expansion scenarios in relation to loss of self-incompatibility.

Linkage disequilibrium and recombination rate estimates in the self-incompatibility region of Arabidopsis lyrata

Genetic diversity is unusually high at loci in the S-locus region of the self-incompatible species of the flowering plant, Arabidopsis lyrata, not just in the S loci themselves, but also at two nearby loci. In a previous study of a single natural population from Iceland, we attributed this elevated polymorphism to linkage disequilibrium (LD) between variants at loci close to the S locus and the S alleles, which are maintained in the population by balancing selection. With the four S-flanking loci whose diversity we previously studied, we could not determine the extent of the region linked to the S loci in which neutral sites are affected. We also could not exclude the possibility of a population bottleneck, or of admixture, as causes of the LD. We have now studied four more distant loci flanking the S-locus region, and more populations, and we analyze the results using a theoretical model of the effect of balancing selection on diversity at linked neutral sites within and between different functional S-allelic classes. In the model, diversity is a function of the number of selectively maintained alleles and the recombination distances from the selectively maintained sites. We use the model to estimate the number of different functional S alleles, their turnover rate, and recombination rates between the S-locus region and other loci. Our estimates suggest that there is a small region of very low recombination surrounding the S-locus region.

Population structure at the S-locus of Sorbus aucuparia L. (Rosaceae: Maloideae)

Low sequence divergence within functional alleles is predicted for the self-incompatibility locus because of strong negative frequency-dependent selection. Nevertheless, sequence variation within functional alleles is essential for current models of the evolution of new mating types. We genotyped the stylar self-incompatibility RNase of 20 Sorbus aucuparia from a population in the Pyrenees mountains of France in order to compare alleles found there to those previously sampled in a Belgian population. Both populations returned 20 different alleles from samples of 20 individuals, providing maximum-likelihood estimates of 24.4 (95% CI 20-34) alleles in each. Ten alleles occurred in both samples. The maximum likelihood (ML) estimate of the overlap in the alleles present in both populations was 16, meaning that an estimated eight alleles are private to each population, and a total of 32 alleles occur across the two populations examined. We used Fisher's (1961) missing plot method to estimate that 40 alleles occur in the species. In accord with population genetics theory, we observed at most one synonymous sequence difference between copies of alleles sampled from the different populations and no variation within populations. Phylogenetic analysis shows that nearly every allele in S. aucuparia arose prior to divergence of this species from members of three different genera of the Rosaceae subfamily, Maloideae. Lack of observable sequence variation within alleles, coupled with the slow pace of allelic relative to taxonomic diversification, implies that finding intermediate stages in the process of new allele creation will be difficult in this group.

A general model to explore complex dominance patterns in plant sporophytic self-incompatibility systems

We developed a general model of sporophytic self-incompatibility under negative frequency-dependent selection allowing complex patterns of dominance among alleles. We used this model deterministically to investigate the effects on equilibrium allelic frequencies of the number of dominance classes, the number of alleles per dominance class, the asymmetry in dominance expression between pollen and pistil, and whether selection acts on male fitness only or both on male and on female fitnesses. We show that the so-called "recessive effect" occurs under a wide variety of situations. We found emerging properties of finite population models with several alleles per dominance class such as that higher numbers of alleles are maintained in more dominant classes and that the number of dominance classes can evolve. We also investigated the occurrence of homozygous genotypes and found that substantial proportions of those can occur for the most recessive alleles. We used the model for two species with complex dominance patterns to test whether allelic frequencies in natural populations are in agreement with the distribution predicted by our model. We suggest that the model can be used to test explicitly for additional, allele-specific, selective forces.

Distribution of S-alleles in island populations of flowering cherry, Prunus lannesiana var. speciosa

We surveyed the distribution of S-alleles in natural island populations of Prunus lannesiana var. speciosa sampled from seven sites on the Izu Peninsula and six Izu islands, Japan. The S-genotypes of sampled individuals were determined by Southern analysis of RFLPs generated by restriction enzyme digestion of genomic DNA, using cDNA of the S-RNase gene as a probe. All individuals were heterozygous, as expected under gametophytic self-incompatibility (GSI). Sixty-three S-alleles were observed in the species, but 12 private to the Izu Peninsula population seemed to be derived from related species, giving a total of 75. The estimated number of S-alleles in each population ranged from 26 to 62, and was inversely correlated with the respective population's distance from the Izu Peninsula, the closest point in the mainland to the islands. This geographical cline in the estimated numbers of S-alleles suggests that gene flow to and from the distant island populations was less frequent, and that the studied species has migrated from the mainland to the Izu islands. The genetic relationship at the S-locus among populations also gave an "isolation by distance" pattern. The genetic differentiation at the S-locus among the populations was very low (F(ST) = 0.014, p < 0.001). The number of S-alleles in the species did not seem to depend on genetic differences associated with population subdivisions. This might be due to the greater effective migration rates of S-alleles, as expected under balancing selection in GSI.

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 low recombination affect diversity near the self-incompatibility loci of the plant Arabidopsis lyrata

The self-incompatibility (S-) locus region of plants in the Brassica family is a small genome region. In Arabidopsis lyrata, the S-genes, SRK and SCR, encode the functional female and pollen recognition proteins, which must be coadapted to maintain correct associations between the two component genes, and thus self-incompatibility (SI). Recombinants would be self-compatible and thus probably disadvantageous in self-incompatible species. Therefore, tight linkage between the two genes in incompatibility systems is predicted to evolve to avoid producing such recombinant haplotypes. The evolution of low recombination in S-locus regions has not been rigorously tested. To test whether these regions' per-nucleotide recombination rates differ from those elsewhere in the genome, and to investigate whether the A. lyrata S-loci have the predicted effect on diversity in their immediate genome region, we studied diversity in genes that are linked to the S-loci but are not involved in incompatibility and are not under balancing selection. Compared with other A. lyrata loci, genes linked to the S-loci have extraordinarily high polymorphism. Our estimated recombination in this region, from fitting a model of the effects of S-allele polymorphism on linked neutral sites, supports the hypothesis of locally suppressed recombination around the S-locus.

Linkage disequilibrium between incompatibility locus region genes in the plant Arabidopsis lyrata

We have studied diversity in Arabidopsis lyrata of sequences orthologous to the ARK3 gene of A. thaliana. Our main goal was to test for recombination in the S-locus region. In A. thaliana, the single-copy ARK3 gene is closely linked to the non-functional copies of the self-incompatibility loci, and the ortholog in A. lyrata (a self-incompatible species) is in the homologous genome region and is known as Aly8. It is thus of interest to test whether Aly8 sequence diversity is elevated due to close linkage to the highly polymorphic incompatibility locus, as is theoretically predicted. However, Aly8 is not a single-copy gene, and the presence of paralogs could also lead to the appearance of elevated diversity. We established a typing approach based on different lengths of Aly8 PCR products and show that most A. lyrata haplotypes have a single copy, but some have two gene copies, both closely linked to the incompatibility locus, one being a pseudogene. We determined the phase of multiple haplotypes in families of plants from Icelandic and other populations. Different Aly8 sequence types are associated with different SRK alleles, while haplotypes with the same SRK sequences tend to have the same Aly8 sequence. There is evidence of some exchange of sequences between different Aly8 sequences, making it difficult to determine which ones are allelic or to estimate the diversity. However, the homogeneity of the Aly8 sequences of each S-haplotype suggests that recombination between the loci has been very infrequent over the evolutionary history of these populations. Overall, the results suggest that recombination rarely occurs in the interval between the S-loci and Aly8 and that linkage to the S-loci can probably account for the observed high Aly8 diversity.

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.

The evolutionary ecology of the major histocompatibility complex

The major histocompatibility complex (MHC) has become a paradigm for how selection can act to maintain adaptively important genetic diversity in natural populations. Here, we review the contribution of studies on the MHC in non-model species to our understanding of how selection affects MHC diversity, emphasising how ecological and ethological processes influence the tempo and mode of evolution at the MHC, and conversely, how variability at the MHC affects individual fitness, population dynamics and viability. We focus on three main areas: the types of information that have been used to detect the action of selection on MHC genes; the relative contributions of parasite-mediated and sexual selection on the maintenance of MHC diversity; and possible future lines of research that may help resolve some of the unanswered issues associated with MHC evolution.

Balancing selection in the wild: testing population genetics theory of self-incompatibility in the rare species Brassica insularis

Self-incompatibility (SI) systems are widespread mechanisms that prevent self-fertilization in angiosperms. They are generally encoded by one genome region containing several multiallelic genes, usually called the S-locus. They involve a recognition step between the pollen and the pistil component and pollen is rejected when it shares alleles with the pistil. The direct consequence is that rare alleles are favored, such that the S-alleles are subject to negative frequency-dependent selection. Several theoretical articles have predicted the specific patterns of polymorphism, compared to neutral loci, expected for such genes under balancing selection. For instance, many more alleles should be maintained and populations should be less differentiated than for neutral loci. However, empirical tests of these predictions in natural populations have remained scarce. Here, we compare the genetic structure at the S-locus and microsatellite markers for five natural populations of the rare species Brassica insularis. As in other Brassica species, B. insularis has a sporophytic SI system for which molecular markers are available. Our results match well the theoretical predictions and constitute the first general comparison of S-allele and neutral polymorphism.

Phenotypic and genotypic expression of self-incompatibility haplotypes in Arabidopsis lyrata suggests unique origin of alleles in different dominance classes

The highly divergent alleles of the SRK gene in outcrossing Arabidopsis lyrata have provided important insights into the evolutionary history of self-incompatibility (SI) alleles and serve as an ideal model for studies of the evolutionary and molecular interactions between alleles in cell-cell recognition systems in general. One tantalizing question is how new specificities arise in systems that require coordination between male and female components. Allelic recruitment via gene conversion has been proposed as one possibility, based on the division of DNA sequences at the SRK locus into two distinctive groups: (1) sequences whose relationships are not well resolved and display the long branch lengths expected for a gene under balancing selection (Class A); and (2) sequences falling into a well-supported group with shorter branch lengths (Class B) that are closely related to an unlinked paralogous locus. The purpose of this study was to determine if differences in phenotype (site of expression assayed using allele-specific reverse transcription-polymerase chain reaction) or function (dominance relationships assayed through controlled pollinations) accompany the sequence-based classification. Expression of Class A alleles was restricted to floral tissues, as predicted for genes involved in the SI response. In contrast, Class B alleles, despite being tightly linked to the SI phenotype, were unexpectedly expressed in both leaves and floral tissues; the same pattern found for a related unlinked paralogous sequence. Whereas Class A included haplotypes in three different dominance classes, all Class B haplotypes were found to be recessive to all except one Class A haplotype. In addition, mapping of expression and dominance patterns onto an S-domain-based genealogy suggested that allelic dominance may be determined more by evolutionary history than by frequency-dependent selection for lowered dominance as some theories suggest. The possibility that interlocus gene conversion might have contributed to allelic diversity is discussed.

Single-locus gametophytic incompatibility in autotetraploids

It is known that a single-locus gametophytic self-incompatibility (GSI) system can persist with just two distinct alleles in an autotetraploid population, in contrast to diploid GSI systems, assuming "competitive interaction" in which heteroallelic pollen is universally compatible. The steady-state population structure of a GSI system in autotetraploids was investigated in an undivided population assuming "competitive interaction." A deterministic model was developed to predict the frequencies of genotypes with two, three, or four distinct S alleles, assuming no mutation or population subdivision. The model showed that unlike in diploid GSI systems, the limiting values of the frequencies of genotype classes do not minimize pollen wastage.