THE DISTRIBUTION OF SELF-STERILITY ALLELES IN POPULATIONS

Selection at work in self-incompatible Arabidopsis lyrata: mating patterns in a natural population

Identification and characterization of the self-incompatibility genes in Brassicaceae species now allow typing of self-incompatibility haplotypes in natural populations. In this study we sampled and mapped all 88 individuals in a small population of Arabidopsis lyrata from Iceland. The self-incompatibility haplotypes at the SRK gene were typed for all the plants and some of their progeny and used to investigate the realized mating patterns in the population. The observed frequencies of haplotypes were found to change considerably from the parent generation to the offspring generation around their deterministic equilibria as determined from the known dominance relations among haplotypes. We provide direct evidence that the incompatibility system discriminates against matings among adjacent individuals. Multiple paternity is very common, causing mate availability among progeny of a single mother to be much larger than expected for single paternity.

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.

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.

Maximum-likelihood estimation of rates of recombination within mating-type regions

Features common to many mating-type regions include recombination suppression over large genomic tracts and cosegregation of genes of various functions, not necessarily related to reproduction. Model systems for homomorphic self-incompatibility (SI) in flowering plants share these characteristics. We introduce a method for the exact computation of the joint probability of numbers of neutral mutations segregating at the determinant of mating type and at a linked marker locus. The underlying Markov model incorporates strong balancing selection into a two-locus coalescent. We apply the method to obtain a maximum-likelihood estimate of the rate of recombination between a marker locus, 48A, and S-RNase, the determinant of SI specificity in pistils of Nicotiana alata. Even though the sampled haplotypes show complete allelic linkage disequilibrium and recombinants have never been detected, a highly significant deficiency of synonymous substitutions at 48A compared to S-RNase suggests a history of recombination. Our maximum-likelihood estimate indicates a rate of recombination of perhaps 3 orders of magnitude greater than the rate of synonymous mutation. This approach may facilitate the construction of genetic maps of regions tightly linked to targets of strong balancing selection.

Evolution under tight linkage to mating type

Recent large-scale sequencing studies of mating type loci in a number of organisms offer insight into the origin and evolution of these genomic regions. Extensive tracts containing genes with a wide diversity of functions typically cosegregate with mating type. Cases in which mating type determination entails complementarity between distinct transcription units may descend from systems in which close physical linkage facilitated the coordinated expression and cosegregation of the interacting genes. In response to the particular selection pressures associated with the maintenance of more than one mating type, this nucleus of low recombination may expand over evolutionary time, engulfing neighboring tracts bearing genes with no direct role in reproduction. This scenario is consistent with the present-day structure of some mating type loci, including regulators of homomorphic self-incompatibility in angiosperms (S-loci). Recombination suppression and enforced S-locus heterozygosity accelerate the accumulation of genetic load and promote genetic associations between S-alleles and degenerating genes in cosegregating tracts. This S-allele-specific load may influence the evolution of self-incompatibility systems.

ON THE EVOLUTIONARY COSTS OF SELF-INCOMPATIBILITY: INCOMPLETE REPRODUCTIVE COMPENSATION DUE TO POLLEN LIMITATION

Pollen limitation affects plants with diverse reproductive systems and ecologies. In self-incompatible (SI) species, pollen limitation may preclude full reproductive compensation for prezygotic rejection of pollen. We present a model designed to explore the effects of incomplete reproductive compensation on evolutionary changes at a modifier locus that regulates the level of SI expression. Our results indicate that incomplete reproductive compensation greatly increases the evolutionary costs of SI, particularly in populations with low S-allele diversity. The evolutionary fate of modifiers of SI expression depends on the rate at which they are transmitted to future generations as well as the effects of SI on offspring number and quality. Partial SI expression can represent a stable condition rather than an evolutionarily transient state between full expression and full suppression. This unanticipated result provides the first theoretical support for the evolutionary stability of such mixed mating systems, the existence of which has recently been documented.

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.

Uneven segregation of sporophytic self-incompatibility alleles in Arabidopsis lyrata

Self-incompatibility in Arabidopsis lyrata is sporophytically controlled by the multi-allelic S-locus. Self-incompatibility alleles (S-alleles) are under strong negative frequency dependent selection because pollen carrying common S-alleles have fewer mating opportunities. Population genetics theory predicts that deleterious alleles can accumulate if linked to the S-locus. This was tested by studying segregation of S-alleles in 11 large full sib families in A. lyrata. Significant segregation distortion leading to an up to fourfold difference in transmission rates was found in six families. Differences in transmission rates were not significantly different in reciprocal crosses and the distortions observed were compatible with selection acting at the gametic stage alone. The S-allele with the largest segregation advantage is also the most recessive, and is very common in natural populations concordant with its apparent segregation advantage. These results imply that frequencies of S-alleles in populations of A. lyrata cannot be predicted based on simple models of frequency-dependent selection alone.

