Recombination and selection at Brassica self-incompatibility loci

Creating novel specificities in a fungal nonself recognition system by single step homologous recombination events

In many organisms, two component systems have evolved to discriminate self from nonself. While the molecular function of the two components has been elucidated in several systems, the evolutionary events leading to the large number of different specificities for self-nonself recognition found in most systems remain obscure. We have investigated the variation within a multiallelic nonself recognition system in the phytopathogenic basidiomycete Ustilago maydis by means of sequence analysis and functional studies. The multiallelic b mating type locus of U. maydis ensures outbreeding during sexual development. Nonself recognition is specified by the two homeodomain proteins, bE and bW, encoded by the b locus. While bE-bW combinations from the same allele do not dimerize, bE and bW proteins originating from different alleles form a heterodimeric complex that functions as master regulator for sexual and pathogenic development. We show that novel specificities of the b mating type locus have arisen by single homologous recombination events between distinct b alleles that lead to a simultaneous exchange of subdomains involved in dimerization in both bE and bW, altering the specificity of both proteins in a single step.

Amino acid repeats cause extraordinary coding sequence variation in the social amoeba Dictyostelium discoideum

Protein sequences are normally the most conserved elements of genomes owing to purifying selection to maintain their functions. We document an extraordinary amount of within-species protein sequence variation in the model eukaryote Dictyostelium discoideum stemming from triplet DNA repeats coding for long strings of single amino acids. D. discoideum has a very large number of such strings, many of which are polyglutamine repeats, the same sequence that causes various human neurological disorders in humans, like Huntington’s disease. We show here that D. discoideum coding repeat loci are highly variable among individuals, making D. discoideum a candidate for the most variable proteome. The coding repeat loci are not significantly less variable than similar non-coding triplet repeats. This pattern is consistent with these amino-acid repeats being largely non-functional sequences evolving primarily by mutation and drift.

Inbreeding depression in self-incompatible North-American Arabidopsis lyrata: disentangling genomic and S-locus-specific genetic load

Newly formed selfing lineages may express recessive genetic load and suffer inbreeding depression. This can have a genome-wide genetic basis, or be due to loci linked to genes under balancing selection. Understanding the genetic architecture of inbreeding depression is important in the context of the maintenance of self-incompatibility and understanding the evolutionary dynamics of S-alleles. We addressed this using North-American subspecies of Arabidopsis lyrata. This species is normally self-incompatible and outcrossing, but some populations have undergone a transition to selfing. The goals of this study were to: (1) quantify the strength of inbreeding depression in North-American populations of A. lyrata; and (2) disentangle the relative contribution of S-linked genetic load compared with overall inbreeding depression. We enforced selfing in self-incompatible plants with known S-locus genotype by treatment with CO(2), and compared the performance of selfed vs outcrossed progeny. We found significant inbreeding depression for germination rate (δ=0.33), survival rate to 4 weeks (δ=0.45) and early growth (δ=0.07), but not for flowering rate. For two out of four S-alleles in our design, we detected significant S-linked load reflected by an under-representation of S-locus homozygotes in selfed progeny. The presence or absence of S-linked load could not be explained by the dominance level of S-alleles. Instead, the random nature of the mutation process may explain differences in the recessive deleterious load among lineages.

Creation and maintenance of variation in allorecognition Loci: molecular analysis in various model systems

Allorecognition is the ability of an organism to differentiate self or close relatives from unrelated conspecifics. Effective allorecognition systems are critical to the survival of organisms; they prevent inbreeding and facilitate fusions between close relatives. Where the loci governing allorecognition outcomes have been identified, the corresponding proteins often exhibit exceptional polymorphism. Two important questions about this polymorphism remain unresolved: how is it created, and how is it maintained. Because the genetic bases of several allorecognition systems have now been identified, including alr1 and alr2 in Hydractinia, fusion histocompatibility in Botryllus, the het (vic) loci in fungi, tgrB1 and tgrC1 in Dictyostelium, and self-incompatibility (SI) loci in several plant families, we are now poised to achieve a clearer understanding of how these loci evolve. In this review, we summarize what is currently known about the evolution of allorecognition loci, highlight open questions, and suggest future directions.

Genetic architecture of life history traits and environment-specific trade-offs

