THE DISTRIBUTION OF SELF-STERILITY ALLELES IN POPULATIONS

Seeing selection in S allele sequences

New data on allelic sequence diversity in natural populations provide evidence for natural selection acting on the self-incompatibility loci of two plant species; there are interesting parallels with, and differences from, other polymorphic systems such as mammalian MHC loci.

Correlation of genetic and physical maps at the A mating-type locus of Coprinus cinereus

The A mating type locus of Coprinus cinereus is remarkable for its extreme diversity, with over 100 different alleles in natural populations. Classical genetic studies have demonstrated that this hypervariability arises in part from recombination between two subloci of A, alpha and beta, although more recent population genetic data have indicated a third segregating sublocus. In this study, we characterized the molecular basis by which recombination generates nonparental A mating types. We mapped the frequency and location of all recombination events in two crosses and correlated the genetic and physical maps of A. We found that all recombination events were located in 6 kb of noncoding DNA between the alpha and beta subloci and that the rate of recombination in this noncoding region matched that generally observed for this genome. No recombination within gene clusters or within coding regions was observed, and the two alpha and beta subloci described in genetic analyses correlated with the previously characterized alpha and beta gene clusters. We propose that pairs of genes constitute both the sex determining and the hereditary unit of A.

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

The self-incompatibility (S) locus of flowering plants offers an example of extreme polymorphism maintained by balancing selection. Estimates of recent and long-term effective population size (Ne) were determined for two solanaceous species by examination of S-allele diversity. Estimates of recent Ne in two solanaceous species differed by an order of magnitude, consistent with differences in the species' ecology. In one species, the evidence was consistent with historical population restriction despite a large recent Ne. In the other, no severe bottleneck was indicated over millions of years. Bottlenecks are integral to founder-event speciation, and loci that are subject to balancing selection can be used to evaluate the frequency of this mode of speciation.

Origin of allelic diversity in antirrhinum S locus RNases

In many plant species, self-incompatibility (SI) is genetically controlled by a single multiallelic S locus. Previous analysis of S alleles in the Solanaceae, in which S locus ribonucleases (S RNases) are responsible for stylar expression of SI, has demonstrated that allelic diversity predated speciation within this family. To understand how allelic diversity has evolved, we investigated the molecular basis of gametophytic SI in Antirrhinum, a member of the Scrophulariaceae, which is closely related to the Solanaceae. We have characterized three Antirrhinum cDNAs encoding polypeptides homologous to S RNases and shown that they are encoded by genes at the S locus. RNA in situ hybridization revealed that the Antirrhinum S RNase are primarily expressed in the stylar transmitting tissue. This expression is consistent with their proposed role in arresting the growth of self-pollen tubes. S alleles from the Scrophulariaceae form a separate group from those of the Solanaceae, indicating that new S alleles have been generated since these families separated (approximately 40 million years). We propose that the recruitment of an ancestral RNase gene into SI occurred during an early stage of angiosperm evolution and that, since that time, new alleles subsequently have arisen at a low rate.

S-allele sequence diversity in natural populations of Solanum carolinense (Horsenettle)

S-allele diversity in Solanum carolinense was surveyed in two natural populations, located in Tennessee and North Carolina, with a molecular assay to determine the genotype of individual plants. A total of 13 different S-alleles were identified and sequenced. There is high overlap between the two populations sampled, with 10 alleles shared in common, one allele found only in Tennessee, and two found only in North Carolina. The number of alleles in this species appears to be extremely low compared with other species with gametophytic self-incompatibility. Sequence comparisons show that most alleles are extremely different one from another in their primary sequence and a phylogenetic analysis indicates extensive trans-specific evolution of S-lineages. In addition, some alleles appear to be derived much more recently. The implications of these observations are discussed in the light of recent theoretical results on S-allele population diversity and persistence.

A generalized least-squares estimate for the origin of sporophytic self-incompatibility