A simple method for computing exact probabilities of mutation numbers

We describe a method for the recursive computation of exact probability distributions for the number of neutral mutations segregating in samples of arbitrary size and configuration. Construction of the recursions requires only characterization of evolutionary changes as a Markov process and determination of one-step transition matrices. We address the pattern of nucleotide diversity at a neutral marker locus linked to a determinant of mating type. Under a reformulation of parameters, the method also applies directly to metapopulation models with island migration among demes. Characterization of complete probability distributions facilitates parameter estimation and hypothesis testing by likelihood- as well as moment-based approaches.

EVOLUTION OF HIERARCHICAL CYTOPLASMIC INHERITANCE IN THE PLASMODIAL SLIME MOLD PHYSARUM POLYCEPHALUM

A striking linear dominance relationship for uniparental mitochondrial transmission is known between many mating types of plasmodial slime mold Physarum polycephalum. We herein examine how such hierarchical cytoplasmic inheritance evolves in isogamous organisms with many self-incompatible mating types. We assume that a nuclear locus determines the mating type of gametes and that another nuclear locus controls the digestion of mitochondria DNAs (mtDNAs) of the recipient gamete after fusion. We then examine the coupled genetic dynamics for the evolution of self-incompatible mating types and biased mitochondrial transmission between them. In Physarum, a multiallelic nuclear locus matA controls both the mating type of the gametes and the selective elimination of the mtDNA in the zygotes. We theoretically examine two potential mechanisms that might be responsible for the preferential digestion of mitochondria in the zygote. In the first model, the preferential digestion of mitochondria is assumed to be the outcome of differential expression levels of a suppressor gene carried by each gamete (suppression-power model). In the second model (site-specific nuclease model), the digestion of mtDNAs is assumed to be due to their cleavage by a site-specific nuclease that cuts the mtDNA at unmethylated recognition sites. Also assumed is that the mtDNAs are methylated at the same recognition site prior to the fusion, thereby being protected against the nuclease of the same gamete, and that the suppressor alleles convey information for the recognition sequences of nuclease and methylase. In both models, we found that a linear dominance hierarchy evolves as a consequence of the buildup of a strong linkage disequilibrium between the mating-type locus and the suppressor locus, though it fails to evolve if the recombination rate between the two loci is larger than a threshold. This threshold recombination rate depends on the number of mating types and the degree of fitness reduction in the heteroplasmic zygotes. If the recombination rate is above the threshold, suppressor alleles are equally distributed in each mating type at evolutionary equilibrium. Based on the theoretical results of the site-specific nuclease model, we propose that a nested subsequence structure in the recognition sequence should underlie the linear dominance hierarchy of mitochondrial transmission.

Signatures of selection among sex-determining alleles of the honey bee

Patterns of DNA polymorphisms are a primary tool for dissecting signatures of selection; however, the underlying selective forces are poorly understood for most genes. A classical example of diversifying selection is the complementary sex-determining locus that is found in the very large insect order Hymenoptera (bees, wasps, ants, and sawflies). The gene responsible for sex determination, the complementary sex determiner (csd), has been most recently identified in the honey bee. Females are heterozygous at this locus. Males result when there is only one functional allele present, as a result of either homozygosity (fertilized eggs) or, more commonly, hemizygosity (unfertilized eggs). The homozygotes, diploid males, do not reproduce and have zero fitness, which implies positive selection in favor of rare alleles. Large differences in csd cDNA sequences within and between four populations were found that fall into two major groups, types I and II. Type I consists of several allelic lineages that were maintained over an extended period, an indication of balancing selection. Diversifying selection has operated on several confined parts of the protein, as shown by an excess of nonsynonymous differences. Elevated sequence differences indicate another selected part near a repeat region. These findings have general implications about the understanding of both the function of the multiallelic mechanism and the adaptive processes on the level of nucleotide sequences. Moreover, the first csd sequence data are a notable basis for the avoidance of diploid males in bee selection programs by allele-assisted breeding.

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.

Parasitic exploitation as an engine of diversity

Parasitic exploitation occurs within and between a wide variety of taxa in a plethora of diverse contexts. Theoretical and empirical analyses indicate that parasitic exploitation can generate substantial genetic and phenotypic polymorphism within species. Under some circumstances, parasitic exploitation may also be an important factor causing reproductive isolation and promoting speciation. Here we review research relevant to the relationship between parasitic exploitation, within species-polymorphism, and speciation in some of the major arenas in which such exploitation has been studied. This includes research on the vertebrate major histocompatibility loci, plant-pathogen interactions, the evolution of sexual reproduction, intragenomic conflict, sexual conflict, kin mimicry and social parasitism, tropical forest diversity and the evolution of language. We conclude by discussing some of the issues raised by comparing the effect of parasitic exploitation on polymorphism and speciation in different contexts.

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.