Life history theory predicts the evolution of trait combinations that enhance fitness, and the occurrence of trade-offs depends in part on the magnitude of variation in growth rate or acquisition. Using recombinant inbred lines, we examined the genetic architecture of age and size at reproduction across abiotic conditions encountered by cultivars and naturalized populations of Brassica rapa. We found that genotypes are plastic to seasonal setting, such that reproduction was accelerated under conditions encountered by summer annual populations and genetic variances for age at reproduction varied across simulated seasonal settings. Using an acquisition-allocation model, we predicted the likelihood of trade-offs. Consistent with predicted relationships, we observed a trade-off where early maturity is associated with small size at maturity under simulated summer and fall annual conditions but not under winter annual conditions. The trade-off in the summer annual setting was observed despite significant genotypic variation in growth rate, which is often expected to decouple age and size at reproduction because rapidly growing genotypes could mature early and attain a larger size relative to slowly growing genotypes that mature later. The absence of a trade-off in the winter setting is presumably attributable to the absence of genotypic differences in age at reproduction. We observed QTL for age at reproduction that jointly regulated size at reproduction in both the summer and fall annual settings, but these QTL were environment-specific (i.e. different QTL contributed to the trade-off in the fall vs. summer annual settings). Thus, at least some of the genetic mechanisms underlying observed trade-offs differed across environments.

Sporophytic self-incompatibility genes and mating system variation in Arabis alpina

BACKGROUND AND AIMS

Sporophytic self-incompatibility (SI) prevents inbreeding in many members of the Brassicaceae, and has been well documented in a variety of high-profile species. Arabis alpina is currently being developed as a model system for studying the ecological genetics of arctic-alpine environments, and is the focus of numerous studies on population structure and alpine phylogeography. Although it is highly inbreeding throughout most of its range, populations in central Italy have been identified that show inbreeding coefficients (F(IS)) more typical of self-incompatible relatives. The purpose of this study was to establish whether this variation is due to a functioning SI system.

METHODS

Outcrossing rate estimates were calculated based on 16 allozyme loci and self-compatibility assessed based on controlled pollinations for six Italian populations that have previously been shown to vary in F(IS) values. Putative SRK alleles (the gene controlling the female component of SI in other Brassicaceae) amplified from A. alpina were compared with those published for other species. Linkage of putative SRK alleles and SI phenotypes was assessed using a diallel cross.

KEY RESULTS

Functional avoidance of inbreeding is demonstrated in three populations of A. alpina, corresponding with previous F(IS) values. The presence is described of 15 putative SRK-like alleles, which show high sequence identity to known alleles from Brassica and Arabidopsis and the high levels of synonymous and nonsynonymous variation typical of genes under balancing selection. Also, orthologues of two other members of the S-receptor kinase gene family, Aly8 (ARK3) and Aly9 (AtS1) are identified. Further to this, co-segregation between some of the putative S-alleles and compatibility phenotypes was demonstrated using a full-sibling cross design.

CONCLUSIONS

The results strongly suggest that, as with other species in the Brassicaceae, A. alpina has a sporophytic SI system but shows variation in the strength of SI within and between populations.

Genetic linkage map for Amylostereum areolatum reveals an association between vegetative growth and sexual and self-recognition

Amylostereum areolatum is a filamentous fungus that grows through tip extension, branching and hyphal fusion. In the homokaryotic phase, the hyphae of different individuals are capable of fusing followed by heterokaryon formation, only if they have dissimilar allelic specificities at their mating-type (mat) loci. In turn, hyphal fusion between heterokaryons persists only when they share the same alleles at all of their heterokaryon incompatibility (het) loci. In this study we present the first genetic linkage map for A. areolatum, onto which the mat and het loci, as well as quantitative trait loci (QTLs) for mycelial growth rate are mapped. The recognition loci (mat-A and het-A) are positioned near QTLs associated with mycelial growth, suggesting that the genetic determinants influencing recognition and growth rate in A. areolatum are closely associated. This was confirmed when isolates associated with specific mat and het loci displayed significantly different mycelial growth rates. Although the link between growth and sexual recognition has previously been observed in other fungi, this is the first time that an association between growth and self-recognition has been shown.

Nonallelic interactions between het-c and a polymorphic locus, pin-c, are essential for nonself recognition and programmed cell death in Neurospora crassa

Nonself recognition in filamentous fungi is conferred by genetic differences at het (heterokaryon incompatibility) loci. When individuals that differ in het specificity undergo hyphal fusion, the heterokaryon undergoes a programmed cell death reaction or is highly unstable. In Neurospora crassa, three allelic specificities at the het-c locus are conferred by a highly polymorphic domain. This domain shows trans-species polymorphisms indicative of balancing selection, consistent with the role of het loci in nonself recognition. We determined that a locus closely linked to het-c, called pin-c (partner for incompatibility with het-c) was required for het-c nonself recognition and heterokaryon incompatibility (HI). The pin-c alleles in isolates that differ in het-c specificity were extremely polymorphic. Heterokaryon and transformation tests showed that nonself recognition was mediated by synergistic nonallelic interactions between het-c and pin-c, while allelic interactions at het-c increased the severity of the HI phenotype. The pin-c locus encodes a protein containing a HET domain; predicted proteins containing HET domains are frequent in filamentous ascomycete genomes. These data suggest that nonallelic interactions may be important in nonself recognition in filamentous fungi and that proteins containing a HET domain may be a key factor in these interactions.