Analysis of nucleotide sequences that regulate the expression of self-incompatibility in flowering plants affords a direct means of examining classical hypotheses for the origin and evolution of this major feature of mating systems. Departing from the classical view of monophyly of all forms of self-incompatibility, the current paradigm for the origin of self-incompatibility postulates multiple episodes of recruitment and modification of preexisting genes. In Brassica, the S locus, which regulates sporophytic self-incompatibility, shows homology to a multigene family present both in self-compatible congeners and in groups for which this form of self-incompatibility is atypical. A phylogenetic analysis of S-allele sequences together with homologous sequences that do not cosegregate with self-incompatibility permits dating the change of function that marked the origin of self-incompatibility. A generalized least-squares method is introduced that provides closed-form expressions for estimates and standard errors for function-specific divergence rates and times of divergence among sequences. This analysis suggests that the age of the sporophytic self-incompatibility system expressed in Brassica exceeds species divergence within the genus by four- to fivefold. The extraordinarily high levels of sequence diversity exhibited by S alleles appears to reflect their ancient derivation, with the alternative hypothesis of hypermutability rejected by the analysis.

Linkage analysis of sex determination in the honey bee (Apis mellifera)

A colony-level phenotype was used to map the major sex determination locus (designated X) in the honey bee (Apis mellifera). Individual queen bees (reproductive females) were mated to single drones (fertile males) by instrumental insemination. Haploid drone progeny of an F1 queen were each backcrossed to daughter queens from one of the parental lines. Ninety-eight of the resulting colonies containing back-cross progeny were evaluated for the trait 'low brood-viability' resulting from the production of diploid drones that were homozygous at X. DNA samples from the haploid drone fathers of these colonies were used individually in polymerase chain reactions (PCR) with 10-base primers. These reactions generated random amplified polymorphic DNA (RAPD) markers that were analyzed for cosegregation with the colony-level phenotype. One RAPD marker allele was shared by 22 of 25 drones that fathered low brood-viability colonies. The RAPD marker fragment was cloned and partially sequenced. Two primers were designed that define a sequence-tagged site (STS) for this locus. The primers amplified DNA marker fragments that cosegregated with the original RAPD marker. In order to more precisely estimate the linkage between X and the STS locus, another group of bees consisting of progeny from one of the low-brood viability colonies was used in segregation analysis. Four diploid drones and 181 of their diploid sisters (workers, nonfertile females) were tested for segregation of the RAPD and STS markers. The cosegregating RAPD and STS markers were codominant due to the occurrence of fragment-length alleles. The four diploid drones were homozygous for these markers but only three of the 181 workers were homozygotes (recombinants). Therefore the distance between X and the STS locus was estimated at 1.6 cM. An additional linked marker was found that was 6.6 cM from the STS locus.

Gene and allelic genealogies at a gametophytic self-incompatibility locus

The properties of gene and allelic genealogies at a gametophytic self-incompatibility locus in plants have been investigated analytically and checked against extensive numerical simulations. It is found that, as with overdominant loci, there are two genealogical processes with markedly different time scales. First, functionally distinct allelic lines diverge on an extremely long time scale which is inversely related to the mutation rate to new alleles. These alleles show a genealogical structure which is similar, after an appropriate rescaling of time, to that described by the coalescent process for genes at a neutral locus. Second, gene copies sampled within the same functional allelic line show genealogical relationships similar to neutral gene genealogies but on a much shorter time scale, which is on the same order of magnitude as the harmonic mean of the number of gene copies within an allelic line. These results are discussed in relation to data showing trans-specific polymorphisms for alleles at the gametophytic self-incompatibility locus in the Solanaceae. It is shown that population sizes on the order of 4 x 10(5) and a mutation rate per locus per generation as high as 10(-6) could account for estimated allelic divergence times in this family.

The S11 and S13 self incompatibility alleles in Solanum chacoense Bitt. are remarkably similar

A genomic clone of the S11 allele from the self-incompatibility locus (S locus) in Solanum chacoense Bitt. has been isolated by cross-hybridization to the S. chacoense S13 allele and sequenced. The sequence of the S11 allele contains all the features expected for S genes of the Solanaceae, and S11 expression, as assessed by northern blots and RNA-PCR, was similar to that of other S. chacoense S alleles. The S11 protein sequence shares 95% identity with the phenotypically distinct S13 protein of S. chacoense and is the gametophytic S allele with the highest similarity to an existing allele so far discovered. Only 10 amino acid changes differentiate the mature proteins from these two alleles, which sets a new lower limit to the number of changes that can produce an altered S allele specificity. The amino acid substitutions are not clustered, suggesting that an accumulation of random point mutations can generate S allele diversity. The S11 intron is unusual in that it could be translated in frame with the coding sequence, thus suggesting an additional mechanism for the generation of new S alleles.