Diversity and linkage of genes in the self-incompatibility gene family in Arabidopsis lyrata

We report studies of seven members of the S-domain gene family in Arabidopsis lyrata, a member of the Brassicaceae that has a sporophytic self-incompatibility (SI) system. Orthologs for five loci are identifiable in the self-compatible relative A. thaliana. Like the Brassica stigmatic incompatibility protein locus (SRK), some of these genes have kinase domains. We show that several of these genes are unlinked to the putative A. lyrata SRK, Aly13. These genes have much lower nonsynonymous and synonymous polymorphism than Aly13 in the S-domains within natural populations, and differentiation between populations is higher, consistent with balancing selection at the Aly13 locus. One gene (Aly8) is linked to Aly13 and has high diversity. No departures from neutrality were detected for any of the loci. Comparing different loci within A. lyrata, sites corresponding to hypervariable regions in the Brassica S-loci (SLG and SRK) and in comparable regions of Aly13 have greater replacement site divergence than the rest of the S-domain. This suggests that the high polymorphism in these regions of incompatibility loci is due to balancing selection acting on sites within or near these regions, combined with low selective constraints.

Effects of inbreeding on the genetic diversity of populations

The study of variability within species is important to all biologists who use genetic markers. Since the discovery of molecular variability among normal individuals, data have been collected from a wide range of organisms, and it is important to understand the major factors affecting diversity levels and patterns. Comparisons of inbreeding and outcrossing populations can contribute to this understanding, and therefore studying plant populations is important, because related species often have different breeding systems. DNA sequence data are now starting to become available from suitable plant and animal populations, to measure and compare variability levels and test predictions.

The population genetics of sporophytic self-incompatibility in Senecio squalidus L. (Asteraceae): avoidance of mating constraints imposed by low S-allele number

Senecio squalidus L. (Asteraceae) has been the subject of several ecological and population genetic studies due to its well-documented history of introduction, establishment and spread throughout Britain in the past 300 years. Our recent studies have focused on identifying and quantifying factors associated with the sporophytic self-incompatibility (SSI) system of S. squalidus that may have contributed to its success as a colonist. These findings are of general biological interest because they provide important insights into the short-term evolutionary dynamics of a plant mating system. The number of S-alleles in populations and their dominance interactions were investigated in eight wild British populations using cross-diallel studies. The numbers of S-alleles in British S. squalidus populations are typically low (average of 5.3 S-alleles) and the entire British population is estimated to possess no more than 7-11 S-alleles. Such low numbers of S-alleles are most probably a consequence of population bottlenecks associated with introduction and colonization. Potential evolutionary impacts on SSI caused by a paucity of S-alleles, such as restricted mate availability, are discussed, and we suggest that increased dominance interactions between S-alleles may be an important short-term means of increasing mate availability when S-allele numbers are low.

Genealogy-dependent variation in viability among self-incompatibility genotypes

Many hermaphroditic plants avoid self-fertilization by rejecting pollen that express genetically determined specificities in common with the pistil. The S-locus, comprising the determinants of pistil and pollen specificity, typically shows extremely high polymorphism, with dozens to hundreds of specificities maintained within species. This article explores a conjecture, motivated by empirical findings, that the expression of recessive deleterious factors at sites closely linked to the S-locus may cause greater declines in the viability of zygotes constituted from more closely related S-alleles. Diffusion approximation models incorporating variation in viability among S-locus genotypes and antagonistic interactions between a new specificity and its immediate parent specificity are constructed and analyzed. Results indicate that variation in viability tends to reduce the number of specificities maintained in a population at stochastic steady state, and that genealogy-based antagonism reduces the rate of bifurcation of S-allele lineages. These effects may account for some of the unusual features observed in empirical studies of S-allele genealogies.

Estimating the number, frequency, and dominance of S-alleles in a natural population of Arabidopsis lyrata(Brassicaceae) with sporophytic control of self-incompatibility

In homomorphic plant self-incompatibility (SI) systems, large numbers of alleles may be maintained at a single Mendelian locus. Most estimators of the number of alleles present in natural populations are designed for gametophytic self-incompatibility systems (GSI) in which the recognition phenotype of the pollen is determined by its own haploid genotype. In sporophytic systems (SSI), the recognition phenotype of the pollen is determined by the diploid genotype of its parent, and dominance differs among alleles. We describe research aimed at estimates of S-allele numbers in a natural population of Arabidopsis lyrata (Brassicaceae), whose SSI system has recently been described. Using a combination of pollination studies and PCR-based identification of alleles at a locus equivalent to the Brassica SRK gene, we identified and sequenced 11 putative alleles in a sample of 20 individuals from different maternal seed sets. The pollination results indicate that we have not amplified all alleles that must be present. Extensive partial incompatibility, nonreciprocal compatibility differences, and evidence of weakened expression of SI in some genotypes, prevent us from determining the exact number of missing alleles based only on cross-pollination data. Although we show that none of the theoretical models currently proposed is completely appropriate for estimating the number of alleles in this system, we estimate that there are between 13 and 16 different S-alleles in our sample, probably between 16 and 25 alleles in the population, and discuss the relative frequency of alleles in relation to dominance.