RPW8 and resistance to powdery mildew pathogens in natural populations of Arabidopsis lyrata

It is not clear to what extent the orthologues of genes that are adaptively important in one species also contribute to adaptive variation in others. Here, we examine Arabidopsis lyrata to assess the functional and evolutionary significance of natural variation in an orthologue of the gene RPW8 known to be a major determinant of powdery mildew resistance in Arabidopsis thaliana. We assessed the sequence variation at RPW8 and the associated resistance reaction in populations of A. lyrata ssp. petraea. Neutrality tests were performed to understand the importance of local adaptation in maintaining variation at the locus. Highly truncated RPW8 proteins were frequent in all populations and were associated with an increased risk of susceptibility. Haplotypes encoding full-length proteins were highly significantly associated with resistance. There were no signatures of selection at the species-wide level, but some evidence for positive selection in two populations. RPW8 in A. lyrata appears to have a role in powdery mildew resistance, similar to its orthologue in A. thaliana. Unlike A. thaliana, A. lyrata contains a genetic component that can act independently of RPW8 to confer resistance to powdery mildew pathogens. Infrequent local selective sweeps may favour different alleles in different populations, and thereby contribute to the maintenance of species-wide variation at the locus.

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.

In vivo detection of residues required for ligand-selective activation of the S-locus receptor in Arabidopsis

The self-incompatibility response of crucifers is a barrier to fertilization in which arrest of pollen tube development is mediated by allele-specific interactions between polymorphic receptors and ligands encoded by the S-locus haplotype. Activation of stigma-expressed S-locus receptor kinase (SRK) [1] by pollen coat-localized S-locus cysteine-rich (SCR) ligand [2-5] and the resulting rejection of pollen occurs only if receptor and ligand are encoded by the same S haplotype [4, 6-8]. To identify residues within the SRK extracellular domain (eSRK) that are required for its ligand-selective activation, we assayed chimeric receptors and receptor variants containing substitutions at polymorphic sites in Arabidopsis thaliana[9, 10]. We show that only a small number of the approximately 100 polymorphic residues in eSRK are required for ligand-specific activation of self-incompatibility in vivo. These essential residues occur in two noncontiguous clusters located at equivalent positions in the two variants tested. They also correspond to sites showing elevated levels of substitutions in other SRKs, suggesting that these residues could define self-incompatibility specificity in most SRKs. The results demonstrate that the majority of eSRK residues that show signals of positive selection and previously surmised to function as specificity determinants are not essential for specificity in the SRK-SCR interaction.

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.

The evolution and diversification of S-locus haplotypes in the Brassicaceae family

Self-incompatibility (SI) in the Brassicaceae plant family is controlled by the SRK and SCR genes situated at the S locus. A large number of S haplotypes have been identified, mainly in cultivated species of the Brassica and Raphanus genera, but recently also in wild Arabidopsis species. Here, we used DNA sequences from the SRK and SCR genes of the wild Brassica species Brassica cretica, together with publicly available sequence data from other Brassicaceae species, to investigate the evolutionary relationships among S haplotypes in the Brassicaceae family. The results reveal that wild and cultivated Brassica species have similar levels of SRK diversity, indicating that domestication has had but a minor effect on S-locus diversity in Brassica. Our results also show that a common set of S haplotypes was present in the ancestor of the Brassica and Arabidopsis genera, that only a small number of haplotypes survived in the Brassica lineage after its separation from Arabidopsis, and that diversification within the two Brassica dominance classes occurred after the split between the two lineages. We also found indications that recombination may have occurred between the kinase domain of SRK and the SCR gene in Brassica.

Recombination and nucleotide diversity in the sex chromosomal pseudoautosomal region of the emu, Dromaius novaehollandiae