Genetic mapping and protein product diversity of the self-incompatibility locus in wild tomato (Lycopersicon peruvianum)

Phenotypic diversity of self-incompatibility (S) alleles within nine natural populations of Lycopersicon peruvianum was investigated. Only 7 incompatible responses were observed of a total of 276 unique combinations tested, on the basis of controlled pollinations, indicating the large number of alleles that exist within these populations. Molecular weight polymorphism for specific major stylar proteins observed on SDS-PAGE was also evident in two of the populations examined. Five proteins were shown to map to the S locus and to be associated with different S alleles through controlled pollinations and segregation of the proteins. Two of these S related proteins had been described previously in terms of spatial and temporal expression consistent with their involvement in self-incompatibility (Mau et al., Planta 169, 184-191, 1986). A mapping population derived from a fully compatible cross was used to establish linkage of the S locus to two DNA markers, CD15 and TG184, that lie on chromosome 1. The order of the markers and estimates of map distances are given.

Molecular diversity at the self-incompatibility locus is a salient feature in natural populations of wild tomato (Lycopersicon peruvianum)

A cDNA encoding a stylar protein was cloned from flowers of self-incompatible wild tomato (Lycopersicon peruvianum). The corresponding gene was mapped to the S locus, which is responsible for self-incompatibility. The nucleotide sequence was determined for this allele, and compared to other S-related sequences in the Solanaceae. The S allele was used to probe DNA from 92 plants comprising 10 natural populations of Lycopersicon peruvianum. Hybridization was conducted under moderate and permissive stringencies in order to detect homologous sequences. Few alleles were detected, even under permissive conditions, underscoring the great sequence diversity at this locus. Those alleles that were detected are highly homologous. Sequences could not be detected in self-incompatible Nicotiana alata, self-compatible L. esculentum (cultivated tomato) or self-compatible L. hirsutum. However, hybridization to an individual of self-incompatible L. hirsutum revealed a closely related sequence that maps to the S locus in this reproductively isolated species. This supports the finding that S locus polymorphism predates speciation. The extraordinarily high degree of sequence diversity present in the gametophytic self-incompatibility system is discussed in the context of other highly divergent systems representing several kingdoms.

Polymorphism and balancing selection at major histocompatibility complex loci

Amino acid replacements in the peptide-binding region (PBR) of the functional major histocompatibility complex (Mhc) genes appear to be driven by balancing selection. Of the various types of balancing selection, we have examined a model equivalent to overdominance that confers heterozygote advantage. As discussed by A. Robertson, overdominance selection tends to maintain alleles that have more or less the same degree of heterozygote advantage. Because of this symmetry, the model makes various testable predictions about the genealogical relationships among different alleles and provides ways of analyzing DNA sequences of Mhc alleles. In this paper, we analyze DNA sequences of 85 alleles at the HLA-A, -B, -C, -DRB1 and -DQB1 loci with respect to the number of alleles and extent of nucleotide differences at the PBR, as well as at the synonymous (presumably neutral) sites. Theory suggests that the number of alleles that differ at the sites targeted by selection (presumably the nonsynonymous sites in the PBR) should be equal to the mean number of nucleotide substitutions among pairs of alleles. We also demonstrate that the nucleotide substitution rate at the targeted sites relative to that of neutral sites may be much larger than 1. The predictions of the presented model are in surprisingly good agreement with the actual data and thus provide means for inferring certain population parameters. For overdominance selection in a finite population at equilibrium, the product of selection intensity (s) against homozygotes and the effective population size (N) is estimated to be 350-3000, being largest at the B locus and smallest at the C locus. We argue that N is of the order of 10(5) and s is several percent at most, if the mutation rate per site per generation is 10(-8).