The population genetics of sporophytic self-incompatibility in Senecio squalidus L. (Asteraceae) I: S allele diversity in a natural population

Twenty-six individuals of the sporophytic self-incompatible (SSI) weed, Senecio squalidus were crossed in a full diallel to determine the number and frequency of S alleles in an Oxford population. Incompatibility phenotypes were determined by fruit-set results and the mating patterns observed fitted a SSI model that allowed us to identify six S alleles. Standard population S allele number estimators were modified to deal with S allele data from a species with SSI. These modified estimators predicted a total number of approximately six S alleles for the entire Oxford population of S. squalidus. This estimate of S allele number is low compared to other estimates of S allele diversity in species with SSI. Low S allele diversity in S. squalidus is expected to have arisen as a consequence of a disturbed population history since its introduction and subsequent colonisation of the British Isles. Other features of the SSI system in S. squalidus were also investigated: (a) the strength of self-incompatibility response; (b) the nature of S allele dominance interactions; and (c) the relative frequencies of S phenotypes. These are discussed in view of the low S allele diversity estimates and the known population history of S. squalidus.

Multilocus self-recognition systems in fungi as a cause of trans-species polymorphism

Trans-species polymorphism, meaning the presence of alleles in different species that are more similar to each other than they are to alleles in the same species, has been found at loci associated with vegetative incompatibility in filamentous fungi. If individuals differ at one or more of these loci (termed het for heterokaryon), they cannot form stable heterokaryons after vegetative fusion. At the het-c locus in Neurospora crassa and related species there is clear evidence of trans-species polymorphism: three alleles have persisted for approximately 30 million years. We analyze a population genetic model of multilocus vegetative incompatibility and find the conditions under which trans-species polymorphism will occur. In the model, several unlinked loci determine the vegetative compatibility group (VCG) of an individual. Individuals of different VCGs fail to form productive heterokaryons, while those of the same VCG form viable heterokaryons. However, viable heterokaryon formation between individuals of the same VCG results in a loss in fitness, presumably via transfer of infectious agents by hyphal fusion or exploitation by aggressive genotypes. The result is a form of balancing selection on all loci affecting an individual's VCG. We analyze this model by making use of a Markov chain/strong selection, weak mutation (SSWM) approximation. We find that trans-species polymorphism of the type that has been found at the het-c locus is expected to occur only when the appearance of new incompatibility alleles is strongly constrained, because the rate of mutation to such alleles is very low, because the number of possible incompatibility alleles at each locus is restricted, or because the number of incompatibility loci is limited.

S-allele diversity in Sorbus aucuparia and Crataegus monogyna (Rosaceae: Maloideae)

RT-PCR was used to obtain the first estimates from natural populations of allelic diversity at the RNase-based gametophytic self-incompatibility locus in the Rosaceae. A total of 20 alleles were retrieved from 20 Sorbus aucuparia individuals, whereas 17 alleles were found in 13 Crataegus monogyna samples. Estimates of population-level allele numbers fall within the range observed in the Solanaceae, the only other family with RNase-based incompatibility for which estimates are available. The nucleotide diversity of S-allele sequences was found to be much lower in the two Rosaceae species as compared with the Solanaceae. This was not due to a lower sequence divergence among most closely related alleles. Rather, it is the depth of the entire genealogy that differs markedly in the two families, with Rosaceae S-alleles exhibiting more recent apparent coalescence. We also investigated patterns of selection at the molecular level by comparing nucleotide diversity at synonymous and nonsynonymous sites. Stabilizing selection was inferred for the 5' region of the molecule, while evidence of diversifying selection was present elsewhere.

Molecular mechanisms underlying the breakdown of gametophytic self-incompatibility

The breakdown of self-incompatibility has occurred repeatedly throughout the evolution of flowering plants and has profound impacts on the genetic structure of populations. Recent advances in understanding of the molecular basis of self-incompatibility have provided insights into the mechanisms of its loss in natural populations, especially in the tomato family, the Solanaceae. In the Solanaceae, the gene that controls self-incompatibility in the style codes for a ribonuclease that causes the degradation of RNA in pollen tubes bearing an allele at the S-locus that matches either of the two alleles held by the maternal plant. The pollen component of the S-locus has yet to be identified. Loss of self-incompatibility can be attributed to three types of causes: duplication of the S-locus, mutations that cause loss of S-RNase activity, and mutations that do not cause loss of S-RNase activity. Duplication of the S-locus has been well studied in radiation-induced mutants but may be a relatively rare cause of the breakdown of self-incompatibility in nature. Point mutations within the S-locus that disrupt the production of S-RNase have been documented in natural populations. There are also a number of mutants in which S-RNase production is unimpaired, yet self-incompatibility is disrupted. The identity and function of these mutations is not well understood. Careful work on a handful of model organisms will enable population biologists to better understand the breakdown of self-incompatibility in nature.