Pseudoautosomal regions (PARs) shared by avian Z and W sex chromosomes are typically small homologous regions within which recombination still occurs and are hypothesized to share the properties of autosomes. We capitalized on the unusual structure of the sex chromosomes of emus, Dromaius novaehollandiae, which consist almost entirely of PAR shared by both sex chromosomes, to test this hypothesis. We compared recombination, linkage disequilibrium (LD), GC content, and nucleotide diversity between pseudoautosomal and autosomal loci derived from 11 emu bacterial artificial chromosome (BAC) clones that were mapped to chromosomes by fluorescent in situ hybridization. Nucleotide diversity (pi = 4N(e)mu) was not significantly lower in pseudoautosomal loci (14 loci, 1.9 +/- 2.4 x 10(-3)) than autosomal loci (8 loci, 4.2 +/- 6.1 x 10(-3)). By contrast, recombination per site within BAC-end sequences (rho = 4Nc) (pseudoautosomal, 3.9 +/- 6.9 x 10(-2); autosomal, 2.3 +/- 3.7 x 10(-2)) was higher and average LD (D') (pseudoautosomal, 4.2 +/- 0.2 x 10(-1); autosomal, 4.7 +/- 0.5 x 10(-1)) slightly lower in pseudoautosomal sequences. We also report evidence of deviation from a simple neutral model in the PAR and in autosomal loci, possibly caused by departures from demographic equilibrium, such as population growth. This study provides a snapshot of the population genetics of avian sex chromosomes at an early stage of differentiation.

Patterns of divergence among conifer ESTs and polymorphism in Pinus sylvestris identify putative selective sweeps

Finding genes that are under positive selection is a difficult task, especially in non-model organisms. Here, we have analyzed expressed sequence tag (EST) data from 4 species (Pinus pinaster, Pinus taeda, Picea glauca, and Pseudotsuga menziesii) to investigate selection patterns during their evolution and to identify genes likely to be under positive selection. To confirm selection, population samples of these genes have been sequenced in Pinus sylvestris, a species that was not included in the EST data set. The estimates of branch-specific Ka/Ks (nonsynonymous/synonymous substitution rates) across all genes in the EST data set were similar or smaller than estimates from other higher plant species. There was no evidence for the traditional indication of positive selection, Ka/Ks above 1. However, several lines of evidence based on polymorphism patterns suggest that genes with high Ka/Ks (0.20-0.52) in the EST data set are in fact more affected by positive selection in P. sylvestris than genes with low Ka/Ks (0.01-0.04). The high Ka/Ks genes have a lower level of polymorphism and more negative Tajima's D than the low Ka/Ks genes. Further, in the high Ka/Ks group, the Hudson-Kreitman-Aguade test is significant. This suggests that the EST data set is a good starting point for finding genes under positive selection in conifers and that even moderate Ka/Ks values could be indicative of selection. A group of 5 genes with high Ka/Ks collectively show evidence for positive selection within P. sylvestris.

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.

Gene organization of the mating type regions in the ectomycorrhizal fungus Laccaria bicolor reveals distinct evolution between the two mating type loci

In natural conditions, basidiomycete ectomycorrhizal fungi such as Laccaria bicolor are typically in the dikaryotic state when forming symbioses with trees, meaning that two genetically different individuals have to fuse or 'mate'. Nevertheless, nothing is known about the molecular mechanisms of mating in these ecologically important fungi. Here, advantage was taken of the first sequenced genome of the ectomycorrhizal fungus, Laccaria bicolor, to determine the genes that govern the establishment of cell-type identity and orchestrate mating. The L. bicolor mating type loci were identified through genomic screening. The evolutionary history of the genomic regions that contained them was determined by genome-wide comparison of L. bicolor sequences with those of known tetrapolar and bipolar basidiomycete species, and by phylogenetic reconstruction of gene family history. It is shown that the genes of the two mating type loci, A and B, are conserved across the Agaricales, but they are contained in regions of the genome with different evolutionary histories. The A locus is in a region where the gene order is under strong selection across the Agaricales. By contrast, the B locus is in a region where the gene order is likely under a low selection pressure but where gene duplication, translocation and transposon insertion are frequent.

Genetic causes and consequences of the breakdown of self-incompatibility: case studies in the Brassicaceae

The genetic consequences of inbreeding is a subject that has received thorough theoretical attention and has been of interest to empirical biologists since the time of Darwin. Particularly for species with genetically controlled mechanisms to promote outcrossing (self-incompatibility or SI systems), it is expected that high levels of genetic load should accumulate through sheltering of deleterious recessive mutations. Nevertheless, transitions to selfing are common across angiosperms, which suggests that the potentially negative consequences of reduced heterozygosity and genetic diversity are balanced by other factors, such as reproductive assurance. This mini-review focuses on empirical research in the Brassicaceae to emphasize some of the genetic consequences of shifts to inbreeding in terms of mechanisms for loss of SI, changes in genetic diversity following loss of SI, and inbreeding depression in relation to outcrossing history. Despite the long history of theoretical attention, there are still some surprisingly large gaps in our understanding in each of these areas. Rather than providing a complete overview, examples are drawn predominantly from published and emerging data from Arabidopsis thaliana and its relatives to highlight recent progress and remaining questions. We are currently on the brink of major breakthroughs in understanding due both to advances in sequencing technology and a shift in focus from crop plants to natural populations, where critical factors such as population structure, phylogeography, demographic history, partial compatibility and individual variation can be taken into account when investigating the nature of the selective forces regulating mating system evolution.

Hitch-hiking to a locus under balancing selection: high sequence diversity and low population subdivision at the S-locus genomic region inArabidopsis halleri

Hitch-hiking to a site under balancing selection is expected to produce a local increase in nucleotide polymorphism and a decrease in population differentiation compared with the background genomic level, but empirical evidence supporting these predictions is scarce. We surveyed molecular diversity at four genes flanking the region controlling self-incompatibility (the S-locus) in samples from six populations of the herbaceous plantArabidopsis halleri, and compared their polymorphism with sequences from five control genes unlinked to the S-locus. As a preliminary verification, the S-locus flanking genes were shown to co-segregate withSRK, the gene involved in the self-incompatibility reaction at the pistil level. In agreement with theory, our results demonstrated a significant peak of nucleotide diversity around the S-locus as well as a significant decrease in population genetic structure in the S-locus region compared with both control genes and a set of seven unlinked microsatellite markers. This is consistent with the theoretical expectation that balancing selection is increasing the effective migration rate in subdivided populations. Although only four S-locus flanking genes were investigated, our results suggest that these two signatures of the hitch-hiking effect are localized in a very narrow genomic region.

Detecting groups of coevolving positions in a molecule: a clustering approach

Background

Although the patterns of co-substitutions in RNA is now well characterized, detection of coevolving positions in proteins remains a difficult task. It has been recognized that the signal is typically weak, due to the fact that (i) amino-acid are characterized by various biochemical properties, so that distinct amino acids changes are not functionally equivalent, and (ii) a given mutation can be compensated by more than one mutation, at more than one position.

Results

We present a new method based on phylogenetic substitution mapping. The two above-mentioned problems are addressed by (i) the introduction of a weighted mapping, which accounts for the biochemical effects (volume, polarity, charge) of amino-acid changes, (ii) the use of a clustering approach to detect groups of coevolving sites of virtually any size, and (iii) the distinction between biochemical compensation and other coevolutionary mechanisms. We apply this methodology to a previously studied data set of bacterial ribosomal RNA, and to three protein data sets (myoglobin of vertebrates, S-locus Receptor Kinase and Methionine Amino-Peptidase).

Conclusion

We succeed in detecting groups of sites which significantly depart the null hypothesis of independence. Group sizes range from pairs to groups of size ≃ 10, depending on the substitution weights used. The structural and functional relevance of these groups of sites are assessed, and the various evolutionary processes potentially generating correlated substitution patterns are discussed.

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.

Effects of recombination on hitchhiking diversity in the Brassica self-incompatibility locus complex

In self-incompatibility, a number of S haplotypes are maintained by frequency-dependent selection, which results in trans-specific S haplotypes. The region of several kilobases (approximately 40-60 kb) from SP6 to SP2, including self-incompatibility-related genes and some adjacent genes in Brassica rapa, has high nucleotide diversity due to the hitchhiking effect, and therefore we call this region the "S-locus complex." Recombination in the S-locus complex is considered to be suppressed. We sequenced regions of >50 kb of the S-locus complex of three S haplotypes in B. rapa and found higher nucleotide diversity in intergenic regions than in coding regions. Two highly similar regions of >10 kb were found between BrS-8 and BrS-46. Phylogenetic analysis using trans-specific S haplotypes (called interspecific pairs) of B. rapa and B. oleracea suggested that recombination reduced the nucleotide diversity in these two regions and that the genes not involved in self-incompatibility in the S-locus complex and the kinase domain, but not the S domain, of SRK have also experienced recombination. Recombination may reduce hitchhiking diversity in the S-locus complex, whereas the region from the S domain to SP11 would disfavor recombination.

Pronounced differences of recombination activity at the sex determination locus of the honeybee, a locus under strong balancing selection

Recombination decreases the association of linked nucleotide sites and can influence levels of polymorphism in natural populations. When coupled with selection, recombination may relax potential conflict among linked genes, a concept that has played a central role in research on the evolution of recombination. The sex determination locus (SDL) of the honeybee is an informative example for exploring the combined forces of recombination, selection, and linkage on sequence evolution. Balancing selection at SDL is very strong and homozygous individuals at SDL are eliminated by worker bees. The recombination rate is increased up to four times that of the genomewide average in the region surrounding SDL. Analysis of nucleotide diversity (pi) reveals a sevenfold increase of polymorphism within the sex determination gene complementary sex determiner (csd) that rapidly declines within 45 kb to levels of genomewide estimates. Although no recombination was observed within SDL, which contains csd, analyses of heterogeneity, shared polymorphic sites, and linkage disequilibrium (LD) show that recombination has contributed to the evolution of the 5' part of some csd sequences. Gene conversion, however, has not obviously contributed to the evolution of csd sequences. The local control of recombination appears to be related to SDL function and mode of selection. The homogenizing force of recombination is reduced within SDL, which preserves allelic differences and specificity, while the increase of recombination activity around SDL relaxes conflict between SDL and linked genes.

Comparison of the genome structure of the self-incompatibility (S) locus in interspecific pairs of S haplotypes

The determinants of recognition specificity of self-incompatibility in Brassica are SRK in the stigma and SP11/SCR in the pollen, both of which are encoded in the S locus. The nucleotide sequence analyses of many SRK and SP11/SCR alleles have identified several interspecific pairs of S haplotypes having highly similar sequences between B. oleracea and B. rapa. These interspecific pairs of S haplotypes are considered to be derived from common ancestors and to have maintained the same recognition specificity after speciation. In this study, the genome structures of three interspecific pairs of S haplotypes were compared by sequencing SRK, SP11/SCR, and their flanking regions. Regions between SRK and SP11/SCR in B. oleracea were demonstrated to be much longer than those of B. rapa and several retrotransposon-like sequences were identified in the S locus in B. oleracea. Among the seven retrotransposon-like sequences, six sequences were found to belong to the ty3 gypsy group. The gag sequences of the retrotransposon-like sequences were phylogenetically different from each other. In Southern blot analysis using retrotransposon-like sequences as probes, the B. oleracea genome showed more signals than the B. rapa genome did. These findings suggest a role for the S locus and genome evolution in self-incompatible plant species.

Phenotypic variation and natural selection at catsup, a pleiotropic quantitative trait gene in Drosophila

Quantitative traits are shaped by networks of pleiotropic genes . To understand the mechanisms that maintain genetic variation for quantitative traits in natural populations and to predict responses to artificial and natural selection, we must evaluate pleiotropic effects of underlying quantitative trait genes and define functional allelic variation at the level of quantitative trait nucleotides (QTNs). Catecholamines up (Catsup), which encodes a negative regulator of tyrosine hydroxylase , the rate-limiting step in the synthesis of the neurotransmitter dopamine, is a pleiotropic quantitative trait gene in Drosophila melanogaster. We used association mapping to determine whether the same or different QTNs at Catsup are associated with naturally occurring variation in multiple quantitative traits. We sequenced 169 Catsup alleles from a single population and detected 33 polymorphisms with little linkage disequilibrium (LD). Different molecular polymorphisms in Catsup are independently associated with variation in longevity, locomotor behavior, and sensory bristle number. Most of these polymorphisms are potentially functional variants in protein coding regions, have large effects, and are not common. Thus, Catsup is a pleiotropic quantitative trait gene, but individual QTNs do not have pleiotropic effects. Molecular population genetic analyses of Catsup sequences are consistent with balancing selection maintaining multiple functional polymorphisms.

Duplication, balancing selection and trans-species evolution explain the high levels of polymorphism of the DQA MHC class II gene in voles (Arvicolinae)

Major histocompatibility complex (MHC) genes play important role in host-parasite interactions and parasites are crucial factors influencing the population dynamics of hosts. We described the structure and diversity of exon 2 of the MHC class II DQA gene in three species of voles (Arvicolinae) exhibiting regular multi-annual fluctuations of population density and analysed the processes leading to the observed MHC polymorphism. By using cloning-sequencing methodology and capillary electrophoresis-single strand conformation polymorphism, we described seven sequences in the water, eight in the common, and seven in the bank voles coming from an area of 70 km(2) around the Nozeroy canton in the Jura Mountains (Franche Comté, France). All exon 2 sequences translate to give unique amino acid sequences and positive selection was found to act very intensively on antigen binding sites. We documented the presence of recombination at vole DQA region but its importance in generating allelic polymorphism seems to be relatively limited. For the first time within rodents, we documented the duplication of the DQA gene in all three species with both copies being transcriptionally active. Phylogenetic analysis of allelic sequences revealed extensive trans-species polymorphism within the subfamily although no alleles were shared between species in our data set. We discuss possible role of parasites in forming the recent polymorphism pattern of the DQA locus in voles.

The impact of intragenic recombination on phylogenetic reconstruction at the sectional level in Eucalyptus when using a single copy nuclear gene (cinnamoyl CoA reductase)

Low copy number nuclear genes have been found to be useful for phylogenetic reconstruction at different taxonomic levels. This study investigated the utility of a single copy gene, cinnamoyl CoA reductase (CCR), for resolving phylogenetic relationships at the sectional level within Eucalyptus. The monophyly of sections Exsertaria and Latoangulatae was explored, using section Maidenaria as an outgroup, and the impact of intragenic recombination on phylogenetic reconstruction examined. Phylogenetic analysis did not resolve monophyletic groups. Latoangulatae and Maidenaria were polyphyletic or paraphyletic. Exsertaria species formed a clade but included a single Latoangulatae species (E. major). Recombination analysis identified two intragenic recombination events that involved species from different sections, which have probably been facilitated by inter-sectional hybridisation. One of the events most likely occurred prior to speciation, with several Latoangulatae species having the recombinant allele. The other event may have occurred after speciation, since only one of two E. globulus samples possessed the recombinant allele. This is the first detailed report of intragenic recombination in both CCR and Eucalyptus, and between species of different sections of a plant genus. The occurrence of intragenic recombination may explain the anomalous positions of some species within the phylogenetic tree, and indicates that phylogenetic analysis of Eucalyptus using nuclear genes will be problematic unless recombination is taken into account.

Evidence for recombination in scorpion mitochondrial DNA (Scorpiones: Buthidae)

There has been very little undisputed evidence for recombination in animal mitochondrial DNA (mtDNA) provided so far. Previous unpublished results suggestive of mtDNA recombination in the scorpion family Buthidae, together with cytological evidence for a unique mechanism of mitochondrial fusion in that family, prompted us to investigate this group in more details. First, we sequenced the complete mtDNA genome of Mesobuthus gibbosus, and chose two genes opposing each other (16S and coxI). We then sequenced 150 individuals from the natural populations of four species of Buthidae (Old World genera Buthus and Mesobuthus). We observed strong evidence for widespread recombination through highly significant negative correlations between linkage disequilibrium and physical distance in three out of four species. The evidence is further confirmed when using five other tests for recombination and by the presence of a high amount of homoplasy in phylogenetic trees.

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.

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.

Molecular Population Genetics and the Search for Adaptive Evolution in Plants

The first papers on plant molecular population genetics were published approximately 10 years ago. Since that time, well over 50 additional studies of plant nucleotide polymorphism have been published, and many of these studies focused on detecting the signature of balancing or positive selection at a locus. In this review, we discuss some of the theoretical and statistical issues surrounding the detection of selection, with focus on plant populations, and we also summarize the empirical plant molecular population genetics literature. At face value, the literature suggests that a history of balancing or positive selection in plant genes is rampant. In two well-studied taxa (maize and Arabidopsis) over 20% of studied genes have been interpreted as containing the signature of selection. We argue that this is probably an overstatement of the prevalence of natural selection in plant genomes, for two reasons. First, demographic effects are difficult to incorporate and have generally not been well integrated into the plant population genetics literature. Second, the genes studied to date are not a random sample, so selected genes may be overrepresented. The next generation of studies in plant molecular population genetics requires additional sampling of local populations, explicit comparisons among loci, and improved theoretical methods to control for demography. Eventually, candidate loci should be confirmed by explicit consideration of phenotypic effects.

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.

Low linkage disequilibrium indicative of recombination in foot-and-mouth disease virus gene sequence alignments

We have applied tests for detecting recombination to genes of foot-and-mouth disease virus (FMDV). Our approach estimated summary statistics of linkage disequilibrium (LD), which are sensitive to recombination. Using the genealogical relationships, rate heterogeneity and mutation parameters estimated from individual sets of aligned gene sequences, we simulated matching RNA sequence datasets without recombination. These simulated datasets allowed for recurrent mutations at any site to mimic homoplasy in virus sequence data and allow construction of null distributions for LD parameters expected in the absence of recombination. We tested for recombination in two ways: by comparing LD in observed data with corresponding null distributions obtained from simulated data; and by testing for a negative relationship between observed LD between pairs of polymorphic nucleotide sites and inter-site distance. We applied these tests to six FMDV datasets from four serotypes and found some evidence for recombination in all of them.

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.

Molecular evolution of the s locus controlling mating in the brassicaceae

Flowering plants possess self-incompatibility (SI) mechanisms that promote outbreeding and thereby increase their genetic diversity. In the self-incompatible Brassicaceae, recognition and rejection of self-pollen is based on a receptor-ligand interaction between male and female SI determinants. A transmembrane receptor kinase (S locus Receptor Kinase, SRK) determines the SI specificity in stigmatic cells, whereas a pollen coat-localized ligand (S locus Cysteine-Rich, SCR) determines the SI specificity in pollen. During recent years, major advances have been made in the understanding of the molecular basis of self-pollen recognition by stigmatic cells. In this review, we will focus on evolutionary aspects of the SI system in Brassicaceae. We will describe how the study of the molecular aspect of SI, not only in the historical Brassica model but also in Arabidopsis species, has contributed to highlight certain aspects of evolution of SI in the Brassicaceae.

Evidence for rare recombination at the gametophytic self-incompatibility locus

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

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.

The evolutionary genomics of pathogen recombination

A pressing problem in studying the evolution of microbial pathogens is to determine the extent to which these genomes recombine. This information is essential for locating pathogenicity loci by using association studies or population genetic approaches. Recombination also complicates the use of phylogenetic approaches to estimate evolutionary parameters such as selection pressures. Reliable methods that detect and estimate the rate of recombination are, therefore, vital. This article reviews the approaches that are available for detecting and estimating recombination in microbial pathogens and how they can be used to understand pathogen evolution and to identify medically relevant loci.

Coevolution of the S-locus genes SRK, SLG and SP11/SCR in Brassica oleracea and B. rapa

Brassica self-incompatibility (SI) is controlled by SLG and SRK expressed in the stigma and by SP11/SCR expressed in the anther. We determined the sequences of the S domains of 36 SRK alleles, 13 SLG alleles, and 14 SP11 alleles from Brassica oleracea and B. rapa. We found three S haplotypes lacking SLG genes in B. rapa, confirming that SLG is not essential for the SI recognition system. Together with reported sequences, the nucleotide diversities per synonymous and nonsynonymous site (pi(S) and pi(N)) at the SRK, SLG, and SP11 loci within B. oleracea were computed. The ratios of pi(N):pi(S) for SP11 and the hypervariable region of SRK were significantly >1, suggesting operation of diversifying selection to maintain the diversity of these regions. In the phylogenetic trees of 12 SP11 sequences and their linked SRK alleles, the tree topology was not significantly different between SP11 and SRK, suggesting a tight linkage of male and female SI determinants during the evolutionary course of these haplotypes. Genetic exchanges between SLG and SRK seem to be frequent; three such recent exchanges were detected. The evolution of S haplotypes and the effect of gene conversion on self-incompatibility are discussed.

Self-incompatibility: how to stay incompatible

Plant self-incompatibility is controlled by different genes for the recognition reactions of pollen and stigmas, yet correct association of the two genes have been maintained in two Brassica species.

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.

Recombination, balancing selection and phylogenies in MHC and self-incompatibility genes

Using a coalescent model of multiallelic balancing selection with recombination, the genealogical process as a function of recombinational distance from a site under selection is investigated. We find that the shape of the phylogenetic tree is independent of the distance to the site under selection. Only the timescale changes from the value predicted by Takahata's allelic genealogy at the site under selection, converging with increasing recombination to the timescale of the neutral coalescent. However, if nucleotide sequences are simulated over a recombining region containing a site under balancing selection, a phylogenetic tree constructed while ignoring such recombination is strongly affected. This is true even for small rates of recombination. Published studies of multiallelic balancing selection, i.e., the major histocompatibility complex (MHC) of vertebrates, gametophytic and sporophytic self-incompatibility of plants, and incompatibility of fungi, all observe allelic genealogies with unexpected shapes. We conclude that small absolute levels of recombination are compatible with these observed distortions of the shape of the allelic genealogy, suggesting a possible cause of these observations. Furthermore, we illustrate that the variance in the coalescent with recombination process makes it difficult to locate sites under selection and to estimate the selection coefficient from levels of variability.

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.

Intrahaplotype polymorphism at the Brassica S locus

The S locus receptor kinase and the S locus glycoproteins are encoded by genes located at the S locus, which controls the self-incompatibility response in Brassica. In class II self-incompatibility haplotypes, S locus glycoproteins can be encoded by two different genes, SLGA and SLGB. In this study, we analyzed the sequences of these genes in several independently isolated plants, all of which carry the same S haplotype (S(2)). Two groups of S(2) haplotypes could be distinguished depending on whether SRK was associated with SLGA or SLGB. Surprisingly, SRK alleles from the two groups could be distinguished at the sequence level, suggesting that recombination rarely occurs between haplotypes of the two groups. An analysis of the distribution of polymorphisms along the S domain of SRK showed that hypervariable domains I and II tend to be conserved within haplotypes but to be highly variable between haplotypes. This is consistent with these domains playing a role in the determination of haplotype specificity.

Identification and characterization of a polymorphic receptor kinase gene linked to the self-incompatibility locus of Arabidopsis lyrata

We study the segregation of variants of a putative self-incompatibility gene in Arabidopsis lyrata. This gene encodes a sequence that is homologous to the protein encoded by the SRK gene involved in self-incompatibility in Brassica species. We show by diallel pollinations of plants in several full-sib families that seven different sequences of the gene in A. lyrata are linked to different S-alleles, and segregation analysis in further sibships shows that four other sequences behave as allelic to these. The family data on incompatibility provide evidence for dominance classes among the S-alleles, as expected for a sporophytic SI system. We observe no division into pollen-dominant and pollen-recessive classes of alleles as has been found in Brassica, but our alleles fall into at least three dominance classes in both pollen and stigma expression. The diversity among sequences of the A. lyrata putative S-alleles is greater than among the published Brassica SRK sequences, and, unlike Brassica, the alleles do not cluster into groups with similar dominance.