Excess nonsynonymous substitution of shared polymorphic sites among self-incompatibility alleles of Solanaceae

The function of the self-incompatibility locus (S locus) of many plant species dictates that natural selection will favor high levels of protein diversity. Pairwise sequence comparisons between S alleles from four species of Solanaceae reveal remarkably high sequence diversity and evidence for shared polymorphism. The level of amino acid constraint was found to be significantly heterogeneous among different regions of the gene, with some regions being highly constrained and others appearing to be virtually unconstrained. In some regions of the protein, there was an excess of nonsynonymous over synonymous substitution, consistent with the strong diversifying selection that must operate on this locus. These hypervariable regions are candidates for the sites that determine functional allelic identity. Simple contingency table tests show that sites that have polymorphism shared between species have more nonsynonymous substitution than polymorphic sites that do not exhibit shared polymorphism. This is consistent with the idea that adaptive evolution favoring amino acid replacement is occurring at sites with shared polymorphism. Tests of clustered polymorphism reveal that an unusually low rate of recombination must be occurring in this locus, allowing very ancient alleles to preserve their identity.

Polymorphism at the self-incompatibility locus in Solanaceae predates speciation

Sequences of 11 alleles of the gametophytic self-incompatibility locus (S locus) from three species of the Solanaceae family have recently been determined. Pairwise comparisons of these alleles reveal two unexpected observations: (i) amino acid sequence similarity can be as low as 40% within species and (ii) some interspecific similarities are higher than intraspecific similarities. The gene genealogy clearly illustrates this unusual pattern of relationships. The data suggest that some of the polymorphism at the S locus existed prior to the divergence of these species and has been maintained to the present. In support of this hypothesis, the number of shared polymorphic sites was found to exceed the number found in simulations with independent accumulation of mutations. Strictly neutral evolution is exceedingly unlikely to maintain the polymorphism for such a long time. The allele multiplicity and extreme age of the alleles is consistent with Wright's classic one-locus population genetic model of gametophytic self-incompatibility. Similarities between the plant S locus and the mammalian major histocompatibility complex are discussed.

Allelic genealogy under overdominant and frequency-dependent selection and polymorphism of major histocompatibility complex loci

To explain the long-term persistence of polymorphic alleles (trans-specific polymorphism) at the major histocompatibility complex (MHC) loci in rodents and primates, a computer simulation study was conducted about the coalescence time of different alleles sampled under various forms of selection. At the same time, average heterozygosity, the number of alleles in a sample, and the rate of codon substitution were examined to explain the mechanism of maintenance of polymorphism at the MHC loci. The results obtained are as follows. (1) The coalescence time for neutral alleles is too short to explain the trans-specific polymorphism at the MHC loci. (2) Under overdominant selection, the coalescence time can be tens of millions of years, depending on the parameter values used. The average heterozygosity and the number of alleles observed are also high enough to explain MHC polymorphism. (3) The pathogen adaptation model proposed by Snell is incapable of explaining MHC polymorphism, since the coalescence time for this model is too short and the expected heterozygosity and the expected number of alleles are too small. (4) From the mathematical point of view, the minority advantage model of frequency-dependent selection is capable of explaining a high degree of polymorphism and trans-specific polymorphism. (5) The molecular mimicry hypothesis also gives a sufficiently long coalescence time when the mutation rate is low in the host but very high in the parasite. However, the expected heterozygosity and the expected number of alleles tend to be too small. (6) Consideration of the molecular mechanism of the function of MHC molecules and other biological observations suggest that the most important factor for the maintenance of MHC polymorphism is overdominant selection. However, some experiments are necessary to distinguish between the overdominance and frequency-dependent selection hypotheses.

Numbers of sporophytic self-incompatibility alleles in populations of wild radish

To estimate the numbers of sporophytic S-alleles in two adjacent populations of wild radish, we performed 701 reciprocal crosses among 50 individuals. Each cross was replicated five times in each direction. Sixteen plants were fully intercompatible, indicating the presence of at least 32 S-alleles in the two populations. A minimum of 22 S-alleles occur in a single population. The frequency of incompatibility was significantly higher for within-population crosses (14.5%) than for between-population crosses (7.8%). This suggests that the two populations differ in the composition and frequency of alleles at the S-locus.