Molecular variation at the self-incompatibility locus in natural populations of the genera Antirrhinum and Misopates

The self-incompatibility system of flowering plants is a classic example of extreme allelic polymorphism maintained by frequency-dependent selection. We used primers designed from three published Antirrhinum hispanicum S-allele sequences in PCR reactions with genomic DNA of plants sampled from natural populations of Antirrhinum and Misopates species. Not surprisingly, given the polymorphism of S-alleles, only a minority of individuals yielded PCR products of the expected size. These yielded 35 genomic sequences, of nine different sequence types of which eight are highly similar to the A. hispanicum S-allele sequences, and one to a very similar unpublished Antirrhinum S-like RNase sequence. The sequence types are well separated from the S-RNase sequences from Solanaceae and Rosaceae, and also from most known "S-like" RNase sequences (which encode proteins not involved in self-incompatibility). An association with incompatibility types has so far been established for only one of the putative S-alleles, but we describe evidence that the other sequences are also S-alleles. Variability in these sequences follows the pattern of conserved and hypervariable regions seen in other S-RNases, but no regions have higher replacement than silent diversity, unlike the results in some other species.

Inbreeding depression in small populations of self-incompatible plants

Self-incompatibility (SI) is a widespread mechanism that prevents inbreeding in flowering plants. In many species, SI is controlled by a single locus (the S locus) where numerous alleles are maintained by negative frequency-dependent selection. Inbreeding depression, the decline in fitness of selfed individuals compared to outcrossed ones, is an essential factor in the evolution of SI systems. Conversely, breeding systems influence levels of inbreeding depression. Little is known about the joint effect of SI and drift on inbreeding depression. Here we studied, using a two-locus model, the effect of SI (frequency-dependent selection) on a locus subject to recurrent deleterious mutations causing inbreeding depression. Simulations were performed to assess the effect of population size and linkage between the two loci on the level of inbreeding depression and genetic load. We show that the sheltering of deleterious alleles linked to the S locus strengthens inbreeding depression in small populations. We discuss the implications of our results for the evolution of SI systems.

Consequences of population structure on genes under balancing selection

This paper describes a new approach to modeling population structure for genes under strong balancing selection of the type seen in plant self-incompatibility systems and the major histocompatibility complex (MHC) system of vertebrates. Simple analytic solutions for the number of alleles maintained at equilibrium and the expected proportion of alleles shared between demes at various levels are derived and checked against simulation results. The theory accurately captures the dynamics of allele number in a subdivided population and identifies important values of m (migration rate) at which allele number and distribution change qualitatively. Starting from a panmictic population, as migration among demes decreases a qualitative change in dynamics is seen at approximately m(crit) approximately equal to the square root of(s/4piNT) where NT is the total population size and s is a measure of the strength of selection. At this point, demes can no longer maintain their panmictic allele number, due to increasing isolation from the total population. Another qualitative change occurs at a migration rate on the same order of magnitude as the mutation rate, mu. At this point, the demes are highly differentiated for allele complement, and the total number of alleles in the population is increased. Because in general u << m<(crit) at intermediate migration rates slightly fewer alleles may be maintained in the total population than are maintained at panmixia. Within this range, total allele number may not be the best indicator of whether a population is effectively panmictic, and some caution should be used when interpreting samples from such populations. The theory presented here can help to analyze data from genes under balancing selection in subdivided populations.

On the origin of self-incompatibility haplotypes: transition through self-compatible intermediates

Self-incompatibility (SI) in flowering plants entails the inhibition of fertilization by pollen that express specificities in common with the pistil. In species of the Solanaceae, Rosaceae, and Scrophulariaceae, the inhibiting factor is an extracellular ribonuclease (S-RNase) secreted by stylar tissue. A distinct but as yet unknown gene (provisionally called pollen-S) appears to determine the specific S-RNase from which a pollen tube accepts inhibition. The S-RNase gene and pollen-S segregate with the classically defined S-locus. The origin of a new specificity appears to require, at minimum, mutations in both genes. We explore the conditions under which new specificities may arise from an intermediate state of loss of self-recognition. Our evolutionary analysis of mutations that affect either pistil or pollen specificity indicates that natural selection favors mutations in pollen-S that reduce the set of pistils from which the pollen accepts inhibition and disfavors mutations in the S-RNase gene that cause the nonreciprocal acceptance of pollen specificities. We describe the range of parameters (rate of receipt of self-pollen and relative viability of inbred offspring) that permits the generation of a succession of new specificities. This evolutionary pathway begins with the partial breakdown of SI upon the appearance of a mutation in pollen-S that frees pollen from inhibition by any S-RNase presently in the population and ends with the restoration of SI by a mutation in the S-RNase gene that enables pistils to reject the new pollen type.

Roles of lineage sorting and phylogenetic relationship in the genetic diversity at the self-incompatibility locus of Solanaceae

Allelic polymorphism at the S locus that determines the gametophytic self-incompatibility (GSI) system in the pistil predates speciation. Understanding the evolution of a GSI system therefore requires knowledge of how lineage sorting and interspecific phylogenetic relationship affect S allele polymorphism. In searching for patterns of lineage sorting among species of various phylogenetic relationships, 22 S-alleles from 34 genets randomly taken at three Tennessee sites from a newly known GSI species Physalis longifolia were sequenced. Analyses of these data along with the previous sequences of three solanaceous species indicate that much of the combined allelic genealogy may be explained by lineage sorting and phylogenetic relationship. Using the mean terminal branch lengths of trans-specific alleles on the allelic genealogy to infer phylogenetic relationship among species, P. longifolia was found to be more closely related to P. cinerascens than to P. crassifolia. Nonetheless, the distribution of terminal branch lengths of P. longifolia was more similar to that of P. crassifolia than to that of P. cinerascens, suggesting phylogenetic relationship may have little effect on species-specific polymorphism. Similar habitat and growth characters, yet contrasting S-polymorphism, between P. longifolia and P. cinerascens also reject previous hypotheses that habitat and growth characters are the major factors responsible for interspecific differences in S-polymorphism. A likely scenario is that species-specific S-polymorphism is based on lineage sorting whose effect is further modified by species age and historical changes in population parameters.

Polymorphism at the ribosomal DNA spacers and its relation to breeding structure of the widespread mushroom Schizophyllum commune

The common split-gilled mushroom Schizophyllum commune is found throughout the world on woody substrates. This study addresses the dispersal and population structure of this fungal species by studying the phylogeny and evolutionary dynamics of ribosomal DNA (rDNA) spacer regions. Extensive sampling (n = 195) of sequences of the intergenic spacer region (IGS1) revealed a large number of unique haplotypes (n = 143). The phylogeny of these IGS1 sequences revealed strong geographic patterns and supported three evolutionarily distinct lineages within the global population. The same three geographic lineages were found in phylogenetic analysis of both other rDNA spacer regions (IGS2 and ITS). However, nested clade analysis of the IGS1 phylogeny suggested the population structure of S. commune has undergone recent changes, such as a long distance colonization of western North America from Europe as well as a recent range expansion in the Caribbean. Among all spacer regions, variation in length and nucleotide sequence was observed between but not within the tandem rDNA repeats (arrays). This pattern is consistent with strong within-array and weak among-array homogenizing forces. We present evidence for the suppression of recombination between rDNA arrays on homologous chromosomes that may account for this pattern of concerted evolution.

Mutational origin of new mating type specificities in flowering plants

Many hermaphroditic plants avoid self-fertilization by rejecting pollen that express genetically-determined specificities in common with the pistil. Self-incompatibility systems typically show extremely high genetic diversity, some maintaining hundreds of specificities. This article addresses the genetic and evolutionary mechanisms through which new mating specificities arise. Recent investigations of the genetic and physiological basis of self-incompatibility are reviewed. Two evolutionary pathways are considered: one which requires full expression of self-incompatibility in all intermediates and one in which new mating specificities arise through episodes of partial breakdown and restoration of self-incompatibility.

A model for the evolution of high frequencies of males in an androdioecious plant based on a cross-compatibility advantage of males

Lloyd's (1975) and Charlesworth & Charlesworth's (1978) phenotypic selection models for the maintenance of androdioecy predict that males (female-sterile individuals) must have an advantage in fertility (K) of at least two in order to invade a hermaphroditic population, and that their equilibrium frequency (x(eq)=(K - 2)/2(K - 1)) is always less than 0.5. In this paper, we develop a model in which male fertility is frequency-dependent, a situation not investigated in the previous models, to explore the conditions under which a high frequency of males (i.e. more than 50%) could be maintained at equilibrium. We demonstrate that a gametophytic self-incompatibility (GSI) locus linked to a nuclear sex determination locus can favour rare alleles through male function, by causing frequency-dependent selection. Thus, the spread of a female-sterility allele in a hermaphroditic population may be induced. In contrast with the previous models, our model can explain male frequencies greater than 50% in a functionally androdioecious species, as long as there is (i) dominance of female-sterility at the sex locus, and (ii) a few alleles at the self-incompatibility locus, even if the advantage in fertility of male phenotype is lower than two.