On the evolution of genetic incompatibility systems. V. Origin of sporophytic self-incompatibility in response to overdominance in viability

Conditions for the origin of partial sporophytic self-incompatibility (SSI) are obtained from two quantitative models, which differ with respect to the determination of offspring viability. Offspring viability depends solely on the source (self or nonself) of the fertilizing pollen in the first model, which describes changes only at a primitive S-locus itself. Two loci evolve in the second model: overdominant viability selection maintains an arbitrary number of alleles at one locus, with SSI under the control of a separate locus. In both cases, the origin of SSI requires that the relative change in the numbers of offspring derived by the two reproductive modes compensate for the twofold cost of outcrossing. In the first model studied, the viability of inbred offspring fully determines the relative change in the numbers of inbred and outbred offspring produced. In the second model, the relative change in offspring numbers depends in addition on associations between the S-locus and the viability locus. Because these two-locus associations are comparable in magnitude to the differences between the viabilities of inbred and outbred offspring, SSI can arise under less restrictive conditions than expected from the one-locus model. Greater allelic multiplicity at the viability locus facilitates the origin of SSI by reducing the relative viability of inbred offspring. Tight linkage between the S-locus and the viability locus and high rates of receipt of self-pollen promote the generation and maintenance of associations between the S-locus and the viability locus. In populations in which more than two viability alleles are maintained, the active S-allele can invade even in the absence of linkage with the viability locus. The present study establishes that incompatibility systems can arise in response to identity disequilibrium between a modifier of incompatibility and a locus subject to overdominant viability selection; in particular, compensation for the twofold cost of outcrossing does not require preexisting gametic level disequilibria.

Coevolution of the major histocompatibility complex and the t-complex in the mouse. I. Generation and maintenance of high complementarity associations

A quantitative model is developed to explore the effects of prezygotic and postzygotic incompatibility on the origin and maintenance of associations between the major histocompatibility complex (MHC) and the t-complex in the mouse. Incompatibility is represented by a reduction in the rate of conception or gestation of offspring derived from sperm bearing MHC antigens in common with the mother. Incompatibility encourages the evolution of associations from a state of complete independence between the two complexes by promoting the invasion of all novel antigens, including those that exhibit associations with the t-complex. Incompatibility can modify the relative numbers of antigens associated with each haplotype by actively promoting the exclusion or invasion of recombinants that bear formerly +-specific or t-specific antigens on the alternative haplotype. The results of the analysis indicate that the state of complete independence between the MHC and the t-complex is not preserved over evolutionary time in the presence of incompatibility. Further, the expression of incompatibility maintains fully associated states that include a single antigen associated with the t-haplotype and up to three to five antigens associated with the +-haplotype within a single population.

On the evolution of genetic incompatibility systems. IV. Modification of response to an existing antigen polymorphism under partial selfing

A 2-locus model of the evolution of self-incompatibility in a population practicing partial selfing is presented. An allele is introduced at a modifier locus which influences the strength of the rejection reaction expressed by the style in response to antigens recognized in pollen. Two causes of inbreeding depression are investigated. First, offspring viability depends solely on the source (self or non-self) of the fertilizing pollen. Second, offspring viability declines with the expression of recessive deleterious alleles, segregating at a third (disease) locus, which exhibit an imperfect association with antigen alleles. Evolutionary changes occurring at the disease locus are not considered in this study. The condition under which a modifier allele that intensifies the incompatibility reaction increases when rare depends upon the number of antigens, the frequency of recessive deleterious alleles at the disease locus, and the level of association between the antigen locus and the disease locus. It is the improvement of viability among offspring derived by outcrossing, rather than the prevention of self-fertilization, that may represent the primary evolutionary function of genetic incompatibility systems.

On the evolution of genetic incompatibility systems. III. Introduction of weak gametophytic self-incompatibility under partial inbreeding