Evolutionary dynamics of self-incompatibility alleles in Brassica

Self-incompatibility in Brassica entails the rejection of pollen grains that express specificities held in common with the seed parent. In Brassica, pollen specificity is encoded at the multipartite S-locus, a complex region comprising many expressed genes. A number of species within the Brassicaceae express sporophytic self-incompatibility, under which individual pollen grains bear specificities determined by one or both S-haplotypes of the pollen parent. Classical genetic and nucleotide-level analyses of the S-locus have revealed a dichotomy in sequence and function among S-haplotypes; in particular, all class I haplotypes show dominance over all class II haplotypes in determination of pollen specificity. Analysis of an evolutionary model that explicitly incorporates features of the Brassica system, including the class dichotomy, indicates that class II haplotypes may invade populations at lower rates and decline to extinction at higher rates than class I haplotypes. This analysis suggests convergence to an evolutionarily persistent state characterized by the maintenance in high frequency of a single class II haplotype together with many class I haplotypes, each in low frequency. This expectation appears to be consistent with empirical observations of high frequencies of relatively few distinct recessive haplotypes.

The ancestral graph and gene genealogy under frequency-dependent selection

Minority-advantage frequency-dependent selection has been proposed as the cause for the high level of observed polymorphism in some self/nonself-recognition systems. We present a mathematically rigorous derivation of the ancestral graph for a sample of genes that evolved according to a haploid infinite-alleles model of minority-advantage frequency-dependent selection. In the case of sufficiently weak selection, the gene genealogy can be extracted from the ancestral graph. We demonstrate that the gene genealogy under this model is identical to that obtained for a diploid model with heterozygote advantage. The case of strong selection is exemplified by a one-locus haploid self-incompatibility system; in this context, we investigate the number of alleles that can be maintained in a spatial versus a non-spatial habitat. Finally, we compare gametophytic self-incompatibility to the haploid self-incompatibility model.

Cellular and molecular mechanisms of sexual incompatibility in plants and fungi

Plants and fungi show an astonishing diversity of mechanisms to promote outbreeding, the most widespread of which is sexual incompatibility. Sexual incompatibility involves molecular recognition between mating partners. In fungi and algae, highly polymorphic mating-type loci mediate mating through complementary interactions between molecules encoded or regulated by different mating-type haplotypes, whereas in flowering plants polymorphic self-incompatibility loci regulate mate recognition through oppositional interactions between molecules encoded by the same self-incompatibility haplotypes. This subtle mechanistic difference is a consequence of the different life cycles of fungi, algae, and flowering plants. Recent molecular and biochemical studies have provided fascinating insights into the mechanisms of mate recognition and are beginning to shed light on evolution and population genetics of these extraordinarily polymorphic genetic systems of incompatibility.

Overdominant alleles in a population of variable size

An approximate method is developed to predict the number of strongly overdominant alleles in a population of which the size varies with time. The approximation relies on the strong-selection weak-mutation (SSWM) method introduced by J. H. Gillespie and leads to a Markov chain model that describes the number of common alleles in the population. The parameters of the transition matrix of the Markov chain depend in a simple way on the population size. For a population of constant size, the Markov chain leads to results that are nearly the same as those of N. Takahata. The Markov chain allows the prediction of the numbers of common alleles during and after a population bottleneck and the numbers of alleles surviving from before a bottleneck. This method is also adapted to modeling the case in which there are two classes of alleles, with one class causing a reduction in fitness relative to the other class. Very slight selection against one class can strongly affect the relative frequencies of the two classes and the relative ages of alleles in each class.

Identification of regions in which positive selection may operate in S-RNase of Rosaceae: implication for S-allele-specific recognition sites in S-RNase

A stylar S-RNase is associated with gametophytic self-incompatibility in the Rosaceae, Solanaceae, and Scrophulariaceae. This S-RNase is responsible for S-allele-specific recognition in the self-incompatible reaction, but how it functions in specific discrimination is not clear. Window analysis of the numbers of synonymous (dS) and non-synonymous (dN) substitutions in rosaceous S-RNases detected four regions with an excess of dN over dS in which positive selection may operate (PS regions). The topology of the secondary structure of the S-RNases predicted by the PHD method is very similar to that of fungal RNase Rh whose tertiary structure is known. When the sequences of S-RNases are aligned with the sequence of RNase Rh based on the predicted secondary structures, the four PS regions correspond to two surface sites on the tertiary structure of RNase Rh. These findings suggest that in S-RNases the PS regions also form two sites and are candidates for the recognition sites for S-allele-specific discrimination.

Allelic genealogies in sporophytic self-incompatibility systems in plants

Expectations for the time scale and structure of allelic genealogies in finite populations are formed under three models of sporophytic self-incompatibility. The models differ in the dominance interactions among the alleles that determine the self-incompatibility phenotype: In the SSIcod model, alleles act codominantly in both pollen and style, in the SSIdom model, alleles form a dominance hierarchy, and in SSIdomcod, alleles are codominant in the style and show a dominance hierarchy in the pollen. Coalescence times of alleles rarely differ more than threefold from those under gametophytic self-incompatibility, and transspecific polymorphism is therefore expected to be equally common. The previously reported directional turnover process of alleles in the SSIdomcod model results in coalescence times lower and substitution rates higher than those in the other models. The SSIdom model assumes strong asymmetries in allelic action, and the most recessive extant allele is likely to be the most recent common ancestor. Despite these asymmetries, the expected shape of the allele genealogies does not deviate markedly from the shape of a neutral gene genealogy. The application of the results to sequence surveys of alleles, including interspecific comparisons, is discussed.

Flowering plant self-incompatibility: the molecular population genetics of Brassica S-loci

Self-incompatibility systems in different angiosperm families are reviewed, and the evidence that incompatibility has arisen several times is outlined. New data on the sequence polymorphism of self-incompatibility loci from two different angiosperm families are compared with results from other highly polymorphic loci, particularly MHC loci. We discuss what molecular genetic analyses of these sequences can tell us about the nature and maintenance of the polymorphism at self-incompatibility loci. We suggest that there is evidence for recombination at the Brassica self-incompatibility loci, so that it may be possible to discern regions that are particularly functionally important in the recognition reaction, even though the long-term maintenance of polymorphisms in these amino acid residues has caused the evolution of many other sequence differences between alleles.

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

The actual and effective number of gametophytic self-incompatibility alleles maintained at mutation-drift-selection equilibrium in a finite population subdivided as in the island model is investigated by stochastic simulations. The existing theory founded by Wright predicts that for a given population size the number of alleles maintained increases monotonically with decreasing migration as is the case for neutral alleles. The simulation results here show that this is not true. At migration rates above Nm = 0.01-0.1, the actual and effective number of alleles is lower than for an undivided population with the same number of individuals, and, contrary to Wright's theoretical expectation, the number of alleles is not much higher than for an undivided population unless Nm < 0.001. The same pattern is observed in a model where the alleles display symmetrical overdominant selection. This broadens the applicability of the results to include proposed models for the major histocompatibility (MHC) loci. For a subdivided population over a large range of migration rates, it appears that the number of self-incompatibility alleles (or MHC-alleles) observed can provide a rough estimate of the total number of individuals in the population but it underestimates the neutral effective size of the subdivided population.

Variation at the b Mating Type Locus of Ustilago maydis

ABSTRACT Population level diversity at the Ustilago maydis b mating type locus was determined in samples from four Minnesota locations using a combination of plate mating techniques and a polymerase chain reaction (PCR)-based assay. The PCR method allows rapid identification of b types from samples of natural populations and utilizes the hypervariable regions of the b locus that determine mating type specificity. Results demonstrated high levels of b diversity within populations, with one population yielding 17 of the total 18 b types found in the study. Pairwise G(ST) values were in the range of 0.02 to 0.05, and common b mating types were found across broad geographic distances. These data demonstrated that very low levels of differentiation among U. maydis populations occur with respect to b locus variation. Consistent with frequency-dependent selection models, b types were represented at approximately equal frequencies within the entire Minnesota population. However, neutral evolutionary models for patterns of geographic distribution and variation at b cannot be entirely excluded. The importance to agricultural practices of understanding population genetic processes is discussed.

Evolutionary dynamics of sporophytic self-incompatibility alleles in plants

The stationary frequency distribution and allelic dynamics in finite populations are analyzed through stochastic simulations in three models of single-locus, multi-allelic sporophytic self-incompatibility. The models differ in the dominance relationships among alleles. In one model, alleles act codominantly in both pollen and style (SSIcod), in the second, alleles form a dominance hierarchy in pollen and style (SSIdom). In the third model, alleles interact codominantly in the style and form a dominance hierarchy in the pollen (SSIdomcod). The SSIcod model behaves similarly to the model of gametophytic self-incompatibility, but the selection intensity is stronger. With dominance, dominant alleles invade the population more easily than recessive alleles and have a lower frequency at equilibrium. In the SSIdom model, recessive alleles have both a higher allele frequency and higher expected life span. In the SSIdomcod model, however, loss due to drift occurs more easily for pollen-recessive than for pollen-dominant alleles, and therefore, dominant alleles have a higher expected life span than the more recessive alleles. The process of allelic turnover in the SSIdomcod and SSIdom models is closely approximated by a random walk on a dominance ladder. Implications of the results for experimental studies of sporophytic self-incompatibility in natural populations are discussed.