I explore the proposition that genetic incompatibility systems serve as a means for parents to evaluate and discriminate among their own offspring. Conditions for the initial increase of gametophytic self-incompatibility in a self-compatible population undergoing selfing, sibmating, and random outcrossing are reported. The adaptive value of reducing the concordance between offspring and maternal genotypes depends upon the relative changes in the numbers of offspring derived by the three modes, parent-offspring relatedness, and the magnitude of distortion of transmission ratios through pollen. Recessivity of stylar expression and low rates of receipt of pollen from related individuals facilitate the evolution of self-incompatibility. Viewed as a means of preferential maternal investment in offspring of high quality, self-incompatibility may be regarded as serving a function in common with diverse phenomena, including sexual selection, brood reduction, and other forms of prezygotic and postzygotic incompatibility. Associations between incompatibility loci and loci expressing inbreeding depression are expected to improve the reliability of the level of concordance at incompatibility loci as a measure of genomic homozygosity and offspring quality.

Relative roles of mutation and selection in the maintenance of genetic variability

The extent and pattern of protein and DNA polymorphisms are discussed with emphasis on the mechanism of maintenance of the polymorphisms. Statistical studies suggest that a large proportion of genetic variability at the molecular level is maintained by a mutation-drift balance. At some loci, such as those for histocompatibility in mammals, however, a form of overdominant selection seems to be involved. In the presence of overdominant selection, polymorphic alleles may be maintained for tens of millions of years, so that the number of nucleotide differences between alleles is often very large, as in the case of self-incompatibility alleles in plants. There are also an increasing number of examples in which an adaptive change of a morphological or physiological character is caused by a single nucleotide substitution. Nevertheless, these mutations seem to be a small proportion of the total nucleotide changes that contribute to genetic variability and evolution. Although there are many examples of frequency-dependent selection, this form of selection is apparently unimportant for the maintenance of genetic variability except in some special cases. Observations on the evolutionary change of DNA suggest that the driving force of evolution is mutation rather than selection.

Linkage disequilibrium and gametophytic self-incompatibility

The approach to linkage equilibrium of a locus linked to the locus determining gametophytic self-incompatibility (S) is considered. For the simplest case of three alleles at the S locus and two at the linked locus it is necessary to consider 3 measures of linkage disequilibrium. These are found to approach their equilibrium value of zero in one of three ways: 1) steadily declining to zero; 2) oscillating as decline proceeds; 3) a combination: 2) followed by 1). Linkage equilibrium may be established before genotype frequencies reach their expectation under random crossing. Earlier studies (Li 1951; Moran 1962) of the approach to S allele equilibrium have been based on the assumption that all types of pollen take part in fertilizations equally frequently. Such an assumption leads to simpler expressions for changes in S gene frequencies but is extremely unrealistic and, in particular, leads to a different rate of approach to equilibrium from the more comprehensive model. It is shown that even in the absence of selection it is not possible to predict the equilibrium gene frequency of a linked locus until S allele equilibrium is reached. This frequency may be either higher or lower than that calculated from a gene count in the starting genotype pool. However, these two gene frequencies may stabilize long before linkage equilibrium is achieved. An examination of selection against one genotype at the linked locus is undertaken. If linkage is complete, lethality can be less effective at reducing the gene frequency than is less intense selection (in only a few generations of selection). Here too linkage equilibrium may be established with selection still effective in bringing about a decline in gene frequency. An examination of the analysis and conclusions of Rasmuson (1980) shows that because these were based on the inadequate formulae previously discussed and exclude phenomena discussed above, they are misleading. The possibility of a gametophytic self-incompatibility system providing a sufficient condition for the sheltering of lethals in the absence of the condition of complete linkage to the S locus (r=0) is shown to be unlikely.

Gametophytic Self-Incompatibility Reexamined

The conventional hypothesis of gametophytic self-incompatibility in the angiosperms involves one to four multiallelic incompatibility loci and the positive inhibition of incompatible pollen tubes. However, this concept does not accommodate recent experimental data indicating that there may be many loci. An alternative hypothesis which incorporates many loci and complementary pollen-style interactions suggests that there may be no S gene, as previously thought, and that gametophytic self-incompatibility is perhaps merely one aspect of extensive pollen-style interactions.